circonvt.m

信号分析与处理的一些源代码

function y=circonvt(x1,x2,N)
%To compute y=x1@x2
if length(x1)>N
   error('Length(x1)is not greater than N');
end

if length(x2)>N
   error('Length(x2)is not greater than N');
end

x1=[x1,zeros(1,N-length(x1))];
x2=[x2,zeros(1,N-length(x2))];
m=[0:N-1];
x2=x2(mod(-m,N)+1);
H=zeros(N,N);
for n=1:N;
   H(n,:)=cirshift(x2,n-1,N);
end
y=x1*H';