Uncertainty Constraint on Headphone Secondary Path Function for Designing Cascade Biquad Feedback Controller with Improved Noise Reduction Performance

: The uncertainty in the secondary path of active noise control (ANC) headphones affects the waterbed effect and stability of the feedback system. This study focuses on the uncertainty of the secondary path when real users wear headphones and proposes a new uncertainty constraint based on the measured results of the secondary path transfer function under different wearing conditions of a dummy head and limited subjects. This constraint and a cascaded second-order infinite impulse response filter with fixed coefficients are used to formulate a control strategic function, which is optimized using the Improved Grey Wolf Optimizer (IGWO) algorithm to obtain the optimal controller with better noise reduction performance. The proposed method and simulation model are validated based on the experimental test results. The results demonstrate that the safety factor and waterbed suppressing factor contained in the proposed uncertainty constraint ensure more stable noise reduction and effective suppression of the waterbed effect for new subjects without a priori data.


Introduction
In the field of commercial headphones, ANC has become a common application [1].Controllers in this realm are classified into two categories based on the time variability of their filtering coefficients: fixed controllers and adaptive controllers.Fixed controllers contribute to reduced computational demands and extended battery life, thereby improving the standby time of headphones.Recently, Bai et al. [2] found that compared to an adaptive controller, a fixed controller achieved satisfactory noise reduction performance and signal tracking quality, which are highly valuable merits in practical applications.
ANC controllers are mainly divided into feedforward, feedback, and hybrid controllers in terms of their control structures.Notably, the feedback control strategies have been paid particular attention in a long history of research [3].The introduction of feedback mechanisms leads to two primary issues: firstly, according to the Bode sensitivity integral theorem, reducing low-frequency noise often results in the amplification of noise at other frequencies, a phenomenon commonly known as the "waterbed effect" [4]; secondly, the stability of the system varies across different usage environments, primarily owing to the impact of secondary path uncertainty [5,6] on system stability.Therefore, the basic principle of designing a feedback controller for ANC headphones is to maximize the low-frequency noise reduction performance while robustly maintaining the stability of the system and limit the noise enhancement at high frequencies.
Rafaely [7] pointed out that owing to the additional delays in a digital system, the bandwidth often becomes limited.The controllers for ANC headphones were initially implemented using analog circuits, with various design methods proposed by researchers such as Bai and Lee [8], Pawełczyk [9], and Hu et al. [10,11].However, analog circuits have drawbacks such as significant discrete component errors and temperature drift, deteriorating the noise reduction performance and even system stability, so precise and low-power controller designs are needed in practical applications.
The uncertainty of the secondary path exceeds the deviation of the controller, often leading to stability issues.The modeling errors of the secondary path may deteriorate the stability bound of the system [12].Eriksson and Allie [13] first attempted to use random noise for online identification of a system's secondary path, while Morgan [14] proposed a prediction method applicable for both online and offline identification.However, the method of online identification with added noise is not suitable for ANC headphones as the identification error signal may contaminate the error signal of the noise reduction system itself.Given the close proximity of the headphone system to the human ear, even minor noise can be perceived as discomforting, and the presence of identification noise is unacceptable.Gan discussed using music or speech signals present in the headphone system as identification signals for the secondary path [15].More recently, Yang et al. [16] proposed a method to control the identification error energy, and the system stability was significantly improved, although further research is required for headphone applications.
In practice, the transfer function of a headphone's secondary path is usually obtained through offline identification and is then constrained by disk [17].The uncertainty in the secondary path manifests mainly in the variance between a single headphone's use by a single subject or a dummy head, in terms of wearing and removing the headphones, and the differences in multiple wearings by one or multiple subjects.Zhao et al. [18] noted that putting on and taking off headphones introduces large variance in the transfer functions of the secondary path, affecting the convergence of the adaptive system, and this variance can be reduced using the internal model control (IMC) feedback control implemented with analog circuits.Yu and Hu [10] used the results of four wearings to establish the controller's objective function and constraints and then created an analog circuit controller.Guldenschuh and Callafon [19] found that the low-frequency drastic changes in the secondary path when putting on and removing headphones could be identified using an adaptive filter, based on which a feasible weighting coefficient update rule was established through experiments.Leading companies like BOSE, SONY and Apple have adopted proximity sensors, offering quick and accurate responses in their ANC headphones to monitor the states of wearing and removal [20][21][22], and this approach can replace algorithmic detection from an engineering application perspective, reducing the algorithm's consumption of system power.For the latter case, Ti's engineer Krishnamurthy et al. [23] used the average result of 10 real subjects to obtain the secondary path transfer function and employed a smaller update coefficient to handle the greater variance.Benois and Zolzer [24] constructed different weighting coefficient constraints using real secondary path data from two dummy heads and six subjects wearing the headphones at three different angles and used the sequential quadratic programming (SQP) algorithm [25] to derive the FIR filter and compared the results.Hilgemann and Jay [26] built an uncertainty model using the real part of 78 secondary path transfer functions tested on a single subject and also used the SQP algorithm to calculate the FIR filter.
However, the implementation methods mentioned above use either analog circuits or FIR filters.Analog circuits have significant errors, and the FIR controllers will introduce additional delay and cannot meet the requirements of current ANC headphones.Wang et al. [27] and An et al. [28] implemented controller designs using the IIR filters and optimized them with the genetic optimization algorithm (GA) combined with the Nelder-Mead (NM) algorithm and differential evolution (DE) algorithm, respectively.In the former controller, the application of dual optimization increased the system complexity, while the latter controller based on the DE algorithm operates directly in continuous space, facilitating the global minimum search but bringing high computational demand [29].The GWO [30] algorithm and its improved version, IGWO [31], as typical representatives of swarm intelligence algorithms, are characterized by their simple structure, minor parameter adjustment, and ease of implementation.The IGWO algorithm contains an information feedback mechanism, and a balance is achieved between local optimization and global search, thus offering good performance in both solution accuracy and convergence speed.
In this paper, a method is proposed to design a feedback ANC headphone controller.A new secondary path uncertainty constraint is proposed based on the results of multiple wearing tests on real subjects and a dummy head.This constraint is applied in the process of optimizing the controller, so that new subjects can obtain stable noise reduction and moderate waterbed lift when using the headphone.In Section 2, a safety factor and a waterbed suppression factor are introduced into a proposed uncertainty constraint, which is integrated into the H2/H∞ framework to obtain a feedback control strategic function.In Section 3, the control strategic function is optimized using the IGWO algorithm to obtain the optimal controller composed of cascaded biquad IIR filters, and the noise reduction performance is simulated in two cases with or without the proposed uncertainty constraint.The reliability of the proposed method is validated based on the measured results of the experimental test.Finally, conclusions are drawn in Section 4.

Parametric Filter
A biquad IIR filter has the basic form as follows: The form of each filter coefficient uses the peaking or high-shelf prototypes, whose coefficients [32] are shown in Table 1 below.
In the table, 20) . ( where Fs is the predefined sampling frequency.Among the other variables, Fc denotes the center frequency in Hz; Q represents the steepness of the central frequency and is a dimensionless number; and g denotes the gain of each filter, in dB.These three parameters can uniquely characterize a biquad IIR filter.A positive g means the filter shape is peaking or a high-shelf shape, and a negative g means notching or a low-shelf shape.The feedback controller is composed of cascaded biquad IIR filters: where H n (z) is the nth biquad filter formed by parameter Fc n Q n and g n , and n = 1, 2, . .., N; g total is a real number representing the total gain of the cascade controller.This feedback controller is determined by using the following parameter vector:

Feedback Controller
In Figure 1a, the secondary path G defined in the acoustic domain goes from the headphone speaker (SPK) to the feedback microphone (FB MIC) in the front cavity of the headphone, corresponding to the controller circuit H consisting of a microphone preamplifier, feedback filter, and power amplifier in the electronic domain.The external noise passing through the headphone shell into the chamber is the primary noise d.After the controller is turned on, the remaining error noise is e.In Figure 1b, where d(t) and e(t) are the primary noise signal and error noise signal, respectively, and G(jω) and H(jω) are the transfer functions of the secondary path and controller circuit in the complex frequency domain, respectively.is a real number representing the total gain of the cascade c feedback controller is determined by using the following parameter vector:

Feedback Controller
In Figure 1a, the secondary path G defined in the acoustic domain headphone speaker (SPK) to the feedback microphone (FB MIC) in the fron headphone, corresponding to the controller circuit H consisting of a mi amplifier, feedback filter, and power amplifier in the electronic domain.The passing through the headphone shell into the chamber is the primary noise controller is turned on, the remaining error noise is e .In Figure 1b, where d are the primary noise signal and error noise signal, respectively, and (j G  are the transfer functions of the secondary path and controller circuit in th quency domain, respectively.The system sensitivity function is defined as the transfer function bet mary noise signal and the error noise signal, as follows: is the open-loop transfer function of the system the complementary sensitivity function of the system is defined as follows: The system sensitivity function is defined as the transfer function between the primary noise signal and the error noise signal, as follows: where is the open-loop transfer function of the system.In addition, the complementary sensitivity function of the system is defined as follows: From the point of view of noise suppression, when the primary noise exists, the error signal picked up by the error microphone should be the smallest in terms of mean square. where E[ ] represents the mathematical expectation, and Φ e (ω) and Φ d (ω) are the power spectral densities of e(t) and d(t), respectively.
Appl.Sci.2024, 14, 2236 5 of 13 The objective function and constraint of the H2/H∞ framework is set as follows: where w 1 (ω) is the root mean square value of the power spectral density of primary noise, and w 2 (ω) is the limit of noise lift at each frequency.
From the perspective of system stability, an open-loop transfer function should have a certain amplitude margin, defined as L AM , and phase margin, defined as L PM .The stability constraint of the feedback system is set as follows: For an invariable specified secondary path, the above objective function and constraints, defined as method 1, can ensure that the system can achieve a good balance between noise reduction and waterbed lift, but the secondary path in a real system has uncertainties, so it is necessary to add an uncertainty constraint.

Secondary Path Uncertainty Constraint
The uncertainty discussed in this paper is described by multiplicative uncertainty.The multiplicative form of the secondary path can be expressed as where G 0 (jω) is the initial secondary path transfer function, and the expectation of ∆(jω) represents the degree of multiplicative uncertainty: From a statistical point of view, the degree of uncertainty is minimal when the initial secondary path G 0 (jω) is exactly the expected value of the secondary path E[G(jω)].In practical applications, we should test as many secondary path transfer functions as possible to improve the scope of system stability.In the case of limited numbers of measurement, using the mean value of the secondary path transfer function obtained from all tests helps to reduce the degree of system uncertainty.
In order to ensure the stability of the system in each actual scenario, set w 3 (ω) as the greatest degree of uncertainty ∥∆(jω)∥ ∞ , and because the measured data are always limited, a safety factor δ (where 0 < δ < 1) is set to improve the robustness of the system stability to uncertainty.The larger the number of test samples, the closer δ is to 1.The weighting coefficient of uncertainty constraint is denoted as w 3 (ω)/δ.
As shown in Figure 2, the open-loop transfer function is restricted to the disk whose center is the open-loop transfer function L 0 (jω) and whose radius is |L 0 (jω)w 3 (ω)/δ|, resulting from the fluctuation of the secondary path transfer function.The Nyquist stability condition requires that all points in the disk do not pass through the point (1, j0), and the radius |L 0 (jω)w 3 (ω)/δ| is less than the distance |L 0 (jω) − 1| between the center L 0 (jω) and (1, j0).

Control Strategic Function
Use a specified measured secondary path transfer function and a feedback co with parameter vector 1 X to build a feedback control strategic function.This proc denoted as method 1: Considering the uncertainty of the secondary path, use the mean of multipl ured secondary path transfer functions and a feedback controller with paramete 2 X to rebuild the feedback control strategic function with the uncertainty constrai procedure is denoted as method 2: where () Ux is the unit step function In the Nyquist plot, noise rising is more than 1/χ in the red disk with center (1, j0) and radius χ, which is defined as the waterbed suppressing factor.The open-loop transfer function is restricted in the black disk.When these two disks are non-intersected, the noise rising caused by uncertainty is less than 1/χ; hence, an uncertainty constraint is proposed as below: |L 0 (jω) Equation ( 14) considers the constraint conditions of secondary path uncertainty.In order to prevent the bandwidth of noise reduction being too narrow, χ should be less than w 2 (ω).When χ = 0, the constraint has only met the Nyquist stability but not limited the waterbed lift for the secondary path uncertainty.

Control Strategic Function
Use a specified measured secondary path transfer function and a feedback controller with parameter vector X 1 to build a feedback control strategic function.This procedure is denoted as method 1: Considering the uncertainty of the secondary path, use the mean of multiple measured secondary path transfer functions and a feedback controller with parameter vector X 2 to rebuild the feedback control strategic function with the uncertainty constraint.This procedure is denoted as method 2: where U(x) is the unit step function In Equation ( 16), the first term on the right side represents the amplitude of the sensitivity function, indicating the noise reduction performance.It reaches the minimum when the noise reduction is equal to the frequency characteristic of primary noise w 1 (ω), and the best feedback reduction is achieved.
The second term on the right side of Equation ( 16) represents the constraint on the Nyquist stability.If the open-loop function does not meet the conditions regarding amplitude margin L AM or phase margin L PM , this term is multiplied by a penalty factor γ 1 , which is a positive real number.
The third term pertains to the constraint on waterbed lift of the mean of secondary paths.If the noise enhancement exceeds a preset upper limit 1/w 2 (ω), this term is also multiplied by another positive real number penalty factor γ 2 .
The fourth term relates to the uncertainty constraint of the secondary path.If the constraint is broken through, this term is multiplied by a positive real number penalty factor γ 3 .If a single specified secondary path design is used, this term of the strategic function is omitted.
In this study, the IGWO algorithm is used to solve the control strategic function.The initial parameter range can be customized based on experience, thus offering good performance in terms of solution accuracy and convergence speed.It is suggested that the optimization process should be repeated multiple times to find the global optimized filter coefficients.

Test Setting
In this section, experiments are conducted using ANC headphones designed by Goertek (Qingdao, China) to validate the proposed method.The secondary path transfer function test is shown in Figure 3a, using Audio Precision 525 (AP525) (manufactured by Audio Precision Inc., Beaverton, OR, USA) as the measuring device.The ANC test experimental environment shown in Figure 3b.The experiment is conducted in an anechoic room with a low-frequency cut-off frequency of 80 Hz and free field radius of 1 m.The noise source is generated by a Tannoy's Precision 6 speaker (Tannoy Brand is owned by Music Tribe.Manufactured in Coatbridge, UK), which has an effective frequency response range of 60 Hz to 50 kHz.A B&K's Type 4128C dummy Head and Torso Simulator (HATS) (manufactured in Teknikerbyen 28, Virum, Denmark) is used as the headphone wearing and sound recording device.The headphones are worn on the HATS facing the test speaker, and the distance between the test speaker and the HATS is 0.5 m.The real person ANC test is also conducted, and a microphone is placed at the ear canal to collect the residual noise, as shown in Figure 3c.A controller with the ADAU1777 (manufactured by Analog Devices, Inc., Norwood, MA, USA) as its core component is adopted to implement the feedback controller.The sampling frequency of ADAU1777 is configured to be 768 kHz, and the system latency from ADC to DSP and then to DAC is 4.7 µs, which can be ignored.
Appl.Sci.2024, 14, x FOR PEER REVIEW 8 of 14 wearing states, so it is necessary to conduct the headphone squeeze scenario test.During the test, the headphones are worn in normal and simulated squeeze state alternatively.In the latter state, the headphones are tightened with a rubber band and the strength of the rubber band is adjusted until the headphone cushion has been squished.In order to increase the samples of secondary path, four real subjects and HATS substitute test are conducted under the same experimental conditions.The four subjects each have a head circumference of 52.8 to 59.1 cm.The subjects are two men and two women: one man and woman wear glasses, and the subjects with long hair have their hair pressed when wearing headphones.The parameters of the subjects are shown in Table 2.During the test, the headphones are adjusted to fit the wearing state of each subject's (both HATS and human subjects) headphones in the anechoic room for the first round of test, and then the strength of the rubber band is adjusted to make the headphones squished for the second round of test.In each round of test, the headphones are kept still for three times on each subject and the average value is recorded.The test results are shown in Figure 4a,b, displaying the amplitude and phase response of the secondary path transfer function, respectively, under the normal wearing state of the first round.In gen- During the secondary path transfer function test, the AP525's left output is connected to the speaker of the headphone, and the right input is connected to the error microphone of the headphone.The left input is directly connected to the right output to offset the delay caused by the device's own ADC and DAC.The secondary path transfer function is tested using the steady-state noise signals.
As current commercial ANC headphones use a proximity sensor to detect if the headphones are being worn, this experiment only focuses on the wearing state of headphones.However, the proximity sensor cannot distinguish between squeezed wearing and normal wearing states, so it is necessary to conduct the headphone squeeze scenario test.During the test, the headphones are worn in normal and simulated squeeze state alternatively.In the latter state, the headphones are tightened with a rubber band and the strength of the rubber band is adjusted until the headphone cushion has been squished.In order to increase the samples of secondary path, four real subjects and HATS substitute test are conducted under the same experimental conditions.The four subjects each have a head circumference of 52.8 to 59.1 cm.The subjects are two men and two women: one man and woman wear glasses, and the subjects with long hair have their hair pressed when wearing headphones.The parameters of the subjects are shown in Table 2.During the test, the headphones are adjusted to fit the wearing state of each subject's (both HATS and human subjects) headphones in the anechoic room for the first round of test, and then the strength of the rubber band is adjusted to make the headphones squished for the round of test.In each round of test, the headphones are kept still for three times on each subject and the average value is recorded.The test results are shown in Figure 4a,b, displaying the amplitude and phase response of the secondary path transfer function, respectively, under the normal wearing state of the first round.In general, large differences are present in the low frequency range, and the results of male subjects have larger low-frequency amplitude than those of the females, since the male head circumference is larger and the headphones are more easily squeezed and sealed against the head.Air leakage may be created when the headphones press glasses legs or long hair, resulting in low-frequency attenuation.In Figure 4c,d, results of the second round with squished wearing of headphones are shown.In general, the headphones are almost sealed, so the amplitude difference is smaller in the low frequency range.The test results of the HATS in both rounds (green line) demonstrate a more pronounced rise and dip caused by high-frequency reflections because the hardness of HATS is different from that of human subjects.The measured secondary path results show fluctuation across different su likely causing a system stability risk if designed with a specified secondary path.To this, two optimization methods are proposed: the first method uses a normal wearin The measured secondary path results show fluctuation across different subjects, likely causing a system stability risk if designed with a specified secondary path.To verify this, two optimization methods are proposed: the first method uses a normal wearing test result on the HATS as the specified secondary path transfer function; in the second method, the mean of all test results is used as the secondary path transfer function.The secondary path transfer functions used in both methods are shown below in Figure 5.The measured secondary path results show fluctuation across different su likely causing a system stability risk if designed with a specified secondary path.To this, two optimization methods are proposed: the first method uses a normal wearin result on the HATS as the specified secondary path transfer function; in the s method, the mean of all test results is used as the secondary path transfer function secondary path transfer functions used in both methods are shown below in Figure

Controller Optimization and Validation
Firstly, for frequency discretization, the sampling rate is set to be 40 kHz, corres ing to 512 frequency points in the frequency domain.The weighting coefficient 1 ( w set to 15 from 117.2 Hz to 390.6 Hz and 1 for 39.1 Hz and 1.21 kHz, and the values remaining frequency points are obtained using discretized frequency interpolation is shown in dB by the blue curve in Figure 6a.The weighting coefficient 2 () w  is a fixed value 10 − from 1.21 kHz to 20 kHz.For method 2, the weighting coef 3 () w  is calculated by all the measured secondary path transfer functions and is s

Controller Optimization and Validation
Firstly, for frequency discretization, the sampling rate is set to be 40 kHz, corresponding to 512 frequency points in the frequency domain.The weighting coefficient w 1 (ω) is set to 15 from 117.2 Hz to 390.6 Hz and 1 for 39.1 Hz and 1.21 kHz, and the values at the remaining frequency points are obtained using discretized frequency interpolation.This is shown in dB by the blue curve in Figure 6a.The weighting coefficient w 2 (ω) is given a fixed value 10 (−3/20) from 1.21 kHz to 20 kHz.For method 2, the weighting coefficient w 3 (ω) is calculated by all the measured secondary path transfer functions and is shown in Figure 6b.The safety factor δ is set to be 0.9 and the waterbed suppressing factor χ is set to be 10 (−5 /20) .The penalty factors are set to be the same: γ 1 = γ 2 = γ 3 = 10,000, for simplifying the optimization.
Appl.Sci.2024, 14, x FOR PEER REVIEW 10 of in Figure 6b.The safety factor  is set to be 0.9 and the waterbed suppressing factor is set to be  For a cascaded IIR controller, the function of H are used to process the waterbed effect in the mi dle frequency band; and 5 H ensures the stability of the ultra-high frequency band abo 10 kHz.The upper and lower limits of the filter parameters ub and lb are configured Table 3 according to the author's experience.In the IGWO algorithm, the population si is set to be 50, the dimension is 16, there are 10,000 iteration times, and the optimizati process is repeated 100 times to obtain the global optimal solution.The filter paramete For a cascaded IIR controller, the function of H 1 is to achieve more noise reduction at low frequencies; H 2 , H 3 and H 4 are used to process the waterbed effect in the middle frequency band; and H 5 ensures the stability of the ultra-high frequency band above 10 kHz.The upper and lower limits of the filter parameters ub and lb are configured in Table 3 according to the author's experience.In the IGWO algorithm, the population size is set to be 50, the dimension is 16, there are 10,000 iteration times, and the optimization process is repeated 100 times to obtain the global optimal solution.The filter parameters calculated using the two methods are shown in Table 3.The controller responses of the two methods are shown in Figure 7.  Method 1′s secondary path transfer function's amplitude being less pronoun the frequency range of several hundred to three thousand Hertz compared to met necessitates a higher amplitude for the controller to compensate.Method 2 restric amplitude beyond 10 kHz to guarantee system stability across various user scenario preserves a larger phase margin in the 2 to 5 kHz range to effectively manage the wat effect.Next, in addition to the HATS and four subjects on whom the secondary path fer functions have been tested, four additional subjects are added for the secondary and noise reduction test, and their parameters are shown in Table 4.  Method 1 ′ s secondary path transfer function's amplitude being less pronounced in the frequency range of several hundred to three thousand Hertz compared to method 2 necessitates a higher amplitude for the controller to compensate.Method 2 restricts the amplitude beyond 10 kHz to guarantee system stability across various user scenarios and preserves a larger phase margin in the 2 to 5 kHz range to effectively manage the waterbed effect.Next, in addition to the HATS and four subjects on whom the secondary path transfer functions have been tested, four additional subjects are added for the secondary path and noise reduction test, and their parameters are shown in Table 4. Noise reduction is simulated for all the measured secondary path transfer functions, with the controllers optimized using method 1 and method 2, and the results are shown in Figures 8a and 8c, respectively.In Figure 8a, the green line denotes the noise reduction simulation result for the secondary path transfer function of the headphones worn normally on the HATS.The gray and blue curves depict the simulated results for other secondary paths under normal and squeezed wearing states, respectively.Figure 8c illustrates the noise reduction simulation results for the secondary path under normal wearing state with the gray and black curves representing these outcomes.In contrast, the blue and cyan curves correspond to the simulated results under a squeezed wearing state, where the black and cyan curves specifically indicate the results for the four newly included subjects.The simulated results shown in Figure 8a suggest that method 1 will have a relatively serious waterbed effect and stability risk at high frequencies (as indicated by the arrows), while the results in Figure 8c suggest that method 2 will have a relatively mild waterbed lift (controlled within 5 dB) under all measured secondary paths.This comparison indeed underscores the advantages of method 2 in noise reduction and system stability.It also showcases method 2's robustness and adaptability, as it maintains consistent effectiveness even with the addition of new subjects who were not included in the initial controller design.Subsequently, the noise reduction is also measured with the controllers optimized using method 1 and method 2, and the results are shown in Figure 8b,d, respectively.The overall measured results are generally in agreement with the simulation predictions, indicating a strong correlation between theoretical predictions and practical outcomes.Across all subjects, method 2 has demonstrated stable performance with a moderate waterbed lift; this suggests that method 2 offers a reliable approach to commercial ANC headphones.
Appl.Sci.2024, 14, x FOR PEER REVIEW 12 consistent effectiveness even with the addition of new subjects who were not include the initial controller design.Subsequently, the noise reduction is also measured with controllers optimized using method 1 and method 2, and the results are shown in Fig 8b,d, respectively.The overall measured results are generally in agreement with the s ulation predictions, indicating a strong correlation between theoretical predictions practical outcomes.Across all subjects, method 2 has demonstrated stable performa with a moderate waterbed lift; this suggests that method 2 offers a reliable approac commercial ANC headphones.

Conclusions
In this paper, an uncertainty constraint is proposed based on multiple measured ondary path transfer functions obtained from different subjects and a dummy hea different wearing states.It includes a safety factor, which is used to make the second path transfer function obtained from a limited number of subjects compatible with m subjects, and a waterbed suppressing factor, which is introduced to limit the waterbed With the aid of this constraint, the optimal parameters of feedback controller compo of cascade biquad IIR filters are found with the IGWO algorithm, and experiments

Conclusions
In this paper, an uncertainty constraint is proposed based on multiple measured secondary path transfer functions obtained from different subjects and a dummy head in different wearing states.It includes a safety factor, which is used to make the secondary path transfer function obtained from a limited number of subjects compatible with more subjects, and a waterbed suppressing factor,

Figure 1 .
Figure 1.(a) Schematic drawing of feedback control headphone; (b) block diagram o tem.

Figure 1 .
Figure 1.(a) Schematic drawing of feedback control headphone; (b) block diagram of feedback system.

13 0
Appl.Sci.2024, 14, 2236 6 of Equation (14) considers the constraint conditions of secondary path uncerta order to prevent the bandwidth of noise reduction being too narrow,  should than 2 () w  .When 0  = , the constraint has only met the Nyquist stability but n ited the waterbed lift for the secondary path uncertainty.

Figure 3 .
Figure 3. (a) Test setting for secondary path; (b) ANC test environment; (c) microphone in the ear of subject for ANC test.

Figure 3 .
Figure 3. (a) Test setting for secondary path; (b) ANC test environment; (c) microphone in the ear of subject for ANC test.

Figure 4 .
Figure 4. Amplitude and phase response of tested secondary path transfer function (a,b) und mal wearing state and (c,d) under squeezed wearing state.

Figure 4 .
Figure 4. Amplitude and phase response of tested secondary path transfer function (a,b) under normal wearing state and (c,d) under squeezed wearing state.

Figure 4 .
Figure 4. Amplitude and phase response of tested secondary path transfer function (a,b) und mal wearing state and (c,d) under squeezed wearing state.

Figure 5 .
Figure 5. Secondary path transfer functions used in the two methods, (a) amplitude response phase response.

Figure 5 .
Figure 5. Secondary path transfer functions used in the two methods, (a) amplitude response, (b) phase response.

( 5/ 20 ) 10 −
. The penalty factors are set to be the same:

Figure 7 .
Figure 7. (a) Amplitude response of controller; (b) phase response of controller.

Table 2 .
Parameter of subjects.

Table 2 .
Parameter of subjects.

Table 3 .
Parameters of the controller.

Table 4 .
Parameter of new subjects.

Table 4 .
Parameter of new subjects.