求解一阶非线性微分方程:特定方程解法与经典方法咨询
嘿,咱们一步步来解这个微分方程吧!你之前试了(u = yx)的代换没成功很正常,换个视角重新整理方程就能找到更清晰的解法啦。
原方程为:
$$y dx+(xy+2x-ye^y) dy=0$$
先把方程整理成关于(\frac{dx}{dy})的形式(两边除以(dy),假设(y \neq 0)):
$$y \frac{dx}{dy} + xy + 2x - y e^y = 0$$
把含(x)的项合并,移项后得到标准的一阶线性非齐次微分方程(以(x)为因变量,(y)为自变量):
$$\frac{dx}{dy} + x\left(1 + \frac{2}{y}\right) = e^y$$
这个形式符合(\frac{dx}{dy} + P(y)x = Q(y)),其中(P(y) = 1 + \frac{2}{y}),(Q(y) = e^y),非常适合用积分因子法求解。
积分因子(\mu(y))的计算公式为:
$$\mu(y) = \exp\left(\int P(y) dy\right)$$
先计算积分部分:
$$\int \left(1 + \frac{2}{y}\right) dy = y + 2\ln|y| + C$$
忽略常数项,取最简形式的积分因子:
$$\mu(y) = \exp\left(y + 2\ln|y|\right) = e^y \cdot y^2$$
将原方程两边乘以积分因子(\mu(y)),左边会变成乘积的导数形式:
$$y^2 e^y \frac{dx}{dy} + y^2 e^y \left(1 + \frac{2}{y}\right)x = y^2 e^y \cdot e^y$$
$$\frac{d}{dy}\left(y^2 e^y x\right) = y^2 e^{2y}$$
接下来对两边积分:
$$\int \frac{d}{dy}\left(y^2 e^y x\right) dy = \int y^2 e^{2y} dy$$
左边积分结果直接是(y^2 e^y x),右边需要用分部积分法计算:
设(u = y^2),(dv = e^{2y} dy),则(du = 2y dy),(v = \frac{1}{2}e^{2y}),第一次分部得:
$$\frac{1}{2}y^2 e^{2y} - \int y e^{2y} dy$$
对(\int y e^{2y} dy)再次分部:设(u = y),(dv = e^{2y} dy),则(du = dy),(v = \frac{1}{2}e^{2y}),得到:
$$\frac{1}{2}y e^{2y} - \frac{1}{4}e^{2y} + C$$
把结果合并,右边积分最终为:
$$\frac{1}{2}y^2 e^{2y} - \frac{1}{2}y e^{2y} + \frac{1}{4}e^{2y} + C$$
最后解出(x):
$$x = \frac{\frac{1}{2}y^2 e^{2y} - \frac{1}{2}y e^{2y} + \frac{1}{4}e^{2y} + C}{y^2 e^y}$$
化简后得到通解:
$$x = \frac{ey(2y2 - 2y + 1)}{4y^2} + \frac{C}{y^2 e^y}$$
你之前尝试的(u=yx)(即(x = \frac{u}{y}))代入原方程后,其实会得到另一个一阶线性方程:
$$\frac{du}{dy} + u\left(1 + \frac{1}{y}\right) = y e^y$$
这个方程也是可解的,但整理过程相对繁琐,容易绕弯路。而直接将(x)视为(y)的函数,能更快识别出线性方程的形式,解法更直接。
另外,也可以先检查原方程是否为全微分方程:计算(\frac{\partial M}{\partial y})((M=y))和(\frac{\partial N}{\partial x})((N=xy+2x-ye^y)),得到(\frac{\partial M}{\partial y}=1),(\frac{\partial N}{\partial x}=y+2),两者不相等,所以不是全微分方程,无需往这个方向尝试。
内容的提问来源于stack exchange,提问作者Rabih Assaf

