Issue 
Acta Acust.
Volume 7, 2023



Article Number  2  
Number of page(s)  10  
Section  Active Control  
DOI  https://doi.org/10.1051/aacus/2022058  
Published online  09 January 2023 
Scientific Article
Robust direct acoustic impedance control using two microphones for mixed feedforwardfeedback controller
Signal Processing Laboratory 2, EPFL, 1015 Lausanne, Switzerland
^{*} Corresponding author: maxime.volery@epfl.ch
Received:
14
July
2022
Accepted:
24
November
2022
This paper presents an acoustic impedance control architecture for an electroacoustic absorber combining both feedforward and feedback microphonebased strategies on a currentdriven loudspeaker. Feedforward systems enable good performance for direct impedance control. However, inaccuracies in the required actuator model can lead to a loss of passivity, which can cause unstable behaviour. The feedback contribution allows the absorber to better handle model errors and still achieve an accurate impedance, preserving passivity. Numerical and experimental studies were conducted to compare this new architecture against a stateoftheart feedforward control method.
Key words: Active sound absorption / Electrodynamic loudspeaker / Feedback control / Feedforward control / Model uncertainty / Passivity / Pressure control
© The Author(s), Published by EDP Sciences, 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Electroacoustic absorption consists in controlling the acoustic impedance presented by an electroacoustic actuator, typically an electrodynamic loudspeaker [1]. The control of this impedance can be done passively, by loading the voice coil of the loudspeaker with an appropriate electrical impedance [2, 3], or actively, using one or more sensors controlling the voltage or current applied to the actuator. Active electroacoustic absorbers have a wide range of applications, spanning from room acoustics [4] to aircraft engine noise reduction [5] thanks to their advantage of being tuneable, broadband and of subwavelength dimensions. Most of the stateofthe art active absorber designs are either not tuneable, such as in the hybrid passive/active absorption concept [6] or require both a pressure and velocity sensor for a feedback implementation. The sensing of the velocity can, for instance, be achieved using an accelerometer placed on the loudspeaker cone [7] (not acceptable for small loudspeakers), two closely placed microphones [8] (not practical because upstream from the impedance plane) or a Wheatstone bridge [9] (requires fine resistors and inductance tuning).
However, should the model of the actuator be known, a feedforward architecture [4] can be used where only a single sensor is needed. Also, thanks to the model inversion, direct impedance control can be achieved accurately, whereas other methods only approach the target impedance. Nevertheless, due to some inevitable inaccuracies in the estimation of the model parameters and the delay of the numerical controller, a mismatch between the target impedance and the achieved one will eventually occur. This mismatch can cause a loss of acoustic passivity of the absorber, meaning that it is injecting energy into the acoustic environment instead of absorbing it. Such behaviour is unwelcome, even if it occurs outside of the frequency band of interest, because it can result in an unstable positive acoustic feedback. In other words, if at a given frequency, the absorber injects more energy than the acoustic environment dissipates, energy will buildup, leading to an instability [10].
Combining both a feedforward and a feedback loop can help reduce the inaccuracies while keeping the same performances, enabling a better fit with the analytical target impedance. The membrane velocity estimation needed for the feedback implementation can be obtained via a microphone placed inside the cavity of the loudspeaker [11, 12]. Indeed, for wavelengths smaller than the cabinet dimensions, the acoustic pressure behind the actuator is proportional to its membrane displacement and can be used to control it. With this configuration, the size and complexity of the proposed mixed feedforwardfeedback strategy does not fundamentally change from the former feedforwardonly architecture and can be directly compared.
This paper is organized as follows. In Section 2, a model of the electrodynamic loudspeaker is introduced before the description of the twoinput control architecture. Section 3 presents a MonteCarlo analysis of the sensitivity of the achieved absorption to the model estimation errors. Experimental validation of the proposed architecture is given in Section 4 for three different control configurations, and Section 5 provides conclusion and opens some future perspectives for the presented concept.
2 Robust electroacoustic absorber design
2.1 Model of the electrodynamic loudspeaker
An electrodynamic loudspeaker can be modelled as a massspringdamper system, of mass M_{ms}, mechanical compliance C_{ms} and mechanical resistance R_{ms} It is thus a second order resonator [13]. Three forces act on its membrane: the pressure in front of the membrane p_{f}, the pressure behind the membrane p_{b} and the Lorentz force due to the current i flowing in the voice coil. When mounted on an enclosure, the contribution from the rear pressure can be modelled as a specific compliance C_{sb} for wavelengths smaller than the cabinet dimensions. This compliance is the ratio between membrane displacement and the pressure in the cavity, and is linked to the volume of the cavity V_{b} as follows:
(1)where ρ_{0} is the mass density of air, c_{0} the speed of sound in the air, and S_{d} the effective piston area of the loudspeaker. The membrane motion is described by Newton’s second law of motion
(2)where v is the membrane inwards velocity, Bl the coil force factor, and C_{mc} the combined mechanical compliance of the loudspeaker and the cabinet. The pressure in the cabinet p_{b} is directly proportional to the membrane displacement
In the Laplace domain, with Laplace variable s, equations (2) and (3) are written
(7)is the specific impedance of the loudspeaker, its specific resistance, its natural resonance angular frequency and its (passive) quality factor. From the representation of the impedance of equation (7), it is straightforward to notice that the passive loudspeaker (i = 0) mounted on a cabinet is indeed a second order resonator.
Because an accurate model of the electrical impedance of the loudspeaker is complex to develop and to estimate [14, 15], and that the electrical force applied on the membrane is directly proportional to the current flowing in the coil, as shown in equation (4), it is interesting to drive the loudspeaker using a current source rather than a voltage source, as has been done in [4]. In the following, the loudspeaker is driven in current. An implementation of such a current source is given in Appendix A.
2.2 Formulation of the TwoInput SingleOutput controller
Direct impedance control allows to reach a desired target impedance on the membrane of the loudspeaker instead of the passive one . A feedforwardcontroller [4] measures the pressure in front of the membrane and relies on the model of the actuator to find the current to inject in the voice coil to get the appropriate membrane velocity such that the desired target impedance is met. It is therefore capable of reaching a wide range of target impedances. However, this also implies that an accurate model of the loudspeaker must be given to the controller, and that any inaccuracy in this model can have an important impact on the obtained results (i.e., the achieved impedance will deviate from the target one). Adding a feedback loop along with the feedforward architecture can help reduce this problem. To implement feedback on top of the feedforward architecture, a measure of the velocity of the membrane is needed in addition to the pressure in front of it. This can be achieved by sensing the pressure in the cavity closing the rear face of the actuator because the pressure in it is proportional to the displacement of the membrane at low frequencies, as shown in equation (5).
It appears now that the controller has two inputs: the pressure in front of the membrane p_{f} and the pressure behind it p_{b} and has a single output: the current i injected in the moving coil of the loudspeaker. This output current can therefore be expressed as
(8)where both H_{1} and H_{2} are linear timeinvariant systems. An illustration of such a controller is shown in Figure 1, and its detailed block diagram in Figure 2. In the latter, it is clearly visible that H_{1}(s) is the feedforward part of the controller and H_{2}(s) the feedback part.
Figure 1 Controlled absorber. The twoinput controller is depicted on the right in the dashed rectangle. 
Figure 2 Block diagram of the mixed feedforwardfeedback controlled absorber. 
In order to achieve a target impedance Z_{st} (s), it follows from equations (4), (5) and (8) that H_{1} and H_{2} must satisfy the relation
There is an infinite number of realizations that satisfy equation (9), but feedback from the membrane velocity is desired. This feedback in velocity G(s) is the combination of the controller H_{2}, the compliance of the enclosure and the force factor. And because the modelling of the box as a constant compliance is only valid for wavelengths smaller than the dimension of the box, G(s) should have a lowpass behaviour. A first order lowpass filter is chosen for G(s) such that the controller is of the smallest degree possible:
(10)where k_{g} ≥ 0 is a dimensionless tuneable feedback gain and ω_{g} is the cutoff angular frequency of the lowpass filter G(s). The two control transfer functions are thus
In equations (11) and (12), it can be observed that the controller is proper, and that by setting G = 0, only H_{1}(s) is left, and is equal to the stateofthe art feedforward controller from [4] without any feedback. Furthermore, equations (11) and (12) can also be interpreted as the superposition of the pure feedforward implementation and a pure feedback implementation where the error between target velocity and achieved velocity is fed as a current to the loudspeaker with feedback gain G(s), as in [16].
However, not any arbitrary impedance can be achieved: to avoid divergence of the control transfer function H_{1}(s) for low and high frequencies, the asymptotes of the target impedance should behave as a compliance for low frequency, and a mass for high frequencies, as it is the case for the passive impedance. In this article, the considered target impedance is a multidegreeoffreedom resonator, which is the result of N second order resonators connected in parallel, as used in [17]
(13)where R_{st,n}, ω_{t,n} and Q_{t,n} are respectively the specific resistance, the resonance angular frequency and the quality factor of the nth resonator. Different realizations of the target impedance could also be considered, but the following derivation will consider the form of equation (13) without loss of generality.
2.3 Proof of stability
A pole analysis of the feedback loop created by H_{2}(s) is required to show the stability properties of the absorber. Each transfer functions H_{1}(s) and H_{2}(s) are individually (open loop) proper and stable. There is one feedforward loop, which is stable if its components are stable, and a feedback loop which is stable if the real part of all its poles is negative. These poles are the solutions of
(14)where T(s) is the closed loop transfer function between (1 − FH_{1})p_{f} and v. This is equivalent to solving
(18)and it is interesting to notice that equation (15) does not depend on the target impedance. The closed loop T(s) is stable if and only if the Hurwitz matrix
(19)corresponding to the polynomial of equation (15) has all its three leading principal minors which are positive [18]:
This means that k_{g} must satisfy
(23)which is always true for nonnegative values of k_{g}.
2.4 Sensitivity to parameter variations
To analyse the robustness of the proposed method to parameter estimation accuracy, the sensitivity functions of the achieved impedance are calculated. When the estimated values , and of the parameters , and respectively are used in the controller transfer functions from equations (11) and (12), the achieved impedance is
The sensitivity function of this achieved impedance with respect to a parameter x is defined as the ratio between the percentage of change in the achieved impedance Z_{sa} to the percentage of change in the parameter x [19]:
(28)for parameters , and , respectively. The limit when G(s) → ∞ of , and are respectively 0, 0 and 1. It can therefore be concluded that any variation in the estimation and will be less significant when the magnitude of G(s) is larger. This is however not true for , for which the error on the achieved impedance becomes proportional to the error in when the magnitude of G(s) is large.
3 Numerical sensitivity analysis
In this section, a numerical sensitivity analysis is presented for three different control targets: a singledegreeoffreedom resonator whose resonance is shifted with respect to the passive one, a broadband absorption centred at the passive resonance and a twodegreeoffreedom impedance with two distinct shifted resonances. The target impedances and the control parameters are defined according to equation (13) and are reported for each case in Table 1.
Target impedances and control parameters for the three considered configurations.
The numerical sensitivity analysis consists in evaluating the achieved normal incidence absorption coefficient α_{a} 10^{5} times, with random Gaussian deviations of 5% on the estimated parameters , , , and . This absorption coefficient is defined as the ratio between absorbed and incident power. It lies between 0 and 1 for acoustically passive systems, whereas it is smaller than one if the system is acoustically active (for which energy is injected in the acoustic domain instead of being absorbed). It is calculated from the achieved impedance Z_{sa}(s) as
(29)where Z_{sa}, the achieved impedance is evaluated according to equation (24). At every simulated frequency, the values of the first and the third quartiles of the absorption coefficient are reported in Figures 3–5 for each considered target. In these figures, it is observable that the absorption coefficient with only feedforward deviates further away from the target than with the mixed feedforwardfeedback control. It can even reach negative values around the passive resonance of the actuator. With feedback however, it is much better controlled around this resonance, but at the price of lower accuracy for other frequencies.
Figure 3 First and third quartiles of the achieved absorption for the singledegreeoffreedom absorber with 10^{5} random relative errors of 5% standard deviation on the five estimated parameters. 
Figure 4 First and third quartiles of the achieved absorption for the broadband absorber with 10^{5} random relative errors of 5% standard deviation on the five estimated parameters. 
Figure 5 First and third quartiles of the achieved absorption for the twodegreeoffreedom absorber with 10^{5} random relative errors of 5% standard deviation on the five estimated parameters. 
Although the feedback does not bring much improvement for the broadband absorption shown in Figure 4, it does for the two other cases. In an Ultra High Bypass Ratio aircraft engine application, the sound to absorb is typically tonal, and an absorber with multiples rays of absorptions would be convenient [5, 20]. Also, in this application, the optimal impedance would not be ρ_{0}c_{0} but rather consists of a given resistive part and a reactive part, as explained in [21], for which this new architecture can bring interesting improvements.
4 Experimental results
4.1 Experimental setup
The measurement setup used to experimentally assess this new control architecture is shown in Figure 6, and schematised in Figure 7. The two microphones used to control the electroacoustic absorber are connected to the fieldprogrammable gate array (FPGA) controller through a signal conditioner. The digital filter running on the FPGA is the bilinear transform of equations (11) and (12) with a sampling frequency of 50 kHz. For better numerical stability, the digital filter is realized as a cascade of secondorder sections [22]. The output voltage of the controller is converted into a current by a homemade voltagecontrolled current source whose schematic is described in Appendix A. A short study on the impact of the position of the rear microphone is available in Appendix B.
Figure 6 Experimental setup used to measure the impedance presented by the absorber. 1) Electroacoustic resonator; 2) Measurement microphones; 3) Frequency analyser; 4) Power amplifier; 5) Sound source; 6) IEPE signal conditioner; 7) FPGA controller; 8) current pump. 
Figure 7 Schematic of the experimental setup used to measure the impedance presented by the absorber. 
Experimental setup used to measure the impedance presented by the absorber. 1) Electroacoustic resonator; 2) measurement microphones; 3) frequency analyser; 4) power amplifier; 5) sound source; 6) IEPE signal conditioner; 7) FPGA controller; 8) current pump.
The achieved impedance presented by the absorber is measured using a Kundt’s tube after ISO 105342 [23]. A multichannel frequency analyser feeds white noise to the amplified external source during 60 (resulting in a sound pressure level up to 105 dB at the absorber position) while measuring the signals from the two measurement microphones p_{1} and p_{2}. From the transfer function p_{2}(s)/p_{1}(s) and the waveguide dimensions Δx and x_{1}, the reflection coefficient of the termination of the waveguide, and thus its impedance too, can be recovered [23]. The estimation of the transfer function is done with a linear averaging of 1 length Hann windows overlapping by 66.67%, with a 1 Hz frequency resolution. All the hardware equipment used is listed in Table 1.
4.2 Transducer parameters identification
To implement the filters from equations (11) and (12), five parameters of the electrodynamic loudspeaker are needed: R_{ss}, ω_{0}, Q_{ms}, F and C_{sb}. The estimation of the specific mass M_{ss} = M_{ms}/S_{d}, resistance R_{ss} and stiffness K_{sc} = 1/(S_{d} C_{mc}) are obtained by a polynomial fitting of the measured passive (i = 0) impedance curve:
(30)where + denotes the Moore–Penrose pseudoinverse, ω is a vector containing the measured angular frequencies, Z_{ss}(s) is the measured specific impedance and and are vectors of respectively zeros or ones of the same size as ω. The parameters ω_{0} and Q_{ms} are straightforward to derive from the result of equation (30). Then, can be estimated as presented in [24], using the proportional controller :
(31)where Z_{1}(s) is the specific impedance measured with the constant feedforward controller of gain K_{1} and ω_{n} is the nth element of ω. Finally, the box specific compliance can be found using the proportional controller i = K_{2}p_{b}:
(32)where Z_{2}(s) is the specific impedance measured with the constant feedback controller of gain K_{2}.
All these measured parameters of the electrodynamic absorber are reported in Table 2. The frequency band considered in equations (30), (31) and (32) is from 170 Hz to 250 Hz with steps of 1 Hz. Note that these parameters describe the termination of the Kundt’s tube. To get the loudspeaker parameters, they must be scaled by S_{d}/S_{duct}, where S_{duct} is the cross section of the duct. However, this is not necessary if one is interested in controlling the impedance of the whole termination instead of only the loudspeaker. Indeed, using the cross section S_{duct} instead of S_{d} is equivalent to a scaling of v, and thus a scaling of the impedances and the box compliance. It therefore has no impact on the equations if all the measured impedances as well as the target one are considered with the same crosssection. It is also interesting to notice that the calibration of the two control microphones is not necessary. Indeed, in both equations (31) and (32) the errors in the microphone sensitivities are embedded in the estimation of F and C_{sb}.
Experimental setup equipment list.
4.3 Impedance measurements
The three considered target impedances are described by the parameters from Table 3. To highlight the advantage of the mixed feedforwardfeedback controller, a 5% error was purposely included in the model of the loudspeaker, needed to build the controller transfer functions, such that . In Figures 8–10, the passive, the target and the achieved absorption coefficients with and without the feedback contribution are drawn.
Figure 8 Experimentally obtained absorption coefficients for the singledegreeoffreedom absorber, with = 0.95F. 
Figure 9 Experimentally obtained absorption coefficients for the broadband absorber, with = 0.95F. 
Figure 10 Experimentally obtained absorption coefficients for the twodegreeoffreedom absorber, with = 0.95F. 
Measured ThieleSmall parameters of the Monacor SPX30M loudspeaker mounted on a cabinet.
Like for the numerical study, it is observed that the passive resonant behaviour is still present in the achieved impedances without feedback, reaching in some cases a negative value of absorption and adding a degree of freedom to the achieved impedance. The mixed feedforwardfeedback controller is capable to overcome this issue, does truly behave as the target and is more accurate, especially around the passive resonance of the loudspeaker. Note that the lack of precision at lower frequencies (i.e., lower than 100 Hz) for both controllers is inherent to the Kundt’s tube measurement. Indeed, the termination reflection coefficient Γ(s) is given in [23] as
(33)where H_{12}(s) = p_{2}(s)/p_{1}(s) is the transfer function between the two measurement microphones, k is the wave number and Δx and x_{1} are dimensions visible in Figure 7. When the frequency tends to zero, equation (33) becomes ill conditioned because both H_{12}(s) and e^{±jkΔx} tend to one. Equation (33) is therefore very sensitive to the measurement errors in H_{12} for low frequencies.
5 Conclusions
This article presented a new direct impedance control architecture providing a more accurate and robust control on the actual impedance than previously reported in the literature. The concept of mixed feedforwardfeedback control is based on an already existing feedforward implementation, but to achieve a better accuracy, it is combined with a feedback loop that relies on the sensing of the displacement of the actuator to adjust the driving current. Displacement sensing is done through a microphone placed in the enclosure of the loudspeaker, effective at low frequencies. Even if it is not a noticeable improvement for broadband absorption, as targeted by the feedforward architecture [4], it does significantly improve the passivity, and thus the stability, of a multidegreeoffreedom absorber, as formerly used in aircraft engine noise reduction applications. Additionally, in such an environment, the estimated parameters of the absorber might change significantly with the static pressure, surrounding temperature or humidity. With the feedback contribution, the sensitivity to errors is lowered, and is therefore more adapted to drifting parameters.
This design could be further improved, typically by investigating different relations between H_{1} and H_{2} in equation (9). Also, a more sophisticated model of the relationship between the membrane velocity and the pressure in the cavity could be considered to extend the feedback contribution to higher frequencies or larger loudspeaker enclosures. For this, a more elaborated fitting should be used in equation (32) rather than a constant real value. Furthermore, the mixed feedforwardfeedback control could also be used to linearize actuators at high sound pressure levels, at which their stiffness is no longer linear and typically depends on the membrane position.
Acknowledgments
This project has received funding from the Clean Sky 2 Joint Undertaking under the European Union’s Horizon 2020 research and innovation program under grant Agreement No. 821093. This publication reflects only the authors’ view, and the JU is not responsible for any use that may be made of the information it contains.
Appendix A
A.1 Current source
The voltage controlled current source used to drive the loudspeaker for the experimental measurements is depicted in Figure 11 and is inspired from the application report [25]. The chosen operational amplifier is a TL288CP from Texas Instruments. The output current can be shown to be
(34)When R_{1} = R_{2} and R_{3} = R_{4} + R_{5}, it simplifies to a proportional relation between input voltage and output current, regardless of the load impedance Z_{L}:
(35)With the values from Figure 11, a suitable voltage controlled current source for driving a loudspeaker is obtained:
Figure 11 Voltage controlled current source schematic. R_{1} = R_{2} = 92 kΩ, R_{3} = R_{4} = 1:1 kΩ and R_{5} = 1.2 Ω. 
The current delivered by the operational amplifier is
(37)which is approximately the output current of the current pump since R_{5} and the Z_{L} are both much smaller than R_{3}. For the 2 DOF case from Table 3, the highest current is required when all the incident pressure is concentrated at 100.8 Hz. The maximal output current for the TL288CP is of 80 mA, which is reached when the incident pressure is 117 dB SPL at 100.8 Hz.
Appendix B
B.1 Microphone position in the cavity
For wavelengths much smaller than the dimension of the enclosure of the loudspeaker, the pressure in the cavity is proportional to the displacement of the membrane. However, as the frequency increases, the model of the box is becoming worse, and cavity modes appear. The position of the microphone in the cavity can help mitigate this effect.
Frequencydomain simulations have been conducted using the finite element simulation software COMSOL Multiphysics to find an optimal microphone position. The obtained relationships from the membrane displacement to the pressure at the position of the microphone p_{b}/ξ are reported in Figure 12 for the two geometries shown in Figure 13. In this graph, it is visible that the first cavity mode happens at 2.2 kHz. The response of the microphone at position 1 has the flattest response up to this frequency and is therefore chosen in the experimental absorber prototype. However, to avoid instabilities at high frequencies, some melamine foam was added in the enclosure, which will damp higher frequencies and remove the undesired spikes.
Figure 12 Simulated transfer function between rear microphone pressure and membrane displacement . 
Figure 13 Simulated geometry, with the two microphones positions. Membrane is drawn in a thin line and the magnet is hatched. Units in mm. 
References
 H.F. Olson, E.G. May: Electronic sound absorber. The Journal of the Acoustical Society of America 25 (1953) 1130–1136. https://doi.org/10.1121/1.1907249. [Google Scholar]
 A.J. Fleming, D. Niederberger, S.O.R. Moheimani, M. Morari: Control of resonant acoustic sound fields by electrical shunting of a loudspeaker. IEEE Transactions on Control Systems Technology 15 (2007) 689–703. https://doi.org/10.1109/TCST.2006.890276. [CrossRef] [Google Scholar]
 R. Boulandet, E. Rivet, H. Lissek: Sensorless electroacoustic absorbers through synthesized impedance control for damping lowfrequency modes in cavities. Acta Acustica united with Acustica 102 (2016) 696–704. https://doi.org/10.3813/AAA.918986. [CrossRef] [Google Scholar]
 E. Rivet, S. Karkar, H. Lissek: Broadband lowfrequency electroacoustic absorbers through hybrid sensor/shuntbased impedance control. IEEE Transactions on Control Systems Technology 25 (2017) 63–72. https://doi.org/10.1109/TCST.2016.2547981. [CrossRef] [Google Scholar]
 R. Boulandet, H. Lissek, S. Karkar, M. Collet, G. Matten, M. Ouisse, M. Versaevel: Duct modes damping through an adjustable electroacoustic liner under grazing incidence. Journal of Sound and Vibration 426 (2018) 19–33. https://doi.org/10.1016/j.jsv.2018.04.009. [Google Scholar]
 M.A. Galland, B. Mazeaud, N. Sellen: Hybrid passive/active absorbers for flow ducts. Applied Acoustics 66 (2005) 691–708. https://doi.org/10.1016/j.apacoust.2004.09.007. [CrossRef] [Google Scholar]
 T.J. Cox, P. D’Antonio: Acoustic absorbers and diffusers. CRC Press, London, 2004. https://doi.org/10.4324/9781482266412. [Google Scholar]
 F. OrduñaBustamante, P. Nelson: An adaptive controller for the active absorption of sound. The Journal of the Acoustical Society of America 91 (1992) 2740–2747. https://doi.org/10.1121/1.403779. [Google Scholar]
 X. Meynial, H. Lissek: Active reflectors for room acoustics. Proceedings of the Institute of Acoustics 21 (1999) 99–106. http://infoscience.epfl.ch/record/148493. [Google Scholar]
 E. De Bono, M. Collet, G. Matten, S. Karkar, H. Lissek, M. Ouisse, K. Billon, T. Laurence, M. Volery: Effect of time delay on the impedance control of a pressurebased, currentdriven electroacoustic absorber. Journal of Sound and Vibration 537 (2022) 117201. https://doi.org/10.1016/j.jsv.2022.117201. [Google Scholar]
 X. Guo, R. Fleury, H. Lissek: Improving sound absorption through nonlinear active electroacoustic resonators. Physical Review Applied 13 (2020) 014018. https://doi.org/10.1103/PhysRevApplied.13.014018. [Google Scholar]
 X. Guo, M. Volery, H. Lissek: PIDlike active impedance control for electroacoustic resonators to design tunable singledegreeoffreedom sound absorbers. Journal of Sound and Vibration 525 (2022) 116784. https://doi.org/10.1016/j.jsv.2022.116784. [Google Scholar]
 M. Rossi: Audio. Presses Polytechniques et Universitaires Romandes, Lausanne, 2007. [Google Scholar]
 W.M. Leach Jr: Loudspeaker voicecoil inductance losses: circuit models, parameter estimation, and effect on frequency response. Journal of the Audio Engineering Society 50 (2002) 442–450. [Google Scholar]
 K. Thorborg, A.D. Unruh: Electrical equivalent circuit model for dynamic movingcoil transducers incorporating a semiinductor. Journal of the Audio Engineering Society 56 (2008) 696–709. [Google Scholar]
 M. Volery, H. Lissek: Achieving direct acoustic impedance control with only two microphones, in Forum Acusticum, Lyon, France, December, 2020, 1573–1574. https://doi.org/10.48465/fa.2020.0586. [Google Scholar]
 E. Rivet, S. Karkar, H. Lissek: On the optimisation of multidegreeoffreedom acoustic impedances of lowfrequency electroacoustic absorbers for room modal equalisation. Acta Acustica united with Acustica 103 (2017) 1025–1036. https://doi.org/10.3813/AAA.919132. [CrossRef] [Google Scholar]
 A. Hurwitz: Ueber die bedingungen, unter welchen eine gleichung nur wurzeln mit negativen reellen theilen besitzt. Mathematische Annalen 46 (1895) 273–284. https://doi.org/10.1007/bf01446812. [Google Scholar]
 S. Shinners: Modern control system theory and design. 2nd ed., John Wiley & Sons, Nashville, TN, 1998. [Google Scholar]
 E. Salze, A. Pereira, P. Souchotte, J. Regnard, F. GeaAguilera, M. Gruber: New modular fan rig for advanced aeroacoustic tests – acoustic characterization of the facility, in 25th AIAA/CEAS Aeroacoustics Conference, May 2019 in Delft. https://doi.org/10.2514/6.20192603. [Google Scholar]
 B. Tester: The propagation and attenuation of sound in lined ducts containing uniform or “plug” flow. Journal of Sound and Vibration 28, 2 (1973) 151–203. https://doi.org/10.1016/S0022460X(73)801026. [Google Scholar]
 S.K. Mitra: Digital signal processing: a computerbased approach. McGrawHill, New York, NY, 1998. [Google Scholar]
 ISO 10534–2:1998: Determination of sound absorption coefficient and impedance in impedance tubes – Part 2: Transferfunction method. International Organization for Standardization, Geneva, CH, 1998. [Google Scholar]
 E. De Bono: Electroactive boundary control for noise mitigation: local and advective strategies. PhD Thesis. Ecole Centrale de Lyon, 2021. 2021LYSEC024. [Google Scholar]
 Texas Instruments: A comprehensive study of the Howland current pump. Application Report, January 2008. [Google Scholar]
Cite this article as: Volery M. Guo X. & Lissek H. 2023. Robust direct acoustic impedance control using two microphones for mixed feedforwardfeedback controller. Acta Acustica, 7, 2.
All Tables
Target impedances and control parameters for the three considered configurations.
Measured ThieleSmall parameters of the Monacor SPX30M loudspeaker mounted on a cabinet.
All Figures
Figure 1 Controlled absorber. The twoinput controller is depicted on the right in the dashed rectangle. 

In the text 
Figure 2 Block diagram of the mixed feedforwardfeedback controlled absorber. 

In the text 
Figure 3 First and third quartiles of the achieved absorption for the singledegreeoffreedom absorber with 10^{5} random relative errors of 5% standard deviation on the five estimated parameters. 

In the text 
Figure 4 First and third quartiles of the achieved absorption for the broadband absorber with 10^{5} random relative errors of 5% standard deviation on the five estimated parameters. 

In the text 
Figure 5 First and third quartiles of the achieved absorption for the twodegreeoffreedom absorber with 10^{5} random relative errors of 5% standard deviation on the five estimated parameters. 

In the text 
Figure 6 Experimental setup used to measure the impedance presented by the absorber. 1) Electroacoustic resonator; 2) Measurement microphones; 3) Frequency analyser; 4) Power amplifier; 5) Sound source; 6) IEPE signal conditioner; 7) FPGA controller; 8) current pump. 

In the text 
Figure 7 Schematic of the experimental setup used to measure the impedance presented by the absorber. 

In the text 
Figure 8 Experimentally obtained absorption coefficients for the singledegreeoffreedom absorber, with = 0.95F. 

In the text 
Figure 9 Experimentally obtained absorption coefficients for the broadband absorber, with = 0.95F. 

In the text 
Figure 10 Experimentally obtained absorption coefficients for the twodegreeoffreedom absorber, with = 0.95F. 

In the text 
Figure 11 Voltage controlled current source schematic. R_{1} = R_{2} = 92 kΩ, R_{3} = R_{4} = 1:1 kΩ and R_{5} = 1.2 Ω. 

In the text 
Figure 12 Simulated transfer function between rear microphone pressure and membrane displacement . 

In the text 
Figure 13 Simulated geometry, with the two microphones positions. Membrane is drawn in a thin line and the magnet is hatched. Units in mm. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.