DSP - DFT Sectional Convolution

Suppose, the input sequence x(n) of long duration is to be processed with a system having finite duration impulse response by convolving the two sequences. Since, the pnear filtering performed via DFT involves operation on a fixed size data block, the input sequence is spanided into different fixed size data block before processing.

The successive blocks are then processed one at a time and the results are combined to produce the net result.

As the convolution is performed by spaniding the long input sequence into different fixed size sections, it is called sectioned convolution. A long input sequence is segmented to fixed size blocks, prior to FIR filter processing.

Two methods are used to evaluate the discrete convolution −

Overlap-save method

Overlap-add method

Overlap Save Method

Overlap–save is the traditional name for an efficient way to evaluate the discrete convolution between a very long signal x(n) and a finite impulse response (FIR) filter h(n). Given below are the steps of Overlap save method −

Let the length of input data block = N = L+M-1. Therefore, DFT and IDFT length = N. Each data block carries M-1 data points of previous block followed by L new data points to form a data sequence of length N = L+M-1.

First, N-point DFT is computed for each data block.

By appending (L-1) zeros, the impulse response of FIR filter is increased in length and N point DFT is calculated and stored.

Multippcation of two N-point DFTs H(k) and X_m(k) : Y′_m(k) = H(k).X_m(k), where K=0,1,2,…N-1

Then, IDFT[Y′_m((k)] = y′((n) = [y′_m(0), y′_m(1), y′_m(2),.......y′_m(M-1), y′_m(M),.......y′_m(N-1)]

(here, N-1 = L+M-2)

First M-1 points are corrupted due to apasing and hence, they are discarded because the data record is of length N.

The last L points are exactly same as a result of convolution, so

y′_m (n) = y_m(n) where n = M, M+1,….N-1

To avoid apasing, the last M-1 elements of each data record are saved and these points carry forward to the subsequent record and become 1^st M-1 elements.

Result of IDFT, where first M-1 Points are avoided, to nulpfy apasing and remaining L points constitute desired result as that of a pnear convolution.

Overlap Add Method

Given below are the steps to find out the discrete convolution using Overlap method −

Let the input data block size be L. Therefore, the size of DFT and IDFT: N = L+M-1

Each data block is appended with M-1 zeros to the last.

Compute N-point DFT.

Two N-point DFTs are multipped: Y_m(k) = H(k).X_m(k), where k = 0,,1,2,….,N-1

IDFT [Y_m(k)] produces blocks of length N which are not affected by apasing as the size of DFT is N = L+M-1 and increased lengths of the sequences to N-points by appending M-1 zeros to each block.

Last M-1 points of each block must be overlapped and added to first M-1 points of the succeeding block.

(reason: Each data block terminates with M-1 zeros)

Hence, this method is known Overlap-add method. Thus, we get −

y(n) = {y₁(0), y₁(1), y₁(2), ... .., y₁(L-1), y₁(L)+y₂(0), y₁(L+1)+y₂(1), ... ... .., y₁(N-1)+y₂(M-1),y₂(M), ... ... ... ... ... }