矩形 波 フーリエ 変換。 3. フーリエ変換 (やる夫で学ぶディジタル信号処理)

MTFの測定-チャート法-

やらない夫 あそこで をくくり出さなかったどうなるか,計算してみればわかるんだが,フーリエ変換の式の方の先頭に がついて,逆変換の式の方には何もつかなくなるんだ.だから,「フーリエ変換」の式をきれいに見せたかったら今回みたいに でくくればいいし,「フーリエ逆変換」の式をきれいに見せたいならば,くくらずに導出したもので定義すればよかった.どっちにしろ,フーリエ変換して,またフーリエ逆変換すればちゃんと元に戻るからな.どっちでもよいんだけど,我々は前者を採用したってことだ.教科書によっては両方に をつけているのもあるしな. やる夫 どっちでもいいってのはあまり納得いかないお.定義が変わったら周波数成分の値が定数倍だけ変わってしまうお. やらない夫 変わってもいいんだよ.例えばもとの時間信号の振幅が,そうだな,電圧だったとしようか.じゃあそれをフーリエ変換したときの周波数成分の単位はどうなる? ダメだお. で にしたときの値が定まらないお. やらない夫 ああ,それがさっき言った根本的な問題だ.よくよく考えてみると,デルタ関数なんて普通の「関数」ではなかったわけだろ.普通の意味の積分の結果としてさくっと出てくるようなものではないんだ.特殊な取扱いが必要になる. sinc 関数をフーリエ逆変換するときのように「計算が難しい」のとはちょっと話が一味違う. やる夫 じゃあどうすればいいのかお? 連続時間 で定義された関数 のうち実用上重要なものの多く に対して,式 で計算される を のフーリエ変換と呼ぶ. あるいはこの計算をすること自体をフーリエ変換と呼ぶ•。 4 フーリエ級数とフーリエ変換の関係 やる夫 うーん やらない夫 どうした? やらない夫 あきらめないためには超関数論にまともに踏み込まないといけないからな.そこで「フーリエ変換対をセットにする」考え方の出番だ. をフーリエ変換したら になると知っているんだから, のフーリエ逆変換が になることもわかっていると考えてしまう. やる夫 さっきも言ったけど,ずるいお. やらない夫 さっきも言ったが,ずるくて構わない.ともかく,複素指数関数をフーリエ変換しなくちゃならない状況になったときに,あー,これはデルタ関数になるなと思い出して,逆から計算できるようになれれば勝ちだ.差し詰め,試合に負けて勝負に勝つといったところか. やる夫 あ,やっぱり試合は終了なのかお. やらない夫 同様に,周波数領域のデルタ関数 とか,それを定数 だけシフトした とかのフーリエ逆変換も計算しておこうか.やり方はほとんど同じだ. やる夫 同じように計算すると 3. やる夫 そりゃ,やっぱり時刻 ずつシフトしてやればいいはずだお.周波数 のサイン波は 1 周期が だから,位相で考えると,比例計算で 3. DFTの練習問題としてよく出てくるのが, 1 方形波 2 正弦波 の2つです. 正直フーリ変換と違ってDFTは手計算できるものがかなり限られてくるので この2つを押さえておくだけでもだいぶ違うかなと思います. 今回は 1 方形波を扱います. 2 正弦波は次回扱おうと思います. DFTを実際に手計算で計算し,それをPythonでFFTしたものと一致するかを確かめてみたいと思います. 動作環境• , , をそれぞれ, の振幅スペクトル,位相スペクトル,パワースペクトルと呼ぶ. やる夫 あれ? やる夫 えっと,基本角周波数が で,その整数倍の周波数成分だけがでてくるんだったお. やらない夫 そうだな.だからスペクトルは飛び飛びに値を持つことになる.図でかくとこんな感じだったな. やる夫 なんか強引な気がするお.そんなんでいいんかお. やらない夫 やや乱暴かな.まあ気にするな.ともかく,周期を長くしていったときに,周波数領域がどういう風に変化していくかを考えていこう.で,出発地点に戻ると,周期 のときは,周波数領域では おきに飛び飛びに値を持つんだったわけだろ. やる夫 そうだお.さっきのグラフの通りだお. やらない夫 周期が になったらどうなる? やる夫 うーん, は反比例のグラフだお.反比例と sin をかけたグラフだから, が正のときは,sin なんだけど振幅が に反比例して減っていくようなグラフになるお. が負のときは…反比例の部分が負だから,sin 関数の正負がひっくり返ったものになって,その振幅はやっぱり の絶対値に反比例して減っていくわけだお.だから左右対称なグラフになりそうだお.よくわからないのは の近辺だお.反比例は無限大に,sin はゼロに近づいていくから,かけ合わせた結果どうなるのか,すぐにはわからんお. やらない夫 のときの値が なのは計算の結果わかっていただろう.で,実はちゃんと連続につながったグラフになるんだ. の場合をプロットしてみるとこうなる. なので,単に sinc 関数と言われた場合は,実際にはどっちを指しているかちょっと注意が必要だ. やる夫 面倒くさいお. やらない夫 まあとにかく,定数倍はさておくとして,矩形関数と sinc 関数がフーリエ変換対の関係になっていることを,しっかり把握しておいてくれ. やる夫 ということは,sinc 関数に対してフーリエ変換の計算をすれば矩形関数が出てくるのかお? 何か納得行かない顔をしてるな. やる夫 今日の話って元々,周期信号しか扱えなかったフーリエ級数展開を,周期的じゃない信号に適用できるようにするって流れだったお.それで出てきたのがフーリエ変換なわけだお. やらない夫 そうだったな. やる夫 でも,さっき計算した とか とかも周期的な信号なんだお.周期的な信号もフーリエ変換できるなら,フーリエ級数はもう要らないのかお? やらない夫 そうだな.「重ね合わせ」という言葉であれば,総和の場合も積分の場合も,まあそんなに違和感無く表現できてる気がするが,どうだろう.まあ語感は人それぞれかも知れないけどな. ともかく,一般の時間信号は,あらゆる実数を周波数とする複素指数関数の重ね合わせで表すことができる,ということだ.これがフーリエ逆変換の意味だ. やる夫 逆? やらない夫 出てくる.出てくるんだが,その計算は割とややこしい.何しろ sinc 関数は不定積分が初等関数の組み合わせで書けないんだ.なので積分の計算にいろいろと技巧が必要だ.というわけで「sinc と矩形はフーリエ変換対」と覚えておいて,例えば時間領域の sinc 関数のフーリエ変換が必要になったときには,周波数領域の矩形関数から考えて逆算するようにする方が楽ちんだ. やる夫 ふーん,なんだかずるいお. やらない夫 いいんだよ,ずるくても.この「フーリエ変換対をセットにして考える」という戦略はとても重要だ.sinc 関数の場合は単に計算が面倒なのを回避するだけだが,もっと根本的な問題を回避する場合にも有用なんだ. やる夫 もっと根本的って,どんな場合だお? ああ,周期信号だから,1周期分のコピーが延々に続くわけだお.だから無限に足し合わせることになるんだお. やらない夫 そう,だから の整数倍のところでは,フーリエ変換の値は無限大に発散する.それが,周期関数のフーリエ変換がデルタ関数の並んだものになる理由だ.三角関数や複素指数関数をフーリエ変換したときにデルタ関数が出てくるのも,その特殊な場合になっているだけだな. やる夫 なるほど,辻褄は合ってるお. やらない夫 デルタ関数が「普通の意味での関数」ではなかったことを思い出してくれ.言い換えると,周期信号は,普通の意味ではフーリエ変換が存在しないってことだ.それでは不便なので,関数の意味を拡張して考えている.本来はもっと厳密な扱いが必要だ.もし避けることができるなら避けておきたい. やる夫 正直,どっちみち厳密性にはあまり興味ないから,今のままのゆるーい理解でいいなら,別に避けなくてもいいんじゃないかお? 積分っていっても飛び飛びにしか値がないんだお.こんなもの積分しても何も出てこないお. やらない夫 ああ.まず,高校で習ったような,等幅の短冊の面積の総和で近似して極限を取るような考え方では,こういう飛び飛びにしか値がないような関数の積分は定義できない.だからもう少しうまく定義された積分を導入する必要があるんだが,いずれにせよ,こういう「面積のない」関数の積分は 0 になる. やる夫 なんかよくわからないけど,面積がないから積分が 0 って話は抵抗なく受け入れられるお. やらない夫 フーリエ逆変換したら全部 0 になっちゃうようではお話にならないわけだ.だから,各 に比例したデルタ関数を考えてやることにして「面積を持つ」ようにしてやるわけだ.そうすればフーリエ逆変換の公式で,うまく時間信号に戻るようになる. やる夫 うーん,雰囲気はわからないでもないけど,いまいち理解しにくいお. やらない夫 雰囲気だけでもつかんでおくといい.というのは,この考え方は後々いろんなところで出てくるんだ. 座標軸上の飛び飛びの点でしか値をもたないような,つまり離散的な信号に対して,連続信号用の処理を適用したいときに,デルタ関数をかけておくことで各点の値に「面積」を持たせるという考え方だ.離散信号と連続信号が同じ密度を持つようにするためのトリックだと思ってもいい. やる夫 ピンと来ないお. やらない夫 実は同じような考え方は今回の最初の方で既に使っているんだ.フーリエ級数からフーリエ変換に移行するときに,だろ.短冊の面積に置き換える代わりに,「デルタ関数の面積」に置き換えているのが今の話だ. やる夫 うーん,わかったような,わからないような感じだお. やらない夫 まあ,またすぐに別の具体例が出てくるので,そのとき話そう.. やらない夫 ほら,フーリエ変換の公式を導くときにだろう.その分の辻褄を合わせるために出てきたものだ. やる夫 あー,そういえばそうだったお. やらない夫 今まで見てきたような話を,フーリエ逆変換の視点から見ておくことも重要だ.各周波数成分が有限値 で,それ以外が 0 になっているようなスペクトル が与えられたとしよう.そのままフーリエ逆変換の公式に入れるとどうなる? やらない夫 数学的な厳密性に興味がない場合でも,例えば実際にコンピュータで計算しようと思ったら,無限大を扱うのは厄介だろ. やる夫 ああ,それはそうかも知れないお.紙の上に描くなら矢印にすればいいだけだけど,コンピュータではそうもいかないお. やらない夫 そういう面倒さを避けるために,1 周期だけ積分することにしたのがフーリエ係数だ,と考えてもいいかな.どうせ周期的な信号なんだから,1 周期分だけ考えれば各周波数成分がどういう「割合」で含まれているかを知るには十分だ. やる夫 結局,周期信号をフーリエ変換すると,フーリエ級数展開したときの に比例した高さのデルタ関数が並ぶことになるってことでいいのかお? 離散フーリエ変換 DFT の仕組みを完全に理解する 2 前提知識 高校2年生程度でも分かるぐらいを目指して書いています.しかし高度な高校数学は使わないので,数学の断片的な知識があれば中学生でも分かると思います. 以下,知っていると嬉しいことをリストします.• 連続時間 で定義された関数 のうち実用上重要なものの多く に対して,式 で計算される を のフーリエ変換と呼ぶ. あるいはこの計算をすること自体をフーリエ変換と呼ぶ• あ,そうか,フーリエ級数展開に対応するのは,フーリエ変換じゃなくてフーリエ逆変換の方なんだお.なんか混乱しそうだお. やらない夫 そう,フーリエ変換は,フーリエ係数の計算の方に対応している.それぞれの周波数成分がどのくらい含まれているかを知るための計算になっているということだな. は一般に複素数になるから,振幅と位相を持っている.フーリエ係数 と同様に,周波数 の成分の振幅と初期位相を表しているわけだ. やる夫 フーリエ級数展開やフーリエ係数の計算を「変換」と呼んじゃダメなのかお? 38 やる夫 積分したら 1 になるわけだお.これがどう重要なのかお? 1のフーリエ変換 ってのが今回の話になる.結論からいうと,それがフーリエ変換だ. やる夫 「級数」が「変換」に変わるんかお.なんか「周期的」かどうかとは全く異質な話に聞こえるお. やらない夫 そうかもな.まあその辺は追々理解してもらえばいい.ともかく出発地点はフーリエ級数だ.周期 の時間信号を周波数成分に分解するんだった.どんな周波数成分が出てくる? やる夫 うーん,スペクトルの線の間隔がどんどん狭くなっていくお.だから,飛び飛びじゃないスペクトルになるのかお. やらない夫 そういうことだ. から の連続時間上で定義された時間関数は,周波数領域で見ると, から の連続周波数上で定義されたスペクトルになる.ちょっと議論は乱暴だったけど,ああ何かそうなりそうだな,と納得してもらえればとりあえず OK としよう. やる夫 ふーん,まあ言ってることの雰囲気はわかるお. やらない夫 さて,実際にそういう極限を考えたときに,数式としてはどんな形になるのかっていうのが次の話だ.ところがちょっと問題があって,今の話の流れで考えていても,実は答えにはたどり着けないんだ. やる夫 ちゃぶ台返しかお.じゃあ今までの話はなんだったんだお. やらない夫 まあそう言うな.飛び飛びの離散周波数から連続周波数になっていくイメージを持ってもらいたかっただけだ.でも,どんなに間隔が細かくなっても線は線のままだからな.そのままじゃ連続にはならない.なのでそこはちょっと連続化のための手続きを踏んでやる必要がある. やる夫 どういうことかお. やらない夫 フーリエ級数展開の式から出発しよう.前回の式 ,つまりこれだ. やらない夫 そういうことだ.これで,この短冊の面積をすべて足し合わせると になるようにできたわけだ.こうやって「総和を計算する問題」を「面積を計算する問題」に書き換えておいてから,分割をどんどん細かくしていけば,「面積を積分で求める問題」に持って行くことができる. やる夫 うーん,なんか微妙にしっくり来ないけど,そんなもんなのかお. やらない夫 同じ無限でも,「整数が無限にある」というときの無限と「実数が無限にある」というときの無限との間には大きなギャップがあるんだ.だから「線」のまま間隔を狭くしていっても連続にはならない.そのギャップを,面積を持つ短冊を考えることで埋めていると思ってくれ. 今の話を数式で書くとこうなる.まずフーリエ級数の式を,面積の総和だと思って書き換える. そしてさっきの式 の方をフーリエ逆変換と呼ぶ. やる夫 いつの間にか「級数展開」が「変換」になったお. やらない夫 いつの間にかというか,いつ「変換」になったかと敢えて答えるなら,無限に飛ばして連続化したときだな.その時点で「連続時間上の関数」と「連続周波数上の関数」の相互間の「変換」になったと考えている. フーリエ変換の計算式の右辺には時間変数 と周波数変数 が含まれているが, で積分するから, だけが残る.連続時間上の関数から連続周波数上の関数への変換になるわけだ.フーリエ逆変換の方は,右辺を で積分しているから, だけが残るんだな.時間関数への変換になる. やる夫 結局,周波数が連続になっただけで,フーリエ級数と同じようなものだと思っていいのかお? やる夫 どうなるって言われても,まあ公式につっこんでみるお. やらない夫 の解釈としてはその通りだな.時間領域で見てきた複素指数関数を,時間と周波数を入れ替えて考えればいい.さて, のスペクトルがこういう螺旋になるというのは,どう捉えればいいだろう? やらない夫 そうだな,基本的な考え方は同じだ.フーリエ級数は,周期的な時間信号を無限個の複素指数関数の足し合わせで表現したわけだ.ただし無限といっても高々「整数の個数」の無限だ.周波数成分は飛び飛びにしか存在しないが,それで元の時間関数が十分に再現できた. これに対して,周期的とは限らない一般の時間信号を表現しようと思うと.周波数としてはあらゆる実数を考えなくてはならなくなる.数式で表現すると複素指数関数の「総和」ではなくて「積分」で表現しなくてはならないわけだ. やる夫 フーリエ級数の「複素指数関数の足し合わせで表す」っていう考え方は直観的にわかりやすかったお.でも総和じゃなくて積分になるとどうもピンと来ないお. やらない夫 そうかもしれないが,本質的には全く同じことなんだ.同じイメージを持っていて構わない.ただし「足し合わせ」という言葉を使うのはさすがに違和感があるので,「重ね合わせ」という言葉を使うことが多い. やる夫 「重ね合わせの原理」とかいう場合の重ね合わせと同じかお? , , をそれぞれ, の振幅スペクトル,位相スペクトル,パワースペクトルと呼ぶ. やる夫 あれ? やる夫 うーん, は反比例のグラフだお.反比例と sin をかけたグラフだから, が正のときは,sin なんだけど振幅が に反比例して減っていくようなグラフになるお. が負のときは…反比例の部分が負だから,sin 関数の正負がひっくり返ったものになって,その振幅はやっぱり の絶対値に反比例して減っていくわけだお.だから左右対称なグラフになりそうだお.よくわからないのは の近辺だお.反比例は無限大に,sin はゼロに近づいていくから,かけ合わせた結果どうなるのか,すぐにはわからんお. やらない夫 のときの値が なのは計算の結果わかっていただろう.で,実はちゃんと連続につながったグラフになるんだ. の場合をプロットしてみるとこうなる. なので,単に sinc 関数と言われた場合は,実際にはどっちを指しているかちょっと注意が必要だ. やる夫 面倒くさいお. やらない夫 まあとにかく,定数倍はさておくとして,矩形関数と sinc 関数がフーリエ変換対の関係になっていることを,しっかり把握しておいてくれ. やる夫 ということは,sinc 関数に対してフーリエ変換の計算をすれば矩形関数が出てくるのかお? は角周波数を表す連続変数である. は に含まれる角周波数 の振動成分の量 振幅・位相 を表す.• やらない夫 まず注意してほしいのは,この式では の全区間を積分しているが, が 0 でない値を持つのは を跨ぐ瞬間のみだということだ.だから, を含むような積分範囲を取れば,積分値は必ず 1 になる. それが何を意味するかというと,デルタ関数は,高さは無限大だけど面積は有限で 1 だということだ. やる夫 高さが無限大で,幅が0で,かけたら1になるような短冊だってことかお. やらない夫 そう考える手もあるかな.ともかく,単に「無限大です」ってんじゃなく,何らかの意味で「大きさを考えられる」という点が重要だ.単に無限大だと言われた場合は,その2倍とか3倍とかを考えることに意味がない.でもデルタ関数の場合は, とか とかがちゃんと意味を持っている. やる夫 高さはどれも無限大だけど,面積はそれぞれ 2 と 3 だってことかお. やらない夫 そういうことだ.この性質は,他の関数とデルタ関数をかけ合わせて積分するときに重要だ.有限の値を持つ関数 とかけ合わせて積分すると やる夫 のときの瞬時値が取り出されるわけだお. やらない夫 そう.どの場合も,積分範囲はインパルスの立っているところを含んでさえいれば OK だってことに注意しておこう. やる夫 しかし,奇妙な関数というか不思議な関数だお.こんなもの実在するのかお? やる夫 なんか,ある一点でだけ無限大の値を持って,それ以外の点では 0 になるようなやつだったお.ピーンとインパルスが立ってる感じだお. やらない夫 ああ,イメージとしてはそれで OK だ.無理やり数式で書くとすると 3. 〒673-0036 兵庫県明石市松江62-14 有 ゴッドフット企画 感想をお聞かせ下さい 脱線事故シミュレーションもEXCELで (Yahoo Japan 掲載) Copyright C GODFOOT 2007. やらない夫 そういうことになるな. やる夫 どうして式 は, じゃなくて が並んだものになるのかお? フーリエ変換の定義と性質 やらない夫 そうそう,それを説明してなかった.これまではずっと連続時間信号について話をしてきたけど,次回から離散時間信号についての話に入るんだ.つまり,時間軸上で飛び飛びの時刻にしか値をもたないような信号だな.いよいよ「ディジタル」信号処理の世界に入っていくわけだ. 離散時間信号について考えるとき,いわゆる普通の角周波数とは別に「正規化角周波数」という概念が出てくる.小文字の はそっちの方で使おうと思うんだ.だからそれと区別するために,普通の角周波数は と書くことにする.ちょっと戸惑うかもしれないが,まあ我慢してくれ. やらない夫 互いにフーリエ変換と逆変換の関係になっているものを「フーリエ変換対」と呼ぶことがある.後々の説明で必要になるものをいくつか計算しておこうと思う. やる夫 あまり計算好きじゃないお. やらない夫 まあ数学の演習じゃないので,必要最低限に留めようと思う. そうそう,以下では が のフーリエ変換であることをこんな風に表すことにする.これらは割と標準的な記法だ. 26 単にこういうことかお. やらない夫 いいだろう.結局,どんな形のスペクトルになるかわかるか? 2 デルタ関数と複素指数関数 やる夫 デルタ関数…数学の授業で習った気はするお. やらない夫 正確にはディラックのデルタ関数とか,あるいは単位インパルス関数と呼ばれることもあるが,どんなものだったか覚えているか? ってのが今回の話になる.結論からいうと,それがフーリエ変換だ. やる夫 「級数」が「変換」に変わるんかお.なんか「周期的」かどうかとは全く異質な話に聞こえるお. やらない夫 そうかもな.まあその辺は追々理解してもらえばいい.ともかく出発地点はフーリエ級数だ.周期 の時間信号を周波数成分に分解するんだった.どんな周波数成分が出てくる? やらない夫 うーん,実在という言葉の意味によるな.例えば音声信号とか電気信号として物理的に存在するかというと,振幅が無限大なんてのは無理だから,存在はしない.数学的にも,普通の「関数」としては存在しないと考えた方がいいだろうな.しかし,超関数という概念を導入することで,ちゃんと定義することができる.だから数学的にはちゃんと実在しているともいえる. やる夫 超関数って,その中二病っぽい響きの単語は何なんだお.そういう難しそうなのは勘弁してほしいお. やらない夫 まあそこに深入りする気はないので安心してくれ.とにかく信号処理を考える上ではものすごく重要な概念なので,やや天下りだが,こういうような性質をもった「関数っぽいもの」が存在すると考えればいい. やる夫 そうしますお. やらない夫 というわけで,デルタ関数をフーリエ変換するとどうなるか,というのがここからの主題だ.早速だが, をフーリエ変換してみたらどうなる? やらない夫 いい質問だな.じゃあ,三角関数や複素指数関数に限らず一般論として,周期 で周期的な信号 をフーリエ変換したときにどうなるかを考えてみようか. 3. import numpy as np import matplotlib. から,式 によって元の が復元できる.この計算をフーリエ逆変換と呼ぶ. あるいは「 は のフーリエ逆変換である」という言い方もする• は角周波数を表す連続変数である. は に含まれる角周波数 の振動成分の量 振幅・位相 を表す.• やる夫 は に関わらず 1 だお.あ,そうか,だからあらゆる周波数が等しい振幅で含まれているっていう点では と同じなんだお.問題は位相だお. だから,周波数 の成分は, だけ位相が遅れていることになるお.これってどういうことだお…. やらない夫 もう一度時間領域に戻って考えてみるといいぞ. ってのは, から時刻 だけシフトしているわけだろ.そういう関数を合成するには,各周波数のサイン波をどうすればいいと思う? エクセルを用いたフーリエ変換 FFT エクセルを用いたフーリエ変換 FFT Excelを用いた科学技術計算が第2版になりました 30年10月 ! 大学の教科書にも採択されています amazon: 楽天: 図-1 フーリエ変換を行う例題の波形グラフ 音波、電磁波、地震波などの波は大きさ(振幅)、周波数、位相が異なる三角関数波 sin,cos の組み合わせで表すことができる。

20
やらない夫 そうだな.「重ね合わせ」という言葉であれば,総和の場合も積分の場合も,まあそんなに違和感無く表現できてる気がするが,どうだろう.まあ語感は人それぞれかも知れないけどな. ともかく,一般の時間信号は,あらゆる実数を周波数とする複素指数関数の重ね合わせで表すことができる,ということだ.これがフーリエ逆変換の意味だ. やる夫 逆? やらない夫 そうそう,それを説明してなかった.これまではずっと連続時間信号について話をしてきたけど,次回から離散時間信号についての話に入るんだ.つまり,時間軸上で飛び飛びの時刻にしか値をもたないような信号だな.いよいよ「ディジタル」信号処理の世界に入っていくわけだ. 離散時間信号について考えるとき,いわゆる普通の角周波数とは別に「正規化角周波数」という概念が出てくる.小文字の はそっちの方で使おうと思うんだ.だからそれと区別するために,普通の角周波数は と書くことにする.ちょっと戸惑うかもしれないが,まあ我慢してくれ. やらない夫 互いにフーリエ変換と逆変換の関係になっているものを「フーリエ変換対」と呼ぶことがある.後々の説明で必要になるものをいくつか計算しておこうと思う. やる夫 あまり計算好きじゃないお. やらない夫 まあ数学の演習じゃないので,必要最低限に留めようと思う. そうそう,以下では が のフーリエ変換であることをこんな風に表すことにする.これらは割と標準的な記法だ. 3. あ,そうか,フーリエ級数展開に対応するのは,フーリエ変換じゃなくてフーリエ逆変換の方なんだお.なんか混乱しそうだお. やらない夫 そう,フーリエ変換は,フーリエ係数の計算の方に対応している.それぞれの周波数成分がどのくらい含まれているかを知るための計算になっているということだな. は一般に複素数になるから,振幅と位相を持っている.フーリエ係数 と同様に,周波数 の成分の振幅と初期位相を表しているわけだ. やる夫 フーリエ級数展開やフーリエ係数の計算を「変換」と呼んじゃダメなのかお? 図-6 フーリエ変換結果 今回、波形のデータ数 n を128、時間刻み dt を0. 61 やる夫 ほとんど同じだお.フーリエ係数を計算するときの 倍をやめて,積分範囲が 1 周期分だったのを から までに変えたのがフーリエ変換になるお. やらない夫 そうだな.フーリエ係数を計算するときは,周期関数の 1 周期分だけを積分することで有限の値 を得ていたわけだ.フーリエ変換するときは,まず 倍が必要だがそれはさておいても,この 1 周期分の積分を無限個足し合わせなきゃならない. やる夫 んー,無限…? フーリエ変換は時間 t の関数である波形 f t を周波数 k の分布関数F k に変換し、その逆がフーリエ逆変換である。

MTFの測定-チャート法-

11

3. フーリエ変換 (やる夫で学ぶディジタル信号処理)

。 。 。

矩形 波 フーリエ 変換

15
。 。

3. フーリエ変換 (やる夫で学ぶディジタル信号処理)

。 。

。 。

矩形 波 フーリエ 変換

3. フーリエ変換 (やる夫で学ぶディジタル信号処理)

12

MTFの測定-チャート法-

17
。 。

MTFの測定-チャート法-

。 。 。

5

MTFの測定-チャート法-

。 。

。 。

エクセルを用いたフーリエ変換(FFT)

。 。

11
。 。