7-微分方程

2225 字
11 分钟
7-微分方程

各种结构的微分方程通过各种方式来分离变量

总览#

graph LR a["可分离变量"] b["有理分式"] b --->|"换元"| a

可分离变量#

f(y)dy=g(x)dxf(y) dy = g(x) dx

有理分式#

{y=φ(u),u=ay+f1(x)by+f2(x)\left\{ \begin{aligned} & y' = \varphi(u), \\ & u = \frac{ay + f_1(x)}{by + f_2(x)} \end{aligned} \right.

当且仅当f1(x)f_1(x)f2(x)f_2(x)是线性函数时有

y=f1f2uabuduF(u)=dxaf2bf1F(u)=(abu)2(f1f2uabu+φ(u))\begin{aligned} & y = -\frac{f_1 - f_2u}{a-bu} \\ & \frac{du}{F(u)} = \frac{dx}{af_2 - bf_1} \\ & F(u) = (a-bu)^2(\frac{f_1'-f_2'u}{a-bu} + \varphi(u)) \end{aligned}

一阶线性微分方程#

变系数齐次#

y+P(x)y=0y=±eP(x)dx=CeP(x)dx\begin{aligned} & y' + P(x)y = 0 \\ & y = \pm e^{\int -P(x) dx} = C \cdot e^{\int{-P(x)}dx} \end{aligned}

变系数非齐次#

y+P(x)y=Q(x)y' + P(x)y = Q(x)

使用常数变易法,设y=ueP(x)dxy = u e^{\int{-P(x)}dx}

y=ueP(x)dxP(x)ueP(x)dx=ueP(x)dxP(x)yy+P(x)y=ueP(x)dx=Q(x)du=Q(x)eP(x)dxdxy=Q(x)eP(x)dxdxeP(x)dx\begin{aligned} & y' = u'e^{\int{-P(x)}dx} - P(x) ue^{\int{-P(x)}dx} = u'e^{\int{-P(x)}dx} - P(x)y \\ & y' + P(x)y = u'e^{\int{-P(x)}dx} = Q(x) \\ & du = Q(x)e^{\int{P(x)}dx} dx \\ & y = \int{Q(x)e^{\int{P(x)}dx} dx}\cdot e^{\int{-P(x)}dx} \end{aligned}

高阶线性微分方程#

二阶线性微分方程#

变系数齐次#

y+P(x)y+Q(x)y=0y'' + P(x)y' + Q(x)y = 0

1. 已知一特解(降阶法)#

已知一解 y1y_1,设:

y2=u(x)y1y_2 = u(x) y_1

计算:

y2=uy1+uy1y2=uy1+2uy1+uy1\begin{aligned} y_2' &= u'y_1 + uy_1' \\ y_2'' &= u''y_1 + 2u'y_1' + uy_1'' \end{aligned}

代入原方程,并利用 y1y_1 满足原方程消去项,得到:

uy1+u(2y1+P(x)y1)=0u''y_1 + u'(2y_1' + P(x)y_1) = 0

令:

v=uv = u'

得到一阶方程:

v+(2y1y1+P(x))v=0v' + \left(2\frac{y_1'}{y_1} + P(x)\right)v = 0

解得:

v=Cy12eP(x)dxv = \frac{C}{y_1^2} e^{-\int P(x)dx}

再积分:

u=eP(x)dxy12dxu = \int \frac{e^{-\int P(x)dx}}{y_1^2} dx

最终:

y2=y1eP(x)dxy12dxy_2 = y_1 \int \frac{e^{-\int P(x)dx}}{y_1^2} dx

通解#

y=C1y1+C2y2y = C_1 y_1 + C_2 y_2

变系数非齐次#

y+P(x)y+Q(x)y=R(x)y'' + P(x)y' + Q(x)y = R(x)

1. 已知一特解#

若已知一个特解 y1y_1,设:

y=uy1y = u y_1

代入可得关于 uu 的一阶非齐次微分方程:

降阶后再求解\Rightarrow \text{降阶后再求解}

2. 已知齐次解#

已知齐次解:

y1,  y2y_1, \; y_2

设:

y=C1(x)y1+C2(x)y2y = C_1(x)y_1 + C_2(x)y_2

附加条件:

{C1y1+C2y2=0C1y1+C2y2=R(x)\begin{cases} C_1' y_1 + C_2' y_2 = 0 \\ C_1' y_1' + C_2' y_2' = R(x) \end{cases}

设朗斯基行列式:

W=y1y2y1y20W = \begin{vmatrix} y_1 & y_2 \\ y_1' & y_2' \end{vmatrix} \neq 0

解得:

C1=y2R(x)W,C2=y1R(x)WC_1' = -\frac{y_2 R(x)}{W}, \quad C_2' = \frac{y_1 R(x)}{W}

积分:

C1=y2R(x)Wdx,C2=y1R(x)WdxC_1 = \int -\frac{y_2 R(x)}{W} dx, \quad C_2 = \int \frac{y_1 R(x)}{W} dx

通解#

y=C1y1+C2y2+ypy = C_1 y_1 + C_2 y_2 + y_p

欧拉方程(Cauchy–Euler)#


1. 标准形式#

x2y+axy+by=0x^2 y'' + a x y' + b y = 0

2. 基本解法#

设:

x=etx = e^ty=xry = x^r

代入得到特征方程:

r(r1)+ar+b=0r(r-1)+ar+b=0

3. 解的分类#

(1) 不同实根 r1r2r_1 \ne r_2#

y=C1xr1+C2xr2y = C_1 x^{r_1} + C_2 x^{r_2}

(2) 重根 rr#

y=(C1+C2lnx)xry = (C_1 + C_2 \ln x)x^r

(3) 共轭复根 r=α±βir = \alpha \pm \beta i#

y=xα[C1cos(βlnx)+C2sin(βlnx)]y = x^\alpha \left[ C_1 \cos(\beta \ln x) + C_2 \sin(\beta \ln x) \right]

4. 非齐次欧拉方程#

x2y+axy+by=f(x)x^2 y'' + a x y' + b y = f(x)先解齐次解,再用常数变易法求特解\Rightarrow \text{先解齐次解,再用常数变易法求特解}

5. 本质#

x=et转化为常系数线性微分方程x = e^t \Rightarrow \text{转化为常系数线性微分方程}

其他补充#


1. 降阶公式#

y2=y1eP(x)dxy12dxy_2 = y_1 \int \frac{e^{-\int P(x)dx}}{y_1^2} dx

2. 朗斯基行列式#

W(y1,y2)=y1y2y1y2W(y_1,y_2) = \begin{vmatrix} y_1 & y_2 \\ y_1' & y_2' \end{vmatrix}W0线性无关W \ne 0 \Rightarrow \text{线性无关}

3. 方法选择#

已知一个解降阶法非齐次 + 已知两个解常数变易法欧拉型y=xr\begin{aligned} & \text{已知一个解} \Rightarrow \text{降阶法} \\ & \text{非齐次 + 已知两个解} \Rightarrow \text{常数变易法} \\ & \text{欧拉型} \Rightarrow y=x^r \\ \end{aligned}

伯努利方程#

y+P(x)y=yn11ndy1ndx+P(x)y1n=Q(x)换元一阶线性微分方程\begin{aligned} & y' + P(x)y = y^n \\ & \frac{1}{1-n} \cdot \frac{dy^{1-n}}{dx} + P(x)y^{1-n} = Q(x) \overset{\text{换元}}\Rightarrow \text{一阶线性微分方程} \end{aligned}

高阶微分方程#

n阶常系数齐次线性方程#

1. 标准形式#

any(n)+an1y(n1)++a1y+a0y=0a_n y^{(n)} + a_{n-1} y^{(n-1)} + \cdots + a_1 y' + a_0 y = 0

2. 特征方程#

设:

y=erxy = e^{rx}

代入得特征方程:

anrn+an1rn1++a1r+a0=0a_n r^n + a_{n-1} r^{n-1} + \cdots + a_1 r + a_0 = 0

3. 解的结构#

设特征方程根为 rr


(1) 互异实根 r1,r2,,rkr_1, r_2, \dots, r_k#

y=C1er1x+C2er2x++Ckerkxy = C_1 e^{r_1 x} + C_2 e^{r_2 x} + \cdots + C_k e^{r_k x}

(2) 重根 rr(重数为 mm#

y=(C1+C2x++Cmxm1)erxy = (C_1 + C_2 x + \cdots + C_m x^{m-1}) e^{rx}

(3) 共轭复根 r=α±βir = \alpha \pm \beta i(重数为1)#

y=eαx(C1cosβx+C2sinβx)y = e^{\alpha x} \left( C_1 \cos \beta x + C_2 \sin \beta x \right)

(4) 复根重数为 mm#

y=eαx[(C1+C2x++Cmxm1)cosβx+(D1+D2x++Dmxm1)sinβx]y = e^{\alpha x} \left[ (C_1 + C_2 x + \cdots + C_m x^{m-1}) \cos \beta x + (D_1 + D_2 x + \cdots + D_m x^{m-1}) \sin \beta x \right]

4. 通解#

y=(所有线性无关解)y = \sum \text{(所有线性无关解)}

n阶常系数非齐次线性方程#

1. 标准形式#

any(n)+an1y(n1)++a0y=f(x)a_n y^{(n)} + a_{n-1} y^{(n-1)} + \cdots + a_0 y = f(x)

2. 通解结构#

y=yh+ypy = y_h + y_p

其中:

  • yhy_h:对应齐次方程通解
  • ypy_p:任一特解

3. 特解方法一:待定系数法(核心)#

传统的#

适用:f(x)f(x) 为以下类型:

  • 多项式
  • 指数函数 eλxe^{\lambda x}
  • 三角函数 sinωx, cosωx\sin \omega x,\ \cos \omega x
  • 它们的乘积和线性组合

(1) 基本假设形式#

若:

f(x)=eλxPn(x)f(x) = e^{\lambda x} P_n(x)

设:

yp=xseλxQn(x)y_p = x^s e^{\lambda x} Q_n(x)

其中:

  • Qn(x)Q_n(x) 为同阶多项式
  • ss 为“调整次数”

(2) 调整规则(关键)#

λ\lambda 是特征方程的根,且重数为 kk

s=k\Rightarrow s = k

否则:

s=0\Rightarrow s = 0

(3) 三角函数情况#

若:

f(x)=eαx(Acosβx+Bsinβx)f(x) = e^{\alpha x} (A \cos \beta x + B \sin \beta x)

设:

yp=xseαx(Ccosβx+Dsinβx)y_p = x^s e^{\alpha x} (C \cos \beta x + D \sin \beta x)

结合算子的#

将待定系数法用微分算子统一表达,可以把“猜形式 + 调整次数”转化为一个结构化计算过程。


(1) 指数平移(核心变换)#

设:

L(D)y=f(x),f(x)=eλxPn(x)L(D)y = f(x), \quad f(x)=e^{\lambda x}P_n(x)

令:

yp=eλxR(x)y_p = e^{\lambda x} R(x)

则有算子恒等式:

L(D)(eλxR(x))=eλxL(D+λ)R(x)L(D)\big(e^{\lambda x}R(x)\big) = e^{\lambda x} L(D+\lambda)R(x)

因此原方程等价于:

L(D+λ)R(x)=Pn(x)L(D+\lambda)R(x) = P_n(x)

(2) 点积形式(算子展开)#

将算子展开:

L(D+λ)=j=0nL(j)(λ)j!DjL(D+\lambda) = \sum_{j=0}^n \frac{L^{(j)}(\lambda)}{j!} D^j

于是得到:

j=0nL(j)(λ)j!R(j)(x)=Pn(x)\sum_{j=0}^n \frac{L^{(j)}(\lambda)}{j!} R^{(j)}(x) = P_n(x)

可写成“向量点积形式”:

(R, R, , R(n))(L(λ), L(λ), L(λ)2!, )=Pn(x)(R,\ R',\ \dots,\ R^{(n)}) \cdot \left( L(\lambda),\ L'(\lambda),\ \frac{L''(\lambda)}{2!},\ \dots \right) = P_n(x)

这一步将微分算子的作用转化为有限维代数计算。


(3) 调整次数的算子解释#

传统规则:

λ 是特征根,重数 ks=k\lambda \text{ 是特征根,重数 } k \Rightarrow s=k

在算子中表现为:

若:

L(λ)=L(λ)==L(k1)(λ)=0L(\lambda)=L'(\lambda)=\cdots=L^{(k-1)}(\lambda)=0

则算子退化为:

L(D+λ)R(x)=L(k)(λ)k!R(k)(x)+L(D+\lambda)R(x) = \frac{L^{(k)}(\lambda)}{k!} R^{(k)}(x) + \cdots

即低阶导数项全部消失。


这意味着:

  • 方程阶数下降
  • 原多项式空间不再足够
  • 必须提高 R(x)R(x) 的次数

等价于传统中的:

yp=xkeλxQn(x)y_p = x^k e^{\lambda x} Q_n(x)

(4) 三角函数的统一处理#

若:

f(x)=eαx(Acosβx+Bsinβx)f(x)=e^{\alpha x}(A\cos\beta x + B\sin\beta x)

设:

λ=α+iβ\lambda = \alpha + i\beta

转化为:

f(x)=(eλxP(x))f(x)=\Re\big(e^{\lambda x}P(x)\big)

按指数情形计算:

L(D+λ)R(x)=P(x)L(D+\lambda)R(x)=P(x)

求得复解后取实部。


(5) 方法本质总结#

待定系数法可以统一理解为:

  1. 用指数平移:

    DD+λD \to D+\lambda
  2. 将微分算子展开为导数的线性组合

  3. 在多项式空间中解方程

  4. 若算子在 λ\lambda 处退化(零点),则通过提高次数(乘 xkx^k)补偿


(6) 一句话理解#

待定系数法本质上是:

在指数因子 eλxe^{\lambda x} 下,将微分算子 L(D)L(D) 转化为 L(D+λ)L(D+\lambda)
并在多项式空间中通过有限维代数(点积形式)求解特解。

4. 特解方法二:常数变易法(一般理论)#

设齐次解为:

y1,y2,,yny_1, y_2, \dots, y_n

设:

yp=i=1nCi(x)yiy_p = \sum_{i=1}^n C_i(x) y_i

满足条件:

{Ciyi=0Ciyi=0Ciyi(n2)=0Ciyi(n1)=f(x)an\begin{cases} \sum C_i' y_i = 0 \\ \sum C_i' y_i' = 0 \\ \cdots \\ \sum C_i' y_i^{(n-2)} = 0 \\ \sum C_i' y_i^{(n-1)} = \frac{f(x)}{a_n} \end{cases}

解线性方程组:#

设朗斯基行列式:

W=y1y2yny1y2yny1(n1)y2(n1)yn(n1)0W = \begin{vmatrix} y_1 & y_2 & \cdots & y_n \\ y_1' & y_2' & \cdots & y_n' \\ \vdots & \vdots & & \vdots \\ y_1^{(n-1)} & y_2^{(n-1)} & \cdots & y_n^{(n-1)} \end{vmatrix} \neq 0可解出 Ci(x)\Rightarrow \text{可解出 } C_i'(x)

再积分得到 Ci(x)C_i(x)


5. 解题结构总结#

齐次方程特征方程非齐次方程y=yh+ypf(x) 规则型待定系数法f(x) 一般型常数变易法\begin{aligned} & \text{齐次方程} \Rightarrow \text{特征方程} \\ & \text{非齐次方程} \Rightarrow y = y_h + y_p \\ & f(x)\ \text{规则型} \Rightarrow \text{待定系数法} \\ & f(x)\ \text{一般型} \Rightarrow \text{常数变易法} \end{aligned}

微分算子#

来源于该文章的总结

从线性变换的角度详解求微分方程特解的微分算子法 - 知乎


一、基本思想#

设常系数线性微分方程:

L(D)y=f(x)L(D)y = f(x)

其中:

L(D)=anDn++a1D+a0L(D) = a_n D^n + \cdots + a_1 D + a_0

核心思路是把 DD 当作“代数变量”,类比多项式运算。


二、指数型右端 ekxe^{kx}#

情况 1:非共振(L(k)0L(k) \neq 0#

yp=1L(k)ekxy_p = \frac{1}{L(k)} e^{kx}

情况 2:一次共振(L(k)=0, L(k)0L(k)=0,\ L'(k)\neq 0#

yp=1L(k)xekxy_p = \frac{1}{L'(k)} x e^{kx}

情况 3:高重根#

若:

L(k)=L(k)==L(m1)(k)=0,L(m)(k)0L(k)=L'(k)=\cdots=L^{(m-1)}(k)=0,\quad L^{(m)}(k)\neq 0

则:

yp=1L(m)(k)xmekxy_p = \frac{1}{L^{(m)}(k)} x^m e^{kx}

本质:乘 (x^m) 对应补齐算子零空间(广义特征向量 / Jordan 链)。


三、三角函数((\sin\omega x,\ \cos\omega x))#

利用欧拉公式:

sinωx,cosωxe±iωx\sin\omega x,\cos\omega x \longrightarrow e^{\pm i\omega x}

步骤:

  1. 转为复指数形式
  2. 按指数规则计算
  3. 取实部或虚部

四、多项式右端 (x^n)#

设:

yp=1L(D)xny_p = \frac{1}{L(D)} x^n

通过微分算子“除法”(类似多项式除法)求解,即在多项式空间中求算子逆作用。


五、指数 × 多项式 xnekxx^n e^{kx}#

利用指数平移:

DD+kD \to D + k

将方程转化为:

L(D+k)R(x)=xnL(D+k)R(x) = x^n

求得 R(x)\后:

yp=ekxR(x)y_p = e^{kx} R(x)

六、统一结构#

1. 算子代数化#

L(D)L(r)L(D) \leftrightarrow L(r)

2. 指数平移#

L(D)(ekxR(x))=ekxL(D+k)R(x)L(D)\big(e^{kx}R(x)\big) = e^{kx}L(D+k)R(x)

3. 共振处理#

L(k)=0乘 xL(k)=0 \Rightarrow \text{乘 } x

对应:

  • 特征根重数
  • 算子不可逆
  • Jordan 结构

七、本质理解#

微分算子法本质上是将微分方程转化为算子多项式在特定函数空间上的代数运算,并通过指数平移将问题降到多项式空间求解。


八、与待定系数法的关系#

  • 待定系数法:直接设解的形式并代入求系数
  • 微分算子法:通过算子结构推导解的形式

两者本质一致,但微分算子法更结构化,适合推广与统一理解。

文章分享

如果这篇文章对你有帮助,欢迎分享给更多人!

7-微分方程
https://skaco2.com/posts/02-math/7-微分方程/
作者
SKACO2
发布于
2026-05-03
许可协议
CC BY-NC-SA 4.0

评论区

Profile Image of the Author
SKACO2
Hello……
公告
欢迎来到我的博客!
音乐
封面

音乐

暂未播放

0:00 0:00
暂无歌词
分类
标签
站点统计
文章
53
分类
8
标签
54
总字数
58,255
运行时长
0
最后活动
0 天前

目录