我在试着做一个带通滤波器。在我工作的时候,我注意到了一些奇怪的东西。
滤波器的频率响应正常(scipy.signal.freqz公司). 但是当我把它应用到正弦信号上时,滤波器的形状会改变。在
图3=信号频谱
图4=根据freqz()函数的滤波器频率响应
图1=滤波器应用信号的频谱
最后,为什么滤波器的形状(也是幅度)在信号上应用时会发生变化?在
from scipy.signal import butter, lfilter, buttord, freqz
import numpy as np
import matplotlib.pyplot as plt
import math
t = np.arange(0, 8192, 1)
data = 1+np.sin(.15*np.pi*t + np.pi/2)
b, a = butter(6, [0.35, 0.55], btype='band')
y = lfilter(b, a, data)
y = np.fft.rfft(y)
y = np.abs(y)
x = np.fft.rfft(data)
x = np.abs(x)
w, h = freqz(b, a, worN=4096)
plt.figure(1)
plt.plot(np.arange(0,4097,1),y)
plt.figure(2)
plt.plot(np.arange(0,8192,1), data)
plt.figure(3)
plt.plot(np.arange(0,4097,1), x)
plt.figure(4)
plt.plot(np.arange(0,4096,1),np.abs(h))
plt.show()