Pythonでのハイパスフィルタが実装できない

Question

pythonを使って音声解析をしています。
現在1秒ごとに波形の表示とその波形のフーリエ変換を行っています。
以下のコードが該当するのですが、300ヘルツ以下をカットした状態でプロットしたいです。
ハイパスフィルタというのですがPython初心者でありどのようにコードを変更すればよいのかわかりません。
わかる方がいらっしゃいましたら変えるコードを教えていただきたいです。
何卒宜しくお願い致します。

以下のサイトのコードを参考に実行しようとしたのですが、このコードを実行すると300ヘルツ以下をカットしたいのに300以上がカットされてしまいます。
https://algorithm.joho.info/programming/python/numpy-fast-fourier-transform/#toc3

# -*- coding: utf-8 -*-
import wave

from matplotlib.colors import same_color
import numpy as np
import matplotlib.pyplot as plt 
 
def main():
    filename = "04-agarikata-kara-01-p6TAskMJUgg.wav"
    wf = wave.open(filename, "r" ) 
     
    # WAVファイルの情報を表示
    print ("オーディオチャンネル数(モノラル: 1 ステレオ:2 ) : ", wf.getnchannels())
    print ("サンプルサイズ(バイト数) : ", wf.getsampwidth())
    print ("サンプリングレート : ", wf.getframerate())
    print ("オーディオフレーム数 : ", wf.getnframes())
    print ("圧縮形式 : ", wf.getcomptype())
    print ("（圧縮形式） : ", wf.getcompname())
    print ("パラメータ : ", wf.getparams())
    print ("記録時間(Sec) : ", float(wf.getnframes()) / wf.getframerate() )       
     
    fs = wf.getframerate()                         # サンプリング周波数

    fc = 300                                     # カットオフ周波数

    g = wf.readframes(wf.getnframes())
    g = np.frombuffer(g, dtype= "int16")/32768.0    # -1～1に正規化

    mado = 1#s
    N = mado*fs                                 # サンプル数＝窓（1秒）×サンプリング周波数

    freq = np.linspace(0, 1.0/mado, N)
    
    wf.close()
    for ki in range(round(wf.getnframes()/fs)-1):
        n0 = ki                                     # サンプリング開始位置
        gz = g[::2]                                 # サンプル数片方にわける
        gg = gz[ki*N:ki*N+N]
        # print(fs,ki,N,gz,gg)
        G = np.fft.fft(gg)                          # 高速フーリエ変換
        amp = [np.sqrt(c.real ** 2 + c.imag ** 2) for c in G]       # 振幅スペクトル
        phase = [np.arctan2(int(c.imag), int(c.real)) for c in G]   # 位相スペクトル
        flist = np.fft.fftfreq(N, d=1.0/fs)             # 周波数リスト
        flist[(freq < fc)] = 0
        flist[(freq > 1/(mado*2))] = 0
        
        # 波形サンプルを描画
        fig = plt.figure()
        plt.subplot(211)
        plt.plot(np.linspace(n0,n0+mado,N),gg)
        plt.axis([n0, n0+mado, -0.2, 0.2])
        plt.xlabel("時間 [sample]",fontname="MS Gothic")
        plt.ylabel("振幅",fontname="MS Gothic")

        # 振幅スペクトルを描画
        plt.subplot(212)
        plt.plot(flist, amp, marker='o', linestyle='-')
        plt.axis([0, fs/2, 0, int(max(amp))])
        plt.xlabel("周波数 [Hz]",fontname="MS Gothic")
        plt.ylabel("振幅スペクトル",fontname="MS Gothic")
        output ="C:\Users\hiro2\Desktop\sound\kara"
        fig.savefig(output+"\"+str(ki)+"_"+str(mado)+".jpg")
        plt.close()
        fig = plt.figure()
        plt.close()
 
if __name__ == '__main__':
    main()

Comments

[y_waiwai]

このままではコードが読めないので、質問を編集し、<code>ボタンを押し、出てくる’’’の枠の中にコードを貼り付けてください

[hirohiro1999]

ありがとうございます。
直します！

[1T2R3M4]

質問テンプレートにあったと思いますが、
調べたこと、試したことを追記していただけませんか。

[hirohiro1999]

はい！

[jbpb0]

pythonのコードの一番最初の行のすぐ上に
```python
だけの行を追加してください
また、pythonのコードの一番最後の行のすぐ下に
```
だけの行を追加してください

現状、コードがとても読み辛いです
質問にコードを載せる際に上記をやってくれたら、他人がコードを読みやすくなり、コードの実行による現象確認もやりやすくなるので、回答されやすくなります

[hirohiro1999]

すみません。
投稿が不慣れで。。。
このような形でよろしでしょうか？

[jbpb0]

大丈夫です

[jbpb0]

flist[(freq < fc)] = 0
flist[(freq > 1/(mado*2))] = 0
↓ 変更
amp[((flist < fc) & (flist > -fc))] = 0

で、どうでしょうか？

[hirohiro1999]

替えてみます！

[jbpb0]

どうでしたか？

[hirohiro1999]

うまくいきました！
ありがとうございました！！！

Answer 1

flist[(freq < fc)] = 0
flist[(freq > 1/(mado*2))] = 0

↓ 変更

amp[((flist < fc) & (flist > -fc))] = 0