伯努利微分方程求解器

伯努利微分方程求解器

伯努利微分方程求解器 能够解决最著名的一阶非线性微分方程之一 —— 伯努利方程 y' + P(x)y = Q(x)yⁿ —— 并将经典的教科书推导过程转化为交互式的分步演示。它通过代换 v = y¹⁻ⁿ 将方程线性化，构建积分因子 μ(x)，并将生成的闭式解曲线叠加在 RK4 数值解和斜率场上，让你能同时看到所有细节。

什么是伯努利微分方程？

伯努利方程由雅各布·伯努利于 1695 年提出，是形如下式的一阶 ODE：

当 n = 0 时，方程已经是线性的；当 n = 1 时，它是可分离变量的。对于任何其他实数 n，方程是非线性的，但经典的代换 v = y¹⁻ⁿ 可将其转化为关于 v 的一阶线性 ODE，从而可以使用标准的积分因子技巧求解。

伯努利方法的六个步骤

从 y' + P(x)y = Q(x)yⁿ 开始：

1. 除以 yⁿ： \( y^{-n}y' + P(x)\,y^{1-n} = Q(x) \)。
2. 代换 v = y¹⁻ⁿ： 注意 \( v' = (1-n)y^{-n}y' \)，因此 \( y^{-n}y' = v'/(1-n) \)。
3. 线性化： \( v' + (1-n)P(x)\,v = (1-n)Q(x) \) —— 一个关于 v 的一阶线性 ODE。
4. 积分因子： \( \mu(x) = \exp\!\left(\int (1-n)P(x)\,dx\right) \)，因此 \( (\mu v)' = \mu(1-n)Q(x) \)。
5. 求解 v(x)： \( v(x) = \frac{1}{\mu(x)}\left[\mu(x_0)v_0 + \int_{x_0}^{x}\mu(t)(1-n)Q(t)\,dt\right] \).
6. 回代： \( y(x) = v(x)^{1/(1-n)} \)。

当涉及的积分是初等函数时，你可以得到干净的闭式解；否则，计算器会使用辛普森法则进行数值评估，以绘制解曲线。

自动处理的特殊情况

计算示例 — n = 2，逻辑斯谛型

考虑 y' + y/x = x·y²，初始条件为 y(1) = 1。此处 P(x) = 1/x，Q(x) = x，n = 2，因此 1 − n = −1。

1. 代换 v = y⁻¹ = 1/y。则 v' = −y⁻²y'，方程变为 v' − (1/x)v = −x。
2. 积分因子：μ(x) = exp(∫−1/x dx) = 1/x。
3. (μ·v)' = μ·(−x) = −1。积分：(1/x)·v = −x + C，即 v = −x² + Cx。
4. 代入初始条件：在 x = 1 时，v = 1/1 = 1，因此 1 = −1 + C ⇒ C = 2。因此 v(x) = −x² + 2x。
5. 回代：y(x) = 1/v(x) = 1/(2x − x²) = 1/(x(2 − x))。

闭式解 y = 1/(x(2−x)) 在 x = 0 和 x = 2 处具有垂直渐近线 —— 这正是斜率场一眼就能看出的特征。

如何使用此计算器

1. 填写方程构建器。 在蓝色方框中输入 P(x) 和 Q(x)，在右上角的小框中输入指数 n。布局镜像了标准形式 y' + P(x)y = Q(x)yⁿ。
2. 设置初始条件 (x₀, y₀) 和绘图范围 [x min, x max]。范围应包含 x₀。
3. 点击求解。 计算器会检测是否属于特殊情况（n = 0 或 n = 1）并显示匹配的推导过程。否则，它将运行完整的六步伯努利代换，并使用 MathJax 渲染方程。
4. 查看图表。 橙色曲线是 RK4 数值解。蓝色虚线是通过积分因子评估的闭式解。箭头场显示了各处的 y'，因此你也可以直观观察其他解。
5. 复制采样点的 CSV 格式（如果你想将轨迹导入其他程序）。

技巧、陷阱与边缘情况

• 非整数 n 要求 y > 0。 当 n = 1/2 且 y₀ ≤ 0 时，求解器会标记此组合，因为 yⁿ 会变为复数。
• y₀ = 0 通常是奇异的。 对于 Q ≠ 0 且 n > 0 的伯努利方程，具有平凡解 y ≡ 0，但这通常不是你想要的分支。
• 避免 P(x) 在 x₀ 附近爆炸。 像 1/x 这样的表达式要求 x₀ ≠ 0；求解器在运行前会进行验证。
• 大指数 (|n| > 20) 会被拒绝以防止溢出。在实际应用中，很少出现 n 如此大的伯努利方程。
• 垂直渐近线。 如果 RK4 发散，请尝试将 x 范围缩小到 x₀ 的一侧，即解保持有限的一侧。

伯努利方程的应用场景

• 人口动态 — 逻辑斯谛方程 y' = ry(1 − y/K) 本质上是伯努利方程（整理后 n = 2）。
• 化学动力学 — 自催化反应通常遵循 y' ∝ y − y²。
• 电路分析 — 某些含有非线性电阻的 RL 电路会产生伯努利形式。
• 流体力学 — 相似性简化后的边界层方程。
• 流行病模型 — SIR 模型中的易感人群比例可以简化为伯努利形式。
• 经济增长 — 储蓄率恒定的索洛-斯旺模型（Solow–Swan model）是 n = α 的伯努利方程。

常见问题解答

什么是伯努利微分方程？

伯努利方程是一种形如 y' + P(x)y = Q(x)yⁿ 的一阶 ODE，其中 P 和 Q 是连续函数，n 是任何实数。它是非线性 ODE 的经典例子，可以通过代换 v = y¹⁻ⁿ 转化为线性方程。

代换 v = y¹⁻ⁿ 是如何工作的？

将原方程乘以 y⁻ⁿ，使每个 y 项变为 y¹⁻ⁿ 或 y⁻ⁿy'。设置 v = y¹⁻ⁿ 得到 v' = (1−n)y⁻ⁿy'。代换后将伯努利方程转变为 v' + (1−n)P(x)v = (1−n)Q(x)，该方程关于 v 是线性的，可以用积分因子求解。

当 n = 0 或 n = 1 时会发生什么？

当 n = 0 时，方程已经是一阶线性的，因此不需要代换。当 n = 1 时，伯努利公式会除以 1 − n = 0，因此我们单独处理：方程简化为 y' = (Q(x) − P(x))·y，这是一个可分离变量方程，其闭式解为 y = y₀·exp(∫(Q−P) dx)。

伯努利方程总能得到闭式解吗？

原则上是可以的，但涉及积分因子的最终积分可能没有初等原函数。在这种情况下，计算器会使用辛普森法则进行数值评估并绘制解曲线。该方法本身总是能将伯努利 ODE 简化为求积问题。

为什么负数 y 和非整数 n 会导致问题？

如果 n 不是整数，yⁿ 定义为 exp(n·ln y)，且仅在 y > 0 时为实数。输入负数 y 会产生复数。求解器会标记这种情况，并要求 y₀ > 0 或使用整数指数，以保持解为实值。

斜率场显示了什么？

斜率场是由微小切线段组成的网格，其角度等于该 (x, y) 点处的 y'。任何解曲线都必须遵循这些切线，因此斜率场让你能一次性看到所有解的定性形状，而初始条件则确定了那条特定的曲线。

延伸阅读

• 伯努利微分方程 — 维基百科
• 积分因子 — 维基百科
• 逻辑斯谛函数 — 维基百科
• 斜率场 — 维基百科