滤波器相关知识

Posted 2022-09-12 windyz

滤波器相关

理想滤波器幅频特性中通带和阻带的增益均为0

实际滤波器通带与阻带之间存在过渡带（迁移带）

数字滤波器进行滤波：需要对模拟信号通过ADC进行A/D转换称数字信号后在进行数字滤波，滤波后在进行D/A转换回模拟信号

但ADC之前要通过抗混叠滤波器，一般采用模拟滤波器

频率响应：幅频响应+相频响应

单位冲激响应--傅里叶变换-->频率响应

当输出信号幅值降低至最大值的1/sqrt(2)=0.707倍时，对应的频率称截止频率

当输出下降速度趋于稳定时，频率每增大10倍，幅值对应降低10倍，即减小20dB

伯德图就是坐标轴取了对数的幅频特性曲线

其优点在于横轴频率范围大，纵轴频率范围大，将频率响应不规则曲线转变成规则曲线

n阶巴特沃斯低通滤波器幅频响应图：

$ |H(w)|=\frac1\sqrt1+(\fracww_c)^2n$

w_c是滤波器截止频率，n是滤波器阶数

# 连载481： 1、2、3阶巴特沃斯低通滤波器波特图
omega = np.arange(1,500000)
A = 1/np.sqrt(1+(omega/5000)**2)
plt.loglog(omega,A,label='First-order Butterworth Filter')
A = 1/np.sqrt(1+(omega/5000)**4)
plt.loglog(omega,A,label='Second-order Butterworth Filter')
A = 1/np.sqrt(1+(omega/5000)**6)
plt.loglog(omega,A,label='Third-order Butterworth Filter')
plt.legend(loc="best",fontsize=10)
plt.xlabel('rad/s')
plt.ylabel('dB')
plt.show()

n阶巴特沃斯高通滤波器幅频响应图：

$ |H(w)|=\frac1\sqrt1+(\fracw_cw)^2n$

w_c是滤波器截止频率，n是滤波器阶数

# high-pass filter
omega = np.arange(1,500000)
H = 1/(1+1j*5000/omega)
A = abs(H)
plt.loglog(omega,A)

n阶巴特沃斯带通滤波器幅频响应图：

$ |H(w)|=\frac1\sqrt1+(\fracw^2-w_0Bw)^2n$

w_0是滤波器中心频率，B是滤波器带宽，n是滤波器阶数

截止频率\(\omega_c = \frac\sqrtB^2+4w_0^2\pm B2\)

# band-pass filter
omega = np.arange(0.1,100000)
H = 1/(1+1j*(omega**2-100**2)/(200*omega))
A = abs(H)
plt.loglog(omega,A)
plt.xlabel('$\omega$(rad/s)')
plt.ylabel('|H(w)| (dB)')

连载483，484，485待补充