CN104063610A - Simulation noise source and original noise sensor transmitting model identifying method - Google Patents

Simulation noise source and original noise sensor transmitting model identifying method Download PDF

Info

Publication number
CN104063610A
CN104063610A CN201410311149.1A CN201410311149A CN104063610A CN 104063610 A CN104063610 A CN 104063610A CN 201410311149 A CN201410311149 A CN 201410311149A CN 104063610 A CN104063610 A CN 104063610A
Authority
CN
China
Prior art keywords
noise
original
filter
white noise
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201410311149.1A
Other languages
Chinese (zh)
Inventor
张华军
赵菁
蔡炜
周登科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wisdri Engineering and Research Incorporation Ltd
Original Assignee
Wisdri Engineering and Research Incorporation Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wisdri Engineering and Research Incorporation Ltd filed Critical Wisdri Engineering and Research Incorporation Ltd
Priority to CN201410311149.1A priority Critical patent/CN104063610A/en
Publication of CN104063610A publication Critical patent/CN104063610A/en
Pending legal-status Critical Current

Links

Landscapes

  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

The invention provides a simulation noise source and original noise sensor transmitting model identifying method. The method comprises the following steps that white noise is simulated and generated according to a standard white noise generation method, the white noise passes through a loud speaker to spread in air, and an original noise sensor detects the signal that the sound wave of the loud speaker passes the air and attenuates and the signal that the white noise passes a filter and filters; assuming that original values of parameters of the filter are zero, the signal that the white noise passes the filter and filters is calculated, and the error between the signal detected by the original noise sensor and the signal that the white noise passes the filter and filters is calculated; the efficient of the filter is updated, and the above steps are repeated; when the error satisfies the requirement, the calculation is stopped, and the efficient of the filter is output. The method can obtain a more accurate noise transmitting model, is suitable for identifying the transmitting model between the loud speaker and the original noise sensor during the active noise controlling design, and provides the basis for actively controlling the noise and obtaining good de-noising effects.

Description

Analogue noise source and original noise transducer TRANSFER MODEL discrimination method
Technical field
The present invention relates to iron and steel metallurgical industry Noise Control field, be specifically related to analogue noise source and original noise transducer TRANSFER MODEL discrimination method.
Background technology
In the industrial circle such as iron and steel, metallurgy, due to the operation of various big machineries, particularly operation high-energy noise sources such as pressure reducing valve of blast furnace group, causes the head of district of the steel plant phase among serious noise pollution.Owing to can causing staff's health and psychology to sustain damage for a long time, therefore in engineering, often adopt initiatively denoising and two kinds of method noise decrease energy of passive denoising in this high intensity noise pollutes.In active noise control technique, because original noise transducer has not only received original noise, it has also received the analogue noise signal that analogue noise loudspeaker sends simultaneously, in order to reach good denoising effect, need identification loudspeaker to original noise transducer noise transfer model before, classic method is often difficult to identification and obtains noise transfer model accurately.
Summary of the invention
The technical problem to be solved in the present invention is: a kind of analogue noise source and original noise transducer TRANSFER MODEL discrimination method are provided, can obtain noise transfer model more accurately.
The present invention for solving the problems of the technologies described above taked technical scheme is: a kind of analogue noise source and original noise transducer TRANSFER MODEL discrimination method, is characterized in that: it comprises the following steps:
1) making the sampling period is T s, generate white noise y (k) according to the white noise production method simulation of standard, y (k) to be propagated in air through loudspeaker, original noise transducer detects the signal x ' of loudspeaker sound wave after attenuation of air (k); White noise y (k) is the filtered signal x of device (k)=D (z) y (k) after filtering, D (z) represents the transport function between analogue noise source and original noise transducer, and expression formula is D (z)=d 0+ d 1z -1+ d 2z -2+ ... + d nz -n, wherein z represents shifting function, the order that n is wave filter, d 0, d 1..., d nfor filter coefficient;
2) suppose that the initial value of filter parameter is 0;
3) calculate the white noise filtered signal x of device (k) after filtering, wherein d i(k) represent filter D (z) coefficient d ivalue in the time of the k time sampling;
4) calculate signal x ' that original noise transducer detects (k) and the white noise error e (k) between the filtered signal x of device (k) after filtering, e (k)=x ' (k)-x (k);
5) judge the size of e (k):
In the time that e (k) does not reach preset requirement, upgrade filter coefficient, make d i(k)=d i(k-1)+ε × e (k) × y (k-i), wherein i=0,1,2 ..., n-1, ε is iteration step length, ε satisfies condition wherein P yfor the power of white noise y (k); Repeat above-mentioned steps 3)-4);
In the time that reaching preset requirement, error stops calculating output filter coefficient.
Press such scheme, step 5) in, when e (k)>=0.1, upgrade filter coefficient, make d i(k)=d i(k-1)+ε × e (k) × y (k-i), wherein i=0,1,2 ..., n-1, ε is iteration step length, ε satisfies condition wherein P yfor the power of white noise y (k); Repeat above-mentioned steps 3)-4);
In the time of e (k) <0.1, stop calculating output filter coefficient.
Beneficial effect of the present invention is: can obtain noise transfer model more accurately than classic method, while being applicable to very much active noise controlling design, the TRANSFER MODEL between loudspeaker and original noise transducer is carried out to identification, provide the foundation for active noise controlling obtains good denoising effect.
Brief description of the drawings
Fig. 1 is the identification principle figure of one embodiment of the invention.
Fig. 2 is one embodiment of the invention median filter parameter d 0convergence curve.
Fig. 3 is one embodiment of the invention median filter parameter d 1convergence curve.
Fig. 4 is one embodiment of the invention median filter parameter d 2convergence curve.
Fig. 5 is one embodiment of the invention median filter parameter d 3convergence curve.
Fig. 6 is one embodiment of the invention median filter parameter d 4convergence curve.
Fig. 7 is one embodiment of the invention median filter parameter d 5convergence curve.
Fig. 8 is one embodiment of the invention median filter parameter d 6convergence curve.
Fig. 9 is error e in one embodiment of the invention (k) change curve.
Embodiment
Below in conjunction with instantiation and accompanying drawing, the present invention will be further described.
On some directions of propagation, sound is along straight line to propagating at a distance, and first sensor of foremost is original noise transducer, and centre is loudspeaker, the error pick-up of finally placing.Wherein original noise transducer be for check noise source along the direction of propagation noise signal on original noise transducer position, and loudspeaker is the sound wave offseting with noisiness sending by active noise control technique, last error pick-up is the device for detecting original noise sound wave and loudspeaker sound wave, because the sound wave that original signal and loudspeaker send is cancelled each other in air, therefore actual what receive is the acoustic signals that original noise and loudspeaker sound wave superpose mutually to error pick-up, it can check the denoising performance quality of noise active control technology.Because noise is to the so-called first passage of travel path of error pick-up, so the sound wave just loudspeaker being sent is called second channel to the travel path of error pick-up.
Analogue noise provided by the invention source and original noise transducer TRANSFER MODEL discrimination method, comprise the following steps:
1) set up noise transfer passage according to Fig. 1, making the sampling period is T s, generate white noise y (k) according to the white noise production method simulation of standard, y (k) to be propagated in air through loudspeaker, original noise transducer detects the signal x ' of loudspeaker sound wave after attenuation of air (k); White noise y (k) is the filtered signal x of device (k)=D (z) y (k) after filtering, D (z) represents the transport function between analogue noise source and original noise transducer, and expression formula is D (z)=d 0+ d 1z -1+ d 2z -2+ ... + d nz -n, wherein z represents shifting function, the order that n is wave filter, d 0, d 1..., d nfor filter coefficient;
2) suppose that the initial value of filter parameter is 0;
3) calculate the white noise filtered signal x of device (k) after filtering, wherein d i(k) represent filter D (z) coefficient d ivalue in the time of the k time sampling;
4) calculate signal x ' that original noise transducer detects (k) and the white noise error e (k) between the filtered signal x of device (k) after filtering, e (k)=x ' (k)-x (k);
5) judge the size of e (k):
When e (k)>=0.1, upgrade filter coefficient, make d i(k)=d i(k-1)+ε × e (k) × y (k-i), wherein i=0,1,2 ..., n-1, ε is iteration step length, ε satisfies condition wherein P yfor the power of white noise y (k); Repeat above-mentioned steps 3)-4);
In the time of error e (k) <0.1, stop calculating output filter coefficient.
Principle of work of the present invention is: when hypothetical simulation noise is propagated in air through loudspeaker, be H (z) from loudspeaker to the transport function original noise transducer, so in accompanying drawing 1 white noise y (k) after the attenuation of air original noise transducer detect the signal obtaining become x ' (k)=H (z) y (k), and signal x (k)=D (z) y (k), if represent that the error of x ' (k) and between x (k) is at 1 o'clock, D (z) is the transport function between analogue noise source and original noise transducer.The present invention regards the transport function between analogue noise source and original noise transducer as a finite impulse wave filter, and its expression-form is D (z)=d 0+ d 1z -1+ d 2z -2+ ... + d nz -n, wherein symbol " z " represents shifting function, n is filter order, d 0, d 1..., d nfor filter coefficient.According to the form of matrix operation, signal x (k) can be expressed as wherein vector θ=[d 0d 1d 2d n] t.
Make the accurately TRANSFER MODEL between approaching simulation noise source and original noise transducer of D (z), can be by solving following minimization problem:
Wherein N is sampled data summation.
Meet the minimized parameter θ of formula (1) and be according to Optimum Theory is known:
Due to
Therefore formula (2) can be converted to:
Due to matrix inversion operation in the time that unknown parameter is more, calculated amount is very large, therefore can utilize a less constant to replace, and rule of thumb formula (5) can be reduced to following form here:
Wherein constant p yfor the power of white noise y (k).
Due to therefore formula (6) can be reduced to
Again because therefore θ ncan be expressed as:
d 0 d 1 &CenterDot; &CenterDot; &CenterDot; d n N = &theta; N - 1 + &epsiv; &times; e ( N ) &times; y ( N ) y ( N - 1 ) &CenterDot; &CenterDot; &CenterDot; y ( N - n ) - - - ( 8 ) ,
Each parameter calculation formula that can obtain transport function D (z) according to formula (8) is:
d i(N)＝d i(N-1)+ε×e(N)×y(N-i) (9)，
Wherein i=0,1,2 ..., n-1.
Taking a certain active noise controlling system as example, wherein 40 meters of the original noise transducers of loudspeaker distance, for the TRANSFER MODEL between identification loudspeaker and original noise transducer, it is T that present case makes the sampling period s=0.00001 second, according to step 1) generate white noise sequence y (k), propagate in air through loudspeaker.Establishing TRANSFER MODEL according to step is 7 rank, and the initial parameter of designing filter is 0, i.e. d 0=0, d 1=0, d 2=0 ..., d 6=0.According to step 3)-4) iterative computation filter parameter progressively, after the calculating in 0.2 second, error e (k) <0.1, the filter parameter that identification obtains is respectively d 0=0.321862672903827, d 1=0.0657865193590414, d 2=0.0689975523138478, d 3=0.0701804274787997, d 4=0.0689975523138478, d 5=0.0657865193590414, d 6=0.311862672903827, each parameter iteration convergence process curve is as shown in accompanying drawing 2-8, and error e (k) change curve as shown in Figure 9.
Above embodiment is only for calculating thought of the present invention and feature are described, its object is to make those skilled in the art can understand content of the present invention and implement according to this, and protection scope of the present invention is not limited to above-described embodiment.So the disclosed principle of all foundations, equivalent variations or the modification that mentality of designing is done, all within protection scope of the present invention.

Claims (2)

1. analogue noise source and an original noise transducer TRANSFER MODEL discrimination method, is characterized in that: it comprises the following steps:
1) making the sampling period is T s, generate white noise y (k) according to the white noise production method simulation of standard, y (k) to be propagated in air through loudspeaker, original noise transducer detects the signal x ' of loudspeaker sound wave after attenuation of air (k); White noise y (k) is the filtered signal x of device (k)=D (z) y (k) after filtering, D (z) represents the transport function between analogue noise source and original noise transducer, and expression formula is D (z)=d 0+ d 1z -1+ d 2z -2+ ... + d nz -n, wherein z represents shifting function, the order that n is wave filter, d 0, d 1..., d nfor filter coefficient;
2) suppose that the initial value of filter parameter is 0;
3) calculate the white noise filtered signal x of device (k) after filtering, wherein d i(k) represent filter D (z) coefficient d ivalue in the time of the k time sampling;
4) calculate signal x ' that original noise transducer detects (k) and the white noise error e (k) between the filtered signal x of device (k) after filtering, e (k)=x ' (k)-x (k);
5) judge the size of e (k):
In the time that e (k) does not reach preset requirement, upgrade filter coefficient, make d i(k)=d i(k-1)+ε × e (k) × y (k-i), wherein i=0,1,2 ..., n-1, ε is iteration step length, ε satisfies condition wherein P yfor the power of white noise y (k); Repeat above-mentioned steps 3)-4);
In the time that reaching preset requirement, error stops calculating output filter coefficient.
2. analogue noise according to claim 1 source and original noise transducer TRANSFER MODEL discrimination method, is characterized in that: step 5) in, when e (k)>=0.1, upgrade filter coefficient, make d i(k)=d i(k-1)+ε × e (k) × y (k-i), wherein i=0,1,2 ..., n-1, ε is iteration step length, ε satisfies condition wherein P yfor the power of white noise y (k); Repeat above-mentioned steps 3)-4);
In the time of e (k) <0.1, stop calculating output filter coefficient.
CN201410311149.1A 2014-07-01 2014-07-01 Simulation noise source and original noise sensor transmitting model identifying method Pending CN104063610A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410311149.1A CN104063610A (en) 2014-07-01 2014-07-01 Simulation noise source and original noise sensor transmitting model identifying method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410311149.1A CN104063610A (en) 2014-07-01 2014-07-01 Simulation noise source and original noise sensor transmitting model identifying method

Publications (1)

Publication Number Publication Date
CN104063610A true CN104063610A (en) 2014-09-24

Family

ID=51551320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410311149.1A Pending CN104063610A (en) 2014-07-01 2014-07-01 Simulation noise source and original noise sensor transmitting model identifying method

Country Status (1)

Country Link
CN (1) CN104063610A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105848056A (en) * 2016-04-01 2016-08-10 张俊斌 Novel active noise reduction method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0568282A2 (en) * 1992-04-28 1993-11-03 Westinghouse Electric Corporation Noise control system
JPH06259085A (en) * 1993-03-09 1994-09-16 Fujitsu Ltd Estimating method for transfer characteristic of active noise control system
JPH07129250A (en) * 1993-10-28 1995-05-19 Nissan Motor Co Ltd Passive vibration controller and passive noise controller
CN1573917A (en) * 2003-06-17 2005-02-02 本田技研工业株式会社 Active vibratory noise control apparatus
CN1651726A (en) * 2005-02-18 2005-08-10 中国船舶重工集团公司第七一一研究所 Adaptive active control silencer
CN101354885A (en) * 2007-01-16 2009-01-28 哈曼贝克自动系统股份有限公司 Active noise control system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0568282A2 (en) * 1992-04-28 1993-11-03 Westinghouse Electric Corporation Noise control system
JPH06259085A (en) * 1993-03-09 1994-09-16 Fujitsu Ltd Estimating method for transfer characteristic of active noise control system
JPH07129250A (en) * 1993-10-28 1995-05-19 Nissan Motor Co Ltd Passive vibration controller and passive noise controller
CN1573917A (en) * 2003-06-17 2005-02-02 本田技研工业株式会社 Active vibratory noise control apparatus
CN1651726A (en) * 2005-02-18 2005-08-10 中国船舶重工集团公司第七一一研究所 Adaptive active control silencer
CN101354885A (en) * 2007-01-16 2009-01-28 哈曼贝克自动系统股份有限公司 Active noise control system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SEN M. KUO 等: "Active Noise Control: A Tutorial Review", 《PROCEEDING OF THE IEEE》 *
SEN M.KUO ET AL: "Design of active noise control systems with the TMS320 family", 《TEXAS INSTRUMENTS》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105848056A (en) * 2016-04-01 2016-08-10 张俊斌 Novel active noise reduction method

Similar Documents

Publication Publication Date Title
CN102843633B (en) The control of loudspeaker output
CN106068007A (en) Nonlinear loudspeaker system identifying method
CN100578622C (en) A kind of adaptive microphone array system and audio signal processing method thereof
CN105659080B (en) System and method for crossing noise filtering from from the acoustic energy of valve
CN103685111A (en) Calculating method of digital pre-distortion parameters and pre-distortion system
CN105916079A (en) Loudspeaker nonlinearity compensation method and apparatus thereof
CN106531145B (en) Recurrence active noise control method based on M estimator
CN103105773A (en) Sound parametric array control method based on neural network opposite identification and self-adaptive piping and instrument diagram (PID)
CN112233644A (en) Filtering-X least mean square active noise control method based on quaternion adaptive filter
CN103632009B (en) A kind of method for designing of analog feedback active noise reduction earphone
CN109379652B (en) Earphone active noise control secondary channel off-line identification method
CN109556797A (en) The pipeline leakage detection and location method with convolutional neural networks is decomposed based on spline local mean value
CN104864274A (en) Enhancement method for pipeline leakage acoustic signals
CN104716928B (en) A kind of digital filter of online zero phase-shift iir digital filter
CN103399589A (en) Random vibration control method for electro-hydraulic acceleration servo system
Shi et al. Prediction method of ball valve internal leakage rate based on acoustic emission technology
CN112865721B (en) Signal processing method, device, equipment, storage medium, chip and module equipment
CN104063610A (en) Simulation noise source and original noise sensor transmitting model identifying method
CN107241070B (en) Method for determining DPD coefficient, FPGA and DPD processing system
CN104123438A (en) Method for recognizing second noise transmission channel model
TW201414991A (en) Method for measuring electroacoustic parameters of transducer
Liu et al. Fuzzy pattern recognition of impact acoustic signals for nondestructive evaluation
CN106782594A (en) One kind suppresses the unmatched method of frequency in arrowband active noise controlling
Akingbade et al. Separation of digital audio signals using least-mean-square (LMS) adaptive algorithm
CN208025008U (en) A kind of duct noise active control device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20140924

RJ01 Rejection of invention patent application after publication