SST Lab Dokuwiki Header header picture

ユーザ用ツール

サイト用ツール


amath2:楔型_三角波

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


楔形関数のフーリエ変換/三角波のフーリエ展開

$$ f(x) = \int_{-\infty}^{\infty}F(u)e^{iux}du $$ とすると$F(u)$は $$ F(u)=\frac{1}{2\pi}\int_{-\infty}^{\infty}f(t)e^{-iut}dt $$ これに$-a<x<a$で $$ f(x) = b\left(1-\frac{|x|}{a}\right) $$ を代入すると $$ F(u)=\frac{b}{2\pi}\int_{-a}^{a}\left(1-\frac{|t|}{a}\right)e^{-iut}dt $$ ただし$a>0$, $b>0$とする。

$f(x)$が偶関数であることを用いると $$ F(u)=\frac{b}{\pi}\int_{0}^{a}\left(1-\frac{t}{a}\right)\cos(ut)dt$$

ここで $$ \frac{d}{dt}\left\{\left(1-\frac{t}{a}\right)\sin(ut)\right\}=-\frac{1}{a}\sin(ut)+\left(1-\frac{t}{a}\right)u\cos(ut)$$

この両辺を$t$で不定積分すると、 $$ \left(1-\frac{t}{a}\right)\sin(ut)=\frac{1}{au}\cos(ut)+u\int\left(1-\frac{t}{a}\right)\cos(ut)dt$$ ただし積分定数は省略した。これより$[0,a]$の定積分は $$ \left[\left(1-\frac{t}{a}\right)\sin(ut)\right]_{t=0}^{t=a}=0=\frac{1}{au}(\cos(au)-1)+u\int_0^a\left(1-\frac{t}{a}\right)\cos(ut)dt$$ したがって

したがって $$ F(u)=\frac{b}{\pi}\int_{0}^{a}\left(1-\frac{t}{a}\right)\cos(ut)dt=\frac{b}{\pi}\frac{1}{au^2}(1-\cos(au))$$

これはF(u)に関する常微分方程式であるので、これを解くと $$\frac{1}{F(u)}\frac{d}{du}F(u)=-\frac{1}{2a}u$$ $$\int\frac{1}{F(u)}\frac{d}{du}F(u)du=-\frac{1}{2a}\int udu$$ $$\int\frac{1}{F(u)}d(F(u))=-\frac{1}{2a}\frac{u^2}{2}+C$$ ここで$C$は積分常数である。 $$\ln|F(u)|=-\frac{u^2}{4a}+C$$ $$F(u)=Ce^{-\frac{u^2}{4a}}$$ ただし、積分常数は出るたびに記号$C$にまとめなおしている。この$C$は$u=0$と置くことにより $$F(0)=C$$ したがって$F(u)$は $$F(u)=F(0)e^{-\frac{u^2}{4a}}$$

一方$F(0)$は $$ F(0)=\frac{1}{2\pi}\int_{-\infty}^{\infty}e^{-at^2}dt $$ であるから $$ I=\int_{-\infty}^{\infty}e^{-at^2}dt $$ とおくと $$ F(0)=\frac{1}{2\pi}I $$

積分変数はなんでもいいので $$ I=\int_{-\infty}^{\infty}e^{-ax^2}dx $$ $$ I=\int_{-\infty}^{\infty}e^{-ay^2}dy $$ これより $$ I^2=\int_{-\infty}^{\infty}e^{-ax^2}dx\int_{-\infty}^{\infty}e^{-ay^2}dy $$ これは$x,y$の逐次積分とも、二重積分ともとることができるので $$ I^2=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-a(x^2+y^2)}dxdy $$

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

$$ I^2=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-a(x^2+y^2)}dxdy=\int_{0}^{2\pi}\int_{0}^{\infty}e^{-ar^2}rdrd\theta $$ これは計算できて、 $$ I^2=\int_{0}^{2\pi}\int_{0}^{\infty}e^{-ar^2}rdrd\theta =2\pi\int_{0}^{\infty}e^{-ar^2}rdr=\frac{\pi}{a}\int_{0}^{\infty}e^{-ar^2}(2ar)dr $$

$$ \int_{0}^{\infty}e^{-ar^2}(2ar)dr=\int_{0}^{\infty}e^{-ar^2}\frac{d(ar^2)}{dr}dr=\int_{0}^{\infty}e^{-ar^2}d(ar^2)=[-e^{-ar^2}]_0^{\infty}=-(0-1)=1 $$

よって $$ I^2=\frac{\pi}{a} $$ $$ I=\sqrt{\frac{\pi}{a}} $$ $$ F(0)=\frac{1}{2\pi}I=\frac{1}{2\pi}\sqrt{\frac{\pi}{a}}=\frac{1}{\sqrt{4a\pi}} $$

まとめると $$F(u)=\frac{1}{\sqrt{4a\pi}}e^{-\frac{u^2}{4a}}$$

したがって $$ f(x) = \int_{-\infty}^{\infty}F(u)e^{iux}du $$ を具体的に書き下すと $$ e^{-ax^2} =\frac{1}{\sqrt{4a\pi}}\int_{-\infty}^{\infty}e^{-\frac{u^2}{4a}}e^{iux}du $$

amath2/楔型_三角波.1595912464.txt.gz · 最終更新: 2022/08/23 13:34 (外部編集)

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki