CN105575382A - Complete parallel narrow-band active noise control method with rapid and stable convergence - Google Patents

Abstract

Provided is a complete parallel narrow-band active noise control method with rapid and stable convergence. The invention relates to a complete parallel narrow-band active noise control method with rapid and stable convergence. Problems that the coefficient of a main controller of each conventional frequency channel is still updated by overall residual errors of a system, the controllers of the channels are mutually affected, and the convergence speed of the system is low are solved. The specific process includes: 1) a primary noise signal p (n) with the sinusoidal characteristic is established; 2) a secondary noise signal y (n) with the equal amplitude and opposite phase of the primary noise signal p (n) in step 1) is synthesized; 3) after the primary noise signal p (n) in step 1) and the secondary noise signal y (n) in step 2) are cancelled, a general residual error signal e (n) of the system is detected via a microphone; and 4) according to the general residual error signal e (n) of the system detected by step 3), a residual error signal ei (n) of the ith frequency channel separated from a residual error separating subsystem is obtained. The complete parallel narrow-band active noise control method is applied to the field of noise processing.

Description

A kind of complete parallel arrowband active noise controlling method of fast and stable convergence

Technical field

The present invention relates to the complete parallel arrowband active noise controlling method of fast and stable convergence.

Background technology

In the productive life of reality, there is a large amount of hazardous noises, such as, the noise that the medium-and-large-sized cutting machine of factory produces, can cause very large harm to operator, is reduced or to eliminate these noise signals significant by effective measures.From early 1970s, active noise controlling (ANC) system is applied to reducing these hazardous noises.ANC is the destructive interference principle based on sound, produces the secondary noise signal identical with primary noise frequency, amplitude is close, phase place is contrary by loudspeaker, thus effectively reduces or eliminate primary noise.

The hazardous noise based on low frequency is produced by rotating machineries such as engine, cutting machine, exhaust fans, it has periodically or nearly periodicity, due in arrowband active noise control system, usually non-acoustic sensor (as accelerometer) is utilized to obtain this kind of harmonic noise signals frequency, and then produce reference signal by signal generator, thus avoid acoustic feedback problem, therefore, in this kind of low-frequency noise of elimination, traditional arrowband active noise control system obtains research and apply widely.

For multi-frequency narrowband noise signals, traditional narrow active noise control system is that multiple second order self-adaptive filter parallel is connected to form multi-channel structure, achieve each narrow band noise component by the independent processing that walks abreast, and the master controller coefficient of each frequency channel is still upgraded by the overall residual error of system, cause influencing each other between each channel controller, the speed of convergence of system is low.

Summary of the invention

The master controller coefficient that the object of the invention is to solve existing each frequency channel is but still upgraded by the overall residual error of system, cause influencing each other between each channel controller, the problem that the speed of convergence of system is low, and the complete parallel arrowband active noise controlling method proposing the convergence of a kind of fast and stable.

Above-mentioned goal of the invention is achieved through the following technical solutions:

Step one, foundation have primary noise signal p (n) of sinuso sine protractor;

Step 2, synthesis and secondary noise signal y (n) that primary noise signal p (n) amplitude is equal, phase place is contrary in step one;

Step 3, after primary noise signal p (n) in step one disappears mutually with secondary noise signal y (n) in step 2, record overall system residual error signal e (n) by microphone;

Step 4, according to overall system residual error signal e (n) recorded in step 3, draw the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem _i(n).

Invention effect

Adopt the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence of the present invention, in system of the present invention, comprise residual error isolated subsystem, achieve and overall system residual error is separated according to frequency self-adaption, make the renewal of each channel controller all from the residual error that himself frequency is corresponding, and no longer by the impact of other frequency residuals errors, constitute the arrowband active noise control system of input frequency signal and residual error frequency signal complete parallel structure one to one, in guarantee steady-state behaviour situation, improve the speed of convergence of system, solve in existing arrowband active noise control system, each narrow band noise component is by the independent processing that walks abreast, and each controller coefficient is still upgraded by the overall residual error of system, cause influencing each other between each channel controller, the problem that the speed of convergence of system is low, the speed of convergence of system is made to improve nearly 1 times.

Accompanying drawing explanation

Fig. 1 is process flow diagram of the present invention;

Fig. 2 is Traditional parallel arrowband active noise control system figure, cos (ω _in) be cosine signal, sin (ω _in) be sinusoidal signal, for cosine reference input, for sinusoidal reference input, be the coefficient of i-th cosine component, be the coefficient of i-th sinusoidal component, ∑ is all summations, y _in () is that the i-th tunnel exports, y ₁n () is that the 1st tunnel exports, y _qn () is that q road exports, y (n) is total output ,+for being added, S (z) is secondary channel, y _pn () is the output through primary channel, p (n) is primary channel, and e (n) is overall residual error ,-for subtracting each other, for the estimation of secondary channel, for cosine reference input after filtering, for sinusoidal reference input after filtering, FXLMS is the least mean square algorithm of filtering X passage;

Fig. 3 is residual error isolated subsystem figure, e _in () is the output after separation, for the coefficient of i-th cosine component of isolated subsystem, for the coefficient of i-th sinusoidal component of isolated subsystem, LMS is least mean square algorithm;

Fig. 4 is the arrowband active noise control system structural drawing of the complete parallel of band residual error isolated subsystem,

Fig. 5 is arrowband active noise control system and the legacy system convergence process comparison diagram of complete parallel.

Embodiment

Embodiment one: composition graphs 1, Fig. 2, Fig. 3, Fig. 4 and Fig. 5 illustrate present embodiment, the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence of present embodiment, specifically prepare according to following steps:

Step one, foundation have primary noise signal p (n) of sinuso sine protractor;

Step 2, synthesis and secondary noise signal y (n) that primary noise signal p (n) amplitude is equal, phase place is contrary in step one;

Step 3, after primary noise signal p (n) in step one disappears mutually with secondary noise signal y (n) in step 2, record overall system residual error signal e (n) by microphone;

Step 4, according to overall system residual error signal e (n) recorded in step 3, draw the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem _i(n).

Embodiment two: present embodiment and embodiment one unlike: described p (n) initial value is 0, e (n) initial value is 0, e _in () initial value is 0, the equal assignment 0 of initial value of all the other all variablees.

Other step and parameter identical with embodiment one.

Embodiment three: present embodiment and embodiment one or two unlike: set up primary noise signal p (n) with sinuso sine protractor in described step one, detailed process is:

$p\left(n\right)=\underset{i=1}{\overset{q}{\Σ}}\[a_{i}cos\left({\mathrm{\ω}}_{i}n\right)+b_{i}sin\left({\mathrm{\ω}}_{i}n\right)\]+v_p\left(n\right)$

Wherein, ω _ibe i-th frequency, i is frequency number, and i span is 1≤i≤q; Q is frequency number, and the span of q is positive integer; for the discrete fourier coefficient (DFC) of primary noise signal p (n), a _ibe the coefficient of i-th cosine component, b _ibe the coefficient of i-th sinusoidal component, v _p(n) for average be zero, variance is white Gaussian noise, n be algorithm iteration run point, n represents the time point number of interative computation, value be more than or equal to 1 any positive integer; for representing variance value size.

Other step and parameter identical with embodiment one or two.

Embodiment four: present embodiment and one of embodiment one to three unlike: synthesize in described step 2 and secondary noise signal y (n) that primary noise signal p (n) amplitude is equal, phase place is contrary in step one, detailed process is:

First, i-th frequencies omega of primary noise signal is obtained by accelerometer non-acoustic sensor _i(i=1,2 ... q),

Then, i-th frequencies omega of primary noise signal p (n) _i(i=1,2 ... q) sine and cosine reference signal is produced through signal generator wherein, for sinusoidal reference signal, for cosine reference signal, T is the cycle;

By the discrete fourier coefficient of secondary noise signal y (n) synthon system adjustment sine and cosine reference signal amplitude and phase place, synthesize secondary noise signal y (n) equal with primary noise signal p (n) amplitude, phase place is contrary, (sound wave has mutual interference effect in transmitting procedure to utilize the destructive interference principle of sound wave, two frequency is identical, direction of vibration is identical and the sound wave that the sound source of acting in agreement sends superposes mutually time just there will be interference, if their phase place is identical, two ripple superposition amplitude increase acoustic pressure and strengthen; Otherwise their phase place is contrary, two ripple superposition amplitude reduce acoustic pressure and weaken, if two wave amplitude are the same, will offset completely) eliminate primary noise signal p (n).

Other step and parameter identical with one of embodiment one to three.

Embodiment five: one of present embodiment and embodiment one to four unlike: described secondary noise signal y (n) is provided by following formula:

$\begin{array}{c}y\left(n\right)=\underset{i=1}{\overset{q}{\Σ}}{y}_i\left(n\right)\\ =\underset{i=1}{\overset{q}{\Σ}}{\hat{a}}_i\left(n\right)x_{a_i}\left(n\right)+{\hat{b}}_i\left(n\right)x_{b_i}\left(n\right)\end{array}$

In formula, $x_{a_i}\left(n\right)=cos\left({\mathrm{\ω}}_{i}n\right),x_{b_i}\left(n\right)=sin\left({\mathrm{\ω}}_{i}n\right),{\{{\hat{a}}_i\left(n\right),{\hat{b}}_i\left(n\right)\}}_{i=1}^q$ For the discrete fourier coefficient of secondary noise signal y (n) synthon system, be the coefficient of the cosine component of the n-th time point of i-th frequency, be the coefficient of the sinusoidal component of the n-th time point of i-th frequency, q is frequency number, and q span is positive integer, y _in () is the i-th tunnel output signal, span is that the noise to be offset with target is relevant.

Other step and parameter identical with one of embodiment one to four.

Embodiment six: one of present embodiment and embodiment one to five are unlike the discrete fourier coefficient of described secondary noise signal y (n) synthon system detailed process is:

Residual error isolated subsystem adopts least mean square algorithm (Filter-XLeastMeanSquare, FXLMS) to receive the reference signal after secondary channel filters and the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem in step 4 _i(n); The discrete fourier coefficient of secondary noise signal y (n) synthon system is upgraded again by least mean square algorithm module the more new formula of described discrete fourier coefficient is:

${\hat{a}}_i(n+1)={\hat{a}}_i\left(n\right)+{\mathrm{\μ}}_i{e}_i\left(n\right){\hat{x}}_{a_i}\left(n\right);$

${\hat{b}}_i(n+1)={\hat{b}}_i\left(n\right)+{\mathrm{\μ}}_i{e}_i\left(n\right){\hat{x}}_{b_i}\left(n\right);$

In formula, μ _ifor the step-length of secondary noise signal syntheses subsystem 1;

Sinusoidal reference signal with cosine reference signal estimate through secondary channel filtering obtains filtered reference signal:

${\hat{x}}_{a_i}\left(n\right)=\underset{m=0}{\overset{\hat{M}-1}{\Σ}}{\hat{s}}_m{x}_{a_i}(n-m);$

${\hat{x}}_{b_i}\left(n\right)=\underset{m=0}{\overset{\hat{M}-1}{\Σ}}{\hat{s}}_m{x}_{b_i}(n-m);$

In formula, for secondary channel model coefficient, for secondary channel model order, m is the sequence number of secondary channel, and span is 0 to arrive with for the reference-input signal of pair of orthogonal, with for the reference signal pair after secondary channel filters, it is the coefficient of m secondary channel.

Other step and parameter identical with one of embodiment one to four.

Embodiment seven: present embodiment and one of embodiment one to six unlike: in described step 3 after secondary noise signal y (n) produced in primary noise signal p (n) in step one and step 2 disappears mutually, record overall system residual error signal e (n) by microphone; Detailed process is:

$\begin{array}{c}e\left(n\right)=p\left(n\right)-y_p\left(n\right)\\ =p\left(n\right)-\underset{j=0}{\overset{M-1}{\Σ}}{s}_{j}y(n-j)\end{array}$

Wherein, for the coefficient of secondary channel S (z), M-1 is the exponent number of S (z), and M is positive integer; y _pn () is the output signal of y (n) after primary channel, y (n) through primary channel be y (n) through air borne, through electronic equipment; P is the english abbreviation of primary, and j is a jth filter loop variable, the span of n to be 1≤n≤N, N be more than or equal to 1 any positive integer, the span of j is that the length of secondary channel is from 0 to M-1, s _jfor a jth secondary channel coefficient.

Other step and parameter identical with one of embodiment one to six.

Embodiment eight: one of present embodiment and embodiment one to seven unlike: according to overall system residual error signal e (n) recorded in step 3 in described step 4, draw the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem 2 _i(n), detailed process is:

Containing q frequency component in overall system residual error signal e (n) recorded in step 3, utilize sine and cosine reference signal by the discrete fourier coefficient of residual error isolated subsystem 2 adjusting amplitude and phase place, thus the residual error of each frequency separated from overall system residual error e (n), is the discrete fourier coefficient of secondary noise signal y (n) synthon system renewal the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem 2 is provided _i(n), it is calculated by following formula:

$e_i\left(n\right)={\hat{a}}_{e,i}\left(n\right)x_{a_i}\left(n\right)+{\hat{b}}_{e,i}\left(n\right)x_{b_i}\left(n\right)$

Wherein, for i-th cosine fourier coefficient component of error separate device output subsystem, for i-th sinusoidal Fourier coefficients component of error separate device output subsystem, for sinusoidal reference signal, for cosine reference signal.

Other step and parameter identical with one of embodiment one to seven.

Embodiment nine: one of present embodiment and embodiment one to eight unlike: described by residual error signal e _in () brings step 2 into, draw step 4, until complete N time of circular flow, N be more than or equal to 1 any positive integer, N is artificial setting.

Other step and parameter identical with one of embodiment one to eight.

Embodiment ten: one of present embodiment and embodiment one to nine are unlike the discrete fourier coefficient of described residual error isolated subsystem 2 adjust amplitude and phase place, detailed process is:

By leak least mean square algorithm (LeakyLeastMeanSquare, LLMS) module carry out real-time Adjustable calculation (namely in cost function, introduce leakage factor α, its cost function is specially $J_e\left(n\right)=\frac{1}{2}{e}^2\left(n\right)+\frac{1}{2}\underset{j=1}{\overset{q}{\Σ}}{\mathrm{\α}}_j\[{\hat{a}}_{e_j}^2\left(n\right)+{\hat{b}}_{e_j}^2\left(n\right)\]),$ Concrete more new formula is:

${\hat{a}}_{e,i}(n+1)=(1-{\mathrm{\μ}}_e{\mathrm{\α}}_i){\hat{a}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{a_i}\left(n\right)x_{a_i}\left(n\right)=\mathrm{\ξ}{\hat{a}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{a_i}\left(n\right);$

${\hat{b}}_{e,i}(n+1)=(1-{\mathrm{\μ}}_e{\mathrm{\α}}_i){\hat{b}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{b_i}\left(n\right)=\mathrm{\ξ}{\hat{b}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{b_i}\left(n\right);$

In formula, e (n) is overall residual error, for sinusoidal reference signal, for cosine reference signal, α _ifor leakage factor, μ _efor the step-length of residual error isolated subsystem, ξ=1-μ _eα _i, its span is (0,1), and ξ is system leak coefficient.

Other step and parameter identical with one of embodiment one to nine.

Embodiment 1

In order to investigate the performance of native system, we and traditional arrowband active noise control system carry out simulation comparison, here representative simulation result is only provided, as shown below, give when leakage factor α=0.12, the all square residual error simulation comparison curve (green line represents traditional arrowband active noise control system, and red line represents native system) of two systems.Under the prerequisite ensureing identical steady-state error, can find that native system is 500 everywhere convergents, and legacy system is 1000 convergences later, new system improves nearly 1 times compared with the speed of convergence of legacy system.

Claims (10)

1. a complete parallel arrowband active noise controlling method for fast and stable convergence, is characterized in that, the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence specifically carries out work according to following steps:

Step one, foundation have primary noise signal p (n) of sinuso sine protractor;

Step 2, synthesis and secondary noise signal y (n) that primary noise signal p (n) amplitude is equal, phase place is contrary in step one;

Step 3, after primary noise signal p (n) in step one disappears mutually with secondary noise signal y (n) in step 2, record overall system residual error signal e (n) by microphone;

Step 4, according to overall system residual error signal e (n) recorded in step 3, draw the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem _i(n).

2. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 1, it is characterized in that, described p (n) initial value is 0, e (n) initial value is 0, e _in () initial value is 0, the equal assignment 0 of initial value of all the other all variablees.

3. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 2, it is characterized in that, set up primary noise signal p (n) with sinuso sine protractor in described step one, detailed process is:

$p\left(n\right)=\underset{i=1}{\overset{q}{\Σ}}\[a_{i}cos\left({\mathrm{\ω}}_{i}n\right)+b_{i}sin\left({\mathrm{\ω}}_{i}n\right)\]+v_p\left(n\right)$

Wherein, ω _ibe i-th frequency, i is frequency number, and i span is 1≤i≤q; Q is frequency number, and the span of q is positive integer; for the discrete fourier coefficient of primary noise signal p (n), a _ibe the coefficient of i-th cosine component, b _ibe the coefficient of i-th sinusoidal component, v _p(n) for average be zero, variance is white Gaussian noise, n be algorithm iteration run point, n represents the time point number of interative computation, value be more than or equal to 1 any positive integer; for representing variance value size.

4. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 3, it is characterized in that, synthesis and secondary noise signal y (n) that primary noise signal p (n) amplitude is equal, phase place is contrary in step one in described step 2, detailed process is:

First, i-th frequencies omega of primary noise signal is obtained by accelerometer non-acoustic sensor _i(i=1,2...q),

Then, i-th frequencies omega of primary noise signal p (n) _isine and cosine reference signal is produced through signal generator ${\left[\begin{array}{cc}{x}_{a_i}\left(n\right)& x_{b_i}\left(n\right)\end{array}\right]}^T,$ Wherein, for sinusoidal reference signal, for cosine reference signal, T is the cycle;

By the discrete fourier coefficient of secondary noise signal y (n) synthon system adjustment sine and cosine reference signal ${\left[\begin{array}{cc}{x}_{a_i}\left(n\right)& x_{b_i}\left(n\right)\end{array}\right]}^T$ Amplitude and phase place, synthesize secondary noise signal y (n) equal with primary noise signal p (n) amplitude, phase place is contrary, utilize sound wave destructive interference principle eliminate primary noise signal p (n).

5. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 4, it is characterized in that, described secondary noise signal y (n) is provided by following formula:

$\begin{array}{c}y\left(n\right)=\underset{i=1}{\overset{q}{\Σ}}{y}_i\left(n\right)\\ =\underset{i=1}{\overset{q}{\Σ}}{\hat{a}}_i\left(n\right)x_{a_i}\left(n\right)+{\hat{b}}_i\left(n\right)x_{b_i}\left(n\right)\end{array}$

In formula, $x_{a_i}\left(n\right)=cos\left({\mathrm{\ω}}_{i}n\right),$ $x_{b_i}\left(n\right)=sin\left({\mathrm{\ω}}_{i}n\right),$ ${\{{\hat{a}}_i\left(n\right),{\hat{b}}_i\left(n\right)\}}_{i=1}^q$ For the discrete fourier coefficient of secondary noise signal y (n) synthon system, be the coefficient of the cosine component of the n-th time point of i-th frequency, be the coefficient of the sinusoidal component of the n-th time point of i-th frequency, q is frequency number, and q span is positive integer, y _in () is the i-th tunnel output signal.

6. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 5, is characterized in that, the discrete fourier coefficient of described secondary noise signal y (n) synthon system detailed process is:

Residual error isolated subsystem adopts least mean square algorithm to receive the reference signal after secondary channel filters ${\left[\begin{array}{cc}{\hat{x}}_{a_i}\left(n\right)& {\hat{x}}_{b_i}\left(n\right)\end{array}\right]}^T,$ And the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem in step 4 _i(n); The discrete fourier coefficient of secondary noise signal y (n) synthon system is upgraded again by least mean square algorithm module the more new formula of described discrete fourier coefficient is:

${\hat{a}}_i(n+1)={\hat{a}}_i\left(n\right)+{\mathrm{\μ}}_i{e}_i\left(n\right){\hat{x}}_{a_i}\left(n\right);$

${\hat{b}}_i(n+1)={\hat{b}}_i\left(n\right)+{\mathrm{\μ}}_i{e}_i\left(n\right){\hat{x}}_{b_i}\left(n\right);$

In formula, μ _ifor the step-length of secondary noise signal syntheses subsystem 1;

Sinusoidal reference signal with cosine reference signal estimate through secondary channel filtering obtains filtered reference signal:

${\hat{x}}_{a_i}\left(n\right)=\underset{m=0}{\overset{\hat{M}-1}{\Σ}}{\hat{s}}_m{x}_{a_i}(n-m);$

${\hat{x}}_{b_i}\left(n\right)=\underset{m=0}{\overset{\hat{M}-1}{\Σ}}{\hat{s}}_m{x}_{b_i}(n-m);$

In formula, for secondary channel model coefficient, for secondary channel model order, m is the sequence number of secondary channel, and span is 0 to arrive with for the reference-input signal of pair of orthogonal, with for the reference signal pair after secondary channel filters, it is the coefficient of m secondary channel.

7. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 6, it is characterized in that, in described step 3 after secondary noise signal y (n) produced in primary noise signal p (n) in step one with step 2 disappears mutually, record overall system residual error signal e (n) by microphone; Detailed process is:

$\begin{array}{c}e\left(n\right)=p\left(n\right)-y_p\left(n\right)\\ =p\left(n\right)-\underset{j=0}{\overset{M-1}{\Σ}}{s}_{j}y(n-j)\end{array}$

Wherein, for the coefficient of secondary channel S (z), M-1 is the exponent number of S (z), and M is positive integer; y _pn () is the output signal of y (n) after primary channel, j is a jth filter loop variable, the span of n to be 1≤n≤N, N be more than or equal to 1 any positive integer, the span of j is that the length of secondary channel is from 0 to M-1, s _jfor a jth secondary channel coefficient.

8. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 7, it is characterized in that, according to overall system residual error signal e (n) recorded in step 3 in described step 4, draw the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem _i(n), detailed process is:

Containing q frequency component in overall system residual error signal e (n) recorded in step 3, utilize sine and cosine reference signal ${\left[\begin{array}{cc}{x}_{a_i}\left(n\right),& x_{b_i}\left(n\right)\end{array}\right]}^T,$ By the discrete fourier coefficient of residual error isolated subsystem adjusting amplitude and phase place, thus the residual error of each frequency separated from overall system residual error e (n), is the discrete fourier coefficient of secondary noise signal y (n) synthon system renewal the residual error signal e of isolated i-th frequency channel of residual error isolated subsystem is provided _i(n), it is calculated by following formula:

$e_i\left(n\right)={\hat{a}}_{e,i}\left(n\right)x_{a_i}\left(n\right)+{\hat{b}}_{e,i}\left(n\right)x_{b_i}\left(n\right)$

Wherein, for i-th cosine fourier coefficient component of error separate device output subsystem, for i-th sinusoidal Fourier coefficients component of error separate device output subsystem, for sinusoidal reference signal, for cosine reference signal.

9. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 8, is characterized in that, described by residual error signal e _in () brings step 2 into, draw step 4, until complete N time of circular flow, N be more than or equal to 1 any positive integer, N is artificial setting.

10. the complete parallel arrowband active noise controlling method of a kind of fast and stable convergence according to claim 9, is characterized in that, the discrete fourier coefficient of described residual error isolated subsystem adjust amplitude and phase place, detailed process is:

Carry out real-time Adjustable calculation by leaking least mean square algorithm module, specifically more new formula is:

${\hat{a}}_{e,i}(n+1)=(1-{\mathrm{\μ}}_e{\mathrm{\α}}_i){\hat{a}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{a_i}\left(n\right)=\mathrm{\ξ}{\hat{a}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{a_i}\left(n\right);$

${\hat{b}}_{e,i}(n+1)=(1-{\mathrm{\μ}}_e{\mathrm{\α}}_i){\hat{b}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{b_i}\left(n\right)=\mathrm{\ξ}{\hat{b}}_{e,i}\left(n\right)+{\mathrm{\μ}}_{e}e\left(n\right)x_{b_i}\left(n\right);$

In formula, e (n) is overall residual error, for sinusoidal reference signal, for cosine reference signal, α _ifor leakage factor, μ _efor the step-length of residual error isolated subsystem, ξ=1-μ _eα _i, its span is (0,1), and ξ is system leak coefficient.