Multidelay block frequency domain adaptive filter explained
The multidelay block frequency domain adaptive filter (MDF) algorithm is a block-based frequency domain implementation of the (normalised) Least mean squares filter (LMS) algorithm.
Introduction
The MDF algorithm is based on the fact that convolutions may be efficiently computed in the frequency domain (thanks to the fast Fourier transform). However, the algorithm differs from the fast LMS algorithm in that block size it uses may be smaller than the filter length. If both are equal, then MDF reduces to the FLMS algorithm.
The advantages of MDF over the (N)LMS algorithm are:
- Lower algorithmic complexity
- Partial de-correlation of the input (which 'may' lead to faster convergence)
Variable definitions
Let
be the length of the processing blocks,
be the number of blocks and
denote the 2Nx2N Fourier transform matrix. The variables are defined as:
}(\ell) = \mathbf\left[\mathbf{0}_{1xN}, e(\ell N),\dots,e(\ell N-N-1) \right]^T
}_k(\ell) = \mathrm \left\
}(\ell) = \left[\underline{\mathbf{x}}_0(\ell), \underline{\mathbf{x}}_1(\ell), \dots, \underline{\mathbf{x}}_{K-1}(\ell) \right]
}(\ell) = \mathbf\left[\mathbf{0}_{1xN}, d(\ell N),\dots,d(\ell N-N-1) \right]^TWith normalisation matrices
and
:
G1=F\begin{bmatrix}
0N x &0N x \\
0N x &IN x \\
\end{bmatrix}FH
}_2 = \mathbf\begin\mathbf_ & \mathbf_ \\\mathbf_ & \mathbf_ \\\end\mathbf^H
G2=\operatorname{diag}\left\{\tilde{G
}_2, \tilde_2, \dots, \tilde_2 \right\}
In practice, when multiplying a column vector
by
, we take the inverse FFT of
, set the first
values in the result to zero and then take the FFT. This is meant to remove the effects of the circular convolution.
Algorithm description
For each block, the MDF algorithm is computed as:
}}(\ell) = \mathbf_1 \underline(\ell) \underline(\ell-1)
}(\ell) = \underline(\ell) - \underline(\ell)
\Phixx(\ell)=\underline{X
}^H(\ell)\underline(\ell)
}}(\ell) = \underline(\ell-1) + \mu\mathbf_2\mathbf_\mathbf^(\ell) \underline^H(\ell) \underline(\ell)
It is worth noting that, while the algorithm is more easily expressed in matrix form, the actual implementation requires no matrix multiplications. For instance the normalisation matrix computation
}^H(\ell)\underline(\ell) reduces to an element-wise vector multiplication because
}(\ell) is block-diagonal. The same goes for other multiplications.
References
- J.-S. Soo and K. Pang, “Multidelay block frequency domain adaptive filter,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 38, no. 2, pp. 373–376, 1990.
- H. Buchner, J. Benesty, W. Kellermann, "An Extended Multidelay Filter: Fast Low-Delay Algorithms for Very High-Order Adaptive Systems". Proc. IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), 2003.
- A free implementation of the MDF algorithm is available in Speex (main source file)
See also
