常系数非齐次线性常微分方程解法之一pdf
- 格式:pdf
- 大小:187.27 KB
- 文档页数:4
关于非齐次线性常系数微分方程特解的微分算子解法的若干示例一、表示符号把某函数对于自变量x 的导数写成D ,即D=dxd 。
例如,函数y 对x 的一阶导数为y dxdy '=,可以表示成Dy ,同理,y ''可以写成2D y ,三阶、四阶….以此类推D1则代表着求积分,如D1x ,就是⎰xdx ,参看复习指导二、 微分方程的表示如果非齐次方程按降阶写成:)x (f y a y a ya y a n 1n )1n (1)n (0=+'+++-- (1)当然,你也可以写成:)x (f y p y p y p y n 1n )1n (1)n (=+'+++-- ,本质都一样,这种形式相当于(1)方程两边同时除以a 0(0≠)。
这里我们以(1)式为准。
用微分子形式表示方程(1):)x (f y a Dy a y D a y D a n 1n 1n 1n 0=++++-- 方程左边把公因子y 提出来:f(x))y a D a D a D (a n 1n 1n 1n 0=++++--上式中,把)a D a Da D (a n 1n 1n 1n0++++-- 看作关于D 的一个函数表达式,表示成F (D )即F (D )=)a D a Da D (a n 1n 1n 1n 0++++--则方程(1)最终可以写成:F (D )y=f (x )三、 相关结论 F (D )kxe=kxe·F (k )甲也可以写成:)F(k ee )D (F 1kxkx=,(分母不为零时),若分母为零,参见指导书表格内的公式证明:F (D )kxe =kxn 1n 1n 1n0)ea D a Da D (a ++++--=)(ea )(ea )(ea )(ea kxn kx1n )1n (kx1)n (kx0+'+++--=kxn kx1n kx1-n 1kxn 0ea kea eka e k a ++++-kxn 1n 1-n 1n0-kx=F (k )kxe甲注意此处方程左右两端的写法,表达的意义是不一样的,左边F (D )是求导,具体来说左边是kxn 1n 1n 1n0)ea D a D a D (a ++++-- ,即)(ea )(e a )(ea )(ea kxn kx1n )1n (kx1)n (kx0+'+++-- ,而方程右边则是)(ekx乘于多项式F (k )其中,左边的带下划线的部分的函数形式与F (D )一样,因此写成F (k )形式,只是字母 是常数k ,而不是求导了,意义也就不同了,它只是个关于k 的多项式了。
n阶常系数非齐次线性微分方程特解的简便解法
n阶常系数非齐次线性微分方程特解的简便解法是一种简单有效的求解n阶常系数非齐次线性微分方程特解的数值解法。
首先,根据给定的n阶非齐次线性微分方程,确定它的一组特权根以及其置换的相应特权向量。
其次,利用以上n项特权向量构造n阶特权伴随矩阵,然后解出该伴随矩阵的方程组,就可以确定该特解的系数基向量和整体解。
最后,使用前面求得的系数基向量和特权根构造出特解,即可得到n阶常系数非齐次线性微分方程特解要求的解。
另外,关于n阶常系数非齐次线性微分方程特解的简便解法有一个重要的常用结论,即当方程组有多个特权根时,特解就是由各自特权向量的乘积组成的。
这一定理可以使解决非齐次线性微分方程特解简便许多,算法的复杂度也降低了很多。
总的来说,n阶常系数非齐次线性微分方程特解的简便解法是一种非常有效、简单易操作的数值求解方法,可以帮助我们更加因材施教、快速有效地确定并获得满足特解要求的解。
常系数非齐次线性微分方程的特解简单解法
经常系数非齐次线性微分方程(Nonhomogeneous Linear Equation with Constant Coefficient,简称NLCC)是数学分析中一类数学模型,应用广泛,有
着丰富的实际应用价值。
其特解的简单解法尤为重要,为解决NLCC特解提供了一
种有效的方法。
特解是NLCC的一种解法思想,即采用分析相应特征根以确定其特解解析式,
同时应用解析法去求解。
首先,根据模型特征系数确定特征方程,而特征方程常可先用单根定理求得特征根;其次,将特征函数形式成一系列有联系的特征方程;最后根据求出的特征根,用解析法求解特征函数的线性组合,即可解出特解。
NLCC的解决思路,在解析类型或者数值类型上均可有效应用,可以存在多种
解决方案和思路。
特解的简单解法是其中的一种,使用概率化的手段,可以有效地减少 NLCC的解法过程、加快实际解答。
特解的简单解法在求解特征方程是,重点
考虑特征根的具体内容,从复数空间分布和解析特征方程开始,可以为求解NLCC
提供有益的支持。
特解的简单解法不仅能够快速准确地求出特解,而且同时避免了NLCC解法过
程中不必要的工作,从而大大加快了解题速度,同时也可以减少解题复杂程度。
其思想由古老,但仍将在当下的数学应用中发挥有效的支撑作用,是值得积极发展的一门技术。
回顾至今,特解的简单解法甚至成为处理NLCC方程的重要组成部分,为后世
学者提供了有用的解决框架。
未来,将有更多的研究面向进一步发展特解简单解法,使其在快速精准地求解NLCC特解方面拥有更强的能力。
常系数线性微分方程复习
一、常系数线性微分方程的形式和名词解释
1. n 阶常系数线性微分方程的标准形式为:
)
(1)1(1)(t f y a y a y a y n n n n =+′+++−−L
其中 a 1,a 2,L ,a n 是常数,f (t )为连续函数
2. n 阶微分方程的含有n 个独立的任意常数的解,叫做一般解(通解)。
3. 微分方程不含任意常数的解,叫做特解。
4. 把微分方程与初始条件合在一起叫做微分方程的初值问题。
初值问题的解是即满足
微分方程又满足初始条件的特解。
二、常系数线性齐次微分方程的解法
01)1(1)(=+′+++−−y a y a y a y n n n n L
其中a 1,a 2,L ,a n 是常数,等号右端自由项为零
1. 求齐次线性微分方程的特征方程(只要将齐次线性微分方程式中的 y (k )换写成 λk ,
k = 0,1,L ,n ,即得其特征方程)。
011
1=++++−−n n n n a a a λλ
λL
2. 求特征方程的根(称为微分方程的特征根)。
3. 求得了方程的 n 个特征根,就可得到微分方程的n 个线性无关的一般解(根的形
式不同,解的形式也不同)。
(1) 特征方程有n 个互异的实根 λ1, λ2 ,L ,λn 。
方程的通解为 t n t t
c c c y n 21e e e
21λλλ+++=L
例 求齐次微分方程032=−′−′′y y y 的通解
特征方程
0322=−−λλ 求出特征方程的根3121=−=λλ
方程的通解 t t
c c y −+=e e
231
(2) 特征方程有n 个实根,但存在重根(设λ0是方程的k 重根)。
方程的通解为
t n t k t k k c c t c t c c y k n 10e e )e (1121λλλ++++++=++−L L
例 求齐次微分方程043=−′′+′′′y y y 的通解
特征方程0432
3
=−+λλ 求出特征方程的根21
321−===λλλ
方程的通解为 t t
t t c c c y 23221e e
e −−++=
(3) n 个特征根中存在复数根的情况(举例说明)
a. 存在1对不重复的复数根 a ± j β ,n -2个互异的实根。
方程的通解为 t n t t t
c c t c t c y n 3e e sin e cos e
321λλ++++=L ββαα
例 求齐次微分方程05832=−′+′′+′′′y y y y 的通解
特征方程058322
3
=−++λλλ
求出特征方程的根2j 12/1321±−===λλλ
方程的通解为 t c t c c y t t t 2sin e 2cos e e
322
/1−−++=
b. 存在2对重复的复数根 a ± j β ,n -4个互异的实根。
方程的通解为
t
n t
t t t t c c t t c t t c t c t c y n 5e
e
sin e cos e sin e cos e 54321λλ++++++=L ββββαααα
例 求齐次微分方程的通解04444)2()3()4()
5(=+′++++y y y y y y
特征方程
0)2)(1(0
4444222
345=++=+++++λλλλλλλ
求出特征方程的根 (二重根)2j 1321±==−=λλλ
方程的通解为
)2sin 2cos (2sin 2cos e 54321t c t c t t c t c c y t ++++=−
三、常系数线性非齐次微分方程的解法
)
(1)1(1)(t f y a y a y a y n n n n =+′+++−−L
其中a 1,a 2,L ,a n 是常数,f (t )为连续函数 解的形式为 )()(t Y t y y +=
其中:
0)(1)1(1)(=+′+++−−y a y a y a y t y n n n n L 性微分方程是方程所对应的齐次线
的通解。
的任一特解。
是非齐次线性微分方程)(t Y
求解步骤:
第1步:求方程对应的常系数线性齐次微分方程的通解(称作自由分量); 第2步:求常系数线性非齐次微分方程的任一个特解(称作强制分量); 第3步:将自由分量与强制分量相加,得到待求微分方程方程的一般解;
第4步:根据初始条件确定一般解中的待定系数,从而得到方程初值问题的解(最终解答)。
求常系数线性非齐次微分方程的一个特解(强制分量),可用待定系数法。
待定系数法:
根据方程等式右端自由项f (t )的函数类型,猜想它的特解是何种函数类型(包括常数),然后将其代入方程来确定所猜的函数中的系数。
例 求方程66)2()3(2
+−=−′−+′′−t t y y t y t 的一个特解。
通过观察可知,c bt at y ++=2
可能是上述方程的一个特解,将其代入方程得
6
6)()2)(2()2)(3(22+−=++−+−+−t t c bt at b at t a t 66)26(622+−=−−+−t t c b a at at
b c a 2,1−==⇒
取 b = 0,则 c = 0,于是 y = t 2 是方程的一个特解
常见函数 f (t ) 所对应的特解函数类型
f (t )(自由项) 特解的函数类型 C (常数) C 1(常数) at
e 特征根)≠α(e at
C at sin ,at cos
特征根)
(≠±+αj at
C at C cos sin 21
k t 11
21+−++++k k k k
C t C t
C t C L
求特解也可用常数变易法,可参考线性微分方程的相关资料。
在求解电路问题时,电路的稳态响应是描述该电路动态响应所对应的微分方程的一个特解。
例 已知:i L (0)=2A ,u C (0)=0,R =50Ω ,L =0.5H , C =100μF 。
求:i L (t ) 。
解
以 i L 为变量列出微分方程
442
2
102102d d 200d d ×=×++L L
L i t i t
i (1) 求通解(自由分量)
020*******=++P λ特征方程
特征根 λ = -100 ± j 100
)100sin(100cos 100sin )(10010021001θ+=+=−−−t Ke t e C t e C t i t t t L 通解
(2) 求特解(强制分量,稳态解)A 1=L i (3) 全解
)100sin(1)(100θ++=−t Ke t i t L 全解
(4) 由初值定积分常数
i L (0+
)=2A , u C (0+
)=0 (已知)
0)0(1)0(1d d 0===++
+C L
L u L u L t i )100cos(100)100sin(100d d 100100θθ+++−=−−t Ke t Ke t
i t t L
⎪⎩⎪⎨⎧=+−→==+→=++ 0cos 100sin 1000
2sin 12)0(0θθθK K dt
di K i L
L o 452==θK 得
0A )45100sin(21)(100≥++=∴−t t e t i t L o
50 V
u C。