數(shù)學(xué)物理方程第二章分離變量法word版_第1頁
數(shù)學(xué)物理方程第二章分離變量法word版_第2頁
數(shù)學(xué)物理方程第二章分離變量法word版_第3頁
數(shù)學(xué)物理方程第二章分離變量法word版_第4頁
數(shù)學(xué)物理方程第二章分離變量法word版_第5頁
已閱讀5頁,還剩22頁未讀 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、第五講補(bǔ)充常微分方程求解相關(guān)知識。第二章 分離變量法l 偏微分方程定解問題常用解法,分離變量法。l 解常微分方程定解問題時,通??偸窍惹蟪鑫⒎址匠痰奶亟?,由線性無關(guān)的特解疊加出通解,而后用定解條件定出疊加系數(shù)l 一階線性偏微分方程的求解問題,基本方法也是轉(zhuǎn)化為一階線性常微分方程組的求解問題l 對于二階以及更高階的偏微分方程定解問題,情況有些不同:即使可以先求出通解,由于通解中含有待定函數(shù),一般來說,很難直接根據(jù)定解條件定出,因此,通常的辦法就是把它轉(zhuǎn)化為常微分方程問題(第六講)§2.1 有界弦的自由振動什么是分離變量法?使用分離變量法應(yīng)具備那些條件? 下面通過兩端固定的弦的自由振動問

2、題來說明。定解問題:考慮長為l ,兩端固定的弦的自由振動,其數(shù)理方程及定解條件為分析:1. 方程和邊界條件都是齊次的,求這樣的問題可用疊加原理。2. 我們知道,在解常微分方程定解問題時,通常總是先求出微分方程的特解,由線性無關(guān)的特解疊加出通解,而后用定解條件定出疊加系數(shù)。啟發(fā):能否運(yùn)用類似求常微分方程定解問題的方法求偏微分方程?也既是能否先找出滿足齊次方程及齊次邊界條件的足夠多的特解,再用其作線性組合使其滿足初始條件。由分析,我們現(xiàn)在試求方程的變量分離形式:的非零解。l 將代入方程,可得此式中,左端是關(guān)于的函數(shù),右端是關(guān)于的函數(shù)。因此,左端和右端相等,就必須等于一個與無關(guān)的常數(shù)。設(shè)為,則有l(wèi)

3、將邊界條件代入得此時,必有這就完成了用分離變量法求解偏微分方程定解問題的第一步:分離變量l 目標(biāo):分離變量形式的解l 結(jié)果:得到函數(shù)滿足的常微分方程和邊界條件以及滿足的常微分方程,l 條件:偏微分方程和邊界條件都是齊次的現(xiàn)在我們求解函數(shù)滿足的常微分方程定解問題。我們發(fā)現(xiàn):方程中含有待定常數(shù),定解條件是齊次邊界條件,與一般的常微分方程的初值問題不同:l 并非對任一的,都有既滿足齊次方程有滿足邊界條件的非零解;l 只有當(dāng)取某些特定值時,才有既滿足方程又滿足邊界條件的非零解。有非零解的稱為該問題的特征值相應(yīng)的非零解稱特征函數(shù)而滿足的常微分方程的定解問題稱特征值問題。第二步:求解特征值問題1) 若,方

4、程的通解形式為由定解條件知,從而,不符合要求。2) 若,方程的通解形式為由邊界條件知,從而,不符合要求。3) 若,方程的通解形式為代入邊界條件得從而得特征值問題的一系列特征值及相應(yīng)的特征函數(shù)第三步:求特解,并疊加出一般解求解了特征值問題后,將每特征值代入函數(shù)滿足的方程可得出相應(yīng)的解因此,也就得到滿足偏微分方程和邊界條件的特解注:l 這樣的特解有無窮多個l 每個特解都滿足齊次方程和齊次邊界條件l 一般來說,單獨(dú)任何一個特解不可能恰好滿足定解問題的初始條件,即無法找到滿足。l 把全部特解疊加起來我們知,只要級數(shù)收斂,并且二次可微,則也滿足齊次邊值問題。下面選擇合適的使?jié)M足初始條件,即第四步:運(yùn)用特

5、征值函數(shù)的正交性定疊加系數(shù)事實(shí)上,我們知道補(bǔ)充內(nèi)容:f(x)的傅里葉級數(shù)其中總結(jié):利用分離變量法求解偏微分方程定解問題的基本步驟:第一步:分離變量這一步所以能夠?qū)崿F(xiàn),先決條件使偏微分方程和邊界條件都是齊次的。而分離變量的結(jié)果是得到含有待定常數(shù)的齊次常微分方程和齊次邊界條件,即特征值問題;第二步:求解特征值問題;第三步:求出全部特解,并進(jìn)一步疊加出一般解(形式解);第四步:利用特征函數(shù)的正交性確定疊加系數(shù)。(第七講)嚴(yán)格來說,上面得到的還是形式解,對于具體問題,還必須驗(yàn)證:1) 這樣得到的是否滿足偏微分方程,換句話說,級數(shù)解是否可以逐項(xiàng)求二階導(dǎo)數(shù);2) 是否滿足邊界條件,即是否連續(xù);3) 確定系

6、數(shù)時,逐項(xiàng)積分是否合理。關(guān)于上三個問題,都涉及到級數(shù)解的收斂性,由于系數(shù)都是由決定的,因而的性質(zhì)就決定了上三個問題的回答??梢宰C明若三次可微,二次可微,則問題解存在,且此解可用上面的級數(shù)形式給出(見復(fù)旦大學(xué)數(shù)學(xué)物理方程)。從理論上講,分離變量法之所以成功,要取決于下列幾個條件:1) 特征值問題有解;2) 定解問題的解一定可以按照特征值函數(shù)展開,也既是說,特征值函數(shù)是完備的;3) 特征值函數(shù)一定具有正交性。以后適當(dāng)回答這些問題。解的物理意義先看特解其中。l 代表一個駐波駐波:頻率和振幅均相同、振動方向一致、傳播方向相反的兩列波疊加后形成的波。波在介質(zhì)中傳播時其波形不斷向前推進(jìn),故稱行波;上述兩列

7、波疊加后波形并不向前推進(jìn),故稱駐波。l 表示弦上各點(diǎn)的振幅分布l 是振動的固有頻率,稱為弦的固有頻率或特征(本征)頻率l 為初相位,由初始條件決定l 在的各點(diǎn)上,振動振幅恒為零,稱為波節(jié)(節(jié)點(diǎn)),包含兩個端點(diǎn)共有個節(jié)點(diǎn)l 在的各點(diǎn)上,振幅絕對值恒為最大,稱為波峰(腹點(diǎn)),共有n個l 滿足定解問題的級數(shù)解則是這些駐波的疊加,因此也稱分離變量法為駐波法l 就兩端點(diǎn)固定的弦來說,固有頻率中有一最小值,即,稱為基頻,其他頻率都是其倍數(shù),稱為倍頻。實(shí)際例子:l 弦的基頻決定了聲音的音調(diào)。在弦樂器中,當(dāng)弦的質(zhì)料一定時,可以通過改變弦的繃緊程度,調(diào)解的大小。l 的相對大小,決定了聲音的頻譜分布,即決定了音色

8、。l 和數(shù)與弦的能量成正比,決定了聲音的強(qiáng)度。分離變量法舉例例題1. 有一根長為10個單位的弦,兩端固定,初速為零,初位移為,求弦作微小橫振動的位移。解:設(shè)位移為,它的定解問題的解。給定,顯然,這個問題的傅立葉級數(shù)解可由給出,其系數(shù)為因此,所求的解為例題2. 解定解問題解:運(yùn)用分離變量可得將邊界條件代入可得相應(yīng)的特征值問題重復(fù)前面的解法,知當(dāng)時,特征值問題有解,此時通解形式為代入邊界條件得從而求得一系列特征值和特征函數(shù)與這些特征值相對應(yīng)得的通解表示為于是,所求定解問題的形式解可表示為利用初始條件確定其中的系數(shù)得故所求的解為(第八講)§2.2 有限桿上的熱傳導(dǎo)定解問題:一均勻細(xì)桿,長為

9、,兩端坐標(biāo)為。桿的側(cè)面絕熱,且在端點(diǎn)處溫度為零,而在處桿的熱量自由發(fā)散到周圍溫度為0的介質(zhì)中。初始溫度為,求桿上的溫度變化情況,即考慮下定解問題:仍用分離變量法求解。此定解問題的邊界條件為第三類邊界條件。類似§2.1中步驟,設(shè),代入上面的方程可得從而可得通解由邊界條件知從而令 上方程的解可以看作曲線,交點(diǎn)的橫坐標(biāo),顯然他們有無窮多個,于是方程有無窮多個根。用下符號表示其無窮多個正根于是得到特征值問題的無窮個特征值及相應(yīng)的特征函數(shù)再由方程, 可得,從而我們得到滿足邊界條件的一組特解由于方程和邊界條件是齊次的,所以仍滿足此方程和邊界條件。下面研究一下其是否滿足初始條件??梢宰C明在區(qū)域0,

10、l上具有正交性,即證明:完成。令于是,從而得到定解問題得解。§2.3 圓域內(nèi)的二維Laplace方程的定解問題平面極坐標(biāo)和直角坐標(biāo)的關(guān)系是由此可得即是由復(fù)合函數(shù)求導(dǎo)法則,可得進(jìn)一步,可得在此基礎(chǔ)上,還可以得到柱坐標(biāo)系下的Laplace算符(第九講)考慮圓域內(nèi)的穩(wěn)定問題:其在極坐標(biāo)下的表示形式:因圓域內(nèi)溫度不可能為無限,尤其是在圓盤中心點(diǎn)的溫度應(yīng)該有限,并且表示同一點(diǎn),故而我們有下約束下面用分離變量法求解該問題。令 代入極坐標(biāo)下方程可得:從而可得常微分方程由有限性及周期邊界條件知,從而得定解問題求解:時,通解為由周期邊界條件可得 從而,不可取。時,通解為由周期邊界條件可得 B任意,說明

11、為一特征值,相應(yīng)得特征函數(shù)為。時,通解為因以為周期,所以有 從而可得特征值特征函數(shù)為接下來,求特解,并疊加出一般解。由Euler方程若令,即,則上方程可寫為故時,通解時,通解為為保證,所以可得,即從而,滿足齊次方程和周期條件及有限性的解可以表示為級數(shù)最后,為了確定系數(shù),我們利用邊界條件可得運(yùn)用性質(zhì)從而可得因而,我們有利用下面的求和公式所以,稱此表達(dá)式為圓域內(nèi)的Poisson公式,它的作用是把解寫成積分形式,便于作理論上的研究。例題 解下列定解問題:解:利用公式可知, 所以。(第十講)Laplace變換法定義:函數(shù)當(dāng)時有定義,當(dāng)s屬于某區(qū)間內(nèi)時廣義積分收斂,則由此積分確定的函數(shù)稱為的Laplac

12、e變換。記,同時把稱為的逆Laplace變換,記為。性質(zhì):1 線性性質(zhì):a,b為常數(shù),則對,的拉氏變換同時存在的s有2 若則3 微分性質(zhì):設(shè)在區(qū)間上連續(xù),其導(dǎo)數(shù)是分段連續(xù)的,且存在常數(shù)M,使得 ,對于。則對任,拉氏變換存在,且 。更一般的形式 4 積分性質(zhì):設(shè)在區(qū)間上連續(xù),其導(dǎo)數(shù)是分段連續(xù)的,且存在常數(shù)M,使得 ,對于。則對任,則有5 設(shè)在區(qū)間上連續(xù),其導(dǎo)數(shù)是分段連續(xù)的,且存在常數(shù)M,使得 ,對于,。則對任,則有,更一般的是6 設(shè)在區(qū)間上分段連續(xù)的,且存在常數(shù)M,使得 ,對于,且存在,則對任,7 卷積性質(zhì):設(shè),在區(qū)間上分段連續(xù)的,且存在常數(shù)M,使得 ,則對任,卷積的拉氏變換存在且8 平移性質(zhì):

13、若對存在,則對于,則其中 為Heavivide函數(shù)。一些常見函數(shù)的Laplace變換:1);2) ;3) ;4) ;5) ;6) 7) 8) 9) §2.4 非齊次方程的解法分離變量法成功的關(guān)鍵是:方程和邊界條件都是齊次的。若方程非齊次的,邊界條件為齊次的,能否運(yùn)用分離變量法?若能,如何求解?下面以弦的強(qiáng)迫振動為例,來討論非齊次方程的解法。問題模型:兩端固定的弦,受強(qiáng)迫力作用產(chǎn)生振動現(xiàn)象。其定解問題如下: (1)此情況下,弦的振動是由兩部分干擾引起的,一是強(qiáng)迫力,一是初始狀態(tài),因此,振動可以看作為僅由強(qiáng)迫力引起的振動和僅由初始狀態(tài)引起的振動。啟發(fā):可設(shè)解為其中表示僅強(qiáng)迫振動的位移,滿

14、足 (2)而表示純初始狀態(tài)引起的振動位移,滿足 (3)注:不難驗(yàn)證,只要為(2)的解,為(3)的解,則必為(1)的解。問題(3)由分離變量法很容易求解,因此(2)的求解是關(guān)鍵。下面討論如何求解純強(qiáng)迫振動問題的解。思路:類似于線性非齊次常微分方程中常用的常數(shù)變易法(補(bǔ)充線性非齊次常微分方程的求解),希望問題(2)的解可分解為無窮多個駐波的疊加,而每個駐波的波形由相應(yīng)的齊次方程通過分離變量所得的特征值問題的特征函數(shù)所決定,既解有形式:(運(yùn)用對應(yīng)齊次問題的特解的線性組合去構(gòu)造非齊次問題的特解)(4)其中為待定系數(shù)。為確定,我們將自由項(xiàng)也按特征函數(shù)展開 (5)將(4)和(5)代入問題(2)的方程得由此

15、得再將(4)代入初始條件得因此確定只需求解如下定解問題 (6)用Laplace變換法(或參數(shù)變易法)解出(6),得所以,將此解與問題(3)的解加起來就得定解問題(1)的解。注:l 求解上面的非齊次方程的方法稱為特征函數(shù)法,其實(shí)質(zhì)就是將方程自由項(xiàng)及解都按齊次方程所對應(yīng)的一族特征函數(shù)展開。l 該方法也適用于其他類的齊次邊界問題。(第十一講)例:在環(huán)型域內(nèi)求解下定解問題 (7)解:因解域?yàn)榄h(huán)形區(qū)域,故可選平面極坐標(biāo)系,利用平面極坐標(biāo)和直角坐標(biāo)的關(guān)系則上問題可表示為 (8)此方程為非齊次方程附有齊次邊界條件。用特征函數(shù)法求解。由§2.3節(jié)中得到關(guān)于圓域內(nèi)laplace方程對應(yīng)的特征函數(shù),可令

16、問題(8)的解為代入(8)中方程可得比較兩端系數(shù)可得(9)再由(8)的條件可得 (10)(9)中第二、三個方程都是齊次的Euler方程,其通解為由邊界條件(10)可得 下面求。因(9)第一方程為非齊次Euler方程,首先用待定系數(shù)法可求的其一特解,從而,他的通解為,由條件(10)第一式,可得所以故。(第十二講)§2.5 非齊次邊界條件的處理前面討論的定解問題,無論方程是齊次還是非齊次的,邊界條件都是齊次的,對于邊界條件為非齊次的如何處理?原理:將邊界齊次化。下以波動方程的定解問題為例說明。考慮如下定解問題: (1)為應(yīng)用分離變量法,設(shè)法作一代換將邊界條件齊次化,為此令 (2)適當(dāng)選擇,使得的邊界條件化為齊次的,即 (3)因此,由(1)和(2)可知,滿足 (4)所以只要找到滿足(4)就達(dá)到目的。對于這樣的有很多選法,例如取其為一次式,即要滿足(4),則可得從而這樣只要作代換 (5)就能使得滿足齊次邊界條件,即有 (6)顯然(6)可以用特征函數(shù)法求解。注:1 上面選擇一次式,是因?yàn)楹唵我浊?,可以有別的形式;2 若f,u1,u2與t無關(guān),則可選適當(dāng)?shù)氖沟脻M足的方程和邊界條件都是齊次的,減少求解的工作量。(見例1)3 此方法對其他類邊界條件依然成立,只不過是表達(dá)式不同而

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論