MATLAB program for the design and implementation of Chebyshev low pass filter

Program code

%Chebyshev low pass filter

clc;

clear all;

close all;

fp=input('Please enter the pass band frequency: ');

fs=input('Please enter the stop band frequency: ');

rp=input('Please enter the pass band attenuation: ');

rs=input('Please enter the stop band attenuation: ');

sf=input('Please enter the sampling frequency: ');

wp=2*fp/sf;

ws=2*fs/sf;

[n,wn]=cheb1ord(wp,ws,rp,rs);

[b,a]=cheby1(n,rp,wn,'low');

freqz(b,a);

title('Chebyshev low pass filter');

Example of output

Please enter the pass band frequency: 2000

Please enter the stop band frequency: 1600

Please enter the pass band attenuation: .4

Please enter the stop band attenuation: 50

Please enter the sampling frequency: 8000 

_________________________________________

MATLAB program for the design and implementation of Chebyshev high pass filter

Program code

%Chebyshev high pass filter

clc;

clear all;

close all;

fp=input('Please enter the pass band frequency: ');

fs=input('Please enter the stop band frequency: ');

rp=input('Please enter the pass band attenuation: ');

rs=input('Please enter the stop band attenuation: ');

sf=input('Please enter the sampling frequency: ');

wp=2*fp/sf;

ws=2*fs/sf;

[n,wn]=cheb1ord(wp,ws,rp,rs);

[b,a]=cheby1(n,rp,wn,'high');

freqz(b,a);

title('Chebyshev high pass filter');

Example of output

Please enter the pass band frequency: 2000

Please enter the stop band frequency: 1600

Please enter the pass band attenuation: 0.3

Please enter the stop band attenuation: 40

Please enter the sampling frequency: 8000

__________________________

MATLAB program for the design and implementation of Chebyshev band pass filter

Program code

%Chebyshev band pass filter

clc;

clear all;

close all;

fp1=input ('Please input the first pass band frequency = ');

fs1=input ('Please input the first stop band frequency = ');

fp2=input ('Please input the second pass band frequency = ');

fs2=input ('Please input the second stop band frequency = ');

rp=input ('Please input the pass band attenuation = ');

rs=input ('Please input the stop band attenuation = ');

sf=input ('Please input the sampling frequency = ');

p1=2*fp1/sf;

s1=2*fs1/sf;

p2=2*fp2/sf;

s2=2*fs2/sf;

p=[p1,p2];

s=[s1,s2];

[n,w]=cheb1ord(p,s,rp,rs);

[f1,f2]=cheby1(n,rp,w,'bandpass');

freqz(f1,f2);

title ('Chebyshev Band Pass Filter');

Example of output

Please input the first pass band frequency = 2000

Please input the first stop band frequency = 1600

Please input the second pass band frequency = 3000

Please input the second stop band frequency = 3200

Please input the pass band attenuation = 0.5

Please input the stop band attenuation = 50

Please input the sampling frequency = 9000

_______________________

MATLAB program for the design and implementation of Butterworth band pass filter

Program code

clc;

clear all;

close all;

fp1=input ('Please input the first pass band frequency = ');

fs1=input ('Please input the first stop band frequency = ');

fp2=input ('Please input the second pass band frequency = ');

fs2=input ('Please input the second stop band frequency = ');

rp=input ('Please input the pass band attenuation = ');

rs=input ('Please input the stop band attenuation = ');

sf=input ('Please input the sampling frequency = ');

p1=2*fp1/sf;

s1=2*fs1/sf;

p2=2*fp2/sf;

s2=2*fs2/sf;

p=[p1,p2];

s=[s1,s2];

[n,w]=buttord(p,s,rp,rs);

[f1,f2]=butter (n,w,'bandpass');

freqz(f1,f2);

title ('Butterworth Band Pass Filter');

Example of output

Please input the first pass band frequency = 2000

Please input the first stop band frequency = 1700

Please input the second pass band frequency = 3000

Please input the second stop band frequency = 3200

Please input the pass band attenuation = 0.5

Please input the stop band attenuation = 60

Please input the sampling frequency = 9000

_____________________________

MATLAB program to perform linear convolution using circular convolution

Program code

%linear convolution using circular convolution

clc;

clear all;

close all;

x=input('Please enter the first sequence x[n] = ');

h=input('Please enter the second sequence h[n] = ');

len_x=length(x);

len_h=length(h);

h=[h zeros(1,len_x-1)];

x=[x zeros(1,len_h-1)];

L=len_x+len_h-1;

for i=0:L-1

    for j=1:L-i;

        a(j+i,i+1)=x(j);

    end

end

for m=1:L-1

    for n=m+1:L

        a(m,n)=x(L-n+m+1);

    end

end

b=transpose(h);

y=a*b;

y=transpose(y);

subplot(311);

stem(x);

disp('The sequence obtained after linear convolution of x[n] with h[n], is given below. ');

disp(y);

grid on;

title('First sequence x[n] ');

subplot(312);

stem(h);

grid on;

title('Second sequence h[n] ');

subplot(313);

stem(y);

Example of output

Please enter the first sequence x[n] = [1 2 3 4 5 6]

Please enter the second sequence h[n] = [1 2 3 4 5 6]

The sequence obtained after linear convolution of x[n] with h[n], is given below.

     1     4    10    20    35    56    70    76    73    60    36

________________________________________

MATLAB program for the design and implementation of Butterworth high pass filter

Program code

clc;
clear all;
close all;
fp=input('Please enter the first pass band frequency: ');
fs=input('Please enter the first stop band frequency: ');
rp=input('Please enter the pass band attenuation: ');
rs=input('Please enter the stop band attenuation: ');
sf=input('Please enter the sampling frequency: ');
wp=2*fp/sf;
ws=2*fs/sf;
[n,wn]=buttord(wp,ws,rp,rs);
[b,a]=butter(n,wn,'high');
freqz(b,a);
title('Butterworth High Pass Filter');

Example of output

Please enter the first pass band frequency: 2000

Please enter the first stop band frequency: 1300

Please enter the first pass band attenuation: 0.35

Please enter the first stop band attenuation: 50

Please enter the sampling frequency: 8000

_____________________________

MATLAB program for the design and implementation of Butterworth low pass filter

Program code

clc;
clear all;
close all;
fp=input('Please enter the first pass band frequency: ');
fs=input('Please enter the first stop band frequency: ');
rp=input('Please enter the pass band attenuation: ');
rs=input('Please enter the stop band attenuation: ');
sf=input('Please enter the sampling frequency: ');
wp=2*fp/sf;
ws=2*fs/sf;
[n,wn]=buttord(wp,ws,rp,rs);
[b,a]=butter(n,wn,'low');
freqz(b,a);
title('Butterworth Low Pass Filter');

Example of output

Please enter the first pass band frequency: 2000

Please enter the first stop band frequency: 1600

Please enter the pass band attenuation: .5

Please enter the stop band attenuation: 53

Please enter the sampling frequency: 9500

__________________________________

MATLAB program to perform circular convolution of two signals

Program Code

%circular convolution

clc;

close all;

clear all;

x=input('Please enter the first sequence x[n] = ');

h=input('Please enter the second sequence h[n] = ');

L1=length(x);

L2=length(h);

Len=max(L1,L2);

x=[x zeros(1,Len-L1)];

h=[h zeros(1,Len-L2)];

for i=0:(Len-1)

    y(i+1)=0;

    for j=0:(Len-1)

        m=mod(i-j,Len);

        y(i+1)=y(i+1)+(x(j+1)*h(m+1));

    end

end

disp('Convoluted sequence of x[n] and h[n] is given below: ');

disp(y);

subplot(311);

stem(x);

title('First sequence x[n]');

grid on;

subplot(312);

stem(h);

title('Second sequence h[n]');

grid on;

subplot(313);

stem(y);

title('Convoluted sequence ');

grid on;

Example of Output

Please enter the first sequence x[n] = [1 3 2 1]
Please enter the second sequence h[n] = [2 1 2]
Convoluted sequence of x[n] and h[n] is given below: 
     7     9     9    10

________________________________________

MATLAB program to perform linear convolution of two signals ( without using MATLAB function)

Program Code

clc;

close all;

clear all;

xt=input('Please enter the input sequence in time domain x[n]= ');

lxt=length(xt);

ht=input('Please enter the impulse sequence h[n]= ');

lht=length(ht);

ext=[xt zeros(1,(lht-1))];

eht=[ht zeros(1,(lxt-1))];

xdft=fft(ext);

hdft=fft(eht);

freqm=xdft.*hdft;

yt=[ifft(freqm)];

display('The convoluted sequence is given below');

disp(yt);

subplot(311);

stem(xt);

xlabel('Time');

ylabel('Magnitude');

title('Input sequence x[n]');

subplot(312);

stem(xt);

xlabel('Time');

ylabel('Magnitude');

title('Input sequence 1');

subplot(313);

stem(yt);

xlabel('Time');

ylabel('Magnitude');

title('Convoluted sequence');

Example of Output

Please enter the input sequence in time domain x[n]= [1 2 3 4 5 6]

Please enter the impulse sequence h[n]= [1 2 3 4 5 6]

The convoluted sequence is given below

  Columns 1 through 5 

    1.0000    4.0000   10.0000   20.0000   35.0000

  Columns 6 through 10 

   56.0000   70.0000   76.0000   73.0000   60.0000

  Column 11 

   36.0000

Matlab program to find the linear convolution of two signals (using matlab functions)

Program Code

%linear convolution (using matlab functions)
clc;
close all;
clear all;
x1=input('Please enter the input sequence x[n]= ');
x2=input('Please enter the starting time index of x[n]= ');
h1=input('Please enter the impulse response h[n]= ');
h2=input('Please enter the starting time index of h[n]= ');
y=conv(x1,h1);
n=x2+h2:length(y)+x2+h2-1;
display('The convoluted sequence is given below:');
y
stem(n,y);
xlabel('Time');
ylabel('Amplitude');
title('Linear Convolution');

Example of Output

Please enter the input sequence x[n]= [4 3 1 2]
Please enter the starting time index of x[n]= -2
Please enter the impulse response h[n]= [1 4 3]
Please enter the starting time index of h[n]= -1
The convoluted sequence is given below:

y =


     4    19    25    15    11     6

_________________________________