MATLABでFFTしたwavフィルタにハイパスフィルタをかけたい

Question

前提・実現したいこと

今回下の画像のようなFFTを作成し、低周波に出ているピークを取り除きたくハイパスフィルタを使おうと思いこちらのサイトを見つけ、ローパスフィルタのプログラムではありますが、参考にさせていただくことにしました。wavファイルをcsvにファイルにした際に、そのサイズを小さくするためにハイパスフィルタを用いることにしました。

https://qiita.com/tomochiii/items/2145ee7a23f5dd4cff9a

このサイトには３つのプログラムが掲載されており１つ目はクリアできているものとします。lowpass.mとaudio_filter.mに関しては、どちらも実行してみましたが、どちらも正規化周波数のグラフがプロットされ、目標となる下の画像のハイパスフィルタでピークを取り除いた画像になりません。FFTした画像の下がプロットされた正規化周波数の画像です。

power_Spectrum.mのプログラムでのsignalは、wavファイルをcsvファイルにしてそちらを配列にして入力しました。
サンプリング周波数は１６０００です。

教えていただきたいこと

・lowpass.mでのプログラムの流れを詳しく教えていただきたいです
・なぜ、lowpass.mとaudio_filter.mを実行した際に正規化周波数のグラフがプロットされてしまうのでしょうか

補足情報（FW/ツールのバージョンなど）

バージョン
MATLB R2020b

Comments

[jbpb0]

> 低周波に出ているピークを取り除きたくローパスフィルタを使おうと思い

「低周波」のピークを下げたいんですよね?
ローパスフィルタは、「低周波」はそのまま残して、「高周波」の振幅を減らすフィルタですよ

下記Webページに、ローパスフィルタ適用前後のf特のグラフがあるので、それ見てください
四つ並んでるグラフの、右上が元のf特、右下がローパスフィルタ適用後のf特です
https://algorithm.joho.info/programming/python/numpy-ifft-lowpass-denoise/

[ilbisonte]

ご指摘ありがとうございます。勘違いしていました。
すぐに質問を編集いたします。

[yuki23]

> どちらも正規化周波数のグラフがプロットされ、目標となる下の画像のハイパスフィルタでピークを取り除いた画像になりません
具体的にどういうものが表示されるんでしょうか？勝手に要約するのではなく質問文に貼ってください

[ilbisonte]

申し訳ありません
具体的な画像を提示しました、

[ozwk]

audio_filter.mを読んでいるといきなりFs_sineという未定義の変数がpowerSpectrum()に与えられているように見えますが、動いてますか?

Answer 1

audio_filter関数にFs_sineという未定義の変数が出てくるのでエラーで強制終了しているからです。
Fsに変えて実行してみてください。

lowpass.mでのプログラムの流れを詳しく教えていただきたいです

たった3行の関数の何を説明すればいいんでしょうか？ 質問は具体的にお願いします。

Comments

[ilbisonte]

回答ありがとうございます。
Fsに変えて実行してみましたが、ここでもう一つ質問させてください。
audio_filter関数でのfilenameとはwaveファイルで良いのでしょうか。

曖昧質問で申し訳ありません。lowpass.mのの説明というより、なぜ正規化周波数のグラフがプロットされたのかが知りたいです。

[yuki23]

> audio_filter関数でのfilenameとはwaveファイルで良いのでしょうか。
はい

> なぜ正規化周波数のグラフがプロットされたのかが知りたいです。
freqz(b,a)
は、出力引数なしなのでグラフをプロットします https://jp.mathworks.com/help/signal/ref/freqz.html