整理文档很辛苦,赏杯茶钱您下走!
还剩 ... 页未读,继续阅读 >>
免费阅读已结束,点击下载阅读编辑剩下 ... 页
阅读已结束,您可以下载文档离线阅读编辑
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;