hsolpsav.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 15 行

function [y]=hsolpsav(x,h,N)
%用FFT实现分段卷积的高速重叠保留法
%y=输出序列；x=输入序列；h=脉冲响应；N=段长（基2）
N=2^(ceil(log10(N)/log10(2)));
Lenx=length(x);M=length(h);
M1=M-1;L=N-M1;
h=fft(h,N);
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=fft(x(k*L+1:k*L+N));
    Y(k+1,:)=real(ifft(xk.*h));
end
Y=Y(:,M:N)';y=(Y(:))';