circonv.m

来自「是一个完成循环卷积的MATLAB程序 实现方法实在时域完成的」· M 代码 · 共 10 行

function y = circonv(x1,x2,N)
%求序列的循环卷积，方法y(n) = sum(x1(m)*x2(mod((n-m),N)+1))
%检查x1,x2的长度的语句略
x1 = [x1,zeros(1,N-length(x1))];
x2 = [x2,zeros(1,N-length(x2))];
m = [0:N-1];
x2m = x2(mod(-m,N)+1);
H = toeplitz(x2m,[0,x2(2:N)]);
y = x1*H;