每日一题 #002
题目:
已知微分方程 y′′−2(1+tan2x)y=0 有一个解 y1=tanx,求其通解。
[!NOTE]
关键点:刘维尔公式 常数变易法
[!TIP]
本文中验证两特解线性无关的方法为朗斯基行列式(Wronskian Determinant),更加直观与明显。
在日常计算时,也可选用比例法,即检查是否存在一个常数 k,使得 y2(x)y1(x)=k:
- 如果比值是常数,则两个函数线性相关。
- 如果比值依赖于 x(不是常数),则两个函数线性无关。
示例:exe2x=ex 结果随 x 变化,因此两个函数线性无关。
PROOF1:(使用刘维尔公式)
步骤1:将方程化为标准形式
二阶线性齐次微分方程的通用标准形式为:
y′′+P(x)y′+Q(x)y=0
对比原方程 y′′−2(1+tan2x)y=0,无 y′ 项,因此对应系数为:
P(x)=0,Q(x)=−2(1+tan2x)
步骤2:刘维尔降阶公式(二阶方程专用)
对于二阶线性齐次方程,若已知一个非零特解 y1(x),则第二个与它线性无关的特解,可由刘维尔公式直接推导:
y2(x)=y1(x)⋅∫y12(x)e−∫P(x)dxdx
公式原理:基于刘维尔定理,朗斯基行列式满足 W(x)=W(x0)e−∫x0xP(t)dt,以此实现降阶,快速求解第二个特解。
步骤3:代入特解计算第二个线性无关解
已知特解 y1(x)=tanx,且 P(x)=0,分步计算:
-
计算指数项
∫P(x)dx=∫0dx=0,e−∫P(x)dx=e0=1
-
代入刘维尔降阶公式
y2(x)=tanx⋅∫tan2x1dx=tanx⋅∫cot2xdx
-
计算三角积分(利用恒等式 cot2x=csc2x−1)
∫cot2xdx=∫(csc2x−1)dx=−cotx−x+C
-
取积分常数 C=0(仅需一个特解,常数可吸收到通解的任意常数中),代入化简
y2(x)=tanx⋅(−cotx−x)=−1−xtanx
步骤4:验证线性无关性(朗斯基行列式)
两个函数线性无关的充要条件:其朗斯基行列式在定义域内不恒为0。
朗斯基行列式定义:
W(y1,y2)=y1y1′y2y2′
先求导:
y1′=sec2x=1+tan2x,y2′=−tanx−xsec2x
代入行列式展开计算:
W=tanx⋅(−tanx−xsec2x)−sec2x⋅(−1−xtanx)=−tan2x−xtanxsec2x+sec2x+xtanxsec2x=sec2x−tan2x=1=0
朗斯基行列式在定义域内恒不为0,因此 y1 与 y2 线性无关。
步骤5:写出方程的通解
二阶线性齐次微分方程的通解,为两个线性无关特解的线性组合:
y=C1tanx+C2(1+xtanx)
注:C1,C2 为任意常数,原 y2 的负号已吸收进任意常数 C2 中,不影响通解的一般性。
定义域补充说明
由于 tanx、sec2x 的定义域为 x=kπ+2π,k∈Z,因此通解在每个区间 (kπ−2π,kπ+2π) 内成立。
利用"刘维尔定理"解题的思路:
刘维尔公式的核心价值,是解决「二阶线性齐次微分方程,已知一个非零特解,求第二个线性无关特解,最终得到通解」的问题。
二阶线性齐次微分方程的通解结构定理,是所有操作的前提:
对于标准形式的二阶线性齐次方程:
y′′+P(x)y′+Q(x)y=0
若 y1(x) 和 y2(x) 是方程的两个线性无关的特解(即不成比例),则方程的通解为:
y=C1y1(x)+C2y2(x)(C1,C2为任意常数)
回到本题,题目已经给出一个特解 y1=tanx,因此我们的唯一核心任务,就是找到第二个和 y1 线性无关的特解 y2,而刘维尔公式,就是专门用来快速求这个 y2 的固化公式,不用每次都手动降阶推导。
刘维尔降阶公式:
y2(x)=y1(x)⋅∫y12(x)e−∫P(x)dxdx
一句话总结:刘维尔公式把「求二阶方程特解」的复杂问题,直接简化成了「求一个不定积分」的简单问题,跳过了手动设函数、求导、化简、降阶的全过程。
解题思路:
将通解中与 y1 搭配的常数,替换为关于 x 的未知函数 u(x),即设第二个线性无关特解为 y2=u(x)⋅y1(x)。
将其代入原方程后,利用 y1 是解的性质,可消去含 u(x) 的项,将二阶方程降为关于 u′(x) 的一阶微分方程,求解后即可得到 u(x),最终得到第二个特解与通解。
原方程:
y′′−2(1+tan2x)y=0
对比标准形式 y′′+P(x)y′+Q(x)y=0,得:
P(x)=0,Q(x)=−2(1+tan2x)
已知方程的一个非零特解:y1=tanx ,代入刘维尔公式得:
y2=−1−xtanx
接下来验证无关性:
W(y1,y2)=tanxsec2x−1−xtanx−tanx−xsec2x=tanx(−tanx−xsec2x)−sec2x(−1−xtanx)=−tan2x−xtanxsec2x+sec2x+xtanxsec2x=sec2x−tan2x=1=0
朗斯基行列式恒不为0,故 y1 与 y2 线性无关。
二阶线性齐次方程的通解为两个线性无关特解的线性组合,将 y2 的负号吸收到任意常数 C2 中,最终通解为:
y=C1tanx+C2(1+xtanx)(C1,C2 为任意常数)
PROOF2:(使用常数变易法)
对于标准形式的二阶线性齐次微分方程:
y′′+P(x)y′+Q(x)y=0且已知一个非零特解 y1(x)
常数变易法的核心是:
将通解中与 y1 搭配的常数,替换为关于 x 的未知函数 u(x),即设第二个线性无关特解为 y2=u(x)⋅y1(x)。
将其代入原方程后,利用 y1 是解的性质,可消去含 u(x) 的项,将二阶方程降为关于 u′(x) 的一阶微分方程,求解后即可得到 u(x),最终得到第二个特解与通解。
原方程:
y′′−2(1+tan2x)y=0
对比标准形式 y′′+P(x)y′+Q(x)y=0,得:
P(x)=0,Q(x)=−2(1+tan2x)
已知方程的一个非零特解:y1=tanx (补充三角恒等式:1+tan2x=sec2x,后续化简用)
根据常数变易法,设第二个与 y1 线性无关的特解为:
y2=u(x)⋅y1=u(x)⋅tanx
其中 u(x) 是待求的非常数函数(非常数才能保证与 y1 线性无关)。
对 y2 逐次求导(乘积求导法则):
-
一阶导数:
y2′=u′(x)tanx+u(x)sec2x
-
二阶导数:
y2′′=u′′(x)tanx+u′(x)sec2x+u′(x)sec2x+u(x)⋅2sec2xtanx=u′′(x)tanx+2u′(x)sec2x+2u(x)sec2xtanx
将 y2、y2′′ 代入原方程 y′′−2sec2x⋅y=0:
[u′′tanx+2u′sec2x+2usec2xtanx]−2sec2x⋅[utanx]=0
展开后观察项的抵消:含 u(x) 的项 2usec2xtanx 与 −2usec2xtanx 完全抵消(这是常数变易法降阶的核心,因为 y1 是原方程的解,含 u 的项必然抵消),最终化简为:
u′′tanx+2u′sec2x=0
此时二阶方程已降为关于 u′ 的一阶微分方程。
令 v=u′(x),则 u′′=v′,代入上式得:
v′tanx+2vsec2x=0
这是可分离变量的一阶微分方程,分离变量:
vdv=−2⋅tanxsec2xdx
两边同时积分(仅需一个特解,积分常数取0即可):
∫vdv=−2∫tanxsec2xdx
左边积分:
∫vdv=ln∣v∣
右边积分:令 t=tanx,则 dt=sec2xdx,因此
∫tanxsec2xdx=∫tdt=ln∣tanx∣
代入后整理:
ln∣v∣=−2ln∣tanx∣=ln∣cot2x∣
两边取指数,得:
v=u′(x)=cot2x
对 u′(x) 积分求 u(x)(积分常数仍取0):
u(x)=∫u′(x)dx=∫cot2xdx
利用三角恒等式 cot2x=csc2x−1 化简积分:
u(x)=∫(csc2x−1)dx=−cotx−x
即:
u(x)=−cotx−x
将 u(x) 代入 y2=u(x)⋅tanx,化简得:
y2=(−cotx−x)⋅tanx=−cotx⋅tanx−xtanx=−1−xtanx
接下来验证无关性:
W(y1,y2)=tanxsec2x−1−xtanx−tanx−xsec2x=tanx(−tanx−xsec2x)−sec2x(−1−xtanx)=−tan2x−xtanxsec2x+sec2x+xtanxsec2x=sec2x−tan2x=1=0
朗斯基行列式恒不为0,故 y1 与 y2 线性无关。
二阶线性齐次方程的通解为两个线性无关特解的线性组合,将 y2 的负号吸收到任意常数 C2 中,最终通解为:
y=C1tanx+C2(1+xtanx)(C1,C2 为任意常数)
*刘维尔公式的证明
看不懂没关系,会用公式就可以。
从二阶线性齐次微分方程的标准形式出发(保证推导的通用性,要求 y′′ 的系数为1):
y′′+P(x)y′+Q(x)y=0(1)
其中 P(x),Q(x) 是区间 I 上的连续函数;已知方程 (1) 在区间 I 上的一个非零特解 y1(x)(即 y1(x)≡0,x∈I)。
根据常数变易的核心思想,设与 y1(x) 线性无关的第二个特解为:
y2(x)=u(x)⋅y1(x)(2)
其中 u(x) 是待求的非常数函数(非常数才能保证 y2 与 y1 不成比例,即线性无关)。
利用乘积求导法则,对 (2) 式逐次求导:
-
一阶导数:
y2′=u′(x)y1(x)+u(x)y1′(x)(3)
-
二阶导数:
y2′′=u′′(x)y1(x)+2u′(x)y1′(x)+u(x)y1′′(x)(4)
将 (2)(3)(4) 代入原方程 (1),得:
[u′′y1+2u′y1′+uy1′′]+P(x)[u′y1+uy1′]+Q(x)[uy1]=0
将式子按 u′′,u′,u 分组整理:
u′′y1+u′[2y1′+P(x)y1]+u[y1′′+P(x)y1′+Q(x)y1]=0(5)
因为 y1(x) 是原方程 (1) 的特解,所以它必然满足原方程,即:
y1′′+P(x)y1′+Q(x)y1=0
因此 (5) 式中含 u 的项整体为0,直接抵消。此时二阶方程被降为关于 u′ 的一阶微分方程:
u′′y1+u′[2y1′+P(x)y1]=0(6)
令 v=u′(x),则 u′′=v′,代入 (6) 式得一阶齐次微分方程:
v′y1+v[2y1′+P(x)y1]=0
由于 y1(x)≡0,两边同时除以 y1,整理得:
v′+v(2⋅y1y1′+P(x))=0
这是可分离变量的微分方程,分离变量(v=0,否则 u 为常数,无意义):
vdv=−(2⋅y1y1′+P(x))dx(7)
对 (7) 式两边同时积分,我们仅需一个特解,因此积分常数取0(常数最终会被通解的任意常数吸收,不影响结果):
∫vdv=−∫(2⋅y1y1′+P(x))dx
分别计算左右两边积分:
-
左边:
∫vdv=ln∣v∣
-
右边:利用积分的线性性质拆分,其中
∫y1y1′dx=ln∣y1∣(复合函数积分法则)
因此:
−∫(2⋅y1y1′+P(x))dx=−2ln∣y1∣−∫P(x)dx
合并后得:
ln∣v∣=−2ln∣y1∣−∫P(x)dx
利用对数运算法则化简右边:
−2ln∣y1∣=ln∣y1−2∣=lny121
因此:
ln∣v∣=lny121−∫P(x)dx
两边同时取指数,消去对数:
v=y121⋅e−∫P(x)dx
由于 v=u′(x),因此我们得到了 u′(x) 的通用表达式:
u′(x)=y12(x)e−∫P(x)dx(8)
对 (8) 式两边再次积分,仍取积分常数为0(仅需一个特解),得:
u(x)=∫y12(x)e−∫P(x)dxdx
将 u(x) 代回我们最初设的特解形式 (2),最终得到第二个线性无关特解的通用计算公式,即刘维尔降阶公式:
y2(x)=y1(x)⋅∫y12(x)e−∫P(x)dxdx(9)
根据线性无关的充要条件,两个解的朗斯基行列式在区间内不恒为0,则线性无关。将 y1,y2 代入朗斯基行列式:
W(y1,y2)=y1y1′y2y2′=y1y2′−y1′y2=y1(u′y1+uy1′)−y1′⋅uy1=u′⋅y12
将 (8) 式的 u′ 代入,得:
W(y1,y2)=y12e−∫P(x)dx⋅y12=e−∫P(x)dx
由于指数函数恒不为0,因此 W(y1,y2)=0,即 y2 与 y1 必然线性无关,公式的有效性得证。
Q.E.D
希望这篇推文能帮到你:)
新年快乐~
LOVE FROM YUNRAN♡