circonvf.m

matlab例子

function y=circonvf(x1,x2,N)
%在x1和x2：（时域）之间的N点循环卷积
%------------------------------
%[y]=circonvf(x1,x2,N)
%y=包含循环卷积的输出序列
%x1=长度N1<=N的输入序列
%x2=长度N2<=N的输入序列
%N=循环缓冲器的大小
%方法y(n)=sum(x1(m)*x2(m)mod N))

%检查x1的长度
if length(x1)>N
   error('N必须>=x1的长度')
end
%检查x2的长度
if length(x2)>N
   error('N必须>=x2的长度')
end
x1=[x1 zeros(1,N-length(x1))];
x2=[x2 zeros(1,N-length(x2))];
X1=dft(x1,N);
X2=dft(x2,N);
X=X1.*X2;
y=idft(X,N);