circonvt.m

来自「卷积实验代码 实现多种功能 包括线性卷积和圆周卷积的动态演示」· M 代码 · 共 24 行

function y=circonvt(x1,x2,N)
%y(n)=x1(n)(N)x2(n)
%序列的循环卷积
if length(x1)>N
    error('N must be >= the x1')
end

if length(x2)>N
    error('N must be >= the x2')
end

x1=[x1 zeros(1,N-length(x1))];
x2=[x2 zeros(1,N-length(x2))];

m=[0:1:N-1];

x2=x2(mod(-m,N)+1);
H=zeros(N,N);
for n=1:1:N
    H(n,:)=cirshftt(x2,n-1,N);
end

y=x1*H';