相关代码 |
|
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(:))'; %装成输出
相关资源 |
|