Matlab Program for Circular Convolution

Matlab Program-Circular Convolution:                             

clc

close all

clear all

n=0:3;

xn=[1,2,3,4];

subplot(2,2,1);

stem(n,xn);

grid on

N=4;

for k=0:1:N-1

for n=0:1:N-1

p=exp(-j*pi*2*n*k/N);

x1(k+1,n+1)=p;

end

end

xk=xn*x1;

title('x(n)sequence');

xlabel('n');

ylabel('Amplitude');

disp(xn);

hn=[1 2 1 2];

N=4;

for k=0:1:N-1

for n=0:1:N-1

q=exp(-j*pi*2*n*k/N);

h1(k+1,n+1)=q;

end

end

hk=hn*h1;

n=0:3;

h=[1 2 1 2];

subplot(2,2,2);

stem(n,h);

grid on

title('h(n)sequence');

xlabel('n');

ylabel('Amplitude');

disp(hk);

yk=xk.*hk;

N=4;

for k=0:1:N-1

for n=0:1:N-1

t=exp(j*pi*2*n*k/N);

y1(k+1,n+1)=t/N;

end

end

yn=yk*y1;

s=0:3;

subplot(2,2,3);

stem(s,yn);

grid on

title('y(n)sequence');

xlabel('n');

ylabel('Amplitude');

disp(yk);

OUTPUT: