hsolpadd / 快速重叠相加法

 

• 空间

function [ y ] = hsolpadd( x, h, N )
% High-speed Overlap-Add method of block convolutions using FFT
% -------------------------------------------------------------------------------
% [ y ] = hsolpadd( x, h, N )
% y = output sequence
% x = input real sequence
% h = impulse response
% N = block length (must be a power of two)
%

Lenx = length(x); L = length(h);

% 块长度N要求
if N < L
    error( 'N must be >= length(h)' );
end

%将x分为K块
K = ceil ( Lenx / N ); % 每块N点,向正无穷取整

% x 补足长度
x = [ x zeros(1, K*N-Lenx) ];

% 循环卷积长度M >= N+L