求解带双初始条件的微分方程$y'=y^2-x$及代换后求解疑问

嘿，你的变量替换思路其实是可行的，不过第一步推导时不小心漏了个系数2，导致后续方程出错了。我帮你纠正这个小失误，然后一步步完成求解，再应用给定的初始条件：

第一步：纠正变量替换后的方程 #

原方程是 $y' = y^2 - x$，你设 $t = y^2$，这个替换没问题。对 $t$ 求导得 $2yy' = t'$，所以 $y' = \frac{t'}{2\sqrt{t}}$。代入原方程时，应该得到：

\frac{t'}{2\sqrt{t}} = t - x

两边乘以 $2\sqrt{t}$，正确的方程是：

t' = 2t\sqrt{t} - 2x\sqrt{t} \quad \Rightarrow \quad t' = 2t^{3/2} - 2x t^{1/2}

你之前的推导里少了右边第二项的系数2，这是关键的小疏忽~

第二步：转化为线性微分方程 #

接下来我们把这个方程整理成更容易求解的形式：

• 两边除以 $t^{3/2}$（因为两个初始条件下y都不为0，所以t>0，这个操作合法）：

  \frac{t'}{t^{3/2}} = 2 - \frac{2x}{t^{1/2}}
  

• 设 $u = \frac{1}{t^{1/2}} = \frac{1}{y}$（因为 $t=y^2$，所以 $u=1/\sqrt{t}=1/y$），对u求导：

  u' = -\frac{1}{2}t^{-3/2}t' \quad \Rightarrow \quad \frac{t'}{t^{3/2}} = -2u'
  

• 代入上面的方程，整理后得到一阶线性非齐次微分方程：

  u' - x u = -1
  

第三步：用积分因子法求解线性方程 #

对于线性方程 $u' + P(x)u = Q(x)$，我们用积分因子法：

1. 计算积分因子：

   \mu(x) = \exp\left(\int -x dx\right) = \exp\left(-\frac{x^2}{2}\right)
   

2. 方程两边乘以积分因子，左边会变成乘积的导数：

   \exp\left(-\frac{x^2}{2}\right)u' - x\exp\left(-\frac{x^2}{2}\right)u = -\exp\left(-\frac{x^2}{2}\right)
   

   左边等价于 $\frac{d}{dx}\left[u \cdot \exp\left(-\frac{x^2}{2}\right)\right]$
3. 两边积分：

   u \cdot \exp\left(-\frac{x^2}{2}\right) = -\int \exp\left(-\frac{x^2}{2}\right)dx + C
   

   右边的积分无法用初等函数表示，我们用**误差函数（erf）**来表示：
   $\int \exp\left(-\frac{x^2}{2}\right)dx = \sqrt{\frac{\pi}{2}} \text{erf}\left(\frac{x}{\sqrt{2}}\right) + C$，其中误差函数定义为 $\text{erf}(z) = \frac{2}{\sqrt{\pi}}\int_0^z \exp(-t^2)dt$
4. 解出u：

   u = \exp\left(\frac{x^2}{2}\right)\left[ C - \sqrt{\frac{\pi}{2}} \text{erf}\left(\frac{x}{\sqrt{2}}\right) \right]
   

第四步：换回y并应用初始条件 #

因为 $u = 1/y$，所以：

y(x) = \frac{1}{\exp\left(\frac{x^2}{2}\right)\left[ C - \sqrt{\frac{\pi}{2}} \text{erf}\left(\frac{x}{\sqrt{2}}\right) \right]}

初始条件1：$y(0)=1$ #

当x=0时，$\exp(0)=1$，$\text{erf}(0)=0$，代入得：

1 = \frac{1}{C - 0} \quad \Rightarrow \quad C=1

对应的解为：

y(x) = \frac{1}{\exp\left(\frac{x^2}{2}\right)\left[ 1 - \sqrt{\frac{\pi}{2}} \text{erf}\left(\frac{x}{\sqrt{2}}\right) \right]}

初始条件2：$y(0)=0.5$ #

当x=0时，代入得：

0.5 = \frac{1}{C - 0} \quad \Rightarrow \quad C=2

对应的解为：

y(x) = \frac{1}{\exp\left(\frac{x^2}{2}\right)\left[ 2 - \sqrt{\frac{\pi}{2}} \text{erf}\left(\frac{x}{\sqrt{2}}\right) \right]}

补充说明 #

• 这个方程属于Riccati型微分方程，由于无法找到初等函数形式的特解，所以最终解需要用误差函数这种特殊函数来表示。如果需要数值结果，可以用数值方法（比如龙格-库塔法）结合初始条件计算近似值。
• 注意解的定义域：当分母为0时，解会出现奇点，比如对于y(0)=1的解，当 $1 - \sqrt{\frac{\pi}{2}} \text{erf}\left(\frac{x}{\sqrt{2}}\right)=0$ 时，对应的x值就是解的奇点，超过这个值解就不存在了。

内容的提问来源于stack exchange，提问作者Inverse Problem