循环卷积功能函数MATLAB仿真程序代码

functionfn=circonvt(x1,x2,N)if(length(x1)N||length(x2)N)error('N的长度必须大于输入数据的长度');endx1=[x1,zeros(1,N-length(x1))];x2=[x2,zeros(1,N-length(x2))];m=0:N-1;x=zeros(N,N);forn=0:N-1x(:,n+1)=x2(mod((n-m),N)+1)';endfn=x1*x;