第二章分離變量3ppt課件

1、2.4 非齊次方程的解法非齊次方程的解法 通過前面課程的學(xué)習(xí)，我們已經(jīng)了解，用分離變量法求通過前面課程的學(xué)習(xí)，我們已經(jīng)了解，用分離變量法求解偏微分方程定解問題，這個(gè)定解問題必須是線性、齊解偏微分方程定解問題，這個(gè)定解問題必須是線性、齊次方程、齊次邊界條件。那么對(duì)于非齊次方程和非齊次次方程、齊次邊界條件。那么對(duì)于非齊次方程和非齊次邊界條件如何進(jìn)行處理？邊界條件如何進(jìn)行處理？非齊次方程、齊次邊界條件非齊次方程、齊次邊界條件考慮如下定解問題：考慮如下定解問題：22222( , ),0,0 (2.37)uuaf x tx l ttx 00,0 (2.38)xx luut00( ),( ),0 (2.3

2、9)ttuuxxxlt 從物理上看：從物理上看：在現(xiàn)在的情況，弦的振動(dòng)是由兩部分干擾引起的，一是在現(xiàn)在的情況，弦的振動(dòng)是由兩部分干擾引起的，一是強(qiáng)迫力，一是初始狀態(tài)，所以由物理意義可知，此時(shí)的強(qiáng)迫力，一是初始狀態(tài)，所以由物理意義可知，此時(shí)的振動(dòng)可以看作為僅由強(qiáng)迫力引起的振動(dòng)和僅由初始狀態(tài)振動(dòng)可以看作為僅由強(qiáng)迫力引起的振動(dòng)和僅由初始狀態(tài)引起的振動(dòng)的合成。引起的振動(dòng)的合成。從數(shù)學(xué)上看：從數(shù)學(xué)上看：就是將將弦的振動(dòng)位移分解為強(qiáng)迫力引起振動(dòng)的位移與就是將將弦的振動(dòng)位移分解為強(qiáng)迫力引起振動(dòng)的位移與初始狀態(tài)引起振動(dòng)的位移的和。即初始狀態(tài)引起振動(dòng)的位移的和。即( , )V x t設(shè)設(shè)強(qiáng)迫力引起振動(dòng)的位移；強(qiáng)

3、迫力引起振動(dòng)的位移；( , )W x t初始狀態(tài)引起振動(dòng)的位移。初始狀態(tài)引起振動(dòng)的位移。即即( , )( , )( , )u x tV x tW x t（2.40） ( , )V x t滿足如下定解問題：滿足如下定解問題：22222( , ),0,0VVaf x txl ttx00,0 xx lVVt000,0,0ttVVxlt（2.41） 22222,0,0WWaxl ttx00,0 xx lWWt00( ),( ),0ttWWxxxlt（2.42） ( , )W x t滿足如下定解問題：滿足如下定解問題：22222000+ + ( , ),0,0+ + 0,0+ ,+ ,0 xx lttV

4、Vaf x txl ttxVVtVVxlt22222000 ( ) ( )xx lttWWatxWWWWxxt22222000( , ),0,00,00,0,0 xx lttVVaf x txl ttxVVtVVxlt22222000,0,00,0( ),( ),0 xx lttWWaxl ttxWWtWWxxxlt22222000( , ),0,00,00,0,0 xx lttVVaf x txl ttxVVtVVxlt22222000,0,00,0( ),( ),0 xx lttWWaxl ttxWWtWWxxxlt( , )( , )22222()()( , ),0,0u x tu x

5、tVWVWaf x txl ttx ( , )( , )0()()0,0u x tu x txx lVWVWt ( , )( , )00()()( ),( ),0u x tu x tttVWVWxxxlt 不難驗(yàn)證，若不難驗(yàn)證，若V是是2.41的解，的解，W是是2.42的解，那么的解，那么uV+W一定就是原定解問題的解。一定就是原定解問題的解。00,0 (2.38)xx luut00( ),( ),0 (2.39)ttuuxxxlt 22222( , ),0,0 (2.37)uuaf x tx l ttx 問題問題2.42可以直接用分離變量法求解，因此現(xiàn)在的可以直接用分離變量法求解，因此現(xiàn)在的

6、問題只要討論如何解問題問題只要討論如何解問題2.41就行了。就行了。22222,0,0WWaxl ttx00,0 xx lWWt00( ),( ),0ttWWxxxlt（2.42） 22222( , ),0,0VVaf x txl ttx00,0 xx lVVt000,0,0ttVVxlt（2.41） 復(fù)習(xí)：參數(shù)常數(shù)變易法復(fù)習(xí)：參數(shù)常數(shù)變易法下列形式的一階常微分方程下列形式的一階常微分方程假設(shè)假設(shè)( )0Q x ，則方程成為，則方程成為( )( )yP x yQ x一階線性非齊次一階線性非齊次( )0yP x y一階線性齊次一階線性齊次一階線性齊次常微分方程的解法一階線性齊次常微分方程的解法(

7、 )0yP x y是可分離變量方程。分離變量，得是可分離變量方程。分離變量，得( )dyP x dxy ( )dyP x dxy 兩邊積分，得兩邊積分，得ln( )lnyP x dxC所以，方程的通解公式為所以，方程的通解公式為( )P x dxyCe一階線性非齊次常微分方程的解法一階線性非齊次常微分方程的解法齊次方程與非齊次方程的差異，在于齊次方程與非齊次方程的差異，在于( )0Q x 。因而，我們可以設(shè)想它們的通解之間會(huì)有一定的聯(lián)系。因而，我們可以設(shè)想它們的通解之間會(huì)有一定的聯(lián)系。設(shè)設(shè)1( )yy x，是齊次方程的一個(gè)解，則當(dāng)，是齊次方程的一個(gè)解，則當(dāng)C為常數(shù)時(shí)，為常數(shù)時(shí)，1( )yCy

8、x1( )yy x仍然是仍然是( )0yP x y的一個(gè)解。它不可能滿足線性非齊次方程。如果我們把的一個(gè)解。它不可能滿足線性非齊次方程。如果我們把C看作看作x的函數(shù)，并將的函數(shù)，并將1( )( )yC x y x代入線性非齊次方程中，得代入線性非齊次方程中，得設(shè)設(shè)1( )( )yC x y x是非齊次方程的解，將是非齊次方程的解，將1( )( )yC x y x及其導(dǎo)數(shù)及其導(dǎo)數(shù)11( )( )( )yC x y xC x y代入方程代入方程( )( )yP x yQ x常數(shù)變易常數(shù)變易得得111( )( ) ( ) ( ) ( )yyC x yC x yP x C x yQ x 111( )(

9、 ) ( ) ( ) ( )yyC x yC x yP x C x yQ x 即即111( )( )( ) ( )C x yC x yP x yQ x因因1( )yy x是對(duì)應(yīng)線性齊次方程的解，所以是對(duì)應(yīng)線性齊次方程的解，所以11( )0yP x y因而，有因而，有1( )( )C x yQ x1( )( )C xQ x y兩邊積分，得兩邊積分，得1( )( )C xQ x ydx C1( )( )C xQ x ydx C1( )( )yC x y x將將代入到代入到得得111( )( )yCy xyQ xy dx在實(shí)際運(yùn)算中，我們?nèi)≡趯?shí)際運(yùn)算中，我們?nèi)? )1P x dxye于是，得到一階線

10、性非齊次常微分方程的通解公式：于是，得到一階線性非齊次常微分方程的通解公式：( )( )1( )P x dxP x dxyeCyQ x edx22222000,0,0 (2.1)0,0,0 (2.2)( ),( ),0 (2.3)xx lttuuaxl ttxuutuuxxxlt齊次方程、齊次邊界條件下，定解問題齊次方程、齊次邊界條件下，定解問題的解為的解為11( , )( , ) (cossin)sinnnnnnu x tu x tn an anCtDtxlll11( )1( , )( , ) (cossin)sin ( )sinnnnnnnutnnu x tu x tn an anCtDt

11、xlllnu txl 22222000( , ),0,00,00,0,0 xx lttVVaf x txl ttxVVtVVxlt而定解問題而定解問題（2.41） 只比齊次方程多一個(gè)自由項(xiàng)只比齊次方程多一個(gè)自由項(xiàng)f(x,t)，所以設(shè)想，所以設(shè)想2.41的解的解22222000,0,0 (2.1)0,0,0 (2.2)( ),( ),0 (2.3)xx lttuuaxl ttxuutuuxxxlt有如下形式：有如下形式：1( , )( )sinnnnV x tv txl（2.43） ( )( )yP x yQ x( )0yP x y一階線性非齊次一階線性非齊次一階線性齊次一階線性齊次1( )yy

12、 x是一個(gè)特解是一個(gè)特解1( )yCy x也滿足方程也滿足方程1( )( )yC x y x是非齊次方程的解是非齊次方程的解1( )( )C xQ x y1( )( )C xQ x ydx C其中其中( )nv t是待定的函數(shù)。是待定的函數(shù)。也按特征函數(shù)系展開成如下的也按特征函數(shù)系展開成如下的級(jí)數(shù)：級(jí)數(shù)： 1( , )( )sinnnnf x tf txl（2.44） 其中其中 為了確定為了確定( )nv t，將自由項(xiàng)，將自由項(xiàng)f(x,t)02( )( , )sinlnnf tf x txdxll1( , )( )sinnnnV x tv txl（2.43） 22222( , ),0,0VVa

13、f x txl ttx將將2.43及及2.44代入代入 得得22221( )( )( )sin0nnnna nnv tv tf txll（2.45） （2.44） 1( , )( )sinnnnf x tf txl22221( )( )( )sin0nnnna nnv tv tf txll（2.45） 2222( )( )( )0nnna nv tv tf tl即即2222( )( )( )nnna nv tv tf tl我們又得到一個(gè)常微分方程。我們又得到一個(gè)常微分方程。因?yàn)橐驗(yàn)閟innxl是線性無關(guān)的，所以只有是線性無關(guān)的，所以只有再將再將2.43代入代入2.41中的初始條件得中的初始條件

14、得 1( , )( )sinnnnV x tv txl（2.43） 22222( , ),0,0VVaf x txl ttx00,0 xx lVVt000,0,0ttVVxlt（2.41） (0)0,(0)0nnvv同樣可得同樣可得因而，只需解如下的常微分方程初值問題：因而，只需解如下的常微分方程初值問題：2222( )( )( )nnna nv tv tf tl(0)0,(0)0,(1,2,)nnvvn（2.46） 用拉普拉斯變換法解這個(gè)二階常系數(shù)非齊次常微分方程用拉普拉斯變換法解這個(gè)二階常系數(shù)非齊次常微分方程 拉普拉斯變換理論又稱為運(yùn)算微積分，或稱為算子微積分是在19世紀(jì)末發(fā)展起來的首先是

15、英國工程師亥維賽德(O.Heaviside)發(fā)明了用運(yùn)算法解決當(dāng)時(shí)電工計(jì)算中出現(xiàn)的一些問題，但是缺乏嚴(yán)密的數(shù)學(xué)論證后來由法國數(shù)學(xué)家拉普拉斯(P.S.Laplace)給出了嚴(yán)密的數(shù)學(xué)定義，稱之為拉普拉斯變換方法 解：解：在方程在方程2.46）2222( )( )( )nnna nv tv tf tl(0)0,(0)0,(1,2,)nnvvn得兩端取關(guān)于得兩端取關(guān)于t的拉普拉斯變換，得的拉普拉斯變換，得22222( )( )( )nnna np UpUpFpl( )nUp其中其中和和( )nFp分別是分別是( ),( )nnv tf t得拉普拉斯變換得拉普拉斯變換（2.46）222221( )(

16、)nnUpFpa npl由于由于222221a npl的逆拉普拉斯變化為的逆拉普拉斯變化為1sinn atn al利用拉普拉斯變化的卷積性質(zhì)，得利用拉普拉斯變化的卷積性質(zhì)，得01()( )( )sintnnn a tv tfdn al所以，所以，011()( , )( )sinsintnnn a tnV x tfdxn all22222,0,0WWaxl ttx00,0 xx lWWt00( ),( ),0ttWWxxxlt（2.42） 將這個(gè)解與將這個(gè)解與2.42）22222( , ),0,0 (2.37)uuaf x tx l ttx 00,0 (2.38)xx luut00( ),( )

17、,0 (2.39)ttuuxxxlt 的解加起來，就得到原定解問題的解加起來，就得到原定解問題2.372.39的解。的解。還可以用沖量定理法求解非齊次振動(dòng)方程定解問題。還可以用沖量定理法求解非齊次振動(dòng)方程定解問題。22222( , ),0,0VVaf x txl ttx00,0 xx lVVt000,0,0ttVVxlt（2.41） 沖量定理法的前提是初始條件均取零值。沖量定理法的前提是初始條件均取零值。這里所給的求解問題這里所給的求解問題2.41的方法，其實(shí)質(zhì)是將方程的方法，其實(shí)質(zhì)是將方程的自由項(xiàng)及解都按齊次方程所對(duì)應(yīng)的一族特征函數(shù)展開的自由項(xiàng)及解都按齊次方程所對(duì)應(yīng)的一族特征函數(shù)展開。隨著方

18、程與邊界條件的不同，特征函數(shù)族也就不同，。隨著方程與邊界條件的不同，特征函數(shù)族也就不同，但總是把非齊次方程的解按相應(yīng)的特征函數(shù)展開。所但總是把非齊次方程的解按相應(yīng)的特征函數(shù)展開。所以這種方法也叫做特征函數(shù)法。以這種方法也叫做特征函數(shù)法。 22222( , ),0,0VVaf x txl ttx00,0 xx lVVt000,0,0ttVVxlt（2.41） 解解 由于求解區(qū)域是環(huán)形區(qū)域，所以我們選用平面極坐由于求解區(qū)域是環(huán)形區(qū)域，所以我們選用平面極坐標(biāo)系，利用直角坐標(biāo)系與極坐標(biāo)系之間的關(guān)系標(biāo)系，利用直角坐標(biāo)系與極坐標(biāo)系之間的關(guān)系cossinxy可將上述定解問題用極坐標(biāo)可將上述定解問題用極坐標(biāo)

19、, 形式：形式： 222211()12 cos2 ,02 (2.47)0,0,02 (2.48)abuuabuu 這是一個(gè)非齊次方程附有齊次邊界條件的定解問題。采這是一個(gè)非齊次方程附有齊次邊界條件的定解問題。采用特征函數(shù)法，并注意到在用特征函數(shù)法，并注意到在2.3中得到的關(guān)于圓域內(nèi)拉中得到的關(guān)于圓域內(nèi)拉普拉斯方程所對(duì)應(yīng)的特征函數(shù)，可令問題普拉斯方程所對(duì)應(yīng)的特征函數(shù)，可令問題2.47）、）、（2.48的解為的解為0( , )( )cos( )sinnnnuAnBn 代入代入2.47）2202221( )( )( )cos1( )( )( )sin12cos2nnnnnnnnAAAnnBBBn2222214( )( )( )12AAA（2.49） 221( )( )( )0(2)nnnnAAAn（2.50） 222211()12cos2uu并整理得到并整理得到0( , )( )cos( )sinnnnuAnBn 221(