RORO

ふつうの日記(移転したい)

1次元波動方程式の一般解をラプラス変換とフーリエ変換で求める

免責事項: 私は電磁気学も流体力学も微分方程式も学んだことはありません。デタラメです。

ラプラス変換とフーリエ変換だけでやってみる試みです。デタラメです。


1次元の波動方程式は次の通り:

22xu(x,t)=1c222tu(x,t)\frac{\partial^2}{\partial^2 x} u(x, t) = \frac{1}{c^2} \frac{\partial^2}{\partial^2 t} u(x, t)

初期条件は u(x,0)=g(x)u(x, 0) = g(x), u(x,0)=h(x)u’(x, 0) = h(x) とする。

ttについては因果的なので、tt についてラプラス変換すると:

22xU(x,s)=1c2(s2U(x,s)su(x,0)u(x,0))=1c2(s2U(x,s)sg(x)h(x))\begin{aligned}\frac{\partial^2}{\partial^2 x} U(x, s) &= \frac{1}{c^2} \left( s^2 U(x, s) - s u(x, 0) - u’(x, 0)\right) \\ &= \frac{1}{c^2} \left( s^2 U(x, s) - s g(x) - h(x)\right) \\ \end{aligned}

さらに xx についてフーリエ変換すると:

(jω)2U(ω,s)=1c2(s2U(ω,s)sG(ω)H(ω))\begin{aligned}(-j\omega)^2 U(\omega, s) &= \frac{1}{c^2} \left( s^2 U(\omega, s) - s G(\omega) - H(\omega)\right)\end{aligned}

これを整理すると:

U(ω,s)=12c2(s/c)ω2+(s/c)2G(ω)+12cs2(s/c)ω2+(s/c)2H(ω)\begin{aligned}U(\omega, s) &= \frac{1}{2c} \cdot \frac{2 (s / c)}{\omega^2 + (s / c)^2} G(\omega) + \frac{1}{2cs} \cdot \frac{2 (s / c)}{\omega^2 + (s / c)^2} H(\omega)\end{aligned}

これを逆フーリエ変換すると(\circledastは畳み込みを表すこととする):

U(x,s)=[12ce(s/c)xxg(x)]+[12cse(s/c)xxh(x)]U(x, s) = \left[ \frac{1}{2c} e^{-(s/c)|x|} \circledast_x g(x) \right] + \left[ \frac{1}{2cs} e^{-(s/c)|x|} \circledast_x h(x) \right]

さらに逆ラプラス変換すると:

u(x,t)=12cδ(txc)xg(x)+12c0t(δ(τxc)xh(x))dτ=12δ(ctx)xg(x)+120th(xcτ)+h(x+cτ)dτ=12(g(xct)+g(x+ct))+12cxctx+cth(ξ)dξ\begin{aligned} u(x, t) &= \frac{1}{2c} \delta(t - \frac{|x|}{c}) \circledast_x g(x) + \frac{1}{2c}\int^{t}_{0} \left( \delta(\tau - \frac{|x|}{c}) \circledast_x h(x) \right) d\tau \\ &= \frac{1}{2} \delta(ct - |x|) \circledast_x g(x) + \frac{1}{2} \int^{t}_{0} h(x - c\tau) + h(x + c\tau) d\tau \\ &= \frac{1}{2} \left( g(x - ct) + g(x + ct) \right) + \frac{1}{2c}\int^{x+ct}_{x-ct}h(\xi) d\xi \\ \end{aligned}

Wikipediaに載ってるダランベールの式と同じ解が得られた。


合ってるのかな?