mobile wallpaper 1mobile wallpaper 2mobile wallpaper 3
2327 字
6 分钟
微积分下 06 - 常系数线性方程

常系数微分方程

F(y)=i=1naiy(i)F(y)=\sum_{i=1}^{n}a_iy^{(i)}

有一些独特的性质

齐次#

我们考虑二阶齐次线性常系数微分方程

y+py+qy=0y''+py'+qy=0

有没有办法像处理多项式一样处理这个微分方程呢?为了做到这一点,我们需要将其变为类似于a2+pa+qa^2+pa+q的形式,于是我们自然想到了erxe^rx,因为它一旦求导,主干不会改变,只会多出一个常数rr,于是我们令y=erxy=e^{rx},上式变为

erx(r2+pr+q)=0e^{rx}(r^2+pr+q)=0

显然erx0e^{rx}\neq0,就有

r2+pr+q=0r^2+pr+q=0

这时只需算出r1,r2r_1,r_2,显然e1re2rC\frac{e^r_1}{e^r_2}\neq C,就可以得出通解为y=C1er1x+C2er2xy=C_1e^{r_1x}+C_2e^{r_2x}

不过存在有重根 r=r1=r2r=r_1=r_2 的情况, 这时候

y=(C1+C2x)erxy=(C_1+C_2x)e^{rx}

r1,r2r_1,r_2 也叫微分方程的特征根

非齐次#

那么,对于等式右侧不为0,而是Q(x)的情况呢?

这时候我们引入一个微分算子,表示对x求导这一操作本身

D=ddxD=\frac{d}{dx}

于是我们就得到一个特殊的多项式

(D2+pD+q)y=Q(x)y=1D2+pD+qQ(x)\begin {aligned} (D^2+pD+q)y=Q(x)\\ y=\frac{1}{D^2+pD+q}Q(x) \end {aligned}

我们令

P(D)=D2+pD+qP(D)=D^2+pD+q

这时候,根据Q(x)Q(x)的不同,会出现解法的差异

Q(x)=erxQ(x)=e^{rx}#

由于指数函数求导就是

D(erx)=rerxD(e^{rx})=re^{rx}

所以积分也就是

1Derx=1kerx\frac{1}{D}e^{rx}=\frac{1}{k}e^{rx}

因此

y=1P(D)erx=1P(k)erx=1k2+pk+qerxy=\frac{1}{P(D)}e^{rx}=\frac{1}{P(k)}e^{rx}=\frac{1}{k^2+pk+q}e^{rx}

也就是特解了

Q(x)=erxPm(x)Q(x)=e^{rx}P_m(x)#

其中 Pm(x)P_m(x) 为多项式

类似上面一个情况,我们猜测一个y的形式,令一个特解为

yp=erxQ(x)y_p=e^{rx}Q(x)

于是

yp=rerxQ(x)+erxQ(x)=erx(rQ(x)+Q(x))y_p'=re^{rx}Q(x)+e^{rx}Q'(x)=e^{rx}(rQ(x)+Q'(x))yp=erx(Q(x)+2rQ(x)+r2Q(x))y_p''=e^{rx}(Q''(x)+2rQ'(x)+r^2Q(x))

代回到原方程,约掉 erxe^{rx} 就有

Q+(2r+p)Q+(r2+pr+q)Q=Pm(x)Q''+(2r+p)Q'+(r^2+pr+q)Q=P_m(x)

就得到了新的方程

这时候我们把原方程的特征方程设为

C(λ)=λ2+pλ+q=0C(\lambda)=\lambda^2+p\lambda+q=0

于是就有

C(λ)=2λ+pC'(\lambda)=2\lambda+p

所以就可以把Q的方程写为

Q(x)+C(r)Q(x)+C(r)Q(x)=Pm(x)Q''(x)+C'(r)Q'(x)+C(r)Q(x)=P_m(x)

然后我们根据r和 C(λ)=0C(\lambda)=0 的关系分情况讨论

r不是 C(λ)=0C(\lambda)=0 的根#

这时候,方程左右两边相等,而右边是m次多项式,而左侧最高次数项是Q(x)

这时候 C(r)0C(r)\neq 0 ,因此Q(x)本身就是一个m次多项式,令 Q(x)=Rm(x)Q(x)=R_m(x)

这时候我们只需待定系数

Rm(x)=C0+C1x+...+CmxmR_m(x)=C_0+C_1x+...+C_mx^m

再代入 Q(x)+C(r)Q(x)+C(r)Q(x)=Pm(x)Q''(x)+C'(r)Q'(x)+C(r)Q(x)=P_m(x) 计算出系数,就可以得出y的特解为

yp=erxRm(x)y_p=e^{rx}R_m(x)

r是 C(λ)=0C(\lambda)=0 的根#

这时候 C(r)=0,C(r)0C(r)=0, C'(r)\neq 0 ,就有

Q(x)+C(r)Q(x)=Pm(x)Q''(x)+C'(r)Q'(x)=P_m(x)

这时候,要求 Q(x)Q'(x) 是一个m+1次常数项,而且由于常数项求导就会消失,我们直接假定没有常数项为0,再从 Q(x)Q'(x) 中提出一个x,就有

你可能想问,卧槽凭啥,你怎么想假定什么就假定什么。但是不要忘记我们是在求「特解」,只需是一个解即可,而非通解

Q(x)=xRm(x)Q(x)=xR_m(x)

同理代入 Q(x)+C(r)Q(x)+C(r)Q(x)=Pm(x)Q''(x)+C'(r)Q'(x)+C(r)Q(x)=P_m(x) 然后待定系数地求解即可有特解

yp=xerxRm(x)y_p=xe^{rx}R_m(x)

r是 C(λ)=0C(\lambda) = 0 的二重根#

这时候 C(λ)=(λr)2C(\lambda)=(\lambda-r)^2C(r)=C(r)=0C(r)=C'(r)=0 所以同上我们设

Q(x)=x2Rm(x)Q(x)=x^2R_m(x)

待定系数法求解即有

yp=x2Rm(x)y_p=x^2R_m(x)

Q(x)=sinkxcoskxQ(x)=\sin{kx}或\cos{kx}#

我们根据三角函数的特性知道,对正弦或余弦函数求导两次就会有

D2y(kx)=k2y(kx),y(x)=sinxcosxD^2y(kx)=-k^2y(kx), y(x)=\sin x或\cos x

于是我们就i可以直接用 k2-k^2 替换 D2D^2

y=1k2+pD+qsinkxy=\frac{1}{-k^2+pD+q}\sin{kx}

这时候我们乘以分母的共轭项,对分母升次来处理单独的 DD

y=pDq+k2(pD+qk2)(pDq+k2)sinkx=pDq+k2p2(k2)(qk2)2sinkx(p2(k2)(qk2)20)y=\frac{pD-q+k^2}{(pD+q-k^2)(pD-q+k^2)}\sin{kx}=\frac{pD-q+k^2}{p^2(-k^2) - (q-k^2)^2}\sin{kx} \quad (p^2(-k^2) - (q-k^2)^2\neq0)

这样一来,分母变为常数,就很好解决了

y=pkp2(k2)(qk2)2coskx+q+k2p2(k2)(qk2)2sinkxy=\frac{pk}{p^2(-k^2) - (q-k^2)^2}\cos{kx}+\frac{-q+k^2}{p^2(-k^2) - (q-k^2)^2}sin{kx}

我们就有了方程的特解

分母为零时,p=0,q=k2p=0,q=k^2y+k2y=sinkxy''+k^2y=\sin{kx},实际上就是共振方程

在这一情况下

yp=x2kcoskxy_p=-\frac{x}{2k}coskx

Q(x)=erx[Pl(x)cosωx+Qn(x)sinωx]Q(x)=e^{rx}[P_l(x)\cos{\omega x}+Q_n(x)\sin{\omega x}]#

这时候,我们不希望处理三角函数,因此直接使用欧拉公式来化解三角函数

欧拉公式: eiθ=cosθ+isinθe^{i\theta}=\cos\theta+i\sin\theta

就有

cosθ=12(eiθ+eiθ)\cos\theta=\frac{1}{2}(e^{i\theta}+e^{-i\theta})sinθ=12i(eiθeiθ)\sin\theta=\frac{1}{2i}(e^{i\theta}-e^{-i\theta})

这时候我们构造一个复数多项式

Hm(x)=Pl(x)iQn(x)H_m(x)=P_l(x)-iQ_n(x)

然后构造

Hm(x)e(r+iω)x=Hm(x)erx(cosωx+isinωx),其中m=max(l,n)H_m(x)e^{(r+i\omega)x}=H_m(x)e^{rx}(\cos{\omega x}+i\sin{\omega x}), 其中m=max(l,n)

展开就有

erx[Pl(x)cosωx+Qn(x)sinωx]+ierx[Pl(x)sinωxQn(x)cosωx]e^{rx}[P_l(x)\cos{\omega x +Q_n(x)\sin{\omega x}}]+ie^{rx}[P_l(x)\sin{\omega x}-Q_n(x)\cos{\omega x}]

不难发现,这个式子左侧就是 Q(x)Q(x)Q(x)Q(x) 就是一个复变函数的实部

于是我们构造一个复数函数 z(x)z(x) 为上述式子,它的特解的实部就是 ypy_p 了,再令 λ=r+iω\lambda = r+i\omega,它满足

z+pz+qz=Hm(x)eλxz''+pz'+qz=H_m(x)e^{\lambda x}

这时候,就得到了我们朝思暮想的多项式形式的方程,成功拿走了三角函数,那么我们运用之前的结论,就可以设特解

zp=xkRm(x)eλx,kλ作为特征方程解的重数,不是解则为0z_p=x^kR_m(x)e^{\lambda x},k为\lambda作为特征方程解的重数,不是解则为0

那么显然 Rm(x)R_m(x) 是复数多项式,于是我们拆开, Rm(x)=Am(x)+iBm(x)R_m(x)=A_m(x)+iB_m(x)eλx=erx(cosωx+isinωx)e^{\lambda x}=e^{rx}(\cos{\omega x}+i\sin{\omega x}) (欧拉公式),那么

zp=xk[Am(x)+iBm(x)]erx(cosωx+isinωx)z_p = x^k [A_m(x) + i B_m(x)] e^{rx} (\cos{\omega x} + i\sin{\omega x})

像刚刚那样,我们提取出实数部分,就有了特解形式的公式

yp=Re[z(x)]=xkerx[Am(x)cosωxBm(x)sinωx]y_p= \text{Re}[z^*(x)] = x^k e^{rx} [A_m(x)\cos{\omega x} - B_m(x)\sin{\omega x}]

同样使用待定系数的方法回代,求解出 Am(x)A_m(x)Bm(x)B_m(x) ,我们就得到了特解

通解#

你可能发现了刚刚我们只求了特解,那通解呢?现在让我们假定特解是 ypy_p ,它满足了

P(D)yp=Q(x)P(D)y_p=Q(x)

但是它不通用,于是我们想,可以加入一个等价于0的项,而且和原有方程相关,然后我们就发现

P(D)y=0P(D)y=0

的通解,也就是这一非齐次方程对应的其次方程的通解 yhy_h 就可以很好地充当0了!

好糟糕的表述

于是我们把这两组方程加在一起,就有了通解的形式

y=yh+ypy=y_h+y_p

给出特解求通解和反推微分方程#

对于常系数齐次线性微分方程,当我们需要反推的时候,就需要用到特征根了,特征根对应通解中的项的方式如下:

特征根通解中的对应项
实数rCerxCe^{rx}
一对共轭复数根 r=α±βir=\alpha \pm \beta i两项 eαx(C1cosβx+C2sinβxe^{\alpha x}(C_1\cos{\beta x}+C_2\sin{\beta x})
k重实根rk项 erx(C1+C2x+...+Ckxk1)e^{rx}(C_1+C_2x+...+C_k x^{k-1})
一对k重复数根 r=α±βir=\alpha \pm \beta i2k项 ekx[(C1+C2x+...+Ckxk1)cosβx+(D1+D2x+...+Dkxk1)sinβx]e^{kx}[(C_1+C_2x+...+C_k x^{k-1})\cos{\beta x} + (D_1+D_2x+...+D_k x^{k-1})\sin{\beta x}]

例题#

y1(x)=sinx+ex,y2(x)=cosx+ex,y3(x)=exy_1(x)=\sin{x}+e^x, y_2(x)=\cos{x}+e^x,y_3(x)=e^x为一个二阶线性微分方程的三个特解。求此微分方程的通解和微分方程的具体形式

对于线性非齐次微分方程,任意两个特解的,必定是对应齐次方程的解。

已知三个特解为: y1(x)=sinx+exy_1(x) = \sin x + e^x y2(x)=cosx+exy_2(x) = \cos x + e^x y3(x)=exy_3(x) = e^x

我们用它们作差,来寻找齐次方程的解: y1y3=(sinx+ex)ex=sinxy_1 - y_3 = (\sin x + e^x) - e^x = \sin x y2y3=(cosx+ex)ex=cosxy_2 - y_3 = (\cos x + e^x) - e^x = \cos x

因为 sinx\sin xcosx\cos x 线性无关(二者不成比例),且原方程是二阶的,所以我们已经找齐了齐次方程的基础解系,即 yh1=sinxy_{h1} = \sin x, yh2=cosxy_{h2} = \cos x

二阶线性非齐次微分方程的通解公式为: 通解 = 齐次方程的通解 + 非齐次方程的任意一个特解

所以,该方程的通解为: y(x)=C1sinx+C2cosx+exy(x) = C_1 \sin x + C_2 \cos x + e^x (注:特解选 y1,y2,y3y_1, y_2, y_3 中的哪一个都可以,通常选最简单的 y3=exy_3=e^xC1,C2C_1, C_2 为任意常数)

2. 求微分方程的具体形式#

既然齐次方程的基础解系是 sinx\sin xcosx\cos x,这对应着特征根为 r1,2=±ir_{1,2} = \pm i。 由特征根反推特征方程:r2+1=0r^2 + 1 = 0。 因此,对应的齐次方程左端形式为:y+y=0y'' + y = 0

这意味着原非齐次方程的形式一定为: y+y=f(x)y'' + y = f(x)

现在我们需要求出右端的非齐次项 f(x)f(x)。只需将任意一个特解(比如最简单的 y3=exy_3 = e^x)代入方程左边即可: f(x)=(ex)+(ex)=ex+ex=2exf(x) = (e^x)'' + (e^x) = e^x + e^x = 2e^x

所以,该微分方程的具体形式为: y+y=2exy'' + y = 2e^x

(你可以把 y1y_1y2y_2 代入检验,等式两边一定也是相等的)


第二部分:方法教学 —— 如何求解这类问题?#

这类问题通常的题设是:“已知某 nn 阶线性微分方程的 n+1n+1 个特解,求方程通解及原方程”。 掌握以下几个核心步骤,这类题就是送分题:

第一步:利用“特解之差”找齐次方程的基础解系#

核心原理:非齐次方程的解 - 非齐次方程的解 == 齐次方程的解。

  • 做法:通常选其中最简单的一个特解(记为 y0y_0),让其他的特解都去减它,得到 y1y0,y2y0y_1-y_0, y_2-y_0 \dots
  • 目的:对于二阶方程,你需要找到两个线性无关的齐次解;对于三阶,需要找到三个。

第二步:写出通解#

核心原理:解的结构定理。

  • 做法:把第一步得到的齐次解乘上任意常数 C1,C2C_1, C_2 拼起来,再加上原题目中给的任意一个特解。
  • 公式y=C1(y1y0)+C2(y2y0)+y0y = C_1(y_1-y_0) + C_2(y_2-y_0) + y_0

第三步:反推微分方程的形式#

这一步有两种情况,我教你两招:

情况A:基础解系很漂亮(常系数方程)—— 像本题一样

  1. 观察基础解系(如 eax,sin(bx),cos(bx),xeaxe^{ax}, \sin(bx), \cos(bx), x e^{ax} 等)。
  2. 写出对应的特征根 rr
    • 基础解系是 e2x,e3xe^{2x}, e^{3x} \Rightarrow 特征根是 r=2,r=3r=2, r=3 \Rightarrow (r2)(r3)=0r25r+6=0y5y+6y=0(r-2)(r-3)=0 \Rightarrow r^2-5r+6=0 \Rightarrow y''-5y'+6y=0
    • 基础解系是 sinx,cosx\sin x, \cos x \Rightarrow 特征根是 r=±ir=\pm i \Rightarrow r2+1=0y+y=0r^2+1=0 \Rightarrow y''+y=0
  3. 得到左边后,把已知的特解代入左边,算出来的结果就是方程右边的 f(x)f(x)

情况B:基础解系比较乱(变系数方程)—— 万能消元法 如果减出来的齐次解是 xxx2x^2 这种,不能用特征方程(因为不是常系数),怎么反推呢? 设通解为 y=C1yh1+C2yh2+yy = C_1 y_{h1} + C_2 y_{h2} + y^* 移项得到齐次部分:Y=yy=C1yh1+C2yh2Y = y - y^* = C_1 y_{h1} + C_2 y_{h2}

  1. YY 求一阶导数:Y=C1yh1+C2yh2Y' = C_1 y_{h1}' + C_2 y_{h2}'
  2. YY 求二阶导数:Y=C1yh1+C2yh2Y'' = C_1 y_{h1}'' + C_2 y_{h2}''
  3. C1,C2C_1, C_2 当作未知数,利用这三个式子消去 C1C_1C2C_2,找出一个只包含 Y,Y,YY, Y', Y'' 的等式。
  4. 最后把 Y=yyY = y - y^* 代回去化简,就能得到最终的微分方程。

(注:大部分考研或期末考试题,用情况A的“特征根反推法”就能秒杀,情况B作为进阶防身即可。)

特征方程式とは?特征方程即为微分方程 C1y(n)+...+Cny=Q(x)C_1y^{(n)}+...+C_ny=Q(x) 中,将每个y的i阶导数都替换为 rir^i 后,得到的方程,计算出来的值即为特征值

那么,特征根是共轭复根,即r=a±bir=a\pm bi的的情况呢?我们在上述的表格中已经说明过了,这时候通解中含有三角函数

n阶线性常系数微分方程拥有一样的解法,不再赘述了

分享

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

微积分下 06 - 常系数线性方程
https://hatoya-doublepigeonblog.pages.dev/posts/微积分下_06_-_常系数线性方程/
作者
两只鸽子
发布于
2026-03-15
许可协议
CC BY-NC-SA 4.0

部分信息可能已经过时