Processing math: 100%

応用数学 第11回 (4) 熱伝導方程式再び

問題設定

 前回扱った熱伝導方程式 を、今日は 境界条件を外して 解きます。 変数 t はフーリエ変換の変数に使いますので、今日は時刻を s で表します。
 u=u(x,s) についての熱伝導方程式
c22ux2=us
を初期条件
u(x,0)=f(x)
のもとで解け。( ただし f(x) は与えられた関数で、時刻 s=0 での熱分布を表す。)
 境界条件をはずしたということは、棒が無限の長さを持っているということです。

Step 1 ( フーリエ変換 )

 u の、x についてのフーリエ変換ˆu と表すと、 Th.4 (2) より
^(2ux2)=(it)2ˆu=t2ˆu
ˆus=s(12πu(x,s)eitxdx)=12π(su(x,s))eitxdx=12π(c22ux2)eitxdx=c2^(2ux2)=c2t2ˆu s の関数として解くと
ˆu=Cec2t2s
初期条件より
ˆu(t,0)=^u(x,0)=ˆf(t)
ゆえ
C=Cec2t2×0=ˆu(t,0)=ˆf(t)
よって
ˆu(t,s)=ˆf(t)ec2st2

Step 2 ( フーリエ逆変換 )

 この式をフーリエ逆変換すると、Th.7 より
u(x,s)=12πf(x)(ec2st2 のフーリエ逆変換 )
Ex.8 を使うために 14α=c2s とおくと α=14c2s
(ec2st2 のフーリエ逆変換 )=2αeαx2=1c2sex2/(4c2s)
よって
u(x,s)=12cπsf(xy)ey2/(4c2s)dy.
積分表示の形で解が得られました。


 次回から勉強するラプラス変換は、 この「変換して、簡単に解いて、逆変換」というパターンで微分方程式を解くための 強力な道具です。