Rでフーリエ変換した後に特定周波数をカットしたい

質問

現在Rを用いてフーリエ変換の勉強をしております。下記のサイトなどを参考に、フーリエ変換は出来るようになったのですが、特定の周波数成分だけカットしてフーリエ逆変換を行い、その周波数成分がカットされた信号を作りたいのですが、どのようにするのがRでは良いでしょうか？
　下記のサイトなどを参考にしました。
　http://www.abstractnew.com/2014/04/the-fast-fourier-transform-fft-without.html

　ご存じの方おられましたらご教示願えれば幸いです。

回答No.1

DCI4

特定の周波数成分だけカットしてフーリエ逆変換を行い、
その周波数成分がカットされた信号を作りたいのですが、
どのようにするのが良いでしょうか？

★回答

Rは知らん　でもみんな同じ　いくらでも便利計算ツールはある　メーカーもいっぱい

フーリエ変換は出来るようになったのですが→では　逆変換も出来るということです

アルゴリズムは逆も同じ　対称変換対だ

離散フーリエ変換　数列が出てくるだけ　

ｈ（ｎ）⇔H(ｍ)

ｔの時系列　の式が　ｈ（ｎ）なら

FFTすれば　H(ｍ）のωの　複素数　系列が出てくる

0～ωs/2　が帯域だから　それの繰り返し

好きなωのデーター範囲をゼロにして　逆変換でOK

H(ｍ）で　　ｍ＝10～100　まで　ゼロ　　とか　すきにやるだけ

ω＝2πｆ　の式

ωs＝サンプリング角周波数　　ｆｓ；サンプリング周波数

ωは　系列　0・・・・・・・ωs　まで

　　　　　　ｍ0・・・・・ｍωs

　　　　　　ｍ＝1　　2　　3　・・・・・・　FFTの個数まで　　（複素周波数数列）

　　　　　　ｎ＝　1　　2　　3　・・・・・・　FFTの個数まで　　（時間数列）

とか出てくるだけ

それが　ｎ　の系列で出てくるだけ

時間域　ｈ（ｎ）の数列が出てくるだけ

系列の間隔は　時間域　Tｓ

周波数域間隔は　1/Ts＝Fs　　　　（ωs＝2πFs）

※注意　最初にFFTする前に　対象の時系列の信号が1/2×ωsで帯域制限されてないとだめよ意味なし

・・・・・・・・・・・ざっくり回答　おしまい・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・

グラフの書き方　御参考以下アンサー
http://okwave.jp/qa/q9069692.html
http://okwave.jp/qa/q9080187.html

お礼

ありがとうございます。しかし、Rでやるにはどうすると良いのかという質問なのでその点を踏まえてお願い致します。