陕西邮电职业技术学院-调查显示

`
第十二章：微分方程
教学目的：
1．了解微分方程及其解、阶、通解，初始条件和特等概念。
2．熟练掌握变量可分离的微分方程及一阶线性微分方程的解法。
3．会解齐次微分方程、伯努利方程和全微分方程，会用简单的变量代换解某些微分方程。
4.会用降阶法解下列微分方程：
y
(n)
f(x)
，
y

f(x,y

)
和
y

f( y,y

)

5.理解线性微分方程解的性质及解的结构定理。
6 ．掌握二阶常系数齐次线性微分方程的解法，并会解某些高于二阶的常系数齐次线性微
分方程。
7.求自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性
微分方程的 特解和通解。
8.会解欧拉方程，会解包含两个未知函数的一阶常系数线性微分方程组。
9．会解微分方程组（或方程组）解决一些简单的应用问题。
教学重点：
1、可分离的微分方程及一阶线性微分方程的解法
2、可降阶的高阶微分方程
y
(n)
f(x)
，
y
f(x,y

)
和
y

f(y,y

)

3、二阶常系数齐次线性微分方程；
4、自由项为多项式、 指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性
微分方程；
教学难点：
1、齐次微分方程、伯努利方程和全微分方程；
2、线性微分方程解的性质及解的结构定理；
3、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方
程的特解。
4、欧拉方程
Word文档

`
§12 1 微分方程的基本概念

函数是客观事物的部联系在数量方面的反映 利用函数关系又可以对客观事物的规律性进
行研究 因此如何寻找出所需要的函数关系 在实践中具有重要意义 在许多问题中 往往不能
直接找出所需要的函数关系 但是根据问题所提供的情况 有时可以列出含有要找的函数及其
导数的关系式 这样的关系就是所谓微分方程 微分方程建立以后 对它进行研究 找出未知函
数来 这就是解微分方程
几个概念
微分方程 表示未知函数、未知函数的导数与自变量之间的关系的方程 叫微分方程
常微分方程 未知函数是一元函数的微分方程 叫常微分方程
偏微分方程 未知函数是多元函数的微分方程 叫偏微分方程
微分方程的阶 微分方程中所出现的未知函数的最高阶导数的阶数 叫微分方程的阶
x
3
yx
2
y4xy3x
2

y
(4)
4y10y12y5ysin2x
y
(
n
)
10
一般n阶微分方程
F(x y y     y
(
n
)
)0
y
(
n
)
f(x y y     y
(
n
1)
) 
微分方程的解 满足微分方程的函数(把函数代入微分方程能使该方程成为恒等式)叫做该
微分方程的解 确切地说 设函数y

(x)在区间I上有n阶连续导数 如果在区间I上
F[x

(x)

(x)   

(
n
)
(x)]0
那么函数y

(x)就叫做微分方程F(x y y    y
(
n
)
)0在区间I上的解
通解 如果微分方程的解中含有任意常数 且任意常数的个数与微分方程的阶数相同 这样
的解叫做微分方程的通解
初始条件 用于确定通解中任意常数的条件 称为初始条件 如
xx
0
时 yy
0
 y y
0

一般写成


y
xx
0
y
0

y

xx
0
y
0
特解 确定了通解中的任意常数以后 就得到微分方程的特解 即不含任意常数的解
初值问题 求微分方程满足初始条件的解的问题称为初值问题
Word文档

`
如求微分方程yf(x y)满足初始条件
y
xx
0
y
0
的解的问题 记为

y

f(x,y)



y
xx
0
y
0

积分曲线 微分方程的解的图形是一条曲线 叫做微分方程的积分曲线
例1 一曲线通过点(1 2) 且在该曲线上任一点M(x y)处的切线的斜率为2x 求这曲线的方
程
解 设所求曲线的方程为yy(x) 根据导数的几何意义 可知未知函数yy(x)应满足关系式
(称为微分方程)

dy
2x
 (1)
dx
此外 未知函数yy(x)还应满足下列条件
x1时 y2 简记为y|
x1
2 (2)
把(1)式两端积分 得(称为微分方程的通解)

y2xdx
 即yx
2
C (3)
其中C是任意常数
把条件“x1时 y2”代入(3)式 得
21
2
C
由此定出C1 把C1代入(3)式 得所求曲线方程(称为微分方程满足条件y|
x

1
2的解)
yx
2
1
例2 列车在平直线路上以20ms(相当于72kmh)的速度行驶 当制动时列车获得加速度
04ms
2
 问开始制动后多少时间列车才能停住 以及列车在这段时间里行驶了多少路程?
解 设列车在开始制动后t秒时行驶了s米 根据题意 反映制动阶段列车运动规律的函数
ss(t)应满足关系式
2
ds
0.4
 (4)
dt
2

此外 未知函数ss(t)还应满足下列条件
t0时 s0
v
ds
20
 简记为s|=0 s|=20 (5)
t0t0
dt
把(4)式两端积分一次 得
Word文档

`

v
ds
0.4tC
 (6)
1
dt
再积分一次 得
s02t
2
C
1
t C
2
 (7)
这里C
1
 C
2
都是任意常数
把条件v|
t0
20代入(6)得
20C
1

把条件s|
t0
0代入(7)得0C
2

把C
1
 C
2
的值代入(6)及(7)式得
v04t 20 (8)
s02t
2
20t (9)
在(8)式中令v0 得到列车从开始制动到完全停住所需的时间

t
20
50
(s)
0.4
再把t50代入(9) 得到列车在制动阶段行驶的路程
s0250
2
2050500(m)

解 设列车在开始制动后t秒时行驶了s米
s04 并且s|
t0
=0 s|
t0
=20
把等式s04两端积分一次 得
s04tC
1
 即v04tC
1
(C
1
是任意常数)
再积分一次 得
s02t
2
C
1
t C
2
(C
1
 C
2
都C
1
是任意常数)
由v|
t0
20得20C
1
 于是v04t 20
由s|
t0
0得0C
2
 于是s02t
2
20t
令v0 得t50(s) 于是列车在制动阶段行驶的路程
s0250
2
2050500(m)
例3 验证 函数
xC
1
cos ktC
2
sin kt
是微分方程
Word文档

`
2
dx
k
2
x0

dt
2
的解
解 求所给函数的导数

dx
kCsinktkCcoskt

12
dt
2
dx
k
2
Ccosktk
2
Csinktk< br>2
(CcosktCsinkt)

1212
2
dt
2
dx
将
2
及x的表达式代入所给方程 得
dt
k
2
(C
1
cos ktC
2
sin kt) k
2
(C
1
cos ktC
2
sin kt)0
2
dx
这表明函数xC
1
cosktC
2
sinkt 满足方程
2
k
2
x0
 因此所给函数是所给方程的解
dt
2
dx
例4 已知函数xC
1
coskt C
2
sinkt(k0)是微分方程
2
k
2
x0的通解 求满足初始条件
dt
x|
t0
A x|
t0
0
的特解
解 由条件x|
t0
A及xC
1
cos ktC
2
sin kt 得
C
1
A
再由条件x|
t0
0 及x(t) kC
1
sin ktkC
2
cos kt 得
C
2
0
把C
1
、C
2
的值代入xC
1
cos ktC
2
sin kt中 得
xAcos kt

§12 2 可分离变量的微分方程
观察与分析
1 求微分方程y2x的通解 为此把方程两边积分 得
yx
2
C
一般地 方程yf(x)的通解为
yf(x)dxC
(此处积分后不再加任意常数)
2 求微分方程y2xy
2
的通解

Word文档

`
因为y是未知的 所以积分
2xy
2
dx
无法进行 方程两边直

接积分不能求出通解
为求通解可将方程变为


1
dy2xdx
 两边积分 得
y
2
1
x
2
C
 或
y
2
1

y
xC
可以验证函数
y
1
是原方程的通解
x
2
C
一般地 如果一阶微分方程y

(x, y)能写成
g(y)dyf(x)dx
形式 则两边积分可得一个不含未知函数的导数的方程
G(y)F(x)C
由方程G(y)F(x)C所确定的隐函数就是原方程的通解
对称形式的一阶微分方程
一阶微分方程有时也写成如下对称形式
P(x y)dxQ(x y)dy0
在这种方程中 变量x与y 是对称的
若把x看作自变量、y看作未知函数 则当Q(x,y)0时 有

dyP(x,y)


dxQ(x,y)
dx

Q(x,y)

dyP(x,y)
若把y看作自变量、x看作未知函数 则当P(x,y)0时 有

可分离变量的微分方程
如果一个一阶微分方程能写成
g(y)dyf(x)dx (或写成y

(x)

(y))
的形式 就是说 能把微分方程写成一端只含y的函数和dy 另一端只含x的函数和dx 那么原
方程就称为可分离变量的微分方程
讨论 下列方程中哪些是可分离变量的微分方程？
(1) y2xy 是 y
1
dy2xdx 
(2)3x
2
5xy0 是 dy(3x
2
5x)dx
Word文档

`
(3)(x
2
y
2
)dxxydy=0 不是
(4)y1xy
2
xy
2
 是 y(1x)(1y
2
)
(5)y10
x

y
 是 10

y
dy10
x
dx
(6)
y


x

y
 不是
yx
可分离变量的微分方程的解法
第一步 分离变量 将方程写成g(y)dy f(x)dx的形式
第二步 两端积分
g(y)dyf(x)dx
 设积分后得G(y)F(x)C
第三步 求出由G(y)F(x)C所确定的隐函数y

(x)或x

(y)
G(y)F(x)C  y

(x)或x

(y)都是方程的通解 其中G(y)F(x)C称为隐式(通)解
例1 求微分方程

dy
2xy
的通解
dx
解 此方程为可分离变量方程 分离变量后得

1
dy2xdx

y
1
两边积分得


y
dy

2xdx

2
即 ln|y|x
2
C
1

从而
yexC
1
e
C
1
e
x

2
因为
e
C
1
仍是任意常数 把它记作C 便得所给方程的通解

yCe
x

解 此方程为可分离变量方程 分离变量后得

2
1
dy2xdx

y
1
dy2xdx


y

两边积分得

即 ln|y|x
2
lnC
Word文档

`
从而
yCe
x


例2 铀的衰变速度与当时未衰变的原子的含量M成正比 已知t0时铀的含量为M
0
 求在
衰变过程中铀含量M(t)随时间t变化的规律
解 铀的衰变速度就是M(t)对时间t的导数
2
dM

dt
由于铀的衰变速度与其含量成正比 故得微分方程

dM


M

dt
dM
0

dt
其中

(

>0)是常数

前的曲面号表示当t增加时M单调减少 即
由题意 初始条件为
M|
t0
M
0

将方程分离变量得

dM


dt

M
dM
(

)dt


M

两边积分 得
即 lnM

tlnC 也即MCe


t

由初始条件 得M
0
Ce
0
C
所以铀含量M( t)随时间t变化的规律MM
0
e


t


例3 设降落伞从跳伞塔下落后 所受空气阻力与速度成正比 并设降落伞离开跳伞塔时速
度为零 求降落伞下落速度与时间的函数关系
解 设降落伞下落速度为v(t) 降落伞所受外力为Fmgkv( k为比例系数) 根据牛顿第二
运动定律Fma 得函数v(t)应满足的方程为

m
dv
mgkv

dt
初始条件为
v|
t0
0
方程分离变量 得
Word文档

`

dv

dt

mgkvm
dv

dt


mgkv

m
两边积分 得

ln(mgkv)
1
k
t
C

m1
kC
1

k
t
mg
e
m
即
v
(
C
)
Ce
k
k
mg
将初始条件v|
t0
0代入通解得
C

k

k
t
mg
(1e
m
)
 于是降落伞下落速度与时间的函数关系为
v
k
dy
例4 求微分方程
1xy
2
xy
2
的通解
dx
解 方程可化为

dy
(1x)(1y
2
)

dx
1
dy(1x)dx

1y
2
分离变量得

两边积分得

1
dy(1x)dx
 即
1
x
2
xC

arctany
1y
2

2
于是原方程的通解为
ytan(x
2< br>xC)

例5有高为1m的半球形容器 水从它的底部小孔流出 小孔横截面面积为1cm
2
 开始时容
器盛满了水 求水从小孔流出过程中容器里水面高度h随时间t变化的规律
解 由水力学知道 水从孔口流出的流量Q可用下列公式计算

Q
1
2
dV
0.62S2gh

dt
其中0 62为流量系数 S为孔口横截面面积 g为重力加速度 现在孔口横截面面积S1cm
2
 故

dV
0.622gh
 或
dV0.622ghdt

dt
另一方面 设在微小时间间隔[t tdt] 水面高度由h降至hdh(dh0) 则又可得到
Word文档

`
dV

r
2
dh
其中r是时刻t的水面半径 右端置负号是由于dh0而dV0的缘故 又因

r100
2
(100h)
2
200hh
2

所以 dV

(200hh
2
)dh
通过比较得到

0.622ghdt

(200hh
2
)dh

这就是未知函数hh(t)应满足的微分方程
此外 开始时容器的水是满的 所以未知函数hh(t)还应满足下列初始条件
h|
t0
100
将方程
0.622ghdt

(200hh
2
)dh
分离变量后得

dt
两端积分 得

t

0.6 22g
13
(200h
2
h
2
)dh
 
0.622g

13
(200h
2
h
2< br>)dh

即
t(
400
h
2
2
h
2
)C

5
0.622g
3
其中C是任意常数
由初始条件得

t(
400
100
2

2< br>100
2
)C

5
0.622g
3

C

35

35

(
400000
200000
)
14
10
5

35
0.622g0.622g
15

因此
t
Word文档

0.622g
(710
5
35
3
2
10h3h
2
)


`
上式表达了水从小孔流出的过程中容器水面高度h与时间t之间的函数关系
§12 3 齐次方程
齐次方程
如果一阶微分方程
dy
f(x,y)
中的函数f(x, y)可写成
dx
y
y
的函数 即
f(x,y)

()
 则称这方程为齐次方程
x
x
下列方程哪些是齐次方程？
dy
yy
2< br>x
2
dyyy
()
2
1
 (1)
xy

yyx0
是齐次方程

dxxdx xx
22
(2)
1xy

1y
2
不 是齐次方程

2
dy1y
2


dx
1x
2
(3)(x
2
y
2
)dxxydy0
dyx
2
y
2
dy
x
y

 是齐次方程

dxxydxyx
(4) (2xy4)dx(xy1)dy0不是齐次方程

dy2xy4


dxxy1
(5)
(2xsh3ych)dx3xchdy0
是齐次方程
y
x
y
x
y
x
yy
2xsh3ych
dy
xx

dy

2
th
y

y




y
dxdx3xx
3xch
x

齐次方程的解法
y
dyy


()
中 令
u
 即yux 有
x
dxx
du


(u)

ux
dx
在齐次方程
分离变量 得

du

dx


(u)ux
两端积分 得
Word文档

`

du

dx



(u)u

x
y
代替u 便得所给齐次方程的通解
x
dydy
例1 解方程
y
2
x
2
xy

dxdx
求出积分后 再用
解 原方程可写成
y
2
()
dyy
x

 < br>dx
xyx
2
y
1
x
2
因此原方程是齐 次方程 令
yux
于是原方程变为
2
duu

ux


dxu1
du

u
 即
x
dxu1
y
u
 则
x
dy
ux
du

dxdx
分离变量 得

(1)du
1
u
dx

x
两边积分 得uln|u|Cln|x|
或写成ln|xu|uC
以
y
代上式中的u 便得所给方程的通解
x

ln|y|
y
C

x
例2 有旋转曲面形状的凹镜 假设由旋转轴上一点O发出的一切光线经此凹镜反射后都与
旋转轴平行 求这旋转曲面的方程
解 设此凹镜是由xOy面上曲线L yy(x)(y>0)绕x轴旋转而成 光源在原点 在L上任取一点
M(x, y) 作L的切线交x轴于A 点O发出的光线经点M反射后是一条平行于x轴射线 由光学及几
何原理可以证明OAOM
Word文档

`
因为
OAAPOPPMcot

OP
而
OMx
2
y
2

于是得微分方程
y
x

y

y
xx
2
y
2

y

整理得
dx

x
(
x
)
2< br>1
 这是齐次方程
dyyy
dx

x
(
x
)
2
1

dyyy
问题归结为解齐次方程
令
即
y
x
v
 即xyv 得
vy
dv
vv
2
1

y
dy
dv
v
2
1

dy
dv

dy

v
2
1
y
分离变量 得
两边积分 得
ln(vv
2
1)lnylnC
,
vv
2
1
y
y
,
(v)
2
v
2
1
,
CC
y
2
2yv
1

C
2
C
以yvx代入上式 得
y
2
2C(x
C
)

2
这是以x轴为轴、焦点在原点的抛物线 它绕x轴旋转所得旋转曲面的方程为

y
2
z
2
2C(x
C
)

2
这就是所求的旋转曲面方程 
例3 设一条河的两岸为平行直线 水流速度为a 有一鸭子从岸边点A游向正对岸点O 设
鸭子的游速为b(b>a) 且鸭子游动方向始终朝着点O 已知OAh 求鸭子游过的迹线的方程
解 取O为坐标原点 河岸朝顺水方向为x轴 y 轴指向对岸 设在时刻t鸭子位于点P(x, y)
则鸭子运动速度
Word文档

`

v(v
x
, v
y
)(
dx
,
dy
)
 故有
dx

v
x

dyv
y
dtdt
x
,
y
)

v(a
bx
, 
by
)

x
2
y
2
x
2
y
2
x
2
y2
x
2
y
2
另一方面
vab(a, 0) b(
因此
dx

v
x

a
(
x
)
2
1
x
dx

a
(
x
)
2
1
x
  即
dyv
y
by y
dybyy
dx

a
(
x
)
2
1
x

dybyy
问题归结为解齐次方程
令

y
x
u
 即xyu 得
y
du

a
u
2
1

dyb
du

a
dy
 分离变量 得
by
u
2
1
两边积分 得
arshu(lnylnC)

b
a
x
1
[(Cy)
1
b
(Cy)
1
b
]
 将
u
代入上式并整理 得
x
y
2C
以x|
yh
0代入上式 得
C
aa
1
 故鸭子游过的轨迹方程为
h
aa
y
1
b
y
1
b
h
()]
 0yh
x[()
2hh
b
将
ux
代入
arshu(lnylnC)
后的整理过程
y
a
arsh
x

b
(lnylnC)

ya



x
shln(Cy)
a

x
1
[(Cy)
a
(Cy)
a
]

y
y2
b
bb
b
bb

b
y
1 1
a
x[(Cy)(Cy)
a
]
x
1
[(Cy)
a
(Cy)
a
]

2C
2
Word文档

`
§12.4 线性微分方程
一、 线性方程
线性方程
方程
dy
P(x)yQ(x)
叫做一阶线性微分方程 dx
dydy
P(x)y0
叫做对应于非齐次线性方程
P(x)y Q(x)
的齐次线性方程
dxdx
dy
dy
y


1
y0
是齐次线性方程
dx
dxx2
如果Q(x)0  则方程称为齐次线性方程 否则方程称为非齐次线性方程
方程
下列方程各是什么类型方程？
(1)
(x2)
(2) 3x
2
5x5y0y3x
2
5x  是非齐次线性方程
(3) yy cos xe
sin
x

 是非齐次线性方程
(4)
dy
10
xy
 不是线性方程
dx
2
3
(y1)
2
dy
3< br>dy
dx
x
0
或

(5)
(y1)
 不是线性方程
x0


dy
dx
(y1)
2
dx
x
3
齐次线性方程的解法
齐次线性方程

dy
P(x)y0
是变量可分离方程 分离变量后得
dx
dy
P(x)dx

y
两边积分 得

ln|y|P(x)dxC
1

P(x)dx
或
yCe

(Ce
C
1
)


这就是齐次线性方程的通解（积分中不再加任意常数）
例1 求方程
(x2)
dy
y
的通解
dx
解 这是齐次线性方程 分离变量得

Word文档
dy
dx


yx2

`
两边积分得
ln|y|ln|x2|lnC
方程的通解为
yC(x2)
非齐次线性方程的解法
将齐次线性方程通解中的常数换成x的未知函数u(x) 把
P(x)dx

yu(x)e


设想成非齐次线性方程的通解 代入非齐次线性方程求得
P(x)dxP(x)dxP(x)dx

u

(x)e

u(x)e

P(x)P(x)u(x)e
Q(x)

化简得
u

(x)Q(x)e

P(x)dx


u(x)Q(x)e


P(x)dx
dxC

于是非齐次线性方程的通解为
P(x)dxP(x)dx
[Q(x)e

dxC]

ye


P(x)dxP(x)dxP(x)dx
e

Q(x)e

dx
 或
yCe
< br>
非齐次线性方程的通解等于对应的齐次线性方程通解与非齐次线性方程的一个特解之和
5
dy2y
(x1)
2
的通解 例2 求方程
dxx1
解 这是一个非齐次线性方程
先求对应的齐次线性方程
分离变量得

dy2y
0
的通解
dxx1
dy
2dx


yx1
两边积分得
ln y2ln (x1)ln C
齐次线性方程的通解为
yC(x1)
2

Word文档

`
用常数变易法 把C换成u 即令yu(x1)
2
 代入所给非齐次线性方程 得
2
u(x1)
2
(x1)
2
< br>u

(x1)2u(x1)
x1
2
5

1
u

(x1)
2

两边积分 得
2

u(x1)
2
C

3
再把上式代入yu(x1)
2
中 即得所求方程的通解为
3
2

y(x1)[(x1)
2
C]

3
2
3
2

Q(x)(x1)
2
 解 这里
P(x)
x1
2
)dx2ln(x1)
 因为

P(x)dx

(
x1
P(x)dx

e

e
2ln(x1)
(x1)
2
< br>5

Q(x)e
所以通解为

ye


P(x)dx
dx

5
(x1)
2
(x1)
2
dx

13
2
(x 1)
2
dx(x1)
2
3



P (x)dxP(x)dx
[

Q(x)e

dxC](x1)
2
[
2
(x1)
2
C]

3
3
例3 有一个电路如图所示 其中电源电动势为EE
msin

t(E
m
、

都是常数) 电阻R和电感
L都是常量 求电流i(t)
解 由电学知道 当电流变化时 L上有感应电动势
L

EL
即
di
 由回路电压定律得出
dt
di
iR0

dt
di

R
i
E

dtLLdi

R
i
E
m
sin

t

dtLL
把EE
m
sin

t代入上式 得

初始条件为
Word文档

`
i|
t0
0
di

R
i
E
m
sin

t
为非齐次线性方程 其中
dtLL
E
R

P(t)

Q(t)
m
sin

t

L
L
方程
由通解公式 得

i(t)e


P(t)dt
[

Q(t)e< br>
P(t)dt
dtC]


R
dt
e
L
(
R
dt
E
m


L
sin

te
L
dtC)

R
tt
E
m

R
L

e(

sin

te
L
dtC)

L

R
t
E
m


2
(Rsin

t

Lcos

t)Ce
L

22
R

L
其中C为任意常数
将初始条件i|
t0
0代入通解 得
C
因此 所求函数i(t)为
t

LE
m

R
E
m
L


i(t)
2
e(Rsin

t

Lcos

t)

R

2
L
2< br>R
2


2
L
2

LE
m

R
2


2
L
2
二、伯努利方程
伯努利方程 方程

dy
P(x)yQ(x)y
n

(n0 1)
dx
叫做伯努利方程
下列方程是什么类型方程？
dy
11
y(12x)y
4
 是伯努利方程
dx33
dy
dy
(2)
yxy
5
 
yxy
5
 是伯努利方程
dx
dx
x
y
1
(3)
y


 
y

yxy
1
 是伯努利方程
yx
x
(1)
(4)
dy
2xy4x
 是线性方程 不是伯努利方程
dx
伯努利方程的解法 以y
n
除方程的两边 得
Word文档

`

y
n
dy
P(x)y
1n
Q(x)

dx
令z y
1
n

 得线性方程
dz
(1n)P(x)z(1n)Q(x)

dx
dyy
例4 求方程
a(lnx)y
2
的通解
dxx

解 以y
2
除方程的两端 得

y
2
dy
1
1
yalnx

dxx
d(y
1
)
1
1
即
yalnx

dxx
令zy
1
 则上述方程成为

dz

1
zalnx

dxx
a
2
这是一个线性方程 它的通解为

zx[C(lnx)
2
]

以y
1
代z  得所求方程的通解为

yx[C(lnx)
2
]1

经过变量代换 某些方程可以化为变量可分离的方程 或化为已知其求解方法的方程
例5 解方程
a
2
dy

1

dxxy
解 若把所给方程变形为

dx
xy

dy
即为一阶线性方程 则按一阶线性方程的解法可求得通解 但这里用变量代换来解所给方程
令xyu 则原方程化为

du
1
1
 即
du

u1

dxudxu
u
dudx

u1
分离变量 得

Word文档

`
两端积分得
uln|u1|xln|C|
以uxy代入上式 得
yln|xy1|ln|C| 或xCe
y
y1


§12 5 全微分方程
全微分方程
一个一阶微分方程写成
P(x, y)dxQ(x, y)dy0
形式后 如果它的左端恰好是某一个函数uu(x, y)的全微分
du(x, y)P(x, y)dxQ(x, y)dy
那么方程P(x, y)dxQ(x, y)dy0就叫做全微分方程 这里

u
P(x,y)

u
Q(x,y)

y
x
而方程可写为
du(x, y)0
全微分方程的判定
若P(x, y)、Q(x, y)在单连通域G具有一阶连续偏导数 且

P

Q

yx
则方程P(x, y)dxQ(x, y)dy0是全微分方程
全微分方程的通解
若方程P(x, y)dxQ(x, y)dy0是全微分方程 且
du(x, y)P(x, y)dxQ(x, y)dy
则 u(x, y)C
即

x
x
0
P(x,y)dx

Q(x
0
,y)dxC ((x
0
,y
0
)G)

y
0
y
是方程P(x, y)dxQ(x, y)dy0的通解
例1 求解(5x
4
3xy
2
y
3
) dx(3x
2
y3xy
2
y
2
)dy0
解 这里
Word文档

`

P
6xy3y
2

Q

yx
xy
所以这是全微分方程 取(x
0
, y
0
)(0, 0) 有

u(x,y)

0
(5x
4
3xy
2
y
3
)dx< br>
y
2
dy

0

x
5
x
2
y
2
xy
3
y
3

于是 方程的通解为

x
5
x
2
y
2
xy
3
y
3
C

积分因子
若方程P(x, y)dxQ(x, y)dy0不是全微分方程 但存在一函数




(x, y) (

(x, y)0) 使方程


(x, y)P(x, y)dx

(x, y)Q(x, y)dy0
是全微分方程 则函数

(x, y)叫做方程P(x, y)dxQ(x, y)dy0的积分因子
例2 通过观察求方程的积分因子并求其通解:
(1)ydxxdy0
(2)(1xy)ydx(1xy)xdy0
解 (1)方程ydxxdy0不是全微分方程
因为
d()3
2
1
3
3
2
1
3
x
yydxxdy

2
y
所以
1
是方程ydxxdy0的积分因子 于是
y
2
ydxxdy
x
C

0
是全微分方程 所给方程的通解为
y
y
2
(2)方程 (1xy)ydx(1xy)xdy0不是全微分方程
将方程的各项重新合并 得
(ydxxdy)xy(ydxxdy)0
再把它改写成

d(xy)x
2
y
2
(
Word文档
dx

dy
)0

xy

`
这时容易看出
1
为积分因子 乘以该积分因子后 方程就变为
(xy)
2

d(xy)
dx
dy
0

(xy)
2
xy
积分得通解
1x
x

ln||lnC
 即
Ce
xy

xyy
y
我们也可用积分因子的方法来解一阶线性方程yP(x)yQ(x)
可以验证
< br>(x)e

乘以

(x)e


y

e

即
y

e

P(x)dx
1
是一阶线性方程yP(x)yQ(x)的一个积分因子 在一阶线性方程的两边
P(x)dx
得
P(x)dx
yP(x)e

y[e

P(x)dx
Q(x)e

P(x)dx

P(x)dxP(x)dxP(x)dx

]

Q(x)e

亦即
[ye

P(x)dxP(x)dx

]

Q(x)e

两边积分 便得通解

ye

P(x)dx


Q(x)e

P (x)dx
dxC

P(x)dxP(x)dx
[Q(x)e

dxC]
 或
ye


例3用积分因子求
dy
2xy4x
的通解
dx
解 方程的积分因子为


(x)e

2xdx
e
x

2
2
方程两边乘以
e
x
得

y

e
x
2xe
x
y4xe
x
 即
(e
x
y)

4xe
x

于是
e
x
y4xe
x
dx2e
x
C

2
22222

22
Word文档

`
因此原方程的通解为
y4xe
x
dx2Ce
x


22
§12 6 可降阶的高阶微分方程
一
、y
(
n
)
f (x)型的微分方程
解法 积分n 次

y
(n1)
f(x)dxC
1



y
(n2)
[f(x)dxC
1
]dxC
2


  
例1 求微分方程ye
2
x
cos x 的通解
解 对所给方程接连积分三次 得

y

e
2x
sinxC
1


y

e
2x
cosxC
1
xC
2


ye
2x
sinxC
1x
2
C
2
xC
3

这就是所给方程的通解

或
y

e
2x
sinx2C
1


y

e
2x
cosx2C
1
xC
2


ye
2x
 sinxC
1
x
2
C
2
xC
3

这就是所给方程的通解
例2 质量为m的质点受力F的作用沿Ox轴作直线运动 设力F仅是时间t的函数FF(t)
在开始时刻t0时F(0)F
0
 随着时间t的增大 此力F均匀地减小 直到tT时 F(T)0 如果开
始时质点位于原点 且初速度为零 求这质点的运动规律
解 设xx(t)表示在时刻t时质点的位置 根据牛顿第二定律 质点运动的微分方程为

m
1
2
1
41
8
1
2
1
2
1
4
1
8d
2
x
F(t)

dt
2
由题设 力F(t)随t增大而均匀地减小 且t0时 F(0)F
0
 所以F(t)F
0
kt 又当tT时 F(T)0
从而
Word文档

`

F(t)F
0
(1)

于是质点运动的微分方程又写为
t
T
d
2
x
< br>F
0
(1
t
)

T
dt
2
m
dx
|0
 其初始条件为
x|
t0
0

dt
t0

把微分方程两边积分 得
2
F
0
dxt

(t)C
1

dtm2T
再积分一次 得
F
0
1
2
t
3

x(t)C
1
tC
2

m26T
dx
|0
 由初始条件x|
t0
0
dt
t0
得C
1
C
2
0
于是所求质点的运动规律为
F
0
1
2
t
3

x(t)
 0tT
m26T
解 设xx(t)表示在时刻t时质点的位置
根据牛顿第二定律 质点运动的微分方程为
mxF(t)
由题设 F(t)是线性函数 且过点(0 F
0
)和(T 0)
故
F(t)
t
1
 即
F(t)F
0
(1
t
)

F
0
T
T
于是质点运动的微分方程又写为

x


F
0
(1
t
)

mT
其初始条件为x|
t0
0 x|
t0
0
把微分方程两边积分 得

x

2
F
0
(t
t
)C
1

m2T
再积分一次 得
F
0
1
2
t
3

x(t)C
2

m26T
由初始条件x|
t0
0 x|
t0
0
Word文档

`
得C
1
C
2
0
于是所求质点的运动规律为
F
0
1
2
t
3

x(t)
 0tT
m26T
二、y f(x y)型的微分方程
解法
设yp则方程化为
pf(x p)
设pf(x p)的通解为p

(xC
1
) 则

dy


(x,C
1
)

dx
原方程的通解为

y


(x,C
1
)dxC
2

例3 求微分方程
1xy''2xy'
满足初始条件 y|
x0
1 y|
x0
3的特解
解 所给方程是yf(x y)型的 设yp 代入方程并分离变量后 有


2

dp
2x
dx

p
1x
2
两边积分 得
ln|p|ln(1x
2
)C
即 pyC
1
(1x
2
) (C
1
e
C
)
由条件y|
x0
3 得C
1
3
所以 y3(1x
2
)
两边再积分 得 yx
3
3xC
2

又由条件y|
x0
1 得C
2
1
于是所求的特解为
yx
3
3x1
例4 设有一均匀、柔软的绳索 两端固定 绳索仅受重力的作用而下垂 试问该绳索在平衡
状态时是怎样的曲线?
三、yf(y y)型的微分方程
解法 设yp有

y


原方程化为
dpdpdydp
p

dxdydxdy
Word文档

`
dp
f(y,p)

dy
dp
f(y,p)
的通解为yp

(y C
1
) 则原方程的通解为 设方程
p
dy

p

dy


(y,C
1
)
xC
2

例5 求微分yyy
2
0的通解
解 设yp 则
y

p
代入方程 得

yp
dp

dy
dp
2
p0

dy
在y0、p0时 约去p并分离变量 得

dpdy


py
两边积分得
ln|p|ln|y|lnc
即 pCy或yCy(Cc)
再分离变量并两边积分 便得原方程的通解为
ln|y|Cxlnc
1

或 yC
1
e
Cx

(C
1
c
1
)
例6 一个离地面很高的物体受地球引力的作用由静止开始落向地面 求它落到地面时的速
度和所需的时间（不计空气阻力）
§12 7 高阶线性微分方程
一、二阶线性微分方程举例
例1 设有一个弹簧 上端固定 下端挂一个质量为m 的物体 取x 轴铅直向下 并取物体
的平衡位置为坐标原点
给物体一个初始速度v
0
0后 物体在平衡位置附近作上下振动 在振动过程中 物体的位置
x是t的函数 xx(t)
设弹簧的弹性系数为c 则恢复力fcx
又设物体在运动过程中受到的阻力的大小与速度成正比 比例系数为

 则

R

Word文档
dx

dt

`
由牛顿第二定律得
2
dxdx

m
2
cx

dt
dt
移项 并记
2n

m

k
2

c

m
2
dx
2n
dx
k
2
x0
 则上式化为
dt
dt
2
这就是在有阻尼的情况下 物体自由振动的微分方程
如果振动物体还受到铅直扰力
FHsin pt
的作用 则有
2
dx
2n
dx
k
2
xhsinpt

dt
dt
2
H
其中
h
 这就是强迫振动的微分方程
m
例2 设有一个由电阻R、自感L、电容C和电源E串联组成的电路 其中R、L、及C为
常数 电源电动势是时间t的函数 EE
m
sin

t 这里E
m
及

也是常数
设电路中的电流为i(t) 电容器极板上的电量为q(t) 两极板间的电压为u
c
 自感电动势为E
L
 由
电学知道

i
q
dq
di

u
c


E
L
L

C
dt
dt
根据回路电压定律 得
di

q
Ri0

dtC
d
2
u
c
du
c
即
LCRCu
c
E
m
sin

t

dt
dt
2

EL
或写成
d
2
u
c
du
c
E
m
2

2



usin

t

0c
dtLC
dt
2
R



1
 这就是串联电路的振荡方程 其中


0
2L
LC
如果电容器经充电后撤去外电源(E0) 则上述成为
Word文档

`
d
2
u
c
du
c
2

2



0
u
c
0

2
dt
dt
二阶线性微分方程 二阶线性微分方程的一般形式为
yP(x)yQ(x)yf(x)
若方程右端f(x)0时 方程称为齐次的 否则称为非齐次的
二、线性微分方程的解的结构
先讨论二阶齐次线性方程
d
2
ydy
yP(x)yQ(x)y0 即
2
P(x)Q(x)y0

dx
dx
定理1 如果函数y
1
(x)与y
2
(x)是方程
yP(x)yQ(x)y0
的两个解 那么
yC
1
y
1
(x)C
2
y
2
(x)
也是方程的解 其中C
1
、C
2
是任意常数
齐次线性方程的这个性质表明它的解符合叠加原理
证明 [C
1
y1
C
2
y
2
]C
1
y
1
C
2
y
2

[C
1
y
1
C
2
y
2
]C
1
y
1
C
2
y
2

因为y
1
与y
2
是方程yP(x)yQ(x)y0 所以有
y
1
P(x)y
1
Q(x)y
1
0及y
2
P(x)y
2
Q(x)y
2
0
从而 [C
1
y
1
C
2
y
2
]P(x)[ C
1
y
1
C
2
y
2
]Q(x)[ C
1
y
1
C
2
y
2
]
C
1
[y
1
P(x)y
1
Q(x)y
1
]C
2
[y
2
P(x)y
2
Q(x )y
2
]000
这就证明了yC
1
y
1(x)C
2
y
2
(x)也是方程yP(x)yQ(x)y 0的解
函数的线性相关与线性无关
设y
1
(x) y
2
(x)     y
n
(x)为定义在区间I上的n个函数 如果存在n个不全为零的常数k
1
 k
2

    k
n
 使得当xI 时有恒等式
k
1
y
1
(x)k
2
y
2
(x)     k
n
y
n
(x)0
成立 那么称这n个函数在区间I上线性相关 否则称为线性无关
判别两个函数线性相关性的方法
对于两个函数 它们线性相关与否 只要看它们的比是否为常数 如果比为常数 那么它们
就线性相关 否则就线性无关
Word文档

`
例如 1 cos
2
x  sin
2
x 在整个数轴上是线性相关的 函数1 x x
2
在任何区间(a, b)是线性无关
的
定理2 如果如果函数y
1
(x)与y
2
(x)是方程
yP(x)yQ(x)y0
的两个线性无关的解 那么
yC
1
y
1
(x)C
2
y
2
(x) (C
1
、C
2
是任意常数)
是方程的通解
例3 验证y
1
cos x与y
2
sin x是方程yy0的线性无关解 并写出其通解
解 因为
y
1
y
1
cos xcos x0
y
2
y
2
sin xsin x0
所以y
1
cos x与y
2
sin x都是方程的解
因为对于任意两个常数k
1
、k
2
 要使
k
1
cos xk
2
sin x0
只有k
1
k
2
0 所以cos x与sin x在(, )是线性无关的
因此y
1
cos x与y
2
sin x是方程yy0的线性无关解
方程的通解为yC
1
cos xC
2
sin x
例4 验证y
1
x与y
2
e
x
是方程(x1)y xyy0的线性无关解 并写出其通解
解 因为
(x1)y
1
xy
1
y
1
0xx0
(x1)y
2
xy
2
y
2
(x1)e
x
xe
x
e
x
0
所以y
1
x与y
2
e
x
都是方程的解
因为比值e

x
x 不恒为常数 所以y
1
x与y
2
e
x
在(, )是线性无关的
因此y
1
x 与y
2
e
x
是方程(x1)yxyy0的线性无关解
方程的通解为yC
1
xC
2
e

x

推论 如果y
1
(x) y
2
(x)    y
n
(x)是方程
y
(
n
)
a
1
(x)y
(
n
 1)
    a
n1
(x)y a
n
(x)y0
的n个线性无关的解 那么 此方程的通解为
yC
1
y
1
(x)C
2
y
2
(x)     C
n
y
n
(x)
Word文档

`
其中C
1
 C
2
    C
n
为任意常数
二阶非齐次线性方程解的结构
我们把方程
yP(x)yQ(x)y0
叫做与非齐次方程
yP(x)yQ(x)yf(x)
对应的齐次方程
定理3 设y*(x)是二阶非齐次线性方程
yP(x)yQ(x)yf(x)
的一个特解 Y(x)是对应的齐次方程的通解 那么
yY(x)y*(x)
是二阶非齐次线性微分方程的通解
证明提示 [Y(x)y*(x)]P(x)[ Y(x)y*(x)]Q(x)[ Y(x)y*(x)]
 [Y P(x)Y Q(x)Y ][ y* P(x)y* Q(x)y*]
0 f(x) f(x)
例如 YC
1
cos xC
2
sin x 是齐次方程yy0的通解 y*x
2
2是yyx
2
的一个特解 因此
yC
1
cos xC
2
sin xx
2
2
是方程yyx
2
的通解
定理4 设非齐次线性微分方程 yP(x)yQ(x)yf(x)的右端f(x)几个函数之和
如 yP(x)yQ(x)yf
1
(x) f
2
(x)
而y
1
*(x)与y
2
*(x)分别是方程
yP(x)yQ(x)yf
1
(x)与yP(x)yQ(x)yf< br>2
(x)
的特解 那么y
1
*(x)y
2
*(x)就是原方程的特解
证明提示
[y
1
y
2
*]P(x)[ y
1
*y
2
*]Q(x)[ y
1
*y
2
*]
[ y
1
*P(x) y
1
*Q(x) y
1
*][ y
2
*P(x) y
2
*Q(x) y
2
*]
f
1
(x)f
2
(x)
§12 8 二阶常系数齐次线性微分方程
二阶常系数齐次线性微分方程
Word文档

`
方程 ypyqy0
称为二阶常系数齐次线性微分方程 其中p、q均为常数
如果y
1
、y
2
是二阶常系数齐次线性微分方程的两个线性无关解 那么y C
1
y
1
C
2
y
2
就是它的
通解
我们看看 能否适当选取r 使ye
rx

满足二阶常系数齐次线性微分方程 为此将ye
rx
代入
方程
ypyqy0
得
(r
2
prq)e
rx

0
由此可见 只要r满足代数方程r
2
prq0 函数ye
rx
就是微分方程的解
特征方程 方程r
2
prq0叫做微分方程ypyqy0的特征方程 特征方程的两 个根r
1
、r
2
pp
2
4q
可用公式
r
求出
1,2

2
特征方程的根与通解的关系
(1)特征方程有两个不相等的实根r
1
、r
2
时
函 数
y
1
e
r
1
x
、
y
2
e
r
2
x
是方程的两个线性无关的解
这是因为
函数
y
1
e
r
1
x
、
y
2
e
r
2
x
y
1
e
r
1
x
(r
1
r
2
)x
e
是方程的 解 又不是常数
y
2
e
r
2
x
因此方程的通解为

yC
1
e
r
1
x
C
2
er
2
x

(2)特征方程有两个相等的实根r
1
r
2
时
函数
y
1
e
r
1
x
、
y
2
xe
r
1
x
是二阶常系数齐次线性微分方程的两个线性无关的解
这是因为
y
1
e
r
1
x
是方程的解 又
r
1
xr
1
x
2
r
1
x

(xe
r
1
x
)

p(xe
r
1
x
)

q(xe
r
1
x
)(2r
1
xr
1
)ep(1xr
1
)eqxe

Word文档

`
r
1
x
2

e
r
1
x
(2r
1
p)xe(r
1
pr
1
q)0

y
2
xe
r
1
x
x
不是常数 所以
y
2
xe
也是方程的解 且
y
1
e
r
1
x
r
1
x
因此方程的通解为

yC
1
e
r
1x
C
2
xe
r
1
x

(3)特征方程有一对共轭复根r
1, 2


i

时
函数ye
(
< br>
i

)
x
、ye
(


i

)
x
是微分方程的两个线性无关的复数形式的解
函数y e

x
cos

x、ye

x
sin

x是微分方程的两个线性无关的实数形式的解
函数y
1e
(


i

)
x
和y
2
e
(


i

)
x
都是方程的 解 而由欧拉公式 得
y
1
e
(

i

)
x
e

x
(cos

xisin

x)
y
2
e
(


i

)
x
e

x
(cos

xisin

x)
1
y
1
y
2
2e

x
cos

x
e< br>
x
cos

x(y
1
y
2
)

2
1
y
1
y
2
2ie< br>
x
sin

x
e

x
sin

x(y
1
y
2
)

2i
故e

x
cos

x、y
2
e
x
sin

x也是方程解
可以验证 y
1
e

x
cos

x、y
2
e

x
sin

x是方程的线性无关解
因此方程的通解为
ye

x
(C
1
cos

xC
2sin

x )
求二阶常系数齐次线性微分方程ypyqy0的通解的步骤为
第一步 写出微分方程的特征方程
r
2
prq0
第二步 求出特征方程的两个根r
1
、r
2

第三步 根据特征方程的两个根的不同情况 写出微分方程的通解
例1 求微分方程y2y3y0的通解
解 所给微分方程的特征方程为
r
2
2r30 即(r1)(r3)0
其根r
1
1 r
2
3是两个不相等的实根 因此所求通解为
Word文档

`
yC
1
e

x
C
2e
3
x

例2 求方程y2yy0满足初始条件y|
x0
4、y|
x0
2的特解
解 所给方程的特征方程为
r
2
2r10 即(r1)
2
0
其根r
1
r
2
1是两个相等的实根 因此所给微分方程的通解为
y(C
1
C
2
x)e

x

将条件y|
x0
4代入通解 得C
1
4 从而
y(4C
2
x)e

x

将上式对x求导 得
y(C
2
4C
2
x)e

x

再把条件y|
x0
2代入上式 得C
2
2 于是所求特解为
x(42x)e

x

例 3 求微分方程y2y5y 0的通解
解 所给方程的特征方程为
r
2
2r50
特征方程的根为r
1
12i r
2
12i 是一对共轭复根
因此所求通解为
ye
x
(C
1
cos2xC
2
sin2x)
n 阶常系数齐次线性微分方程 方程
y
(
n
)
p
1
y
(
n
1)
p
2
y
(
n
2)
     p
n1
yp
n
y0
称为n 阶常系数齐次线性微分方程 其中 p
1
 p
2
     p
n1
 p
n
都是常数
二阶常系数齐次线性微分方程所用的方法以及方程的通解形式 可推广到n 阶常系数齐次
线性微分方程上去
引入微分算子D 及微分算子的n次多项式
L(D)=D
n

p
1
D
n
1
p
2
D
n
2
     p
n1
Dp
n

则n阶常系数齐次线性微分方程可记作
(D
n

p
1
D
n
1
p
2
D
n
2
     p
n1
Dp
n
)y0或L(D)y0
注 D叫做微分算子D
0
yy Dyy D
2
yy D
3
yy   D
n
yy
(
n
)

分析 令ye
rx
 则
Word文档

`
L(D)yL(D)e
rx
(r
n

p
1
r
n
1
p
2
r
n
2
     p
n1
rp
n
)e
rx
L(r)e
rx

因此如果r是多项式L(r)的根 则ye
rx
是微分方程L(D)y0的解
n 阶常系数齐次线性微分方程的特征方程
L(r)r
n

p
1
r
n
1
p
2
r
n
2
     p
n1
rp
n
0
称为微分方程L(D)y0的特征方程
特征方程的根与通解中项的对应
单实根r 对应于一项 Ce
rx


一对单复根r
1

2


i

对应于两项 e

x
(C
1
cos

xC
2
sin

x)
k重实根r对应于k项 e
rx
(C
1
C
2
x    C
k
x
k
1
)

一
对k 重复根r
1

2


i

对应于2k项
e

x
[(C
1
C
2
x    C
k
x
k
1
)cos

x( D
1
D
2
x    D
k
x
k
1
)sin

x]
例4 求方程y
(4)
2y5y0 的通解
解 这里的特征方程为
r
4
2r
3
5r
2
0 即r
2
(r
2
2r5)0
它的根是r
1
r
2
0和r
3

4
12i
因此所给微分方程的通解为
yC< br>1
C
2
xe
x
(C
3
cos2xC< br>4
sin2x)
例5 求方程y
(4)



4
y0的通解 其中

0
解 这里的特征方程为
r
4


4
0
它的根为
r
1,2


2
(1i)

r
3,4


2
(1i)

因此所给微分方程的通解为

ye
2
x
(C
1
cos

2
xC
2
sin

2
x)e


2
x
(C
3
cos
2
xC
4
sin

2
x)

§12 9 二阶常系数非齐次线性微分方程
二阶常系数非齐次线性微分方程
方程 ypyqyf(x)
称为二阶常系数非齐次线性微分方程 其中p、q是常数
Word文档

`
二阶常系数非齐次线性微分方程的通解是对应的齐次方程
的通解yY(x)与非齐次方程本身的一个特解yy*(x)之和
yY(x) y*(x)
当f(x)为两种特殊形式时 方程的特解的求法
一、 f(x)P
m
(x)e

x

型
当f(x)P
m
(x)e

x
时 可以猜想 方程的特解也应具有这种形式 因此 设特解形式为
y*Q(x)e

x
 将其代入方程 得等式
Q(x)(2

p)Q(x)(

2
p
< br>q)Q(x)P
m
(x)
(1)如果

不是特征方程r
2
prq0 的根 则

2
p

q0
要使上式成立 Q(x)应设为m 次多项式
Q
m
(x)b
0x
m
b
1
x
m
1
    b
m1
xb
m

通过比较等式两边同次项系数 可确定b
0
 b
1
     b
m
 并得所求特解
y*Q
m
(x)e

x

(2)如果

是特征方程 r
2
prq0 的单根 则

2
p

q0 但2

p0 
要使等式
Q(x)(2

p)Q(x)(< br>
2
p

q)Q(x)P
m
(x)
成立 Q(x)应设为（m1） 次多项式
Q(x)xQ
m
(x)
Q
m
(x)b
0
x
m

b
1
x
m
1
    b
m1
xb
m

通过比较等式两边同次项系数 可确定b
0
 b
1
     b
m
 并得所求特解
y*xQ
m
(x)e

x

(3)如果

是特征方程 r
2
prq0的二重根 则

2
p

q02

p0。
要使等式
Q(x)(2

p)Q(x)(< br>
2
p

q)Q(x)P
m
(x)
成立 Q(x)应设为（m2）次多项式
Q(x)x
2
Q
m
(x)
Q
m< br>(x)b
0
x
m
b
1
x
m
1
    b
m1
xb
m

通过比较等式两边同次项系数 可确定b
0
 b
1
     b
m
 并得所求特解
y*x
2
Q
m
(x)e

x

Word文档

`
综上所述 我们有如下结论
如果f(x)P
m
(x)e

x
 则二阶常系数非齐次线性微分方程ypyqy f(x)有形如
y*x
k

Q
m
(x)e

x

的特解 其中Q
m
(x)是与P
m
(x)同次的多项式 而k 按

不是特征方程的根、是特征方程的单根
或是特征方程的的重根依次取为0、1或2 
例1 求微分方程y2y3y3x1的一个特解
解 这是二阶常系数非齐次线性微分方程 且函数f(x)是P
m
(x)e

x
型(其中P
m
(x)3x1

0)
与所给方程对应的齐次方程为
y2y3y0
它的特征方程为
r
2
2r30
由于这里

0不是特征方程的根 所以应设特解为
y*b
0
xb
1

把它代入所给方程 得
3b
0
x2b
0
3b
1
3x1
比较两端x同次幂的系数 得
0


 3b
0
3 2b
0
3b
1
1
 2b3b1
1

0

3b3
由此求得b
0
1
b
1

 于是求得所给方程的一个特解为

y*x

例2 求微分方程y5y6yxe
2
x
的通解
解 所给方程是二阶常系数非齐次线性微分方程 且f(x)是P
m
(x)e

x
型(其中P
m
(x)x

2)
与所给方程对应的齐次方程为
y5y6y0
它的特征方程为
r
2
5r

60
特征方程有两个实根r
1
2 r
2
3 于是所给方程对应的齐次方程的通解为
YC
1
e
2
x
C
2
e
3
x


Word文档
1
3
1
3

`
由于

2是特征方程的单根 所以应设方程的特解为
y*x(b
0
xb
1
)e
2
x

把它代入所给方程 得
2b
0
x2b
0
b
1
x
比较两端x同次幂的系数 得
0


 2b
0
1 2b
0
b
1
0
2b b0

01

2b1
由此求得
b
0

1
 b1 于是求得所给方程的一个特解为
1
2

y*x(x1)e
2x

从而所给方程的通解为

yC
1
e
2x
C
2
e
3x
 (x
2
2x)e
2x

提示
y*x(b
0
xb
1
)e
2
x
(b
0
x
2
b
1
x)e
2
x

[(b
0< br>x
2
b
1
x)e
2
x
][(2b0
xb
1
)(b
0
x
2
b
1< br>x)2]e
2
x

[(b
0
x
2b
1
x)e
2
x
][2b
0
2(2 b
0
xb
1
)2(b
0
x
2
b< br>1
x)2
2
]e
2
x


y *5y*6y*[(b
0
x
2
b
1
x)e< br>2
x
]5[(b
0
x
2
b
1
x)e
2
x
]6[(b
0
x
2
b
1
x)e
2
x
]
[2b
0
2(2b
0
xb
1
)2(b
0
x
2
b
1< br>x)2
2
]e
2
x
5[(2b
0
xb
1
)(b
0
x
2
b
1
x)2]e< br>2
x
6(b
0
x
2
b
1
x)e
2
x

[2b
0
4(2b
0
xb< br>1
)5(2b
0
xb
1
)]e
2
x[2b
0
x2b
0
b
1
]e
2
x

二、ypyqye

x
[P
l
(x)cos

xP
n
(x)sin

x]
特解形式应用欧拉公式可得
e

x
[P
l
(x) cos

xP
n
(x)sin

x]
i

xi

xi

x
eee

P
n
(x)]

l
22i
1
(x)iP(x)]e
(

i

) x

1
[P(x)iP(x)]e
(

i
< br>)x

[P
nn
2
l
2
l1
2
1
2
e

x
[P(x)
ei

x

P(x)e
(

i
)x
P(x)e
(

i

)x

Word文档

`
其中
P(x)(P
l
P
n
i)

P(x)(P
l
P
n
i)
 而mmax{l n}
设方程ypyqyP(x)e
(


i

)
x
的特解为y
1
*x
k
Qm
(x)e
(


i

)
x

则
y
1
*x
k
Q
m
(x)e(

i

)
必是方程
y

py

qyP(x)e
(

i

)
的特 解
其中k按

i

不是特征方程的根或是特征方程的根依次 取0或1
于是方程ypyqye

x
[P
l
(x)cos

xP
n
(x)sin

x] 的特解为

y*x
k
Q
m
(x)e(

i

)x
x
k
Q
m
(x)e
(

i

)x


x
k
e

x
[Q
m
(x)(cos
xisin

x)Q
m
(x)(cos

xisin

x)

x
k
< br>e

x
[R
(1)
m
(x)cos

xR
(2)
m
(x)sin

x]
综上所述 我们有如下结论
如果f(x)e

x

[P
l
(x)cos

xP
n
(x)sin

x] 则二阶常系数非齐次线性微分方程
ypyqyf(x)
的特解可设为
y*x
k

e

x
[R
(1)
m
(x)cos

xR
(2)
m
(x)sin

x]
其中R
(1)
m
(x)、R
(2)
m
( x)是m次多项式 mmax{l n} 而k 按

i

(或

i

)不是特征方程的根或
是特征方程的单根依次取0或1
例3 求微分方程yyxcos2x的一个特解
解 所给方程是二阶常系数非齐次线性微分方程
且f(x)属于e

x
[P
l
(x)cos

xP
n
(x)sin

x]型(其中

0

2 P
l
(x)x P
n
(x)0）
与所给方程对应的齐次方程为
yy0
它的特征方程为
r
2
10
由于这里

i

2i 不是特征方程的根 所以应设特解为
y*(axb)cos2x(cxd )sin2x
把它代入所给方程 得
(3ax3b4c)cos2x(3cx3d4a)sin2xxcos2x
Word文档
1
2
1
2

`
比较两端同类项的系数 得
a
 b0 c0
d
于是求得一个特解为
y*xcos2xsin2x

提示
y*(axb)cos2x(cxd)sin2x
y*acos2x2(axb)sin2xcsin2x2(cxd)cos2x
(2cxa2d)cos2x(2ax2bc)sin2x
1
3
4

9
1
3
4
9
y*2ccos2x2(2cxa2d)sin2x2asin2x2(2ax2bc )cos2x
(4ax4b4c)cos2x(4cx4a4d)sin2x
y* y*(3ax3b4c)cos2x(3cx4a3d)sin2x


3a1

3b4c0
14
由

 得
a
 b0 c0
d

3c0
39


4a3d0
§12 10 微分方程的幂级数解法
当微分方程的解不能用初等函数或其积分表达时 我们就要寻求其它解法 常用的有幂级
数解法和数值解法 本节我们简单地介绍微分方程的幂级数解法
求一阶微分方程
(yy
0
)的多项式
f(x y )a
00
a
10
(xx
0
)a
01
(yy
0
)    a
im
(xx
0
)
l
(yy
0
)
m

这时我们可以设所求特解可展开为xx
0
的幂级数
yy
0
a
1
(xx
0
)a
2
(x x
0
)
2
    a
n
(xx
0
)
n
    
其中a
1
 a
2
     a
n
     是待定的系数 把所设特解代入微分方程中 便得一恒等式 比较这恒
等式两端xx
0
的同次幂的系数 就可定出常数a
1
 a
2
     从而得到所求的特解
例1 求方程
dy
f(x,y)
满足初始条件
y|
xx
0
y
0
的特解 其中函数f(x y)是(xx
0
)、
dx
dy
xy
2
满足y|
x0
0的特解
dx
解 这时x
0
0 y
0
0 故设
ya
1
xa
2
x
2
a
3
x
3a
4
x
4
    
把y及y的幂级数展开式代入原方程 得
Word文档

`
a
1
2a
2
x3a
3
x
2
4a
4
x
3
5a
5
x
4< br>   
x(a
1
xa
2
x
2
a
3
x
3
a
4
x
4
    )
2

xa
1
2
x
2
2a
1
a
2
x
3
(a
2
2
2a
1
a
3
)x
4
    
由此 比较恒等式两端x的同次幂的系数 得
a
1
0
a
2

 a
3
0 a
4
0
a
5

1
2
1
    
20
于是所求解的幂级数展开式的开始几项为

yx
2

1
2
1
x
5
   

20
定理 如果方程
yP(x)yQ(x)y0
中的系数P(x)与Q(x)可在R
y

a
n
x
n

n0

的解
例2 求微分方程yxy 0的满足初始条件y|
x0
0 y|
x0
1的特解
解 这里P(x)0 Q(x)x在整个数轴上满足定理的条件 因此所求的解可在整个数轴上展开
成x的幂级数
ya
0a
1
xa
2
x
2
a
34
3xa
4
x   


a
n
x
n

n0

由条件y|
x0
0 得a
0
0 由ya
1
2a
2
x3a
3
x
2
 4a
4
x
3
   及y|
x0
1 得a
1
1 于是

yxa
2
xa
3
xa
4
x    x

a
n
x
n

n2
234


y

12a
2
x3a
3
x4a
4
x    1

na
n
x
n1

n2
23


y

2a2
32a
3
x43a
4
x
2
    

n(n1)a
n
x
n2

n2

yxa
2
x
2
a
3
x
3
a
4
x
4
   
x

a
n
x
n

n2

y12a
2
x3a
3

x
2
4a
4
x
3
   
1

na
n
x
n1

n2

y2a
2
x32a
3
x43a
4
x
2
   


n(n1)a
n
x
n2

n2

Word文档

`

把y及y代入方程yxy 0 得
2a2
32a
3
x43a
4
x
2
    n(n1)a
n
x
n
2
  
x(xa
2
x
2
a
3
x
3
a4
x
4
  a
n
x
n
  )0

即 2a
2
32a
3
x (43a
4
1)x
2
(54a
5
a
2
)x
3

(65a
6
a
3
)x
4
    [(n2)(n1)a
n2
a
n1
]x
n
    0
于是有

a
2
0, a
3
0, a
4

一般地
a
n2

1
, a0, a0,   

6
43
5
a
n1
(n3 4   )
(n2)(n1)
由递推公式可得
a
a
4
1

1
, a
8
0, a
9
0, a
10

7
,   

7676431091097643
1
一般地
a
3m1

(m1 2   )
(3m1)(3m)    7643

a
7

所求的特解为

yx


1
x
4

1
x
7

1
x
10
   

4376431097643
Word文档