Untitled

clc; clear all; close all;

%% Linear Convolution
x=[1 2 3 4 5];
h=[-2 -1 -3 -5];

N=length(x)+length(h)-1;

x1=[x zeros(1,N-length(x))];
h1=[h zeros(1,N-length(h))];


for i=1:N
    for j=1:N
        H(i,j)=h1(mod((i-j),N)+1);
    end
end

disp ('Linear Convolution Output')
y=(H*x1')'
y1=conv(x,h)

%% Circular Convolution
x=[1 2 3 4 5];
h=[-2 -1 -3 -5];

N=max(length(x),length(h));

x1=[x zeros(1,N-length(x))];
h1=[h zeros(1,N-length(h))];


for i=1:N
    for j=1:N
        H1(i,j)=h1(mod((i-j),N)+1);
    end
end

disp('Circular Convolution Output')
y=(H1*x1')'
y=cconv(x,h,N)

%% Cross Correlation
x=[1 2 3 4 5];
h=[-2 -1 -3 -5];
h2=fliplr(h);

N=length(x)+length(h)-1;

x1=[x zeros(1,N-length(x))];
h1=[h2 zeros(1,N-length(h))];


for i=1:N
    for j=1:N
        H2(i,j)=h1(mod((i-j),N)+1);
    end
end

disp ('Cross Correlation Output')
y=(H2*x1')'
y=xcorr(x,h)

%% Auto Correlation

x=[1 2 3 4 5];
h=fliplr(x);

N=length(x)+length(h)-1;

x1=[x zeros(1,N-length(x))];
h1=[h zeros(1,N-length(h))];


for i=1:N
    for j=1:N
        H3(i,j)=h1(mod((i-j),N)+1);
    end
end

disp ('Auto Correlation Output')
y=(H3*x1')'
y=xcorr(x,x)