« TULLY'S COFFEE 海岸竹橋通り店 | トップページ | デニーズでウバ茶 »

Rでフィルター

Octaveでフィルターをやったので、次にRを使ってみました。
元データは周波数解析 Fast Viewerのdata5000.xlsを使わせていただきました。
さらにフーリエ変換の実際も参考にさせていただいています。

#データ読み込み
DATA <- scan("data.prn")
#FFT計算
FFT <- fft(DATA)
#パワー計算
MOD <- Mod(FFT)
#FFT結果のプロット(横軸がきちんと周波数にしていません)
plot(MOD, xlim=c(0, 2500))
#フィルター
FFT[51] <- 0+0i
FFT[4951] <- 0+0i
#逆FFT計算
INFFT <- (fft(FFT, inverse = TRUE)/5000)
#結果のプロット
par(mfcol = c(2,1))
plot(DATA, ylim=c(-1,1), type="line")
plot(Re(INFFT), ylim=c(-1,1), type="line")

最後のプロットはこんな感じです。→View image
ちゃんと除去されていますね。

|

« TULLY'S COFFEE 海岸竹橋通り店 | トップページ | デニーズでウバ茶 »

コメント

コメントを書く



(ウェブ上には掲載しません)


コメントは記事投稿者が公開するまで表示されません。



トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/44043/1702557

この記事へのトラックバック一覧です: Rでフィルター:

« TULLY'S COFFEE 海岸竹橋通り店 | トップページ | デニーズでウバ茶 »