CN108229030A - A kind of design method of active noise reduction system controller parameter - Google Patents
A kind of design method of active noise reduction system controller parameter Download PDFInfo
- Publication number
- CN108229030A CN108229030A CN201810010686.0A CN201810010686A CN108229030A CN 108229030 A CN108229030 A CN 108229030A CN 201810010686 A CN201810010686 A CN 201810010686A CN 108229030 A CN108229030 A CN 108229030A
- Authority
- CN
- China
- Prior art keywords
- noise reduction
- space
- probability
- energy
- noise
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
The present invention has developed a kind of design method of active noise reduction system controller parameter, the probability for going out each position energy according to the probability calculation of head part's movement each position is averaged, the space gross energy in noise reduction region is obtained with reference to the probability that each position head part for measuring statistics occurs and the integral energy of each position acquired, in the case of this integral energy minimum, calculate design of the controller parameter for controller in practical application, so as to make the noise reduction of system optimal, it can be understood as in the position noise reduction for needing most noise reduction, to realize the optimal noise reduction of overall space, it is not needing to add additional noise, it does not need in the case that sensor is set in human ear road, the noise level in space near human ear is reduced well.
Description
Technical field
The present invention relates to noise reduction field more particularly to a kind of design methods of active noise reduction system controller parameter.
Background technology
When designing active noise reduction system for the spatial field near human ear, generally require and first use the mold number of people or the person
Carry out test early period, placed specifically at the duct of mold sensor and measurement noise field human ear position characteristic with
And transmission function of the secondary sound source at human ear, these measurement data can be stored in signal processor simultaneously after testing early period
The systematic parameter of active noise reduction controller is designed according to these measurement results, but when system is applied to true environment, noise reduction
Target is the human ear position of user in actual use, however the head position of user is uncertain, meeting in use
There is a certain range of random movement near the mold number of people placement location of test early period, this can lead to swashing for sub-loudspeaker
It encourages the voice signal that source generates in practical human ear position and the data obtained in being tested in early period is not consistent, therefore according to early period
The controller noise reduction in actual use of test data design can be affected, then in practical applications can not possibly be
Sensor is placed at user's human ear and measures transmission function of the secondary sound source in practical human ear position in real time.Even if use can be allowed
Family goes out wearable sensors in ear, and this method measured in real time can not reach ideal effect, because measuring loud speaker
In addition transmission function needs send out an additional noise pumping signal, then by ear while active noise reduction with loud speaker
The corresponding response signal of the additional excitation is extracted in sensor could complete real-time testing, but this additional excitation sound source can make
Active noise reduction effect is deteriorated.
Invention content
The present invention is intended to provide a kind of controller parameter design method of active noise reduction system, can not need to plus add make an uproar
Sound is not needed in the case of setting sensor at human ear, and for the random fluctuation of user's head position, to reduce human ear well attached
The noise level of near space.
To realize the above-mentioned technical purpose, the present invention is using following technical scheme, a kind of active noise reduction system controller parameter
Design method, include the following steps：
S1, it determines noise reduction region, multiple incoming devices and output device is installed in the space for needing noise reduction；
S2, the probability that statistics head part appears in each random site is measured；
The probability of each position energy that S3, calculating head part occur is averaged, due to having measured statistics in previous experiments
Next probability existing for the head part of spatial points position needs the energy of measurement space each point, within this space error wheat
The quantity of gram wind may be 2,4 or n, depending on the feature with specific space：
α represents some determining position of space in formula, and l represents the number of some error microphone, thenHead is worked as in expression
Signal of the portion in α at l-th of error microphone,Represent original at l-th of error microphone when head is in α
Noise source signal, the reversed acoustic signals sent out for system behind plus sige, then the formula represents, misses for l-th when head is in α
The noise signal after noise reduction is carried out at poor microphone, Ben Shi represents the noise signal near a microphone, if practical
There are several microphones in system, need to calculate the signal progress phase adduction of several microphones into the whole energy of each position accordingly
Amount, the integral energy of alpha position are represented with the following formula：
Each position that S4, the probability that each position head part appearance counted is measured with reference to step S2 and step S3 are acquired
Integral energy space gross energy J in noise reduction region is obtained：
S5, the space gross energy that step S4 is acquired is minimized, derives the optimal parameter w of controller.
Preferably, step S2 is by recording the characteristics of motion of each tested person personnel in noise reduction regional space so as to remember
Record the probability residing for spatial points head part.
Preferably, obtained during gross energy minimum at human ear in the parameter of controller, that is, noise reduction experiment early period measurement process
Parameter value, noise reduction is best during ensureing noise reduction system practical application.
The probability that the present invention goes out each position energy according to the probability calculation of head part's movement each position is averaged, whole so as to obtain
The energy in body space in the case of this integral energy minimum, calculates controller parameter for controller in practical application
Design, so as to make the noise reduction of system optimal, it can be understood as in the position noise reduction for needing most noise reduction, to realize overall space
Optimal noise reduction is not being needed to add additional noise, not needed in human ear road in the case of setting sensor, reduction well
The noise level in space near human ear.
Description of the drawings
Fig. 1 is the flow chart of the present invention.
Specific embodiment
The embodiment of the present invention is described below in detail, the example of the embodiment is shown in the drawings, wherein from beginning to end
Same or similar label represents same or similar element or the element with same or like function.Below with reference to attached
The embodiment of figure description is exemplary, and is only used for explaining the present invention, and is not considered as limiting the invention.
In the description of the present invention, unless otherwise prescribed and limit, it should be noted that term " installation ", " connected ",
" connection " should be interpreted broadly, for example, it may be the connection inside mechanical connection or electrical connection or two elements, it can
To be to be connected directly, can also be indirectly connected by intermediary.It for the ordinary skill in the art, can basis
Concrete condition understands the concrete meaning of above-mentioned term.
Referring to a kind of design side of Fig. 1 descriptions active noise reduction system controller parameter according to embodiments of the present invention
Method includes the following steps：
S1, it determines noise reduction region, multiple incoming devices and output device is installed in the space for needing noise reduction；
S2, the probability that statistics head part appears in each random site is measured；
The probability of each position energy that S3, calculating head part occur is averaged, due to having measured statistics in previous experiments
Next probability existing for the head part of spatial points position needs the energy of measurement space each point, within this space error wheat
The quantity of gram wind may be 2,4 or n, depending on the feature with specific space：
α represents some determining position of space in formula, and l represents the number of some error microphone, thenHead is worked as in expression
Signal of the portion in α at l-th of error microphone,Represent original at l-th of error microphone when head is in α
Noise source signal, the reversed acoustic signals sent out for system behind plus sige, then the formula represents, misses for l-th when head is in α
The noise signal after noise reduction is carried out at poor microphone, Ben Shi represents the noise signal near a microphone, if practical
There are several microphones in system, need to calculate the signal progress phase adduction of several microphones into the whole energy of each position accordingly
Amount, the integral energy of alpha position are represented with the following formula：
Each position that S4, the probability that each position head part appearance counted is measured with reference to step S2 and step S3 are acquired
Integral energy space gross energy J in noise reduction region is obtained：
S5, the space gross energy that step S4 is acquired is minimized, derives the optimal parameter w of controller.
Preferably, step S2 is by recording the characteristics of motion of each tested person personnel in noise reduction regional space so as to remember
Record the probability residing for spatial points head part.
Preferably, obtained during gross energy minimum at human ear in the parameter of controller, that is, noise reduction experiment early period measurement process
Parameter value, noise reduction is best during ensureing noise reduction system practical application.
The probability that the present invention goes out each position energy according to the probability calculation of head part's movement each position is averaged, whole so as to obtain
The energy in body space in the case of this integral energy minimum, calculates controller parameter for controller in practical application
Design, so as to make the noise reduction of system optimal, it can be understood as in the position noise reduction for needing most noise reduction, to realize overall space
Optimal noise reduction is not being needed to add additional noise, not needed in human ear road in the case of setting sensor, reduction well
The noise level in space near human ear.
In the description of this specification, the description meaning of reference term " one embodiment ", " example " or " some examples " etc.
Refer at least one reality that the present invention is contained in reference to the embodiment or example particular features, structures, materials, or characteristics described
It applies in example or example.In the present specification, schematic expression of the above terms are not necessarily referring to identical embodiment or show
Example.Moreover, particular features, structures, materials, or characteristics described can be in any one or more embodiments or example
It combines in an appropriate manner.
Although an embodiment of the present invention has been shown and described, it will be understood by those skilled in the art that：Not
In the case of being detached from the principle of the present invention and objective a variety of change, modification, replacement and modification can be carried out to these embodiments, this
The range of invention is limited by claim and its equivalent.
Claims (3)
1. a kind of design method of active noise reduction system controller parameter, which is characterized in that include the following steps：
S1, it determines noise reduction region, multiple incoming devices and output device is installed in the space for needing noise reduction；
S2, the probability that statistics head part appears in each random site is measured；
The probability of each position energy that S3, calculating head part occur is averaged, and sky has been counted due to having been measured in previous experiments
Between probability existing for various point locations head part, next need the energy of measurement space each point, within this space error microphone
Quantity may be 2,4 or n, depending on the feature with specific space：
α represents some determining position of space in formula, and l represents the number of some error microphone, and m is the number of noise reduction horn, n
It it is the time, j is last moment, thenRepresent the signal at l-th of error microphone when head is in α,It represents to work as
Raw noise source signal of the head in α at l-th of error microphone, the reversed acoustic signals sent out for system behind plus sige,
Wherein g (α) represents the coefficient for the transmission function that loud speaker is responded in alpha position, which is measured by measurement early period, and u, which is represented, to be driven
Dynamic loud speaker sends out the signal of noise reduction sound wave, and u (α) represents the drive signal in alpha position, then the formula represents, works as head
The noise signal after noise reduction is carried out at l-th of error microphone in α, Ben Shi represents the noise letter near a microphone
Number, if there are several microphones in real system, need the signal progress phase adduction of several microphones calculating everybody accordingly
The integral energy put, the integral energy of alpha position are represented with the following formula：
T in formula refers to the transposition of matrix, and what E was represented is expectation of a random variable, and d (α) represents noise signal, and R passes through formula
Conversion is got；
S4, the whole of probability that each position head part of statistics occurs and each positions for acquiring of step S3 is measured with reference to step S2
The space gross energy J in noise reduction region is obtained in physical efficiency amount：
F (α) represents the probability density function at α in formula, is got by the measurement statistics random to each location probability early period,
A, b, c are the formula transfer of gross energy to be asked to get the parameter come in return with reference to step S2 and S3；
S5, the space gross energy that step S4 is acquired is minimized, derives the optimal parameter w of controller.
2. design method according to claim 1, which is characterized in that step S2 is being dropped by recording each tested person personnel
The characteristics of motion made an uproar in regional space is so as to record the probability residing for spatial points head part.
3. design method according to claim 1, which is characterized in that the parameter of controller, that is, noise reduction experiment early period measured
The parameter value obtained during gross energy minimum at human ear in journey, noise reduction is best during ensureing noise reduction system practical application.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810010686.0A CN108229030B (en) | 2018-01-05 | 2018-01-05 | Design method of controller parameters of active noise reduction system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810010686.0A CN108229030B (en) | 2018-01-05 | 2018-01-05 | Design method of controller parameters of active noise reduction system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108229030A true CN108229030A (en) | 2018-06-29 |
| CN108229030B CN108229030B (en) | 2021-07-02 |
Family
ID=62643065
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810010686.0A Active CN108229030B (en) | 2018-01-05 | 2018-01-05 | Design method of controller parameters of active noise reduction system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108229030B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108877758A (en) * | 2018-07-04 | 2018-11-23 | 北京安声科技有限公司 | A kind of spatial field active denoising method |
| CN109714668A (en) * | 2019-01-07 | 2019-05-03 | 哈尔滨工业大学（深圳） | Indoor active noise reduction device, noise-reduction method and storage medium |
| CN111083591A (en) * | 2019-12-13 | 2020-04-28 | 恒玄科技(北京)有限公司 | Configuration method and device of noise reduction earphone and noise reduction earphone |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103244465A (en) * | 2013-05-03 | 2013-08-14 | 湖北文理学院 | Active noise control system for fan of portable computer |
| CN103901401A (en) * | 2014-04-10 | 2014-07-02 | 北京大学深圳研究生院 | Binaural sound source positioning method based on binaural matching filter |
| WO2015100430A1 (en) * | 2013-12-24 | 2015-07-02 | Digimarc Corporation | Methods and system for cue detection from audio input, low-power data processing and related arrangements |
| US9619755B2 (en) * | 2013-10-23 | 2017-04-11 | Mitsubishi Electric Research Laboratories, Inc. | Method for estimating parameters of a graph spectral filter using training data |
| CN107464552A (en) * | 2017-08-24 | 2017-12-12 | 徐银海 | A kind of distributed locomotive active noise reduction system and method |
-
2018
- 2018-01-05 CN CN201810010686.0A patent/CN108229030B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103244465A (en) * | 2013-05-03 | 2013-08-14 | 湖北文理学院 | Active noise control system for fan of portable computer |
| US9619755B2 (en) * | 2013-10-23 | 2017-04-11 | Mitsubishi Electric Research Laboratories, Inc. | Method for estimating parameters of a graph spectral filter using training data |
| WO2015100430A1 (en) * | 2013-12-24 | 2015-07-02 | Digimarc Corporation | Methods and system for cue detection from audio input, low-power data processing and related arrangements |
| CN103901401A (en) * | 2014-04-10 | 2014-07-02 | 北京大学深圳研究生院 | Binaural sound source positioning method based on binaural matching filter |
| CN107464552A (en) * | 2017-08-24 | 2017-12-12 | 徐银海 | A kind of distributed locomotive active noise reduction system and method |
Non-Patent Citations (3)
| Title |
|---|
| ANDRZEJ STANKIEWICZ: "Magnetic Reader Instability as a Source of Low-Frequency Noise", 《IEEE TRANSACTIONS ON MAGNETICS》 * |
| LIU, LIANG 等: "Optimized Uplink Transmission in Multi-Antenna C-RAN With Spatial Compression and Forward", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
| 魏庆朝 等: "公路隧道内主动降噪声场研究", 《中国公路学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108877758A (en) * | 2018-07-04 | 2018-11-23 | 北京安声科技有限公司 | A kind of spatial field active denoising method |
| CN109714668A (en) * | 2019-01-07 | 2019-05-03 | 哈尔滨工业大学（深圳） | Indoor active noise reduction device, noise-reduction method and storage medium |
| CN111083591A (en) * | 2019-12-13 | 2020-04-28 | 恒玄科技(北京)有限公司 | Configuration method and device of noise reduction earphone and noise reduction earphone |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108229030B (en) | 2021-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108877758A (en) | A kind of spatial field active denoising method | |
| CN108229030A (en) | A kind of design method of active noise reduction system controller parameter | |
| Ishizaka et al. | Direct determination of vocal tract wall impedance | |
| CN101378607B (en) | Sound processing apparatus and method for correcting phase difference | |
| CN104239736A (en) | Structure damage diagnosis method based on power spectrum and intelligent algorithms | |
| EP3163914A3 (en) | Sound level estimation | |
| WO2009070476A1 (en) | Active noise cancellation using a predictive model approach | |
| Pfiffner et al. | A MEMS condenser microphone-based intracochlear acoustic receiver | |
| WO2009144721A3 (en) | Method and apparatus for examining subjects for particular physiological conditions utilizing acoustic information | |
| MY147063A (en) | Ultrasonic surface monitoring | |
| Hudde et al. | Measuring and modeling basic properties of the human middle ear and ear canal. Part III: Eardrum impedances, transfer functions and model calculations | |
| CN114580071B (en) | Adjacent floor shock insulation damping comfort evaluation method based on computer aided design | |
| KR20160148057A (en) | Ear-shaped part, artificial head part, and measurement system and measurement method using these | |
| CN108644965A (en) | A kind of air-conditioning and the distributed active denoising method of central air-conditioning | |
| WO2004092700A3 (en) | A method and device for determining acoustical transfer impedance | |
| JP5442215B2 (en) | Ultrasonic distance measurement system | |
| WO2005043949A3 (en) | Method for obtaining diagnostic information relating to a patient having an implanted transducer | |
| CN109871824A (en) | The multi-modal separation method of supersonic guide-wave and its system based on management loading | |
| CN109900451A (en) | The method for correcting the distortion of wind tunnel experiment pressure measuring model wind pressure signal | |
| WO2016192277A1 (en) | Bone conduction sound transmission device and method | |
| CN103813258B (en) | Method and system for acquiring diaphragm compliance | |
| CN102736064A (en) | Compression sensor-based positioning method of sound source of hearing aid | |
| Rosowski et al. | The radiation impedance of the external ear of cat: Measurements and applications | |
| EP2133665A3 (en) | System and method for objective self-diagnosis of measurement device calibration condition | |
| CN102370485B (en) | Digital automatic acoustic-control otoacoustic immittance detection system and digital signal compensation method thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |