主題
Search

波動方程——一維

下載 Mathematica 筆記本下載 Wolfram 筆記本

一維波動方程由下式給出

(partial^2psi)/(partialx^2)=1/(v^2)(partial^2psi)/(partialt^2).
(1)

為了指定一個波，該方程受制於邊界條件

psi(0,t)=0
(2)
psi(L,t)=0,
(3)

和初始條件

psi(x,0)=f(x)
(4)
(partialpsi)/(partialt)(x,0)=g(x).
(5)

一維波動方程可以透過達朗貝爾解法、使用傅立葉變換方法或透過分離變數法精確求解。

達朗貝爾在 1746 年設計了他的解法，尤拉隨後在 1748 年擴充套件了該方法。令

xi=x-vt
(6)
eta=x+vt.
(7)

根據鏈式法則，

(partial^2psi)/(partialx^2)=(partial^2psi)/(partialxi^2)+2(partial^2psi)/(partialxipartialeta)+(partial^2psi)/(partialeta^2)
(8)
1/(v^2)(partial^2psi)/(partialt^2)=(partial^2psi)/(partialxi^2)-2(partial^2psi)/(partialxipartialeta)+(partial^2psi)/(partialeta^2).
(9)

則波動方程變為

(partial^2psi)/(partialxipartialeta)=0.
(10)

該方程的任何解都具有以下形式

psi(xi,eta)=f(eta)+g(xi)=f(x+vt)+g(x-vt),
(11)

其中 f 和 g 是任意函式。它們表示兩個沿相反方向傳播的波形，f 沿負x方向傳播，g 沿正x方向傳播。

一維波動方程也可以透過對每一側應用傅立葉變換來求解，

int_(-infty)^infty(partial^2psi(x,t))/(partialx^2)e^(-2piikx)dx=1/(v^2)int_(-infty)^infty(partial^2psi(x,t))/(partialt^2)e^(-2piikx)dk,
(12)

它由傅立葉變換導數恆等式給出，如下所示

(2piik)^2Psi(k,t)=1/(v^2)(partial^2Psi(k,t))/(partialt^2),
(13)

其中

Psi(k,t)=F_x[psi(x,t)](k)
(14)
=int_(-infty)^inftypsi(x,t)e^(-2piikx)dx.
(15)

這有解

Psi(k,t)=A(k)e^(2piikvt)+B(k)e^(-2piikvt).
(16)

取逆傅立葉變換得到

psi(x,t)=int_(-infty)^inftyPsi(k,t)e^(2piikx)dk
(17)
=int_(-infty)^infty[A(k)e^(2piikvt)+B(k)e^(-2piikvt)]e^(-2piikx)dk
(18)
=int_(-infty)^inftyA(k)e^(-2piik(x-vt))dk+int_(-infty)^inftyB(k)e^(-2piik(x+vt))dk
(19)
=f_1(x-vt)+f_2(x+vt),
(20)

其中

f_1(u)=F_k[A(k)](u)=int_(-infty)^inftyA(k)e^(-2piiku)dk
(21)
f_2(u)=F_k[B(k)](u)=int_(-infty)^inftyB(k)e^(-2piiku)dk.
(22)

此解仍然受所有其他初始條件和邊界條件的約束。

一維波動方程可以透過使用試解分離變數法來求解

psi(x,t)=X(x)T(t).
(23)

這給出

T(d^2X)/(dx^2)=1/(v^2)X(d^2T)/(dt^2)
(24)
1/X(d^2X)/(dx^2)=1/(v^2)1/T(d^2T)/(dt^2)=-k^2.
(25)

因此，X 的解是

X(x)=Ccos(kx)+Dsin(kx).
(26)

重寫 (25) 得到

1/T(d^2T)/(dt^2)=-v^2k^2=-omega^2,
(27)

因此，T 的解是

T(t)=Ecos(omegat)+Fsin(omegat),
(28)

其中 v=omega/k。將邊界條件 psi(0,t)=psi(L,t)=0 應用於 (◇) 得到

C=0 kL=mpi,
(29)

其中 m 是一個整數。將 (◇)、(◇) 和 (29) 代回 (◇) 中的 psi 得到，對於 m 的特定值，

psi_m(x,t)=[E_msin(omega_mt)+F_mcos(omega_mt)]D_msin((mpix)/L)
(30)
=[A_mcos(omega_mt)+B_msin(omega_mt)]sin((mpix)/L).
(31)

初始條件 psi^.(x,0)=0 然後給出 B_m=0，因此 (31) 變為

psi_m(x,t)=A_mcos(omega_mt)sin((mpix)/L).
(32)

通解是所有可能的 m 值的總和，因此

psi(x,t)=sum_(m=1)^inftyA_mcos(omega_mt)sin((mpix)/L).
(33)

再次使用正弦函式的正交性，

int_0^Lsin((lpix)/L)sin((mpix)/L)dx=1/2Ldelta_(lm),
(34)

其中 delta_(lm) 是由下式定義的克羅內克 delta

delta_(mn)={1 m=n; 0 m!=n,
(35)

給出

int_0^Lpsi(x,0)sin((mpix)/L)dx=sum_(l=1)^(infty)A_lsin((lpix)/L)sin((mpix)/L)dx
(36)
=sum_(l=1)^(infty)A_l1/2Ldelta_(lm)
(37)
=1/2LA_m,
(38)

所以我們有

A_m=2/Lint_0^Lpsi(x,0)sin((mpix)/L)dx.
(39)

對於特定的初始變形，A_ms 的計算在傅立葉正弦級數部分匯出。我們已經發現 B_m=0，因此弦的運動方程 (◇) 為，

omega_m=vk_m=(vmpi)/L,
(40)

是

psi(x,t)=sum_(m=1)^inftyA_mcos((vmpit)/L)sin((mpix)/L),
(41)

其中 A_m 係數 由 (◇) 給出。

一個阻尼一維波

(partial^2psi)/(partialx^2)=1/(v^2)(partial^2psi)/(partialt^2)+b(partialpsi)/(partialt),
(42)

給定邊界條件

psi(0,t)=0
(43)
psi(L,t)=0,
(44)

初始條件

psi(x,0)=f(x)
(45)
(partialpsi)/(partialt)(x,0)=g(x),
(46)

和附加約束

0<b<(2pi)/(Lv),
(47)

也可以作為傅立葉級數求解。

psi(x,t)=sum_(n=1)^inftysin((npix)/L)e^(-v^2bt/2)[a_nsin(mu_nt)+b_ncos(mu_nt)],
(48)

其中

mu_n=(sqrt(4v^2n^2pi^2-b^2L^2v^4))/(2L)=(vsqrt(4n^2pi^2-b^2L^2v^2))/(2L)
(49)
b_n=2/Lint_0^Lsin((npix)/L)f(x)dx
(50)
a_n=2/(Lmu_n){int_0^Lsin((npix)/L)[g(x)+(v^2b)/2f(x)]dx}.
(51)

另請參閱

達朗貝爾算符, 達朗貝爾解法, Korteweg-de Vries 方程, 拉普拉斯算符, 電報方程, 波動方程, 波動方程——圓盤, 波動方程——矩形, 波動方程——三角形

使用 探索

參考文獻

Abramowitz, M. 和 Stegun, I. A. (Eds.). "拋物線和扁球座標系中的波動方程。" §21.5 in 數學函式手冊，包含公式、圖表和數學表格，第 9 次印刷。 New York: Dover, pp. 752-753, 1972.Morse, P. M. 和 Feshbach, H. 理論物理方法，第一部分。 New York: McGraw-Hill, pp. 124-125 和 271, 1953.Zwillinger, D. (Ed.). CRC 標準數學表格和公式手冊。 Boca Raton, FL: CRC Press, p. 417, 1995.Zwillinger, D. 微分方程手冊，第 3 版。 Boston, MA: Academic Press, p. 130, 1997.

請引用本文為

Weisstein, Eric W. "波動方程——一維。" 來自 —— 資源。 https://mathworld.tw/WaveEquation1-Dimensional.html

主題分類