第一篇：ode微分方程函数的书写格式

x0=[0;0]

f=@(t,x)[-41.86*x(1)^5 + 113.155*x(1)^40.58945*x(1)^4*x(2)+ 4.931*x(1)^43.6245*x(1)^30.41191*x(1)^2*x(2)0.27278*x(1)*x(2)^4 + 0.034887*x(1)*x(2)^3-0.91755*x(1)*x(2)^2 + 0.031226*x(1)*x(2)0.48587*x(2)^5 + 5.8169e-08*x(2)^41.9321e-07*x(2)^2-0.90217*x(2)];

opt=odeset;

opt.RelTol=1e-6;

tic,[t,x]=ode15s(f,[-1,10],x0,opt);toc

length(t)

figure(1)

plot(t,x(:,1))

第二篇：微分方程

高等数学教案

§12 微分方程

第十二章

微分方程

教学目的：

1．了解微分方程及其解、阶、通解，初始条件和特等概念。2．熟练掌握变量可分离的微分方程及一阶线性微分方程的解法。

3．会解齐次微分方程、伯努利方程和全微分方程，会用简单的变量代换解某些微分方程。4． 会用降阶法解下列微分方程：y(n)f(x)，yf(x,y)和yf(y,y)

5． 理解线性微分方程解的性质及解的结构定理。

6．掌握二阶常系数齐次线性微分方程的解法，并会解某些高于二阶的常系数齐次线性微分方程。

7.求自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程的特解和通解。

8.会解欧拉方程，会解包含两个未知函数的一阶常系数线性微分方程组。9．会解微分方程组（或方程组）解决一些简单的应用问题。教学重点：

1、可分离的微分方程及一阶线性微分方程的解法

(n)

2、可降阶的高阶微分方程yf(x)，yf(x,y)和yf(y,y)

3、二阶常系数齐次线性微分方程；

4、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程；

教学难点：

1、齐次微分方程、伯努利方程和全微分方程；

2、线性微分方程解的性质及解的结构定理；

3、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程的特解。

高等数学教案

§12 微分方程

4、欧拉方程

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

函数是客观事物的内部联系在数量方面的反映 利用函数关系又可以对客观事物的规律性进行研究 因此如何寻找出所需要的函数关系 在实践中具有重要意义 在许多问题中 往往不能直接找出所需要的函数关系 但是根据问题所提供的情况 有时可以列出含有要找的函数及其导数的关系式 这样的关系就是所谓微分方程 微分方程建立以后 对它进行研究 找出未知函数来 这就是解微分方程

例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x2C

(3)其中C是任意常数

把条件“x1时 y2”代入(3)式 得

212C

由此定出C1 把C1代入(3)式 得所求曲线方程(称为微分方程满足条件y|x12的解)

yx21

例2 列车在平直线路上以20m/s(相当于72km/h)的速度行驶 当制动时列车获得加速度04m/s2 问开始制动后多少时间列车才能停住 以及列车在这段时间里行驶了多少路程?

解 设列车在开始制动后t秒时行驶了s米 根据题意 反映制动阶段列车运动规律的函数ss(t)应满足关系式

d2s0.4

(4)dt22 此外 未知函数ss(t)还应满足下列条件

高等数学教案

§12 微分方程

t0时 s0 vds20 简记为s|=0 s|=20

(5)

t0t0dt

把(4)式两端积分一次 得

vds0.4tC

(6)1dt再积分一次 得

s02t2 C1t C2

(7)这里C1 C2都是任意常数

把条件v|t020代入(6)得

20C1

把条件s|t00代入(7)得0C2

把C1 C2的值代入(6)及(7)式得

v04t 20

(8)

s02t220t

(9)在(8)式中令v0 得到列车从开始制动到完全停住所需的时间

t2050(s)

0.4再把t50代入(9) 得到列车在制动阶段行驶的路程

s025022050500(m)

解 设列车在开始制动后t秒时行驶了s米

s04 并且s|t0=0 s|t0=20

把等式s04两端积分一次 得

s04tC1 即v04tC1(C1是任意常数)

再积分一次 得

s02t2 C1t C2(C1 C2都C1是任意常数)

由v|t020得20C1 于是v04t 20

由s|t00得0C2 于是s02t220t

令v0 得t50(s) 于是列车在制动阶段行驶的路程

s025022050500(m) 高等数学教案

§12 微分方程

几个概念

微分方程 表示未知函数、未知函数的导数与自变量之间的关系的方程 叫微分方程

常微分方程 未知函数是一元函数的微分方程 叫常微分方程

偏微分方程 未知函数是多元函数的微分方程 叫偏微分方程

微分方程的阶 微分方程中所出现的未知函数的最高阶导数的阶数 叫微分方程的阶

x3 yx2 y4xy3x2 

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x0 时 yy0  y y0 

一般写成



yxx0y0 yxx0y0

特解 确定了通解中的任意常数以后 就得到微分方程的特解 即不含任意常数的解

初值问题 求微分方程满足初始条件的解的问题称为初值问题

如求微分方程yf(x

y)满足初始条件yxx0y0的解的问题 记为

yf(x,y)

 yxx0y0

积分曲线 微分方程的解的图形是一条曲线 叫做微分方程的积分曲线

例3 验证 函数 高等数学教案

§12 微分方程

xC1cos ktC2 sin kt 是微分方程

2dxk2x0

dt2的解

解 求所给函数的导数

dxkCsinktkCcoskt 12dtd2xk2Ccosktk2Csinktk2(CcosktCsinkt) 1212dt22dx将2及x的表达式代入所给方程 得 dt

k2(C1cos ktC2sin kt) k2(C1cos ktC2sin kt)0

2dx

这表明函数xC1cosktC2sinkt 满足方程2k2x0 因此所给函数是所给方程的解

dt2dx

例4 已知函数xC1cosktC2sinkt(k0)是微分方程2k2x0的通解 求满足初始条件 dt

x| t0 A x| t0 0 的特解

解

由条件x| t0 A及xC1 cos ktC2 sin kt 得

C1A

再由条件x| t0 0 及x(t)kC1sin ktkC2cos kt 得

C20

把C1、C2的值代入xC1cos ktC2sin kt中 得

xAcos kt

§12 2 可分离变量的微分方程

观察与分析

1 求微分方程y2x的通解 为此把方程两边积分 得 yx2C

一般地 方程yf(x)的通解为yf(x)dxC(此处积分后不再加任意常数)

高等数学教案

§12 微分方程

2 求微分方程y2xy2 的通解

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

为求通解可将方程变为

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

y21x2C 或y21

yxC可以验证函数y1是原方程的通解

x2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)若把y看作自变量、x看作未知函数 则当P(x,y)0时 有

dxQ(x,y)

dyP(x,y)

可分离变量的微分方程

如果一个一阶微分方程能写成

g(y)dyf(x)dx(或写成y(x)(y))的形式 就是说 能把微分方程写成一端只含y的函数和dy 另一端只含x的函数和dx 那么原方程就称为可分离变量的微分方程

讨论 下列方程中哪些是可分离变量的微分方程？(1)y2xy

是 y1dy2xdx  高等数学教案

§12 微分方程

(2)3x25xy0

是 dy(3x25x)dx(3)(x2y2)dxxydy=0

不是

(4)y1xy2xy2 是 y(1x)(1y2)(5)y10xy

是 10ydy10xdx(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

解

此方程为可分离变量方程 分离变量后得

1dy2xdx

y1dy2xdx

y两边积分得

即

ln|y|x2C1

从而

yex2C1eC1ex 2因为eC1仍是任意常数 把它记作C 便得所给方程的通解

yCex

解

此方程为可分离变量方程 分离变量后得

21dy2xdx

y1dy2xdx

y7 两边积分得

即

ln|y|x2lnC

高等数学教案

§12 微分方程

从而

yCex

例2 铀的衰变速度与当时未衰变的原子的含量M成正比 已知t0时铀的含量为M0 求在衰变过程中铀含量M(t)随时间t变化的规律

解 铀的衰变速度就是M(t)对时间t的导数2dM

dt

由于铀的衰变速度与其含量成正比 故得微分方程

dMM

dtdM0

dt其中(>0)是常数 前的曲面号表示当t增加时M单调减少 即由题意 初始条件为

M|t0M0

将方程分离变量得

dMdt

M两边积分 得M()dt dM即

lnMtlnC 也即MCet

由初始条件 得M0Ce0C

所以铀含量M(t)随时间t变化的规律MM0et 

例3 设降落伞从跳伞塔下落后 所受空气阻力与速度成正比 并设降落伞离开跳伞塔时速度为零 求降落伞下落速度与时间的函数关系

解

设降落伞下落速度为v(t) 降落伞所受外力为Fmgkv(k为比例系数) 根据牛顿第二运动定律Fma 得函数v(t)应满足的方程为

mdvmgkv

dt初始条件为

v|t00

方程分离变量 得

高等数学教案

§12 微分方程

dvdt

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

mgkvm

ln(mgkv)1ktC

m1kC1ktmgem即

v()

CeCkkmg将初始条件v|t00代入通解得C

kktmg于是降落伞下落速度与时间的函数关系为v(1em)

kdy

例4 求微分方程1xy2xy2的通解

dx

解 方程可化为

dy(1x)(1y2)

dx1dy(1x)dx

1y2分离变量得

两边积分得

1dy(1x)dx 即arctany1x2xC

1y22于是原方程的通解为ytan(x2xC)

例4 有高为1m的半球形容器 水从它的底部小孔流出 小孔横截面面积为1cm2 开始时容器内盛满了水 求水从小孔流出过程中容器里水面高度h随时间t变化的规律

解 由水力学知道 水从孔口流出的流量Q可用下列公式计算

Q12dV0.62S2gh

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

dV0.622gh 或dV0.622ghdt

dt

另一方面 设在微小时间间隔[t tdt]内 水面高度由h降至hdh(dh0) 则又可得到 高等数学教案

§12 微分方程

dVr2dh

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

r1002(100h)2200hh2

所以

dV(200hh2)dh

通过比较得到

0.622ghdt(200hh2)dh

这就是未知函数hh(t)应满足的微分方程

此外 开始时容器内的水是满的 所以未知函数hh(t)还应满足下列初始条件

h|t0100

将方程0.622ghdt(200hh2)dh分离变量后得

dt两端积分 得

t0.622g132(200hh2)dh

0.622g132(200hh2)dh

(400h22h2)C

即

t50.622g3其中C是任意常数

由初始条件得 35(400100221002)C

t50.622gC35(400000200000)14105

350.622g0.622g150.622g(7105353210h3h2) 因此

t 高等数学教案

§12 微分方程

上式表达了水从小孔流出的过程中容器内水面高度h与时间t之间的函数关系

§12 3 齐次方程

齐次方程

如果一阶微分方程dyf(x,y)中的函数f(x, y)可写成 dxyy的函数 即f(x,y)() 则称这方程为齐次方程

xx

下列方程哪些是齐次方程？

dyyy2x2dyyy()21

(1)xyyyx0是齐次方程dxxdxxx222dy1y

(2)1xy1y不是齐次方程

dx1x222dyx2y2dyxy

(3)(xy)dxxydy0是齐次方程 dxxydxyx2

2(4)(2xy4)dx(xy1)dy0不是齐次方程

(5)(2xsh3ych)dx3xchdy2xy4

dxxy1yxyxydy0是齐次方程

xyy2xsh3ychdyxxdy2thyy 

ydxdx3xx3xchx

齐次方程的解法

ydyy()中 令u 即yux 有 dxxxdu(u)

uxdx

在齐次方程分离变量 得 高等数学教案

§12 微分方程

dudx (u)uxdudx

(u)ux两端积分 得

y代替u 便得所给齐次方程的通解

xdydy

例

1解方程y2x2xy

dxdx求出积分后 再用

解

原方程可写成

y2()dyyx



dxxyx2y1x2因此原方程是齐次方程 令

yux 于是原方程变为

2duu

ux

dxu1duu 即

xdxu1yu 则 xdyuxdu

dxdx分离变量 得

(1)du1udx

x两边积分 得uln|u|Cln|x|

或写成ln|xu|uC

以y代上式中的u 便得所给方程的通解 x

ln|y|yC

x

例2 有旋转曲面形状的凹镜 假设由旋转轴上一点O发出的一切光线经此凹镜反射后都与旋转轴平行 求这旋转曲面的方程 高等数学教案

§12 微分方程

解 设此凹镜是由xOy面上曲线L yy(x)(y>0)绕x轴旋转而成 光源在原点 在L上任取一点M(x, y) 作L的切线交x轴于A 点O发出的光线经点M反射后是一条平行于x轴射线 由光学及几何原理可以证明OAOM

因为

OAAPOPPMcotOP而

OMx2y2

于是得微分方程

yx

yyxx2y2 y整理得dxx(x)21 这是齐次方程

dyyydxx(x)21

dyyy

问题归结为解齐次方程

令即

yxv 即xyv 得vydvvv21

dyydvv21

dydvdy

v21y分离变量 得两边积分 得 ln(vv21)lnylnC, vv21yy, (v)2v21, CCy22yv1

C2C以yvx代入上式 得y22C(xC)

2这是以x轴为轴、焦点在原点的抛物线 它绕x轴旋转所得旋转曲面的方程为

y2z22C(xC) 2这就是所求的旋转曲面方程

例3 设河边点O的正对岸为点A 河宽OAh 两岸为平行直线 水流速度为a 有一鸭子从点A游向点O 设鸭子的游速为b(b>a) 且鸭子游动方向始终朝着点 O 求鸭子游过的迹线的方程 高等数学教案

§12 微分方程

例3 设一条河的两岸为平行直线 水流速度为a 有一鸭子从岸边点A游向正对岸点O 设鸭子的游速为b(b>a) 且鸭子游动方向始终朝着点O 已知OAh 求鸭子游过的迹线的方程

解 取O为坐标原点 河岸朝顺水方向为x轴 y 轴指向对岸 设在时刻t鸭子位于点P(x, y) 则鸭子运动速度

v(vx, vy)(dx, dy) 故有dxvx

dyvydtdtx, y) v(abx, by)

x2y2x2y2x2y2x2y2另一方面 vab(a, 0)b(因此dxvxa(x)21x 即dxa(x)21x

dyvybyydybyy

问题归结为解齐次方程dxa(x)21x

dybyy

令

yxu 即xyu 得

yduau21

dybduady 分离变量 得u21by两边积分 得 arshu(lnylnC) bax1[(Cy)1b(Cy)1b]

将u代入上式并整理 得xy2C以x|yh0代入上式 得Caa1 故鸭子游过的轨迹方程为

haay1by1bh

x[()()] 0yh

2hhb将ux代入arshu(lnylnC)后的整理过程

yaarshxb(lnylnC)ya 高等数学教案

§12 微分方程

xshln(Cy)ax1[(Cy)a(Cy)a] yy2bbbbbbby11ax[(Cy)(Cy)a]x1[(Cy)a(Cy)a]

2C2

§12.4 线性微分方程

一、线性方程

线性方程

方程dyP(x)yQ(x)叫做一阶线性微分方程 dxdydyP(x)y0叫做对应于非齐次线性方程P(x)yQ(x)的齐次线性方程

dxdxdydyy1y0是齐次线性方程 dxdxx2如果Q(x)0  则方程称为齐次线性方程 否则方程称为非齐次线性方程

方程

下列方程各是什么类型方程？

(1)(x2)

(2)3x25x5y0y3x25x  是非齐次线性方程

(3)yy cos xesin x  是非齐次线性方程

(4)dy10xy 不是线性方程 dx23(y1)2dy3dydxx0或

(5)(y1) 不是线性方程 x0dydx(y1)2dxx

3齐次线性方程的解法

齐次线性方程

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

y两边积分 得

ln|y|P(x)dxC1  15 高等数学教案

§12 微分方程

P(x)dx(CeC1)

或

yCe这就是齐次线性方程的通解（积分中不再加任意常数）

例

1求方程(x2)dyy的通解

dx

解

这是齐次线性方程 分离变量得

dydx

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P(x)P(x)u(x)eQ(x)

u(x)e化简得

u(x)Q(x)eP(x)dx

u(x)Q(x)eP(x)dxdxC

于是非齐次线性方程的通解为

P(x)dxP(x)dx

ye[Q(x)edxC] P(x)dxP(x)dxP(x)dx或

yCeeQ(x)edx 非齐次线性方程的通解等于对应的齐次线性方程通解与非齐次线性方程的一个特解之和

5dy2y

例2 求方程(x1)2的通解

dxx1

解

这是一个非齐次线性方程

高等数学教案

§12 微分方程

先求对应的齐次线性方程分离变量得

dy2y0的通解

dxx1dy2dx

yx1两边积分得

ln y2ln(x1)ln C

齐次线性方程的通解为

yC(x1)2

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

52u(x1)2(x1)2

u(x1)2u(x1)x1 1u(x1)2

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

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

y(x1)[(x1)2C]

3232 Q(x)(x1)2

解 这里P(x)x12)dx2ln(x1) 因为

P(x)dx(x1P(x)dxe2ln(x1)(x1)2

e Q(x)eP(x)dxdx5(x1)2(x1)2dx132(x1)2dx(x1)23所以通解为

yeP(x)dx[Q(x)eP(x)dxdxC](x1)[2(x1)2C]

323

例3 有一个电路如图所示 其中电源电动势为EEmsint(Em、都是常数) 电阻R和电感L都是常量 求电流i(t)高等数学教案

§12 微分方程

解

由电学知道 当电流变化时 L上有感应电动势L

EL即

di 由回路电压定律得出

dtdiiR0

dtdiRiE

dtLLdiRiEmsin t

dtLL

把EEmsin t代入上式 得

初始条件为

i|t00

diRiEmsin t为非齐次线性方程 其中 dtLLER

P(t) Q(t)msin t

LL

方程由通解公式 得

i(t)eP(t)dt[Q(t)eP(t)dtdtC]RdteL(RdtEmLsin teLdtC)

RttEmR

eL(sinteLdtC)

LRtEm

2(Rsin t Lcos t)CeL

22RL其中C为任意常数

将初始条件i|t00代入通解 得C因此 所求函数i(t)为

t LEmREmL

i(t)2e(Rsin t Lcos t)

R2L2R22L2 LEm

R22L

2二、伯努利方程

伯努利方程 方程

dyP(x)yQ(x)yn(n0 1)dx叫做伯努利方程

下列方程是什么类型方程？ 高等数学教案

§12 微分方程

dy11y(12x)y4 是伯努利方程 dx33dydy

(2)yxy5 yxy5 是伯努利方程

dxdxxy

1(3)y yyxy1 是伯努利方程 yxx

(1)

(4)dy2xy4x 是线性方程 不是伯努利方程 dxdyP(x)y1nQ(x)dx

伯努利方程的解法 以yn除方程的两边 得

yn令z y1n  得线性方程

dz(1n)P(x)z(1n)Q(x)

dxdyya(lnx)y2的通解

dxx

例4 求方程

解 以y2除方程的两端 得

dy11yalnx

dxxd(y1)11即

yalnx

dxx

y2令zy1 则上述方程成为

dz1zalnx

dxxa2这是一个线性方程 它的通解为

zx[C(lnx)2]

以y1代z  得所求方程的通解为

yx[C(lnx)2]1

经过变量代换 某些方程可以化为变量可分离的方程 或化为已知其求解方法的方程

例5 解方程a2dy1

dxxy

解

若把所给方程变形为 高等数学教案

§12 微分方程

dxxy

dy即为一阶线性方程 则按一阶线性方程的解法可求得通解 但这里用变量代换来解所给方程

令xyu 则原方程化为

du11 即duu1

dxudxu分离变量 得

ududx

u1两端积分得

uln|u1|xln|C|

以uxy代入上式 得

yln|xy1|ln|C| 或xCey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是全微分方程 且 高等数学教案

§12 微分方程

du(x, y)P(x, y)dxQ(x, y)dy 则

u(x, y)C

即

xP(x,y)dxyQ(x0,y)dxC((x0,y0)G)

00xy是方程P(x, y)dxQ(x, y)dy0的通解

例1 求解(5x43xy2y3)dx(3x2y3xy2y2)dy0

解 这里

P6xy3y2Q

yxxy所以这是全微分方程 取(x0, y0)(0, 0) 有

u(x,y)0(5x43xy2y3)dxy2dy

0

x5于是 方程的通解为

x53x2y2xy31y3

233x2y2xy31y3C

3积分因子 若方程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()xyydxxdy

2y所以1是方程ydxxdy0的积分因子 于是 y2 21 高等数学教案

§12 微分方程

ydxxdyxC 是全微分方程 所给方程的通解为0yy

2(2)方程(1xy)ydx(1xy)xdy0不是全微分方程

将方程的各项重新合并 得

(ydxxdy)xy(ydxxdy)0

再把它改写成 d(xy)x2y2(这时容易看出dxdy)0

xy1为积分因子 乘以该积分因子后 方程就变为(xy)2

d(xy)dxdy0

2xy(xy)积分得通解

x1x

ln||lnC 即Cexy

yxyy

我们也可用积分因子的方法来解一阶线性方程yP(x)yQ(x)

可以验证(x)e两边乘以(x)e

ye即

ye亦即

[yeP(x)dx1是一阶线性方程yP(x)yQ(x)的一个积分因子 在一阶线性方程的P(x)dx得

P(x)dxP(x)dxyP(x)ey[eQ(x)eP(x)dx

P(x)dxP(x)dxP(x)dx]Q(x)e

P(x)dxP(x)dx]Q(x)e

两边积分 便得通解

yeP(x)dxP(x)dxQ(x)edxC

P(x)dxP(x)dx或

ye[Q(x)edxC]  22 高等数学教案

§12 微分方程

例3用积分因子求dy2xy4x的通解

dx

解 方程的积分因子为

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

yex2xexy4xex 即(exy)4xex

于是

exy4xexdx2exC 22222222因此原方程的通解为y4xexdx2Cex 22

§12 6 可降阶的高阶微分方程

一、y(n)f(x)型的微分方程

解法 积分n 次

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

y(n2)[f(x)dxC1]dxC2 

  

例1 求微分方程ye2xcos x 的通解

解 对所给方程接连积分三次 得

ye2xsinxC1

ye2xcosxC1xC2

ye2xsinxC1x2C2xC3

这就是所给方程的通解

或

ye2xsinx2C1 1214181212高等数学教案

§12 微分方程

ye2xcosx2C1xC2

ye2xsinxC1x2C2xC3

这就是所给方程的通解

例2 质量为m的质点受力F的作用沿Ox轴作直线运动 设力F仅是时间t的函数FF(t) 在开始时刻t0时F(0)F0 随着时间t的增大 此力F均匀地减小 直到tT时 F(T)0 如果开始时质点位于原点 且初速度为零 求这质点的运动规律

解 设xx(t)表示在时刻t时质点的位置 根据牛顿第二定律 质点运动的微分方程为

m1418d2xF(t)

dt2由题设 力F(t)随t增大而均匀地减小 且t0时 F(0)F0 所以F(t)F0kt 又当tT时 F(T)0 从而

F(t)F0(1)

于是质点运动的微分方程又写为 tTd2xF0(1t)

Tdt2mdx|0 其初始条件为x|t00

dtt0

把微分方程两边积分 得

2F0dxt

(t)C1

dtm2T再积分一次 得

F012t x(t)C1tC2

m26Tdx|0 由初始条件x|t00

dtt0得C1C20

于是所求质点的运动规律为

F012t3

x(t) 0tT

m26T

解 设xx(t)表示在时刻t时质点的位置

根据牛顿第二定律 质点运动的微分方程为

mxF(t)

由题设 F(t)是线性函数 且过点(0 F0)和(T 0) 高等数学教案

§12 微分方程

故

F(t)t1 即F(t)F0(1t) F0TT于是质点运动的微分方程又写为

xF0(1t)

mT其初始条件为x|t00 x|t00

把微分方程两边积分 得

2F0t

x(t)C1 m2T再积分一次 得

F012t3

x(t)C2

m26T由初始条件x|t00 x|t00

得C1C20

于是所求质点的运动规律为

F012t3

x(t) 0tT

m26T

二、y f(x y)型的微分方程

解法 设yp则方程化为

pf(x p)

设pf(x p)的通解为p(xC1) 则

dy(x,C1)

dx原方程的通解为

y(x,C1)dxC2

例3 求微分方程

(1x2)y2xy 满足初始条件

y|x01 y|x03 的特解

解 所给方程是yf(x y)型的 设yp 代入方程并分离变量后 有

dp2xdx

p1x2两边积分 得 高等数学教案

§12 微分方程

ln|p|ln(1x2)C

即

pyC1(1x2)(C1eC)

由条件y|x03 得C13

所以

y3(1x2)

两边再积分 得 yx33xC2

又由条件y|x01 得C21

于是所求的特解为

yx33x1

例4 设有一均匀、柔软的绳索 两端固定 绳索仅受重力的作用而下垂 试问该绳索在平衡状态时是怎样的曲线?

三、yf(y y)型的微分方程

解法 设yp有

y原方程化为 dpdpdydp

pdxdydxdydpf(y,p)

dydp设方程pf(y,p)的通解为yp(y C1) 则原方程的通解为

dy

p

dy(y,C1)xC2

例5 求微分yyy20的通解

解 设yp 则yp代入方程 得

ypdp

dydp2p0

dy

在y0、p0时 约去p并分离变量 得

dpdy

py两边积分得

ln|p|ln|y|lnc

即

pCy或yCy(Cc)

再分离变量并两边积分 便得原方程的通解为 高等数学教案

§12 微分方程

ln|y|Cxlnc1

或

yC1eCx(C1c1)

例5 求微分yyy20的通解

解 设yp 则原方程化为

ypdp2p0

dy当y0、p0时 有

dp1p0

dyy于是

peydy1C1y

即

yC1y0

从而原方程的通解为

yC2e

例6 一个离地面很高的物体受地球引力的作用由静止开始落向地面 求它落 到地面时的速度和所需的时间（不计空气阻力）

§12 7 高阶线性微分方程 一、二阶线性微分方程举例

例1 设有一个弹簧 上端固定 下端挂一个质量为m 的物体 取x 轴铅直向下 并取物体的平衡位置为坐标原点

给物体一个初始速度v00后 物体在平衡位置附近作上下振动 在振动过程中 物体的位置x是t的函数 xx(t)

设弹簧的弹性系数为c 则恢复力fcx

又设物体在运动过程中受到的阻力的大小与速度成正比 比例系数为 则

RC1dxC2eC1x

dx

dt

由牛顿第二定律得

2dxdx

m2cxdtdt

移项 并记2n m k2c

m27 高等数学教案

§12 微分方程

则上式化为

d2x2ndxk2x0

dtdt2这就是在有阻尼的情况下 物体自由振动的微分方程

如果振动物体还受到铅直扰力

FHsin pt 的作用 则有

2dx2ndxk2xhsinpt

dtdt2H其中h 这就是强迫振动的微分方程

m

例2 设有一个由电阻R、自感L、电容C和电源E串联组成的电路 其中R、L、及C为常数 电源电动势是时间t的函数 EEmsint 这里Em及也是常数

设电路中的电流为i(t) 电容器极板上的电量为q(t) 两极板间的电压为uc 自感电动势为EL  由电学知道

iqdqdi uc ELL

Cdtdt根据回路电压定律 得

diqRi0

dtCd2ucduc即

LCRCucEmsint

dtdt EL或写成

d2ucducEm2

2usint

0c2dtLCdtR 1 这就是串联电路的振荡方程 其中02LLC

如果电容器经充电后撤去外电源(E0) 则上述成为

d2ucduc2

20uc0

2dtdt

二阶线性微分方程 二阶线性微分方程的一般形式为

yP(x)yQ(x)yf(x)

若方程右端f(x)0时 方程称为齐次的 否则称为非齐次的

二、线性微分方程的解的结构 高等数学教案

§12 微分方程

先讨论二阶齐次线性方程

d2ydy

yP(x)yQ(x)y0 即2P(x)Q(x)y0

dxdx

定理1 如果函数y1(x)与y2(x)是方程

yP(x)yQ(x)y0的两个解 那么

yC1y1(x)C2y2(x)也是方程的解 其中C1、C2是任意常数

齐次线性方程的这个性质表明它的解符合叠加原理

证明 [C1y1C2y2]C1 y1C2 y2

[C1y1C2y2]C1 y1C2 y2

因为y1与y2是方程yP(x)yQ(x)y0 所以有

y1P(x)y1Q(x)y10及y2P(x)y2Q(x)y20

从而

[C1y1C2y2]P(x)[ C1y1C2y2]Q(x)[ C1y1C2y2]

C1[y1P(x)y1Q(x)y1]C2[y2P(x)y2Q(x)y2]000

这就证明了yC1y1(x)C2y2(x)也是方程yP(x)yQ(x)y0的解

函数的线性相关与线性无关

设y1(x) y2(x)     yn(x)为定义在区间I上的n个函数 如果存在n个不全为零的常数k1 k2     kn 使得当xI 时有恒等式

k1y1(x)k2y2(x)

    knyn(x)0 成立 那么称这n个函数在区间I上线性相关 否则称为线性无关

判别两个函数线性相关性的方法

对于两个函数 它们线性相关与否 只要看它们的比是否为常数 如果比为常数 那么它们就线性相关 否则就线性无关

例如 1 cos2x  sin2x 在整个数轴上是线性相关的 函数1 x x2在任何区间(a, b)内是线性无关的

定理2 如果如果函数y1(x)与y2(x)是方程

yP(x)yQ(x)y0 的两个线性无关的解 那么

yC1y1(x)C2y2(x)(C1、C2是任意常数)高等数学教案

§12 微分方程

是方程的通解

例3 验证y1cos x与y2sin x是方程yy0的线性无关解 并写出其通解

解 因为

y1y1cos xcos x0

y2y2sin xsin x0

所以y1cos x与y2sin x都是方程的解

因为对于任意两个常数k1、k2 要使

k1cos xk2sin x0

只有k1k20 所以cos x与sin x在(, )内是线性无关的

因此y1cos x与y2sin x是方程yy0的线性无关解

方程的通解为yC1cos xC2sin x

例4 验证y1x与y2ex是方程(x1)yxyy0的线性无关解 并写出其通解

解 因为

(x1)y1xy1y10xx0

(x1)y2xy2y2(x1)exxexex0

所以y1x与y2ex都是方程的解

因为比值e x/x 不恒为常数 所以y1x与y2ex在(, )内是线性无关的

因此y1x 与y2ex是方程(x1)yxyy0的线性无关解

方程的通解为yC1xC2e x

推论 如果y1(x) y2(x)    yn(x)是方程

y(n)a1(x)y(n1)    an1(x)y an(x)y0 的n个线性无关的解 那么 此方程的通解为

yC1y1(x)C2y2(x)     Cnyn(x)

其中C1 C2    Cn为任意常数

二阶非齐次线性方程解的结构

我们把方程

yP(x)yQ(x)y0 叫做与非齐次方程 高等数学教案

§12 微分方程

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C1cos xC2sin x 是齐次方程yy0的通解 y*x22是yyx2 的一个特解 因此

yC1cos xC2sin xx22 是方程yyx2的通解

定理4 设非齐次线性微分方程 yP(x)yQ(x)yf(x)的右端f(x)几个函数之和 如

yP(x)yQ(x)yf1(x) f2(x)

而y1*(x)与y2*(x)分别是方程

yP(x)yQ(x)yf1(x)与yP(x)yQ(x)yf2(x)的特解 那么y1*(x)y2*(x)就是原方程的特解

证明提示

[y1y2*]P(x)[ y1*y2*]Q(x)[ y1*y2*]

[ y1*P(x)y1*Q(x)y1*][ y2*P(x)y2*Q(x)y2*]

f1(x)f2(x)

§12 9 二阶常系数齐次线性微分方程

二阶常系数齐次线性微分方程 方程

ypyqy0 称为二阶常系数齐次线性微分方程 其中p、q均为常数

如果y1、y2是二阶常系数齐次线性微分方程的两个线性无关解 那么yC1y1C2y2就是它的

高等数学教案

§12 微分方程

通解

我们看看

能否适当选取r 使yerx

满足二阶常系数齐次线性微分方程 为此将yerx代入方程

ypyqy0 得

(r 2prq)erx 0

由此可见 只要r满足代数方程r2prq0 函数yerx就是微分方程的解

特征方程 方程r2prq0叫做微分方程ypyqy0的特征方程 特征方程的两个根r1、r2可用公式

pp24q

r 1,22求出

特征方程的根与通解的关系

(1)特征方程有两个不相等的实根r1、r2时 函数y1er1x、y2er2x是方程的两个线性无关的解

这是因为

函数y1er1x、y2er2xy1er1x(r1r2)x是方程的解 又不是常数

ey2er2x因此方程的通解为

yC1er1xC2er2x

(2)特征方程有两个相等的实根r1r2时 函数y1er1x、y2xer1x是二阶常系数齐次线性微分方程的两个线性无关的解

这是因为 y1er1x是方程的解 又

r1xr1x2r1x

(xer1x)p(xer1x)q(xer1x)(2r1xr1)ep(1xr1)eqxe r1x

2er1x(2r1p)xe(r1pr1q)0

高等数学教案

§12 微分方程

y2xer1x所以y2xe也是方程的解 且x不是常数

y1er1xr1x

因此方程的通解为

yC1er1xC2xer1x

(3)特征方程有一对共轭复根r1, 2i时 函数ye(i)x、ye(i)x是微分方程的两个线性无关的复数形式的解 函数yexcosx、yexsinx是微分方程的两个线性无关的实数形式的解

函数y1e(i)x和y2e(i)x都是方程的解 而由欧拉公式 得

y1e(i)xex(cosxisinx)

y2e(i)xex(cosxisinx)

1y1y22excosx excosx(y1y2)

21y1y22iexsinx exsinx(y1y2)

2i故excosx、y2exsinx也是方程解

可以验证 y1excosx、y2exsinx是方程的线性无关解

因此方程的通解为

yex(C1cosxC2sinx)

求二阶常系数齐次线性微分方程ypyqy0的通解的步骤为

第一步

写出微分方程的特征方程

r2prq0 第二步

求出特征方程的两个根r1、r2

第三步

根据特征方程的两个根的不同情况 写出微分方程的通解

例1 求微分方程y2y3y0的通解

解 所给微分方程的特征方程为

r22r30 即(r1)(r3)0

其根r11 r23是两个不相等的实根 因此所求通解为

yC1exC2e3x

例2 求方程y2yy0满足初始条件y|x0

4、y| x02的特解

解 所给方程的特征方程为

高等数学教案

§12 微分方程

r22r10 即(r1)20

其根r1r21是两个相等的实根 因此所给微分方程的通解为

y(C1C2x)ex

将条件y|x04代入通解 得C14 从而

y(4C2x)ex

将上式对x求导 得

y(C24C2x)ex

再把条件y|x02代入上式 得C22 于是所求特解为

x(42x)ex

例 3 求微分方程y2y5y 0的通解

解 所给方程的特征方程为

r22r50

特征方程的根为r112i r212i 是一对共轭复根

因此所求通解为

yex(C1cos2xC2sin2x)

n 阶常系数齐次线性微分方程 方程

y(n)p1y(n1)p2 y(n2)     pn1ypny0

称为n 阶常系数齐次线性微分方程 其中 p1

p2      pn1 pn都是常数

二阶常系数齐次线性微分方程所用的方法以及方程的通解形式 可推广到n 阶常系数齐次线性微分方程上去

引入微分算子D 及微分算子的n次多项式

L(D)=Dn p1Dn1p2 Dn2      pn1Dpn 则n阶常系数齐次线性微分方程可记作

(Dn p1Dn1p2 Dn2      pn1Dpn)y0或L(D)y0 注 D叫做微分算子D0yy Dyy D2yy D3yy   Dnyy(n)

分析 令yerx 则

L(D)yL(D)erx(rn p1rn1p2 rn2      pn1rpn)erxL(r)erx

因此如果r是多项式L(r)的根 则yerx是微分方程L(D)y0的解

n 阶常系数齐次线性微分方程的特征方程

高等数学教案

§12 微分方程

L(r)rn p1rn1p2 rn2      pn1rpn0 称为微分方程L(D)y0的特征方程

特征方程的根与通解中项的对应

单实根r 对应于一项 Cerx 

一对单复根r1 2 i 对应于两项 ex(C1cosxC2sinx)

k重实根r对应于k项 erx(C1C2x    Ck xk1)

一对k 重复根r1 2 i 对应于2k项

ex[(C1C2x    Ck xk1)cosx(D1D2x    Dk xk1)sinx]

例4 求方程y(4)2y5y0 的通解

解

这里的特征方程为

r42r35r20 即r2(r22r5)0

它的根是r1r20和r3 412i

因此所给微分方程的通解为

yC1C2xex(C3cos2xC4sin2x)

例5 求方程y(4) 4y0的通解 其中0

解

这里的特征方程为

r4 40

它的根为r1,22(1i) r3,42(1i)

因此所给微分方程的通解为

ye

§12 10 二阶常系数非齐次线性微分方程

二阶常系数非齐次线性微分方程 方程

ypyqyf(x)称为二阶常系数非齐次线性微分方程 其中p、q是常数

二阶常系数非齐次线性微分方程的通解是对应的齐次方程 的通解yY(x)与非齐次方程本身的一个特解yy*(x)之和

2x(C1cos2xC2sin2x)e 2x(C3cos2xC4sin2x) 高等数学教案

§12 微分方程

yY(x) y*(x)

当f(x)为两种特殊形式时 方程的特解的求法

一、f(x)Pm(x)ex 型

当f(x)Pm(x)ex时 可以猜想 方程的特解也应具有这种形式 因此 设特解形式为y*Q(x)ex 将其代入方程 得等式

Q(x)(2p)Q(x)(2pq)Q(x)Pm(x)

(1)如果不是特征方程r2prq0 的根 则2pq0 要使上式成立 Q(x)应设为m 次多项式

Qm(x)b0xmb1xm1    bm1xbm 

通过比较等式两边同次项系数 可确定b0 b1     bm 并得所求特解

y*Qm(x)ex

(2)如果是特征方程 r2prq0 的单根 则2pq0 但2p0 要使等式

Q(x)(2p)Q(x)(2pq)Q(x)Pm(x)

成立 Q(x)应设为m1 次多项式

Q(x)xQm(x)

Qm(x)b0xm b1xm1   

bm1xbm 

通过比较等式两边同次项系数 可确定b0 b1   

 bm 并得所求特解

y*xQm(x)ex

(3)如果是特征方程 r2prq0的二重根 则2pq0 2p0 要使等式

Q(x)(2p)Q(x)(2pq)Q(x)Pm(x)

成立 Q(x)应设为m2次多项式

Q(x)x2Qm(x)

Qm(x)b0xmb1xm1    bm1xbm 

通过比较等式两边同次项系数 可确定b0 b1     bm  并得所求特解

y*x2Qm(x)ex

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

y*xk Qm(x)ex 的特解 其中Qm(x)是与Pm(x)同次的多项式 而k 按不是特征方程的根、是特征方程的单根或

高等数学教案

§12 微分方程

是特征方程的的重根依次取为0、1或2

例1 求微分方程y2y3y3x1的一个特解

解 这是二阶常系数非齐次线性微分方程 且函数f(x)是Pm(x)ex型(其中Pm(x)3x1 0)

与所给方程对应的齐次方程为

y2y3y0

它的特征方程为

r22r30

由于这里0不是特征方程的根 所以应设特解为

y*b0xb1

把它代入所给方程 得

3b0x2b03b13x1

比较两端x同次幂的系数 得

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

y*x

例2 求微分方程y5y6yxe2x的通解

解 所给方程是二阶常系数非齐次线性微分方程 且f(x)是Pm(x)ex型(其中Pm(x)x 2)

与所给方程对应的齐次方程为

y5y6y0

它的特征方程为

r25r 60

特征方程有两个实根r12 r23 于是所给方程对应的齐次方程的通解为

YC1e2xC2e3x 

由于2是特征方程的单根 所以应设方程的特解为

y*x(b0xb1)e2x

把它代入所给方程 得

2b0x2b0b1x

1313高等数学教案

§12 微分方程

比较两端x同次幂的系数 得

2b01 2b01 2b0b10 2bb001由此求得b01 b1 于是求得所给方程的一个特解为 121 y*x(x1)e2x

从而所给方程的通解为

yC1e2xC2e3x(x22x)e2x

提示

y*x(b0xb1)e2x(b0x2b1x)e2x

[(b0x2b1x)e2x][(2b0xb1)(b0x2b1x)2]e2x

[(b0x2b1x)e2x][2b02(2b0xb1)2(b0x2b1x)22]e2x

y*5y*6y*[(b0x2b1x)e2x]5[(b0x2b1x)e2x]6[(b0x2b1x)e2x] [2b02(2b0xb1)2(b0x2b1x)22]e2x5[(2b0xb1)(b0x2b1x)2]e2x6(b0x2b1x)e2x [2b04(2b0xb1)5(2b0xb1)]e2x[2b0x2b0b1]e2x

方程ypyqyex[Pl(x)cosxPn(x)sinx]的特解形式

应用欧拉公式可得

ex[Pl(x)cosxPn(x)sinx] 12ei xei xP(x)ei xei x] n22i1(x)iP(x)]e(i)x1[P(x)iP(x)]e(i)x

[Pnn2l2l

ex[Pl(x)

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

其中P(x)(PlPni) P(x)(PlPni) 而mmax{l n}

设方程ypyqyP(x)e(i)x的特解为y1*xkQm(x)e(i)x

1212高等数学教案

§12 微分方程

则y1*xkQm(x)e(i)必是方程ypyqyP(x)e(i)的特解

其中k按i不是特征方程的根或是特征方程的根依次取0或1

于是方程ypyqyex[Pl(x)cosxPn(x)sinx]的特解为

y*xkQm(x)e(i)xxkQm(x)e(i)x

xkex[Qm(x)(cosxisinx)Qm(x)(cosxisinx)

xk ex[R(1)m(x)cosxR(2)m(x)sinx]

综上所述 我们有如下结论

如果f(x)ex [Pl(x)cosxPn(x)sinx] 则二阶常系数非齐次线性微分方程

ypyqyf(x)的特解可设为

y*xk 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[Pl(x)cosxPn(x)sinx]型(其中0 2 Pl(x)x Pn(x)0）

与所给方程对应的齐次方程为

yy0

它的特征方程为

r210

由于这里i2i 不是特征方程的根 所以应设特解为

y*(axb)cos2x(cxd)sin2x

把它代入所给方程 得

(3ax3b4c)cos2x(3cx3d4a)sin2xxcos2x

比较两端同类项的系数 得 a b0 c0 d于是求得一个特解为 y*xcos2xsin2x

134

91349高等数学教案

§12 微分方程

提示

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

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014由 得a b0 c0 d 3c0394a3d0

§12 12 微分方程的幂级数解法

当微分方程的解不能用初等函数或其积分表达时 我们就要寻求其它解法 常用的有幂级数解法和数值解法 本节我们简单地介绍微分方程的幂级数解法

求一阶微分方程的多项式

f(x y)a00a10(xx0)a01(yy0)    aim(xx0)l(yy0)m

这时我们可以设所求特解可展开为xx0的幂级数

yy0a1(xx0)a2(xx0)2    an(xx0)n    

其中a1 a2     an     是待定的系数 把所设特解代入微分方程中 便得一恒等式 比较这恒等式两端xx0的同次幂的系数 就可定出常数a1 a2     从而得到所求的特解

例1 求方程dy(yy0)f(x,y)满足初始条件y|xx0y0的特解 其中函数f(x y)是(xx0)、dxdyxy2满足y|x00的特解

dx

解 这时x00 y00 故设

ya1xa2x2a3x3a4x4    

把y及y的幂级数展开式代入原方程 得

高等数学教案

§12 微分方程

a12a2x3a3x24a4x35a5x4   

x(a1xa2x2a3x3a4x4   )2

xa12x22a1a2x3(a222a1a3)x4    

由此 比较恒等式两端x的同次幂的系数 得

a10 a2 a30 a40 a5121    

20于是所求解的幂级数展开式的开始几项为

yx2121x5    

定理 如果方程

yP(x)yQ(x)y0 中的系数P(x)与Q(x)可在R

yanxn

n0的解

例2 求微分方程yxy 0的满足初始条件y|x00 y|x01的特解

解 这里P(x)0 Q(x)x在整个数轴上满足定理的条件 因此所求的解可在整个数轴上展开成x的幂级数

ya0a1xa2x2a3x3a4x4    anxn

n0由条件y|x00 得a00 由ya12a2x3a3x24a4x3   及y|x01 得a11 于是

yxa2x2a3x3a4x4    xanxn

y12a2x3a3x24a4x3    1nanxn1

n2n2

y2a232a3x43a4x    n(n1)anxn2 

n22

yxa2xa3xa4x   xanxn 234

n2

y12a2x3a3x4a4x   1nanxn1 23

n2

y2a2x32a3x43a4x2    n(n1)anxn2 

n2 41 高等数学教案

§12 微分方程

把y及y代入方程yxy 0 得

2a232a3x43a4x2    n(n1)anxn2  

x(xa2x2a3x3a4x4  anxn  )0

即

2a232a3x(43a41)x2(54a5a 2)x3

(65a6a3)x4    [(n2)(n1)an2an1]xn    0

于是有

a20, a30, a4一般地 an21, a0, a0,    

6435an1(n3 4   )

(n2)(n1)由递推公式可得

aa411, a80, a90, a107,    

76764310910976431一般地 a3m1(m1 2   )

(3m1)(3m)   7643

a7所求的特解为

yx 1x41x71x10    

4376431097643 42

第三篇：微分方程教案

高等数学教案

第七章

微分方程

教学目的：

1．了解微分方程及其解、阶、通解，初始条件和特等概念。2．熟练掌握变量可分离的微分方程及一阶线性微分方程的解法。

3．会解齐次微分方程、伯努利方程和全微分方程，会用简单的变量代换解某些微分方程。4． 会用降阶法解下列微分方程：y(n)f(x)，yf(x,y)和yf(y,y)5． 理解线性微分方程解的性质及解的结构定理。

6．掌握二阶常系数齐次线性微分方程的解法，并会解某些高于二阶的常系数齐次线性微分方程。

7.求自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程的特解和通解。

8.会解欧拉方程，会解包含两个未知函数的一阶常系数线性微分方程组。9．会解微分方程组（或方程组）解决一些简单的应用问题。教学重点：

1、可分离的微分方程及一阶线性微分方程的解法

(n)

2、可降阶的高阶微分方程yf(x)，yf(x,y)和yf(y,y)

3、二阶常系数齐次线性微分方程；

4、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程；

教学难点：

1、齐次微分方程、伯努利方程和全微分方程；

2、线性微分方程解的性质及解的结构定理；

3、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程的特解。

高等数学教案

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

函数是客观事物的内部联系在数量方面的反映 利用函数关系又可以对客观事物的规律性进行研究 因此如何寻找出所需要的函数关系 在实践中具有重要意义 在许多问题中 往往不能直接找出所需要的函数关系 但是根据问题所提供的情况 有时可以列出含有要找的函数及其导数的关系式 这样的关系就是所谓微分方程 微分方程建立以后 对它进行研究 找出未知函数来 这就是解微分方程

例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x2C

(3)其中C是任意常数

把条件“x1时 y2”代入(3)式 得

212C

由此定出C1 把C1代入(3)式 得所求曲线方程(称为微分方程满足条件y|x12的解)

yx21

例2 列车在平直线路上以20m/s(相当于72km/h)的速度行驶 当制动时列车获得加速度04m/s2 问开始制动后多少时间列车才能停住 以及列车在这段时间里行驶了多少路程?

解 设列车在开始制动后t秒时行驶了s米 根据题意 反映制动阶段列车运动规律的函数ss(t)应满足关系式 d2s0.

(4)dt2此外 未知函数ss(t)还应满足下列条件

t0时 s0 vds20 简记为s|=0 s|=20

(5)

t0t0dt高等数学教案

把(4)式两端积分一次 得

vds0.4tC

(6)1dt再积分一次 得

s02t2 C1t C2

(7)这里C1 C2都是任意常数

把条件v|t020代入(6)得

20C1

把条件s|t00代入(7)得0C2

把C1 C2的值代入(6)及(7)式得

v04t 20

(8)

s02t220t

(9)在(8)式中令v0 得到列车从开始制动到完全停住所需的时间

t2050(s)

0.4再把t50代入(9) 得到列车在制动阶段行驶的路程

s025022050500(m)

几个概念

微分方程 表示未知函数、未知函数的导数与自变量之间的关系的方程 叫微分方程

常微分方程 未知函数是一元函数的微分方程 叫常微分方程

偏微分方程 未知函数是多元函数的微分方程 叫偏微分方程

微分方程的阶 微分方程中所出现的未知函数的最高阶导数的阶数 叫微分方程的阶

x3 yx2 y4xy3x2 

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x0 时 yy0  y y0 

一般写成



yxx0y0 yxx0y0

特解 确定了通解中的任意常数以后 就得到微分方程的特解 即不含任意常数的解

初值问题 求微分方程满足初始条件的解的问题称为初值问题

如求微分方程yf(x

y)满足初始条件yxx0y0的解的问题 记为

yf(x,y)

 yxx0y0

积分曲线 微分方程的解的图形是一条曲线 叫做微分方程的积分曲线

d2xk2x0

例3 验证 函数 xC1cos ktC2 sin kt是微分方程

的解

dt

2解 求所给函数的导数

dxkCsinktkCcoskt 12dtd2xk2Ccosktk2Csinktk2(CcosktCsinkt)



1212dt2d2x将2及x的表达式代入所给方程 得 dt

k2(C1cos ktC2sin kt) k2(C1cos ktC2sin kt)0

d2xk2x0

这表明函数xC1cosktC2sinkt 满足方程2 因此所给函数是所给方程的解

dtd2xk2x0

例4 已知函数xC1cosktC2sinkt(k0)是微分方程2的通解 求满足初始条件

dt

x| t0 A x| t0 0 的特解

高等数学教案

解

由条件x| t0 A及xC1 cos ktC2 sin kt 得

C1A

再由条件x| t0 0 及x(t)kC1sin ktkC2cos kt 得

C20

把C1、C2的值代入xC1cos ktC2sin kt中 得

xAcos kt

作业：P298：4

§7 2 可分离变量的微分方程

观察与分析

1 求微分方程y2x的通解 为此把方程两边积分 得 yx2C

一般地 方程yf(x)的通解为yf(x)dxC(此处积分后不再加任意常数)

2 求微分方程y2xy2 的通解

因为y是未知的 所以积分2xy2dx无法进行 方程两边直

接积分不能求出通解

为求通解可将方程变为

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

y21x2C1  或y2yxC可以验证函数y1是原方程的通解

x2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1dy2xdx (2)3x25xy0

是 dy(3x25x)dx(3)(x2y2)dxxydy=0

不是

(4)y1xy2xy2 是 y(1x)(1y2)(5)y10xy

是 10ydy10xdx(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

解

此方程为可分离变量方程 分离变量后得

1dy2xdx

y1dy2xdx

y两边积分得

即

ln|y|x2C1

从而

yex2C1eC1ex 2因为eC1仍是任意常数 把它记作C 便得所给方程的通解

yCex

例2 铀的衰变速度与当时未衰变的原子的含量M成正比 已知t0时铀的含量为M0 求在衰变过程中铀含量M(t)随时间t变化的规律

解 铀的衰变速度就是M(t)对时间t的导数2dM

dtdMM

dtdM0

dt

由于铀的衰变速度与其含量成正比 故得微分方程其中(>0)是常数 前的曲面号表示当t增加时M单调减少 即由题意 初始条件为 M|t0M0

将方程分离变量得

两边积分 得dMdt

MdM()dt

M即

lnMtlnC 也即MCet

由初始条件 得M0Ce0C

所以铀含量M(t)随时间t变化的规律MM0et 

例3 设降落伞从跳伞塔下落后 所受空气阻力与速度成正比 并设降落伞离开跳伞塔时速度为零 求降落伞下落速度与时间的函数关系

解

设降落伞下落速度为v(t) 降落伞所受外力为Fmgkv(k为比例系数) 根据牛顿第二运

高等数学教案

动定律Fma 得函数v(t)应满足的方程为

mdvmgkv

dt初始条件为

v|t00

方程分离变量 得

dvdt

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

ln(mgkv)1ktC

m1kC1ktmgemCe即

v(C)

kkmg将初始条件v|t00代入通解得C

kktmg(1em)

于是降落伞下落速度与时间的函数关系为vkdy1xy2xy2的通解

例4 求微分方程dx

解 方程可化为

dy(1x)(1y2)

dx分离变量得

1dy(1x)dx

1y21dy(1x)dx 即1x2xC

arctany1y22两边积分得

于是原方程的通解为ytan(x2xC)

作业：P304：1（1）（2）（3）（7）（9）（10），2（2）（4），3 12高等数学教案

§7 3 齐次方程

齐次方程

如果一阶微分方程dyf(x,y)中的函数f(x, y)可写成 dxyy的函数 即f(x,y)() 则称这方程为齐次方程

xx

下列方程哪些是齐次方程？

dyyy2x2dyyy

(1)xyyyx0是齐次方程()21

dxxdxxx22dy1y

2(2)1xy1y不是齐次方程

dx1x222dyx2y2dyxy

(3)(xy)dxxydy0是齐次方程 dxxydxyx22

(4)(2xy4)dx(xy1)dy0不是齐次方程

(5)(2xshdy2xy4

dxxy1yyy3ych)dx3xchdy0是齐次方程

xxxyy2xsh3ychdyxxdy2thyy 

ydxdx3xx3xchx

齐次方程的解法

在齐次方程

ux分离变量 得

ydyy()中 令u 即yux 有 dxxxdu(u)

dxdudx (u)uxdudx(u)ux 两端积分 得

高等数学教案

求出积分后 再用y代替u 便得所给齐次方程的通解

xdydyxy

dxdx

例1 解方程y2x2

解

原方程可写成

y2()dyyx

2ydxxyx1x2因此原方程是齐次方程 令

yux 于是原方程变为

2duu

ux

dxu1yu 则 xdyuxdu

dxdx即

xduu

dxu1分离变量 得

(1)du1udx

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

因为

OAAPOPPMcotOPyx

y高等数学教案

而

OMx2y2

于是得微分方程yxx2y2 y整理得dxx(x)21 这是齐次方程

dyyydxx(x)21

dyyy

问题归结为解齐次方程

令即

yxvdvvv21 即xyv 得vy

ydydvv21

dy分离变量 得dvdy

v21yyy, (v)2v21, CC两边积分 得 ln(vv21)lnylnC, vv21y22yv1

C2C以yvx代入上式 得y22C(xC)

2这是以x轴为轴、焦点在原点的抛物线 它绕x轴旋转所得旋转曲面的方程为

y2z22C(xC) 2这就是所求的旋转曲面方程

例3 设一条河的两岸为平行直线 水流速度为a 有一鸭子从岸边点A游向正对岸点O 设鸭子的游速为b(b>a) 且鸭子游动方向始终朝着点O 已知OAh 求鸭子游过的迹线的方程

解 取O为坐标原点 河岸朝顺水方向为x轴 y 轴指向对岸 设在时刻t鸭子位于点P(x, y) 则鸭子运动速度

v(vx, vy)(dx, dy) 故有dxvx

dyvydtdt高等数学教案

另一方面 vab(a, 0)b(x, y) v(abx, by)

x2y2x2y2x2y2x2y2因此dxvxa(x)21x 即dxa(x)21x

dybyydyvybyydxa(x)21x

dybyy

问题归结为解齐次方程

令

yxu 即xyu 得 yduau21

dyb分离变量 得duady

u21by两边积分 得 arshu(lnylnC) bax1[(Cy)1b(Cy)1b]

将u代入上式并整理 得xy2C以x|yh0代入上式 得Caa1 故鸭子游过的轨迹方程为

haay1by1bh()] 0yh

x[()2hhb将ux代入arshu(lnylnC)后的整理过程

yaarshxb(lnylnC)

yaxshln(Cy)ax1[(Cy)a(Cy)a] yy2bbbbyax[(Cy)(Cy)a]x1[(Cy)1a(Cy)1a]

2C2bbb作业：P309：1（1）（3）（5），2

高等数学教案

§7.4 线性微分方程

一、线性方程

线性方程

方程dyP(x)yQ(x)叫做一阶线性微分方程 dxdydyP(x)y0叫做对应于非齐次线性方程P(x)yQ(x)的齐次线性方程

dxdxdydyy1y0是齐次线性方程 dxdxx2如果Q(x)0  则方程称为齐次线性方程 否则方程称为非齐次线性方程

方程

下列方程各是什么类型方程？

(1)(x2)

(2)3x25x5y0y3x25x  是非齐次线性方程

(3)yy cos xesin x  是非齐次线性方程

(4)dy10xy 不是线性方程 dx23dy3(y1)2dydxxx00或

(5)(y1) 不是线性方程

dxdydx(y1)2x

3齐次线性方程的解法

齐次线性方程

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

y两边积分 得

ln|y|P(x)dxC1

P(x)dx(CeC1)

或

yCe这就是齐次线性方程的通解（积分中不再加任意常数）

例

1求方程(x2)dyy的通解

dx

解

这是齐次线性方程 分离变量得

高等数学教案

dydx

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P(x)P(x)u(x)eQ(x)

u(x)e化简得

u(x)Q(x)eP(x)dx

u(x)Q(x)eP(x)dxdxC

于是非齐次线性方程的通解为

P(x)dxP(x)dx

ye[Q(x)edxC] P(x)dxP(x)dxP(x)dx或

yCeeQ(x)edx 非齐次线性方程的通解等于对应的齐次线性方程通解与非齐次线性方程的一个特解之和

5dy2y(x1)2的通解

例2 求方程dxx1

解

这是一个非齐次线性方程

先求对应的齐次线性方程分离变量得

dy2y0的通解

dxx1dy2dx

yx1两边积分得

ln y2ln(x1)ln C

齐次线性方程的通解为

高等数学教案

yC(x1)2

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

52u(x1)2(x1)2

u(x1)2u(x1)x1 1u(x1)2

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

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

y(x1)[(x1)2C]

323

例3 有一个电路如图所示 其中电源电动势为EEmsint(Em、都是常数) 电阻R和电感L都是常量 求电流i(t)

解

由电学知道 当电流变化时 L上有感应电动势L

EL即

di 由回路电压定律得出

dtdiiR0

dtdiRiE

dtLLdiRiEmsin t

dtLL

把EEmsin t代入上式 得

初始条件为

i|t00

diRiEmsin t为非齐次线性方程 其中

dtLLER t

P(t) Q(t)msinLL

方程由通解公式 得

i(t)eP(t)dtdtdtEP(t)dt[Q(t)edtC]eL(msin teLdtC)

LRRRttEmReL(sinteLdtC)

L高等数学教案

RtEm(Rsin t Lcos t)CeL

222RL其中C为任意常数

将初始条件i|t00代入通解 得C因此 所求函数i(t)为

t LEmREmLe(Rsin t Lcos t)

i(t)222222RLRL LEm

R22L

2二、伯努利方程

伯努利方程 方程

dyP(x)yQ(x)yn(n0 1)dx叫做伯努利方程

下列方程是什么类型方程？

(1)

(2)dy1y1(12x)y4 是伯努利方程 dx33dydyyxy5 yxy5 是伯努利方程 dxdxxy

1(3)y yyxy1 是伯努利方程 yxx

(4)dy2xy4x 是线性方程 不是伯努利方程 dxdyP(x)y1nQ(x)dx

伯努利方程的解法 以yn除方程的两边 得

yn令z y1n  得线性方程

dz(1n)P(x)z(1n)Q(x)

dxdyya(lnx)y2的通解

例4 求方程dxx

解 以y2除方程的两端 得

y2dy11yalnx

dxxd(y1)11yalnx

即

dxx高等数学教案

令zy1 则上述方程成为

dz1zalnx

dxxa2这是一个线性方程 它的通解为

zx[C(lnx)2]

以y1代z  得所求方程的通解为

yx[C(lnx)2]1

经过变量代换 某些方程可以化为变量可分离的方程 或化为已知其求解方法的方程

例

5解方程a2dy1

dxxy

解

若把所给方程变形为

dxxy

dy即为一阶线性方程 则按一阶线性方程的解法可求得通解 但这里用变量代换来解所给方程

令xyu 则原方程化为

du11 即duu1

dxudxuududx

u1分离变量 得

两端积分得

uln|u1|xln|C|

以uxy代入上式 得

yln|xy1|ln|C| 或xCeyy1

作业：P315：1（1）（3）（5）（7）（9），2（1）（3）（5），7（1）（2）

§7 5可降阶的高阶微分方程

高等数学教案

一、y(n)f(x)型的微分方程

解法 积分n 次

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

y(n2)[f(x)dxC1]dxC2 

  

例1 求微分方程ye2xcos x 的通解

解 对所给方程接连积分三次 得

ye2xsinxC1

ye2xcosxC1xC2

ye2xsinxC1x2C2xC3

这就是所给方程的通解

或

ye2xsinx2C1

ye2xcosx2C1xC2

ye2xsinxC1x2C2xC3

这就是所给方程的通解

例2 质量为m的质点受力F的作用沿Ox轴作直线运动 设力F仅是时间t的函数FF(t) 在开始时刻t0时F(0)F0 随着时间t的增大 此力F均匀地减小 直到tT时 F(T)0 如果开始时质点位于原点 且初速度为零 求这质点的运动规律

解 设xx(t)表示在时刻t时质点的位置 根据牛顿第二定律 质点运动的微分方程为

m12141812121418d2xF(t)

2dt由题设 力F(t)随t增大而均匀地减小 且t0时 F(0)F0 所以F(t)F0kt 又当tT时 F(T)0 从而

F(t)F0(1)

于是质点运动的微分方程又写为 tTd2xF0(1t)



Tdt2m高等数学教案

其初始条件为x|t00 dx|0

dtt0

把微分方程两边积分 得

dxF0(tt2)C

1

dtm2T再积分一次 得

F012t x(t)C1tC2

m26T由初始条件x|t00 得C1C20

于是所求质点的运动规律为 dx|0

dtt0F012t3

x(t) 0tT

m26T

二、y f(x y)型的微分方程

解法 设yp则方程化为

pf(x p)

设pf(x p)的通解为p(xC1) 则

dy(x,C1)

dx原方程的通解为

y(x,C1)dxC2

例3 求微分方程

(1x2)y2xy 满足初始条件

y|x01 y|x03 的特解

解 所给方程是yf(x y)型的 设yp 代入方程并分离变量后 有

dp2xdx

p1x2两边积分 得

ln|p|ln(1x2)C

即

pyC1(1x2)(C1eC)

由条件y|x03 得C13

所以

y3(1x2)

高等数学教案

两边再积分 得 yx33xC2

又由条件y|x01 得C21

于是所求的特解为

yx33x1

例4 设有一均匀、柔软的绳索 两端固定 绳索仅受重力的作用而下垂 试问该绳索在平衡状态时是怎样的曲线?

三、yf(y y)型的微分方程

解法 设yp有

y原方程化为 dpdpdydpp

dxdydxdydpf(y,p)

dydpf(y,p)的通解为yp(y C1) 则原方程的通解为 设方程pdy

p

dy(y,C1)xC2

dp

dy

例5 求微分yyy20的通解

解 设yp 则yp代入方程 得

ypdp2p0

dy

在y0、p0时 约去p并分离变量 得

dpdy

py两边积分得

ln|p|ln|y|lnc

即

pCy或yCy(Cc)

再分离变量并两边积分 便得原方程的通解为

ln|y|Cxlnc1

或

yC1eCx(C1c1)

作业：P323：1（1）（3）（5）（7）（9），2（1）（3）（5）

高等数学教案

§7 6 高阶线性微分方程 一、二阶线性微分方程举例

例1 设有一个弹簧 上端固定 下端挂一个质量为m 的物体 取x 轴铅直向下 并取物体的平衡位置为坐标原点

给物体一个初始速度v00后 物体在平衡位置附近作上下振动 在振动过程中 物体的位置x是t的函数 xx(t)

设弹簧的弹性系数为c 则恢复力fcx

又设物体在运动过程中受到的阻力的大小与速度成正比 比例系数为 则

Rdx

dt

由牛顿第二定律得

2dxdx

m2cx

dtdt

移项 并记2nc k2

mmd2x2ndxk2x0则上式化为



dtdt2这就是在有阻尼的情况下 物体自由振动的微分方程

如果振动物体还受到铅直扰力

FHsin pt 的作用 则有

d2x2ndxk2xhsinpt



dtdt2H其中h 这就是强迫振动的微分方程

m

例2 设有一个由电阻R、自感L、电容C和电源E串联组成的电路 其中R、L、及C为常

高等数学教案

数 电源电动势是时间t的函数 EEmsint 这里Em及也是常数

设电路中的电流为i(t) 电容器极板上的电量为q(t) 两极板间的电压为uc 自感电动势为EL  由电学知道

iqdqdi uc ELL

CdtdtdiqRi0

dtC根据回路电压定律 得

ELd2ucducRCucEmsint

即

LCdtdt2或写成

d2ucducEm22usint

0c2dtLCdtR 1 这就是串联电路的振荡方程 其中02LLC

如果电容器经充电后撤去外电源(E0) 则上述成为

d2ucduc220uc0

2dtdt

二阶线性微分方程 二阶线性微分方程的一般形式为

yP(x)yQ(x)yf(x)

若方程右端f(x)0时 方程称为齐次的 否则称为非齐次的

二、线性微分方程的解的结构

先讨论二阶齐次线性方程

d2ydyQ(x)y0

yP(x)yQ(x)y0 即2P(x)dxdx

定理

1如果函数y1(x)与y2(x)是方程

yP(x)yQ(x)y0的两个解 那么

yC1y1(x)C2y2(x)也是方程的解 其中C1、C2是任意常数

齐次线性方程的这个性质表明它的解符合叠加原理

证明 [C1y1C2y2]C1 y1C2 y2

高等数学教案

[C1y1C2y2]C1 y1C2 y2

因为y1与y2是方程yP(x)yQ(x)y0 所以有

y1P(x)y1Q(x)y10及y2P(x)y2Q(x)y20

从而

[C1y1C2y2]P(x)[ C1y1C2y2]Q(x)[ C1y1C2y2]

C1[y1P(x)y1Q(x)y1]C2[y2P(x)y2Q(x)y2]000

这就证明了yC1y1(x)C2y2(x)也是方程yP(x)yQ(x)y0的解

函数的线性相关与线性无关

设y1(x) y2(x)     yn(x)为定义在区间I上的n个函数 如果存在n个不全为零的常数k1 k2     kn 使得当xI 时有恒等式

k1y1(x)k2y2(x)

    knyn(x)0 成立 那么称这n个函数在区间I上线性相关 否则称为线性无关

判别两个函数线性相关性的方法

对于两个函数 它们线性相关与否 只要看它们的比是否为常数 如果比为常数 那么它们就线性相关 否则就线性无关

例如 1 cos2x  sin2x 在整个数轴上是线性相关的 函数1 x x2在任何区间(a, b)内是线性无关的

定理2 如果如果函数y1(x)与y2(x)是方程

yP(x)yQ(x)y0 的两个线性无关的解 那么

yC1y1(x)C2y2(x)(C1、C2是任意常数)是方程的通解

例3 验证y1cos x与y2sin x是方程yy0的线性无关解 并写出其通解

解 因为

y1y1cos xcos x0

y2y2sin xsin x0

所以y1cos x与y2sin x都是方程的解

因为对于任意两个常数k1、k2 要使

k1cos xk2sin x0

只有k1k20 所以cos x与sin x在(, )内是线性无关的

因此y1cos x与y2sin x是方程yy0的线性无关解

高等数学教案

方程的通解为yC1cos xC2sin x

例4 验证y1x与y2ex是方程(x1)yxyy0的线性无关解 并写出其通解

解 因为

(x1)y1xy1y10xx0

(x1)y2xy2y2(x1)exxexex0

所以y1x与y2ex都是方程的解

因为比值e x/x 不恒为常数 所以y1x与y2ex在(, )内是线性无关的

因此y1x 与y2ex是方程(x1)yxyy0的线性无关解

方程的通解为yC1xC2e x

推论 如果y1(x) y2(x)    yn(x)是方程

y(n)a1(x)y(n1)    an1(x)y an(x)y0 的n个线性无关的解 那么 此方程的通解为

yC1y1(x)C2y2(x)     Cnyn(x)

其中C1 C2    Cn为任意常数

二阶非齐次线性方程解的结构

我们把方程

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C1cos xC2sin x 是齐次方程yy0的通解 y*x22是yyx2 的一个特解 因此

yC1cos xC2sin xx22

高等数学教案

是方程yyx2的通解

定理4 设非齐次线性微分方程 yP(x)yQ(x)yf(x)的右端f(x)几个函数之和 如

yP(x)yQ(x)yf1(x) f2(x)

而y1*(x)与y2*(x)分别是方程

yP(x)yQ(x)yf1(x)与yP(x)yQ(x)yf2(x)的特解 那么y1*(x)y2*(x)就是原方程的特解

证明提示

[y1y2*]P(x)[ y1*y2*]Q(x)[ y1*y2*]

[ y1*P(x)y1*Q(x)y1*][ y2*P(x)y2*Q(x)y2*]

f1(x)f2(x)

作业：P331：1（1）（3）（5）（7），4（1）（3）（5）

§7 7 二阶常系数齐次线性微分方程

二阶常系数齐次线性微分方程 方程 ypyqy0 称为二阶常系数齐次线性微分方程 其中p、q均为常数

如果y1、y2是二阶常系数齐次线性微分方程的两个线性无关解 那么yC1y1C2y2就是它的通解

我们看看

能否适当选取r 使yerx

满足二阶常系数齐次线性微分方程 为此将yerx代入方程

ypyqy0 得

(r 2prq)erx 0

由此可见 只要r满足代数方程r2prq0 函数yerx就是微分方程的解

特征方程 方程r2prq0叫做微分方程ypyqy0的特征方程 特征方程的两个根r1、r2可用公式

pp24q

r 1,22高等数学教案

求出

特征方程的根与通解的关系

(1)特征方程有两个不相等的实根r1、r2时 函数y1er1x、y2er2x是方程的两个线性无关的解

这是因为

函数y1e因此方程的通解为

yC1er1xC2er2x

(2)特征方程有两个相等的实根r1r2时 函数y1er1x、y2xer1x是二阶常系数齐次线性微分方程的两个线性无关的解

这是因为 y1er1x是方程的解 又

r1xr1x2r1x

(xer1x)p(xer1x)q(xer1x)(2r1xr1xr1)ep(1)eqxe r1x

2er1x(2r1p)xe(r1pr1q)0 r1x、y2er2xy1er1x(r1r2)x是方程的解 又不是常数

ey2er2xy2xer1xx不是常数

所以y2xe也是方程的解 且y1er1xr1x

因此方程的通解为

yC1er1xC2xer1x

(3)特征方程有一对共轭复根r1, 2i时 函数ye(i)x、ye(i)x是微分方程的两个线性无关的复数形式的解 函数yexcosx、yexsinx是微分方程的两个线性无关的实数形式的解

函数y1e(i)x和y2e(i)x都是方程的解 而由欧拉公式 得

y1e(i)xex(cosxisinx)

y2e(i)xex(cosxisinx)

1y1y22excosx excosx(y1y2)

2高等数学教案

1y1y22iexsinx exsinx(y1y2)

2i故excosx、y2exsinx也是方程解

可以验证 y1excosx、y2exsinx是方程的线性无关解

因此方程的通解为

yex(C1cosxC2sinx)

求二阶常系数齐次线性微分方程ypyqy0的通解的步骤为

第一步

写出微分方程的特征方程

r2prq0 第二步

求出特征方程的两个根r1、r2

第三步

根据特征方程的两个根的不同情况 写出微分方程的通解

例1 求微分方程y2y3y0的通解

解 所给微分方程的特征方程为

r22r30 即(r1)(r3)0

其根r11 r23是两个不相等的实根 因此所求通解为

yC1exC2e3x

例2 求方程y2yy0满足初始条件y|x0

4、y| x02的特解

解 所给方程的特征方程为

r22r10 即(r1)20

其根r1r21是两个相等的实根 因此所给微分方程的通解为

y(C1C2x)ex

将条件y|x04代入通解 得C14 从而

y(4C2x)ex

将上式对x求导 得

y(C24C2x)ex

再把条件y|x02代入上式 得C22 于是所求特解为

x(42x)ex

例 3 求微分方程y2y5y 0的通解

解 所给方程的特征方程为

r22r50

高等数学教案

特征方程的根为r112i r212i 是一对共轭复根

因此所求通解为

yex(C1cos2xC2sin2x)

n 阶常系数齐次线性微分方程 方程

y(n)p1y(n1)p2 y(n2)     pn1ypny0

称为n 阶常系数齐次线性微分方程 其中 p1

p2      pn1 pn都是常数

二阶常系数齐次线性微分方程所用的方法以及方程的通解形式 可推广到n 阶常系数齐次线性微分方程上去

引入微分算子D 及微分算子的n次多项式

L(D)=Dn p1Dn1p2 Dn2      pn1Dpn 则n阶常系数齐次线性微分方程可记作

(Dn p1Dn1p2 Dn2      pn1Dpn)y0或L(D)y0 注 D叫做微分算子D0yy Dyy D2yy D3yy   Dnyy(n)

分析 令yerx 则

L(D)yL(D)erx(rn p1rn1p2 rn2      pn1rpn)erxL(r)erx

因此如果r是多项式L(r)的根 则yerx是微分方程L(D)y0的解

n 阶常系数齐次线性微分方程的特征方程

L(r)rn p1rn1p2 rn2      pn1rpn0 称为微分方程L(D)y0的特征方程

特征方程的根与通解中项的对应

单实根r 对应于一项 Cerx 

一对单复根r1 2 i 对应于两项 ex(C1cosxC2sinx)

k重实根r对应于k项 erx(C1C2x    Ck xk1)

一对k 重复根r1 2 i 对应于2k项

ex[(C1C2x    Ck xk1)cosx(D1D2x    Dk xk1)sinx]

例4 求方程y(4)2y5y0 的通解

解

这里的特征方程为

r42r35r20 即r2(r22r5)0

它的根是r1r20和r3 412i

因此所给微分方程的通解为

高等数学教案

yC1C2xex(C3cos2xC4sin2x)

例5 求方程y(4) 4y0的通解 其中0

解

这里的特征方程为

r4 40

它的根为r1,22(1i) r3,42(1i)

因此所给微分方程的通解为

ye2x(C1cos2xC2sin2x)e 2x(C3cos2xC4sin2x)

作业：P340：1（1）（3）（2）（4）（5）（6）（8），2（2）（4）（6）

§7 8 二阶常系数非齐次线性微分方程

二阶常系数非齐次线性微分方程 方程

ypyqyf(x)称为二阶常系数非齐次线性微分方程 其中p、q是常数

二阶常系数非齐次线性微分方程的通解是对应的齐次方程 的通解yY(x)与非齐次方程本身的一个特解yy*(x)之和

yY(x) y*(x)

当f(x)为两种特殊形式时 方程的特解的求法

一、f(x)Pm(x)ex 型

当f(x)Pm(x)ex时 可以猜想 方程的特解也应具有这种形式 因此 设特解形式为y*Q(x)ex 将其代入方程 得等式

Q(x)(2p)Q(x)(2pq)Q(x)Pm(x)

(1)如果不是特征方程r2prq0 的根 则2pq0 要使上式成立 Q(x)应设为m 次多项式

高等数学教案

Qm(x)b0xmb1xm1    bm1xbm 

通过比较等式两边同次项系数 可确定b0 b1     bm 并得所求特解

y*Qm(x)ex

(2)如果是特征方程 r2prq0 的单根 则2pq0 但2p0 要使等式

Q(x)(2p)Q(x)(2pq)Q(x)Pm(x)

成立 Q(x)应设为m1 次多项式

Q(x)xQm(x)

Qm(x)b0xm b1xm1   

bm1xbm 

通过比较等式两边同次项系数 可确定b0 b1   

 bm 并得所求特解

y*xQm(x)ex

(3)如果是特征方程 r2prq0的二重根 则2pq0 2p0 要使等式

Q(x)(2p)Q(x)(2pq)Q(x)Pm(x)

成立 Q(x)应设为m2次多项式

Q(x)x2Qm(x)

Qm(x)b0xmb1xm1    bm1xbm 

通过比较等式两边同次项系数 可确定b0 b1     bm  并得所求特解

y*x2Qm(x)ex

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

y*xk Qm(x)ex 的特解 其中Qm(x)是与Pm(x)同次的多项式 而k 按不是特征方程的根、是特征方程的单根或是特征方程的的重根依次取为0、1或2

例1 求微分方程y2y3y3x1的一个特解

解 这是二阶常系数非齐次线性微分方程 且函数f(x)是Pm(x)ex型(其中Pm(x)3x1 0)

与所给方程对应的齐次方程为

y2y3y0

它的特征方程为

r22r30

由于这里0不是特征方程的根 所以应设特解为

y*b0xb1

高等数学教案

把它代入所给方程 得

3b0x2b03b13x1

比较两端x同次幂的系数 得

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

y*x

例2 求微分方程y5y6yxe2x的通解

解 所给方程是二阶常系数非齐次线性微分方程 且f(x)是Pm(x)ex型(其中Pm(x)x 2)

与所给方程对应的齐次方程为

y5y6y0

它的特征方程为

r25r 60

特征方程有两个实根r12 r23 于是所给方程对应的齐次方程的通解为

YC1e2xC2e3x 

由于2是特征方程的单根 所以应设方程的特解为

y*x(b0xb1)e2x

把它代入所给方程 得

2b0x2b0b1x

比较两端x同次幂的系数 得

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

y*x(x1)e2x

从而所给方程的通解为

yC1e2xC2e3x(x22x)e2x 121212高等数学教案

提示

y*x(b0xb1)e2x(b0x2b1x)e2x

[(b0x2b1x)e2x][(2b0xb1)(b0x2b1x)2]e2x

[(b0x2b1x)e2x][2b02(2b0xb1)2(b0x2b1x)22]e2x

y*5y*6y*[(b0x2b1x)e2x]5[(b0x2b1x)e2x]6[(b0x2b1x)e2x] [2b02(2b0xb1)2(b0x2b1x)22]e2x5[(2b0xb1)(b0x2b1x)2]e2x6(b0x2b1x)e2x [2b04(2b0xb1)5(2b0xb1)]e2x[2b0x2b0b1]e2x

方程ypyqyex[Pl(x)cosxPn(x)sinx]的特解形式

应用欧拉公式可得

ex[Pl(x)cosxPn(x)sinx]

ex[P(x)eli xei xP(x)ei xei x] n22i

[Pe(i)x[Pe(i)x

l(x)iPn(x)]l(x)iPn(x)]

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

其中P(x)(PlPni) P(x)(PlPni) 而mmax{l n}

设方程ypyqyP(x)e(i)x的特解为y1*xkQm(x)e(i)x

则y1*xkQm(x)e(i)必是方程ypyqyP(x)e(i)的特解

其中k按i不是特征方程的根或是特征方程的根依次取0或1

于是方程ypyqyex[Pl(x)cosxPn(x)sinx]的特解为

y*xkQm(x)e(i)xxkQm(x)e(i)x

xkex[Qm(x)(cosxisinx)Qm(x)(cosxisinx)

xk ex[R(1)m(x)cosxR(2)m(x)sinx]

综上所述 我们有如下结论

如果f(x)ex [Pl(x)cosxPn(x)sinx] 则二阶常系数非齐次线性微分方程 12121212高等数学教案

ypyqyf(x)的特解可设为

y*xk 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[Pl(x)cosxPn(x)sinx]型(其中0 2 Pl(x)x Pn(x)0）

与所给方程对应的齐次方程为

yy0

它的特征方程为

r210

由于这里i2i 不是特征方程的根 所以应设特解为

y*(axb)cos2x(cxd)sin2x

把它代入所给方程 得

(3ax3b4c)cos2x(3cx3d4a)sin2xxcos2x

比较两端同类项的系数 得 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

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 134

91349高等数学教案

3a13b4c014由 得a b0 c0 d 3c0394a3d0作业：P347：1（1）（2）（5）（9）2（2）（3）（4）

第四篇：第四章 微分方程讲稿

高等数学C教案

第四章

微分方程

第四章

微分方程

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

导入：(8分钟)函数是客观事物的内部联系在数量方面的反映 利用函数关系又可以对客观事物的规律性进行研究 因此如何寻找出所需要的函数关系 在实践中具有重要意义 在许多问题中 往往不能直接找出所需要的函数关系 但是根据问题所提供的情况 有时可以列出含有要找的函数及其导数的关系式 这样的关系就是所谓微分方程 微分方程建立以后 对它进行研究 找出未知函数来 这就是解微分方程

引例 一曲线通过点(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x2C

(3)其中C是任意常数

把条件“x1时 y2”代入(3)式 得

212C

由此定出C1 把C1代入(3)式 得所求曲线方程(称为微分方程满足条件y|x12的解)

yx21

几个概念

微分方程 表示未知函数、未知函数的导数与自变量之间的关系的方程 叫微分方程

常微分方程 未知函数是一元函数的微分方程 叫常微分方程

偏微分方程 未知函数是多元函数的微分方程 叫偏微分方程

微分方程的阶 微分方程中所出现的未知函数的最高阶导数的阶数 叫微分方程的阶

x3 yx2 y4xy3x2 

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))  高等数学C教案

第四章

微分方程

微分方程的解 满足微分方程的函数(把函数代入微分方程能使该方程成为恒等式)叫做该微分方程的解 确切地说 设函数y(x)在区间I上有n阶连续导数 如果在区间I上

F[x (x) (x)    (n)(x)]0

那么函数y(x)就叫做微分方程F(x y y    y(n))0在区间I上的解

通解 如果微分方程的解中含有任意常数 且任意常数的个数与微分方程的阶数相同 这样的解叫做微分方程的通解

初始条件 用于确定通解中任意常数的条件 称为初始条件 如

xx0 时 yy0  y y0 

一般写成



yxx0y0 yxx0y0

特解 确定了通解中的任意常数以后 就得到微分方程的特解 即不含任意常数的解

初值问题 求微分方程满足初始条件的解的问题称为初值问题

如求微分方程yf(x

y)满足初始条件yxx0y0的解的问题 记为

yf(x,y)

 yxx0y0

积分曲线 微分方程的解的图形是一条曲线 叫做微分方程的积分曲线

§4 2 一阶微分方程

导入：（8分钟）1 求微分方程y2x的通解 为此把方程两边积分 得

yx2C

一般地 方程yf(x)的通解为yf(x)dxC(此处积分后不再加任意常数)

2 求微分方程y2xy2 的通解

因为y是未知的 所以积分2xy2dx无法进行 方程两边直接积分不能求出通解



为求通解可将方程变为

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

y x2C 或y可以验证函数y1y1

x2C1是原方程的通解

x2C

g(y)dyf(x)dx

一般地 如果一阶微分方程y(x, y)能写成

形式 则两边积分可得一个不含未知函数的导数的方程

G(y)F(x)C

由方程G(y)F(x)C所确定的隐函数就是原方程的通解 高等数学C教案

第四章

微分方程

对称形式的一阶微分方程

一阶微分方程有时也写成如下对称形式

P(x y)dxQ(x y)dy0 在这种方程中 变量x与y 是对称的

若把x看作自变量、y看作未知函数 则当Q(x,y)0时 有

dyP(x,y)

dxQ(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1dy2xdx (2)3x25xy0

是 dy(3x25x)dx(3)(x2y2)dxxydy=0

不是

(4)y1xy2xy2 是 y(1x)(1y2)(5)y10xy

是 10ydy10xdx(6)ydxQ(x,y)

dyP(x,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

解

此方程为可分离变量方程 分离变量后得

1dy2xdx

y两边积分得

1dy2xdx

y3 高等数学C教案

第四章

微分方程

即

ln|y|x2C1

从而

yex2C1eC1ex

2因为eC1仍是任意常数 把它记作C 便得所给方程的通解

yCex

例2 铀的衰变速度与当时未衰变的原子的含量M成正比 已知t0时铀的含量为M0 求在衰变过程中铀含量M(t)随时间t变化的规律

解 铀的衰变速度就是M(t)对时间t的导数

2dM

dt

由于铀的衰变速度与其含量成正比 故得微分方程

dMM

dtdM0

dt其中(>0)是常数 前的曲面号表示当t增加时M单调减少 即由题意 初始条件为

M|t0M0

将方程分离变量得

两边积分 得

dMdt

MdM()dt

M即lnMtlnC 也即MCet

由初始条件 得M0Ce0C

所以铀含量M(t)随时间t变化的规律MM0et 

例3 设降落伞从跳伞塔下落后 所受空气阻力与速度成正比 并设降落伞离开跳伞塔时速度为零 求降落伞下落速度与时间的函数关系

解

设降落伞下落速度为v(t) 降落伞所受外力为Fmgkv(k为比例系数) 根据牛顿第二运动定律Fma 得函数v(t)应满足的方程为

m初始条件为

v|t00

方程分离变量 得

两边积分 得

dvmgkv

dtdvdt

mgkvmdvdt

mgkvm 高等数学C教案

第四章

微分方程

ln(mgkv)kC1ktmgemCe(C即v)

kk1ktC

m1将初始条件v|t00代入通解得Cmg

kktmg(1em)

于是降落伞下落速度与时间的函数关系为vk

例4 求微分方程

解 方程可化为 dy1xy2xy2的通解

dx

dy(1x)(1y2)

dx1dy(1x)dx

1y2分离变量得

两边积分得

1dy(1x)dx1x2xC

 即arctany1y22于是原方程的通解为ytan(x2xC)

例5 有高为1m的半球形容器 水从它的底部小孔流出 小孔横截面面积为1cm2 开始时容器内盛满了水 求水从小孔流出过程中容器里水面高度h随时间t变化的规律

解 由水力学知道 水从孔口流出的流量Q可用下列公式计算

Q12dV0.62S2gh

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

dV0.622gh 或dV0.622ghdt

dt

dVr2dh

另一方面 设在微小时间间隔[t tdt]内 水面高度由h降至hdh(dh0) 则又可得到

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

r1002(100h)2200hh2

所以

dV(200hh2)dh

通过比较得到

0.622ghdt(200hh2)dh 高等数学C教案

第四章

微分方程

这就是未知函数hh(t)应满足的微分方程

此外 开始时容器内的水是满的 所以未知函数hh(t)还应满足下列初始条件

h|t0100

将方程0.622ghdt(200hh2)dh分离变量后得

dt两端积分 得

t350.622g13(200h2h2)dh

0.622g13(200h2h2)dh

即 t(400h22h2)C

50.622g3其中C是任意常数

由初始条件得

t(400100221002)C

50.622gC35(400000200000)14105

350.622g0.622g15

因此t0.622g(7105353210h3h2)

上式表达了水从小孔流出的过程中容器内水面高度h与时间t之间的函数关系 二、一阶线性微分方程

方程dyP(x)yQ(x)叫做一阶线性微分方程 dxdydyP(x)y0叫做对应于非齐次线性方程P(x)yQ(x)的齐次线性方程

dxdxdydyy1y0是齐次线性方程

dxx2dx如果Q(x)0  则方程称为齐次线性方程 否则方程称为非齐次线性方程

方程

提问：下列方程各是什么类型方程？

(1)(x2)

(2)3x25x5y0y3x25x  是非齐次线性方程

(3)yy cos xesin x  是非齐次线性方程

(4)dy10xy 不是线性方程 dx6 高等数学C教案

第四章

微分方程

3(y1)2dydy3dxx0或

(5)(y1) 不是线性方程 x032dydxx(y1)dx21、齐次线性方程的解法

齐次线性方程dyP(x)y0是变量可分离方程 分离变量后得 dx

dyP(x)dx

y两边积分 得

ln|y|P(x)dxC1

P(x)dx(CeC1)

或

yCe这就是齐次线性方程的通解（积分中不再加任意常数）

例6 求方程(x2)dyy的通解

dxdydx

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P(x)P(x)u(x)eQ(x)

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

u(x)Q(x)eP(x)dxdxC

于是非齐次线性方程的通解为

P(x)dxP(x)dx

ye[Q(x)edxC]

P(x)dxP(x)dxP(x)dx或

yCeeQ(x)edx

 高等数学C教案

第四章

微分方程

非齐次线性方程的通解等于对应的齐次线性方程通解与非齐次线性方程的一个特解之和

5dy2y(x1)2的通解

例7 求方程dxx

1解

这是一个非齐次线性方程

先求对应的齐次线性方程分离变量得

两边积分得

ln y2ln(x1)ln C

齐次线性方程的通解为

yC(x1)2

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

dy2y0的通解

dxx1dy2dx

yx12u(x1)2(x1)2

u(x1)2u(x1)x12

5两边积分 得

1u(x1)2 u(x1)2C

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

y(x1)[(x1)2C]

3例8 有一个电路如图所示 其中电源电动势为EEmsint(Em、都是常数) 电阻R和电感L都是常量 求电流i(t)

解

由电学知道 当电流变化时 L上有感应电动势L

EL即

di 由回路电压定律得出 dtdiiR0

dtdiRiE

dtLL

把EEmsin t代入上式 得

初始条件为

diRiEmsin t

dtLL8 高等数学C教案

第四章

微分方程

i|t00

方程diRiEmsin t为非齐次线性方程 其中 dtLLER

P(t) Q(t)msin t

LLdtdtEP(t)dt[Q(t)edtC]eL(msin teLdtC)

LRR由通解公式 得

i(t)eP(t)dtRttEmReL(sinteLdtC)

LRtEm(Rsin t Lcos t)CeL

222RL其中C为任意常数

将初始条件i|t00代入通解 得C因此 所求函数i(t)为

t LEmREmLe(Rsin t Lcos t)

i(t)2R2L2R22L2 LEm

R22L2总结：

1、微分方程的相关概念

a、微分方程的阶

b、微分方程的通解与特解

2、可分离变量的微分方程

a、可分离变量的微分方程

b、可转化为可分离变量的微分方程

3、一阶线性微分方程

a、一阶线性齐次微分方程

b、一阶线性非齐次微分方程

c、常数变易法 教学后记：高等数学C教案

第四章

微分方程

作业：

第五篇：微分方程习题答案

微分方程习题答案

习题基本要求：微分方程的阶，判定一阶齐次（非齐次）微分方程，微分方程的通解及特解，可分离变量微分方程及其通解，二阶常系数微分方程的特征根及其三种不同形式的通解，选择题

下列方程哪些是一阶齐次微分方程？ dyyy2x2dyyy(1)xyyyx0是齐次方程(21dxxdxxx2

2dyy2(2)xyy不是齐次方程dx1x22

dyx2y2dyxy(3)(xy)dxxydy0是齐次方程 dxxydxyx22

(4)(2xy4)dx(xy1)dy0不是齐次方程dy2xy4dxxy

1y2()dyydy22dyxy(5)yx是齐次方程dxdxdxxyx21x21、微分方程y“+（yˊ）4-y3=0的阶数是（B）

（A）1（B）2（C）3（D）

42、方程(y-3x)dx –(x+y)dy=0是（B）

（A）可分离变量微分方程（B）齐次方程

（C）一阶非齐次线性微分方程（D）一阶齐次线性微分方程

3、方程xdy+ydx=0的通解为（D）

（A）xy=1（B）xy=3（C）xy=-3（D）xy=C4、方程y”+ yˊ-2 y=0的通解为（C）

----（A）y=e2x+ex（B）y=Ce2x+ex（C）y=C1e2x+C2ex（D）y=e2x+Cex

填空题：

1、方程ydy+xdx=0的通解为22.通解为y=Cex的一阶微分方程为yˊ-y=0.2、满足条件y(0)=3的微分方程dy=2xydx的特解为y=3ex2.3、二阶常系数齐次线性微分方程y“+p yˊ+q y=0的特征方程为r2-

4、微分方程y”-4y=0的通解为2x2x.-

5、微分方程y“-4yˊ-5y=0的通解为x5x6、微分方程y”-4yˊ+13y=0的通解为

7、微分方程y“+2yˊ+y=0的通解解答题

1、求可分离变量微分方程dy=xydx的通解。

解：（1）显然y=0是微分方程的解；

（2）当y≠0时，方程可化为dydyxdx，两边分别积分xdx yy

12x12得方程的解为lnyxC1，即yCe2

212x2由（1）（2）可知微分方程的通解为yCe。

2、求微分方程ex-ydx=dy的通解。

解：方程可化为exdx=eydy，两边积分得∫exdx=∫eydy，于是微分方程的通解为ey = ex+C.3、求微分方程y”-2yˊ-3y=0的通解。

-解：所给微分方程的特征方程为r2-2r-3=0，其根为r1=-1,r2=3，因此所求通解为y=C1ex+C2e3x4、求微分方程y“-5yˊ+6y=0的通解。

解：所给微分方程的特征方程为r2-5r+6=0,其根为r1=2,r2=3.因此所求通解为y=C1e2x+C2e3x。

5、求微分方程y”+2yˊ+y=0的通解。

-解：所给微分方程的特征方程为r2+2r+1=0,其根为r1=r2=-1.因此所求通解为y=(C1+C2x)ex.6、求微分方程y“-4yˊ+4y=0的通解。

解：所给微分方程的特征方程为r2-4r+4=0，其根为r1=r2=2，因此所求通解为y=(C1+C2x)e2x.7、求微分方程y”-2 yˊ+5 y=0的通解。

解：所给方程的特征方程为r2-2r+5=0，其根为r

因此所求通解为y=ex(C1cos2x+C2sin2x)

8、求微分方程y"-4 yˊ+5 y=0的通解。

解：所给方程的特征方程为r2-2r+5=0,其根为r

因此所求通解为y=e2x(C1cosx+C2sinx).12i 2i