巴特沃斯:
1、带阻滤波器设计
带阻滤波器指标:
阻带上边界频率:5Kz;
阻带下边界频率:7Kz;
通带上边界频率:2Kz;
通带下边界频率:9Kz;
通带最大衰减:1dB;
阻带最小衰减:20dB;
设计程序如下:
wp=2*pi*[2000,9000];
ws=2*pi*[5000,7000];
Rp=1;
As=20;
[Nb,wc]=buttord(wp,ws,Rp,As,‘s‘);
[BSB,ASB]=butter(Nb,wc,‘stop‘,‘s‘);
[hk,w]=freqs(BSB,ASB);
subplot(2,2,1)
plot(w,20*log(abs(hk)),‘b‘)
xlabel(‘频率/Hz‘)
ylabel(‘幅度/dB‘)
title(‘幅频特性‘)
subplot(2,2,2)
plot(angle(hk))
xlabel(‘频率/Hz‘)
ylabel(‘相位‘)
title(‘相频特性‘)
程序运行结果:
巴特沃斯模拟带阻滤波器阶数:Nb=4
2、带通滤波器的设计
带通滤波器指标:
通带上边界频率:4Kz;
通带下边界频率:7Kz;
阻带上边界频率:2Kz;
阻带下边界频率:9Kz;
通带最大衰减:1dB;
阻带最小衰减:20dB;
设计程序如下:
wp=2*pi*[4000,7000];
ws=2*pi*[2000,9000];
Rp=1;
As=20;
[N,wc]=buttord(wp,ws,Rp,As,‘s‘);
[BB,AB]=butter(N,wc,‘s‘);
[hk,w]=freqs(BB,AB);
subplot(2,2,1)
plot(w,abs(hk),‘b‘)
xlabel(‘频率/Hz‘)
ylabel(‘幅度/dB‘)
title(‘幅频特性‘)
subplot(2,2,2)
plot(angle(hk))
xlabel(‘频率/Hz‘)
ylabel(‘相位‘)
title(‘相频特性‘)
程序运行结果:
巴特沃斯模拟带通滤波器阶数:N=5
原文:https://www.cnblogs.com/hjj-fighting/p/10324768.html