timecirconv.m

数字信号处理基于matlab 文件内容 /conv(卷积部分) dupconv 重叠保留法 freqcirconv 频域圆周卷积 freqconv 频域卷积 timecirconv

function yn = timecirconv( xn,hn,N )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
len1=length(xn);
len2=length(hn);
if len1>N | len2>N
    error('N must biger than length of xn and hn')
end
xn=[xn zeros(1,N-len1)];
hn=[hn zeros(1,N-len2)];
m=0:N-1;
x=zeros(N,N);
for n=0:N-1
    x(:,n+1)=hn(mod((n-m),N)+1);
end
yn=xn*x;