ovrlpsav.m

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

function [y]=ovrlpsav(x,h,N)
%[y]=ovrlpsav(x,h,N)
%-------------------
%
Lenx=length(x);
M=length(h);
M1=M-1;L=N-M1;
h=[h,zeros(1,N-M)];

x=[zeros(1,M1),x,zeros(1,N-1)];
K=floor((Lenx+M1-1)/L);
Y=zeros(K+1,N);
for k=0:K
   xk=x(k*L+1:k*L+N);
   Y(k+1,:)=circonvt(xk,h,N);
end
Y=Y(:,M:N)';
y=(Y(:))';