微分方程PPT(罗兆富等编)第七章 特征线法、达朗贝尔公式和分离变量法
合集下载
可分离变量的微分方程ppt课件
§.2 可分离变量的微分 方程
1
❖可分离变量的微分方程 如果一个一阶微分方程能写成
g(y)dyf(x)dx (或写成y(x)(y))
的形式 那么原方程就称为可分离变量的微分方程 ❖可分离变量的微分方程的解法
•分离变量 将方程写成g(y)dy f(x)dx的形式
•两端积分 g(y)dy f (x)dx 设积分后得 G(y)F(x)C
)
k
7
结束
m
dv dt
mg
kv
v |t 0 0
提示 降落伞所受外力为Fmgkv(k为比例系数)
牛顿第二运动定律Fma 6 下页
例3 设降落伞从跳伞塔 两边积分得
下落后 所受空气阻力与速度 成正比 并设降落伞离开跳伞
dv mg kv
dt m
塔时速度为零 求降落伞下落 即 速度与时间的函数关系
1 k
ln(
dx 设积分后得 G(y)F(x)C
方程由G(y)F(x)C所确定的函数,称为隐式(通)解 2
讨论
微分方程
分离变量
y2xy
y1dy2xdx
3x25xy0
dy(3x25x)dx
(x2y2)dxxydy=0 ————
y1xy2xy2 y(1x)(1y2)
y10xy y x y
yx
10ydy10xdx ————
是否可分离变量 是 是 不是 是 是
不是
3
例 1 求微分方程 dy 2xy 的通解 dx
解 这是一个可分离变量的微分方程
分离变量得
1 y
dy 2xdx
两边积分得
1 dy y
2xdx
即
ln|y|x2C1 注 加常数的另一方法
1
❖可分离变量的微分方程 如果一个一阶微分方程能写成
g(y)dyf(x)dx (或写成y(x)(y))
的形式 那么原方程就称为可分离变量的微分方程 ❖可分离变量的微分方程的解法
•分离变量 将方程写成g(y)dy f(x)dx的形式
•两端积分 g(y)dy f (x)dx 设积分后得 G(y)F(x)C
)
k
7
结束
m
dv dt
mg
kv
v |t 0 0
提示 降落伞所受外力为Fmgkv(k为比例系数)
牛顿第二运动定律Fma 6 下页
例3 设降落伞从跳伞塔 两边积分得
下落后 所受空气阻力与速度 成正比 并设降落伞离开跳伞
dv mg kv
dt m
塔时速度为零 求降落伞下落 即 速度与时间的函数关系
1 k
ln(
dx 设积分后得 G(y)F(x)C
方程由G(y)F(x)C所确定的函数,称为隐式(通)解 2
讨论
微分方程
分离变量
y2xy
y1dy2xdx
3x25xy0
dy(3x25x)dx
(x2y2)dxxydy=0 ————
y1xy2xy2 y(1x)(1y2)
y10xy y x y
yx
10ydy10xdx ————
是否可分离变量 是 是 不是 是 是
不是
3
例 1 求微分方程 dy 2xy 的通解 dx
解 这是一个可分离变量的微分方程
分离变量得
1 y
dy 2xdx
两边积分得
1 dy y
2xdx
即
ln|y|x2C1 注 加常数的另一方法
微分方程PPT(罗兆富等编)第七章 特征线法、达朗贝尔公式和分离变量法
(7.1.05) (7.1.04)
(7.1.06)
所以不妨设
是特征方程组(7.1.03)的一个首次积分.
因为函数a1, a2, … , an 不同时为零, 准形式的常微分方程组
机动 目录 上页
a an a这样特征方程组 an ( x 等价于下面标 1 ( x1 , x2 ,, xn )u x1 0. 2 ( x1 , x2 ,, xn )ux2 (7.1.03) 1 , x2 ,, xn )u xn 0
a1 ( x1 , x2 ,, xn , u )Vx1 a2 ( x1, x2 ,, xn , u)Vxi an ( x1, x2 ,, xn , u)Vxn
b( x1 , x2 ,, xn , u )Vu 0
(7.1.15)
机动 目录 上页 下页 返回
14
结束
a1 ( x1 , x2 ,, xn , u )u x1 a2 ( x1 , x2 ,, xn , u )u x2 an ( x1, x2 ,, xn , u )u xn
1 ( x1 , x2 , , xn ) C1 2 ( x1 , x2 , , xn ) C2 n 1 ( x1 , x2 , , xn ) Cn 1
(7.1.04)
我们的目标是通过求(7.1.03)的首次积分(7.1.04)来求 一阶齐次线性偏微分方程(7.1.01)的通解. 偏微分方程 (7.1.01) (7.1.03)的首 dx1 dx2的解与它的特征方程 dxn (7.1.03) a1 ( x1, x2 ,, xn ) a2 ( x1, x2 ,, xn ) 次积分之间的关系有如下的定理 . an ( x1, x2 ,, xn )
高等数学-第7章 微分方程
将上式两端积分,并由
中的函数可写成的函数,即
(引进新的未知函数(
代入方程(),便得方程
分离变量,得两端积分,得
代替
解方程
因此是齐次方程。
令,则
两端积分,得
以代入上式中的
方程
离变量后得,两端积分,得
,这是对应的齐次线性方程(
把上式代入(
.
以除)的两端,再通过上述代换得线性方程
型的微分方程
(
..
,那末而方程就成为
但是,因此又得到一个一阶微分方程
)的通解为
(3)
合函数的求导法则把化为对
)就成为
通解为
)的通解为
如果函数均是方程的解,那末
我们所求得的解是不是方程的通解呢?
,那末称此两函数在区间,否则,即
如果
就是该方程的通解,其中
的任一特解,
就是方程的通解。
.如果
的解,那末
(
的系数(
和它的各阶导数都只相差一个常数因子。
将
把代入方程(
(
)的两个根。
特征方程微分方程
(
型,
(是与
不是特征方程的根,
若
型
,
,)其中、
)的重复次数。
高等数学 常微分方程PPT课件
第12页/共35页
【解法】需经过变量代换化为一阶线性微分方程.
除方程两边 , 得
yn d y P( x) y1n Q( x) dx
令 z y1n , 则 dz (1 n) yn d y
dx
dx
dz (1 n) P( x) z (1 n)Q( x) (关于z , x的一阶线性方程) dx
特征方程法
待 定
特征方程的根 及其对应项
系
数
法 f(x)的形式及其
特解形式
高阶方程 可降阶方程
线性方程 解的结构
定理1;定理2 定理3;定理4
欧拉方程
第4页/共35页
微分方程解题思路
一阶方程
作 变 换
降 阶
高阶方程
分离变量法 全微分方程 常数变易法
作变换 积分因子
非非 变全 量微 可分
分方 离程
特征方程法
[提示](1)
原方程化为
令u=xy,得 (2) 将方程改写为
d u u ln u (分离变量方程) dx x
d y 1 y y3 (贝努里方程) d x 2x ln x 2x
令 z y2
第17页/共35页
【例3】 识别下列一阶微分方程的类型,并求解
1)
【解】
y y x
①可分离变量的微分方程
u e P( x)d x P( x) ue P( x)d x P( x) u e P( x)d x Q( x)
即 两端积分得
非齐பைடு நூலகம்方程
dy P(x) y Q(x)
dx
u Q(
对应齐次方程通解
x
)
e
P( x)d
y
x
dx
【解法】需经过变量代换化为一阶线性微分方程.
除方程两边 , 得
yn d y P( x) y1n Q( x) dx
令 z y1n , 则 dz (1 n) yn d y
dx
dx
dz (1 n) P( x) z (1 n)Q( x) (关于z , x的一阶线性方程) dx
特征方程法
待 定
特征方程的根 及其对应项
系
数
法 f(x)的形式及其
特解形式
高阶方程 可降阶方程
线性方程 解的结构
定理1;定理2 定理3;定理4
欧拉方程
第4页/共35页
微分方程解题思路
一阶方程
作 变 换
降 阶
高阶方程
分离变量法 全微分方程 常数变易法
作变换 积分因子
非非 变全 量微 可分
分方 离程
特征方程法
[提示](1)
原方程化为
令u=xy,得 (2) 将方程改写为
d u u ln u (分离变量方程) dx x
d y 1 y y3 (贝努里方程) d x 2x ln x 2x
令 z y2
第17页/共35页
【例3】 识别下列一阶微分方程的类型,并求解
1)
【解】
y y x
①可分离变量的微分方程
u e P( x)d x P( x) ue P( x)d x P( x) u e P( x)d x Q( x)
即 两端积分得
非齐பைடு நூலகம்方程
dy P(x) y Q(x)
dx
u Q(
对应齐次方程通解
x
)
e
P( x)d
y
x
dx
高等数学-第七章 微分方程ppt课件
练习: 求方程 dy ex y 的通解. dx
解法 1 分离变量 e ydy exdx
积分
ey ex C
即
(exC)ey1 0 ( C < 0 )
解法 2 令u x y, 则u 1 y
故有
u 1 eu
积分
1
d
u eu
x
C
(1 eu ) eu 1 eu
du
u ln (1 eu ) x C
解: 设所求曲线方程为 y = y(x) , 则有如下关系式:
dy 2x
①
dx
y x1 2
②
由 ① 得 y 2x dx x2 C (C为任意常数)
由 ② 得 C = 1, 因此所求曲线方程为 y x2 1.
引例2. 列车在平直路上以 20 m s 的速度行驶, 制动时
获得加速度 a 0.4 m s2 , 求制动后列车的运动规律.
解: 设列车在制动后 t 秒行驶了s 米 , 即求 s = s (t) .
已知
d2 dt
s
2
0.4 d
s
s t0 0 , d t
t
0 20
由前一式两次积分, 可得 s 0.2 t 2 C1 t C2
利用后两式可得
C1 20, C2 0
因此所求运动规律为 s 0.2 t 2 20 t
ln y x3 ln C
y Cex3
( C 为任意常数 )
( 此式含分离变量时丢失的解 y = 0 )
x ydx ( x2 1) dy 0
例2. 解初值问题 y(0) 1
解: 分离变量得
dy y
1
x x
2
dx
两边积分得 ln y ln 1 ln C x2 1
求微分方程的解PPT课件
y(0) 1
fun=inline('-2*y+2*x^2+2*x','x','y'); [x,y]=ode23(fun,[0,0.5],1);
注:也可以在 tspan 中指定对求解区间的分割,如:
[x,y]=ode23(fun,[0:0.1:0.5],1); %此时 x=[0:0.1:0.5]
第19页/共23页
solver 为Matlab的ODE求解器(可以是 ode45、ode23、ode113、
ode15s、ode23s、ode23t、ode23tb)
没有一种算法可以有效地解决所有的 ODE 问题,因此MATLAB 提供了多 种ODE求解器,对于不同的ODE,可以调用不同的求解器。
第17页/共23页
Matlab提供的ODE求解器
度均可到 10-3~10-6
ode23t 适度刚性 采用梯形算法
适度刚性情形
ode15s
刚性 多步法;Gear’s 反向数值微 若 ode45 失效时,可
分;精度中等
尝试使用
ode23s 刚性 单步法;2 阶Rosebrock 算 当精度较低时,计算时
法;低精度
间比 ode15s 短
ode23tb 刚性 梯形算法;低精度
x y
|t 0 |t 0
1 0
[x,y]=dsolve('Dx+5*x+y=exp(t)','Dy-x-3*y=0', ... 'x(0)=1', 'y(0)=0', 't')
ezplot(x,y,[0,1.3]);
注:解微分方程组时,如果所给的输出个数与方程个数相同,则方程组的解按词 典顺序输出;如果只给一个输出,则输出的是一个包含解的结构(structure)类型 的数据。
《微分方程 》课件
总结词
需要选择合适的代换变量。
详细描述
在使用变量代换法时,需要选择合适的代换变量,使得微 分方程能够被转化为更简单的形式。这个过程需要一定的 技巧和经验。
积分因子法
总结词
通过寻找积分因子,将微分方程转化为积分方程。
详细描述
积分因子法是通过寻找积分因子,将微分方程转化为积 分方程,从而简化求解过程。这种方法适用于具有特定 形式的一阶非线性微分方程。
总结词
通过引入新的变量代换,简化微分方程的形式。
详细描述
变量代换法是通过引入新的变量代换,将微分方程转化为 更简单的形式,从而简化求解过程。这种方法适用于具有 特定形式的高阶微分方程。
总结词
适用于高阶微分方程。
详细描述
变量代换法主要适用于高阶微分方程,通过引入新的变量 代换,可以将高阶微分方程转化为更简单的形式,从而简 化求解过程。
解法
通常需要使用迭代法、级数法或摄动法等非线性 求解方法。
3
特例
当 p(x,y,y') = 0, q(x,y,y') = a(常数)时,方程 简化为 y'' + ay = f(x),其解法与二阶线性微分 方程类似。
二阶常系数线性微分方程
定义
形如 y'' + ay' + by = f(x) 的微分方程称为二阶常系数线性 微分方程。
《微分方程》PPT课件
目 录
• 微分方程简介 • 一阶微分方程 • 二阶微分方程 • 高阶微分方程 • 微分方程的解法 • 微分方程的应用实例
01
微分方程简介
微分方程的定义
总结词
微分方程是描述数学模型中变量之间 动态关系的方程,通过微分来描述函 数的变化率。
需要选择合适的代换变量。
详细描述
在使用变量代换法时,需要选择合适的代换变量,使得微 分方程能够被转化为更简单的形式。这个过程需要一定的 技巧和经验。
积分因子法
总结词
通过寻找积分因子,将微分方程转化为积分方程。
详细描述
积分因子法是通过寻找积分因子,将微分方程转化为积 分方程,从而简化求解过程。这种方法适用于具有特定 形式的一阶非线性微分方程。
总结词
通过引入新的变量代换,简化微分方程的形式。
详细描述
变量代换法是通过引入新的变量代换,将微分方程转化为 更简单的形式,从而简化求解过程。这种方法适用于具有 特定形式的高阶微分方程。
总结词
适用于高阶微分方程。
详细描述
变量代换法主要适用于高阶微分方程,通过引入新的变量 代换,可以将高阶微分方程转化为更简单的形式,从而简 化求解过程。
解法
通常需要使用迭代法、级数法或摄动法等非线性 求解方法。
3
特例
当 p(x,y,y') = 0, q(x,y,y') = a(常数)时,方程 简化为 y'' + ay = f(x),其解法与二阶线性微分 方程类似。
二阶常系数线性微分方程
定义
形如 y'' + ay' + by = f(x) 的微分方程称为二阶常系数线性 微分方程。
《微分方程》PPT课件
目 录
• 微分方程简介 • 一阶微分方程 • 二阶微分方程 • 高阶微分方程 • 微分方程的解法 • 微分方程的应用实例
01
微分方程简介
微分方程的定义
总结词
微分方程是描述数学模型中变量之间 动态关系的方程,通过微分来描述函 数的变化率。
《微分方程复习大纲》课件
《微分方程复习大纲》 PPT课件
欢迎来到我们的《微分方程复习大纲》PPT课件!在这个课件中,我们将一起 探索微分方程的定义、分类,以及解决这些方程的方法。
微分方程的定义和概念
1 微分方程的含义
解决自然现象中变化和关系的数学方程。包含未知函数及其导数的方程。
2 微分方程的分类
分为常微分方程(ODEs)和偏微分方程(PDEs),根据未知函数的变量类型进行分类。
常微分方程的基本解法
分离变量法
将未知函数分为多个表达式并 分别性代数问 题,利用特征根求解。
常系数线性微分方程 法
利用特征根解常系数线性微分 方程,得到通解。
常微分方程的数值解法
1
欧拉法
使用差分代替微分,逐步逼近微分方程的解。
2
龙格-库塔法
通过多次计算,提高数值解的精确度。
将解函数表示为傅立叶级数,逐步逼近方程的数 值解。
应用实例和习题练习
物理学
模拟物体的运动、热传导、波动等现象。
经济学
预测经济发展、市场价格波动等。
工程学
分析电路、热传导、结构稳定性等问题。
数学建模
挑战各种实际问题,加深对微分方程的理解。
3
改进的欧拉法
控制步长大小,并提供更精确的数值解。
偏微分方程的基本解法
热方程
描述物体温度分布随时间变化的 方程。
波动方程
描述波的传播和震荡的方程。
拉普拉斯方程
描述势场的分布和形状的方程。
偏微分方程的数值解法
有限差分法 有限元法 谱方法
将偏微分方程转化为差分表达式,并逐步计算数 值解。
将解域划分为有限个单元,利用逼近函数计算数 值解。
欢迎来到我们的《微分方程复习大纲》PPT课件!在这个课件中,我们将一起 探索微分方程的定义、分类,以及解决这些方程的方法。
微分方程的定义和概念
1 微分方程的含义
解决自然现象中变化和关系的数学方程。包含未知函数及其导数的方程。
2 微分方程的分类
分为常微分方程(ODEs)和偏微分方程(PDEs),根据未知函数的变量类型进行分类。
常微分方程的基本解法
分离变量法
将未知函数分为多个表达式并 分别性代数问 题,利用特征根求解。
常系数线性微分方程 法
利用特征根解常系数线性微分 方程,得到通解。
常微分方程的数值解法
1
欧拉法
使用差分代替微分,逐步逼近微分方程的解。
2
龙格-库塔法
通过多次计算,提高数值解的精确度。
将解函数表示为傅立叶级数,逐步逼近方程的数 值解。
应用实例和习题练习
物理学
模拟物体的运动、热传导、波动等现象。
经济学
预测经济发展、市场价格波动等。
工程学
分析电路、热传导、结构稳定性等问题。
数学建模
挑战各种实际问题,加深对微分方程的理解。
3
改进的欧拉法
控制步长大小,并提供更精确的数值解。
偏微分方程的基本解法
热方程
描述物体温度分布随时间变化的 方程。
波动方程
描述波的传播和震荡的方程。
拉普拉斯方程
描述势场的分布和形状的方程。
偏微分方程的数值解法
有限差分法 有限元法 谱方法
将偏微分方程转化为差分表达式,并逐步计算数 值解。
将解域划分为有限个单元,利用逼近函数计算数 值解。
微分方程PPT课件
x
ln u 1 3 ln(u 2) u 2 1 ln u ln x ln C ,
2
2
u1 3 Cx.
u(u 2)2
微分方程的解为 ( y x)2 Cy( y 2 x)3 . 20
三. 一阶线性微分方程
一阶线性微分方程的标准形式:
dy P( x) y Q( x) (1) dx
初值问题: 求微分方程满足初始条件的解的问题.
一阶:
y f (x, y),
y(
x0
)
y0
.
过定点的积分曲线;
y f ( x, y, y),
过定点且在定点的切线
二阶:
y(
x0
)
y0 ,
y(x0 )
y0
.
的斜率为定值的积分曲线.
n
阶:
f (x, y, y( x0 )
y, y y0 , y(
x A, dx 0,
t 0
dt t0
C1 A,
而
dx dt
kC1
s in kt
kC2
cos kt,
C2 0.
所求特解为 x Acoskt.
9
注意: 1. 有些方程可能无解.
( y)2 y2 1 0 无实函数解.
2. 方程可能有解而无通解. ( y)2 y2 0 只有特解 y 0 . 3. 通解不一定能包含所有的解.
y 3e2x 中不含任意常数,
故为微分方程的特解.
11
6.2 一阶微分方程
一阶微分方程的一般形式是
F( x, y, y) 0
如果一阶导数可解出,则可写为
dy f ( x, y), dx 或 P( x, y)dx Q( x, y)dy 0
微分方程ppt
VIP时长期间,下载特权不清零。
100W优质文档免费下 载
VIP有效期内的用户可以免费下载VIP免费文档,不消耗下载特权,非会员用户需要消耗下载券/积分获取。
部分付费文档八折起 VIP用户在购买精选付费文档时可享受8折优惠,省上加省;参与折扣的付费文档均会在阅读页标识出折扣价格。
服务特 权
共享文档下载特权
VIP用户有效期内可使用共享文档下载特权下载任意下载券标价的文档(不含付费文档和VIP专享文档),每下载一篇共享文
档消耗一个共享文档下载特权。
年VIP
月VIP
连续包月VIP
享受100次共享文档下载特权,一次 发放,全年内有效
赠每的送次VI的发P类共放型享的决文特定档权。下有载效特期权为自1个VI月P,生发效起放每数量月发由放您一购次买,赠 V不 我I送 清 的P生每 零 设效月 。 置起自 随1每5动 时次月续 取共发费 消享放, 。文一前档次往下,我载持的特续账权有号,效-自
服务特 权
共享文档下载特权
VIP用户有效期内可使用共享文档下载特权下载任意下载券标价的文档(不含付费文档和VIP专享文档),每下载一篇共享文
档消耗一个共享文档下载特权。
年VIP
月VIP
连续包月VIP
享受100次共享文档下载特权,一次 发放,全年内有效
赠每的送次VI的发P类共放型享的决文特定档权。下有载效特期权为自1个V月IP,生发效放起数每量月由发您放购一买次,赠 V不 我IP送 清的生每 零设效月 。置起1自随每5次动时月共续取发享费消放文,。一档前次下往,载我持特的续权账有,号效自-
分 方 程
z z xy z2 x y
zx 5z4 0
常微分方程
偏微分方程
100W优质文档免费下 载
VIP有效期内的用户可以免费下载VIP免费文档,不消耗下载特权,非会员用户需要消耗下载券/积分获取。
部分付费文档八折起 VIP用户在购买精选付费文档时可享受8折优惠,省上加省;参与折扣的付费文档均会在阅读页标识出折扣价格。
服务特 权
共享文档下载特权
VIP用户有效期内可使用共享文档下载特权下载任意下载券标价的文档(不含付费文档和VIP专享文档),每下载一篇共享文
档消耗一个共享文档下载特权。
年VIP
月VIP
连续包月VIP
享受100次共享文档下载特权,一次 发放,全年内有效
赠每的送次VI的发P类共放型享的决文特定档权。下有载效特期权为自1个VI月P,生发效起放每数量月发由放您一购次买,赠 V不 我I送 清 的P生每 零 设效月 。 置起自 随1每5动 时次月续 取共发费 消享放, 。文一前档次往下,我载持的特续账权有号,效-自
服务特 权
共享文档下载特权
VIP用户有效期内可使用共享文档下载特权下载任意下载券标价的文档(不含付费文档和VIP专享文档),每下载一篇共享文
档消耗一个共享文档下载特权。
年VIP
月VIP
连续包月VIP
享受100次共享文档下载特权,一次 发放,全年内有效
赠每的送次VI的发P类共放型享的决文特定档权。下有载效特期权为自1个V月IP,生发效放起数每量月由发您放购一买次,赠 V不 我IP送 清的生每 零设效月 。置起1自随每5次动时月共续取发享费消放文,。一档前次下往,载我持特的续权账有,号效自-
分 方 程
z z xy z2 x y
zx 5z4 0
常微分方程
偏微分方程
微分方程
1)首先建立M-文件 (weif.m) function f = weif(x,y) f=-y+x+1;
2)求解:[x,y]=ode23(‘weif’, [0, 1], 1) 3) 作图形: plot(x, y, ‘r’); 4) 与精确解进行比较
hold on ezplot(‘x+exp(-x)’,[0, 1])
数学实验之 --微分方程
实验目的 应用场景 基本概念 数值解法 软件求解 范例 课堂延伸 布置实验
结束
Matlab软件计算数值解
注意1:
dx(t) dt
f1(t, x(t), y(t))
dy(t) dt
f2 (t, x(t), y(t))
1、建立M文件函数
function xdot = fun(t,z) xdot = [f1(t, z(1), z(2)); f2(t, z(1), z(2))];
函数的 初值
函数文 终值
件
ode23:2阶3级龙格-库塔算法 ode45:4阶5级龙格-库塔算法
数学实验之 --微分方程
实验目的 应用场景 基本概念 数值解法 软件求解 范例 课堂延伸 布置实验
结束
Matlab软件计算数值解
例1 y’= - y+x+1,y(0) = 1
标准形式: y’= f(x , y)
实验目的 应用场景 基本概念 数值解法 软件求解 范例 课堂延伸 布置实验
结束
微分方程的数值解法
“常微分方程初值问题数值解”的提法
而在一系列离散点 x1 x2 xn
求y(xn)的近似值yn(n=1,2,…)
y
通常取等步长h
xn x0 nh
2)求解:[x,y]=ode23(‘weif’, [0, 1], 1) 3) 作图形: plot(x, y, ‘r’); 4) 与精确解进行比较
hold on ezplot(‘x+exp(-x)’,[0, 1])
数学实验之 --微分方程
实验目的 应用场景 基本概念 数值解法 软件求解 范例 课堂延伸 布置实验
结束
Matlab软件计算数值解
注意1:
dx(t) dt
f1(t, x(t), y(t))
dy(t) dt
f2 (t, x(t), y(t))
1、建立M文件函数
function xdot = fun(t,z) xdot = [f1(t, z(1), z(2)); f2(t, z(1), z(2))];
函数的 初值
函数文 终值
件
ode23:2阶3级龙格-库塔算法 ode45:4阶5级龙格-库塔算法
数学实验之 --微分方程
实验目的 应用场景 基本概念 数值解法 软件求解 范例 课堂延伸 布置实验
结束
Matlab软件计算数值解
例1 y’= - y+x+1,y(0) = 1
标准形式: y’= f(x , y)
实验目的 应用场景 基本概念 数值解法 软件求解 范例 课堂延伸 布置实验
结束
微分方程的数值解法
“常微分方程初值问题数值解”的提法
而在一系列离散点 x1 x2 xn
求y(xn)的近似值yn(n=1,2,…)
y
通常取等步长h
xn x0 nh
微分方程ppt课件
F(x, y, y) 0
(1.8)
如果在(1.8)中能将 y 解出,则得到方程
y f (x, y)
(1.9)
或
M (x, y)dx N(x, y)dy 0
(1.10)
(1.8)称为一阶隐式方程,(1.9)称为一阶显式方程,(1.10)称为微 分形式的一阶方程.
14
机动 目录 上页 下页 返回 结束
推得
c1 v0
c2 H
于是,得到满足上述初值条件的特解为
xx(t()t)H12gt122 gt2c1t v0ct 2
(1.14)
22
机动 目录 上页 下页 返回 结束
它描述了初始高度为H,初始速度为v0的自由落体运 动规律.
求微分方程满足初值条件的解的问题称为初值 问题.
于是我们称(1.14)是初值问题
4
机动 目录 上页 下页 返回 结束
目
录
第一章 初等积方法 第二章 基本定理 第三章 一阶线性微分方程组 第四章 n阶线性微分方程 第五章 定性与稳定性理论简介 第六章 一阶偏微分方程初步
5
机动 目录 上页 下页 返回 结束
第一讲
第一章 初等积分法
1.1 微分方程和解
300多年前,由牛顿(Newton,1642-1727)和 莱布尼兹(Leibniz,1646-1716)所创立的微积分学, 是人类科学史上划时代的重大发现,而微积分 的产生和发展,又与求解微分方程问题密切相 关.这是因为,微积分产生的一个重要动因来自 于人们探求物质世界运动规律的需求.
12
机动 目录 上页 下页 返回 结束
例如下面的方程都是常微分方程
dy 2x dx
(1.4)
求解微分方程ppt课件
dM dt,
M
dt
ln M t C1 , M Ce t .
由条件 M(0) M0 得 C M0,所以 M (t ) M0e t . M
M0
铀的衰变规律
t
14
例. 解初值问题
xydx ( x2 1) dy 0 y(0) 1
解: 分离变量得
dy y
1
x x2
dx
两边积分得
即
y x2 1 C ( C 为任意常数 )
Q(x) 0时称为一阶齐次线性微分方程。
dy P( x) y 0 dx
(2)
叫做对应于非齐次线性方程(1)的齐次线性方程。
23
一阶齐次线性方程(2)的解法
分离变量法
dy P( x)dx , y
ln y P( x)dx C1,
得方程(2)的通解 y Ce P( x)dx ,
(3)
这里 P( x)dx 表示P(x)的任一原函数。
第七章 微分方程
§ 1 微分方程的基本概念
1
例 一曲线通过点(1, 2),且曲线上任意点切线
的斜率均等于切点横坐标的2倍 ,求这曲线的 方程。
例 列车在平直线路上以 20m /s 的速度行驶,
制动时列车获得加速度 0.4m /s2 。问开始制 动到停止需多少时间?这段时间列车又走了 多远?
2
微分方程的定义
注:通解(3)包含了方程(2)的全部解。
24
一阶非齐次线性方程(1)的解法
常数变易法,令
y(x) C(x)e P(x)dx.
y e P( x)dx
Q(
x)e
P
(
x
)dx
dx
C
25
常微分方程ppt
例 2.4.5 求解方程
解:仔细观察该方程的特征:
对方程做恒等变形得,
自然做变化
原方程化为:
求解上面的线性方程得:
Riccati方程
定义
形如
的方程称为Riccati方程。 一般情况下,Riccati方程无法用初等积分积分 法求出其解,只是对一些特殊情况,或事先知道了 他的一个特解,才可以求出他的通解。
时,方程可通过适当变化化为变量可分离的方程。
证明: 不妨设 化为
否则可通过变量变化 因此,代替原方程,我们考虑 (2.4.5)
当
时, 上述方程是一个变量可分离的方程
当
时,做变量变化
代入原方程得
这是一个变量可分离的方程。
当
时,做变量变化
代入原方程得: (2.4.6)
其中
再做变换
进一步可把方程变为: (2.4.7) 其中
Riccati方程一些可求解的特殊类型:
1、当 都是常数时, Riccati方程
是变量可分离的方程,可以用分离变量法求解。 2、当 3、当 时, Riccati方程是线性方程。 时, Riccati方程是Bernoulli方程。
4、当Riccati方程的形式为:
时,可利用变量替换 可分离的方程。
引进变量
,则
原方程可化为
这是一个变量可分离的方程。
例 2.4.1 求方程 (2.4.1) 解: 令 则
代
代入原变量得到(2.4.3)的通解为:
其它变化法
利用变量替换法求解微分方程十分灵活,一 般依赖于方程的形式和求导的经验。 例 2.4.2 求方程 (2.4.2) 解:将此方程改写为:
将方程化为变量
5、当Riccati方程有一个特解, 时,可利用变量替换 代入原方程得
《数学分析微分方程》课件
III. 高阶微分方程的解法
特征方程法
将高阶齐次微分方程转化为特征方程,通过解特征方程得到齐次部分的解。
待定系数法
假设解为某些未知函数,代入原方程得到待定系数,通过求导和代入原方程求解未知函数。
常数变易法
假设解为常数的函数,通过求导和代入原方程得到常数的解。
IV. 常系数线性微分方程的解法
特征根法
2
方程得到常数的解。
假设解为某些未知函数,代入原方程得到
待定系数,通过求导和代入原方程求解未
知函数。
3
求解自由项
通过求解无齐次项情况下的特解,再加上 通解,得到非齐次线性微分方程的解。
VI. 傅里叶级数方法
傅里叶级数方法可以将周期函数表示成正弦和余弦函数的无穷级数,通过求解系数得到函数的展开式。
VII. 拉普拉斯变换方法
通过求解特征方程的根,得到齐 次线性微分方程的通解。
待定系数法
假设解为某些未知函数,代入原 方程得到待定系数,通过求导和 代入原方程求解未知函数。
常数变易法
假设解为常数的函数,通过求导 和代入原方程得到常数的解。
V. 变系数线性微分方程的解法
1
常数变易法
假设解为常数的函数,通过求导和代入原
待定系数法
《数学分析微分方程》 PPT课件
欢迎来到《数学分析微分方程》PPT课件。本课件将深入介绍微分方程的基本 概念,并详细讲解一阶、高阶、常系数线性、变系数线性微分方程的解法, 以及傅里叶级数和拉普拉斯变换方法的应用。
I. 介绍微分方程的基本概念
学习微分方程前,我们先了解微分方程的基本概念和意义,掌握微分方程的 分类和形式,并探讨微分方程在实际问题中的应用。
拉普拉斯变换方法是一种将时间域函数转换为复频域函数的方法,通过求解 拉普拉斯变换的积分得到函数的解析表达式。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(7.1.01)
5
下页 返回 结束
dx1 a1 ( x1, x2 ,, xn ) dx a ( x , x ,, x ) n 1 2 n n dx2 a2 ( x1 , x2 ,, xn ) dxn an ( x1 , x2 ,, xn ) dxn 1 an 1 ( x1 , x2 ,, xn ) dx a ( x , x ,, x ) n 1 2 n n
x at C
首次积分!
所以方程的通解为
p( x, t ) ( x at ).
其中 是任意连续可微一元函数.
再注意到初始条件p(x, 0)=f(x), 得 ( x) f ( x), 从而得
到方程的解为
p f ( x at ).
■
机动 目录 上页 下页 返回
11
结束
例3. 用特征线法求解一阶齐是任意连续可微二元函数.
9
机动 目录 上页 下页 返回 结束
例1. 用特征线法求解一阶齐次线性偏微分方程
3 yux 2 xu y 0
解: 根据前面的讨论, 写出特征方程组
dx dy 3 y 2 x
3 2 x y C 首次积分! 2
2
所以方程的通解为
3 2 u( x, y) ( x y ). 2
u [1 ( x1 , x2 ,n , xn ), 2 ( x1 , x2 , , xn ), , n 1 ( x1 , x2 , , xn )]
i 1 i 其中 是任意连续可微n1元函数 .
a ( x , x ,, x ) x
i 1 2 n
0
比较
(7.1.08) ■
a1 ( x1 , x2 ,, xn , u )Vx1 a2 ( x1, x2 ,, xn , u)Vxi an ( x1, x2 ,, xn , u)Vxn
b( x1 , x2 ,, xn , u )Vu 0
(7.1.15)
机动 目录 上页 下页 返回
14
结束
a1 ( x1 , x2 ,, xn , u )u x1 a2 ( x1 , x2 ,, xn , u )u x2 an ( x1, x2 ,, xn , u )u xn
(7.1.07)
因此(7.1.06)也是(7.1.07)的一个首次积分. 再由第三章第 n 1 ai 一节定理3.1知, 有恒等式 0, 两端乘以an, 得
ai ( x1, x2 ,, xn ) 0 xi i 1
n
xn
i 1
an xi
(7.1.08)
这就证明了函数 ( x1 , x2 , , xn )是特征方程组(7.1.03)的一 个首次积分的充要条件为恒等式(7.1.08)成立.
1 ( x1 , x2 , , xn ) C1 2 ( x1 , x2 , , xn ) C2 n 1 ( x1 , x2 , , xn ) Cn 1
(7.1.04)
我们的目标是通过求(7.1.03)的首次积分(7.1.04)来求 一阶齐次线性偏微分方程(7.1.01)的通解. 偏微分方程 (7.1.01) (7.1.03)的首 dx1 dx2的解与它的特征方程 dxn (7.1.03) a1 ( x1, x2 ,, xn ) a2 ( x1, x2 ,, xn ) 次积分之间的关系有如下的定理 . an ( x1, x2 ,, xn )
a1 ( x1 , x2 ,, xn , u )u x1 a2 ( x1 , x2 ,, xn , u )u x2 an ( x1, x2 ,, xn , u )u xn
b( x1 , x2 ,, xn , u )
(7.1.13)
由(7.1.15)可见, 若将V视为关于x1, x2 , … , xn, u的函 数,(7.1.15)就成为关于未知函数V的一阶齐次线性偏微分 方程. 这就证明了,若V(x1, x2 , … , xn, u)=C是一阶非齐次 拟线性偏微分方程(7.1.13)的一个隐函数形式的解,则n+1 元函数 V(x1, x2 , … , xn, u) 是一阶齐次线性偏微分方程 (7.1.15)的解.
一、一阶(拟)线性偏微分方程的通解 1. 一阶齐次线性偏微分方程 考虑形如
a1 ( x1 , x2 ,, xn )ux1 a2 ( x1 , x2 ,, xn )ux2 an ( x1 , x2 ,, xn )uxn 0
(7.1.01) 设u=u(x1, x2, … , xn)是方程(7.1.01)的一个解,则由全微 分法则, 有 du ux dx1 ux dx2 ux dxn (7.1.02)
则方程(7.1.09)的通解为
u ( x, y) [ ( x, y)]
(7.1.10)
其中 是任意连续可微一元函数.
8
机动 目录 上页 下页 返回 结束
注: 当n=3时, 方程(7.1.01)成为
a( x, y, z )ux b( x, y, z)uy c( x, y, z)uz 0
i i
Vu
a1 ( x1 , x2 ,, xn , u )Vx1 a2 ( x1, x2 ,, xn , u)Vxi an ( x1, x2 ,, xn , u)Vxn
b( x1 , x2 ,, xn , u )Vu 0
(7.1.15)
机动 目录 上页 下页 返回
13
结束
u ai ( x1, x2 ,, xn ) 0 xi i 1
n
机动 目录
(7.1.01)
7
上页 下页 返回 结束
注: 当n=2时, 方程(7.1.01)成为
a( x, y )ux b(x , y )u y 0
(7.1.09)
dx dy ,它有一个首次积分 ( x, y ) C , 其特征方程组为 a( x, y) b( x, y)
第七章 特征线法、达朗贝尔公式 和分离变量法
第一节 特征线法 第二节 达朗贝尔公式 反射法 第三节 分离变量法简介
1
机动 目录 上页 下页 返回 结束
第一节 特征线法
一、一阶(拟)线性偏微分方程的通解 1. 一阶齐次线性偏微分方程 考虑形如
a1 ( x1 , x2 ,, xn )ux1 a2 ( x1 , x2 ,, xn )ux2 an ( x1 , x2 ,, xn )uxn 0
b( x1 , x2 ,, xn , u )
(7.1.13)
反过来, 假设n+1元函数V(x1, x2 , … , xn, u)是(7.1.15) 的解, 且Vu≠0, 则由(7.1.15)和(7.1.14)可以推出由方程
V x1 , x2 ,, xn , u 0
b( x1 , x2 ,, xn , u )
(7.1.13)
的一阶齐次拟线性偏微分方程的通解, 其中ai(i=1,2,…,n), b都是n+1个变元x1, x2 , … , xn, u的连续函数,且不全为零.
设V(x1, x2 , … , xn, u)=0是方程(7.1.13)的一个隐函数 形式的解, 注意到u是x1, x2 , … , xn的函数,由隐函数求导 法, 得到 Vx ux ,(i 1,2,, n) (7.1.14)
2
其中 是任意连续可微一元函数.
■
10
机动 目录 上页 下页 返回 结束
例2. 求解交通流线性关系模型
p p a 0, t 0, x x t p ( x,0) f ( x).
解: 根据前面的讨论, 写出特征方程组
dt dx 1 a
6
机动 目录 上页 下页 返回 结束
( x1 , x2 , , xn ) 是特征方程组(7.1.03)的一个首次积分
的充要条件是: u ( x1 , x2 ,, xn ) 是一阶齐次线性偏微分 方程(7.1.01)的解.
因此, 若u ( x1 , x2 ,, xn ) 是一阶齐次线性偏微分方 程(7.1.01)的任意一个解,则它是特征方程组(7.1.03)的一 个首次积分. 再由第三章第一节定理3.5,它可由特征方程组(7.1.03) 的n1个首次积分(7.1.04)来表达
1 2 n
dx1 dx2 dxn a1 ( x1, x2 ,, xn ) a2 ( x1, x2 ,, xn ) an ( x1, x2 ,, xn )
机动 目录 上页
(7.1.03)
3
下页 返回 结束
我们称(7.1.03)为(7.1.01)的特征方程组,由特征方程组 (7.1.03)确定的空间曲线称为特征曲线. 由于特征方程组 (7.1.03)是一个包含n1个方程的常微分方程组, 所以它 有n1个首次积分
(7.1.01) 的一阶齐次线性偏微分方程的通解, 其中ai(i=1,2,…,n)是 自变量x1 , x2 , … , xn的n(n≥2)元连续函数, 且不全为零. 方程(7.1.01)的通解可通过求解一个常微分方程组而 得到, 通常称这种求解方法为特征线法.
2
机动 目录 上页 下页 返回 结束
第一节 特征线法
4
机动 目录 上页 下页 返回 结束
定理7.1 假设已经得到特征方程组(7.1.03)的n1个 首次积分(7.1.04), 则一阶齐次线性偏微分方程(7.1.01) 的通解为
u ( x1 , x2 ,, xn ) [1 ( x1 , x2 ,, xn ), 2 ( x1, x2 ,, xn ),, n 1 ( x1, x2 ,, xn )]
首次积分!
所以方程的通解为
u( x, y) ( x y z, x 2 y 2 z 2 ).