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 find the Fourier transform of an exponential curve

Progam code

clc;

clear all;

close all;

t=0:0.001:1;

a=input('Please enter the multiplication factor of an exponential curve ');

cwtstruct = cwtft(exp(a*t),'plot');

Example of output

Please enter the multiplication factor of an exponential curve  5

_____________________________

MATLAB program to find the Fourier Transform of cosine wave

Program code

clc;

clear all;

close all;

t=0:0.001:1;

cwtstruct = cwtft((cos(2*3.14*1000*t)),'plot');

Output

______________________________________

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 find out the Fourier Transform of sinusoidal waveform

Program code

clc;
clear all;
close all;
t=0:0.001:1;
cwtstruct = cwtft((sin(2*3.14*1000*t)),'plot');

Output

________________________

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

_________________________________

Matlab program to generate logarithmic curve

Program code

clc;
close all;
clear all;
a=input('Please enter the multiplication factor');
t=0:0.001:1;
p=log(a*t);
plot(p);
xlabel('Time');
ylabel('Amplitude');
title('LOGARITHMIC CURVE');
grid on;

Example of output

Please enter the multiplication factor 6

_______________________________