例によってあの本 を読み直しの最中。1回目読んで
わかりにくかった部分もだいぶ読み下せるように
なってきた感じ。
で、ＦＦＴを掛けると計算できるＤＦＴ値とスペクトル
について、ふと素朴な疑問が湧いてきたのでＰＣ上で
シミュレーションを掛けてみたいと。
ある周波数のｓｉｎ波が有った場合に、ＦＦＴを掛ける
際の周波数分解能の兼ね合いから、基本波×整数倍の
周波数にキッチリ乗らない場合もあるわけですが、
1回目に読んだ時の理解では、そういう風に基本波の
整数倍波の狭間に落ちたスペクトルは、両側に
分散して成分が乗る…すなわち、ＬＰＦとＨＰＦを
細かい周波数ごとに掛けていって、結果前後の
周波数成分のところに分散されて片寄せされる
ということになっているっぽいという理解をしました。
そういうフィルターの集まりとしての機能があると…
（今読み直しているところなので、ホントにその
理解で正しいのか自信が無いんだけど）
ってことでこの疑問に関する実験。ある角速度の
ｓｉｎ波を作っておいて、それに近い上下の速度の
波も作ってみて、上下の波を足して２で割ったら
（つまり平均取ったら）元の波に戻るのか？って
話。もし戻れば、上記の話はそのままＯＫなんだけど、
位相が微妙にずれていく波を足し合わせたらその
中間の周波数の波になるんだっけ？怪しい…
ってことで、ＰＣ画面上で実際にグラフに表示して
見ることに。
（１）元の波
（２）（１）よりちょっとだけ遅い波
（３）（１）よりちょっとだけ速い波
（４）（２）と（３）の平均
をグラフに描いてみます。それが↓これ。
（１）と（４）が一致してない。ウネリ？ウナリ？
が出てる。
…よく考えたら当たり前。
http://ja.wikipedia.org/wiki/%E3%81%86%E3%81%AA%E3%82%8A
これだ。
うーん、もうちょっと読み込まないといけませんな…
（追記）
よく考えてみたら、あれかな？
実数部分だけで考えてるからダメなのかな？
今回の件で複素数の計算は少しずつわかって
きたつもりだったけど、アタマで解って
きたってかんじで、まだ体で解ったというレベル
ではないんだろうな。
やっぱドリルとか解いて体で覚えないと
応用するにはまだ一歩って感じなのかもな…
ＦＦＴ自由自在、ＦＦＴ応用自在が欲しいな。
|
おすすめブログ
@goo_blog
