dupconv.m

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

function yn = dupconv( xn,hn )
%DUPCONV Summary of this function goes here
%   duplicate conv of xn and hn ,return yn
M=length(hn);
xn=[xn zeros(1,M-1)];
N=length(xn);
xn1=[zeros(1,M-1),xn(1:floor(N/4))];
xn2=xn((floor(N/4)-M+2):floor(N/2));
xn3=xn((floor(N/2)-M+2):floor(N*3/4));
xn4=xn((floor(N*3/4)-M+2):N);
n1=length(xn1);
n2=length(xn2);
n3=length(xn3);
n4=length(xn4);
yn1=timecirconv(xn1,hn,n1);
yn2=timecirconv(xn2,hn,n2);
yn3=timecirconv(xn3,hn,n3);
yn4=timecirconv(xn4,hn,n4);
yn1=yn1(M:n1);
yn2=yn2(M:n2);
yn3=yn3(M:n3);
yn4=yn4(M:n4);
yn=[yn1 yn2 yn3 yn4];
end