circonvt.m

计DFT

function y=circonvt(x1,x2,N)
%在x1和x2：（时域）之间的N点循环卷积
%--------------------------------
%[y]=circonvt(x1,x2,N)
%y=包括循环卷积的输出序列
%x1=长度N1<=N的输入序列
%x2=长度N2<=N的输入序列
%N=循环缓冲器的大小
%方法y(n)=sum(x1(m)*x2((n-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))];
m=[0:1:N-1];
x2=x2(mod(-m,N)+1);
H=zeros(N,N);
for n=1:1:N
H(n,:)=cirshift(x2,n-1,N);
end
y=x1*H';