Research Article | Open Access
An Alternative Approach to the Balanced Model Truncation Algorithm for Acoustic Minimum-Phase Inverse Filters Order Reduction
We propose an alternative approach to the Balanced Model Truncation method (standard method). This approach reduces substantially the order of minimum-phase inverse filters for equalizing room acoustics. This method is based on a property of the filter z transform function, which modifies the corresponding FIR coefficients before the application of the standard technique to the modified FIR coefficients filter version. In the standard technique, the Hankel singular values plot is the chief guide for a user for the selection of a reduced filter order. Results for minimum-phase inverse filter corresponding to partial equalization of measured acoustic impulse response show the superiority of the proposed method over the standard technique, in terms of reduced filters order selection.
Although partial equalization in room acoustics allows shorter length FIR digital inverse filters than that required in case of complete equalization (perfect dereverberation) [1–3], these are always of very high order (significantly 1000 coefficients or more for an audio-conferencing room), which makes their DSP implementation more complex and inefficient. Such FIR filters have a minimum-phase character, which makes them amenable for possible order reduction when converting them to IIR filters . One very appealing method of FIR to IIR filter conversion is the Balanced Model Truncation technique presented in , which is called here the standard method. In this method, the only allowed way for the selection of a possible reduced filter order depends on the visual Hankel singular values plot, which is the chief guide for a user. The smallest they are the smallest is the filter order. In this paper, we propose an alternative approach to the standard method, which is based on a property of the FIR filter transform function, where corresponding coefficients are modified before the application of the standard method. The goal is to make corresponding Hankel singular values much smaller to substantially reduce the filter order. The standard technique which converts a given single-input/single-output FIR filter to corresponding IIR approximations can be summarized as follows: convert the corresponding transfer function into a matrices system (state-space form), create a Hankel matrix by calculating the Markov parameters, calculate the Hankel singular values, plot and inspect the Hankel singular values and decide upon the reduced system order (equivalently the rank of the Hankel matrix), and finally convert the reduced order state-space form into a transfer function in order to compare the original and the reduced order system against various criteria, and a new order of approximation is attempted if necessary. The paper is organized as follows: Section 2 presents the FIR to IIR filter conversion technique based on Balanced Model Truncation (standard method). In Section 3, an alternative method based on a well-known property of the transform of a digital filter is proposed. Results and comments are given in Section 4, and finally Section 5 concludes the paper.
2. Standard Method
2.1. FIR to IIR Filter Conversion Concept
Let be an coefficients () FIR digital filter where its transfer function in domain is written as follows: This system with the order can be represented in discrete time domain by the following difference equations : where and represent, respectively, the input signal and the output signal of the system and represents the delta function. Notice that is representable by the , , and matrices alone.
Now applying the transform to the system (3), The inverse transform of gives the following infinite impulse response sequence (IIR): represents the Markov parameters for the Hankel matrix which can be defined as follows: Now for the (, , ) system having a finite impulse response (FIR), the rows and columns of zeros are omitted and the following finite Hankel matrix is used: Given such real symmetric matrix (8), the following properties are well known in the literature. (a) can be factorized as follows: where, with being a diagonal matrix containing the real eigenvalues of, an orthogonal matrix containing the corresponding eigenvectors, and a unit matrix, (b)if the matrix contains the eigenvalues of , then contains the eigenvalues of ,(c)the rank of is equal to ,(d)the singular values of are the positive square roots of the eigenvalues of ; then the absolute values of its eigenvalues are given by (10)
2.2. Hankel Singular Values and System Order Reduction
The primary principle used for a system order reduction lies in the elimination of a subsystem associated with the smallest Hankel singular values . Since the rank of is equal to , so the same degree (order) of which represents the system (, , ) then, reducing the rank to an order by applying a partition to the system according to the following partitioning form (11), results in a truncated Hankel matrix with a rank : where with the associated rectangular matrix corresponding to the truncated system () and with the associated rectangular matrix corresponding to the rejected system (), which is associated with the smallest Hankel singular values.
The truncated Hankel matrix with a rank is equivalent to a truncated system (, , ) whose transfer function is given by where with being an extraction of the matrix ’s row from to and its columns from to .
Now having , the transfer function of the truncated system (, , , ), with , can be written as follows: and the corresponding implementable IIR filter coefficients , , , determined after the conversion from the system state-space representation into the transfer function form in domain are given by  with being the reduced filter order.
The original system whose transfer function is and the reduced order system whose transfer function is can then be compared against various criteria, particularly the spectrum magnitude and phase. The plot of Hankel singular values plays the key role in the selection of the filter order reduction.
The standard algorithm can be implemented as follows .
Given coefficients of an FIR filter(1)create the Hankel matrix , (8),(2)decompose the matrix to obtain and , (9),(3)display the Hankel singular values, (10), and choose a desired order of approximation,(4)calculate the matrices , , and , (13),(5)convert the system state-space representation (, , , ), where , into the transfer function form (15) to compare it with that of the original FIR filter (1).
3. Alternative Method
In this method, we use a well-known property of the transform of a filter, being an FIR or IIR, in order to modify the original filter coefficients, that is [6, 7], with being the sequence of the inverse transform of .
The transfer function of an FIR filter can be written as follows: with or
Alternatively, the standard algorithm stated in Section 2 can be applied to the transfer function , ((18), (19)), which represents the transform function of the modified coefficients of the original FIR filter, equivalently the average of the delayed original sequence , , in domain. The corresponding reduced order transfer function, which is the output of the algorithm, will be , with being the reduced order which will be selected from the Hankel singular values plot. The reduced order system whose transfer function is is given by In time domain, this is equivalent to with being the inverse transform of .
The transfer function of the corresponding implementable IIR filter is now given by and the implementation complexity is coefficients, because of the derivative function.
Both techniques were applied to a partial equalization FIR digital filter for an acoustic impulse response, representing the context of an audio-conferencing room (Figure 1) . The corresponding minimum-phase inverse impulse response (FIR prototype) was obtained partially by the Homomorphic method for an optimum . This particular choice of is for avoiding the undesirable effects of the circular deconvolution of the inverse impulse response with the original room impulse response, for some predefined threshold level of magnitude equalization, and which provides some desired sound quality . Figure 6(a) shows its -plane plot. The plot of Hankel singular values of both methods depicted in Figure 2 shows clearly that the alternative method allows the choice of an IIR filter order which is much more reduced than that allowed when using the standard method. When comparing the spectrum magnitude and phase obtained by both methods for the same reduced IIR filter order (with a selected order and a root mean square error of some 0.04) to that of the original FIR prototype, we can see in Figure 3 a better spectral reproduction obtained by the alternative method than that obtained by the standard method, that is, the smoothed form against the affected form (peaks), which is the result of the truncating effect caused by the selection of an underestimated IIR filter order, as displayed by the corresponding Hankel singular values plot in Figure 2.
Both methods were also compared to the direct average of the delayed original FIR sequence (, ) in domain, making sense of noise reduction, that is for the full order. We can see in Figure 3 that despite the fact that the spectral average has completely removed the noise from the original FIR spectrum, which resulted in better smoothed spectral form, the proposed alternative method shows its superiority in tracking nearly closely the form of the original spectrum of the FIR prototype, for both magnitude and phase.
To overcome the complexity implementation problem when using the alternative method, ( coefficients are needed rather than coefficients in the standard method); the standard method has been applied to the transformed version of the reduced order IIR impulse response to an FIR filter (infinite impulse response truncated to a finite number of samples) . The Hankel singular values displayed in Figure 4 allow choosing a new reduced IIR filter order being the same as that delivered by the alternative method, that is, . In Figure 5, we can see a perfect spectral reproduction obtained by the new IIR filter for the reduced order , as compared to that of the transformed IIR to FIR prototype, but with less implementation complexity, that is, coefficients rather than coefficients for the alternative method. This means that the relation (15) is used for DSP implementation rather than (22).
The corresponding magnitude and phase error are both less than −70 dB in the entire frequency band (Figures 5(c) and 5(d)). Figure 6(b) shows the -plane plot for this approximation IIR filter, with only 6 zeros and 6 poles (13 IIR coefficients) rather than more than 1000 zeros (1000 FIR coefficients) for the original FIR prototype (Figure 6(a)).
An alternative approach to the Balanced Model Truncation algorithm (standard method) to substantially reduce minimum-phase inverse filters order for room acoustics equalization is proposed. This method is based on a property of the FIR filter transform function, which modifies the original coefficients of a given FIR sequence, by taking the average of its corresponding transfer function before the conversion to its IIR coefficients approximations. Results obtained by the alternative method show better specifications reproduction of the FIR prototype, with a very low IIR filter order, as compared to that obtained by the standard method.
The authors would like to thank the reviewers for their constructive comments.
- B. D. Radlovic and R. A. Kennedy, “Nonminimum-phase equalization and its subjective importance in room acoustics,” IEEE Transactions on Speech and Audio Processing, vol. 8, no. 6, pp. 728–737, 2000.
- B. D. Radlovic and R. A. Kennedy, “Iterative cepstrum-based approach for speech dereverberation,” in Proceedings of the ISSPA, Brisbane, Australia, August 1999.
- M. Ahfir, I. Kale, and B. Daoud, “Room equalization based on iterative simple complex smoothing of acoustic impulse responses,” in Proceedings of the IEEE ICASSP, vol. 5, pp. 109–112, Toulouse, France, May 2006.
- B. Beliczynski, I. Kale, and G. D. Cain, “Approximation of FIR by IIR digital filters: an algorithm based on balanced model reduction,” IEEE Transactions on Signal Processing, vol. 4, no. 3, pp. 532–542, 1992.
- V. N. Faddeeva, Computational Method of Linear Algebra, Dover, New York, NY, USA, 1959.
- M. Bellanger, Traitement Numerique du Signal, Dunod, Paris, France, 1998.
- M. Kunt, Traitement Numerique des Signaux, Dunod, Paris, France, 1981.
- J. P. Jullien, A. Gilloire, and A. Saliou, “Mesure des réponses impulsionnelles en acoustique,” Note technique NT/LAA/TSS/181, 1984.
- I. Kale, G. D. Cain, and R. C. S. Morling, “Minimum-phase filter design from linear-phase startpoint via balanced model truncation,” Electronics Letters, vol. 31, no. 20, pp. 1728–1729, 1995.
Copyright © 2011 Maamar Ahfir et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.