CN105338462B - A kind of implementation method for reappearing hearing aid insertion gain - Google Patents
A kind of implementation method for reappearing hearing aid insertion gain Download PDFInfo
- Publication number
- CN105338462B CN105338462B CN201510921262.6A CN201510921262A CN105338462B CN 105338462 B CN105338462 B CN 105338462B CN 201510921262 A CN201510921262 A CN 201510921262A CN 105338462 B CN105338462 B CN 105338462B
- Authority
- CN
- China
- Prior art keywords
- hearing aid
- gain
- frequency
- reappearing
- insertion gain
- 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.)
- Expired - Fee Related
Links
Abstract
The invention discloses a kind of implementation methods of reproduction hearing aid insertion gain, and first in specified frequency range, measurement obtains N number of frequency-gain pair of hearing aid to be measured；According to obtained N number of frequency-gain pair, discrete Fu Ye inverse transformation is carried out to it, obtains the unit impulse response h (n) for reappearing hearing aid insertion gain curve；Fourier transformation is carried out to the h (n) again, obtains the frequency response H (e for reappearing hearing aid insertion gain characteristic curvejw)=FFT [h (n)], and thereby determine that coefficient vectorIt willIt is expressed as reappearing the system function of hearing aid insertion gain characteristic curve；According to the system functionObtain the characteristic curve for reappearing hearing aid insertion gain.What this method can objectively reappear hearing aid is actually inserted into gain, to provide advantageous support to objectively evaluating for hearing aid.
Description
Technical field
The present invention relates to hearing aid apparatus technical field more particularly to a kind of realization sides for reappearing hearing aid insertion gain
Method.
Background technique
With being constantly progressive for science and technology, the trend such as intelligent, digitlization, the height that hearing aid has is presented in hearing aid
Hold performance more and more abundant, difference in selling prices is also increasing.And the quality of hearing aid how is evaluated, and it is hearing aid producer, consumption
Person, and test the problem of all paying close attention to very much with doctor.Objective evaluation is carried out to hearing aid at present to be also confined in physical acoustics parameter,
It is international and domestic to have relevant standard criterion.Important performance of the voice quality as hearing aid, assessment mainly are also confined to lead
It sees in evaluation, and subjective evaluating method is not suitable as gauge check means because it expends a large amount of human and material resources and time, it must
Hearing aid voice quality must be evaluated using science objective method, just be avoided that market price confusion, enterprise is unhealthy
Competition, tests and looks at flowers in a fog with doctor, to finally damage consumer's interests.
And the voice quality of hearing aid is objectively evaluated, how hearing aid reappears hearing aid and exists as flight data recorder
Actual measurement insertion gain under various setting states just becomes a major issue.If the insertion that can reappear hearing aid increases
Benefit will provide advantageous support to objectively evaluating for hearing aid, since hearing aid brand is numerous, and every kind of brand hearing aid
Barrier concrete condition is listened test with adjustment according to different patients again during wearing, thus its actual measurement insertion increases
Benefit be it is multifarious, it reappeared, be very difficult thing, in the prior art there is no such solution party
Case.
Summary of the invention
The object of the present invention is to provide a kind of implementation method of reproduction hearing aid insertion gain, this method can be weighed objectively
Existing hearing aid is actually inserted into gain, to provide advantageous support to objectively evaluating for hearing aid.
A kind of implementation method for reappearing hearing aid insertion gain, the method includes：
In specified frequency range, measurement obtains N number of frequency-gain pair of hearing aid to be measured；
According to obtained N number of frequency-gain pair, discrete Fu Ye inverse transformation is carried out to it, is obtained for reappearing hearing aid
It is inserted into the unit impulse response h (n) of gain curve, which is expressed as：
Fourier transformation is carried out to the h (n) again, obtains the frequency response H for reappearing hearing aid insertion gain characteristic curve
(ejw)=FFT [h (n)], and thereby determine that coefficient vectorIt willIt is expressed as reappearing hearing aid insertion gain spy
The system function of linearity curve；
According to the system functionObtain the characteristic curve for reappearing hearing aid insertion gain.
As seen from the above technical solution provided by the invention, the reality that this method can objectively reappear hearing aid is inserted
Enter gain, to provide advantageous support to objectively evaluating for hearing aid.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill in field, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is provided the implementation method flow diagram for reappearing hearing aid insertion gain by the embodiment of the present invention；
Fig. 2 is in example of the present invention to characteristic verification result analogous diagram.
Specific embodiment
With reference to the attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based on this
The embodiment of invention, every other implementation obtained by those of ordinary skill in the art without making creative efforts
Example, belongs to protection scope of the present invention.
The method of the embodiment of the present invention is the insertion gain data using hearing aid actual measurement, designs filter guest
It sees ground and reappears hearing aid insertion gain.The embodiment of the present invention is described in further detail below in conjunction with attached drawing, such as Fig. 1 institute
It is shown as the embodiment of the present invention and the implementation method flow diagram for reappearing hearing aid insertion gain is provided, the method includes：
Step 11：In specified frequency range, measurement obtains N number of frequency-gain pair of hearing aid to be measured；
In this step, frequency-gain pair of hearing aid can be measured at third-octave in [0-8khz] range, but
It is to be not limited at third-octave, in order to obtain more accurate reproduction, the measurement at arbitrary frequency point can be increased, and
And frequency range can any reducing and expansion according to the actual situation.
In the specific implementation, downloading audiogram first against different patients, hearing aid test matching, or by hearing aid
Items setting adjustment is under the specific condition of needs；
Then according to iec 60118-8 canonical measure condition, the insertion gain of hearing aid is tested, in frequency model of concern
In enclosing, measurement obtains the yield value at corresponding one third octave, obtains frequency-gain pair of hearing aid, as wishing to force
The sampled value of the insertion gain actual frequency response of close hearing aid actual measurement, this is to reappear hearing aid insertion gain curve
Foundation.
Step 12：According to obtained N number of frequency-gain pair, discrete Fu Ye inverse transformation is carried out to it, is obtained for reappearing
The unit impulse response h (n) of hearing aid insertion gain curve；
Here, which is expressed as：
W in above formulakIndicate different measurement frequency points；
In the concrete realization, it is assumed that the actual frequency response for the insertion gain that hearing aid actual measurement obtains is Hd(ejw),
The hearing aid then measured in step 11 is considered as in frequency domain 0 in the yield value of each Frequency point and arrives nyquist frequency
(the i.e. half of sample frequency fs：Fs/2 to H between)d(ejw) carry out N point sampling value Hd(k)：
Then with the Hd(k) it as reproduction hearing aid insertion gain-frequency characterisitic sampled value H (k), and enables in above formula：
By discrete Fu Ye inverse transformation, the unit pulse for obtaining reappearing hearing aid insertion gain curve with above-mentioned H (k) is rung
H (n) is answered, i.e.,：
In addition, the reproduction curve of hearing aid insertion gain can be indicated with following formula：
I.e. hearing aid insertion gain can be simulated with a FIR filter, and wherein M is the filter order of reproduction curve,
The selection of order M can preferably be fitted above-mentioned N number of sampled value as standard using the designed hearing aid reproduction curve come out.
Step 13：Fourier transformation is carried out to the h (n) again, obtains the frequency for reappearing hearing aid insertion gain characteristic curve
Respond H (ejw)=FFT [h (n)], and thereby determine that coefficient vectorIt can then incite somebody to actionIt is expressed as reappearing hearing aid
It is inserted into the system function of gain characteristic curve；
Wherein, M is the filter order of reproduction curve, brFor as unit impulse response through determined by Fourier transform filter
Coefficient vectorIn digits.
Specifically, by the interpolation formula in frequency domain sample theorem it is found that can be in the hope of reproduction hearing aid using above-mentioned H (k)
Frequency response H (the e of device insertion gain characteristic curvejw), this H (ejw) genuine property of hearing aid insertion gain will be approached
Frequency response Hd(ejw), H (ejw) interpolation formula be：
In formula：It is interpolating function.
It brings interpolating function into interpolation formula, can be obtained after abbreviation：
Wherein, w is angular frequency, w in above-mentioned formulakIndicate different measurement frequency points.
As can be seen from the above equation, on each stepped-frequency signal, the hearing aid of reproduction insertion gain strictly with hearing aid
The frequency response numerical value of true insertion gain is equal, i.e.,
Step 14：According to the system functionObtain the characteristic curve for reappearing hearing aid insertion gain.
In addition, in the specific implementation, can also be by measuring two side frequency point [wk,wk+1] between transition period sampling
Frequency-gain pair, verify whether acquired hearing aid insertion gain characteristic curve meets the requirements；
If the requirements are not met, then frequency-gain that transition period samples is added to original sampled data to data
In, it reacquires hearing aid and is inserted into gain characteristic curve, until meeting the requirements.
Specifically, by the above-mentioned frequency response H (e for acquiring practical filterjw) process it is found that between sampled point
Frequency response is the Weighted Interpolation function superposition by each sampled point and is formed, thus has certain approximate error, and the error is big
The small shape depending on ideal frequency response, ideal Frequency Response variation is gentler, and interpolated value is closer to ideal value, approximate error
It is smaller；, whereas if the ideal Frequency Response variation between sampled point is steeper, then the error between interpolated value and ideal value is got over
Greatly, thus on the both sides of ideal filter discontinuity point, spike will be generated, and ripple will be generated in passband and stopband, this
When can by increase sampled point, to be improved and be revised.
Above by frequency-gain pair of measurement transition period sampling, acquired hearing aid insertion gain characteristic is verified
The detailed process whether curve meets the requirements is：
An error rate is set first, and measures frequency-gain pair of transition period sampling；
If frequency-gain that measured transition period samples is inserted into gain characteristic with by the hearing aid to data
The error between frequency-insertion gain under the identical frequency that curve is calculated is less than set error rate, then acquired
Hearing aid insertion gain characteristic curve meet the requirements；Otherwise, it is unsatisfactory for requiring.
Above-mentioned implementation method is described in detail with specific example below：
(1) first against the audiogram of certain patient, hearing aid is adjusted to specific setting state.In this example, hearing aid
Volume is adjusted to maximum, closes noise reduction switch, selects omnidirectional's positioning, by hearing aid fits on emulation number of people KAMAR, passes through emulation
The number of people acquires signal to computer.
(2) under IEC 60118-8 standard, that tests hearing aid is actually inserted into gain, obtains corresponding one third frequency multiplication
Yield value at journey obtains N number of frequency-gain pair.Frequency point is freq=[0 0.2e3 0.25e3 0.315e3 0.4e3
0.5e3 0.63e3 0.8e3 1e3 1.25e3 1.6e3 2e3 2.5e3 3.15e3 4e3 5e3 6.3e3 8e3]；0 is
Initial point.The insertion gain of hearing aid at corresponding said frequencies point：[23.8,23.8,24.5,25,24.8,24.7,24.4,
26.2,18.5,13.4,15.9,14.9,11,16,13.3,20.2,17.6,20].
By it as the sampled value of the ideal situation frequency response for the hearing aid insertion gain characteristic for wishing to approach, this is to set
Count weight the foundation of existing hearing aid actual measurement insertion gain.
(3) characteristic curve for reappearing hearing aid insertion gain is obtained, is specifically included：
The first step：Assuming that the frequency response of the insertion gain of hearing aid actual measurement is Hd(ejw), then measurement obtains in (2)
Hearing aid be considered as in frequency domain 0 to 8khz (due to the frequency of hearing aid concern 0-8khz in the yield value of each Frequency point
Signal in range, thus sample frequency is set as 16khz, then nyquist frequency is between 8khz) to Hd(ejw) carry out 18
A point sampling value Hd(k), then with this Hd(k) the frequency characteristic sampled value H (k) as hearing aid insertion gain reproduction curve, i.e.,
Enable N=18：
wk=[0 0.2e3 0.25e3 0.315e3 0.4e3 0.5e3 0.63e3 0.8e3 1e3...1.25e3
1.6e3 2e3 2.5e3 3.15e3 4e3 5e3 6.3e3 8e3]
H (k)=Hd(k)=[23.8,23.8,24.5,25,24.8,24.7,24.4,26.2,18.5,13.4,15.9,
14.9,11,16,13.3,20.2,17.6,20].
Second step：Assuming that the reproduction curve order M=690 of hearing aid insertion gain, i.e.,
Third step：The unit impulse response h of the reproduction curve of hearing aid insertion gain is acquired with inverse discrete Fourier transform
(n)。
It is defined by DFT, hearing aid insertion gain can be uniquely determined with the sampled value H (k) of this 18 frequency domains and reappears song
The unit impulse response h (n) of line, i.e.,：
4th step：Further according to the definition of Fourier transformation, the frequency for acquiring hearing aid insertion gain reproduction characteristic profiles is rung
Answer H (ejw)。
And it thereby determines thatIn coefficient vectorShown in table 1 specific as follows, obtained system in table 1
NumberFor：The first row is followed successively by b from left to right0,b1,b2,b3, the second row is followed successively by b from left to right4,b5,b6,b7, and successively class
It pushes away：
Table 1
Further according to above system functionObtain the characteristic curve for reappearing hearing aid insertion gain.
Further, then H (e is verifiedjw) whether meet the requirement of wave filter technology index, allowable error ε=1% is set,
Increase new sample frequency point 7000hz between former sample frequency [6300hz 8000hz], reproduction curve is read by vernier
Value 18.67db at 7000hz, and the insertion gain at the Frequency point 7000hz that actual measurement obtains is 18.8db, it is seen that
The difference of the two is (18.8-18.67)=0.13db, accounts for the 0.69% of hearing aid insertion gain actual measured value, being less than allows
Error 1%, meets the requirements, and then verifies again and increases new Frequency point between other frequency bands, is also all satisfied requirement, then can recognize
The characteristic curve for being inserted into gain for designed reproduction hearing aid is met the requirements.
It is illustrated in figure 2 in example of the present invention to characteristic verification result analogous diagram, it can be with by the analogous diagram
Find out：Black circle is that the hearing aid that actual measurement obtains is inserted into gain, and solid line is then the spy for reappearing hearing aid insertion gain
Linearity curve, very well, if fitting is bad, needing to add sampled point or continuing to improve terminates N to solve for two parts fitting.
In conclusion being actually inserted into increasing using what method provided by the embodiment of the present invention can objectively reappear hearing aid
Benefit, and the yield value at sampled point is identical with the insertion gain that hearing aid actual measurement obtains for reproduction curve, thus to hearing aid
Objectively evaluating for device provides advantageous support.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Within the technical scope of the present disclosure, any changes or substitutions that can be easily thought of by anyone skilled in the art,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of claims
Subject to enclosing.
Claims (4)
1. a kind of implementation method for reappearing hearing aid insertion gain, which is characterized in that the method includes：
In specified frequency range, measurement obtains N number of frequency-gain pair of hearing aid to be measured；
According to obtained N number of frequency-gain pair, discrete Fu Ye inverse transformation is carried out to it, is obtained for reappearing hearing aid insertion
The unit impulse response h (n) of gain curve, the h (n) are expressed as：
Wherein, wkIndicate different measurement frequency points；
Fourier transformation is carried out to the h (n) again, obtains the frequency response H (e for reappearing hearing aid insertion gain characteristic curvejw)=
FFT [h (n)], and thereby determine that coefficient vectorIt willIt is expressed as reappearing hearing aid insertion gain characteristic song
The system function of line；Wherein, M is the filter order of reproduction curve, brFor identified coefficient vectorIn digits；Its
In,b0、b1、b2…bMAs coefficient vector
In digits；
According to the system functionObtain the characteristic curve for reappearing hearing aid insertion gain.
2. reappearing the implementation method of hearing aid insertion gain according to claim 1, which is characterized in that the method is also wrapped
It includes：
By measuring two side frequency point [wk,wk+1] between transition period sampling frequency-gain pair, verify acquired help
Listen whether device insertion gain characteristic curve meets the requirements；
If the requirements are not met, then is added to frequency-gain that transition period samples in original sampled data to data, weight
The new hearing aid that obtains is inserted into gain characteristic curve, until meeting the requirements.
3. reappearing the implementation method of hearing aid insertion gain according to claim 2, which is characterized in that described by measuring
Frequency-the gain pair sampled during crossing, verifies whether acquired hearing aid insertion gain characteristic curve meets the requirements, specific to wrap
It includes：
An error rate is set first, and measures frequency-gain pair of transition period sampling；
If frequency-gain that measured transition period samples is inserted into gain characteristic curve with by the hearing aid to data
The error between frequency-insertion gain under the identical frequency being calculated is less than set error rate, then acquired to help
Device insertion gain characteristic curve is listened to meet the requirements；Otherwise, it is unsatisfactory for requiring.
4. reappearing the implementation method of hearing aid insertion gain according to claim 1, which is characterized in that described in specified frequency
Within the scope of rate, specifically include：
To measure frequency-gain pair of hearing aid at third-octave in [0-8khz] range.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510921262.6A CN105338462B (en) | 2015-12-12 | 2015-12-12 | A kind of implementation method for reappearing hearing aid insertion gain |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510921262.6A CN105338462B (en) | 2015-12-12 | 2015-12-12 | A kind of implementation method for reappearing hearing aid insertion gain |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105338462A CN105338462A (en) | 2016-02-17 |
| CN105338462B true CN105338462B (en) | 2018-11-27 |
Family
ID=55288679
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510921262.6A Expired - Fee Related CN105338462B (en) | 2015-12-12 | 2015-12-12 | A kind of implementation method for reappearing hearing aid insertion gain |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105338462B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107221339B (en) * | 2017-05-22 | 2020-08-14 | 华北电力大学 | Gain compensation based hearing aid voice quality W-PESQ objective evaluation method |
| CN108040315A (en) * | 2017-10-23 | 2018-05-15 | 广东思派康电子科技有限公司 | A kind of test machine of computer-readable recording medium and the application medium |
| CN114786107B (en) * | 2022-05-10 | 2023-08-22 | 东南大学 | Hearing aid fitting method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4188667A (en) * | 1976-02-23 | 1980-02-12 | Beex Aloysius A | ARMA filter and method for designing the same |
| CN1598927A (en) * | 2004-08-31 | 2005-03-23 | 四川微迪数字技术有限公司 | Chinese voice signal process method for digital deaf-aid |
| CN102549659A (en) * | 2009-10-01 | 2012-07-04 | 高通股份有限公司 | Suppressing noise in an audio signal |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8611554B2 (en) * | 2008-04-22 | 2013-12-17 | Bose Corporation | Hearing assistance apparatus |
-
2015
- 2015-12-12 CN CN201510921262.6A patent/CN105338462B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4188667A (en) * | 1976-02-23 | 1980-02-12 | Beex Aloysius A | ARMA filter and method for designing the same |
| CN1598927A (en) * | 2004-08-31 | 2005-03-23 | 四川微迪数字技术有限公司 | Chinese voice signal process method for digital deaf-aid |
| CN102549659A (en) * | 2009-10-01 | 2012-07-04 | 高通股份有限公司 | Suppressing noise in an audio signal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105338462A (en) | 2016-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9363599B2 (en) | Systems and methods for protecting a speaker | |
| FI102700B (en) | Procedure for locating a single-phase grounding in an electricity distribution network | |
| US9357300B2 (en) | Systems and methods for protecting a speaker | |
| CN105338462B (en) | A kind of implementation method for reappearing hearing aid insertion gain | |
| Wen et al. | Triangular self-convolution window with desirable sidelobe behaviors for harmonic analysis of power system | |
| EP2355542A1 (en) | Control of a loudspeaker output | |
| US20170030958A1 (en) | Transformer parameter estimation using terminal measurements | |
| US11316501B1 (en) | Resampling technique for arbitrary sampling rate conversion | |
| CN105785124A (en) | Method for measuring harmonics and interharmonics of electric power system through spectrum estimation and cross correlation | |
| CN104363554B (en) | A kind of loudspeaker abnormal sound detection method | |
| CN104869519A (en) | Method and system for testing background noise of microphone | |
| CN109506762B (en) | Hydrophone received signal correction method based on filter | |
| CN104729677A (en) | Time domain digit weighting method for non-stable noise signals | |
| CN108736982A (en) | Acoustic communication processing method, device and electronic equipment | |
| US20170146587A1 (en) | Method and device for determining the fault location in the case of a fault on an electric line | |
| CN104840218B (en) | A kind of measurement method and electronic equipment of respiratory rate | |
| CN104568137B (en) | The decision method of specific noise source contribution rate in noise testing | |
| Zamora et al. | Smartphone tuning for accurate ambient noise assessment | |
| JP2016146576A (en) | Measuring method and measuring tool and correction method of reproduction characteristics of earphone and application program of measurement and application program of correction | |
| JP2014157110A (en) | Signal processing device, radar device and signal processing method | |
| Kowalczyk et al. | Virtual room acoustics using finite difference methods. How to model and analyse frequency-dependent boundaries? | |
| Bockstael et al. | Verifying the attenuation of earplugs in situ: Method validation using artificial head and numerical simulations | |
| Wijnings et al. | Characterization of mems microphone sensitivity and phase distributions with applications in array processing | |
| US8917818B2 (en) | Method and system for determining a number of load coils in a transmission line | |
| JP2013137251A (en) | Cable length measurement device, communication system and cable length measurement method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181127 Termination date: 20191212 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |