%FilterCompare.m
%对比相同阶数时四种IIR滤波器的波纹特性


close all;
clear;
clc;


fs = 1000; %Hz


%低通滤波器对比
N1 = 8;
wp = 50 / ( fs / 2);
ws = 100 / ( fs / 2);
alpha_p = 1;
alpha_s = 20;

[ b1,a1] = butter( N1 , ws , 'low' );
[ b2,a2] = cheby1( N1 , alpha_p , ws , 'low' );
[ b3,a3] = cheby2( N1 , alpha_s , ws , 'low' );
[ b4,a4] = ellip ( N1 , alpha_p,alpha_s , ws , 'low' );



[H1,w1] = freqz(b1,a1,1024*4);
[H2,w2] = freqz(b2,a2,1024*4);
[H3,w3] = freqz(b3,a3,1024*4);
[H4,w4] = freqz(b4,a4,1024*4);
w1 = w1/pi *fs/2;
w2 = w2/pi *fs/2;
w3 = w3/pi *fs/2;
w4 = w4/pi *fs/2;
figure(1);
plot(w1,abs(H1),'black','linewidth',0.5);
hold on;
plot(w2,abs(H2),'red','linewidth',0.5);
plot(w3,abs(H3),'green','linewidth',0.5);
plot(w4,abs(H4),'blue','linewidth',0.5);
xlim([0 500]);
ylim([0 1.1]);
legend({ 'butter' , 'cheby1' , 'cheby2' , 'ellip'});
title('同阶数(N=8),同截止频率(f_c=100Hz)的低通滤波器对比');


%高通滤波器
N1 = 8;
wp = 400 / ( fs / 2);
ws = 400 / ( fs / 2);
alpha_p = 1;
alpha_s = 20;
[ b1,a1] = butter( N1 , ws , 'high' );
[ b2,a2] = cheby1( N1 , alpha_p , ws , 'high' );
[ b3,a3] = cheby2( N1 , alpha_s , ws , 'high' );
[ b4,a4] = ellip( N1 , alpha_p,alpha_s , ws , 'high' );
[H1,w1] = freqz(b1,a1,1024*4);
[H2,w2] = freqz(b2,a2,1024*4);
[H3,w3] = freqz(b3,a3,1024*4);
[H4,w4] = freqz(b4,a4,1024*4);
w1 = w1/pi *fs/2;
w2 = w2/pi *fs/2;
w3 = w3/pi *fs/2;
w4 = w4/pi *fs/2;
figure(2);
plot(w1,abs(H1),'black','linewidth',0.5);
hold on;
plot(w2,abs(H2),'red','linewidth',0.5);
plot(w3,abs(H3),'green','linewidth',0.5);
plot(w4,abs(H4),'blue','linewidth',0.5);
xlim([0 500]);
ylim([0 1.1]);
legend({ 'butter' , 'cheby1' , 'cheby2' , 'ellip'});
title('同阶数(N=8),同截止频率(f_c=400Hz)的高通滤波器对比');



%带通滤波器

N1 = 8;
wp = [ 100 , 400 ] / ( fs / 2);
ws = [ 100 , 400 ] / ( fs / 2);
alpha_p = 1;
alpha_s = 20;
[ b1,a1] = butter( N1 , ws , 'bandpass' );
[ b2,a2] = cheby1( N1 , alpha_p , ws , 'bandpass' );
[ b3,a3] = cheby2( N1 , alpha_s , ws , 'bandpass' );
[ b4,a4] = ellip( N1 , alpha_p,alpha_s , ws , 'bandpass' );
[H1,w1] = freqz(b1,a1,1024*4);
[H2,w2] = freqz(b2,a2,1024*4);
[H3,w3] = freqz(b3,a3,1024*4);
[H4,w4] = freqz(b4,a4,1024*4);
w1 = w1/pi *fs/2;
w2 = w2/pi *fs/2;
w3 = w3/pi *fs/2;
w4 = w4/pi *fs/2;
figure(3);
plot(w1,abs(H1),'black','linewidth',0.5);
hold on;
plot(w2,abs(H2),'red','linewidth',0.5);
plot(w3,abs(H3),'green','linewidth',0.5);
plot(w4,abs(H4),'blue','linewidth',0.5);
xlim([0 500]);
ylim([0 1.1]);
legend({ 'butter' , 'cheby1' , 'cheby2' , 'ellip'});
title('同阶数(N=8),同截止频率(f_c_1=100Hz f_c_2 = 400Hz)的带通滤波器对比');



%带阻滤波器
N1 = 8;
wp = [ 100 , 400 ] / ( fs / 2);
ws = [ 100 , 400 ] / ( fs / 2);
alpha_p = 1;
alpha_s = 20;
[ b1,a1] = butter( N1 , ws , 'stop' );
[ b2,a2] = cheby1( N1 , alpha_p , ws , 'stop' );
[ b3,a3] = cheby2( N1 , alpha_s , ws , 'stop' );
[ b4,a4] = ellip( N1 , alpha_p,alpha_s , ws , 'stop' );
[H1,w1] = freqz(b1,a1,1024*4);
[H2,w2] = freqz(b2,a2,1024*4);
[H3,w3] = freqz(b3,a3,1024*4);
[H4,w4] = freqz(b4,a4,1024*4);
w1 = w1/pi *fs/2;
w2 = w2/pi *fs/2;
w3 = w3/pi *fs/2;
w4 = w4/pi *fs/2;
figure(4);
plot(w1,abs(H1),'black','linewidth',0.5);
hold on;
plot(w2,abs(H2),'red','linewidth',0.5);
plot(w3,abs(H3),'green','linewidth',0.5);
plot(w4,abs(H4),'blue','linewidth',0.5);
xlim([0 500]);
ylim([0 1.1]);
legend({ 'butter' , 'cheby1' , 'cheby2' , 'ellip'});
title('同阶数(N=8),同截止频率(f_c_1=100Hz f_c_2 = 400Hz)的带阻滤波器对比');