Scribd
Upload
7 views4 pages

Circular Convolution-1

Circular convolution using MATLAB

Uploaded by

BACHU AKSHITHA
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
7 views4 pages

Circular Convolution-1

Circular convolution using MATLAB

Uploaded by

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

Circular convolution (without built in

function)
% Circular Convolution in MATLAB

clc;

clear;

% Input two sequences

x = input('Enter the first sequence x[n]: ');

h = input('Enter the second sequence h[n]: ');

% Ensure the sequences have the same length by zero-padding the shorter one

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

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

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

% Initialize the output sequence

y = zeros(1, N);

% Perform circular convolution

for n = 1:N

for k = 1:N

index = mod(n - k, N); % Circular index

if index < 0

index = index + N; % Handle negative modulo

end

y(n) = y(n) + x(k) * h(index + 1); % Add 1 for MATLAB 1-based indexing

end

end
% Display the result

disp('The result of the circular convolution y[n] is:');

disp(y);

% Plot x[n]

n=0:N-1;

subplot(3, 1, 1);

stem(n, x, 'filled');

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

xlabel('n');

ylabel('x[n]');

grid on;

% Plot h[n]

subplot(3, 1, 2);

stem(n, h, 'filled');

title('Input Sequence h[n]');

xlabel('n'); ylabel('h[n]');

grid on;

% Plot y[n]

subplot(3, 1, 3);

stem(n, y, 'filled');

title('Circular Convolution Result y[n]');

xlabel('n');

ylabel('y[n]');

grid on;
Circular convolution (with built in function)

% Circular Convolution using cconv

clc;

clear;

% Input two sequences

x = input('Enter the first sequence x[n]: ');

h = input('Enter the second sequence h[n]: ');

% Length of the circular convolution (optional, default is max(length(x), length(h)))

N = input('Enter the length of circular convolution (optional, press Enter for default): ');

if isempty(N) % If no length is specified, use default length

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

end

% Perform circular convolution using cconv

y = cconv(x, h, N);

% Display the result

disp('The result of the circular convolution y[n] is:');

disp(y);

% Plot the input sequences and the result figure;

% Plot x[n]

n=0:length(x)-1;

subplot(3, 1, 1);

stem(n, x, 'filled');

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

xlabel('n'); ylabel('x[n]');
grid on;

% Plot h[n]

n= 0:length(h)-1;

subplot(3, 1, 2);

stem(n, h, 'filled');

title('Input Sequence h[n]');

xlabel('n'); ylabel('h[n]');

grid on;

% Plot y[n]

n= 0:length(y)-1;

subplot(3, 1, 3);

stem(n, y, 'filled');

title('Circular Convolution Result y[n]');

xlabel('n');

ylabel('y[n]');

grid on;

You might also like