(完整版)专题一(二阶常微分方程解法)

合集下载

二阶常系数线性微分方程的解法word版

二阶常系数线性微分方程的解法word版

第八章 讲第四节 二阶常系数线性微分方程一、二阶常系数线形微分方程的概念形如 )(x f qy y p y =+'+'' 1的方程称为二阶常系数线性微分方程.其中p 、q 均为实数,)(x f 为已知的连续函数.如果0)(≡x f ,则方程式 1变成0=+'+''qy y p y 2我们把方程2叫做二阶常系数齐次线性方程,把方程式1叫做二阶常系数非齐次线性方程. 本节我们将讨论其解法.二、二阶常系数齐次线性微分方程1.解的叠加性定理1 如果函数1y 与2y 是式2的两个解, 则2211y C y C y +=也是式2的解,其中21,C C 是任意常数.证明 因为1y 与2y 是方程2的解,所以有0111=+'+''qy y p y 0222=+'+''qy y p y 将2211y C y C y +=代入方程2的左边,得)()()(221122112211y C y C q y C y C p y C y C ++'+'+''+'' =0)()(22221111=+'+''++'+''qy y p y C qy y p y C 所以2211y C y C y +=是方程2的解.定理1说明齐次线性方程的解具有叠加性.叠加起来的解从形式看含有21,C C 两个任意常数,但它不一定是方程式2的通解.2.线性相关、线性无关的概念设,,,,21n y y y 为定义在区间I 内的n 个函数,若存在不全为零的常数,,,,21n k k k 使得当在该区间内有02211≡+++n n y k y k y k , 则称这n个函数在区间I 内线性相关,否则称线性无关.例如 x x 22sin ,cos ,1在实数范围内是线性相关的,因为0sin cos 122≡--x x又如2,,1x x 在任何区间a,b 内是线性无关的,因为在该区间内要使 02321≡++x k x k k必须0321===k k k .对两个函数的情形,若=21y y 常数, 则1y ,2y 线性相关,若≠21y y 常数, 则1y ,2y 线性无关.3.二阶常系数齐次微分方程的解法定理 2 如果1y 与2y 是方程式2的两个线性无关的特解,则212211,(C C y C y C y +=为任意常数是方程式2的通解.例如, 0=+''y y 是二阶齐次线性方程,x y x y cos ,sin 21==是它的两个解,且≠=x y y tan 21常数,即1y ,2y 线性无关, 所以 x C x C y C y C y cos sin 212211+=+=21,C C 是任意常数是方程0=+''y y 的通解.由于指数函数rxe y =r 为常数和它的各阶导数都只差一个常数因子,根据指数函数的这个特点,我们用rx e y =来试着看能否选取适当的常数r ,使rxe y =满足方程2.将rx e y =求导,得 rx rx e r y re y 2,=''='把y y y ''',,代入方程2,得0)(2=++rx eq pr r 因为0≠rx e , 所以只有 02=++q pr r 3只要r 满足方程式3,rx e y =就是方程式2的解.我们把方程式3叫做方程式2的特征方程,特征方程是一个代数方程,其中r r ,2的系数及常数项恰好依次是方程2y y y ,,'''的系数. 特征方程3的两个根为 2422,1q p p r -±-=, 因此方程式2的通解有下列三种不同的情形. (1) 当042>-q p 时,21,r r 是两个不相等的实根. 2421q p p r -+-=,2422q p p r ---= x r x r e y e y 2121,==是方程2的两个特解,并且≠=-x r r e y y )(2121常数,即1y 与2y 线性无关.根据定理2,得方程2的通解为 x r x r e C e C y 2121+=(2) 当042=-q p 时, 21,r r 是两个相等的实根. 221p r r -==,这时只能得到方程2的一个特解x r e y 11=,还需求出另一个解2y ,且≠12y y 常数,设)(12x u y y =, 即 )(12x u e y x r =)2(),(21121211u r u r u e y u r u e y x r x r +'+''=''+'='. 将222,,y y y '''代入方程2, 得 []0)()2(12111=++'++'+''qu u r u p u r u r u e x r 整理,得0])()2([12111=+++'++''u q pr r u p r u e x r由于01≠x r e , 所以 0)()2(1211=+++'++''u q pr r u p r u 因为1r 是特征方程3的二重根, 所以02,01121=+=++p r q pr r从而有 0=''u因为我们只需一个不为常数的解,不妨取x u =,可得到方程2的另一个解 x r xe y 12=.那么,方程2的通解为x r x r xe C e C y 1121+=即 xr e x C C y 1)(21+=.(3) 当042<-q p 时,特征方程3有一对共轭复根 βαβαi r i r -=+=21, 0≠β于是 x i x i e y ey )(2)(1,βαβα-+== 利用欧拉公式 x i x e ix sin cos +=把21,y y 改写为)sin (cos )(1x i x e e e e y x x i x x i ββαβαβα+=⋅==+)sin (cos )(2x i x e e e e y x x i x xi ββαβαβα-=⋅==-- 21,y y 之间成共轭关系,取-1y =x e y y x βαcos )(2121=+, x e y y i y x βαsin )(2121_2=-= 方程2的解具有叠加性,所以-1y ,-2y 还是方程2的解,并且≠==--x x e x e y y x x βββααtan cos sin 12常数,所以方程2的通解为 )sin cos (21x C x C e y x ββα+=综上所述,求二阶常系数线性齐次方程通解的步骤如下:1写出方程2的特征方程02=++q pr r2求特征方程的两个根21,r r3根据21,r r 的不同情形,按下表写出方程2的通解.例1求方程052=+'+''y y y 的通解.解: 所给方程的特征方程为0522=++r ri r i r 21,2121--=+-=所求通解为 )2sin 2cos (21x C x C e y x +=-.例 2 求方程0222=++S dt dS dtS d 满足初始条件2,400-='===t t S S 的特解.解 所给方程的特征方程为0122=++r r121-==r r通解为 te t C C S -+=)(21 将初始条件40==t S 代入,得 41=C ,于是 t e t C S -+=)4(2,对其求导得te t C C S ---=')4(22 将初始条件20-='=t S 代入上式,得 22=C所求特解为t e t S -+=)24(例3求方程032=-'+''y y y 的通解.解 所给方程的特征方程为 0322=-+r r其根为 1,321=-=r r所以原方程的通解为 x x e C e C y 231+=-二、二阶常系数非齐次方程的解法1.解的结构定理3 设*y 是方程1的一个特解,Y 是式1所对应的齐次方程式2的通解,则*+=y Y y 是方程式1的通解.证明 把*+=y Y y 代入方程1的左端:)()()(*++*'+'+*''+''y Y q y Y p y Y=)()(*+*'+*''++'+''qy py y qY Y p Y=)()(0x f x f =+*+=y Y y 使方程1的两端恒等,所以*+=y Y y 是方程1的解. 定理4 设二阶非齐次线性方程1的右端)(x f 是几个函数之和,如 )()(21x f x f qy y p y +=+'+'' 4 而*1y 与*2y 分别是方程 )(1x f qy y p y =+'+''与 )(2x f qy y p y =+'+''的特解,那么**+21y y 就是方程4的特解, 非齐次线性方程1的特解有时可用上述定理来帮助求出.2.)()(x P e x f m x λ=型的解法 )()(x P e x f m x λ=,其中λ为常数,)(x P m 是关于x 的一个m 次多项式. 方程1的右端)(x f 是多项式)(x P m 与指数函数x e λ乘积的导数仍为同一类型函数,因此方程1的特解可能为x e x Q y λ)(=*,其中)(x Q 是某个多项式函数.把 x e x Q y λ)(=*x e x Q x Q y λλ)]()(['+=*'x e x Q x Q x Q y λλλ)]()(2)([2''+'+=*''代入方程1并消去xe λ,得)()()()()2()(2x P x Q q p x Q p x Q m =+++'++''λλλ 5以下分三种不同的情形,分别讨论函数)(x Q 的确定方法:1 若λ不是方程式2的特征方程02=++q pr r 的根, 即02≠++q p λλ,要使式5的两端恒等,可令)(x Q 为另一个m 次多项式)(x Q m :m m m x b x b x b b x Q ++++= 2210)(代入5式,并比较两端关于x 同次幂的系数,就得到关于未知数m b b b ,,,10 的1+m 个方程.联立解方程组可以确定出),,1,0(m i b i =.从而得到所求方程的特解为x m e x Q y λ)(=*2 若λ是特征方程02=++q pr r 的单根, 即02,02≠+=++p q p λλλ,要使式5成立, 则)(x Q '必须要是m 次多项式函数,于是令)()(x xQ x Q m =用同样的方法来确定)(x Q m 的系数),,1,0(m i b i =.3 若λ是特征方程02=++q pr r 的重根,即,02=++q p λλ 02=+p λ.要使5式成立,则)(x Q ''必须是一个m 次多项式,可令)()(2x Q x x Q m =用同样的方法来确定)(x Q m 的系数.综上所述,若方程式1中的x m e x P x f λ)()(=,则式1的特解为x m k e x Q x y λ)(=*其中)(x Q m 是与)(x P m 同次多项式,k 按λ不是特征方程的根,是特征方程的单根或是特征方程的重根依次取0,1或2.例4 求方程x e y y 232-='+''的一个特解.解 )(x f 是x m e x p λ)(型, 且2,3)(-==λx P m对应齐次方程的特征方程为 022=+r r ,特征根根为2,021-==r r . λ=-2是特征方程的单根, 令xe xb y 20-=*,代入原方程解得230-=b故所求特解为 xxe y 223--=* .例5 求方程x e x y y )1(2-='-''的通解.解 先求对应齐次方程02=+'-''y y y 的通解.特征方程为 0122=+-r r , 121==r r齐次方程的通解为 xe x C C Y )(21+=.再求所给方程的特解1)(,1-==x x P m λ由于1=λ是特征方程的二重根,所以x e b ax x y )(2+=*把它代入所给方程,并约去x e 得126-=+x b ax比较系数,得61=a 21-=b于是 xe x x y )216(2-=*所给方程的通解为 x e x x x C C y y y )6121(3221+-+=+=* 3.x B x A x f ϖϖsin cos )(+=型的解法 ,sin cos )(x B x A x f ωω+=其中A 、B 、ω均为常数.此时,方程式1成为x B x A q y p y ωωsin cos +=+'+'' 7这种类型的三角函数的导数,仍属同一类型,因此方程式7的特解*y 也应属同一类型,可以证明式7的特解形式为)sin cos (x b x a x y k ωω+=*其中b a ,为待定常数.k 为一个整数.当ω±i 不是特征方程02=++q pr r 的根, k 取0;当ω±i 不是特征方程02=++q pr r 的根, k 取1;例6 求方程x y y y sin 432=-'+''的一个特解.解 1=ω,ω±i i ±=不是特征方程为0322=-+r r 的根,0=k .因此原方程的特解形式为x b x a y sin cos +=* 于是 x b x a y cos sin +-=*'x b x a y sin cos --=*''将*''*'*y y y ,,代入原方程,得⎩⎨⎧=--=+-442024b a b a 解得 54,52-=-=b a原方程的特解为: x x y sin 54cos 52--=* 例7 求方程x e y y y x sin 32+=-'-''的通解.解 先求对应的齐次方程的通解Y .对应的齐次方程的特征方程为0322=--r r3,121=-=r rx x e C e C Y 321+=-再求非齐次方程的一个特解*y .由于x e x x f -+=2cos 5)(,根据定理4,分别求出方程对应的右端项为,)(1x e x f =x x f sin )(2=的特解*1y 、*2y ,则 **+=*21y y y 是原方程的一个特解.由于1=λ,ω±i i ±=均不是特征方程的根,故特解为)sin cos (21x c x b ae y y y x ++=+=*** 代入原方程,得x e x c b x c b ae x x sin sin )42(cos )24(4=-++--比较系数,得14=-a 024=+c b 142=-c b解之得 51,101,41-==-=c b a . 于是所给方程的一个特解为 x x e y x sin 51cos 10141-+-=* 所以所求方程的通解为x x e e C e C y Y y x x x sin 51cos 10141321-+-+=+=-*.。

(完整版)专题一(二阶常微分方程解法)(可编辑修改word版)

(完整版)专题一(二阶常微分方程解法)(可编辑修改word版)

二阶微分方程:^5+川切纟+0(小%),严)三°时为齐次t/A- dx 2)八丿、时为非齐次二阶常系数齐次线性微分方程及e解法:(*)/+/?/ + qy = 0,其中”,q为常数;求解步骤:1、写出特征方程:⑷F + " + q = 0,其中门r的系数及常数项恰好是(*)式中几Ny的系数2、求出(△)式的两个根m门3、根据八"2的不同1W况,按下表写出(*)式的通解二阶常系数非齐次线性微分方程y + py' + qy =f (x), /?, g为常数/(x) = RPQ)型,为常数:f(X)= e '^lP (x)cos x + P (x)sin jv]型f丿r二阶常系数非齐次线性微分方程的一般形式是)"+卅+ ©=/(%)(1)苴中P是常数。

方程(1)的通解为对应的齐次方程/+p/ + <7y = 0 (2)的通解Y和方程(1)的一个特解y*之和。

即y = r+yt我们已解决了求二阶常系数齐次线性方程通解的问题,所以,我们只需讨论求二阶常系数非齐次线性微分方程的特解的方法Q下而我们只介绍当方程(1)中的/(兀)为如下两种常见形式时求其特解〉'的方法。

由于方程(1)右端函数/(X)是指数函数€ 'X仃加次多项式Pm(x)的乘积,而指数函数与多项式的乘积的导数仍是这类函数,因此,我们推测:方程(1)的特解应为>* = ^以2(x )( 0(尤)是某个次数待;4^的多项式)y* = e^Q(x) + e^Q'(x)ffy* = e ma)+ 200) + 0气x) I 代入方程(1),得厂•[Q"(x) + (2 +p)Q'(x) + (^ + p + q)Q(x)] = e '^ •仇心)消去£ ",得Q"(兀)+ (2 + P )Q©) + (2 + P + q)Q(x) = P…,(x) ⑶讨论2 n1%如果不是特征方程厂+ P 厂+ g = 0 2 + P + q H 0由于P 用(Q 是一个W 次的多项式,欲使(3)的两端恒等,那未2( X )必为一个加 次多项式,设为Qm W = bo 严 + b]*1 +…+ b”一1兀 + bm2 C3\如果是特征方程厂+ pr + q = °的二重根。

二阶常系数线性微分方程的解法word版

二阶常系数线性微分方程的解法word版

第四节 二阶常系数线性微分圆程之阳早格格创做一、二阶常系数线形微分圆程的观念 形如 )(x f qy y p y =+'+'' (1)p 、q 均为真数,)(x f 为已知的连绝函数.如果0)(≡x f ,则圆程式 (1)形成0=+'+''qy y p y (2)咱们把圆程(2)喊干二阶常系数齐次线性圆程,把圆程式(1)喊干二阶常系数非齐次线性圆程. 本节咱们将计划其解法.二、二阶常系数齐次线性微分圆程1.解的叠加性定理1 如果函数1y 与2y 是式(2)的二个解, 则2211y C y C y +=也是式(2)的解,其中21,C C 是任性常数.道明 果为1y 与2y 是圆程(2)的解,所以有 将2211y C y C y +=代进圆程(2)的左边,得=0)()(22221111=+'+''++'+''qy y p y C qy y p y C 所以2211y C y C y +=是圆程(2)的解.定理1道明齐次线性圆程的解具备叠加性. 叠加起去的解从形式瞅含有21,C C 二个任性常数,但是它纷歧定是圆程式(2)的通解.2.线性相闭、线性无闭的观念设,,,,21n y y y 为定义正在区间I 内的n 个函数,若存留没有齐为整的常数,,,,21n k k k 使恰当正在该区间内有02211≡+++n n y k y k y k , 则称那n 个函数正在区间I 内线性相闭,可则称线性无闭.比圆 x x 22sin ,cos ,1正在真数范畴内是线性相闭的,果为又如2,,1x x 正在所有区间(a,b)内是线性无闭的,果为正在该区间内要使必须0321===k k k .对于二个函数的情形,若=21y y 常数, 则1y ,2y 线性相闭,若≠21y y 常数, 则1y ,2y 线性无闭. 3.二阶常系数齐次微分圆程的解法定理2 如果1y 与2y 是圆程式(2)的二个线性无闭的特解,则212211,(C C y C y C y +=为任性常数)是圆程式(2)的通解.比圆,0=+''y y 是二阶齐次线性圆程,x y x y cos ,sin 21==是它的二个解,且≠=x y y tan 21常数,即1y ,2y 线性无闭, 所以 x C x C y C y C y cos sin 212211+=+= (21,C C 是任性常数)是圆程0=+''y y 的通解. 由于指数函数rx e y =(r 为常数)战它的各阶导数皆只好一个常数果子, 根据指数函数的那个特性,咱们用rx e y =去试着瞅是可采用适合的常数r ,使rx e y =谦脚圆程(2).将rx e y =供导,得把y y y ''',,代进圆程(2),得果为0≠rx e , 所以惟有 02=++q pr r (3)只消r 谦脚圆程式(3),rx e y =便是圆程式(2)的解.咱们把圆程式(3)喊干圆程式(2)的特性圆程,特性圆程是一个代数圆程,其中r r ,2的系数及常数项恰佳依次是圆程(2)y y y ,,'''的系数.特性圆程(3)的二个根为 2422,1q p p r -±-=, 果此圆程式(2)的通解有下列三种分歧的情形.(1) 当042>-q p 时,21,r r 是二个没有相等的真根.2421q p p r -+-=,2422q p p r ---= xr x r e y e y 2121,==是圆程(2)的二个特解,而且≠=-x r r e y y )(2121常数,即1y 与2y ,得圆程(2)的通解为x r x r e C e C y 2121+=(2) 当042=-q p 时, 21,r r 是二个相等的真根.221pr r -==,那时只可得到圆程(2)的一个特解x r e y 11=,还需要出另一个解2y ,且≠12y y 常数,设)(12x u y y =, 即 )2(),(21121211u r u r u e y u r u e y x r x r +'+''=''+'='. 将222,,y y y '''代进圆程(2), 得 整治,得由于01≠x r e , 所以 0)()2(1211=+++'++''u q pr r u p r u 果为1r 是特性圆程(3)的二沉根, 所以进而有 0=''u果为咱们只需一个没有为常数的解,无妨与x u =,可得到圆程(2)的另一个解x r xe y 12=.那么,圆程(2)的通解为即 x r e x C C y 1)(21+=. (3) 当042<-q p 时,特性圆程(3)有一对于共轭复根βαβαi r i r -=+=21, (0≠β)于是 x i x i e y e y )(2)(1,βαβα-+==利用欧推公式 x i x e ix sin cos +=把21,y y 改写为 21,y y 之间成共轭闭系,与-1y =x e y y x βαcos )(2121=+, 圆程(2)的解具备叠加性,所以-1y ,-2y 仍旧圆程(2)的解,而且≠==--x x e x e y y x x βββααtan cos sin 12常数,所以圆程(2)的通解为综上所述,供二阶常系数线性齐次圆程通解的步调如下:(1)写出圆程(2)的特性圆程(2)供特性圆程的二个根21,r r(3)根据21,r r 的分歧情形,按下表写出圆程(2)的通解.例1供圆程052=+'+''y y y 的通解.解: 所给圆程的特性圆程为所供通解为 )2sin 2cos (21x C x C e y x +=-.例 2 供圆程0222=++S dt dS dt S d 谦脚初初条件2,400-='===t t S S 的特解.解 所给圆程的特性圆程为通解为 t e t C C S -+=)(21将初初条件40==t S 代进,得 41=C ,于是t e t C S -+=)4(2,对于其供导得将初初条件20-='=t S 代进上式,得所供特解为例3供圆程032=-'+''y y y 的通解.解 所给圆程的特性圆程为 0322=-+r r其根为 1,321=-=r r所以本圆程的通解为 x x e C e C y 231+=-二、二阶常系数非齐次圆程的解法1.解的结构定理3 设*y 是圆程(1)的一个特解,Y 是式(1)所对于应的齐次圆程式(2)的通解,则*+=y Y y 是圆程式(1)的通解.道明 把*+=y Y y 代进圆程(1)的左端:=)()(*+*'+*''++'+''qy py y qY Y p Y=)()(0x f x f =+*+=y Y y 使圆程(1)的二端恒等,所以*+=y Y y 是圆程(1)的解.定理4 设二阶非齐次线性圆程(1)的左端)(x f 是几个函数之战,如)()(21x f x f qy y p y +=+'+''(4)而*1y 与*2y 分别是圆程 )(1x f qy y p y =+'+''与)(2x f qy y p y =+'+''的特解,那么**+21y y 便是圆程(4)的特解, 非齐次线性圆程(1)的特解偶尔可用上述定理去助闲供出.2.)()(x P e x f m x λ=型的解法)()(x P e x f m x λ=,其中λ为常数,)(x P m 是闭于x 的一个m 次多项式.圆程(1)的左端)(x f 是多项式)(x P m 与指数函数x e λ乘积的导数仍为共一典型函数,果此圆程(1)的特解大概为x e x Q y λ)(=*,其中)(x Q 是某个多项式函数. 把 x e x Q y λ)(=*代进圆程(1)并消去x e λ,得)()()()()2()(2x P x Q q p x Q p x Q m =+++'++''λλλ(5) 以下分三种分歧的情形,分别计划函数)(x Q 的决定要领:(1) 若λ没有是圆程式(2)的特性圆程02=++q pr r 的根, 即02≠++q p λλ,要使式(5)的二端恒等,可令)(x Q 为另一个m 次多项式)(x Q m :代进(5)式,并比较二端闭于x 共次幂的系数,便得到闭于已知数m b b b ,,,10 的1+m ),,1,0(m i b i =.进而得到所供圆程的特解为(2) 假如λ特性圆程02=++q pr r 的单根, 即02,02≠+=++p q p λλλ,要使式(5)创造, 则)(x Q '必须假如m 次多项式函数,于是令用共样的要领去决定)(x Q m 的系数),,1,0(m i b i =.(3) 假如λ特性圆程02=++q pr r 的沉根,即,02=++q p λλ02=+p λ.要使(5)式创造,则)(x Q ''必须是一个m 次多项式,可令用共样的要领去决定)(x Q m 的系数.综上所述,若圆程式(1)中的x m e x P x f λ)()(=,则式(1)的特解为其中)(x Q m 是与)(x P m 共次多项式,k 按λ没有是特性圆程的根,是特性圆程的单根或者是特性圆程的沉根依次与0,1或者2.例4 供圆程x e y y 232-='+''的一个特解.解)(x f 是x m e x p λ)(型, 且2,3)(-==λx P m对于应齐次圆程的特性圆程为 022=+r r ,特性根根为2,021-==r r .λ=-2是特性圆程的单根, 令x e xb y 20-=*,代进本圆程解得故所供特解为 x xe y 223--=* . 例5 供圆程x e x y y )1(2-='-''的通解.解 先供对于应齐次圆程02=+'-''y y y 的通解. 特性圆程为 0122=+-r r , 121==r r齐次圆程的通解为 x e x C C Y )(21+=. 再供所给圆程的特解由于1=λ是特性圆程的二沉根,所以把它代进所给圆程,并约去x e 得比较系数,得于是 x e x x y )216(2-=* 所给圆程的通解为 x e x x x C C y y y )6121(3221+-+=+=* 3.x B x A x f ϖϖsin cos )(+=型的解法,sin cos )(x B x A x f ωω+=其中A 、B 、ω均为常数. 此时,圆程式(1)成为x B x A q y p y ωωsin cos +=+'+''(7)那种典型的三角函数的导数,仍属共一典型,果此圆程式(7)的特解*y 也应属共一典型,不妨道明式(7)的特解形式为其中b a ,为待定常数.k 为一个整数. 当ω±i 没有是特性圆程02=++q pr r 的根,k 与0; 当ω±i 没有是特性圆程02=++q pr r 的根,k 与1; 例6 供圆程x y y y sin 432=-'+''的一个特解. 解 1=ω,ω±i i ±=没有是特性圆程为0322=-+r r 的根,0=k .果此本圆程的特解形式为 于是 x b x a y cos sin +-=*'将*''*'*y y y ,,代进本圆程,得解得 54,52-=-=b a本圆程的特解为: x x y sin 54cos 52--=* 例7 供圆程x e y y y x sin 32+=-'-''的通解. 解先供对于应的齐次圆程的通解Y .对于应的齐次圆程的特性圆程为再供非齐次圆程的一个特解*y .由于x e x x f -+=2cos 5)(,根据定理4,分别供出圆程对于应的左端项为,)(1x e x f =x x f sin )(2=的特解*1y 、*2y ,则 **+=*21y y y 是本圆程的一个特解. 由于1=λ,ω±i i ±=均没有是特性圆程的根,故特解为代进本圆程,得比较系数,得解之得 51,101,41-==-=c b a . 于是所给圆程的一个特解为 所以所供圆程的通解为x x e e C e C y Y y x x x sin 51cos 10141321-+-+=+=-*.。

二阶微分方程的解法

二阶微分方程的解法

二阶微分方程的解法二阶微分方程是一种重要的数学工具,使用普通方程难以描述的许多自然现象,可以通过二阶微分方程来描述。

二阶微分方程的解法一般通过分离变量、变量代换、常数变易法、常微分方程定理等多种方法来实现。

1.分离变量法对于形如 y''=f(x)y 的二阶微分方程,可以通过分离变量来解决。

首先将方程转化为 y''/y=f(x),然后对两端同时积分,得到ln|y|=∫f(x)dx+C(常数),则 y=Ae^(∫f(x)dx)或 y=Be^(-∫f(x)dx)。

2.变量代换法当二阶微分方程存在某种特殊的变量代换时,我们可以通过代换来解方程。

例如,对于 y''+p(x)y'+q(x)y=0 的方程,如果我们用y=e^(∫p(x)dx)v(x) 进行代换,则方程转化后的 v(x) 满足 v''+(q(x)-p'(x))v(x)=0,可以进一步使用其他的解法来求解。

3.常数变易法常数变易法主要适用于二阶齐次线性微分方程 y''+p(x)y'+q(x)y=0 的特殊情况。

在解此类方程时,我们常常按照 y=e^(mx) 代入方程,然后解出对应的特征方程。

如果特征方程的根是实数或共轭复数对,那么方程的通解可以表示为y=C1e^(αx)+C2e^(βx),其中 C1,C2 是任意常数,α,β 是特征根;如果特征方程的根是重根,那么方程的通解可以表示为 y=(C1+C2x)e^(mx)。

4.常微分方程定理对于非齐次线性微分方程 y''+p(x)y'+q(x)y=f(x) 的解法,可以利用常微分方程定理(又称为Lagrange公式)来完成。

该定理指出,非齐次线性微分方程的特解可以表示为y*=u(x)y1+v(x)y2,其中 y1,y2分别为解齐次方程 y''+p(x)y'+q(x)y=0,u(x) 和 v(x) 是待定系数函数。

二阶线性常微分方程的级数解法解析课件

二阶线性常微分方程的级数解法解析课件

fn
(s)
sPn
Qn
.
(n 1, 2,
),由于a0 0,必有
f0 (s) s(s 1) sP0 Q0 0 上式为指标方程,其根s1和s2称为正则奇点的指标数.
从而得到方程的一个解w1(z) (z z0 )s1 ak (z z0 )k k 0
求第二个特解
1 s1 s2 整数包括零,则在所设解中取s s2,此时f0 (s2 ) 0,
由于J m
(x)
k 0
k
(1)k !(m
k
1)
( x )m2k,其中m为整数,当 2
k m时, m k 1为负数,函数的值为无穷大,因此对k
求和是从k
m开始,即J m
(x)
k m
k
(1)k !(m
k
1)
( x)m2k 2
令n k m,求和指标从k变到m,则有
Jm (x)
dz2 z dz
z2
在有限远处的奇点为z0 0,且z0 0 是方程的正则奇点.
5.2 方程常点邻域内的解
1.常点邻域内的级数解定理
若p(z)和q(z)在圆形域 | z z0 | R内单值解析,则常微分初值问题
d 2w
dz 2
p(z)
dw dz
q(z)w
0
w(z0 ) a0 , w(z0 ) a1
f0 (s2 k) 0,k 1, 2, 对任选a0 0可唯一确定另外一个解
w2 (z) (z z0 )s2 bk (z z0 )k,w1(z)和w2 (z)线性无关. k 0
2当s1 s2 n 整数,f0 (s2 ) 0,f0 (s2 n) 0,递推到第n步
令a0 a1 an1 0,an 0,可唯一确定ak (k n),从而

二阶常微分方程解

二阶常微分方程解

第七节 二阶常系数线性微分方程的解法在上节我们已经讨论了二阶线性微分方程解的结构,二阶线性微分方程的求解问题,关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解.本节讨论二阶线性方程的一个特殊类型,即二阶常系数线性微分方程及其求解方法.先讨论二阶常系数线性齐次方程的求解方法.§ 二阶常系数线性齐次方程及其求解方法设给定一常系数二阶线性齐次方程为22dx y d +p dxdy +qy =0 其中p 、q 是常数,由上节定理二知,要求方程的通解,只要求出其任意两个线性无关的特解y 1,y 2就可以了,下面讨论这样两个特解的求法.我们先分析方程可能具有什么形式的特解,从方程的形式上来看,它的特点是22dx y d ,dx dy,y 各乘以常数因子后相加等于零,如果能找到一个函数y,其22dx y d ,dxdy ,y 之间只相差一个常数因子,这样的函数有可能是方程的特解,在初等函数中,指数函数e rx,符合上述要求,于是我们令y =e rx其中r 为待定常数来试解将y =e rx,dxdy=re rx,22dx y d =r 2e rx代入方程得 r 2e rx +pre rx +qe rx=0或 e rxr 2+pr +q =0因为e rx≠0,故得r 2+pr +q =0由此可见,若r 是二次方程r 2+pr +q =0的根,那么e rx 就是方程的特解,于是方程的求解问题,就转化为求代数方程的根问题.称式为微分方程的特征方程.特征方程是一个以r 为未知函数的一元二次代数方程.特征方程的两个根r 1,r 2,称为特征根,由代数知识,特征根r 1,r 2有三种可能的情况,下面我们分别进行讨论.1若特证方程有两个不相等的实根r 1,r 2,此时e r 1x ,e r2x 是方程的两个特解.因为 x r xr 21e e =e x)r r (21-≠常数所以e r1x ,e r2x 为线性无关函数,由解的结构定理知,方程的通解为y =C 1e r1x +C 2e r2x2若特征方程有两个相等的实根r 1=r 2,此时p 2-4q =0,即有r 1=r 2=2p-,这样只能得到方程的一个特解y 1=e r 1x,因此,我们还要设法找出另一个满足12y y ≠常数,的特解y 2,故12y y 应是x 的某个函数,设12y y =u,其中u =ux 为待定函数,即 y 2=uy 1=ue r 1x对y 2求一阶,二阶导数得dx dy 2=dxdu e r1x+r 1ue r1x=dx du +r 1uer1x 222dx y d =r 21u +2r 1dx du +22dx ud e r1x将它们代入方程得r 21u +2r 1dx du +22dxu d e r1x+p dxdu +r 1uer1x+que r1x =0或22dx u d +2r 1+p dxdu+r 21+pr 1+que r1x =0因为e r1x ≠0,且因r 1是特征方程的根,故有r 21+pr 1+q =0,又因r 1=-2p故有2r 1+p =0,于是上式成为 22dxu d =0 显然满足22dxud =0的函数很多,我们取其中最简单的一个 ux =x则y 2=xe rx 是方程的另一个特解,且y 1,y 2是两个线性无关的函数,所以方程的通解是y =C 1e r1x +C 2xe r1x =C 1+C 2xe r1x3若特征方程有一对共轭复根 r 1=α+i β,r 2=α-i β此时方程有两个特解y 1=eα+i βxy 2=eα-i βx则通解为y =C 1e α+i βx +C 2e α-i βx其中C 1,C 2为任意常数,但是这种复数形式的解,在应用上不方便.在实际问题中,常常需要实数形式的通解,为此利用欧拉公式e ix =cosx +isinx,e -ix =cosx -isinx有 21e ix+e -ix=cosxi 21e ix-e -ix=sinx21 y 1+y 2=21e αxe i βx+e -i βx=e αxcos βxi 21 y 1-y 2=i21e αxe i βx-e -i βx=e αxsin βx由上节定理一知,21 y 1+y 2,i21y 1-y 2是方程的两个特解,也即eαxcosβx,e αx sin βx 是方程的两个特解:且它们线性无关,由上节定理二知,方程的通解为y =C 1e αx cos βx +C 2e αx sin βx或 y =e αx C 1cos βx +C 2sin βx其中C 1,C 2为任意常数,至此我们已找到了实数形式的通解,其中α,β分别是特征方程复数根的实部和虚部.综上所述,求二阶常系数线性齐次方程的通解,只须先求出其特征方程的根,再根据他的三种情况确定其通解,现列表如下特征方程r 2+pr +q =0的根微分方程22dx y d +p dx dy+qy =0的通解有二个不相等的实根r 1,r 2y =C 1e r1x+C 2e r2x有二重根r 1=r 2y =C 1+C 2xe r1x有一对共轭复根β-α=β+α=i r i r 21y =e αx C 1cos βx +C 2sin βx例1. 求下列二阶常系数线性齐次方程的通解1 22dx y d +3dx dy-10y =0 2 22dx y d -4dx dy +4y =0 3 22dx y d +4dxdy +7y =0 解 1特征方程r 2+3r -10=0有两个不相等的实根r 1=-5,r 2=2所求方程的通解 y =C 1e -5r+C 2e 2x2特征方程r 2-4r +4=0,有两重根 r 1=r 2=2所求方程的通解y =C 1+C 2xe 2x3特征方程r 2+4r +7=0有一对共轭复根r 1=-2+3i r 2=-2-3i所求方程的通解 y =e -2x C 1cos3x +C 2sin 3x§ 二阶常系数线性非齐次方程的解法由上节线性微分方程的结构定理可知,求二阶常系数线性非齐次方程22dx y d +p dxdy +qy =fx 的通解,只要先求出其对应的齐次方程的通解,再求出其一个特解,而后相加就得到非齐次方程的通解,而且对应的齐次方程的通解的解法,前面已经解决,因此下面要解决的问题是求方程的一个特解.方程的特解形式,与方程右边的fx 有关,这里只就fx 的两种常见的形式进行讨论.一、fx =p n xe αx ,其中p n x 是n 次多项式,我们先讨论当α=0时,即当fx =p n x 时方程22dx y d +p dx dy +qy =p nx 的一个特解.1如果q ≠0,我们总可以求得一n 次多项式满足此方程,事实上,可设特解~y =Q nx =a 0x n+a 1xn -1+…+a n,其中a 0,a 1,…a n 是待定常数,将~y 及其导数代入方程,得方程左右两边都是n 次多项式,比较两边x 的同次幂系数,就可确定常数a 0,a 1,…a n .例1. 求22dx y d +dxdy+2y =x 2-3的一个特解. 解 自由项fx =x 2-3是一个二次多项式,又q =2≠0,则可设方程的特解为~y =a 0x 2+a 1x +a 2求导数~'y =2a 0x +a1~"y =2a代入方程有2a 0x 2+2a 0+2a 1x +2a 0+a 1+2a 2=x 2-3比较同次幂系数⎪⎩⎪⎨⎧-=++=+=3a 2a a 20a 2a 21a 2210100 解得 47a 21a 21a 210-=-==所以特解~y =21x 2-21x -472如果q =0,而p ≠0,由于多项式求导一次,其次数要降低一次,此时~y =Q n x 不能满足方程,但它可以被一个n +1次多项式所满足,此时我们可设~y =xQ n x =a 0x n +1+a 1x n +…+a n x代入方程,比较两边系数,就可确定常数a 0,a 1,…a n .例2. 求方程22dx y d +4dxdy=3x 2+2的一个特解. 解 自由项 fx =3x 2+2是一个二次多项式,又q =0,p =4≠0,故设特解~y =a 0x 3+a 1x 2+a 2x求导数~'y =3a 0x 2+2a 1x +a2~"y =6a 0x +2a1代入方程得12a 0x 2+8a 1+6a 0x +2a 1+4a 2=3x 2+2,比较两边同次幂的系数⎪⎩⎪⎨⎧=+=+=2a 4a 20a 6a 83a 1221010 解得 3219a 163a 41a 210=-==所求方程的特解 ~y =41x 3-163x 2+3219x3如果p =0,q =0,则方程变为22dxyd =p nx,此时特解是一个n +2次多项式,可设~y =x 2Q nx,代入方程求得,也可直接通过两次积分求得.下面讨论当α≠0时,即当fx =p n xe αx 时方程22dx y d +p dxdy +qy =p nxe αx的一个特解的求法,方程与方程相比,只是其自由项中多了一个指数函数因子e αx ,如果能通过变量代换将因子e αx 去掉,使得化成式的形式,问题即可解决,为此设y =ue αx ,其中u =ux 是待定函数,对y =ue αx ,求导得dx dy =e αxdxdu+αue αx 求二阶导数 22dx y d =e αx22dx u d +2αe αxdxdu+α2ue αx代入方程得e αx22dx u d +2αdx du +α2u +pe αxdx du +αu +que αx=p n xeαx消去e αx得22dx u d +2α+p dxdu +α2+p α+qu =p nx 由于式与形式一致,于是按的结论有:1如果α2+p α+q ≠0,即α不是特征方程r 2+pr +q =0的根,则可设的特解u =Qn x,从而可设的特解为~y =Q n xe αx2如果α2+p α+q =0,而2α+p ≠0,即α是特征方程r 2+pr +q =0的单根,则可设的特解u =xQ n x,从而可设的特解为~y =xQ n xe αx3如果r 2+p α+q =0,且2α+p =0,此时α是特征方程r 2+pr +q =0的重根,则可设的特解u =x 2Q n x,从而可设的特解为~y =x 2Q n xe αx例3. 求下列方程具有什么样形式的特解122dx y d +5dx dy +6y =e 3x 2 22dx y d +5dx dy +6y =3xe -2x 3 22dx y d +αdxdy +y =-3x 2+1e -x解 1因α=3不是特征方程r 2+5r +6=0的根,故方程具有形如~y =a 0e3x 的特解.2因α=-2是特征方程r 2+5r +6=0的单根,故方程具有形如~y =xa 0x +a 1e -2x的特解.3因α=-1是特征方程r 2+2r +1=0的二重根,所以方程具有形如~y =x 2a 0x 2+a 1x +a 2e -x的特解.例4. 求方程22dxyd +y =x -2e 3x的通解.解 特征方程 r 2+1=0特征根 r =±i 得,对应的齐次方程22dxyd +y =0的通解为 Y =C 1cos x +C 2sin x由于α=3不是特征方程的根,又p n x =x -2为一次多项式,令原方程的特解为~y =a 0x +a 1e 3x此时u =a 0x +a 1,α=3,p =0,q =1,求ux 的导数dxdu =a 0,22dx u d =0,代入22dx u d +2α+p dxdu+α2+αp +qu =x -2得: 10a 0x +10a 1+6a 0=x -2比较两边x 的同次幂的系数有⎩⎨⎧-=+=2a 6a 101a 10010 解得 a 0=101,a 1=-5013于是,得到原方程的一个特解为~y =101x -5013e3x所以原方程的通解是y =Y +~y =C 1cosx +C 2sinx +101x -5013e 3x例5. 求方程22dx y d -2dxdy-3y =x 2+1e -x的通解. 解 特征方程 r 2-2r -3=0特征根 r 1=-1,r 2=3所以原方程对应的齐次方程22dx y d -2dxdy-3y =0的通解Y =C 1e -x +C 2e 3x ,由于α=-1是特征方程的单根,又p n x =x 2+1为二次多项式,令原方程的特解~y =xa 0x 2+a 1x +a 2e -x此时 u =a 0x 3+a 1x 2+a 2x,α=-1,p =-2,q =-3对ux 求导dx du=3a 0x 2+2a 1x +a 222dx ud =6a 0x +2a 1代入22dx u d +2α+p dxdu +α2+pr +qu =x 2+1,得-12a 0x 2+6a 0-8ax +2a 1-4a 2=x 2+1比较x 的同次幂的系数有⎪⎪⎩⎪⎪⎨⎧=--==-0a 8a 6121a 1a 121000 解得 329a 0a 4a 2161a 2011-==--=故所求的非齐次方程的一个特解为~y =-4x 3x 2+4x +89e-x二、fx =p n xe αx cos βx 或p n xe αx sin βx,即求形如22dx y d +p dx dy +qy =p nxe αx cos βx 22dx y d +p dx dy+qy =p nxe αx sin βx 这两种方程的特解.由欧拉公式知道,p n xe αx cos βx,p n xe αx sin x 分别是函数p n xe α+i βx 的实部和虚部.我们先考虑方程22dx y d +p dxdy +qy =p nxe α+i βx方程与方程类型相同,而方程的特解的求法已在前面讨论.由上节定理五知道,方程的特解的实部就是方程的特解,方程的特解的虚部就是方程的特解.因此,只要先求出方程的一个特解,然而取其实部或虚部即可得方程或的一个特解.注意到方程的指数函数e α+i βx 中的α+i ββ≠0是复数,而特征方程是实系数的二次方程,所以α+i β最多只能是它的单根.因此方程的特解形为Q n xeα+i βx或x Qn xeα+i βx.例6. 求方程22dxyd -y =e xcos2x 的通解. 解 特征方程 r 2-1=0特征根 r 1=1,r 2=-1于是原方程对应的齐次方程的通解为Y =C 1e x +C 2e -x为求原方程的一个特解~y .先求方程22dxyd -y =e 1+2ix的一个特解,由于1+2i 不是特征方程的根,且p n x 为零次多项式,故可设u =a 0,此时α=1+2i,p =0,q =-1代入方程22dx u d +2α+p dxdu+α2+αp +qu =1 得1+2i 2-1a 0=1 ,即4i -4a 0=1,得a 0=)1i (41 =-81i +1这样得到22dx y d -y =e 1+2ix的一个特解y =-81i +1e 1+2ix由欧拉公式y =-81i +1e 1+2ix=-81i +1e xcos 2x +isin2x=-81e xcos2x -sin2x +icos2x +sin2x取其实部得原方程的一个特解~y =-81e xcos 2x -sin2x故原方程的通解为y =Y +~y =C 1e x+C 2e-x-81e x cos2x -sin2x 例7. 求方程22dxyd +y =x -2e 3x+xsinx 的通解.解 由上节定理三,定理四,本题的通解只要分别求22dxyd +y =0的特解Y,22dxy d +y =x -2e 3x的一个特解~1y , 22dxy d +y =x sin x 的一个特解~2y 然而相加即可得原方程的通解,由本节例4有Y =C 1cosx +C 2sinx,~1y =101x -5013e3x下面求~2y ,为求~2y 先求方程22dxy d +y =xe ix由于i是特征方程的单根,且pn x=x为一次式,故可设u=xax+a1=a0x2+a1x,此时α=i,p=0,q=1,对u 求导dxdu=2ax+a1,22dxud=2a代入方程22dxud+2α+pdxdu+α2+pα+qu=x得 2a0+2i2ax+a1+0=x即 4iax+2ia1+2a=x比较x的同次幂的系数有:⎩⎨⎧=+=a2ia21ia41得41a41i41a1=-==即方程22dxyd+y=xe ix的一个特解~y=-4ix2+41xe ix=-4ix2+41cosx+isinx=41x2sinx+41xcosx+i-41x2cosx+41xsinx取其虚部,得~2y=-41x2cos x+41x sin x 所以,所求方程的通解y =Y+~1y+~2y=C 1cosx +C 2sinx +101-513e3x-41x 2cosx +41xsinx综上所述,对于二阶常系数线性非齐次方程22dx y d +p dxdy +qy =fx 当自由项fx 为上述所列三种特殊形式时,其特解~y 可用待定系数法求得,其特解形式列表如下:自由项fx 形式特解形式fx =p n x当q ≠0时~y =Q n x当q =0,p ≠0时~y =Q n x当q =0,p =0时~y =x 2Q n xfx =p n xeαx当α不是特征方程根时~y =Q nxeαx当α是特征方程单根时~y =xQ n xe αx当α是特征方程重根时~y =x 2Q n xe αxfx =p n xe αx cos βx 或fx =p n xe αx sin βx利用欧拉公式e i βx =cos βx +isin βx,化为fx =p n xe α+i βx 的形式求特解,再分别取其实部或虚部以上求二阶常系数线性非齐次方程的特解的方法,当然可以用于一阶,也可以推广到高阶的情况.例8. 求y+3y ″+3y ′+y =e x 的通解解 对应的齐次方程的特征方程为r 3+3r 2+3r +1=0 r 1=r 2=r 3=-1所求齐次方程的通解Y =C 1+C 2x +C 3x 2e -x由于α=1不是特征方程的根因此方程的特解~y =a 0e x代入方程可解得a 0=81故所求方程的通解为y =Y +~y =C 1+C 2x +C 3x 2e -x+81e x.§ 欧拉方程下述n 阶线性微分方程a 0xnn n ax y d +a 1x n -11n 1n dxyd --+…+a n -1x dxdy+a ny =fx 称为欧拉方程,其中a 0,a 1,…a n 都是常数,fx 是已知函数.欧拉方程可通过变量替换化为常系数线性方程.下面以二阶为例说明.对于二阶欧拉方程a 0x 222dx y d +a 1x dxdy +a 2y =fx 作变量替换令x =e t,即t =ln x引入新变量t,于是有dx dy =dt dy dx dt =dt dy x 1=x 1dtdy22dx y d =dx d x 1dt dy =x 1dx d dt dy +dt dy dx d x 1 =x 122dt y d dx dt -2x 1dt dy =2x 122dt y d -2x 1dt dy 代入方程得a 022dt y d -dt dy +a 2dtdy+a 1y =fe t即 22dty d +002a a a dt dy +01a a y =0a 1fe t它是yt 的常系数线性微分方程.例9. 求x 222dx y d +x dx dy =6lnx -x1的通解. 解 所求方程是二阶欧拉方程作变换替换,令x =e t ,则dx dy =x 1dxdy22dx y d =2x 122dt y d -2x 1dt dy 代入原方程,可得 22dty d =6t -e -t两次积分,可求得其通解为 y =C 1+C 2t +t 3-e -t代回原来变量,得原方程的通解y =C 1+C 2lnx +lnx3-x1第八节 常系数线性方程组前面讨论的微分方程所含的未知函数及方程的个数都只有一个,但在实际问题中常遇到含有一个自变量的两个或多个未知函数的常微分方程组.本节只讨论常系数线性方程组,并且用代数的方法将其化为常系数线性方程的求解问题.下面以例说明.例1. 求方程组⎪⎩⎪⎨⎧=--=--)2(0y 3x 4dtdy)1(e y 2x dtdx t的通解.解 与解二元线性代数方程组中的消元法相类似,我们设法消去一个未知函数,由1得y =21 dtdx -x -e t3将其代入2得 21 22dt x d -dt dx -e t-4x -23 dtdx -x -e t=0 化简得22dt x d -4dtdx -5x =-2e t它是一个二阶常系数非齐次方程它的通解为 x =C 1e 5t+C 2e -t+41e t代入3得 y =2C 1e 5t-C 2e -t-21e t即所求方程组的通解为⎪⎪⎩⎪⎪⎨⎧--=++=--t t 2t 51t t2t 51e 21e C e C 2y e 41e C e C x例2. 求解方程组⎪⎩⎪⎨⎧++=+-=+)2(t 2y x dtdy dt dx )1(yt dt dydt dx 2的通解解 为消去y,先消去dtdy,为此将1-2得dtdx +x +2y +t =0即有 y =-21 dtdx+x +t 3代入2得dt dx -21dt d dt dx +x +t -x +21 dtdx +x +t -2t =0 即 22dt x d -2dtdx+x =3t -1 这是一个二阶常系数线性非齐次方程,解得x =C 1e t +C 2te t -3t -7代入3得 y =-C 1e t-C 221+te t+t +5 所以原方程组的通解为⎪⎩⎪⎨⎧+++--=--+=5t e )t 21(C e C y 7t 3te C e C x t2t 1t 2t 1。

9.4.3 二阶线性常系数微分方程的解法

9.4.3 二阶线性常系数微分方程的解法
x 2
5 t 2
B. 二阶线性常系数非齐次微分方程的通解 二阶线性常系数非齐次微分方程
y'' p y' q y f ( x )
求通解的步骤:
( 3)
( 1 ) 求对应齐次方程 y' ' p y' q y 0 的通解 yc ( x );
( 2) 求方程 y'' p y' q y f ( x ) 的一个特解 y p ( x );
A. 二阶线性常系数齐次微分方程的通解 二阶常系数齐次线性方程解法 --特征方程法
二阶线性常系数齐次微分方程
y'' p y' q y 0
分析:不妨设 y e
x
(1)
, 代入 (1) 式,
2 e x p e x qe x 0,
由于 e x 0, 则有
( D1 D2 x Dk x k 1 ) sin x ]
例 2 求下列方程通解 ( 5) (4) ( 3) (1) y 4 y 5 y 0;
解 特征方程 5 44 53 0,
特征根: 0 ( 3 重), 2 i.
所求通解为
y e 0 x (C 1 C 2 x C 3 x 2 )
x0
2 得 C 2 6.
x
( 3)
y''4 y'5 y 0.
解 特征方程 2 4 5 0, 特征根: 1,2 2 i . 2x 所以通解为 y e (C1 cos x C 2 sin x ).
利 用特 征值 求 通解 的 法 方, 可 推广 到 注. 一 阶或 高阶 常 系数 齐 线 次 性方 程的 求 解

二阶常系数线性微分方程的解法版共11页文档

二阶常系数线性微分方程的解法版共11页文档

第八章 8.4讲第四节 二阶常系数线性微分方程一、二阶常系数线形微分方程的概念形如 )(x f qy y p y =+'+''(1)的方程称为二阶常系数线性微分方程.其中p 、q 均为实数,)(x f 为已知的连续函数.如果0)(≡x f ,则方程式 (1)变成0=+'+''qy y p y(2)我们把方程(2)叫做二阶常系数齐次线性方程,把方程式(1)叫做二阶常系数非齐次线性方程. 本节我们将讨论其解法.二、二阶常系数齐次线性微分方程1.解的叠加性定理 1 如果函数1y 与2y 是式(2)的两个解, 则2211y C y C y +=也是式(2)的解,其中21,C C 是任意常数. 证明 因为1y 与2y 是方程(2)的解,所以有将2211y C y C y +=代入方程(2)的左边,得所以2211y C y C y +=是方程(2)的解.定理1说明齐次线性方程的解具有叠加性.叠加起来的解从形式看含有21,C C 两个任意常数,但它不一定是方程式(2)的通解.2.线性相关、线性无关的概念设,,,,21n y y y Λ为定义在区间I 内的n 个函数,若存在不全为零的常数,,,,21n k k k Λ使得当在该区间内有02211≡+++n n y k y k y k Λ, 则称这n 个函数在区间I 内线性相关,否则称线性无关.例如 x x 22sin ,cos ,1在实数范围内是线性相关的,因为 又如2,,1x x 在任何区间(a,b)内是线性无关的,因为在该区间内要使必须0321===k k k .对两个函数的情形,若=21y y 常数, 则1y ,2y 线性相关,若≠21y y 常数, 则1y ,2y 线性无关. 3.二阶常系数齐次微分方程的解法定理 2 如果1y 与2y 是方程式(2)的两个线性无关的特解,则212211,(C C y C y C y +=为任意常数)是方程式(2)的通解.例如, 0=+''y y 是二阶齐次线性方程,x y x y cos ,sin 21==是它的两个解,且≠=x y y tan 21常数,即1y ,2y 线性无关, 所以 x C x C y C y C y cos sin 212211+=+= ( 21,C C 是任意常数)是方程0=+''y y 的通解.由于指数函数rx e y =(r 为常数)和它的各阶导数都只差一个常数因子, 根据指数函数的这个特点,我们用rx e y =来试着看能否选取适当的常数r ,使rx e y =满足方程(2).将rx e y =求导,得把y y y ''',,代入方程(2),得因为0≠rx e , 所以只有 02=++q pr r(3)只要r 满足方程式(3),rx e y =就是方程式(2)的解.我们把方程式(3)叫做方程式(2)的特征方程,特征方程是一个代数方程,其中r r ,2的系数及常数项恰好依次是方程(2)y y y ,,'''的系数.特征方程(3)的两个根为 2422,1q p p r -±-=, 因此方程式(2)的通解有下列三种不同的情形.(1) 当042>-q p 时,21,r r 是两个不相等的实根.x r x r e y e y 2121,==是方程(2)的两个特解,并且≠=-x r r e y y )(2121常数,即1y 与2y 线性无关.根据定理2,得方程(2)的通解为x r x r e C e C y 2121+=(2) 当042=-q p 时, 21,r r 是两个相等的实根. 221p r r -==,这时只能得到方程(2)的一个特解x r e y 11=,还需求出另一个解2y ,且≠12y y 常数,设)(12x u y y =, 即 将222,,y y y '''代入方程(2), 得 整理,得由于01≠x r e , 所以 0)()2(1211=+++'++''u q pr r u p r u 因为1r 是特征方程(3)的二重根, 所以从而有 0=''u因为我们只需一个不为常数的解,不妨取x u =,可得到方程(2)的另一个解那么,方程(2)的通解为即 x r e x C C y 1)(21+=. (3) 当042<-q p 时,特征方程(3)有一对共轭复根 于是 x i x i e y e y )(2)(1,βαβα-+==利用欧拉公式 x i x e ix sin cos +=把21,y y 改写为21,y y 之间成共轭关系,取方程(2)的解具有叠加性,所以-1y ,-2y 还是方程(2)的解,并且≠==--x x e x e y y x x βββααtan cos sin 12常数,所以方程(2)的通解为 综上所述,求二阶常系数线性齐次方程通解的步骤如下:(1)写出方程(2)的特征方程(2)求特征方程的两个根21,r r(3)根据21,r r 的不同情形,按下表写出方程(2)的通解.例1求方程052=+'+''y y y 的通解.解: 所给方程的特征方程为所求通解为 )2sin 2cos (21x C x C e y x +=-.例 2 求方程0222=++S dt dS dtS d 满足初始条件2,400-='===t t S S 的特解.解 所给方程的特征方程为通解为 t e t C C S -+=)(21 将初始条件40==t S 代入,得 41=C ,于是t e t C S -+=)4(2,对其求导得 将初始条件20-='=t S 代入上式,得所求特解为例3求方程032=-'+''y y y 的通解.解 所给方程的特征方程为 0322=-+r r其根为 1,321=-=r r所以原方程的通解为 x x e C e C y 231+=-二、二阶常系数非齐次方程的解法1.解的结构定理3 设*y 是方程(1)的一个特解,Y 是式(1)所对应的齐次方程式(2)的通解,则*+=y Y y 是方程式(1)的通解.证明 把*+=y Y y 代入方程(1)的左端:*+=y Y y 使方程(1)的两端恒等,所以*+=y Y y 是方程(1)的解.定理4 设二阶非齐次线性方程(1)的右端)(x f 是几个函数之和,如)()(21x f x f qy y p y +=+'+'' (4)而*1y 与*2y 分别是方程 )(1x f qy y p y =+'+''与 )(2x f qy y p y =+'+''的特解,那么**+21y y 就是方程(4)的特解, 非齐次线性方程(1)的特解有时可用上述定理来帮助求出.2.)()(x P e x f m x λ=型的解法)()(x P e x f m x λ=,其中λ为常数,)(x P m 是关于x 的一个m 次多项式.方程(1)的右端)(x f 是多项式)(x P m 与指数函数x e λ乘积的导数仍为同一类型函数,因此方程(1)的特解可能为x e x Q y λ)(=*,其中)(x Q 是某个多项式函数.把 x e x Q y λ)(=*代入方程(1)并消去x e λ,得)()()()()2()(2x P x Q q p x Q p x Q m =+++'++''λλλ(5)以下分三种不同的情形,分别讨论函数)(x Q 的确定方法:(1) 若λ不是方程式(2)的特征方程02=++q pr r 的根,即02≠++q p λλ,要使式(5)的两端恒等,可令)(x Q 为另一个m 次多项式)(x Q m :代入(5)式,并比较两端关于x 同次幂的系数,就得到关于未知数m b b b ,,,10Λ的1+m 个方程.联立解方程组可以确定出),,1,0(m i b i Λ=.从而得到所求方程的特解为(2) 若λ是特征方程02=++q pr r 的单根, 即02,02≠+=++p q p λλλ,要使式(5)成立, 则)(x Q '必须要是m 次多项式函数,于是令用同样的方法来确定)(x Q m 的系数),,1,0(m i b i Λ=.(3) 若λ是特征方程02=++q pr r 的重根,即,02=++q p λλ 02=+p λ.要使(5)式成立,则)(x Q ''必须是一个m 次多项式,可令 用同样的方法来确定)(x Q m 的系数.综上所述,若方程式(1)中的x m e x P x f λ)()(=,则式(1)的特解为其中)(x Q m 是与)(x P m 同次多项式,k 按λ不是特征方程的根,是特征方程的单根或是特征方程的重根依次取0,1或2.例4 求方程x e y y 232-='+''的一个特解.解 )(x f 是x m e x p λ)(型, 且2,3)(-==λx P m对应齐次方程的特征方程为 022=+r r ,特征根根为2,021-==r r .λ=-2是特征方程的单根, 令x e xb y 20-=*,代入原方程解得故所求特解为 x xe y 223--=* .例5 求方程x e x y y )1(2-='-''的通解.解 先求对应齐次方程02=+'-''y y y 的通解. 特征方程为 0122=+-r r , 121==r r齐次方程的通解为 x e x C C Y )(21+=. 再求所给方程的特解由于1=λ是特征方程的二重根,所以把它代入所给方程,并约去x e 得比较系数,得于是 x e x x y )216(2-=* 所给方程的通解为 x e x x x C C y y y )6121(3221+-+=+=*3.x B x A x f ϖϖsin cos )(+=型的解法,sin cos )(x B x A x f ωω+=其中A 、B 、ω均为常数. 此时,方程式(1)成为x B x A q y p y ωωsin cos +=+'+'' (7)这种类型的三角函数的导数,仍属同一类型,因此方程式(7)的特解*y 也应属同一类型,可以证明式(7)的特解形式为 其中b a ,为待定常数.k 为一个整数.当ω±i 不是特征方程02=++q pr r 的根, k 取0; 当ω±i 不是特征方程02=++q pr r 的根, k 取1; 例6 求方程x y y y sin 432=-'+''的一个特解. 解 1=ω,ω±i i ±=不是特征方程为0322=-+r r 的根,0=k .因此原方程的特解形式为于是 x b x a y cos sin +-=*' 将*''*'*y y y ,,代入原方程,得解得 54,52-=-=b a原方程的特解为: x x y sin 54cos 52--=* 例7 求方程x e y y y x sin 32+=-'-''的通解.解 先求对应的齐次方程的通解Y .对应的齐次方程的特征方程为再求非齐次方程的一个特解*y .由于x e x x f -+=2cos 5)(,根据定理4,分别求出方程对应的右端项为,)(1x e x f =x x f sin )(2=的特解*1y 、*2y ,则**+=*21y y y 是原方程的一个特解.第 11 页 由于1=λ,ω±i i ±=均不是特征方程的根,故特解为 代入原方程,得比较系数,得解之得 51,101,41-==-=c b a . 于是所给方程的一个特解为所以所求方程的通解为希望以上资料对你有所帮助,附励志名言3条:1、理想的路总是为有信心的人预备着。

二阶常微分方程的解法

二阶常微分方程的解法

南京师范大学泰州学院毕业论文(设计)(一六届)题目:二阶常微分方程的解法院(系、部):数学科学与应用学院专业:数学与应用数学姓名:潘陆学号08120146指导教师:刘陆军南京师范大学泰州学院教务处制摘要:本文主要是介绍了二阶常微分方程众多解法中的三种,分别为特征方程法,拉普拉斯变换法和常数变易法,研究并讨论了二阶常微分方程在特征方程法中特征方程根为实根,复根和重根的情形。

我们选用了弹簧振子系统的振子运动,用这三种不同的方法来解决该问题。

关键词:二阶常微分方程;特征根法;常数变易法;拉普拉斯变换Abstract:The main purpose of this paper is the second-order ordinary many differential equation solution of three, respectively as the characteristic equation method, Laplace transform method and variation of constants method, study and discuss the second-order often differential equation in the characteristic equation of the roots of the characteristic equation for real roots, complex roots and root weight. We choose the spring oscillator the oscillator motion, these three different methods to solve the problem.Keywords: second order ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform目录1 绪论 (3)1.1 二阶常微分方程的起源和发展史 (3)1.2 二阶常微分方程的介绍 (3)1.3 研究二阶常微分方程的目的与意义 (4)2 二阶常系数常微分方程的几种解法 (5)2.1 特征方程法 (5)2.1.1 特征根是两个实根的情形 (5)2.1.2 特征根有重根的情形 (6)2.2 常数变易法 (7)2.3 拉普拉斯变换法 (9)3 二阶常微分方程解法的应用(分析例题) (11)3.1 特征方程法 (11)3.2 常数变易法 (13)3.3 拉普拉斯变换法 (14)4 结论和启示 (16)谢辞 (18)参考文献 (19)1 绪论1.1 二阶常微分方程的起源和发展史既然说到了微分方程,就不能不提到海王星的故事,它的发现是人类智慧的硕果,微分方程在其中扮演了重要的角色,并且在其中也包含数学演绎法的作用。

二阶微分方程解法(参考模板)

二阶微分方程解法(参考模板)

第六节 二阶常系数齐次线性微分方程教学目的:使学生掌握二阶常系数齐次线性微分方程的解法,了解二阶常系数非齐次线性微分方程的解法教学重点:二阶常系数齐次线性微分方程的解法 教学过程:一、二阶常系数齐次线性微分方程二阶常系数齐次线性微分方程: 方程 y+py +qy =0称为二阶常系数齐次线性微分方程, 其中p 、q 均为常数.如果y 1、y 2是二阶常系数齐次线性微分方程的两个线性无关解, 那么y =C 1y 1+C 2y 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可用公式2422,1q p p r -±+-= 求出.特征方程的根与通解的关系:(1)特征方程有两个不相等的实根r 1、r 2时, 函数x r e y 11=、x r e y 22=是方程的两个线性无关的解. 这是因为,函数x r e y 11=、x r e y 22=是方程的解, 又xr r xr x r e e e y y )(212121-==不是常数. 因此方程的通解为x r x r e C e C y 2121+=.(2)特征方程有两个相等的实根r 1=r 2时, 函数x r e y 11=、x r xe y 12=是二阶常系数齐次线性微分方程的两个线性无关的解.这是因为, x r e y 11=是方程的解, 又x r x r xr x r x r x r qxe e xr p e xr r xe q xe p xe 111111)1()2()()()(1211++++=+'+'' 0)()2(121111=++++=q pr r xe p r e x r x r ,所以xr xe y 12=也是方程的解, 且x e xe y y xr xr ==1112不是常数. 因此方程的通解为 x r x r xe C e C y 1121+=.(3)特征方程有一对共轭复根r 1, 2=a ib 时, 函数y =e(a +ib )x、y =e(a ib )x是微分方程的两个线性无关的复数形式的解. 函数y =e axcos bx 、y =e axsin bx 是微分方程的两个线性无关的实数形式的解. 函数y 1e(a +ib )x和y 2e(a ib )x都是方程的解 而由欧拉公式 得y 1e (a +ib )x e x (cos x i sin x )y 2e(aib )xe x (cos x i sin x )y 1y 22e x cos x )(21cos 21y y x e x +=βα y 1y 22ie x sin x )(21sin 21y y ix e x -=βα故e ax cos bx 、y 2=e axsin bx 也是方程解.可以验证, y 1=e ax cos bx 、y 2=e axsin bx 是方程的线性无关解. 因此方程的通解为y =e ax(C 1cos bx +C 2sin bx ). 求二阶常系数齐次线性微分方程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是两个不相等的实根, 因此所求通解为 y =C 1e -x+C 2e 3x.例2 求方程y+2y+y=0满足初始条件y|x=0=4、y|x=0=-2的特解.解所给方程的特征方程为r2+2r+1=0, 即(r1)20其根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=12i r2=12i是一对共轭复根因此所求通解为y=e x(C1cos2x+C2sin2x).n阶常系数齐次线性微分方程: 方程y(n) +p1y(n-1)+p2 y(n-2) + + p n-1y+p n y=0,称为n阶常系数齐次线性微分方程, 其中p1, p2 , , p n-1, p n都是常数.二阶常系数齐次线性微分方程所用的方法以及方程的通解形式, 可推广到n阶常系数齐次线性微分方程上去.引入微分算子D及微分算子的n次多项式L(D)=D n+p1D n-1+p2 D n-2 + + p n-1D+p n则n阶常系数齐次线性微分方程可记作(D n+p1D n-1+p2 D n-2 + + p n-1D+p n)y=0或L(D)y0注 D叫做微分算子D0y y D y y D2y y D3y y D n y y(n)分析令y e rx则L(D)y L(D)e rx(r n+p1r n-1+p2 r n-2 + + p n-1r+p n)e rx=L(r)e rx因此如果r是多项式L(r)的根则y e rx是微分方程L(D)y0的解n阶常系数齐次线性微分方程的特征方程L(r)r n+p1r n-1+p2 r n-2 + + p n-1r+p n0称为微分方程L(D)y0的特征方程特征方程的根与通解中项的对应: 单实根r 对应于一项: Ce rx;一对单复根r 1, 2=a ib 对应于两项: e ax(C 1cos bx +C 2sin bx );k 重实根r 对应于k 项: e rx (C 1+C 2x + +C k x k -1); 一对k 重复根r 1, 2=a ib 对应于2k 项:e ax[(C 1+C 2x + +C k x k -1)cos bx +( D 1+D 2x + +D k x k -1)sin bx ]. 例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=12i .因此所给微分方程的通解为y =C 1+C 2x +e x(C 3cos2x +C 4sin2x ). 例5 求方程y (4)+b 4y =0的通解, 其中b 0.解 这里的特征方程为 r 4+b 4=0. 它的根为)1(22,1i r ±=β, )1(24,3i r ±-=β.因此所给微分方程的通解为 )2sin2cos(212x C x C ey xβββ+=)2sin2cos(432x C x C exβββ++-.二、二阶常系数非齐次线性微分方程简介二阶常系数非齐次线性微分方程: 方程y +py +qy =f (x )称为二阶常系数非齐次线性微分方程, 其中p 、q 是常数. 二阶常系数非齐次线性微分方程的通解是对应的齐次方程 的通解y =Y (x )与非齐次方程本身的一个特解y =y *(x )之和:y =Y (x )+ y *(x ).当f (x )为两种特殊形式时, 方程的特解的求法: 一、 f (x )=P m (x )e lx型当f (x )=P m (x )e lx时, 可以猜想, 方程的特解也应具有这种形式. 因此, 设特解形式为y *=Q (x )e lx , 将其代入方程, 得等式 Q(x )+(2l +p )Q(x )+(l 2+pl +q )Q (x )=P m (x ).(1)如果l 不是特征方程r 2+pr +q =0 的根, 则l 2+pl +q 0. 要使上式成立, Q (x )应设为m 次多项式:Q m(x)=b0x m+b1x m-1+ +b m-1x+b m,通过比较等式两边同次项系数, 可确定b0, b1, , b m, 并得所求特解y*=Q m(x)e lx.(2)如果l是特征方程r2+pr+q=0 的单根, 则l2+pl+q=0, 但2l+p0, 要使等式Q(x)+(2l+p)Q(x)+(l2+pl+q)Q(x)=P m(x).成立, Q(x)应设为m+1 次多项式:Q(x)=xQ m(x),Q m(x)=b0x m+b1x m-1+ +b m-1x+b m,通过比较等式两边同次项系数, 可确定b0, b1, , b m, 并得所求特解y*=xQ m(x)e lx.(3)如果l是特征方程r2+pr+q=0的二重根, 则l2+pl+q=0, 2l+p=0, 要使等式Q(x)+(2l+p)Q(x)+(l2+pl+q)Q(x)=P m(x).成立, Q(x)应设为m+2次多项式:Q(x)=x2Q m(x),Q m(x)=b0x m+b1x m-1+ +b m-1x+b m,通过比较等式两边同次项系数, 可确定b0, b1, , b m, 并得所求特解y*=x2Q m(x)e lx.综上所述, 我们有如下结论: 如果f(x)=P m(x)e lx, 则二阶常系数非齐次线性微分方程y+py+qy =f(x)有形如y*=x k Q m(x)e lx的特解, 其中Q m(x)是与P m(x)同次的多项式, 而k按l不是特征方程的根、是特征方程的单根或是特征方程的的重根依次取为0、1或2.例1 求微分方程y-2y-3y=3x+1的一个特解.解这是二阶常系数非齐次线性微分方程, 且函数f(x)是P m(x)e lx型(其中P m(x)=3x+1, l=0).与所给方程对应的齐次方程为y-2y-3y=0,它的特征方程为r2-2r-3=0.由于这里l=0不是特征方程的根, 所以应设特解为y*=b0x+b1.把它代入所给方程, 得-3b0x-2b0-3b1=3x+1,比较两端x同次幂的系数, 得⎩⎨⎧=--=-13233100b b b -3b 0=3, -2b 0-3b 1=1.由此求得b 0=-1, 311=b . 于是求得所给方程的一个特解为 31*+-=x y . 例2 求微分方程y-5y +6y =xe 2x的通解.解 所给方程是二阶常系数非齐次线性微分方程, 且f (x )是P m (x )e lx型(其中P m (x )=x , l =2). 与所给方程对应的齐次方程为y -5y +6y =0,它的特征方程为r 2-5r +6=0.特征方程有两个实根r 1=2, r 2=3. 于是所给方程对应的齐次方程的通解为Y =C 1e 2x +C 2e 3x .由于l =2是特征方程的单根, 所以应设方程的特解为y *=x (b 0x +b 1)e 2x .把它代入所给方程, 得 -2b 0x +2b 0-b 1=x . 比较两端x 同次幂的系数, 得 ⎩⎨⎧=-=-0212100b b b -2b 0=1, 2b 0-b 1=0.由此求得210-=b , b 1=-1. 于是求得所给方程的一个特解为 x e x x y 2)121(*--=. 从而所给方程的通解为x x x e x x e C e C y 223221)2(21+-+=. 提示y *=x (b 0x +b 1)e 2x (b 0x 2+b 1x )e 2x[(b 0x 2+b 1x )e 2x][(2b 0x +b 1)(b 0x 2+b 1x )×2]e2x[(b 0x 2+b 1x )e 2x][2b 02(2b 0x b 1)×2(b 0x 2+b 1x )×22]e 2xy *5y *6y *[(b 0x 2+b 1x )e 2x]5[(b 0x 2+b 1x )e 2x]6[(b 0x 2+b 1x )e 2x][2b 02(2b 0x b 1)×2(b 0x 2+b 1x )×22]e 2x5[(2b 0x +b 1)(b 0x 2+b 1x )×2]e2x6(b 0x 2+b 1x )e 2x[2b 04(2b 0x b 1)5(2b 0x +b 1)]e 2x[2b 0x +2b 0b 1]e 2x方程y+py +qy =e lx[P l (x )cos wx +P n (x )sin wx ]的特解形式应用欧拉公式可得e lx [P l (x )cos wx +P n (x )sin wx ]]2)(2)([ ie e x P e e x P e x i x i nx i xi l x ωωωωλ---++=x i nl x i n l e x iP x P e x iP x P )()()]()([21)]()([21ωλωλ-+++-=x i x i e x P e x P )()()()(ωλωλ-++=,其中)(21)(i P P x P n l -=, )(21)(i P P x P n l +=. 而m =max{l , n }. 设方程y+py+qy =P (x )e(l +iw )x的特解为y 1*=x k Q m (x )e(l +iw )x,则)(1)(*ωλi m k e x Q x y -=必是方程)()(ωλi e x P qy y p y -=+'+''的特解, 其中k 按l iw 不是特征方程的根或是特征方程的根依次取0或1. 于是方程y+py +qy =e lx[P l (x )cos wx +P n (x )sin wx ]的特解为x i m k x i m k e x Q x e x Q x y )()()()(*ωλωλ-++=)sin )(cos ()sin )(cos ([x i x x Q x i x x Q e x m m x k ωωωωλ-++= =x k e lx[R(1)m(x )cos wx +R(2)m(x )sin wx ].综上所述, 我们有如下结论:如果f (x )=e lx[P l (x )cos wx +P n (x )sin wx ], 则二阶常系数非齐次线性微分方程y+py +qy =f (x )的特解可设为y *=x k e lx [R (1)m (x )cos wx +R (2)m (x )sin wx ],其中R(1)m(x )、R(2)m(x )是m 次多项式, m =max{l , n }, 而k 按l +i w (或l -iw )不是特征方程的根或是特征方程的单根依次取0或1. 例3 求微分方程y+y =x cos2x 的一个特解.解 所给方程是二阶常系数非齐次线性微分方程,且f (x )属于e lx[P l (x )cos wx +P n (x )sin wx ]型(其中l =0, w =2, P l (x )=x , P n (x )=0). 与所给方程对应的齐次方程为y +y =0,它的特征方程为r 2+1=0.由于这里l +iw =2i 不是特征方程的根, 所以应设特解为y *=(ax +b )cos2x +(cx +d )sin2x .把它代入所给方程, 得(-3ax -3b +4c )cos2x -(3cx +3d +4a )sin2x =x cos2x . 比较两端同类项的系数, 得 31-=a , b =0, c =0, 94=d . 于是求得一个特解为 x x x y 2sin 942cos 31*+-=. 提示y *=(ax +b )cos2x +(cx +d )sin2x .y *=a cos2x 2(ax +b )sin2x +c sin2x +2(cx +d )cos2x(2cx +a2d )cos2x +(2ax 2b c )sin2xy *=2c cos2x 2(2cx +a 2d )sin2x 2a sin2x +2(2ax 2b c )cos2x(4ax4b4c )cos2x(4cx 4a 4d )sin2xy *y *(3ax 3b 4c )cos2x (3cx 4a 3d )sin2x由⎪⎩⎪⎨⎧=--=-=+-=-0340304313d a c c b a 得31-=a , b =0, c =0, 94=d .(注:文档可能无法思考全面,请浏览后下载,供参考。

二阶常系数线性微分方程的解法

二阶常系数线性微分方程的解法
2.2 二阶常系数线性微分方程 的解法
数学系 贺 丹
2.2 二阶常系数线性微分方程的解法
若二阶线性微分方程为 ay by cy f ( x) ,其中 a, b, c 均为常数,则称该方程为二阶常系数线性微分方程。
二阶常系数线性齐次方程的解法
ay by cy 0 ,

猜想方程①具有 y erx 形式的解,其中 r 为待定常数,
方程②叫做方程①的特征方程。
按特征方程的两个根 r1, r2 的三种可能情况: 1. r1 r2 是两个不相等的实根; 2. r1 r2 是两个相等的实根;
3. r1 i , r2 i 是一对共轭复数。
3
2.2 二阶常系数线性微分方程的解法
1 . 特 征 方 程 的 根 是 两 个 不 相 等 实 数 的 情 形 。
∵函数 y1 和 y2 都是方程①的解,且它们是线性无关的,
∴ 方程①的通解为 y C1 y1 C2 y2 ,即
y e x (C 1 co x s C 2 six n )
(其中 , 为特征方程的复根的实部及虚部)。
7
2.2 二阶常系数线性微分方程的解法
小结:求二阶常系数线性齐次方程通解的步骤
高阶常系数线性齐次方程的解法 n 阶常系数线性齐次方程为
a0 y(n) a1 y(n1) an1 y an y 0 , ③
其特征方程为 a0r n a1r n1 an1r an 0 . ④
方程②是一个一元 n 次方程, 有 n 个根。类似二阶常系
数线性齐次方程,相应地可得到方程①的 n 个线性无关
∵ erx 0 , 2ar b 0, ar 2 br c 0 , ∴ u( x) 0 ,
取 u( x) 0 的一个解 u( x) x ,则 y2 xerx 。

第三节-二阶常系数线性微分方程的解法

第三节-二阶常系数线性微分方程的解法

(*)
情形3
若 r 是特征方程的二重根, 即 r 2 ar b 0 ,
且 2r a 0 , 则令 Q( x) x 2 Qm ( x) , 即
y x Qm ( x ) e
2

rx
15
2 Q ( 2r a )Q ( r ar b)Q Pm ( x )
于是 (2)的通解为
y (C 1 C 2 x ) e
1 x
.
6
情形3 若 0 , 则特征方程(3)有一对共轭复根
1,2 i
x x y e cos x , y e sinx 可以证明, 1 2
是(2)的解,且线性无关, 所以方程(2)的通解为
y e (C1 cos x C 2 sin x )
1 x
,得
2 u (21 a)u (1 a1 b)u 0 ,
因为 1 是方程 2 a b 0 的二重根,
2 故有 1 a1 b 0 , 21 a 0 ,
1 x u x , 即得 y2 x e , u 0 , 取特解

rx
情形2 若 r 是特征方程的单根, 即 r 2 ar b 0 ,
而 2r a 0 , 则令 Q( x ) xQm ( x ) , 即
y xQm ( x ) e r x
14
2 Q ( 2r a )Q ( r ar b)Q Pm ( x )
2
3x
的通解.
特征方程 6 9 0 , 特征根 1, 2 3 ,
对应齐次方程通解 Y (C1 C2 x) e3 x .
因为 r 3 是二重特征根,

二阶线性常微分方程的解法

二阶线性常微分方程的解法

二阶线性常微分方程的解法在数学中,二阶线性常微分方程是一个常见且重要的概念。

本文将介绍二阶线性常微分方程的解法,帮助读者更好地理解和应用这一知识点。

一、二阶线性常微分方程的定义二阶线性常微分方程是指形如下式的微分方程:y''(x) + p(x)y'(x) + q(x)y(x) = g(x)其中y(x)是未知函数,p(x),q(x)和g(x)是已知函数,一般假设其在所考虑的区间上连续。

二、齐次方程的解法首先,我们来研究二阶线性常微分方程的齐次形式,即g(x)为零的情况。

这类方程的解法非常有规律性。

假设y1(x)和y2(x)是二阶线性常微分方程的两个解,那么线性组合c1y1(x) + c2y2(x)也是该方程的解,其中c1和c2是任意常数。

因此,我们可以找到两个解y1(x)和y2(x),并通过线性组合的方式得到方程的通解。

具体的解法有三种情况。

1. 两个不同实数根当方程的特征方程有两个不同的实数根r1和r2时,对应的两个解分别为y1(x) = e^(r1x)和y2(x) = e^(r2x)。

2. 重根当方程的特征方程有一个重根r时,对应的两个解分别为y1(x) =e^(rx)和y2(x) = xe^(rx)。

3. 复数根当方程的特征方程有共轭复数根a±bi时,对应的两个解分别为y1(x) = e^(ax)cos(bx)和y2(x) = e^(ax)sin(bx)。

三、非齐次方程的解法对于非齐次方程,我们需要借助齐次方程的解,通过特解的方法来求解。

假设y1(x)和y2(x)是齐次方程的两个解,我们可以得到非齐次方程的特解为y(x) = u1(x)y1(x) + u2(x)y2(x),其中u1(x)和u2(x)是待定函数。

具体的求解步骤是:1. 将待求特解y(x)代入原方程,消去齐次方程的项,得到u1'(x)y1(x) + u2'(x)y2(x) = g(x)。

二阶常系数线性微分方程的解法-19页PPT文档资料

二阶常系数线性微分方程的解法-19页PPT文档资料
是特征方程的重根
u0
取 u = x , 则得 y2xer1x,因此原方程的通解为 y(C 1C 2x)er1x
机动 目录 上页 下页 返回 结束
3. 当 p24q0时, 特征方程有一对共轭复根
这时原方程有两个复数解:
y1e(i)x e x(cx o isix n ) y2e(i)x e x(cx o isix n )
比较系数, 得 因此特解为 y * x ( 5 c3 o x 3 s s3 i x )n 所求通解为
x (5 c3 o x 3 s3 ix ) n
机动 目录 上页 下页 返回 结束
定理 4.
分别是方程
y p y q f y k ( x )( k 1 ,2 , ,n )
利用解的叠加原理 , 得原方程的线性无关特解:
y11 2(y1y2) excosx y221i(y1y2)exsinx
因此原方程的通解为
y e x ( C 1 co x C s 2 six ) n
机动 目录 上页 下页 返回 结束
小结: ypyqy0(p ,q 为常 ) 数
因此原方程的通解为
例2. 求解初值问题
d2s dt2
2ds dt
s
0
st04,
ds dt
t 0 2
解: 特征方程 r22r10有重根 r1r21,
因此原方程的通解为 s (C 1 C 2t)e t
利用初始条件得
C14, C2 2
于是所求初值问题的解为
机动 目录 上页 下页 返回 结束
(YpYqY)
f(x ) 0 f(x )
复习 目录 上页 下页 返回 结束
故 y Y (x ) y * (x )是非齐次方程的解, 又Y 中含有

二阶常系数常微分方程的初等解法求解技巧

二阶常系数常微分方程的初等解法求解技巧

二阶常系数常微分方程的初等解法求解技巧常微分方程是数学中的一个重要的分支,在数学、物理、工程等领域都有广泛的应用。

其中,二阶常系数常微分方程是最基本的一类常微分方程,其形式如下:$$ a\frac{{d^2y}}{{dt^2}}+b\frac{{dy}}{{dt}}+cy = 0 $$其中,a、b、c是常数,y是未知函数。

一、特征根为实数的情况1.首先,我们将二阶常系数常微分方程变形成特征方程:$$ a\lambda^2+b\lambda+c=0 $$2.求解特征方程得到两个实根,假设为λ1和λ23.根据两个实根求得特解的形式,形式如下:$$ y = C_1e^{\lambda_1 t} + C_2e^{\lambda_2 t} $$其中,C1和C2是待定常数。

二、特征根为复数的情况1.将二阶常系数常微分方程变形成特征方程。

2.求解特征方程得到两个复根,假设为α±βi。

3.根据两个复根求得特解的形式,形式如下:$$ y = e^{\alpha t}(C_1cos(\beta t) + C_2sin(\beta t)) $$其中,C1和C2是待定常数。

三、待定系数法待定系数法是一种适用于二阶常系数常微分方程有特定形式解的求解方法。

1. 如果方程右侧是其中一个函数的线性组合,我们可以假设原方程的特解为该函数的线性组合形式。

例如,如果方程右侧是常数1和指数函数e^kt的线性组合:$$ y_p(t) = A + Be^{kt} $$其中,A、B是待定常数,k是常数。

2.将上述假设代入原方程,得到一个关于A、B和k的代数方程。

3.解代数方程,求得A、B和k的值。

4. 特解为$$ y_p(t) = A + Be^{kt} $$其中,A、B是待定常数,k 是常数。

总结:以上是二阶常系数常微分方程的初等解法求解技巧。

通过找到二阶常系数常微分方程的特征根或使用待定系数法,我们可以求得其通解。

这些技巧在解决实际问题中非常有用,例如在振动、电路等领域的应用中常常会遇到二阶常系数常微分方程的求解。

2.2二阶常系数线性微分方程的解法

2.2二阶常系数线性微分方程的解法
当自由项 f (x) 为两种特殊类型函数时方程②特解的求 法—待定系数法。
13
2.2 二阶常系数线性微分方程的解法
1. f (x) Pm(x)ex ( 其中 pm ( x)是 x 的 m 次多项式 )
这时方程②为 ay by cy Pm ( x)ex

可以设 y Q( x)ex ( 其中Q( x) 是多项式 ) 。
例 1.求方程 y 5 y 6 y 2x 3 的特解。
解: f ( x) 2x 3 (2x 3)e0x ,
属 f ( x) Pm ( x)e x 型( m 1, 0 ),
特征方程为 r2 5r 6 0 , r1 2 , r2 3 ,
∵ 0 不是特征根,
∴设特解为 y Q1( x)e0x Aox A1 ,
得 erx (ar 2 br c) 0 ,但 erx 0 ,故有
ar 2 br c 0 ,

2
2.2 二阶常系数线性微分方程的解法
ar 2 br c 0 ,

若 r 是一元二次方程②的一个根,则 y erx 就是 方程①的一个特解。
方程②叫做方程①的特征方程。
按特征方程的两个根 r1, r2 的三种可能情况: 1. r1 r2 是两个不相等的实根; 2. r1 r2 是两个相等的实根;
9
2.2 二阶常系数线性微分方程的解法
高阶常系数线性齐次方程的解法 n 阶常系数线性齐次方程为
a0 y(n) a1 y(n1) an1 y an y 0 , ③
其特征方程为 a0r n a1r n1 an1r an 0 . ④
方程②是一个一元 n 次方程, 有 n 个根。类似二阶常系
Q( x) 应为 m 次多项式 , Q( x) 应为 m 1 次多项式 ,
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

二阶微分方程:
时为非齐次
时为齐次,0)(0)()()()(22≠≡=++x f x f x f y x Q dx dy x P dx y d 二阶常系数齐次线性微分方程及其解法:
2
122,)(2,,(*)0)(1,0(*)r r y y y r r q pr r q p qy y p y 式的两个根、求出的系数;
式中的系数及常数项恰好是,,其中、写出特征方程:求解步骤:
为常数;
,其中∆'''=++∆=+'+''式的通解:出的不同情况,按下表写、根据(*),321r r
二阶常系数非齐次线性微分方程

为常数;
型,为常数
,]sin )(cos )([)()()(,)(x x P x x P e x f x P e x f q p x f qy y p y n l x m x ωωλλλ+===+'+''
二阶常系数非齐次线性微分方程的一般形式是 ''+'+=y py qy f x () (1)
其中p q ,是常数。

方程(1)的通解为对应的齐次方程
0=+'+''qy y p y (2)
的通解Y 和方程(1)的一个特解*y 之和。

即 *y Y y +=.我们已解决了求二阶常系数齐
次线性方程通解的问题,所以,我们只需讨论求二阶常系数非齐次线性微分方程的特解*
y 的方法。

下面我们只介绍当方程(1)中的)(x f 为如下两种常见形式时求其特解*y 的方法。

一、
f x e P x x m ()()=⋅λ型 由于方程(1)右端函数f x ()是指数函数e x λ⋅与m 次多项式P x m ()的乘积,而指数
函数与多项式的乘积的导数仍是这类函数,因此,我们推测:
方程(1)的特解应为
y e Q x x *⋅=λ()( Q x ()是某个次数待定的多项式 ) y e Q x e Q x x x *⋅⋅'=+'λλλ()() y e Q x Q x Q x x *⋅"=⋅+'+''λλλ[()()()]22 代入方程(1),得
e Q x p Q x p q Q x e P x x x m λλλλλ⋅⋅⋅''++'+++≡⋅[()()()()()]()22 消去e x λ⋅,得
''++'+++≡Q x p Q x p q Q x P x m ()()()()()()22λλλ (3) 讨论
01、如果λ不是特征方程
r pr q 20++=的根。


02≠++q p λλ 由于P x m
()是一个m 次的多项式,欲使(3)的两端恒等,那未Q x ()必为一个m 次多项式,设为
Q x b x b x b x b m m m m m ()=++++--0111Λ
将之代入(3),比较恒等式两端x 的同次幂的系数,就得到以b b b b
m m 01
1,,,,Λ-为未知数的m +1个线性方程的联立方程组,解此方程组可得到这m +1个待定的系数,并得到特解 y e Q x x m *⋅=λ()
02、如果λ是特征方程
r pr q 20++=的单根。


λλ20++=p q ,但 20λ+≠p 欲使(3)式的两端恒等,那么'Q x ()必是一个m 次多项式。

因此,可令 Q x x Q x m ()()=⋅
并且用同样的方法来确定)(x Q 的系数b b b b m m 0
11,,,,Λ-。

03、如果λ是特征方程
r pr q 20++=的二重根。


λλ20++=p q ,且 20λ+=p 。

欲使(3)式的两端恒等,那么''Q x ()必是一个m 次多项式
因此, 可令
Q x x Q x m ()()=⋅2 并且用同样的方法来确定)(x Q 的系数b b b b m m 011,,,,Λ-。

综上所述,我们有结论
如果
f x e P x x m ()()=⋅λ,则方程(1)的特解形式为 y x Q x e k m x *⋅=()λ
其中Q x m ()是与P x m ()同次的多项式,k 的取值应满足条件 k =⎧⎨⎪⎩⎪012λλλ不是特征方程的根是特征方程的单根
是特征方程的二重根
例1求 ''-'+=y y y xe x 562 的通解。

解 特征方程为 0652=+-r r
特征根为 3,221==r r
齐次方程的通解为 x x e C e C Y 3221+=
因为2=λ是特征单根,所以,设非齐次方程的特解为
y x b x b e x *=+()012
则 *'y =+++[()]222020112b x b b x b e x
*''y =++++[()]484240201012b x b b x b b e x
将上述三式代入原方程,得
()-+-≡2200122b x b b e xe x x ,
比较恒等式两端的系数,得
-=-=⎧⎨⎩2120001b b b
解得
21
0-
=b , 11-=b 因此 x
e x x y 2)121(*+-= 所以方程的通解为
y c e c e x x e x x x
=+-+12232121()
二、
f x e P x x P x x x l n ()[()cos ()sin ]=+λωω型 由于方程(1)右端函数为
[]x x p x x p e n l x ωωλsin )(cos )(+,这种形式得到非齐次方程的特解*y 的过程稍微复杂些,所以我们这里就只给出结论
y x e R x x R x x k x m m *=+λωω[()cos ()sin ]()()12
其中,R x m ()()1、R x m ()()2是两个m 次多项式,m l n =max{,},

⎩⎨⎧++=是特征方程的根若不是特征方程的根若ωλωλi i k 10 例2求方程 ''+
=y y x x cos2 的通解。

解 特征方程 r
210+= 特征根 i
r ±=2,1 齐次方程的通解为 x C x C Y sin cos 21+=
这里1,2,0===m ωλ,由于i i 2=+ωλ不是特征方程的根,所以设方程的特解为 y ax b x cx d x *=+++()cos ()sin 22
代入原方程,得
x x x a d Cx x C b ax 2cos 2sin )433(2cos )433(=++-+--
比较两端同类项的系数,得
⎪⎪⎩⎪⎪⎨⎧=--=-=+-=-0430
304313a d C C b a
解得
94,0,0,31===-=d C b a 于是 y x x x *
=-+132492cos sin 所以非齐次方程的通解为
y c x c x x x x =+-+12132492cos sin cos sin。

相关文档
最新文档