帶通濾波 matlab

巴特沃斯：
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