CN108229030B - Design method of controller parameters of active noise reduction system - Google Patents

Design method of controller parameters of active noise reduction system Download PDF

Info

Publication number
CN108229030B
CN108229030B CN201810010686.0A CN201810010686A CN108229030B CN 108229030 B CN108229030 B CN 108229030B CN 201810010686 A CN201810010686 A CN 201810010686A CN 108229030 B CN108229030 B CN 108229030B
Authority
CN
China
Prior art keywords
noise reduction
alpha
head
energy
probability
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.)
Active
Application number
CN201810010686.0A
Other languages
Chinese (zh)
Other versions
CN108229030A (en
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.)
Beijing Ancsonic Technology Co ltd
Original Assignee
Beijing Ancsonic Technology Co 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 Beijing Ancsonic Technology Co ltd filed Critical Beijing Ancsonic Technology Co ltd
Priority to CN201810010686.0A priority Critical patent/CN108229030B/en
Publication of CN108229030A publication Critical patent/CN108229030A/en
Application granted granted Critical
Publication of CN108229030B publication Critical patent/CN108229030B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10Noise analysis or noise optimisation

Abstract

The invention discloses a design method of controller parameters of an active noise reduction system, which is characterized in that the probability average of energy of each position is calculated according to the probability of the head moving each position, the total spatial energy in a noise reduction area is calculated by combining the probability of the head appearing at each position counted by measurement and the calculated total energy of each position, and the calculated controller parameters are used for the design of a controller in practical application under the condition of the minimum total energy, so that the noise reduction effect of the system is optimal.

Description

Design method of controller parameters of active noise reduction system
Technical Field
The invention relates to the field of noise reduction, in particular to a design method of parameters of a controller of an active noise reduction system.
Background
When designing an active noise reduction system aiming at a spatial field near a human ear, a mold human head or a human body is generally required to be used for preliminary testing, specifically, a sensor is placed at an ear canal of the mold, characteristics of a noise field at the position of the human ear and a transfer function of a secondary sound source at the position of the human ear are measured, after the preliminary testing is finished, measured data are stored in a signal processor, and system parameters of an active noise reduction controller are designed according to the measured results, but when the system is applied to a real environment, noise reduction is aimed at the position of the human ear when a user is actually used, however, the position of the head of the user is uncertain in the using process, random movement in a certain range is generated near the position of the mold human head of the previous testing, which causes that a sound signal generated by an excitation source of a secondary speaker at the position of the actual human ear is inconsistent with the data obtained in the previous testing, therefore, the noise reduction effect of the controller designed according to the previous test data can be influenced during actual use, and then in actual application, it is impossible to place a sensor at the position of the human ear of a user to measure the transfer function of the secondary sound source at the position of the actual human ear in real time. Even if the user wears the sensor at the ear, the real-time measurement method cannot achieve the ideal effect because the transfer function of the measuring loudspeaker needs to use the loudspeaker to send an additional noise excitation signal while actively reducing noise, and then the real-time test can be completed by extracting a response signal corresponding to the additional excitation from the ear sensor, but the additional excitation sound source can deteriorate the active noise reduction effect.
Disclosure of Invention
The invention aims to provide a method for designing parameters of a controller of an active noise reduction system, which can well reduce the noise level of a space near the human ear aiming at the random variation of the head position of a user under the conditions of not adding additional noise and not arranging a sensor at the human ear.
In order to achieve the technical purpose, the invention adopts the following technical scheme that the method for designing the parameters of the controller of the active noise reduction system comprises the following steps:
s1, determining a noise reduction area, and installing a plurality of input devices and output devices in a space needing noise reduction;
s2, measuring and counting the probability of the head of the person appearing at each random position;
s3, calculating the probability average of the energy of each position where the head appears, because the probability of the head existing at each point of the space is measured and counted in the previous period of experiment, the energy of each point of the space needs to be measured next, the number of error microphones in the space may be 2, 4 or n, and depends on the characteristics of the specific space:
Figure BDA0001540192320000021
where alpha represents a certain determined position in space and l represents the number of a certain error microphone, then
Figure BDA0001540192320000022
Representing the signal at the ith error microphone when the head is at alpha,
Figure BDA0001540192320000023
the original noise source signal at the ith error microphone when the head is at alpha is shown, and the plus sign is followed by the backward sound wave signal sent by the system, then the formula shows that the noise signal after noise reduction has been carried out at the ith error microphone when the head is at alpha, the formula represents the noise signal near one microphone, if there are several microphones in the actual system, the signals of the several microphones need to be added, and the overall energy of each position is calculated according to the formula, and the overall energy of the alpha position is shown by the following formula:
Figure BDA0001540192320000024
s4, combining the probability of the human head appearing at each position obtained by the measurement statistics of the step S2 and the overall energy of each position obtained by the step S3 to obtain the total spatial energy J in the noise reduction region:
Figure BDA0001540192320000025
and S5, obtaining the minimum value of the total spatial energy obtained in the step S4, and deriving the optimal parameter w of the controller.
Preferably, step S2 records the probability of each human head in the space by recording the motion law of each human subject in the noise reduction region space.
Preferably, the parameters of the controller are parameter values obtained when the total energy at the human ear is minimum in the early noise reduction experiment measurement process, so that the best noise reduction effect is ensured when the noise reduction system is actually applied.
The invention calculates the probability average of the energy of each position according to the probability of the head moving each position, thereby obtaining the energy of the whole space, under the condition that the whole energy is minimum, the parameters of the controller are calculated and used for the design of the controller in practical application, thereby leading the noise reduction effect of the system to be optimal, being understood as noise reduction at the position which needs noise reduction most, realizing the optimal noise reduction effect of the whole space, and well reducing the noise level of the space near the human ear under the condition that additional noise is not required and a sensor is not required to be arranged in the human ear canal.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
In the description of the present invention, unless otherwise specified and limited, it is to be understood that the terms "mounted," "connected," and "connected" are used broadly and can be, for example, mechanically or electrically connected, or can be internal to two elements, directly connected, or indirectly connected through an intermediate medium. The specific meaning of the above terms can be understood by those of ordinary skill in the art as appropriate.
The following describes a design method of parameters of an active noise reduction system controller according to an embodiment of the present invention with reference to fig. 1, including the following steps:
s1, determining a noise reduction area, and installing a plurality of input devices and output devices in a space needing noise reduction;
s2, measuring and counting the probability of the head of the person appearing at each random position;
s3, calculating the probability average of the energy of each position where the head appears, because the probability of the head existing at each point of the space is measured and counted in the previous period of experiment, the energy of each point of the space needs to be measured next, the number of error microphones in the space may be 2, 4 or n, and depends on the characteristics of the specific space:
Figure BDA0001540192320000031
where alpha represents a certain determined position in space and l represents the number of a certain error microphone, then
Figure BDA0001540192320000032
Representing the signal at the ith error microphone when the head is at alpha,
Figure BDA0001540192320000033
the original noise source signal at the ith error microphone when the head is at alpha is shown, and the plus sign is followed by the backward sound wave signal sent by the system, then the formula shows that the noise signal after noise reduction has been carried out at the ith error microphone when the head is at alpha, the formula represents the noise signal near one microphone, if there are several microphones in the actual system, the signals of the several microphones need to be added, and the overall energy of each position is calculated according to the formula, and the overall energy of the alpha position is shown by the following formula:
Figure BDA0001540192320000041
s4, combining the probability of the human head appearing at each position obtained by the measurement statistics of the step S2 and the overall energy of each position obtained by the step S3 to obtain the total spatial energy J in the noise reduction region:
Figure BDA0001540192320000042
and S5, obtaining the minimum value of the total spatial energy obtained in the step S4, and deriving the optimal parameter w of the controller.
Preferably, step S2 records the probability of each human head in the space by recording the motion law of each human subject in the noise reduction region space.
Preferably, the parameters of the controller are parameter values obtained when the total energy at the human ear is minimum in the early noise reduction experiment measurement process, so that the best noise reduction effect is ensured when the noise reduction system is actually applied.
The invention calculates the probability average of the energy of each position according to the probability of the head moving each position, thereby obtaining the energy of the whole space, under the condition that the whole energy is minimum, the parameters of the controller are calculated and used for the design of the controller in practical application, thereby leading the noise reduction effect of the system to be optimal, being understood as noise reduction at the position which needs noise reduction most, realizing the optimal noise reduction effect of the whole space, and well reducing the noise level of the space near the human ear under the condition that additional noise is not required and a sensor is not required to be arranged in the human ear canal.
In the description herein, references to the description of "one embodiment," "an example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the invention have been shown and described, it will be understood by those of ordinary skill in the art that: various changes, modifications, substitutions and alterations can be made to the embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents.

Claims (3)

1. A design method for parameters of a controller of an active noise reduction system is characterized by comprising the following steps:
s1, determining a noise reduction area, and installing a plurality of microphones and loudspeakers in a space needing noise reduction;
s2, measuring and counting the probability of the head of the person appearing at each random position;
s3, calculating the probability average of the energy of each position where the head appears, because the probability of the head existing at each point position in the space is measured and counted in the previous period of experiment, the energy of each point in the space needs to be measured, the number of the error microphones in the space is any number, and the error microphones depend on the characteristics of the specific space:
Figure FDA0002939408200000011
where α represents a certain determined position in space, l represents the number of a certain error microphone, m is the number of a noise reduction horn, n is time, and j is the last time, then
Figure FDA0002939408200000012
Representing the signal at the ith error microphone when the head is at alpha,
Figure FDA0002939408200000013
representing the original noise source signal at the ith error microphone when the head is at alpha, followed by the system-emitted inverse sound wave signal, where g (alpha) represents the coefficient of the transfer function of the loudspeaker response at alpha position, measured by prior measurements, u represents the signal driving the loudspeaker to emit noise-reduced sound waves, and u (alpha) represents the drive signal at alpha position, then
Figure FDA0002939408200000014
It is shown that the noise signal after noise reduction has been performed at the ith error microphone when the head is at alpha,
Figure FDA0002939408200000015
representing noise signals near one microphone, if there are several microphones in the actual system, it is necessary to add the signal energies of several microphones and calculate the overall energy of each position accordingly, and the overall energy of the α position is expressed by the following formula:
Figure FDA0002939408200000016
in the formula
Figure FDA0002939408200000017
For by a plurality of error microphones respectively
Figure FDA0002939408200000018
The column vector is composed, T is the transpose, E represents the expectation of the random variable,
Figure FDA0002939408200000019
for by a plurality of error microphones respectively
Figure FDA00029394082000000110
Constituent column vectors, R(α)(n) is obtained by conversion of a formula;
s4, combining the probability of the human head appearing at each position obtained by the measurement statistics of the step S2 and the overall energy of each position obtained by the step S3 to obtain the total spatial energy J in the noise reduction region:
Figure FDA00029394082000000111
f (alpha) represents a probability density function at alpha, and is obtained through the measurement statistics of random probability of each position in the early stage, and A, b and c are parameters converted from the formula of combining the steps S2 and S3 to obtain total energy;
s5, obtaining the minimum value of the total spatial energy obtained in the step S4, and deducing the optimal parameters of the controller
Figure FDA0002939408200000021
2. The design method of claim 1, wherein step S2 records the probability of each human head in the noise reduction region space by recording the motion law of each human subject in the space.
3. The design method of claim 1, wherein the optimization of the controllerParameter(s)
Figure FDA0002939408200000022
Namely, the parameter value obtained when the total energy at the human ear is minimum in the early noise reduction experiment measurement process S1-S5, so as to ensure that the noise reduction effect is optimal when the noise reduction system is actually applied.
CN201810010686.0A 2018-01-05 2018-01-05 Design method of controller parameters of active noise reduction system Active CN108229030B (en)

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 CN108229030A (en) 2018-06-29
CN108229030B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
Magnetic Reader Instability as a Source of Low-Frequency Noise;Andrzej Stankiewicz;《IEEE Transactions on Magnetics》;20171130;第53卷(第11期);第3302204页 *
Optimized Uplink Transmission in Multi-Antenna C-RAN With Spatial Compression and Forward;Liu, Liang 等;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20151001;第63卷(第19期);第5083-5095页 *
公路隧道内主动降噪声场研究;魏庆朝 等;《中国公路学报》;20170131;第30卷(第1期);第77-82页 *

Also Published As

Publication number Publication date
CN108229030A (en) 2018-06-29

Similar Documents

Publication Publication Date Title
CN108229030B (en) Design method of controller parameters of active noise reduction system
US11638085B2 (en) System, device and method for assessing a fit quality of an earpiece
RU2595636C2 (en) System and method for audio signal generation
US11665488B2 (en) Auditory device assembly
CN111833896B (en) Voice enhancement method, system, device and storage medium for fusing feedback signals
US8045737B2 (en) Method of obtaining settings of a hearing instrument, and a hearing instrument
CN105556260B (en) broadband acoustical holography
US20080064956A1 (en) Method of displaying an elastic image
JP2017503388A5 (en)
KR20080092404A (en) System and method for utilizing inter-microphone level differences for speech enhancement
US11024421B2 (en) Method for automatic determination of an individual function of a DPOAE level
CN108877758A (en) A kind of spatial field active denoising method
CN110715724B (en) Method for measuring and calculating noise in vehicle based on noise of rack and acoustic transmission function
US9414157B2 (en) Method and device for reducing voice reverberation based on double microphones
CN110960224B (en) Hearing threshold and/or hearing status detection systems and methods
JP5737342B2 (en) Sound collecting device and program
US20180041846A1 (en) Method for feedback suppression
JP2000194400A (en) Method and device for processing noisy acoustic signal
EP3644307A1 (en) Tuning method, manufacturing method, computer-readable storage medium and tuning system
TWI429885B (en) Method for visualizing sound source energy distribution in reverberant environment
CN104248439A (en) Eardrum acoustic pressure detection method and eardrum acoustic pressure detection system
CN111968615A (en) Noise reduction processing method and device, terminal equipment and readable storage medium
CN111998934B (en) Sound source sound power testing method
JP4746498B2 (en) Unidirectional microphone
US20140296739A1 (en) Method to measure real-ear-to-coupler difference

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