數(shù)學(xué)物理方法傅里葉變換法

1、積分變換在數(shù)學(xué)物理方程中也有廣泛的用途，變換后，積分變換在數(shù)學(xué)物理方程中也有廣泛的用途，變換后，方程得以化簡，偏微分方程變成常微分方程，求解常微方程得以化簡，偏微分方程變成常微分方程，求解常微方程后，再進(jìn)行逆變換就得到原來偏微分方程的解，方程后，再進(jìn)行逆變換就得到原來偏微分方程的解，同時，積分變換還可能得到有限形式的解，分離變數(shù)同時，積分變換還可能得到有限形式的解，分離變數(shù)法或者傅里葉級數(shù)發(fā)往往不能。法或者傅里葉級數(shù)發(fā)往往不能。本章主要介紹傅里葉變換法在求解偏微分方程中的應(yīng)用。本章主要介紹傅里葉變換法在求解偏微分方程中的應(yīng)用。 itfF tfF)(Fii1 dttftF tfF)(1aFa)(

2、Faxf)(F00Fexxfxi )(F00Fxfexi ,2F2121FFxfxf用分離變數(shù)法求解有界空間的定解問題時，得到的本征值是用分離變數(shù)法求解有界空間的定解問題時，得到的本征值是求解無限長弦的自由振動求解無限長弦的自由振動)(|),(|)(0002xuxuxuautttxxtt應(yīng)用傅里葉變換，即用應(yīng)用傅里葉變換，即用2/ikxe同乘方程和定解條件同乘方程和定解條件中的各項，并對空間變量中的各項，并對空間變量x積分，積分，t看做參數(shù)，則看做參數(shù)，則分分，對于無界空間的定解問題，適用于傅里葉變換法求解。，對于無界空間的定解問題，適用于傅里葉變換法求解。連續(xù)連續(xù)的，所求的解可表示為對連續(xù)本

3、征值求積分的的，所求的解可表示為對連續(xù)本征值求積分的傅里葉積傅里葉積無界空間，分離變數(shù)法求解定解問題時，所得到的本征值是無界空間，分離變數(shù)法求解定解問題時，所得到的本征值是離散的，所求的解可表為對本征值求和的離散的，所求的解可表為對本征值求和的傅里葉級數(shù)傅里葉級數(shù)，對于，對于定解問題變換成：定解問題變換成：)(|),(|00022kUkUUakUtt 其中其中)()(kk，分別是分別是)(),(xx的傅里葉變換，這樣原來的傅里葉變換，這樣原來的定解問題變成了常微分方程及初值條件，通解為：的定解問題變成了常微分方程及初值條件，通解為：ikatikatekBekAktU)()(),(代入初始條件可

4、得：代入初始條件可得：)(121)(21)()(121)(21)(kikakkBkikakkA故故ikatikatikatikatekikaekekikaekktU)(121)(21 )(121)(21),(對對U作逆傅里葉變換，可得最后的結(jié)果如下：作逆傅里葉變換，可得最后的結(jié)果如下： atxatxdaatxatxtxu )(21)()(21),(求解無限長細(xì)桿的熱傳導(dǎo)問題求解無限長細(xì)桿的熱傳導(dǎo)問題)(|)(002xuxuautxxt作傅里葉變換，定解問題變?yōu)椋鹤鞲道锶~變換，定解問題變?yōu)椋?(|0022kUUakUt此常微分方程的初始問題的解為此常微分方程的初始問題的解為takekktU22)

5、(),(進(jìn)行傅里葉逆變換可得：進(jìn)行傅里葉逆變換可得：dkeededkeekktUFtxuikxtakikikxtak2222)(21 )(),(),(1 交換積分次序交換積分次序ddkeetxuxiktak)(21 ),()(22積分公式：積分公式：22224/)/(eadkeekk222222222222222222222224404040)2(4)2(422222eedxeedeedkeedkeedkexkkkk求解無限長細(xì)桿的有源熱傳導(dǎo)問題求解無限長細(xì)桿的有源熱傳導(dǎo)問題0|)(,(02txxtuxtxfuau作傅里葉變換，定解問題變?yōu)榉驱R次常微分方程：作傅里葉變換，定解問題變?yōu)榉驱R次常微

6、分方程：0|);(022tUktFUakU令令)(,xita利用上述公式可得利用上述公式可得detatxutax224)(21)( ),(用用take22同乘方程各項，可得：同乘方程各項，可得：taktakektFektUdtd2222);(),(對對t積分一次，并考慮零初始值可得：積分一次，并考慮零初始值可得： taktakiktaktakddeeefdekFektU0022222222),( );();(進(jìn)行傅里葉逆變換進(jìn)行傅里葉逆變換dkeddeeftxuikxtiktak 0)(22),(21),(交換積分次序可得：交換積分次序可得： txiktakdddkeeftxu0)()(222

7、1),(),( 是單位面積硅片是單位面積硅片表層原有雜質(zhì)總量表層原有雜質(zhì)總量.并利用積分公式可得最后的結(jié)果為：并利用積分公式可得最后的結(jié)果為： ttaxddetaftxu0)(4)(22)(21),(),(在半導(dǎo)體擴散工藝中，雜質(zhì)擴散深度遠(yuǎn)遠(yuǎn)小于硅片厚度，可在半導(dǎo)體擴散工藝中，雜質(zhì)擴散深度遠(yuǎn)遠(yuǎn)小于硅片厚度，可)0( )0(|0|00002xxuuuautxxxxt0硅片，這里求解的是半無界空間硅片，這里求解的是半無界空間x0中的定解問題中的定解問題:有的雜質(zhì)向硅片內(nèi)擴散，但不讓新的雜質(zhì)穿過硅片表面進(jìn)入有的雜質(zhì)向硅片內(nèi)擴散，但不讓新的雜質(zhì)穿過硅片表面進(jìn)入以把硅片看成無限厚，在限定源擴散中，是只讓

8、硅片表層已以把硅片看成無限厚，在限定源擴散中，是只讓硅片表層已沒有雜質(zhì)穿過硅片表面沒有雜質(zhì)穿過硅片表面,即即:0|0 xxu第二類齊次邊界條件第二類齊次邊界條件這種邊界條件意味著偶延拓這種邊界條件意味著偶延拓,即求解以下定解問題即求解以下定解問題)0( )0()0( )0(|0|000002xxxxuuuautxxxxt)(2|00 xut則則)x)(-(2|0002xuuautxxt引用例引用例2結(jié)果可得結(jié)果可得taxtaxetadetatxu22224/04)(022 21)(2),(),(txuxO硅片表面硅片表面123右圖描述了雜質(zhì)濃度右圖描述了雜質(zhì)濃度u(x,t)在硅片中在硅片中0即

9、說明雜質(zhì)總量不變即說明雜質(zhì)總量不變,曲線跟縱軸相交處的切線都是水平的曲線跟縱軸相交處的切線都是水平的,在恒定表面濃度擴散中在恒定表面濃度擴散中,包圍硅片氣體包圍硅片氣體中含有大量的雜質(zhì)原子中含有大量的雜質(zhì)原子,源源不斷穿過硅片表面向內(nèi)部擴散源源不斷穿過硅片表面向內(nèi)部擴散,由由0|00002txxxtuNuuau即硅片表面的濃度梯度為零即硅片表面的濃度梯度為零,表明沒有新的雜質(zhì)進(jìn)入硅片表明沒有新的雜質(zhì)進(jìn)入硅片.度趨于均勻度趨于均勻,曲線下的面積為曲線下的面積為2,3依次對應(yīng)越來越晚的時刻依次對應(yīng)越來越晚的時刻,雜質(zhì)濃雜質(zhì)濃的分布情況的分布情況,曲線曲線1對應(yīng)于較早的時刻對應(yīng)于較早的時刻是半無界空

10、間是半無界空間x0中的定解問題中的定解問題于雜質(zhì)分子充足于雜質(zhì)分子充足,硅片表面雜質(zhì)濃度保持某個常數(shù)硅片表面雜質(zhì)濃度保持某個常數(shù)N0,這里所求這里所求首先把非齊次邊界條件化為齊次邊界條件首先把非齊次邊界條件化為齊次邊界條件,令令),(),(0txwNtxu則化為關(guān)于則化為關(guān)于w的定解問題的定解問題:00000002|0|0NNuwNuwwawttxxxxt這是第一類齊次邊界條件這是第一類齊次邊界條件,意味著奇延拓意味著奇延拓,即即)0()0(|00002xNxNwwawtxxt引用例引用例2結(jié)果可得結(jié)果可得04)(004)(022222121),(detaNdetaNtxwtaxtax第一個積

11、分中令第一個積分中令taddztaxz2/,2/ )(第二個積分中令第二個積分中令taddztaxz2/,2/ )(則有則有taxtaxztaxztaxzdzeNdzeNdzeNtxw2/2/02/02/0222 ),(被積函數(shù)是偶函數(shù)被積函數(shù)是偶函數(shù),故故taxzdzeNtxw2/0022),(記做記做erfx,則則w可寫為可寫為:)2(),(0taxerfNtxw所求的解如下所求的解如下:)2(1),(),(00taxerfNtxwNtxu記做記做erfcx,則有則有taxerfcNtxu2),(0),(txutxO硅片表面硅片表面123右圖描述了雜質(zhì)濃度右圖描述了雜質(zhì)濃度u(x,t)在硅

12、片中在硅片中求解三維無界空間中的波動問題求解三維無界空間中的波動問題)(|),(|00032rUruuauttttt明顯明顯,如果擴散持續(xù)進(jìn)行下去如果擴散持續(xù)進(jìn)行下去,則濃度分布最終將為常數(shù)則濃度分布最終將為常數(shù)N0(虛線虛線)的時刻的時刻,雜質(zhì)濃度趨于均勻的趨勢很雜質(zhì)濃度趨于均勻的趨勢很刻刻,2對應(yīng)于較晚的時刻對應(yīng)于較晚的時刻,3對應(yīng)于更晚對應(yīng)于更晚分布情況分布情況,曲線曲線1對應(yīng)于某個較早的時對應(yīng)于某個較早的時做傅里葉變換做傅里葉變換,問題變換為問題變換為的初始值問題的初始值問題 )(|),(|00022kUkUUakUtt這個方程的解為這個方程的解為)(121)(21),(ikatika

13、tikatikateekikaeekktU再進(jìn)行傅里葉逆變換再進(jìn)行傅里葉逆變換321)(121)( )(21)(),(dkdkdkeeeikakeektrUrikikatikatikatikatVddkdkdkeeeararrikikatikat)(4)(41321)(2Vddkdkdkeeeikrarrikikatikat)(141)(41321)(2Vddkdkdkeeeikrtarrikikatikat)(141)(41321)(2Vddkdkdkeeeikrarrikikatikat)(141)(41321)(2利用利用5.3例例1的結(jié)果的結(jié)果)(41),(rtatrUVddkdkdk

14、eatrrrrik321)()(1)(41raVddkdkdkeatrrrrik321)()(1應(yīng)用延遲定理應(yīng)用延遲定理VdatrrrrrtatrU)| (|)(41),(Vdatrrrrra)| (|)(41出現(xiàn)出現(xiàn))| (atrr對對ratS的積分只要在球面的積分只要在球面上進(jìn)行上進(jìn)行ratS以以r為球心為球心(矢徑矢徑r),半徑為半徑為atratratSSSdatraSdatrtatrU)(41)(41),(Sd 為球面為球面 的面積元的面積元,此即此即ratSr三維無界空間中的波動三維無界空間中的波動,只要知道初始狀況只要知道初始狀況,就可以用泊松公式就可以用泊松公式ratS然后拿初始

15、擾動然后拿初始擾動)(, )(rr按泊松公式在球面按泊松公式在球面 上積分上積分ratS,波動以速度波動以速度a傳播傳播,只有跟點只有跟點r相距相距at的那些點的初始擾動恰好在時刻的那些點的初始擾動恰好在時刻t傳到傳到rrdDT0初始擾動只限于區(qū)域初始擾動只限于區(qū)域T0,如圖如圖,取一定點取一定點r,與與T0ratS跟跟T0不相交不相交,按泊松公式按泊松公式u(r,t)=0,表示擾動的前鋒表示擾動的前鋒沒有到達(dá)沒有到達(dá)r,當(dāng)當(dāng)d/atD/a, 包圍了包圍了T0,但跟但跟T0不相交不相交,u(r,t)=0,表明表明ratS球心球心,以以at為半徑作球面為半徑作球面求以后任一時刻的狀況求以后任一時刻的狀況,具體說具體說,為求時刻為求時刻t在在r的的u(r,t),應(yīng)以應(yīng)以r為為擾動已經(jīng)過去擾動已經(jīng)過去.最小距離為最小距離為d,最大距離為最大距離為D,當(dāng)當(dāng)td/a,yxat,跟跟T0總有重疊總有重疊,積分一般不為零積分一般不為零在點在點(x,y)總有擾動總有擾動,可以看成某種三維波動的剖面可以看成某種三維波動的剖面.計算計算)(1cr