循环卷积函数(Matlab)

function fn=circonvt(x1,x2,N)
%circonvt函数实现输入序列x1和x2的循环卷积,fn为输出序列
%N 为循环卷积长度
%实现方法:fn=sum(x1(m)*x2((n-m) mod N))
if (length(x1)>N|length(x2)>N)%判断输入信号的长度
error('N的长度必须大于输入数据的长度');
end
x1=[x1,zeros(1,N-length(x1))];
x2=[x2,zeros(1,N-length(x2))];
m=0:N-1;
x=zeros(N,N);
for n=0:N-1
    x(:,n+1)=x2(mod((n-m),N)+1)';
end;
fn=x1*x; %循环计算卷积