各种结构的微分方程通过各种方式来分离变量
graph LR
a["可分离变量"]
b["有理分式"]
b --->|"换元"| a
可分离变量#
f(y)dy=g(x)dx
有理分式#
⎩⎨⎧y′=φ(u),u=by+f2(x)ay+f1(x)当且仅当f1(x)和f2(x)是线性函数时有
y=−a−buf1−f2uF(u)du=af2−bf1dxF(u)=(a−bu)2(a−buf1′−f2′u+φ(u))
一阶线性微分方程#
变系数齐次#
y′+P(x)y=0y=±e∫−P(x)dx=C⋅e∫−P(x)dx变系数非齐次#
y′+P(x)y=Q(x)使用常数变易法,设y=ue∫−P(x)dx
y′=u′e∫−P(x)dx−P(x)ue∫−P(x)dx=u′e∫−P(x)dx−P(x)yy′+P(x)y=u′e∫−P(x)dx=Q(x)du=Q(x)e∫P(x)dxdxy=∫Q(x)e∫P(x)dxdx⋅e∫−P(x)dx
高阶线性微分方程#
二阶线性微分方程#
变系数齐次#
y′′+P(x)y′+Q(x)y=0
1. 已知一特解(降阶法)#
已知一解 y1,设:
y2=u(x)y1计算:
y2′y2′′=u′y1+uy1′=u′′y1+2u′y1′+uy1′′代入原方程,并利用 y1 满足原方程消去项,得到:
u′′y1+u′(2y1′+P(x)y1)=0令:
v=u′得到一阶方程:
v′+(2y1y1′+P(x))v=0解得:
v=y12Ce−∫P(x)dx再积分:
u=∫y12e−∫P(x)dxdx最终:
y2=y1∫y12e−∫P(x)dxdx
y=C1y1+C2y2
变系数非齐次#
y′′+P(x)y′+Q(x)y=R(x)
1. 已知一特解#
若已知一个特解 y1,设:
y=uy1代入可得关于 u 的一阶非齐次微分方程:
⇒降阶后再求解
2. 已知齐次解#
已知齐次解:
y1,y2设:
y=C1(x)y1+C2(x)y2附加条件:
{C1′y1+C2′y2=0C1′y1′+C2′y2′=R(x)设朗斯基行列式:
W=y1y1′y2y2′=0解得:
C1′=−Wy2R(x),C2′=Wy1R(x)积分:
C1=∫−Wy2R(x)dx,C2=∫Wy1R(x)dx
y=C1y1+C2y2+yp
欧拉方程(Cauchy–Euler)#
1. 标准形式#
x2y′′+axy′+by=0
2. 基本解法#
设:
x=ety=xr代入得到特征方程:
r(r−1)+ar+b=0
3. 解的分类#
(1) 不同实根 r1=r2#
y=C1xr1+C2xr2
(2) 重根 r#
y=(C1+C2lnx)xr
(3) 共轭复根 r=α±βi#
y=xα[C1cos(βlnx)+C2sin(βlnx)]
4. 非齐次欧拉方程#
x2y′′+axy′+by=f(x)⇒先解齐次解,再用常数变易法求特解
5. 本质#
x=et⇒转化为常系数线性微分方程
其他补充#
1. 降阶公式#
y2=y1∫y12e−∫P(x)dxdx
2. 朗斯基行列式#
W(y1,y2)=y1y1′y2y2′W=0⇒线性无关
3. 方法选择#
已知一个解⇒降阶法非齐次 + 已知两个解⇒常数变易法欧拉型⇒y=xr
伯努利方程#
y′+P(x)y=yn1−n1⋅dxdy1−n+P(x)y1−n=Q(x)⇒换元一阶线性微分方程
高阶微分方程#
n阶常系数齐次线性方程#
1. 标准形式#
any(n)+an−1y(n−1)+⋯+a1y′+a0y=0
2. 特征方程#
设:
y=erx代入得特征方程:
anrn+an−1rn−1+⋯+a1r+a0=0
3. 解的结构#
设特征方程根为 r
(1) 互异实根 r1,r2,…,rk#
y=C1er1x+C2er2x+⋯+Ckerkx
(2) 重根 r(重数为 m)#
y=(C1+C2x+⋯+Cmxm−1)erx
(3) 共轭复根 r=α±βi(重数为1)#
y=eαx(C1cosβx+C2sinβx)
(4) 复根重数为 m#
y=eαx[(C1+C2x+⋯+Cmxm−1)cosβx+(D1+D2x+⋯+Dmxm−1)sinβx]
4. 通解#
y=∑(所有线性无关解)
n阶常系数非齐次线性方程#
1. 标准形式#
any(n)+an−1y(n−1)+⋯+a0y=f(x)
2. 通解结构#
y=yh+yp其中:
- yh:对应齐次方程通解
- yp:任一特解
3. 特解方法一:待定系数法(核心)#
传统的#
适用:f(x) 为以下类型:
- 多项式
- 指数函数 eλx
- 三角函数 sinωx, cosωx
- 它们的乘积和线性组合
(1) 基本假设形式#
若:
f(x)=eλxPn(x)设:
yp=xseλxQn(x)其中:
- Qn(x) 为同阶多项式
- s 为“调整次数”
(2) 调整规则(关键)#
若 λ 是特征方程的根,且重数为 k:
⇒s=k否则:
⇒s=0
(3) 三角函数情况#
若:
f(x)=eαx(Acosβx+Bsinβx)设:
yp=xseαx(Ccosβx+Dsinβx)
结合算子的#
将待定系数法用微分算子统一表达,可以把“猜形式 + 调整次数”转化为一个结构化计算过程。
(1) 指数平移(核心变换)#
设:
L(D)y=f(x),f(x)=eλxPn(x)令:
yp=eλxR(x)则有算子恒等式:
L(D)(eλxR(x))=eλxL(D+λ)R(x)因此原方程等价于:
L(D+λ)R(x)=Pn(x)
(2) 点积形式(算子展开)#
将算子展开:
L(D+λ)=j=0∑nj!L(j)(λ)Dj于是得到:
j=0∑nj!L(j)(λ)R(j)(x)=Pn(x)可写成“向量点积形式”:
(R, R′, …, R(n))⋅(L(λ), L′(λ), 2!L′′(λ), …)=Pn(x)这一步将微分算子的作用转化为有限维代数计算。
(3) 调整次数的算子解释#
传统规则:
λ 是特征根,重数 k⇒s=k在算子中表现为:
若:
L(λ)=L′(λ)=⋯=L(k−1)(λ)=0则算子退化为:
L(D+λ)R(x)=k!L(k)(λ)R(k)(x)+⋯即低阶导数项全部消失。
这意味着:
- 方程阶数下降
- 原多项式空间不再足够
- 必须提高 R(x) 的次数
等价于传统中的:
yp=xkeλxQn(x)
(4) 三角函数的统一处理#
若:
f(x)=eαx(Acosβx+Bsinβx)设:
λ=α+iβ转化为:
f(x)=ℜ(eλxP(x))按指数情形计算:
L(D+λ)R(x)=P(x)求得复解后取实部。
(5) 方法本质总结#
待定系数法可以统一理解为:
-
用指数平移:
D→D+λ
-
将微分算子展开为导数的线性组合
-
在多项式空间中解方程
-
若算子在 λ 处退化(零点),则通过提高次数(乘 xk)补偿
(6) 一句话理解#
待定系数法本质上是:
在指数因子 eλx 下,将微分算子 L(D) 转化为 L(D+λ),
并在多项式空间中通过有限维代数(点积形式)求解特解。
4. 特解方法二:常数变易法(一般理论)#
设齐次解为:
y1,y2,…,yn设:
yp=i=1∑nCi(x)yi满足条件:
⎩⎨⎧∑Ci′yi=0∑Ci′yi′=0⋯∑Ci′yi(n−2)=0∑Ci′yi(n−1)=anf(x)
解线性方程组:#
设朗斯基行列式:
W=y1y1′⋮y1(n−1)y2y2′⋮y2(n−1)⋯⋯⋯ynyn′⋮yn(n−1)=0⇒可解出 Ci′(x)再积分得到 Ci(x)
5. 解题结构总结#
齐次方程⇒特征方程非齐次方程⇒y=yh+ypf(x) 规则型⇒待定系数法f(x) 一般型⇒常数变易法
微分算子#
来源于该文章的总结
从线性变换的角度详解求微分方程特解的微分算子法 - 知乎
一、基本思想#
设常系数线性微分方程:
L(D)y=f(x)其中:
L(D)=anDn+⋯+a1D+a0核心思路是把 D 当作“代数变量”,类比多项式运算。
二、指数型右端 ekx#
情况 1:非共振(L(k)=0)#
yp=L(k)1ekx
情况 2:一次共振(L(k)=0, L′(k)=0)#
yp=L′(k)1xekx
情况 3:高重根#
若:
L(k)=L′(k)=⋯=L(m−1)(k)=0,L(m)(k)=0则:
yp=L(m)(k)1xmekx本质:乘 (x^m) 对应补齐算子零空间(广义特征向量 / Jordan 链)。
三、三角函数((\sin\omega x,\ \cos\omega x))#
利用欧拉公式:
sinωx,cosωx⟶e±iωx步骤:
- 转为复指数形式
- 按指数规则计算
- 取实部或虚部
四、多项式右端 (x^n)#
设:
yp=L(D)1xn通过微分算子“除法”(类似多项式除法)求解,即在多项式空间中求算子逆作用。
五、指数 × 多项式 xnekx#
利用指数平移:
D→D+k将方程转化为:
L(D+k)R(x)=xn求得 R(x)\后:
yp=ekxR(x)
六、统一结构#
1. 算子代数化#
L(D)↔L(r)
2. 指数平移#
L(D)(ekxR(x))=ekxL(D+k)R(x)
3. 共振处理#
L(k)=0⇒乘 x对应:
七、本质理解#
微分算子法本质上是将微分方程转化为算子多项式在特定函数空间上的代数运算,并通过指数平移将问题降到多项式空间求解。
八、与待定系数法的关系#
- 待定系数法:直接设解的形式并代入求系数
- 微分算子法:通过算子结构推导解的形式
两者本质一致,但微分算子法更结构化,适合推广与统一理解。