応用数学 第11回 (4) 熱伝導方程式再び
問題設定
前回扱った熱伝導方程式 を、今日は
境界条件を外して 解きます。
変数
t はフーリエ変換の変数に使いますので、今日は時刻を
s で表します。
問 u=u(x,s) についての熱伝導方程式
c2∂2u∂x2=∂u∂s
を初期条件
u(x,0)=f(x)
のもとで解け。( ただし f(x) は与えられた関数で、時刻 s=0 での熱分布を表す。)
境界条件をはずしたということは、棒が無限の長さを持っているということです。
Step 1 ( フーリエ変換 )
u の、
x についてのフーリエ変換を
ˆu と表すと、
Th.4 (2) より
^(∂2u∂x2)=(it)2ˆu=−t2ˆu
∴ ∂ˆu∂s=∂∂s(1√2π∫∞−∞u(x,s)e−itxdx)=1√2π∫∞−∞(∂∂su(x,s))e−itxdx=1√2π∫∞−∞(c2∂2u∂x2)e−itxdx=c2^(∂2u∂x2)=−c2t2ˆu
s の関数として解くと
ˆu=Ce−c2t2s
初期条件より
ˆu(t,0)=^u(x,0)=ˆf(t)
ゆえ
C=Ce−c2t2×0=ˆu(t,0)=ˆf(t)
よって
ˆu(t,s)=ˆf(t)e−c2st2
Step 2 ( フーリエ逆変換 )
この式をフーリエ逆変換すると、
Th.7 より
u(x,s)=1√2πf(x)∗(e−c2st2 のフーリエ逆変換 )
Ex.8 を使うために
14α=c2s とおくと
α=14c2s で
(e−c2st2 のフーリエ逆変換
)=√2αe−αx2=1c√2se−x2/(4c2s)
よって
u(x,s)=12c√πs∫∞−∞f(x−y)e−y2/(4c2s)dy.
積分表示の形で解が得られました。
次回から勉強するラプラス変換は、
この「変換して、簡単に解いて、逆変換」というパターンで微分方程式を解くための
強力な道具です。