Feedback active noise control system and strategy with online secondary-path modeling

Abstract

The present disclosure presents a feedback active noise control system and strategy with online secondary-path modeling, and belongs to the technical field of active noise control. The linear prediction subsystem takes the residual noise as its input and separates the remaining sinusoidal noise from the broadband noise. The remaining sinusoidal noise is used effectively not only to update the controller but also to scale the auxiliary noise, while the broadband noise serves as a desired input of online secondary-path modeling subsystem. In this way, the coupling between the controller and the online secondary-path modeling subsystem is significantly mitigated, leading to both faster convergence and improved noise reduction performance. A practical scheme for refreshing the entire system is also developed to enhance its robustness against even abrupt changes with the secondary path or the primary noise. The present disclosure enhances the applicability of feedback active noise control in practical applications.

Claims

1. A feedback active noise control system with online secondary-path modeling, the feedback active noise control system comprising: a reference signal synthesis subsystem (1), a secondary sound source synthesis subsystem (2), a linear prediction subsystem (3) and an online secondary-path modeling subsystem (4); wherein the reference signal synthesis subsystem (1) is separately connected to the secondary sound source synthesis subsystem (2) and the linear prediction subsystem (3); the secondary sound source synthesis subsystem (2) is separately connected to the reference signal synthesis subsystem (1) and the online secondary-path modeling subsystem (4); the linear prediction subsystem (3) is separately connected to the reference signal synthesis subsystem (1), the secondary sound source synthesis subsystem (2) and the online secondary-path modeling subsystem (4); and the online secondary-path modeling subsystem (4) is separately connected to the secondary sound source synthesis subsystem (2) and the linear prediction subsystem (3); the reference signal synthesis subsystem (1) is used for synthesizing a reference signal; the secondary sound source synthesis subsystem (2) is used for synthesizing a secondary sound source; the linear prediction subsystem (3) is used for separating a narrowband component and a broadband component from residual noise; the online secondary-path modeling subsystem (4) is used for estimating a time-varying secondary path estimation model on line in real time; the narrowband component separated from the residual noise by the linear prediction subsystem (3) is used for adjusting an amplitude of auxiliary white Gaussian noise, reducing the contribution of injected auxiliary noise to the residual noise, and improving the noise suppression performance of the system; the broadband component and the narrowband component separated from the residual noise by the linear prediction subsystem (3) are respectively used as a desired input of the online secondary-path modeling subsystem (4) and an error signal for the secondary sound source synthesis subsystem (2) so as to improve the independence between a controller and an online secondary-path modeling module, improve the accuracy and speed of online secondary-path modeling, and improve a dynamic performance of the system at the same time; the feedback active noise control system monitoring a possible sudden change of a secondary path or target noise by calculating in real time an energy change of the residual noise after smoothing filtering, and re-initializing a coefficient of the linear prediction subsystem (3), a coefficient of the secondary path estimation model, a coefficient of the secondary sound source synthesis subsystem (2) and a scaling factor of the online secondary-path modeling subsystem (4) so as to improve the capability of the system to deal with a large sudden change of the secondary path or target noise, and improve the robust performance of the feedback active noise control system; the energy of the residual noise after smoothing filtering being:
P.sub.e(n)=λ.sub.mP.sub.e(n−1)+(1−λ.sub.m)e.sup.2(n) wherein n is time instant, n≥0, and λ.sub.m ∈(0,1) is a forgetting factor of smoothing filtering; at time instant n′T.sub.p, the following being obtained by successively performing time averaging and smoothing filtering on the energy P.sub.e(n) of the residual noise after smoothing filtering: $P_{e, T_{p}} (n^{'}) = λ_{m} P_{e, T_{p}} (n^{'} - 1) + (1 - λ_{m}) {.Math.}_{k = 0}^{T_{p} - 1} P_{e} (n - k)$ wherein n′ is a positive integer greater than 1 when n is evenly divided by T.sub.p, and T.sub.p is the length of a time average window; when P.sub.e,T.sub.p(n′)≥αP.sub.e,T.sub.p(n′−1) is satisfied at time instant n, the system performing re-initialization at the moment n+1, wherein α ∈ (1,2) is a threshold parameter.

2. The system according to claim 1, wherein the linear prediction subsystem (3) comprises a D-order delay operator (31) and a linear prediction filter (32), the D-order delay operator (31) and the linear prediction filter (32) are connected in series, a coefficient and a length of the linear prediction filter (32) are respectively {h.sub.j(n)}.sub.j=0.sup.L−1 and L, the coefficient is updated using a least mean square algorithm, and the update formula is:
h.sub.j(n+1)=h.sub.j(n)+μ.sub.he.sub.LP(n)e(n−D−j) wherein μ.sub.h is the update step size of the linear prediction filter and takes a positive value; D is a delay order number; e.sub.LP(n) is the broadband component separated by the linear prediction subsystem (3), and e(n) is the residual noise.

3. The system according to claim 2, wherein the broadband component separated from the residual noise is: $e_{L P} (n) = e (n) - y_{L P} (n) y_{L P} (n) = {.Math.}_{j = 0}^{L - 1} h_{j} (n) e (n - D - j)$ wherein y.sub.LP(n) is the narrowband component separated from the residual noise.

4. The system according to claim 3, wherein the online secondary-path modeling subsystem (4) comprises an online secondary-path modeling module (41) and an auxiliary noise adjustment module (42); the online secondary-path modeling module (41) comprises a secondary path estimation model Ŝ.sub.n(z); the online secondary-path modeling module (41) takes the broadband component as a desired input, takes colored noise v(n) generated by white Gaussian noise after passing through the auxiliary noise adjustment module (42) as a reference input, and uses a least mean square algorithm to estimate and update the time-varying secondary path estimation model on line in real time; a coefficient and length of the secondary path estimation model Ŝ.sub.n(z) of the online secondary-path modeling module (41) are respectively {ŝ.sub.m(n)}.sub.m=0.sup.{circumflex over (M)}−1 and {circumflex over (M)}, and a coefficient update formula is:
Ŝ.sub.m(n+1)=Ŝ.sub.m(n)+μ.sub.se.sub.s(n)v(n−m)
e.sub.s(n)=e.sub.LP(n)−y.sub.s(n) wherein μ.sub.s is the update step size of the secondary path estimation model and takes a positive value; and y.sub.s(n) is an output of the secondary path estimation model of the online secondary-path modeling module (41); the colored noise v(n) is:
v(n)=v.sub.0(n)G.sub.s(n)
G.sub.s(n)=λG.sub.s(n−1)+(1−λ)y.sub.LP.sup.2(n−1) wherein G.sub.s(n) is a scaling factor of the auxiliary noise adjustment module (42); λ is a forgetting factor of the auxiliary noise adjustment module, λ∈(0,1); and v.sub.0(n) is additive white Gaussian noise with a mean value of zero and a variance of σ.sub.0.sup.2.

5. The system according to claim 4, wherein the reference signal synthesis subsystem (1) comprises a secondary path estimation model (11) and a first-order delay operator (12), the secondary path estimation model (11) being provided by the online secondary-path modeling module (41); the reference signal is:
x(n)=e(n−1)+ŷ.sub.0(n−1) wherein e(n−1) is an output of the residual noise e(n) via the first-order delay operator (12); ŷ.sub.0(n) is an output of y.sub.0(n) via the secondary path estimation model (11); and ŷ.sub.0(n−1) is an output of ŷ.sub.0(n) via the first-order delay operator (12).

6. The system according to claim 5, wherein the secondary sound source synthesis subsystem (2) comprises a controller (21) and a filtering-X least mean square algorithm module (22); the filtering-X least mean square algorithm module (22) takes the narrowband component y.sub.LP(n) separated from the residual noise as an error output to be used for updating the coefficient of the controller (21).

7. The system according to claim 6, wherein the controller (21) employs a linear filter, a coefficient and a length of the linear filter being respectively {w.sub.i(n)}.sub.i=0.sup.M.sup.w.sup.−1 and M.sub.w; the coefficient update formula of the controller (21) is:
w.sub.i(n+1)=w.sub.i(n)+μ.sub.wy.sub.LP(n){circumflex over (x)}(n−i) wherein μ.sub.w is the update step size of the controller and takes a positive value; y.sub.LP(n) is the narrowband component separated by the linear prediction subsystem (3); and {circumflex over (x)}(n) is an output of a reference signal x(n) via a secondary path estimation model of the filtering-X least mean square algorithm module (22).

8. The system according to claim 7, wherein a secondary sound source is:
y(n)=y.sub.0(n)−v(n) wherein y.sub.0(n) is an output of the controller (21).

9. An active noise control method, the method being implemented on the basis of the feedback active noise control system with online secondary-path modeling according to claim 8, and the method comprises: step 1: setting system parameters: setting the length and the update step size of the controller (21), the length and the update step size of the linear prediction filter (32), and the length and the update step size of the secondary path estimation model Ŝ.sub.n(z); setting the order number D of a delay operator; setting a forgetting factor of the auxiliary noise adjustment module (42); setting a forgetting factor, a threshold parameter and the length of a time average window required for re-initialization of the system; setting initial values of the coefficients of the controller (21) and the secondary path estimation model Ŝ.sub.n(z), the coefficient of the linear prediction filter (32), and the scaling factor of the auxiliary noise adjustment module (42) all to be zero; step 2: synthesizing the reference signal: adding residual noise e(n) obtained by an error microphone with an output ŷ.sub.0(n) of an output y.sub.0(n) of the controller (21) via the secondary path estimation model (11), and obtaining the reference signal x(n) after an obtained signal is subjected to the first-order delay operator (12);
x(n)=e(n−1)+ŷ.sub.0(n−1) that is, summing the residual noise at time instant n−1 and an output signal of the secondary path estimation model (11) to obtain the reference signal at time instant n by synthesis; step 3: at time instant n, firstly, obtaining y.sub.0(n), by the reference signal x(n), via the controller (21); then, obtaining auxiliary noise v(n) by the auxiliary noise adjustment module (42), and then obtaining a secondary sound source y(n) by synthesis; and finally, separating the residual noise e(n) into a narrowband component y.sub.LP(n) and a broadband component e.sub.LP(n) by the linear prediction subsystem (3); step 4: updating the control system: calculating and updating the coefficient of the controller (21) at time instant n+1 according to the reference signal and the narrowband component y.sub.LP(n); calculating and updating the coefficient of the linear prediction filter (32) at time instant n+1 according to the residual noise e(n) and the narrowband component y.sub.LP(n); calculating and updating the coefficient of the secondary path estimation model Ŝ.sub.n(z) at time instant n+1 according to the auxiliary noise v(n) and the broadband component e.sub.LP(n); and updating a scaling factor of the auxiliary noise adjustment module (42) at time instant n+1 according to the narrowband component y.sub.LP(n); step 5: calculating the energy change of the residual noise after smoothing filtering in real time, i.e., if P.sub.e,T.sub.p(n′)≥αP.sub.e,T.sub.p(n′−1) is satisfied, re-initializing the coefficient of the linear prediction filter (32), the coefficient of the secondary path estimation model Ŝ.sub.n(z), the scaling factor of the auxiliary noise adjustment module (42) and the coefficient of the controller (21) at time instant n+1, and then proceeding to step 6; if P.sub.e,T.sub.p(n′)≥α.sub.e,T.sub.p(n′−1) is not satisfied, then directly proceeding to step 6; and step 6: returning to step 2, and repeating steps 2-5 described above until the system converges and reaches a steady state.

Description

BRIEF DESCRIPTION OF FIGURES

(1) To describe the technical solutions in the examples of the present disclosure more clearly, the following briefly introduces the accompanying drawings required for describing the examples. Apparently, the accompanying drawings in the following description show merely some examples of the present disclosure, and those of ordinary skill in the art may still derive other accompanying drawings from these accompanying drawings without creative efforts.

(2) FIG. 1 is a schematic diagram of a feedback active noise control system with online secondary-path modeling of Example 1;

(3) FIG. 2A is a graph showing the change of residual noise powers of Example 3;

(4) FIG. 2B is a graph showing the change of mean square error (MSE) of secondary path modeling (SPM) of Example 3;

(5) FIG. 2C is a graph showing the change of an auxiliary noise scaling factor of Example 3;

(6) FIG. 3A is a graph showing the change of target noise and residual noise of Example 4; and

(7) FIG. 3B is a graph showing the change of an auxiliary noise scaling factor of Example 4.

DETAILED DESCRIPTION

(8) To make the objectives, technical solutions and advantages of the present disclosure more comprehensible, detailed description is made to specific embodiments of the present disclosure below with reference to the accompanying drawings.

EXAMPLE 1

(9) This example provides a feedback active noise control system with online secondary-path modeling. With reference to FIG. 1, the active noise control system includes a reference signal synthesis subsystem (1), a secondary sound source synthesis subsystem (2), a linear prediction subsystem (3) and an online secondary-path modeling subsystem (4).

(10) The reference signal synthesis subsystem (1) synthesizes a reference signal by superposition of residual noise at the previous moment and a secondary sound source. The secondary sound source synthesis subsystem (2) uses a linear filter as a controller, and the output of the controller is added to the output of an auxiliary noise adjustment module (42) to synthesize the secondary sound source. The linear prediction subsystem (3) is composed of a D-order delay operator (31) and a linear prediction filter (32) in a serial connection manner, so as to separate a narrowband component and a broadband component from the residual noise. With the operation of the feedback active noise control system, the online secondary-path modeling subsystem (4) is used for estimating the time-varying secondary path model on line in real time to improve the stability of the system.

(11) The target noise is:

(12) $\begin{matrix} p (n) = p_{0} (n) + v_{p} (n) = {.Math.}_{i = 1}^{q} [A_{i} \cos (ω_{p, i} n + θ_{i})] + v_{p} (n) \end{matrix}$

(13) where p.sub.0(n) is the narrowband noise component in the target noise; q is the number of narrowband components in the target noise, and {A.sub.i}.sub.i=1.sup.q is the amplitude of the narrowband component; ω.sub.p,i is the frequency of an i.sup.th narrowband component in the target noise; θ.sub.i is an initial phase of the i.sup.th narrowband component; v.sub.p(n) is additive white Gaussian noise with a mean value of zero and a variance of σ.sub.p.sup.2; and n is time instant time, n≥0.

(14) An actual secondary path S(z) represents an acoustic spatial model from a secondary speaker to an error microphone, which may be represented by using a finite impulse response filter or an infinite impulse response filter.

(15) The difference between the target noise p(n) and a signal y.sub.p(n) of a secondary sound source y(n) after passing through the actual secondary path S(z) is the residual noise, i.e., e(n)=p(n)−y.sub.p(n).

(16) The reference signal synthesis subsystem (1) includes a secondary path estimation model and a first-order delay operator. The residual noise e(n) collected by an error microphone and the output ŷ.sub.p(n) of the output y.sub.0(n) of the controller via the secondary path estimation model (11) are added, and an obtained signal can be synthesized into a reference signal after passing through the first-order delay operator (12), i.e.,
x(n)=e(n−1)+{circumflex over (y)}.sub.0(n−1)

(17) where the secondary path estimation model (11) is provided by an online secondary-path modeling module (41).

(18) The secondary sound source synthesis subsystem (2) includes a controller (21), a filtering-X least mean square algorithm module (22) and a linear prediction compensation model (23) with a D-order delay. The controller (21) employs a linear filter, and its coefficient and length are respectively {w.sub.i(n)}.sub.i=0.sup.M.sup.w.sup.−1 and M.sub.w. The filtering-X least mean square algorithm module (22) is used for updating the coefficient of the controller (21), i.e.,
w.sub.i(n+1)=w.sub.i(n)+μ.sub.wy.sub.LP(n){circumflex over (x)}(n−i)

(19) where μ.sub.w is the update step size of the controller and takes a positive value; and y.sub.LP(n) is the narrowband component separated by the linear prediction subsystem (3). A signal {circumflex over (x)}(n) is obtained from the reference signal x(n) via a secondary path estimation model Ŝ.sub.n(z). The output of the controller (21) and the output of the auxiliary noise adjustment module (42) in the online secondary-path modeling subsystem (4) are added, and synthesized to obtain a secondary sound source, i.e., y(n)=y.sub.0(n)—v(n).

(20) The linear prediction subsystem (3) is composed of aD-order delay operator (31) and a linear prediction filter (32) in a serial connection manner. The linear prediction filter (32) is represented by H(z), its coefficient and the length are respectively {h.sub.j(n)}.sub.j=0.sup.L−1 and L, the coefficient is updated with a least mean square algorithm, i.e.,
h.sub.j(n+1)=h.sub.j(n)+μ.sub.he.sub.LP(n)e(n−D−j)

(21) in the formula, μ.sub.h is the update step size of the linear prediction filter and takes a positive value; and e.sub.LP(n) is the broadband component separated by the linear prediction subsystem (3), i.e., the difference between the residual noise and the output of the linear prediction filter (32): e.sub.LP(n)=e(n)−y.sub.LP(n), where

(22) $y_{L P} (n) = {.Math.}_{j = 0}^{L - 1} h_{j} (n) e (n - D - j) .$
The linear prediction subsystem (3) is used for separating the narrowband component y.sub.LP(n) and the broadband component e.sub.LP(n) from the residual noise.

(23) The online secondary-path modeling subsystem (4) includes an online secondary-path modeling module (41) and an auxiliary noise adjustment module (42). The online secondary-path modeling module (41) takes the broadband component separated by the linear prediction subsystem (3) as a desired input e.sub.LP(n), takes colored noise v(n) generated by the auxiliary white Gaussian noise v.sub.0(n) after passing through the auxiliary noise adjustment module (42) as a reference input, and uses a least mean square algorithm to estimate a time-varying secondary path model on line in real time. The coefficient and the length of a corresponding secondary path estimation model are respectively {ŝ.sub.m(n)}.sub.m=0.sup.{circumflex over (M)}−1 and {circumflex over (M)}, and the coefficient update formula is:
Ŝ.sub.m(n+1)=Ŝ.sub.m(n)+μ.sub.se.sub.s(nv(n=m)
e.sub.s(n)=e.sub.LP(n−y.sub.s(n)

(24) in the formula, μ.sub.s is the update step size of the secondary path estimation model and takes a positive value. The stability of the system is improved. The auxiliary noise adjustment module (42) takes the narrowband component y.sub.LP(n) separated by the linear prediction subsystem (3) as an input, and the scaling factor is expressed as:
G.sub.s(n)=λG.sub.s(n−1)+(1−λ)y.sub.LP.sup.2(n−1)

(25) in the formula, a forgetting factor λ∈(0,1) of the auxiliary noise adjustment module usually takes a value close to 1. Then, the colored noise generated by the auxiliary white Gaussian noise v.sub.0(n) after passing through the auxiliary noise adjustment module (42) is v(n)=v.sub.0(n)G.sub.s(n), where v.sub.0(n) is additive white Gaussian noise with a mean value of zero and a variance of σ.sub.0.sup.2.

(26) The system monitors a possible sudden change of the secondary path or target noise by calculating the energy change of the residual noise after smoothing filtering in real time, and re-initializes the coefficient of the linear prediction filter (32), the coefficient of the secondary path estimation model Ŝ.sub.n(z), the coefficient of the controller (21) and the scaling factor of the auxiliary noise adjustment module (42).

(27) The energy of the residual noise after smoothing filtering is:
P.sub.e(n)=λ.sub.mP.sub.e(n−1)+(1−λ.sub.m)e.sup.2(n)

(28) where λ.sub.m ∈(0,1) is a forgetting factor of smoothing filtering.

(29) At time instant n′T.sub.p, the following is obtained by successively performing time averaging and smoothing filtering on the energy P.sub.e(n) of the residual noise after smoothing filtering:
P.sub.e,T.sub.p(n′)=λ.sub.mP.sub.e,T.sub.p(n′−1)+(1−λ.sub.m)Σ.sub.k=0.sup.T.sup.p.sup.−1P.sub.e(n−k)

(30) where n′ is a positive integer greater than 1 when n is evenly divided by T.sub.p, and T.sub.p is the length of a time average window.

(31) When P.sub.e,T.sub.p(n′)≥αP.sub.e,T.sub.p(n′−1) is satisfied at time instant n, the system performs re-initialization at time instant n+1, where α∈(1, 2) is a threshold parameter.

EXAMPLE 2

(32) This Example provides a feedback active noise control strategy with online secondary-path modeling, and the strategy is achieved based on the above active noise control with online secondary-path modeling, and includes

(33) Step 1: Setting System Parameters

(34) setting the length and the update step size of the controller (21), the length and the update step size of the linear prediction filter (32), and the length and the update step size of the secondary path estimation model Ŝ.sub.n(z); setting the order number D of a delay operator; setting a forgetting factor of the auxiliary noise adjustment module (42); setting a forgetting factor, a threshold parameter and the length of a time average window required for re-initialization of the system; setting initial values of the coefficients of the controller (21) and the secondary path estimation model Ŝ.sub.n(z), and the coefficient of the linear prediction filter (32) all to be zero;

(35) Step 2: Synthesizing the Reference Signal

(36) adding residual noise e(n) obtained by an error microphone with an output ŷ.sub.0(n) of an output y.sub.0(n) of the controller (21) via the secondary path estimation model (11), and obtaining the reference signal x(n) after an obtained signal is subjected to the first-order delay operator (12), i.e. x (n)=e(n−1)+ŷ.sub.0(n−1), i.e., summing the residual noise at time instant n−1 and an output signal of the secondary path estimation model (11) to obtain the reference signal at time instant n by synthesis;

(37) Step 3: at time instant n, firstly, obtaining y.sub.0(n), by the reference signal x(n), via the controller (21); then, obtaining auxiliary noise v(n) by the auxiliary noise adjustment module (42), and then obtaining a secondary sound source y(n) by synthesis; and finally, separating the residual noise e(n) into a narrowband component y.sub.LP(n) and a broadband component e.sub.LP(n) by the linear prediction subsystem (3);

(38) Step 4: Updating the Control System

(39) calculating and updating the coefficient of the controller (21) at time instant n+1 according to the reference signal and the narrowband component y.sub.LP(n);

(40) calculating and updating the coefficient of the linear prediction filter (32) at time instant n+1 according to the residual noise e(n) and the narrowband component y.sub.LP(n);

(41) calculating and updating the coefficient of the secondary path estimation model Ŝ.sub.n(z) at time instant n+1 according to the auxiliary noise v(n) and the broadband component e.sub.LP(n); and

(42) updating a scaling factor of the auxiliary noise adjustment module (42) at time instant n+1 according to the narrowband component y.sub.LP(n);

(43) step 5: calculating the energy change of the residual noise after smoothing filtering in real time, i.e., if P.sub.e,T.sub.p(n′)≥αP.sub.e,T.sub.p(n′−1) is satisfied, re-initializing the coefficient of the linear prediction filter (32), the coefficient of the secondary path estimation model Ŝ.sub.n(z), the scaling factor of the auxiliary noise adjustment module (42) and the coefficient of the controller (21) at time instant n+1, and then proceeding to step 6; if P.sub.e,T.sub.p(n′)÷αP.sub.e,T.sub.p(n′−1) is not achieved, then directly proceeding to step 6.

(44) step 6: returning to step 2, and repeating the steps 2-5 described above until the system converges and reaches a steady state to achieve active noise control.

EXAMPLE 3

Verification in Case of Simulated Noise and Simulated Secondary Path

(45) The target noise is composed of five frequency components and additive white Gaussian noise. The normalized angular frequencies of the five frequency components are respectively 0.10π, 0.15π, 0.20π, 0.25π and 0.30π. The amplitudes of the corresponding frequency components are respectively 1.41, 1.00, 0.50, 0.25 and 0.10. The additive white Gaussian noise has a mean value of zero and a variance of 0.10.

(46) In order to simulate the large sudden change of a secondary path, a linear FIR model is adopted for an actual secondary path S(z), its cut-off frequency is 0.5 π, and the model lengths of the first half and the second half are 51 and 31, respectively. The length of a secondary path estimation model is 53, and the update step size of the corresponding coefficient is 0.0005. Auxiliary white Gaussian noise v.sub.0(n) has a mean value of zero and a variance of 0.25. The forgetting factor of an auxiliary noise adjustment module (42) is 0.9995. A delay length of the D-order delay operator (31) is 55; the length of the linear prediction filter (32) is 128, and the update step size of the coefficient is 0.001. The controller (21) adopts a linear filter, its length is 128, and the update step size of the coefficient is 0.000075. λ.sub.m, α and T.sub.p are 0.98, 1.1 and 20, respectively. The number of independent operations is 100. The length of simulation data is 60000.

(47) FIG. 2A is a plot showing the change of target noise and residual noise in the case of simulated noise and a simulated secondary path of Example 3. When the system reaches a steady state, the noise reduction of the first half and the second half are 10.84 dB and 10.46 dB, respectively, and the corresponding residual noise energy of the system are about 0.15 and 0.16, respectively, which are close to the variance of the additive white Gaussian noise in the target noise, i.e. approaching to the ambient noise level, with good target noise suppression performance. FIG. 2B is a plot showing the change of mean square error of secondary-path modeling in this case, and FIG. 2C is a plot showing the change of an auxiliary noise scaling factor in this case, which together show that the system of the present disclosure can not only effectively track the large sudden change of the secondary path, but also has high online secondary-path modeling accuracy.

EXAMPLE 4

Verification in Case of Actual Noise and Actual Secondary Path

(48) Actual noise is from a discharge port of a large cutting machine under working conditions. In order to simulate the larger sudden changes of the target noise, the target noise is divided into the first half and the second half, the first half corresponding to a rotation speed at 1,400 rpm and the second half corresponding to a rotation speed at 1,600 rpm. The actual secondary path is an IIR model that is widely adopted by peers (S. M. Kuo and D. R. Morgan, Active Noise Control Systems-Algorithms and DSP Implementation, New York: Wiley, 1996.). The length of a secondary path estimation model Ŝ(z) is 32 and a corresponding update step size of the coefficient is 0.4. Auxiliary white Gaussian noise v.sub.0(n) has a mean value of zero and a variance of 1.0. The forgetting factor of an auxiliary noise adjustment module (42) is 0.9995. A delay length of a D-order delay operator (31) is 61. The length of the linear prediction filter (32) is 192, and the update step size of the coefficients is 0.5. A controller (21) employs a linear filter with length is 192, and the update step size of the coefficient is 0.040. λ.sub.m, α and T.sub.p are 0.98, 1.8 and 20, respectively. The number of independent operations is 100. The length of simulation data is 120,000.

(49) FIG. 3A is a plot showing the change of target noise and residual noise in the case of actual target noise and an actual secondary path in Example 4. FIG. 3B is a plot showing the change of an auxiliary noise scaling factor in this case. When the system reaches its steady state, the noise reduction of the first half and the second half of the system are 10.55 dB and 12.08 dB, respectively, showing that the system of the present disclosure can not only effectively estimate the actual secondary path in an IIR type, but also has nice performance of suppressing the target noise with large sudden change.

(50) Some of the steps in the examples of the present disclosure may be implemented using software, and the corresponding software program may be stored in a readable storage medium, such as an optical disk or a hard disk.

(51) The above description is only preferred examples of the present disclosure, and is not intended to limit the present disclosure. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present disclosure shall fall within the scope of protection of the present disclosure.

Feedback active noise control system and strategy with online secondary-path modeling

Assignee

Inventors

Cpc classification

Classification Explorer

G10K11/17875

PHYSICS

Classification Explorer

G10K2210/3012

PHYSICS

Classification Explorer

G10K2210/3035

PHYSICS

Classification Explorer

G10K2210/3047

PHYSICS

Classification Explorer

G10K2210/3026

PHYSICS

Classification Explorer

G10K2210/3023

PHYSICS

Classification Explorer

G10K11/17853

PHYSICS

Classification Explorer

G10K11/17855

PHYSICS

Classification Explorer

G10K11/17817

PHYSICS

Classification Explorer

G10K11/17879

PHYSICS

Classification Explorer

G10K2210/121

PHYSICS

Classification Explorer

G10K2210/30232

PHYSICS

Classification Explorer

G10K2210/112

PHYSICS

International classification

Classification Explorer

G10K11/178

PHYSICS

Abstract

Claims

Description