高等数学 · 第八章

常微分方程

基本概念、一阶方程、高阶线性方程与应用

微分方程的基本概念

含未知函数及其导数的方程称为微分方程;出现的最高阶导数的阶数是方程的阶。使方程成立的函数称为解。通解通常含与方程阶数相同个数的独立任意常数,初始条件用于确定特解。解曲线称为积分曲线。

易错 “含 nn 个常数”只是通解的形式特征,还要保证这些常数相互独立。

一阶微分方程

可分离变量

dydx=f(x)g(y)  ⟹  ∫dyg(y)=∫f(x) dx.\frac{\mathrm{d}y}{\mathrm{d}x}=f(x)g(y) \implies \int\frac{\mathrm{d}y}{g(y)}=\int f(x)\,\mathrm{d}x.

条件 g(y)≠0g(y)\neq 0 的区间。使 g(y)=0g(y)=0 的常数解要单独检验。

易错 分离时不要漏掉使分母为零的平衡解。

齐次方程与可化为齐次的方程

形如 dydx=f(yx)\frac{\mathrm{d}y}{\mathrm{d}x}=f\bigl(\frac{y}{x}\bigr)。令 u=y/xu=y/x,则 y=uxy=ux,y′=u+xu′y'=u+xu',化为可分离。

变形 y′=f(ax+by+ca1x+b1y+c1)y'=f\bigl(\frac{ax+by+c}{a_1x+b_1y+c_1}\bigr) 可先平移,化为齐次。

一阶线性方程

y′+p(x)y=q(x).y'+p(x)y=q(x).

通解:

y=e−∫p dx(∫q e∫p dx dx+C).y=\mathrm{e}^{-\int p\,\mathrm{d}x}\left(\int q\,\mathrm{e}^{\int p\,\mathrm{d}x}\,\mathrm{d}x+C\right).

齐次式 q=0q=0 时,y=Ce−∫p dxy=C\mathrm{e}^{-\int p\,\mathrm{d}x}。

用途 积分因子法。标准形必须是 y′+py=qy'+p y=q,先除以 y′y' 的系数。

易错 写成 x′=P(y)x+Q(y)x'=P(y)x+Q(y) 时,应把 yy 当自变量。

伯努利方程

y′+p(x)y=q(x)yn(n≠0,1).y'+p(x)y=q(x)y^n\qquad (n\neq 0,1).

令 z=y1−nz=y^{1-n},化为线性方程。

全微分方程

M(x,y) dx+N(x,y) dy=0M(x,y)\,\mathrm{d}x+N(x,y)\,\mathrm{d}y=0。若 My=NxM_y=N_x,则存在 u(x,y)=Cu(x,y)=C 使 du=Mdx+Ndy\mathrm{d}u=M\mathrm{d}x+N\mathrm{d}y。

u=∫M dx+∫(N−∂∂y∫M dx)dy.u=\int M\,\mathrm{d}x+\int\left(N-\frac{\partial}{\partial y}\int M\,\mathrm{d}x\right)\mathrm{d}y.

若 My≠NxM_y\neq N_x,可尝试只依赖于 xx 或 yy 的积分因子:My−NxN\frac{M_y-N_x}{N} 只含 xx,则 μ=e∫My−NxNdx\mu=\mathrm{e}^{\int\frac{M_y-N_x}{N}\mathrm{d}x}。

高阶微分方程

可降阶方程

不显含 yy:F(x,y′,y′′)=0F(x,y',y'')=0,令 p=y′p=y',y′′=p′y''=p'。

不显含 xx:F(y,y′,y′′)=0F(y,y',y'')=0,令 p=y′p=y',y′′=pdpdyy''=p\dfrac{\mathrm{d}p}{\mathrm{d}y}。

易错 第二种换元后自变量是 yy,积分回来不要漏 yy 与 xx 的关系。

线性方程解的结构

y(n)+pn−1(x)y(n−1)+⋯+p0(x)y=f(x).y^{(n)}+p_{n-1}(x)y^{(n-1)}+\cdots+p_0(x)y=f(x).

齐次通解是 nn 个线性无关特解的线性组合。非齐次通解 == 齐次通解 ++ 一个特解。

叠加:右端 f1+f2f_1+f_2 时特解可相加。

二阶常系数齐次方程

y′′+py′+qy=0.y''+py'+qy=0.

特征方程 r2+pr+q=0r^2+pr+q=0:

根通解
两不等实根 r1,r2r_1,r_2C1er1x+C2er2xC_1\mathrm{e}^{r_1 x}+C_2\mathrm{e}^{r_2 x}
重根 rr(C1+C2x)erx(C_1+C_2 x)\mathrm{e}^{rx}
共轭复根 α+iβ\alpha+\mathrm{i}\beta、α−iβ\alpha-\mathrm{i}\betaeαx(C1cos⁡βx+C2sin⁡βx)\mathrm{e}^{\alpha x}(C_1\cos\beta x+C_2\sin\beta x)

二阶常系数非齐次方程

y′′+py′+qy=f(x).y''+py'+qy=f(x).

待定系数(ff 为多项式、指数、正余弦或其积):

  • f=Pm(x)eλxf=P_m(x)\mathrm{e}^{\lambda x}:若 λ\lambda 不是特征根,特解 Pm(x)eλxP_m(x)\mathrm{e}^{\lambda x};是单根,乘 xx;是重根,乘 x2x^2。
  • f=eαx(Acos⁡βx+Bsin⁡βx)f=\mathrm{e}^{\alpha x}(A\cos\beta x+B\sin\beta x):α+iβ\alpha+\mathrm{i}\beta 不是特征根则同型;是特征根则再乘 xx。

易错 特征根重合必须乘 xx 的相应次幂,否则求不出系数。sin⁡\sin 与 cos⁡\cos 要同时保留,即使右端只有一项。

简单应用

建立模型时先确定自变量、未知函数和变化率关系,再代入初始条件。常见模型包括指数增长与衰减、冷却、混合和振动。求出通解后必须检查单位、初值和解的实际定义域。

搜索

可搜索:洛必达、格林公式、正态分布、等价无穷小…