scipy.signal.firwin()の周波数応答を表示するプログラムが書きたいです

Question

カットオフ周波数が0.75Hz,4.0Hzのバンドパスフィルタを用いて信号解析をしていて、このフィルタの周波数応答をグラフにしたいです。
どうしたらいいでしょうか。

bandpassfilter = signal.firwin(numtaps=65, cutoff= [0.75,4.0], pass_zero=False, nyq=1/(2 * dt))

Answer 1

こんな感じでどうでしょう

import numpy as np
import matplotlib.pyplot as plt

f = np.linspace(0, 10, 1000)
z = np.exp(2j * np.pi * f * dt)
g = (bandpassfilter[None] * z[:, None] ** -np.arange(bandpassfilter.size)[None]).sum(1)

plt.figure()

plt.subplot(211)
plt.plot(f, np.abs(g))
plt.xlabel('Frequency [Hz]')
plt.ylabel('Gain [-]')

plt.subplot(212)
plt.plot(f, np.rad2deg(np.angle(g)))
plt.xlabel('Frequency [Hz]')
plt.ylabel('Phase [rad]')

plt.tight_layout()
plt.show()

Comments

[YukaU]

ありがとうございます！！！！