数字信号处理的比较全面的程序代码

				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(:))';