Processing math: 100%
SST Lab Dokuwiki Header header picture

ユーザ用ツール

サイト用ツール


amath2:ガウス関数

文書の過去の版を表示しています。


ガウス関数のフーリエ変換

f(x)=F(u)eiuxdu とするとF(u)F(u)=12πf(t)eiutdt これに f(x)=eax2 を代入すると F(u)=12πeat2eiutdt ただしa>0とする。

eiz=cosz+isinzを用いると F(u)=12πeat2cos(ut)dt+i12πeat2sin(ut)dt

第二項の被積分関数eat2sin(ut)は奇関数であるから第二項の積分は零になる。したがって F(u)=12πeat2cos(ut)dt

この両辺をuで微分すると、 dduF(u)=u{12πeat2cos(ut)dt}=12πteat2sin(ut)dt

ここでeat2sin(ut)tで微分してみると ddt[eat2sin(ut)]=2ateat2sin(ut)+ueat2cos(ut) この両辺を[,]で定積分すると、 [eat2sin(ut)]=2ateat2sin(ut)dt+ueat2cos(ut)dt

左辺は零になるので teat2sin(ut)dt=u2aeat2cos(ut)dt

したがって dduF(u)=12πu2aeat2cos(ut)dt=u2aF(u)

これはF(u)に関する常微分方程式であるので、これを解くと 1F(u)dduF(u)=12au 1F(u)dduF(u)du=12audu 1F(u)d(F(u))=12au22+C ここでCは積分常数である。 ln|F(u)|=u24a+C F(u)=Ceu24a ただし、積分常数は出るたびに記号Cにまとめなおしている。このCu=0と置くことにより F(0)=C したがってF(u)F(u)=F(0)eu24a

一方F(0)F(0)=12πeat2dt であるから I=eat2dt とおくと F(0)=12πI

積分変数はなんでもいいので I=eax2dx I=eay2dy これより I2=eax2dxeay2dy これはx,yの逐次積分とも、二重積分ともとることができるので I2=ea(x2+y2)dxdy

(x,y)を二次元空間の直交座標とみなしてこれを極座標に変換すると

I2=ea(x2+y2)dxdy=2π00ear2rdrdθ これは計算できて、 I2=2π00ear2rdrdθ=2π0ear2rdr=πa0ear2(2ar)dr

0ear2(2ar)dr=0ear2d(ar2)drdr=0ear2d(ar2)=[ear2]0=(01)=1

よって I2=πa I=πa F(0)=12πI=12ππa=14aπ

amath2/ガウス関数.1578453388.txt.gz · 最終更新: 2022/08/23 13:34 (外部編集)

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki