ovrlpsav.m

来自「数字信号处理的MATLAB仿真加实验指导」· M 代码 · 共 14 行

function [y]=ovrlpsav(x,h,N)       %重叠保留法作分段卷积
Lx=length(x);M=length(h);           %N为段长
M1=M-1;L=N-M1;
h=[h zeros(1,N-M)];
x=[zeros(1,M1),x,zeros(1,N-1)];    %预置（M-1）个零
K=floor((Lx+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)';           %去掉前（M-1）个样本
y=(Y(:))';              %装成输出