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