circonvt2.m

来自「matlab求两个数组的卷积的程序」· M 代码 · 共 10 行

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))];
fn=ifft(fft(x1,N).*fft(x2,N),N);