circonvt.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

function y=circonvt(x1,x2,N)
%[y]=circonvt(x1,x2,N)
%----------------------
% y(n)=sum(x1(m)*x2((n-m)) mod N))
if length(x1)>N
   error('x1<N')
end
if length(x2)>N
   error('x2<N')
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';