根據波動方程式的建模,一個脈衝在一根固定兩端的繩子上的運動
從一個點源發散出的球面波
波動方程式或稱波方程式(英語:wave equation)是一種重要的偏微分方程式,主要描述自然界中的各種的波動現象,包括橫波和縱波,例如聲波、光波和水波。波動方程式抽象自聲學,電磁學,和流體力學等領域。
歷史上許多科學家,如達朗伯特、歐拉、丹尼爾·白努利和拉格朗日等在研究樂器等物體中的弦振動問題時,都對波動方程式理論作出過重要貢獻。
波動方程式是雙曲形偏微分方程式的最典型代表,其最簡形式可表示為:關於位置x 和時間t 的純量函數u(代表各點偏離平衡位置的距離)滿足:
這裡c通常是一個固定常數,代表波的傳播速率。在常壓、20°C的空氣中c為343米/秒(參見音速)。在弦振動問題中,c 依不同弦的密度大小和軸向張力不同可能相差非常大。而在半環螺旋彈簧(一種玩具,英文商標為 Slinky)上,波速可以慢到1米/秒。
在針對實際問題的波動方程式中,一般都將波速表示成可隨波的頻率變化的量,這種處理對應真實物理世界中的色散現象。此時,c 應該用波的相速度代替:
實際問題中對標準波動方程式的另一修正是考慮波速隨振幅的變化,修正後的方程式變成下面的非線性波動方程式:
另需注意的是物體中的波可能是疊加在其他運動(譬如介質的平動,以氣流中傳播的聲波為例)上的。這種情況下,純量u 的表達式將包含一個馬赫因子(對沿流動方向傳播的波為正,對反射波為負)。
三維波動方程式描述了波在均勻各向同性彈性體中的傳播。絕大多數固體都是彈性體,所以波動方程式對地球內部的地震波和用於檢測固體材料中缺陷的超聲波的傳播能給出滿意的描述。在只考慮線性行為時,三維波動方程式的形式比前面更為複雜,它必須同時考慮固體中的縱波和橫波:
式中:
- 和 被稱為彈性體的拉梅常數(也叫「拉梅模量」,英文Lamé constants 或 Lamé moduli),是描述各向同性固體彈性性質的參數;
- 表示密度;
- 是源函數(即外界施加的激振力);
- 表示位移;
注意在上述方程式中,激振力和位移都是向量,所以該方程式也被稱為向量形式的波動方程式。
其他形式的波動方程式還能在量子力學和廣義相對論理論中用到。
[编辑] 純量形式的一維波動方程式
[编辑] 波動方程式的推導
一維波動方程式可用如下的方式推導:一列質量為m的小質點,相鄰質點間用長度h的彈簧連接。彈簧的彈性係數(又稱「倔強係數」)為k:
其中u(x) 表示位於x的質點偏離平衡位置的距離。施加在位於x+h 處的質點m 上的力為:
其中代表根據牛頓第二定律計算的質點慣性力,代表根據虎克定律計算的彈簧作用力。所以根據分析力學中的達朗伯特原理,位於x+h 處質點的運動方程式為:
式中已註明u(x) 是時間t 的顯函數。
若N 個質點間隔均勻地固定在長度L = N h 的彈簧鏈上,總質量M = N m,鏈的總體勁度係數為K = k/N,我們可以將上面的方程式寫為:
取極限 N , h 就得到這個系統的波動方程式:
在這個例子中,波速。
[编辑] 初值問題的解
一維純量形式波動方程式的一般解是由達朗伯特給出的。原方程式可以寫成如下的算子作用形式:
從上面的形式可以看出,若F 和G 為任意函數,那麼它們以下形式的組合
必然滿足原方程式。上面兩項分別對應兩列行波("行"與在"行動"中同音)——F 表示經過該點(x 點)的右行波,G 表示經過該點的左行波。為完全確定F 和G 的最終形式還需考慮如下初始條件:
經帶入運算,就得到了波動方程式著名的達朗伯特行波解,又稱達朗伯特公式:
在古典的意義下,如果並且則。但是,行波函數F和G 也可以是廣義函數,比如狄拉克δ函數。在這種情況下,行波解應被視作左行或右行的一個脈衝。
基本波動方程式是一個線性微分方程式,也就是說同時受到兩列波作用的點的振幅就是兩列波振幅的相加。這意味著可以通過把一列波分解成它的許多分量來分析其行為。傅立葉變換就是將一個波分解成正弦/餘弦分量的方法,因此在波動方程式的求解中很有效。
[编辑] 純量形式的三維波動方程式
三維波動方程式初值問題的解可以通過求解球面波波動方程式得到。求解結果可用於推導二維情況的解。
[编辑] 球面波
球面波方程式的形式不隨空間坐標系統的轉動而變化,所以可以將它寫成僅與距源點距離r 相關的函數。方程式的三維形式為:
將方程式變形為:
此時,因變數ru 滿足一維波動方程式,於是可以利用達朗伯特行波法將解寫成:
其中F 和G 為任意函數,可以理解為以速度c 從中心向外傳播的波和從外面向中心傳播的波。這類從點源傳出的波強度隨距點源距離r 衰減,並且屬於無後效波,可以清晰地搭載信號。這種波僅在奇數維空間中存在(原因將在下一小節中詳細解釋)。幸運的是,我們生活的空間是三維的,所以我們可以清晰地通過聲波和電磁波(都屬於球面波)來互相交流。
[编辑] 時間箭頭的討論
上面方程式的解裡面,分成了兩部分,一部分表示向外傳播的波,一部分則是向內。很明顯,只要將t換成-t,就可以在這兩部分之間轉換。這體現了原始方程式對於時間是對稱的,任意的一個解在時間軸上倒過來看仍然是一個解。
然而,我們所觀察到的實際的波,都是屬於向外傳播的。除非精心地加以調整,我們無法在自然界觀察到向內的波,儘管它們也是波動方程式的合法的解。
關於這個現象,引起了不少討論。有人認爲,實際上它們即使存在,也無法加以觀察。想想如果四周的光向一個物體集中,則因爲沒有光到達我們的眼睛,我們不可能看見這個物體或者發現這個現象(見參考文獻[2] )。
[编辑] 廣義初值問題的解
波動方程式中u 是線性函數,並且不隨時間和空間坐標的平移而改變。所以我們可以通過平移與疊加球面波獲得方程式各種類型的解。令 φ(ξ,η,ζ) 為任意具有三個自變數的函數,球面波形F 為狄拉克δ函數(數學語言是:F 是一個在全空間積分等於1且非零區間收縮至原點的連續函數的弱極限)。設(ξ,η,ζ)位一族球面波的源點,r 為距源點的徑向距離,即:
可定義
稱為三維波動方程式的影響函數,其意義為(ξ,η,ζ)點在t=0 時刻受到短促脈衝δ函數作用後向空間中傳出的波的影響,係數分母 4πc 是為方便後續處理而加上的。
若u 是這一族波函數的加權疊加,且權函數為 φ,則
從δ函數的定義可知,u 還能寫成
式中α、β 和 γ 是單位球面S 上點的坐標,dω 為S 上的面積微元。該結果的意義為:u(t,x,y,z) 是以(x,y,z) 為圓心,ct 為半徑的球面上φ 的平均值的t 倍:
從上式易得
平均值是關於t 的偶函數,所以若
那麼
以上得出的便是波動方程式初值問題的解。從中可以看出,任意點P 在t 時刻受到的波擾動只來自以P 為圓心,ct 為半徑的球面上,而這個球的內部點在這一時刻對P 點的狀態完全沒有影響(因為它們的影響之前就已經傳過P 點了)。換一個角度分析,假設三維空間中任意點P' 在t=0 時刻受到一個脈衝擾動δ,那麼由此發出的球面波在傳過空間中的任意其它點Q 後,便再也不會對Q 的運動狀態產生影響,這就是在物理學中也非常著名的惠更斯原理(Huygens' principle),也稱為無後效現象,表示傳過的球面波不會留下任何後續效應。
下面我們便可以解釋上一小節中留下的問題了。事實上,前面所得到的球面波解僅在奇數維空間中存在。偶數維空間中波動方程式的解是彌散的,也就是說波陣面掠過區域仍然會受其影響。以下面的二維波動方程式(極坐標形式,注意和上一小節三維形式的差別)為例:
可以從三維形式的解通過降維法得到二維波動方程式的影響函數:
其中
設點M(x,y) 到點(ξ,η) 距離為d,那麼從影響函數中可以看出,當t >d /c 即初始擾動已傳過M 點後,M 仍在受到它的影響。二維球面波(柱面波)的這一性質決定了它不能作為傳遞信號的工具,因為這種波(事實上包括所有偶數維空間中的球面波)經過的點受到的是交織在一起的各個不同時刻的擾動。
[编辑] 純量形式的二維波動方程式
二維波動方程式的直角坐標形式為:
如前所述,我們可以從三維波動方程式的解中將u 視為與其中一個自變數無關(降維法)來得到二維形式的解。將初始條件改寫為
則三維形式的解就變成
其中 α 和 β 是單位球面上點的頭兩個坐標分量,dω 是球面上的面積微元。此積分可變換為在(x,y) 為中心, ct 為半徑的圓域D 上的積分:
從這個結果也能得到上一小節最後的結論。
二維波動方程式解的一個例子是緊繃的鼓面的運動。
[编辑] 邊值問題
[编辑] 一維情形
一根自身繃緊,兩端分別固定於x=0 和x=L 的彈性弦在t>0 時刻,0 < x < L 上運動滿足波動方程式。在邊界點處,可以要求u 滿足各種邊界條件。通常遇到的邊界條件都可歸納成下列形式:
其中a、b 非負。若要弦的兩端固定不動,對應上面式子中a、b 趨於無窮大。求解偏微分方程式的分離變數法要求尋找以下形式的解:
將上述假設形式代入原方程式中可以得到:
為使邊值問題有非平凡解,本徵值 λ 須滿足
這是固有值問題的斯圖姆-劉維爾理論的一個特例。若a、b 為正數,則對應的所有本徵值均為正數,方程式的解為三角函數。使u 和ut 滿足平方可積條件的解可以通過適當選取u 和ut 三角級數展開來求得。
[编辑] 多維情形
一維初始值-邊值理論可以拓展至任意維空間中。考慮m 維空間(坐標簡寫為x)中的域D ,B為D的邊界。 當0<t 時,位於D 內的點x 滿足波動方程式。在D 的邊界上,解u 須滿足
其中n 是B 上指向域外的法向向量,a 是定義在B 上的非負函數。要求u 在B 上始終為0的邊界條件相當於令a 趨於無窮。初始條件為
其中f 和g 是定義在D 內的函數。這個問題可以通過將f 和g 展開成域D 內拉普拉斯算子滿足邊界條件的本徵函數系的疊加來求解(這是分離變數法的一般步驟)。也就是求解在域D 內滿足
在邊界B 上滿足
的本徵函數系v 。
在二維情形下,上述本徵函數系可以理解成繃緊地張在邊界B 上的鼓面的自由振動模態。若B 是一個圓,則這些本徵函數是關於極角自變數θ 的三角函數與關於極軸自變數r 的整階貝索函數的乘積。更詳細的說明參見英文版條目亥姆霍茲方程式。
在三維形式下,若邊界是空間中的球面,那麼本徵函數是關於球坐標下兩個極角自變數的球面調和函數,乘以關於徑向自變數ρ 的半奇數階貝索函數。
[编辑] 註釋
[编辑] 參考文獻
- [1] 嚴鎮軍編,《數學物理方程式》,第二版,中國科學技術大學出版社,合肥,2002,第210頁~第224頁,ISBN 7-312-00799-6/O·177
- [2] [英]胡·普賴斯著,肖巍譯,《時間之矢與阿基米德之點—物理學時間的新方向》,上海科學技術出版社,上海,2001, ISBN 7-5323-5737-6
- [3] M. F. Atiyah, R. Bott, L. Garding, Lacunas for hyperbolic differential operators with constant coefficients I, Acta Math., 124 (1970), 109–189.
- [4] M.F. Atiyah, R. Bott, and L. Garding, Lacunas for hyperbolic differential operators with constant coefficients II, Acta Math., 131 (1973), 145–206.
- [5] R. Courant, D. Hilbert, Methods of Mathematical Physics, vol II. Interscience (Wiley) New York, 1962.
[编辑] 參看
[编辑] 外部連結