§4.2常系數(shù)線性微分方程的解法常微分方程課件高教社王高雄教材配套ppt.pps

4.2 常系數(shù)線性微分方程的解法,復值函數(shù)和復值解 常系數(shù)齊次線性特征方程 常系數(shù)非齊次線性方程比較系數(shù)法 拉普拉斯變換法 質點振動,復值函數(shù)和復值解,設實變量t在區(qū)間atb上有實函數(shù)(t)、(t)， 為虛數(shù)，則可在區(qū)間上定義復值函數(shù) z (t)= (t)+i(t) ，稱(t)為其實部, (t)為其虛部。 復值函數(shù)和實函數(shù)有同樣的極限、連續(xù)、導數(shù)、微分、積分及數(shù)乘、函數(shù)和等性質。 滿足微分方程(1) 的復值函數(shù)z(t)稱為方程(1)的復值解。,復值函數(shù)性質,極限 連續(xù) 導數(shù) 微分 積分及函數(shù)和,復指數(shù)函數(shù)性質,設、為實數(shù)，t為實變量，則K= +i 為復數(shù)，復指數(shù)函數(shù)定義為 有 稱復數(shù)K= +i的共軛復數(shù)為 有 復指數(shù)函數(shù)有性質(K、K1 、 K2 為復數(shù)),復指數(shù)函數(shù)性質,證,實系數(shù)線性微分方程的復值解,定理8 對實系數(shù)函數(shù)微分方程(*)的復值解的實部、虛部及其共軛復值函數(shù)均是原方程(*)的解。 定理9 設非齊函數(shù)為復值函數(shù)：f(t)=u(t)+iv(t) 如復值函數(shù)x=U(t)+iV(t)是微分方程(*)的復值解，則實函數(shù)U(t)、V(t)分別是實微分方程的解：,常系數(shù)齊次線性特征方程,系數(shù)為實常數(shù)時的常系數(shù)齊次線性方程 可以尋求指數(shù)函數(shù)形式的解 代入方程右端可得到 其中 是的n次多項式。 稱為微分方程的特征方程。特征方程的根稱為特征根。,(a)實單特征根,(a) 為(實)單根時，齊次線性方程有解 et 設1, 2, n是特征方程的n個相異或相同但非重的實根，則相應的微分方程有n個解 因 i j (i j),行列式 即n個解在區(qū)間上線性無關，構成的基本解組。 方程有通解,(b)復特征根,(b)代數(shù)方程有復根 ,因方程的系數(shù)為實數(shù)， 復根將成對共軛地出現(xiàn)，即 也是特征根 有兩共軛復值解 對應的兩實值解：,(c)(1)重 0 實根,(c) = 1為k重實根時，它表示為 (1) 當1=0時，特征方程有因子k ，于是 即特征方程的形狀為 對應的微分方程為 易見它有k個解 1, t, t2, , tk. 且它們是線性無關的。 因此特征方程的k重零根對應微分方程的k個線性無關的解: 1, t, t2, , tk-1.,(c)(2)重非0實根,(2) = 10為k重實根時，作變量變換 因 有 于是微分方程化為 其中b1,bn仍為常數(shù)。相應的特征方程為 直接計算易得 即 且 因此F()=0的根對應于G()=0的根，且重數(shù)相同。 問題化為前面已討論的情形。,(續(xù)) (c)(2)重非0實根,已知特征方程 G()=0 的k1重根1=0對應于微分方程L1=0的k1個解y=1, t, t2, , tk-1。因此，對應于特征方程F()=0的k1重根1，微分方程L=0有個解 同樣，假設特征方程L=0的其他根2, 3, m的重次依次為k2, k3, km ，且 則微分方程L=0對應地有解,(續(xù)) (c)(2)重非0實根,證 全體個解構成微分方程L=0的基本解組。 反證法。假設這些函數(shù)線性相關，則有 設多項式Pm(t)不恒為零。將上式除于 e1t ，并對t微分k1次，可得 其中 Qr(t)與Pr(t)次數(shù)相同且Qr(t)不恒為零。等式類似，但項數(shù)減少了。,(續(xù)) (c)(2)重非0實根,進一步對施于同樣手續(xù)，除于e(2-1) t 并微分k2次， 則得到項數(shù)更少的類似多項式。繼續(xù)經m-1次后將得到等式 這不可能因 其中Wm是次數(shù)低于Pm(t)的次數(shù)的多項式。 而Rm(t)與Pm(t)次數(shù)相同且Rm(t)不恒為零。 假設函數(shù)線性相關引出矛盾 證明了全部n個解線性無關且構成微分方程的基本解組,(d) k重共軛復根時,= +i 微分方程有k對共軛復值解 或對應的2k對實值解： 類似的，對含有單實根、共軛復根、重實根及重共軛復根的混合情形， 可由前面討論過的分別進行處理，求得對應其特征根的微分方程的解，且這些解線性無關。 當所有根計及其重次總和為微分方程的次數(shù)時，這些解構成微分方程的基本解組。,例1 求方程 的通解,解 特征方程 4-1=0的根為 1=1, 2=-1, 3=i, 4=-i. 有兩個單實根和一對單共軛復根， 故方程的通解為 這里c1, c2, c3, c4,為任意常數(shù)。,例2 求解方程,解 特征方程3+1=0有根 通解為 其中c1, c2, c3為任意常數(shù)。,例3 求方程 的通解,解 特征方程 的根為三重根=1 。 因此方程有通解 其中c1, c2, c3為任意常數(shù)。,例4 求解方程,解 特征方程 的根為=i重根。 方程有4實值解 得通解 其中c1, c2, c3 , c4為任意常數(shù)。,歐拉方程,稱為歐拉方程。其中a1,a2,an 為實常數(shù)。 引進自變量變換 x=et, t=ln x (x=-et同樣，取 t=ln |x| ) 有 代入歐拉方程得常系數(shù)齊次線性微分方程,例5 求解方程,解 令x=et ,方程變?yōu)?特征方程為 有重根=1,通解為 原方程的通解為,常系數(shù)非齊次線性方程比較系數(shù)法,常系數(shù)非齊次線性方程 可按非齊項f(t)的類型用比較系數(shù)法求方程特解。 類型： 特解為 其中k為特征方程F()=0的根的重數(shù)。 當不是特征根時k=0;單根時k=1。 可以將特解代入方程比較t的同冪次項系數(shù)確定Bi。,例6 求方程 的通解,解 先求對應的齊次線性微分方程 的通解。其特征方程為2-2-3=0有根1=3,2=-1. 齊次方程有通解 再求非齊次方程的一個特解，這里f(t)=3t+1, = 0。 因 = 0不是特征方程的根，k=0。 特解形式為 x=A+Bt，其中A,B為待定常數(shù)。 將x代入原方程， 比較t的同冪次系數(shù)得 解為 即 原方程的通解為,例7 求方程 的通解,解 特征方程為3+32+3+1=0，有三重根1,2,3=1 齊次方程有通解 再求非齊次方程 的一個特解， 特解形式為 將特解代入原方程 比較系數(shù)得 從而 故原方程的通解為,b) 類型,其中,為常數(shù)，A(t),B(t)是的最髙次數(shù)為m次的實系數(shù)多項式。 設= +i為特征方程F()=0的k重根。 則方程有特解 這里P(t),Q(t)為待定的的最髙次數(shù)為m次的實系數(shù)多項式。 將特解代入方程，比較t的同冪次系數(shù)， 可得一系列線性方程，解之得特解。 注意：待定系數(shù)法的關鍵是正確寫出特解形式! P(t),Q(t)應為m次完全多項式。,例8 求方程 的通解,解 特征方程為2+4+4=( +2)2=0，有重根1,2=-2 齊次方程有通解 求非齊次方程 的一個特解， 因=2i不是特征方程的根，特解形式為 將特解代入原方程 比較系數(shù)得 即 故原方程的通解為,拉普拉斯變換法,實或復值連續(xù)函數(shù) f(t) 當 t0 滿足 | f(t)|上有定義。 對常系數(shù)非齊線性方程 及初始條件 設f(t)為連續(xù)函數(shù)且滿足原函數(shù)條件。 則方程的解 x(t) 及其各階導數(shù)均是原函數(shù)。,(續(xù)) 拉普拉斯變換法,記 有 如對方程兩邊進行拉普拉斯變換，利用線性性質可得關系式,(續(xù)) 拉普拉斯變換法,即 或 A(s)X(s)=F(s)+B(s) 其中A,B均為已知多項式。因此有 這是方程滿足初始條件的解的象函數(shù)X(s)。 再利用拉普拉斯反變換公式可得方程的解x(t)。 此即為拉普拉斯變換方法。,拉普拉斯變換和反變換表,例9 求方程 滿足初值條件 x(0)=0 的解,解 對方程兩端施行拉普拉斯變換， 得到象函數(shù)方程 由初值條件x(0)=0得 查拉普拉斯變換表知 和 的原函數(shù)分別為 e2t 和 et 利用線性性質，求得X(s)的原函數(shù)為 這就是所要求的解。,例10 求解方程,解 先令 = t-1將方程化為 再對方程兩端實行拉普拉斯變換，利用初值條件得到 即 查拉普拉斯變換表得 從而 此即為所要求的解。,例11 求解方程,解 對方程兩端施行拉普拉斯變換得 即 把上式右端分解為部分分式 各分式分別查拉普拉斯變換表， 利用線性性質，求得的原函數(shù)即方程的解,1-1 例2 數(shù)學擺 阻力與近似方程,近似方程 令,如擺在一個粘性介質中運動，阻力系數(shù)為，則方程為:,則數(shù)學擺近似方程,帶粘性介質的數(shù)學擺近似方程,數(shù)學擺方程,質點線性振動,質點線性振動 (1) 無阻尼 (2) 有阻尼 (a)小阻尼 (2) 大阻尼 (3) 臨界阻尼 (3) 無阻尼強迫振動(= 0) (a) p (b) p= (4) 小阻尼強迫振動 質點非線性振動 第6章討論,質點振動,振動-鐘擺，彈簧，樂器，機械，橋梁，電路. (1) 數(shù)學擺的無阻尼微小自由振動方程為 或 特征方程為2+2 = 0，特征根為共軛復根 1,2= i。 通解為 令 則通解式可改寫為 即,(1) 數(shù)學擺的無阻尼微小自由振動,從通解式 中可看出，不管初狀態(tài) 為何值， 擺的運動總一 個正弦 函數(shù)，且是周期 函數(shù)。 如圖（4.1） 這種運動稱為簡諧振動 稱為圓頻率，A為振幅， 為初相位。 周期為 ，頻率為,(1) 數(shù)學擺的無阻尼微小自由振動,數(shù)學擺的振幅A、初相位均依賴于初值條件。 數(shù)學擺的周期只依賴于擺長l,與初值無關。 當把數(shù)學擺移至位置= 0即 再松開時其初值條件為 將上式代入通解得 初相位 ，振幅A= 0 。 因此，所求的特解為,(2) 有阻尼自由振動,有阻尼自由振動方程為 特征方程為 有特征根 根據(jù) n 大于、小于和等于分為小阻尼、大阻尼和臨界阻尼三種情況進行討論。,(a) 小阻尼n時,小阻尼n時特征方程有共軛復根 通解為 和前面(1)一樣可將上通解改寫為 有小阻尼時擺的運 動已不是周期的， 其振幅逐漸減小。 按=Ae-nt的變化減小， 減小的周期為 當t時(t)擺動地 (振動著)趨于零。如圖,(b) 大阻尼n時,大阻尼時：這時特征方程有兩不同負實根210 通解為 從式中可知擺的運動也不是周期的， 且因方程 最多只有一個解，故擺最多只通過平衡位置一次； 同時由 知當t足夠大時， d/dt與c1的符號相反， 即經過一段時間后， 擺單調地趨于平衡 位置。如圖(4.3),(c) 臨界阻尼 n =時,臨界阻尼n =時：特征方程有重根2=1=-n 方程的通解為 擺的運動也不是周期的，其圖形和圖(4.3)類似，當t時(t)趨于零。擺不具有振動的性質。 值n =稱為阻尼的臨界值，正好可以抑制振動。 實際上，當n時如圖(4.3)，擺不振動； 而當n時如圖(4.2)，擺雖衰減但仍振動。,(3) 無阻尼強迫振動,無阻尼強迫振動方程為 令 其中H為已知常數(shù)，p為外力圓頻率， 此時上式變?yōu)?方程對應的齊次線性微分方程的通解為 這里A,為任意常數(shù)。下面求非齊方程的特解。,(3) (a)無阻尼p,(a) p 有形如 的解，這里M,N是待定常數(shù)。將上式代入方程，比較同類項系數(shù)，可得 最后得方程的通解為 通解中包括無阻尼自由振動的項(固有振動) 和外力引發(fā)的頻率相同、振幅不同的項(強迫振動) 當外力的圓頻率越接近固有圓頻率時，強迫振動項的振幅越大。,(3) (b)無阻尼p=,當p=時, 有形如 的解。將它代入方程，比較同類項系數(shù)，可得 方程的通解為 隨時間的增大，擺的偏離將無限增加， 這種現(xiàn)象稱為共振現(xiàn)象。 但實際上，當擺的偏離將無限增加到一定程度時，方程已不一定能描述擺的運動狀態(tài)了。,(4) 小阻尼n強迫振動,小阻尼n 強迫振動：方程為 由(2)(a)有阻尼自由振動情形知方程對應的齊次線性微分方程的通解為 這里A,為任意常數(shù)， 現(xiàn)求方程的一個特解。此特解有形狀 將上式代入方程，比較同類項系數(shù),(續(xù)) (4) 小阻尼n強迫振動,得到 如令 其中 則特解可寫為 因此程的通解為 其中,(續(xù)) (4) 小阻尼n強迫振動,通解中由有阻尼自由振動的固有振動項和外力引發(fā)的強迫振動項兩部分疊加而成。