cir_conv.m

circular convolution code in matlab

function [result] = cir_conv(x1,x2,n)
m = 4; 
p_n = zeros(n*m,1);

p_n(1,1)=1;
p_n(n+1,1) =1 ;
p_n(2*n+1,1) = 1;
p_n(3*n+1,1) = 1;

tmp = conv(conv(x1,x2),p_n);

result = tmp(2*n + 1:3*n,1);