第 079 篇

考研数学二·零基础讲义 第14讲(五)|常微分方程:非齐次特解(类型 Ⅱ)、欧拉方程与应用

约 11 分钟 3349 字 数学考研数学二零基础讲义高等数学

考研数学二·零基础讲义 · 第 14 讲(五)|常微分方程:非齐次特解(类型 Ⅱ)、欧拉方程与应用

学完本篇,你能为三角函数型自由项设出特解并求解,会用叠加原理处理"几项相加"的右边,了解欧拉方程的换元解法,并能解含变限积分的二阶方程应用题。预计用时:3 小时。

本讲地图

  • 第 14 讲拆成五篇:(一)基本概念、可分离变量与齐次方程;(二)一阶线性方程与应用;(三)可降阶方程与解的结构;(四)常系数齐次方程与非齐次特解(类型 Ⅰ);(五)非齐次特解(类型 Ⅱ)、欧拉方程与应用。
  • 非齐次方程的特解(类型 Ⅱ):f(x)=eαx(Plcos⁡βx+Pnsin⁡βx)f(x)=e^{\alpha x}(P_l\cos\beta x+P_n\sin\beta x)——大纲要求"会解自由项为多项式、指数函数、正弦函数、余弦函数以及它们的和与积"的情形。
  • 叠加原理:右边是几项之和时分别求特解再相加。
  • 欧拉方程(选学):说明:根据近年的数学二考试大纲原文,"会解欧拉方程"只写在数学一的要求里,数学二大纲没有列出;它只是"换元后变成常系数方程",这里作为选学简要介绍。
  • 简单应用:含变限积分的函数方程、初值问题。

前置知识

  • 第 14 讲(三)(四):解的结构、特征方程、类型 Ⅰ 特解。
  • 会对 xcos⁡2xx\cos2x、exsin⁡xe^x\sin x 这类函数求二阶导(见第 05、06 讲);变限积分求导(见第 10 讲)。

一句话心态:类型 Ⅱ 计算量大,但设法才是得分关键:cos⁡\cos、sin⁡\sin 都要设,kk 只有 0 和 1 两种。


一、类型 Ⅱ 特解的设法

1.1 大白话

上一篇的类型 Ⅰ 处理"多项式 × 指数"。这一篇处理右边带 cos⁡\cos、sin⁡\sin 的情形,比如 y′′+y=xcos⁡2xy''+y=x\cos2x、y′′+y′−2y=excos⁡xy''+y'-2y=e^x\cos x。思路完全一样:y′′+py′+qyy''+py'+qy 要等于 f(x)f(x),而 cos⁡βx\cos\beta x、sin⁡βx\sin\beta x 求导时会互相转化,所以特解里要同时放上 cos⁡βx\cos\beta x 和 sin⁡βx\sin\beta x;如果 α±iβ\alpha\pm i\beta 恰好是特征根(即设出的式子本身就是齐次解),就再乘一个 xx "躲开"。

1.2 设法(必背)

类型 Ⅱ:f(x)=eαx[Pl(x)cos⁡βx+Pn(x)sin⁡βx]f(x)=e^{\alpha x}\big[P_l(x)\cos\beta x+P_n(x)\sin\beta x\big]

y∗=xk eαx[Rm(x)cos⁡βx+Sm(x)sin⁡βx],y^*=x^k\,e^{\alpha x}\big[R_m(x)\cos\beta x+S_m(x)\sin\beta x\big],
  • m=max⁡{l,n}m=\max\{l,n\},RmR_m、SmS_m 是两个不同的 mm 次待定多项式;
  • kk:α±iβ\alpha\pm i\beta 不是特征根取 k=0k=0,是特征根取 k=1k=1。

1.3 逐行拆解

  • kk 的含义:如果设出的 y∗y^* 本身就是齐次解,代入左边会得到 0,不可能等于 ff——所以要乘 xx "躲开"它。
  • 类型 Ⅱ 两个三角函数都要设:即使 ff 里只有 cos⁡\cos,y∗y^* 里也要同时有 cos⁡\cos 和 sin⁡\sin,因为求导时它们会互相转化。
  • 类型 Ⅱ 中,sin⁡x\sin x 是 α=0,β=1\alpha=0,\beta=1;excos⁡2xe^x\cos2x 是 α=1,β=2\alpha=1,\beta=2。
  • 设法三问:① ff 属于哪一类?读出 λ\lambda(或 α,β\alpha,\beta)和多项式次数;② 求特征根,看 λ\lambda(或 α±iβ\alpha\pm i\beta)是不是特征根、几重;③ 写出 xkx^k 和待定多项式。
  • 怎么读出 α,β\alpha,\beta:把 f(x)f(x) 写成 eαx[⋯cos⁡βx+⋯sin⁡βx]e^{\alpha x}[\cdots\cos\beta x+\cdots\sin\beta x] 的样子。没有指数因子时 α=0\alpha=0。
  • 为什么 kk 最多是 1:二阶方程的复特征根总是一对单根 α±iβ\alpha\pm i\beta,不会出现二重复根。

二、类型 Ⅱ 例题与叠加原理

2.1 例题

例 1(类型 Ⅱ,不是特征根) 求 y′′+y=xcos⁡2xy''+y=x\cos2x 的一个特解。

解: ① α=0, β=2\alpha=0,\ \beta=2,Pl=xP_l=x(l=1l=1),Pn=0P_n=0,m=1m=1。特征根 ±i\pm i,α±iβ=±2i\alpha\pm i\beta=\pm2i 不是特征根,k=0k=0。 ② 设 y∗=(ax+b)cos⁡2x+(cx+d)sin⁡2xy^*=(ax+b)\cos2x+(cx+d)\sin2x(四个待定系数)。 ③ 分别计算:记 g=(ax+b)cos⁡2xg=(ax+b)\cos2x,h=(cx+d)sin⁡2xh=(cx+d)\sin2x。

g′′+g=−3(ax+b)cos⁡2x−4asin⁡2x,h′′+h=4ccos⁡2x−3(cx+d)sin⁡2x.g''+g=-3(ax+b)\cos2x-4a\sin2x,\qquad h''+h=4c\cos2x-3(cx+d)\sin2x .

(以 gg 为例:g′=acos⁡2x−2(ax+b)sin⁡2xg'=a\cos2x-2(ax+b)\sin2x,g′′=−4asin⁡2x−4(ax+b)cos⁡2xg''=-4a\sin2x-4(ax+b)\cos2x,再加 gg。) ④ 相加并与 xcos⁡2xx\cos2x 比较:

  • cos⁡2x\cos2x 的系数:−3ax−3b+4c=x-3ax-3b+4c=x ⇒ −3a=1-3a=1,−3b+4c=0-3b+4c=0;
  • sin⁡2x\sin2x 的系数:−4a−3cx−3d=0-4a-3cx-3d=0 ⇒ −3c=0-3c=0,−4a−3d=0-4a-3d=0。

解得 a=−13a=-\dfrac13,c=0c=0,b=0b=0,d=49d=\dfrac49。

y∗=−x3cos⁡2x+49sin⁡2x.y^*=-\frac x3\cos2x+\frac49\sin2x .

例 2(类型 Ⅱ,是特征根) 求 y′′+y=sin⁡xy''+y=\sin x 的一个特解。

解:α=0, β=1\alpha=0,\ \beta=1,α±iβ=±i\alpha\pm i\beta=\pm i 正是特征根,k=1k=1,m=0m=0。设 y∗=x(acos⁡x+bsin⁡x)y^*=x(a\cos x+b\sin x)。 记 u=acos⁡x+bsin⁡xu=a\cos x+b\sin x(注意 u′′=−uu''=-u),则 (xu)′′=2u′+xu′′(xu)''=2u'+xu'',所以

(xu)′′+xu=2u′=2(−asin⁡x+bcos⁡x)=sin⁡x.(xu)''+xu=2u'=2(-a\sin x+b\cos x)=\sin x .

比较:−2a=1-2a=1,b=0b=0,a=−12a=-\dfrac12。特解 y∗=−x2cos⁡xy^*=-\dfrac x2\cos x。

例 3(类型 Ⅱ,α≠0\alpha\ne0) 求 y′′+y′−2y=excos⁡xy''+y'-2y=e^x\cos x 的一个特解。

解:α=1, β=1\alpha=1,\ \beta=1。特征根 r=1,−2r=1,-2,1±i1\pm i 不是特征根,k=0k=0,m=0m=0。 设 y∗=ex(acos⁡x+bsin⁡x)y^*=e^x(a\cos x+b\sin x)。代入比较系数(计算较长,思路同例 1),解得 a=−110a=-\dfrac1{10},b=310b=\dfrac{3}{10}:

y∗=ex(310sin⁡x−110cos⁡x).y^*=e^x\left(\frac{3}{10}\sin x-\frac{1}{10}\cos x\right).

例 4(叠加原理) 求 y′′−y=ex+xy''-y=e^x+x 的通解。

解:特征根 ±1\pm1,齐次通解 C1ex+C2e−xC_1e^x+C_2e^{-x}。把右边拆成两部分分别求特解:

  • f1=exf_1=e^x:λ=1\lambda=1 是单根,设 y1∗=axexy_1^*=axe^x。(axex)′′=a(x+2)ex(axe^x)''=a(x+2)e^x,代入 a(x+2)ex−axex=2aex=exa(x+2)e^x-axe^x=2ae^x=e^x,a=12a=\frac12,y1∗=x2exy_1^*=\dfrac x2e^x。
  • f2=xf_2=x:λ=0\lambda=0 不是根,设 y2∗=bx+cy_2^*=bx+c,代入 0−(bx+c)=x0-(bx+c)=x,b=−1, c=0b=-1,\ c=0,y2∗=−xy_2^*=-x。
y=C1ex+C2e−x+x2ex−x.y=C_1e^x+C_2e^{-x}+\frac x2e^x-x .

例 5(只设形式,不求系数) 写出下列方程特解的形式: (1) y′′−3y′+2y=x2e3xy''-3y'+2y=x^2e^{3x};(2) y′′−3y′+2y=xexy''-3y'+2y=xe^{x};(3) y′′−2y′+5y=exsin⁡2xy''-2y'+5y=e^x\sin2x;(4) y′′+4y=cos⁡xy''+4y=\cos x。

解: (1) 特征根 1,21,2;λ=3\lambda=3 不是根,m=2m=2:y∗=(ax2+bx+c)e3xy^*=(ax^2+bx+c)e^{3x}。 (2) λ=1\lambda=1 是单根:y∗=x(ax+b)exy^*=x(ax+b)e^x。 (3) 特征根 1±2i1\pm2i;α±iβ=1±2i\alpha\pm i\beta=1\pm2i 是特征根:y∗=xex(acos⁡2x+bsin⁡2x)y^*=xe^x(a\cos2x+b\sin2x)。(实际解出 y∗=−x4excos⁡2xy^*=-\dfrac x4e^x\cos2x。) (4) 特征根 ±2i\pm2i;±i\pm i 不是特征根:y∗=acos⁡x+bsin⁡xy^*=a\cos x+b\sin x。

2.2 易错点 ⚠️

  • kk 定错是第一大错误:一定先求特征根,再判断 λ\lambda 是不是根、几重。
  • 待定多项式要设完整:f=xf=x 时设 ax+bax+b,不能只设 axax。
  • 类型 Ⅱ 中 cos⁡\cos、sin⁡\sin 都要设,且系数多项式各自独立。
  • 求的是通解时,别忘了加上齐次通解。

三、完整解题流程与设法速查表

3.1 流程与速查

解一道"求 y′′+py′+qy=f(x)y''+py'+qy=f(x) 的通解",按这五步走,一步都不要省:

  1. 写特征方程,求出特征根;
  2. 写齐次通解 YY(三种根对号入座);
  3. 看 f(x)f(x) 的类型,读出 λ\lambda(或 α,β\alpha,\beta)和多项式的次数,判断 kk,设出 y∗y^*;
  4. 代入比较系数,解出待定系数,得到 y∗y^*;
  5. 写通解 y=Y+y∗y=Y+y^*;有初始条件的,最后再定 C1,C2C_1,C_2(注意:初始条件要代入通解,不能只代入齐次通解)。

下面的速查表以 y′′−3y′+2y=f(x)y''-3y'+2y=f(x)(特征根 1,21,2)为例,帮你练"设法":

f(x)f(x) 读出 是否特征根 y∗y^* 的形式
55 λ=0\lambda=0,m=0m=0 否,k=0k=0 aa
x2x^2 λ=0\lambda=0,m=2m=2 否 ax2+bx+cax^2+bx+c
e−xe^{-x} λ=−1\lambda=-1,m=0m=0 否 ae−xae^{-x}
exe^{x} λ=1\lambda=1 单根,k=1k=1 axexaxe^{x}
(x+1)e2x(x+1)e^{2x} λ=2\lambda=2,m=1m=1 单根,k=1k=1 x(ax+b)e2xx(ax+b)e^{2x}
sin⁡x\sin x α=0,β=1\alpha=0,\beta=1 ±i\pm i 否 acos⁡x+bsin⁡xa\cos x+b\sin x
excos⁡xe^x\cos x α=1,β=1\alpha=1,\beta=1 1±i1\pm i 否 ex(acos⁡x+bsin⁡x)e^x(a\cos x+b\sin x)
ex+xe^x+x 拆成两项 分别判断 axex+(bx+c)axe^x+(bx+c)

可以看到,同一个方程,右边换一下,kk 就可能变。所以每次都要拿 λ\lambda 去和特征根比一比。


四、欧拉方程(选学)

4.1 大白话

方程形如 x2y′′+pxy′+qy=f(x)x^2y''+pxy'+qy=f(x):系数里 xx 的次数恰好等于导数的阶数。它不是常系数的,但换元 x=etx=e^t 后就变成常系数方程。数学二大纲没有列出欧拉方程,这里只简单介绍思路。

4.2 方法

令 x=etx=e^t(x>0x>0),记 D=ddtD=\dfrac{d}{dt},则

xy′=Dy,x2y′′=D(D−1)y.xy'=Dy,\qquad x^2y''=D(D-1)y .

代入后得到关于 tt 的常系数线性方程,解出 y(t)y(t) 后把 t=ln⁡xt=\ln x 代回去。

4.3 例题

例 6 求 x2y′′−2xy′+2y=0x^2y''-2xy'+2y=0 的通解(x>0x>0)。

解:代换后 D(D−1)y−2Dy+2y=0D(D-1)y-2Dy+2y=0,即 (D2−3D+2)y=0(D^2-3D+2)y=0,特征方程 r2−3r+2=0r^2-3r+2=0,r=1,2r=1,2。y=C1et+C2e2ty=C_1e^t+C_2e^{2t},代回 et=xe^t=x:

y=C1x+C2x2.y=C_1x+C_2x^2 .

例 7 求 x2y′′+xy′−y=x2x^2y''+xy'-y=x^2 的通解(x>0x>0)。

解:D(D−1)y+Dy−y=e2tD(D-1)y+Dy-y=e^{2t},即 d2ydt2−y=e2t\dfrac{d^2y}{dt^2}-y=e^{2t}。特征根 ±1\pm1;λ=2\lambda=2 不是根,设 y∗=ae2ty^*=ae^{2t}:4a−a=14a-a=1,a=13a=\frac13。 y=C1et+C2e−t+13e2ty=C_1e^t+C_2e^{-t}+\dfrac13e^{2t},代回:

y=C1x+C2x+x23.y=C_1x+\frac{C_2}{x}+\frac{x^2}{3}.

五、简单应用

5.1 大白话

二阶方程的应用题,考研里最常见的是"含变限积分的函数方程":求导两次变成二阶常系数方程,再从原式和一阶导式中代 x=0x=0 找出两个初始条件。

5.2 例题

例 8 设 f(x)f(x) 连续,且 f(x)=sin⁡x−∫0x(x−t)f(t) dtf(x)=\sin x-\displaystyle\int_0^x(x-t)f(t)\,dt,求 f(x)f(x)。

解: ① 先拆开积分(被积函数含 xx,不能直接求导):

f(x)=sin⁡x−x∫0xf(t) dt+∫0xtf(t) dt.f(x)=\sin x-x\int_0^xf(t)\,dt+\int_0^xtf(t)\,dt .

② 求一阶导:−x∫0xf dt-x\int_0^xf\,dt 用乘积法则,导数为 −∫0xf dt−xf(x)-\int_0^xf\,dt-xf(x);∫0xtf(t)dt\int_0^xtf(t)dt 的导数为 xf(x)xf(x)。所以

f′(x)=cos⁡x−∫0xf(t) dt.f'(x)=\cos x-\int_0^xf(t)\,dt .

③ 求二阶导:f′′(x)=−sin⁡x−f(x)f''(x)=-\sin x-f(x),即

f′′+f=−sin⁡x.f''+f=-\sin x .

④ 初始条件:原式代 x=0x=0 得 f(0)=0f(0)=0;②式代 x=0x=0 得 f′(0)=1f'(0)=1。 ⑤ 解方程:齐次通解 C1cos⁡x+C2sin⁡xC_1\cos x+C_2\sin x。由例 2,y′′+y=sin⁡xy''+y=\sin x 的特解是 −x2cos⁡x-\frac x2\cos x,所以 y′′+y=−sin⁡xy''+y=-\sin x 的特解是 x2cos⁡x\frac x2\cos x。

f=C1cos⁡x+C2sin⁡x+x2cos⁡x.f=C_1\cos x+C_2\sin x+\frac x2\cos x .

f(0)=C1=0f(0)=C_1=0;f′=C2cos⁡x+12cos⁡x−x2sin⁡xf'=C_2\cos x+\frac12\cos x-\frac x2\sin x(C1=0C_1=0 后),f′(0)=C2+12=1f'(0)=C_2+\frac12=1,C2=12C_2=\frac12。

f(x)=12sin⁡x+x2cos⁡x.f(x)=\frac12\sin x+\frac x2\cos x .

例 9(初值问题) 求 y′′+y=sin⁡xy''+y=\sin x 满足 y(0)=0, y′(0)=0y(0)=0,\ y'(0)=0 的特解。

解:通解 y=C1cos⁡x+C2sin⁡x−x2cos⁡xy=C_1\cos x+C_2\sin x-\frac x2\cos x。y(0)=C1=0y(0)=C_1=0。y′=C2cos⁡x−12cos⁡x+x2sin⁡xy'=C_2\cos x-\frac12\cos x+\frac x2\sin x,y′(0)=C2−12=0y'(0)=C_2-\frac12=0,C2=12C_2=\frac12。

y=12sin⁡x−x2cos⁡x.y=\frac12\sin x-\frac x2\cos x .

(物理上这叫"共振":外力的频率和系统固有频率相同,振幅 x2\frac x2 随时间越来越大。)

共振:外力频率等于固有频率时,振幅随 x 线性增长 x y 2π 4π 6π 5 −5 O y″ + y = sin x,y(0) = 0,y′(0) = 0 的解 y = (1/2)sin x − (x/2)cos x 橙色虚线 y = ±x/2:振幅越来越大(共振)

图 1:例 9 的解 y = (1/2)sin x − (x/2)cos x:振幅沿 ±x/2 越来越大(动画循环播放)

5.3 易错点 ⚠️

  • 被积函数里含 xx 时,先把 xx 提到积分号外再求导。
  • 两个初始条件分别来自原式和一阶导式,别只找一个。

公式/结论速记卡

主题 核心结论
特解 Ⅱ f=eαx(Plcos⁡βx+Pnsin⁡βx)f=e^{\alpha x}(P_l\cos\beta x+P_n\sin\beta x):y∗=xkeαx(Rmcos⁡βx+Smsin⁡βx)y^*=x^ke^{\alpha x}(R_m\cos\beta x+S_m\sin\beta x),m=max⁡{l,n}m=\max\{l,n\}
kk 的取法 α±iβ\alpha\pm i\beta 不是特征根 k=0k=0,是特征根 k=1k=1
共振小结论 y′′+ω2y=sin⁡ωxy''+\omega^2y=\sin\omega x 的特解 −x2ωcos⁡ωx-\frac{x}{2\omega}\cos\omega x;cos⁡ωx\cos\omega x 对应 x2ωsin⁡ωx\frac{x}{2\omega}\sin\omega x
叠加原理 右边是几项之和,就分别求特解再相加
欧拉(选学) x=etx=e^t:xy′=Dyxy'=Dy,x2y′′=D(D−1)yx^2y''=D(D-1)y
变限积分方程 先把 xx 提出积分号,求导两次,原式和一阶导式各代 x=0x=0

自测练习

1. 写出 y′′+4y=xsin⁡2xy''+4y=x\sin2x 的特解形式,并(选做)求出特解。

答案与提示

α=0,β=2\alpha=0,\beta=2,m=1m=1;特征根 ±2i\pm2i,±2i\pm2i 是特征根,k=1k=1:

y∗=x[(ax+b)cos⁡2x+(cx+d)sin⁡2x].y^*=x\big[(ax+b)\cos2x+(cx+d)\sin2x\big].

代入比较系数可得 a=−18a=-\frac18,b=0b=0,c=0c=0,d=116d=\frac1{16},即 y∗=−x28cos⁡2x+x16sin⁡2xy^*=-\dfrac{x^2}{8}\cos2x+\dfrac{x}{16}\sin2x。

2. 求 y′′+y=cos⁡xy''+y=\cos x 的一个特解。

答案与提示

±i\pm i 是特征根,k=1k=1:设 y∗=x(acos⁡x+bsin⁡x)y^*=x(a\cos x+b\sin x),记 u=acos⁡x+bsin⁡xu=a\cos x+b\sin x,则 (xu)′′+xu=2u′=2(−asin⁡x+bcos⁡x)=cos⁡x(xu)''+xu=2u'=2(-a\sin x+b\cos x)=\cos x,得 a=0a=0,b=12b=\frac12。 y∗=x2sin⁡xy^*=\dfrac x2\sin x。

3. 求 y′′−y=sin⁡xy''-y=\sin x 的一个特解。

答案与提示

特征根 ±1\pm1,±i\pm i 不是特征根,k=0k=0:设 y∗=acos⁡x+bsin⁡xy^*=a\cos x+b\sin x,则 y∗′′−y∗=−2acos⁡x−2bsin⁡x=sin⁡xy^{*\prime\prime}-y^*=-2a\cos x-2b\sin x=\sin x,a=0a=0,b=−12b=-\frac12。 y∗=−12sin⁡xy^*=-\dfrac12\sin x。

4. 求 y′′+y=ex+1y''+y=e^x+1 的一个特解。

答案与提示

叠加:f1=exf_1=e^x,λ=1\lambda=1 不是根,设 aexae^x:2a=12a=1,y1∗=12exy_1^*=\frac12e^x;f2=1f_2=1,λ=0\lambda=0 不是根,设 bb:b=1b=1。 y∗=12ex+1y^*=\dfrac12e^x+1。

5. 设 f(x)f(x) 连续,且 f(x)=x+∫0x(t−x)f(t) dtf(x)=x+\displaystyle\int_0^x(t-x)f(t)\,dt,求 f(x)f(x)。

答案与提示

拆开:f(x)=x+∫0xtf(t)dt−x∫0xf(t)dtf(x)=x+\int_0^xtf(t)dt-x\int_0^xf(t)dt。求导:f′(x)=1+xf(x)−∫0xf dt−xf(x)=1−∫0xf(t)dtf'(x)=1+xf(x)-\int_0^xf\,dt-xf(x)=1-\int_0^xf(t)dt。再求导:f′′=−ff''=-f,即 f′′+f=0f''+f=0。 初始条件:f(0)=0f(0)=0,f′(0)=1f'(0)=1。通解 C1cos⁡x+C2sin⁡xC_1\cos x+C_2\sin x,得 C1=0, C2=1C_1=0,\ C_2=1,f(x)=sin⁡xf(x)=\sin x。

6.(选学)求欧拉方程 x2y′′+3xy′+y=0x^2y''+3xy'+y=0 的通解(x>0x>0)。

答案与提示

D(D−1)y+3Dy+y=0D(D-1)y+3Dy+y=0,即 (D2+2D+1)y=0(D^2+2D+1)y=0,r=−1r=-1(二重)。y=(C1+C2t)e−t=C1+C2ln⁡xxy=(C_1+C_2t)e^{-t}=\dfrac{C_1+C_2\ln x}{x}。


下一讲预告

高等数学部分到此结束。第 15 讲开始线性代数,从行列式讲起:它是什么、怎么算、有哪些让计算变简单的性质。

评论

还没有评论。 来写第一条。

留下一条评论

选填,不会公开

纯文本,链接会自动识别。