circonv.m

来自「自相关函数在噪声检测,识别和有用信号的提取中的应用」· M 代码 · 共 7 行

function y=circonv(x1,x2,N)
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;