Scribd
Upload
4 views

Exp 8

The document outlines the design of various IIR digital filters using Butterworth and Chebyshev methods, including Low Pass, High Pass, Band Pass, and Band Stop filters. It provides MATLAB code for each filter type, detailing the input parameters and plotting the gain and phase responses. The document serves as a practical guide for implementing these filter designs in MATLAB.

Uploaded by

pskiran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
4 views

Exp 8

The document outlines the design of various IIR digital filters using Butterworth and Chebyshev methods, including Low Pass, High Pass, Band Pass, and Band Stop filters. It provides MATLAB code for each filter type, detailing the input parameters and plotting the gain and phase responses. The document serves as a practical guide for implementing these filter designs in MATLAB.

Uploaded by

pskiran
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 7

EXP8: IIR FILTER DESIGN

AIM: Design following IIR Digital Filters using i) Butterworth and ii) Chebyshev designs:

(a) LPF (b) HPF (c) BPF (d) BSF

1. Butterworth
a. LPF
PROGRAM:
clc;
clear all;
rp=input('enter passband ripple');
rs=input('enter stopband ripple');
wp=input('enter passband frequency');
ws=input('enter stopband frequency');
fs=input('enter sampling frequency');
w1=2*wp/fs;
w2=2*ws/fs;
[n,wn]=buttord(w1,w2,rp,rs);
[b,a]=butter(n,wn);
w=0:0.01:pi;
[h,om]=freqz(b,a,w);
m=20*log10(abs(h));
an=angle(h);
subplot(2,1,1);
plot(om/pi,m);
ylabel('gainin db');grid;
xlabel('normalised frequency');
subplot(2,1,2);
plot(om/pi,an);
ylabel('phase in radians');grid;
xlabel('normalisedfrequency');
b.HPF
PROGRAM:
clc;
clear all;
rp=input('enter passband ripple');
rs=input('enter stopband ripple');
wp=input('enter passband frequency');
ws=input('enter stopband frequency');
fs=input('enter sampling frequency');
w1=2*wp/fs;
w2=2*ws/fs;
[n,wn]=buttord(w1,w2,rp,rs);
[b,a]=butter(n,wn,'high');
w=0:0.01:pi;
[h,om]=freqz(b,a,w);
m=20*log10(abs(h));
an=angle(h);
subplot(2,1,1);
plot(om/pi,m);
ylabel('gain in db');grid;
xlabel('normalised frequency');
subplot(2,1,2);
plot(om/pi,an);
ylabel('phase in radians');grid;
xlabel('normalisedfrequency');
c.BPF
PROGRAM:
clc;
close all;
clear all;
alphap=2;
alphas=20;
ws=[0.1,0.5];
wp=[0.2,0.4];
[n,wn]=buttord(wp,ws,alphap,alphas);
[b,a]=butter(n,wn);
w=0:0.01:pi;
[h,ph]=freqz(b,a,w);
m=20*log(abs(h));
an=angle(h);
subplot(2,1,1);plot(ph/pi,m);grid;
ylabel('Gain in dB');xlabel('normalised frequency')
title('Magnitude response of Butterworth Band Pass Filter');
subplot(2,1,2);plot(ph/pi,an);grid;
ylabel('Phase in radians');xlabel('normalised frequency')
title('Phase response of Butterworth Band Pass Filter');

d. BSF
PROGRAM:
clc;
close all;
clear all;
alphap=2;
alphas=20;
ws=[0.2,0.4];
wp=[0.1,0.5];
[n,wn]=buttord(wp,ws,alphap,alphas);
[b,a]=butter(n,wn,'stop');
w=0:0.01:pi;
[h,ph]=freqz(b,a,w);
m=20*log(abs(h));
an=angle(h);
subplot(2,1,1);plot(ph/pi,m);grid;
ylabel('Gain in dB');xlabel('normalised frequency')
title('Magnitude response of Butterworth Band Reject Filter');
subplot(2,1,2);plot(ph/pi,an);grid;
ylabel('Phase in radians');xlabel('normalised frequency')
title('Phase response of Butterworth Band Reject Filter');

TO DESIGN A CHEBYSHEV-1 LOW PASS FILTER:-


PROGRAM:-
clc;
close all;
clear all;
alphap=1;
alphas=15;
wp=0.2;
ws=0.3;
[n,wn]=cheb1ord(wp,ws,alphap,alphas);
[b,a]=cheby1(n,alphap,wn);
w=0:0.01:pi;
[h,ph]=freqz(b,a,w);
m=20*log(abs(h));
an=angle(h);
subplot(2,1,1);plot(ph/pi,m);grid;
ylabel('Gain in dB');xlabel('normalised frequency')
title('Magnitude response of IIR Chebshev-1 Low Pass Filter');
subplot(2,1,2);plot(ph/pi,an);grid;
ylabel('Phase in radians');xlabel('normalised frequency')
title('Phase response of Chebyshev-1 Low Pass Filter');
TO DESIGN A CHEBYSHEV TYPE 1 HIGH PASS FILTER:-
clc;
clear all;
close all;
alphap=0.4;
alphas=30;
wp=0.4;
ws=0.8;
[n,wn]=cheb1ord(wp,ws,alphap,alphas);
[b,a]=cheby1(n,alphap,wn,'high');
w=0:0.01:pi;
[h,ph]=freqz(b,a,w);
m=20*log(abs(h));
an=angle(h);
subplot(2,1,1);
plot(ph/pi,m);grid
ylabel('Gain in dB');xlabel('normalised frequency')
title('Magnitude response of CHEBYSHEV I High Pass Filter ');
subplot(2,1,2);
plot(ph/pi,an);grid
ylabel('Phase in radians');xlabel('normalised frequency')
title('Phase response of CHEBYSHEV I High Pass Filter');

TO DESIGN A CHEBYSHEV-1 BAND PASS FILTER: -


PROGRAM:-
clc;
close all;
clear all;
alphap=2;
alphas=20;
ws=[0.1,0.5];
wp=[0.2,0.4];
[n,wn]=cheb1ord(wp,ws,alphap,alphas);
[b,a]=cheby1(n,alphap,wn);
w=0:0.01:pi;
[h,ph]=freqz(b,a,w);
m=20*log(abs(h));
an=angle(h);
subplot(2,1,1);plot(ph/pi,m);grid;
ylabel('Gain in dB');xlabel('normalised frequency')
title('Magnitude response of Chebyshev-1 Band Pass Filter');
subplot(2,1,2);plot(ph/pi,an);grid;
ylabel('Phase in radians');xlabel('normalised frequency')
title('Phase response of Chebyshev-1 Band Pass Filter');

TO DESIGN A CHEBYSHEV-2 BAND REJECT FILTER:-


PROGRAM:-
clc;
close all;
clear all;
alphap=2;
alphas=20;
ws=[0.2,0.4];
wp=[0.1,0.5];
[n,wn]=cheb2ord(wp,ws,alphap,alphas);
[b,a]=cheby2(n,alphas,wn,'stop');
w=0:0.01:pi;
[h,ph]=freqz(b,a,w);
m=20*log(abs(h));
an=angle(h);
subplot(2,1,1);plot(ph/pi,m);grid;
ylabel('Gain in dB');xlabel('normalised frequency')
title('Magnitude response of Chebyshev-2 Band Reject Filter');
subplot(2,1,2);plot(ph/pi,an);grid;
ylabel('Phase in radians');xlabel('normalised frequency')
title('Phase response of Chebyshev-2 Band Reject Filter');

You might also like