The invention relates to the field of fitting hearing aid devices. More particularly, it relates to a method for fitting a hearing aid device with active occlusion control to a user, said hearing aid device comprising:
Said means for active occlusion control comprising:
Said method comprising the steps of:
A hearing aid device is a device for aiding an individual in regard to its hearing. It may be a hearing aid or hearing prosthesis for compensating a hearing loss of its user. It may also be a hearing protection device which helps individuals to hear without damage in noisy environments. Such a device may transmit speech and attenuate noise by selective amplification. The occlusion effect is an effect experienced by individuals when an ear canal is fully or partially closed by an occluding object. In such a condition, the own voice of the individual and other body conducted sounds are perceived by him- or herself unnaturally loud. The earpiece of a hearing aid device can be such an occluding object. Active occlusion control is a method for reducing the occlusion effect actively. Actively means by destructive interference, i.e. emitting a kind of anti-sound. A passive occlusion control (or passive occlusion reduction) would be the provision of a large vent. However, hearing aids with a large vent are prone to feedback and cannot deliver loud low-frequency sound due to leakage from the canal to the outside and cannot provide good sound cleaning due to leakage from the outside into the canal. Providing hearing protective devices with a large vent renders them useless because low-frequency noise can pass without substantial attenuation through the vent. Occlusion is not to be confused with ampclusion. Users of hearing aid devices may perceive their own voice as being unnatural due to its amplification by the hearing aid device. Ampclusion can be counteracted by reducing the hearing aid device amplification in the frequency range of the users voice. Both occlusion control and ampclusion control aim for providing an own voice perceived as more natural.
U.S. Pat. No. 6,035,050 by Weinfurtner discloses a method for determining optimum parameter sets in a hearing aid. During an optimization phase an optimal user specific parameter set is allocated by selecting one of several trial parameter sets available.
WO 2004/021740 A1 by Rasmussen et al. discloses a method for counteracting the occlusion effect of an electronic device like a hearing aid. Sound conditions in the cavity between the ear piece and the tympanic membrane are determined. The transmission characteristics of the transmission path to the receiver counteracts the occlusion effect.
WO 2006/037156 A1 by Mejia et al. discloses an acoustically transparent occlusion reduction method. An electro-acoustic feedback network produces phase cancelling sounds in the ear. The integration with a hearing aid improves the user's perception of own voice.
WO 2008/017326 A1 by Nordahn discloses a method for in-situ occlusion effect measurement. A hearing aid comprises a microphone for external sounds and a microphone for sounds in the occluded ear. An occlusion effect value is produced from the difference. The user may read a text passage or vocalize a sound such as /iii/ or /uuu/. The hearing aid may be fitted based on the occlusion effect value.
US 2009/238387 by Arndt et al. discloses a method for actively reducing occlusion. A transducer transmission function, which is defined for the transmission path from the input of a receiver via the auditory canal to the output of a microphone, is subjected to an automatic plausibility check.
US 2009/274314 by Arndt et al. discloses a method for determining a degree of closure in hearing devices. Arndt mentions active occlusion reduction. An effective vent diameter specifies the degree of closure. An interpretation of this value is easily possible by a hearing device acoustician.
WO 2010/083888 A1 by Rung et al. discloses a method for in situ occlusion effect measurement. An external sound pressure of an occluded ear is measured by the microphone of a BTE hearing aid. The sound pressure at the eardrum is measured by a hearing aid receiver.
WO 2012/003855 A1 by Rung discloses a method for measuring the occlusion effect of a hearing aid user. The diameter of a ventilation channel may be increased to reduce the occlusion effect. Leakage between bands is regarded in the measurement.
It is an object of the invention to provide a method for fitting active occlusion control means of a hearing aid device in an easy, precise, flexible, robust, sustainable, effective and/or efficient way. This is especially important because active occlusion control does not only reduce occlusion, but also has side effects. A first side effect is a possible instability of the occlusion control loop. A second side effect is the so called waterbed effect according to which there is not only suppression of occlusion sounds but also amplification of sounds at frequencies below and above the suppression. Hence, what is needed is a good trade-off between wanted and unwanted effects suitable for application in practice.
The object can be at least partially achieved by the method of claim 1. Using a complex frequency-dependent plant transfer function and using an objective frequency-dependent occlusion effect function and/or at least one property of it for determining a compensator filter dataset has the advantage that it allows to adapt an active occlusion control means to the needs of a particular individual in an easy, precise and efficient way.
The method of claim 2 can be advantageous in that predefining compensator filter dataset candidates allows to apply audiological expertise prior to the actual fitting, hence a good fitting can be achieved later with less expertise. Candidates can be predefined with regard to the aspects stability and reliability. Selecting between discrete candidates can be easier, more precise, more efficient and less demanding for a fitter and/or a hearing aid device user than adjusting multiple continuous parameters or even curves. There is not even a need for awareness of the multitude of parameters actually applied.
The method of claim 3 can be advantageous in that by scaling the compensator filter the effect of the filter, and thereby the occlusion control strength, can be adjusted in a precise and easy way. It opens up the possibility to provide a user friendly manual adjustability. Good tradeoffs between wanted and unwanted effects may be found. The occlusion control strength may also be maximized up to the bound given by system stability requirements.
The method of claim 4 can be advantageous in that applying selection criteria to a set of compensator filter candidates allows to select a candidate fully automatically or to reduce the number of candidates to be tested by the user and/or the fitter thereby making the choice of an optimum candidate easier and faster.
The method of claim 5 can be advantageous in that actually trying out the hearing aid with different configurations gives a very good indication which fitting is best in the perception of the user. Letting the user actively participate in the fitting improves the acceptance of its results by the user.
The method of claim 6 can be advantageous in that using a complex frequency-dependent vent effect and/or leakage function for determining a compensator filter dataset allows to adapt active occlusion control means to the needs of a particular individual in an especially precise, optimized and efficient way.
The method of claim 7 can be advantageous in that using a fundamental frequency of a voice of the user for determining a compensator filter dataset allows to adapt active occlusion control means to the needs of a particular individual in an especially precise, optimized and efficient way.
The method of claim 9 can be advantageous in that a benefit assessment allows to prevent waste of effort by individuals involved in such a fitting in cases where there is no potential benefit.
The method of claim 10 can be advantageous in that presenting a recorded real life sound stimulus is perceived by the user of the hearing aid as more pleasant than artificially generated stimuli.
Symbols such as “CA”, “P”, “|OE|” or “{C1, C2, C3 . . . }” in the claims are to be regarded as reference signs if they are presented in parentheses and these parentheses are not part of a formula. Reference signs should not be seen as limiting the extent of the matter protected by the claims. Their sole function is to make the claims easier to understand.
Further embodiments and advantages emerge from the claims and the description referring to the figures.
Below, the invention is described in more detail by referring to the drawings showing exemplified embodiments.
The described embodiments are meant as examples and shall not confine the invention.
The term “canal microphone” in the present document is to be interpreted in a broad manner. It is meant to cover all transducers which are suitable for sensing a sound and/or vibration in the residual canal space, for example conventional microphones, but also optical microphones, acceleration sensors and/or strain gauges. The canal microphone 8 may also be integrated or combined with the receiver 7. Both transducers may simply share a common casing and/or wax protection system and be otherwise separate. However, it is also possible that the two transducers share the same membrane or even a common coil. It is also possible to sense the sound in the residual canal space by one or two vent microphones, the sound inlets of which are arranged in the wall of the vent 10. A directional vent microphone or two vent microphones combined with a special processing may allow to determine which sounds in the vent 10 originate from the residual canal space and not from the outside. The canal microphone 8 may also be combined, complemented and/or enhanced with various further sensors.
The sequence and comprehension of measurements and other steps of this flow diagram is purely exemplary and may be composed and varied in various ways. For example the occlusion effect measurement may be carried out before the plant measurement or may be replaced by an estimation based on already existing data. Further, only a magnitude |OE| or a property of the complex objective frequency-dependent occlusion effect function OE may be determined and/or regarded. The frequency-dependent vent effect and/or leakage function VE may be left out completely or only a magnitude |VE| of it or a cutoff frequency fVE of it may be determined and/or regarded. The fundamental frequency F0 may also be left out completely, or instead a fundamental frequency range {F0min, F0max} may be determined and regarded.
The method steps are presented in the claims in particular sequences. These sequences are exemplary and not mandatory, i.e. the claims are to be interpreted such that they cover also carrying out the same steps, but in other sequences, as far as it is feasible. In particular step B and C of claim 1 may be interchanged.
The shown hearing aid 3 is an ITE or in-the-ear hearing aid for compensating a hearing loss. As already indicated, the invention may also be applied for a hearing protection device such as a Serenity DP+ by Phonak™. The hearing aid device fitted according to the invention may also be a distributed or modular hearing aid device. Such a hearing aid device may have a behind-the-ear module as well as an in-the-ear module. The modules are generally electrically connected to each other. The in-the-ear module preferably comprises both the receiver 7 and the canal microphone 8. It is preferable to arrange both transducers in the canal because sound tubes to modules at other locations would introduce delays in the active occlusion control loop which would interfere with its proper functioning. The in-the-ear module may be a custom ear-piece or a one-size-fits-all dome. The vent 10 in an earpiece of a modular hearing aid or in the main body of an ITE hearing aid has preferably a diameter in a range from 0.6 mm to 1.2 mm, in particular 0.8 mm or 1.0 mm. Larger vents may cause feedback problems and impair sound cleaning features. Smaller vents may be prone to plugging and may not provide sufficient pressure equalization and moisture discharge. If the fitting method is carried out in regard to a plurality of users it is advantageous to use the same vent size each time and to accommodate personal preferences by the selection and scaling of the compensator filter dataset C. A hearing protection device has preferably no vent at all to provide maximum noise attenuation. Even though only one hearing aid 3 is shown a typical user will have two hearing aids. Each of them may be fitted as described in this document, in particular one after the other. However, certain steps may be carried out left and right simultaneously and/or in a synergic manner, as for example the measurement of the complex objective frequency-dependent occlusion effect function OE. The same stimulus presentation may be used for measurements at the left and the right hearing aid. Further, results from left and right may be compared for plausibility checks and/or may be combined for obtaining a higher precision. For example the signals of left and right outside microphones may be averaged or be selectively used depending on which signal is best.
The fitting device 12 is represented in
The system is configured for obtaining the complex frequency-dependent plant transfer function P based on a plant measurement and for using it in the determination of the compensator filter dataset C. The plant measurement is carried out with the hearing aid inserted (in-situ) and preferably, if there is a vent, with an open vent. Only if there is substantial environment noise it may be advantageous to close the vent. However, environment noise may also be dealt with by louder plant stimuli. The user 31 is instructed to remain silent during the measurement. The measurement is similar to a feedback measurement. Hence, it may also be advantageously combined with it, in particular such that both measurements are carried out upon a single user or fitter action. The measurement may in particular be started by the fitter 30 by selecting the option “P” on a mode selector control 32, which in turn may switch the system into a plant measurement mode. For the plant measurement, the receiver 7 may be disconnected physically or logically from the hearing aid sound processing means 5, 6 and 9 and may be connected to a signal 28 provided by a plant stimulus generation and/or playback means 15. Different kinds of stimuli may be used, in particular artificially generated stimuli (AGS), recorded real life sound stimuli (RRS), current environment sound stimuli (CES) and/or stimuli generated based on sounds provided by an external device other than the fitting device 12 (EDS). Artificially generated stimuli may include broadband stimuli, such as pink noise and white noise, as well as tonal stimuli, such as stepped or swept sine or complex multi-sine stimuli. An example of a white noise stimulus is a PRBS stimulus (pseudorandom binary sequence) and in particular an MLS (maximum length sequence) stimulus. Recorded real-life stimuli may include music, nature sounds, such as sounds of a waterfall, voice or own voice of the user. Recorded real life stimuli are perceived by the hearing aid user 31 as being more pleasant and entertaining than artificially generated stimuli. The provision of recorded real life stimuli may be carried out by a hearing aid manufacturer and may comprise the steps of picking up environment sounds in the field with a microphone and storing them on a medium such as a hard disk. Recorded real life stimuli may be enhanced by combining them with other stimuli, in particular artificial ones. This allows for example to assure that all frequencies are sufficiently covered by the stimulus. Current environment sound may be used processed or unprocessed as stimulus. The external device may for example be a hi-fi system. Sounds may be transmitted and/or streamed from the external device to the hearing aid 3 by wire or wirelessly, either directly, or indirectly through the fitting device 12 and/or a streaming device such as an iCOM by Phonak™. The sounds may be used processed or unprocessed as stimuli. Finally the plant stimulus may be any result of filtering and/and combining of stimuli such as for example defined by
FCS=α*AGS+β*RRS+γ*CES+δ*EDS
Wherein α, β, γ and δ may be scalars and/or filters. Plant measurement analysis means 18 may calculate a difference of a logarithmic frequency domain representation of the resulting sound and a logarithmic frequency domain representation of the plant stimulus sound. Alternatively a quotient may be calculated of non-logarithmic representations of these sounds. A frequency analysis method may be used, in particular with tonal stimuli. A correlation method may be used, in particular with broadband stimuli. An adaptive algorithm, e.g. a LMS-algorithm (Least-Mean-Squares), may be used if there is no generated stimulus or if a processed or unprocessed environment sound is used as stimulus. More details about such calculations can be found in textbooks about “system identification”. A plausibility check may be carried out for P, in particular for detecting if a wax protection system of receiver 7 and/or microphone 8 is clogged. Preferably the complex frequency-dependent plant transfer function P is measured directly. However, it is also possible to measure only the magnitude |P| of the plant transfer function P and to estimate a phase function φ=arg(P) e.g. by minimum phase considerations, Hilbert transformation and/or application of a sound propagation delay between receiver and microphone. “Complex” may be defined as “including phase information”. It can be advantageous to subdivide the frequency range of the plant measurement, e.g. at 350 Hz, in order to have more low-frequency measurement points at a given FFT (fast Fourier transformation) size for better determining the low frequency overshoot described further down.
The system is further configured for determining the compensator filter dataset C based on an objective frequency-dependent occlusion effect function and/or based on at least one property of it. The function may be a complex function OE or a magnitude function |OE|. The property may be a peak frequency fOE=OEmax at which the occlusion effect magnitude has its maximum or the relevant maximum. It may be also be a substantial occlusion effect frequency range {f1OE, f2OE} in which the occlusion effect is above a threshold and/or in which the occlusion effect is substantially at its maximum. It may also be a root mean square value OERMS of the objective frequency-dependent occlusion effect function. The property may refer to the full frequency range of OE. However, it may also refer to a certain part of the frequency range.
The objective frequency-dependent occlusion effect function and/or the at least one property of it may be obtained based on a measurement while the voice of the user 31 is active and while there are preferably no other outside sounds. The hearing aid 3 is preferably muted, for example by switching off the receiver. The user's voice may be activated by instructing him or her to speak freely, read a text, repeat a word or sentence, ask a question, sweep a vowel and/or speak different vowels. The measurement may be started by the fitter 30 by selecting the option “OE” on a mode selector control 32, which in turn switches the system into an occlusion measurement mode. The voice of the user may be picked up as a canal sound by canal microphone 8 and as a reference sound by a reference microphone, for example the outside microphone 4, an outside microphone of a further not shown hearing aid or any microphone connected to the fitting device 12. The corresponding signals 26 and 27 are transmitted to the fitting device 12. An open ear gain compensation “OEG” may be applied to the reference sound by compensation means 13 thereby obtaining a compensated outside sound.
Alternatively, an inverse open ear gain compensation “1/OEG” may be applied to the canal sound by compensation means 14 thereby obtaining a compensated canal sound. Occlusion measurement analysis means 16 may calculate a difference of a logarithmic frequency domain representation of the canal sound or, as the case may be, the compensated canal sound and a logarithmic frequency domain representation of the reference sound or, as the case may be, the compensated reference sound. Alternatively a quotient may be calculated of non-logarithmic representations of these sounds. If no OEG compensation has been applied yet, it may still be applied to the resulting difference or quotient, or it may not be applied at all since an OEG is usually not much different from 0 dB in the relevant frequency range below 1 kHz.
Instead of activating and measuring the user's voice, an artificial own voice stimulus may be applied in an occlusion effect measurement. The body of the user may be vibrated by vibrating means. Such means may comprise a body stimulus generator and, connected to it, an electromechanical transducer such as a bone conduction headset. A canal sound resulting from such a vibration in the occluded ear canal is picked up by the canal microphone 8. The signal of the outside microphone 4 is ignored. Instead the signal of the body stimulus generator is used as reference sound. In the further processing the sound in the open ear canal can be estimated by applying a compensation to the reference sound similar to the OEG compensation described above. Accordingly, instead, an inverse compensation may be applied to the canal sound or no compensation may be applied at all. Since the vibration stimulus is reproducible, in contrast to the user's voice, a second, subsequent measurement may be carried out with a probe tube in the canal and without hearing aid 3, thereby obtaining a more precise open ear canal sound as reference sound which needs no compensation. Since the probe tube is already in place, the occluded canal sound may be also measured with the probe tube instead of the canal microphone 8.
In embodiments with a vent 10, the objective frequency-dependent occlusion effect function and/or the at least one property of it may refer to the occlusion with open or closed vent. Hence, in the strict sense OE is either OEVented or OEUnnvented. The same applies accordingly for |OE| and the properties of OE. In many cases it is irrelevant which OE is regarded. OEVented is typically only in the low frequencies affected by the vent effect. In a particular embodiment primarily OEVented is used, and is, if necessary derived from OEUnvented by adding the vent effect. For measuring OEUnvented the vent may be temporarily closed.
The objective frequency-dependent occlusion effect function and/or the at least one property of it may also be entered directly by the fitter 30 or user 31. Alternatively fitter 30 or user 31 may enter data from which it can be derived or which can be used in deriving it. The objective frequency-dependent occlusion effect function and/or the at least one property of it may further be obtained by an estimation based on personal and/or hearing aid device data, in particular the size of a residual space between the earpiece of the hearing aid 3 and the eardrum 1, a middle ear compliance and/or an effective leakage. The residual space depends on the penetration depth of the hearing aid earpiece and the ear canal geometry, which can be determined by an impression or scan. The middle ear compliance may be measured by tympanometry. The effective leakage may depend on the weight and/or material of the hearing aid earpiece. If there is no vent, the effective leakage may be determined based on a real ear occluded gain (REOG) measurement. Finally, in a simplified embodiment one average objective frequency-dependent occlusion effect function may be stored in the fitting device 12 and may be used for all fittings.
The system may also be configured for determining the compensator filter dataset C based on a frequency-dependent vent effect and/or leakage function of an earpiece of the hearing aid 3. The function may be specified by a complex function VE, a magnitude function |VE| or simply by a cutoff frequency fVE of a high-pass filter approximation of such a function. The vent effect information can be manually entered. It can also be measured. It can further be derived from the complex frequency-dependent plant transfer function P, in particular by analyzing a roll-off of the complex frequency-dependent plant transfer function P and/or by applying a low-frequency fitting method of a filter, e.g. 2nd order, in regard to the complex frequency-dependent plant transfer function P. The derivation may be carried out by vent effect and/or leakage derivation means 19. Vent effect is caused by the penetration of sound through the vent 10. Leakage occurs when the hearing aid 3 does not exactly fit the ear canal 2, for example because it is not correctly positioned or the canal has changed since the ear impression for manufacturing the earpiece was taken. Vent effect and leakage may be added to each other for defining a so called “effective vent”. The vent effect and/or leakage function may therefore also be called “effective vent function”.
The system may also be configured for determining the compensator filter dataset C based on a fundamental frequency F0, a fundamental frequency range {F0min, F0max} and/or a fundamental spectrum F0Spectrum of the voice of the user 31. This information can be manually entered. It can also be estimated based on data relating to gender and/or age of the user 31. F0 of males is about 125 Hz, F0 of females about 250 Hz and F0 of children about 440 Hz. F0 and the range {F0min, F0max} can further be measured by sensing the voice of the user by outside microphone 4 and/or canal microphone 8. The hearing aid 3 is preferably muted during the measurement. The measurement can be carried out together with the measurement of the objective frequency-dependent occlusion effect function or properties of it, i.e. the same recorded sound data is used for both, determining F0 and/or the range {F0min, F0max} and determining OE, |OE|, fOE=OEmax, the range {f1OE, f2OE} and/or OERMS. The determination of F0 and the range {F0min, F0max} may be carried out by voice measurement analysis means 17. For measuring the range {F0min, F0max} the user may be instructed to speak in pitch and/or loudness varying way, for example a German speaking user may be instructed to ask a question, at the end of which the pitch is generally higher. F0 and the range {F0min, F0max} may also be acquired in a loudness dependent manner, for example by acquiring the values F0soft, F0mid and F0loud or by acquiring a level dependent function F0L(LdB), wherein LdB is a loudness level in decibels or a loudness level class index. F0 is typically higher for louder voice. In a particular embodiment the range {F0min, F0max} is defined such that it accommodates F0soft, F0mid and optionally F0loud.
The above mentioned measurements are preferably carried out during a fitting session, while there is a data connection between the fitting device 12 and the hearing aid 3 and while the user 31 is in a fitting room or a soundproof both. However, it is also possible that these measurements are carried out in the field, during normal use of the hearing aid 3, at particular times, temporarily and/or in fully continuous manner. A sound situation analysis means may determine which parameter can be measured in a particular situation. For example OE and F0 may be measured in quiet environments, while the user is speaking loudly. P may be measured while the user 31 is quiet, the environment is quiet and loud sounds are presented to him or her by the hearing aid 3, as for example when sounds are streamed from a television with muted loudspeakers. Such measurement results may be used instantaneously for automatically readjusting the compensator filter dataset C in the field. However, they may also be stored in the hearing aid 3 for a later, more controlled use during a fitting session. Accordingly, the fitting device 12 may be configured for reading out such measurement results from the hearing aid 3.
The fitting device 12 may comprise a database 22 with a set of raw compensator filter dataset candidates {C1, C2, C3 . . . }. Raw compensator filter dataset candidates may be represented in different ways as described further below. The term “raw” is used because the datasets are usually further processed and in particular scaled before they are applied in the filter 9 as also described further below. However, the term “raw” in this document is not meant to imply that there must be further processing. In addition, the raw datasets may be a result of a preprocessing, hence they may be only “raw” in respect to a certain stage of the fitting method. The raw candidates may in particular have peak magnitude of 0 dB, which guarantees stability if they are applied unprocessed. The set of candidates is generic in that it is not defined for a particular user. The set of candidates is preferably predefined, for example by a hearing aid manufacturer and/or fitting software provider. It may be distributed together with a fitting software or separately, for example on a compact disk or over the internet. Typically the database remains unchanged after the fitting software has been installed or updated and in particular after the fitting in regard to a particular user has started. The set may comprise one or more candidates. For implementing the concept of choosing between candidates a set of two candidates is sufficient. A reasonable number of candidates may be about fifty. However, memory and processing power of a standard computer may allow thousands or millions of candidates. Therefore it is possible to provide candidates even for very rare user profiles. The predefinition of candidates may be based on statistical and/or empirical data. Hypothetical or real fittings or compensator filter datasets may be determined for typical hearing aid device and user profiles and may be evaluated based on criteria as described further below in regard to the candidate selection. The predefinition of candidates may also comprise the steps of providing a set of base filters {CB1, CB2, CB3 . . . } and a set of modification filters {CM1, CM2, CM3 . . . }. Each base filter can then be combined with each subset of modification filters to determine a candidate. For example candidates may be defined as follows:
C
1
=C
B1
C
7
=C
B2
*C
M1
*C
M3
C
15
=C
B4
*C
M2
*C
M4
Such combinations may be calculated in advance and be provided with the fitting software. However, they may also be calculated at runtime. There may also be separate sets of dataset candidates for different user groups, such as for children, females and males. A lookup table may be used to link user groups with sets.
The fitting device 12 may comprise a candidate selection means 24. In a particular embodiment such a selection may result directly in a compensator filter dataset C for use in the hearing aid 3. However, in a preferred embodiment a preferred raw compensator filter dataset CRAW or set of preferred raw compensator dataset candidates {CA, CB, CC, . . . } is obtained by choosing candidates from the set of raw compensator filter dataset candidates {C1, C2, C3 . . . }.
The preferred candidate or candidates are preferably chosen taking into account the complex frequency-dependent plant transfer function P, the objective frequency-dependent occlusion effect function and/or the at least one property of it, i.e. OE, |OE|, fOE=OEmax, {f1OE, f2OE} and/or OERMS, and optionally the frequency-dependent vent effect and/or leakage function VE, |VE| or a cutoff frequency fVE of a high-pass filter approximation of such a function, as well as the fundamental frequency F0 and/or fundamental frequency range {F0min, F0max}.
The quality of a candidate is preferably assessed by applying a selection criterion K or a set of selection criteria {K1, K2, . . . }. The criterion or at least one criterion of the set of criteria is preferably a property of—or is based on one or more properties of—a complex frequency-dependent candidate specific sensitivity function S and/or a complex frequency-dependent candidate specific occlusion modification function OM. S may be defined by
wherein P is the complex frequency-dependent plant transfer function, Cx is the Xth candidate of the set of raw compensator filter dataset candidates {C1, C2, C3 . . . } and gprov is a provisional scalar scaling factor. An example of S is discussed referring to
OM=VE*S
wherein VE is the complex vent effect and/or leakage function.
The provisional scaling factor gprov is provisional in that it is only used for applying the selection criteria, i.e. used for calculating certain values as shown in the criteria table below. It is a purely theoretical value and is not necessarily applied in the actual hearing aid 3. It must therefore not fulfill stability criteria. There are amongst others the following nonexclusive options:
The following table contains examples of selection criteria:
In the specification of the criteria the expression “matching well” is used for describing the relation between a first and a second measure. If both measures are scalars, e.g. decibel values or frequencies, “matching well” means that the absolute value of their difference is small. If both measures are frequency ranges “matching well” means that the lower and upper bounds match well. If both measures are functions “matching well” may in particular mean that an application of the method of least squares indicates a good matching of the two functions.
When carrying out the task of determining one preferred candidate CRAW or a set of preferred candidates {CA, CB, CC, . . . } by applying a criterion K and by choosing from the available compensator filter dataset candidates {C1, C2, C3 . . . }, a quality indicator may be calculated for each candidate thereby obtaining a set of quality indicators {Q1, Q2, Q3 . . . }. A quality indicator may be a numeric representation of a property defined by a criterion K. Depending on the property small or large values may indicated good quality. It may also be a category such as “poor”, “average”, “good” or the like. The quality indicator Q1 for a candidate C1 and a criterion K, namely “Smallness of Smin”, may be defined by:
Q
1=5 min(C1) or Q1=fQ(Smin(C1))
The function fQ(·) allows to derive quality indicators for properties which reflect not directly an extent of quality, for example if values in a certain range indicate good quality. It also allows to normalize the quality indicators of different criteria, for example if one property is a decibel value and another property is a Hertz value. The important feature of the quality indicator is that it provides a basis for comparing the quality of candidates. The following table shows an example:
The preferred raw compensator filter dataset candidate CRAW according to the example would be C2. A set of two preferred raw compensator filter dataset candidates {CA, CB} according to the example would be {C1, C2}.
As already indicated above, not only one criterion K may be applied, but instead a set of criteria {K1, K2, . . . }. In this case a weighting may be provided for each criterion of the set of criteria thereby obtaining a set of weights {W1, W2, . . . }. The weights allow to regard certain criteria more than others. The following table shows an example with three raw compensator filter dataset candidates and three weighted criteria:
Multi-criteria quality indicators Q1, Q2 and Q3 are calculated for the candidates C1, C2 and C3. The multi-criteria quality indicator Qy for a particular Yth candidate Cy is determined by first calculating criterion-specific quality indicators QCYK1, QCYK2 and QCYK3 for the selection criteria K1, K2, and K3 and then combining these criterion-specific quality indicators in a weighted manner by applying a weighting function fW(·):
Q
Y
=f
W({W1, W2, . . . },{QCYK1, QCYK2, . . . })
In a preferred embodiment the weighting function is linear and applies a weighting factor to each criterion-specific quality indicator, as shown by the following formula:
Q
Y
=W
1
*Q
CYK1
+W
2
*Q
CYK2
+W
3
*Q
CYK3
However, the weighting function fW(·) may also be a polynomial and in particular comprise quadratic terms as shown by the following example:
Q
Y
=W
1*(QCYK1)2+W12*(QCYK1*QCYK2)+W22*(QCYK3)2+
The set of weights {W1, W2, . . . } for the set of selection criteria {K1, K2, . . . } can be obtained by carrying out a subjective evaluation of each candidate of the set of raw compensator filter dataset candidates {C1, C2, C3 . . . } by one or more individuals thereby obtaining a set of subjective evaluation results {R1, R2, R3 . . . }. The evaluation may in particular be carried out based on a scaling to a maximum stable active occlusion control strength and/or based on an adjustable scaling. The weights {W1, W2, . . . } are then set such that a set of multi-criteria quality indicators {Q1, Q2, Q3 . . . } calculated based on the set of weights {W1, W2, . . . } substantially best matches the set of subjective evaluation parameters {R1, R2, R3 . . . }. This may comprise carrying out a regression analysis, a stepwise regression analysis, a discriminant analysis and/or a stepwise discriminant analysis.
As already indicated, the compensator filter datasets {C1, C2, C3 . . . } stored in database 22 are “raw”. Before they are actually applied as occlusion filter dataset C the hearing aid 3 they are scaled by a scaling factor g:
C=C
RAW
*g C=C
1
*g
1
C=C
A
*g
A
The scaling factor g influences the strength of the occlusion control. However, if g is chosen too large, the active occlusion control loop may become unstable. Accordingly, there is a maximum allowable scaling factor gmax. This value depends on the raw compensator filter data set such as CRAW or CA and on the complex frequency-dependent plant transfer function P of the particular individual and should therefore be recalculated if any of these parameters changes. In a preferred embodiment g is not manually adjustable but always set to gmax such that the occlusion control is maximized while keeping the system stable. In another embodiment the scaling factor g and therefore the strength of the occlusion control is adjustable manually by the fitter 30 and/or the user 31, in particular by the strength selector user control 33. The adjustment range is preferably limited such that gmax cannot be exceeded. Further, the g may have a particular initial value g0, which can for example be gmax.
The active occlusion control loop is stable and substantially robust against destabilization if the maximum sensitivity Smax does not exceed a predefined value Sbound. The stability of a system with feedback can be assessed based on a Nyquist plot. A distance between the Nyquist plot and the Nyquist point at (−1, i*0) is a stability criterion. The maximum sensitivity Smax is an indicator for this distance and therefor also a stability criterion. The smaller Smax, the more robust is the system against destabilization. Sbound is typically in the range from 4 to 6 dB, in particular at 5 dB. Preferably the system allows to redefine Sbound, since empirical tests may imply other values. gmax may be calculated based on CRAW, P, Sbound and the following equations:
However, since there is no formula for a direct calculation of gmax it may be advantageous to determine gmax by an iterative method. For example g might be increased in one dB-steps and after each increase Smax is calculated and evaluated.
In a particular implementation of the candidate selection means 24 the user and hearing aid specific data such as P, OE, |OE|, VE, |VE|, fVE, F0, F0min, F0max, age, gender, hearing loss, hearing aid coupling and hearing aid type is mapped to a finite number of categories. The preferred raw compensator filter dataset CRAW or the set of preferred raw compensator dataset candidates {CA, CB, Cc, . . . } is then determined without actually calculating criterion data such as Smin. Instead the candidate or candidates for the determined category are looked up in a lookup table. The lookup table may also be combined with a criterion based evaluation. Both, lookup table and criterion based evaluation may be used in an arbitrary sequence to reduce the number of candidates until a target number of candidates has been reached.
As already indicated the candidate selection means 24 may not only provide a preferred raw compensator filter dataset CRAW but instead also a set of preferred raw compensator dataset candidates {CA, CB, Cc, . . . } which is a subset of the set {C1, C2, C3 . . . } stored in the database. The hearing aid 3 is then temporarily and successively configured based on candidates of this subset. Such a demonstration of candidates may be started by the fitter 30 by selecting the option “ABC” on a mode selector 32, which in turn switches the system into a demonstration mode. In a first trial the compensator filter C may be configured with CA*gA, in a second trial with CB*gB and so forth. A particular candidate may also be demonstrated differently scaled. There may be presentations CA*gA1 and CA*gA2. An additional configuration to be evaluated may be “No AOC”, i.e. without active occlusion control. At least two configurations should be presented, wherein one might be the “No AOC” configuration. However, optimally three to five configurations are presented. The user 31 may be instructed to speak, walk, chew, listen to the fitter 30 speaking or listen to a surround sound system. The user 31 and/or the fitter 30 may actively switch between the configurations by actuating a candidate selector user control 34 or the configurations may be presented automatically one after the other for a certain time and/or until a corresponding evaluation result is entered. Eventually, the fitting device 12 obtains an absolute or relative evaluation information in regard to one or more of the demonstrated configurations from the user 31. The user 31 and/or the fitter 30 may enter such information, in particular by a candidate rating user control 35. Based on the information the system determines which of the candidates CA, CB, Cc is the preferred candidate. The result CRAW or the scaled result CRAW*g is then stored in the non-volatile memory of the hearing aid 3, in particular by selecting the option “NVM” on a mode selector control 32. The hearing aid may be then or thereby switched from the fitting mode back to the operation mode.
The compensator filter dataset C may also be determined without the above mentioned candidates, in particular by a calculation based on the equations:
S=(1+P*C)−1 and S=Starget=OE−1
or the equation:
The fitting method of according to the invention may also be used to determine more than one compensator filter data set, for example for different hearing programs or hearing situations such as a CSp for speech, a CSpN speech in noise, a CC for calm situations and a CM for music or for different loudness levels such as a CS for soft, a CM for medium and a CL for loud. Accordingly, more than one compensator filter data set may be stored in the non-volatile memory of the hearing aid 3.
Once a compensator filter dataset C has been determined the occlusion control compensator filter 9 and the pre-equalizer 6 may be configured based on it, such that it becomes part of an active configuration of the signal processor of the hearing aid 3. This may in particular occur during the above mentioned demonstrations, at the end of the fitting session, when the hearing aid is switched on or to another program, when filter data is transmitted by a signal 29 from the fitting device 12 to the hearing aid 3 and/or when filter data is read from the non-volatile memory of the hearing aid 3.
The compensator filter datasets, such as C, CRAW, C1, C2, C3, CA, CB, CC, CSp, CSpN, CC, CM, CS, CM and CL, may be represented in two substantially different ways:
A first way is named here “coefficient format”. It is a representation as a set of scalar filter coefficients. The filter is preferably time-discrete. Such a set may comprise or consist of coefficients of a numerator polynomial in z, for example {b0, b1, b2, . . . }, and coefficients of a denominator polynomial in z, for example {a1, a2, . . . }. A simple implementation would be a “digital biquad filter”. The coefficients may define a filter of nth order. A representation of C in this format is indicated below by the symbol C[cf].
A second way is named here “function format”. It is a representation as a complex frequency-dependent filter function, also referred to as frequency response. Such a function is preferably frequency discrete such that the function can be described by a complex vector of a predefined dimension. A reasonable tradeoff between accuracy and data size can be achieved by a third octave frequency resolution. A higher resolution function may be filtered to obtain a function having such a resolution. Preferably, the frequency resolution applied in measurements, calculations and/or filter definitions is the same. Accordingly, the complex frequency-dependent functions P, OE, OM, VE, S and C have preferably the same frequency resolution and the corresponding vectors have the same dimension. A representation of C in this format is indicated below by the symbol C[ft].
The “coefficient format” has the advantage that it needs less memory and transmission time than the “function format”. “coefficient format” data may be compressed and/or reduced to a data size of about 75 bytes, i.e. less than 100 bytes, per compensator filter dataset C. The “coefficient format” can easily be converted to the “function format”. Vice versa, it is difficult and not very practical to convert the “function format” to the “coefficient format”. The “function format” is much better suited for assessing the filter quality. The formulas comprising “C” in the present document, such as S=1/(1+P*C) are normally calculated based on the “function format”. An exception is the scaling of a “raw” filter compensator filter dataset with a scaling factor, such as C=CRAW*g, and the additive inversion, such as C′=−C, which can be calculated well in both formats.
In the following it is indicated which format is preferably used in which stage of the fitting process:
The predefinition of raw compensator filter dataset candidates {C1, C2, C3 . . . } is preferably at least partially carried out based on the “function format”, because the predefinition involves most likely filter quality assessments.
The storing of raw compensator filter dataset candidates {C1, C2, C3 . . . } in the database 22 is preferably carried out based on the “coefficient format” because of memory and convertibility considerations. However, the candidates may be stored additionally in the “function format”. This allows to save processing time during the fitting session, because it eliminates the format conversion step.
The quality assessment and candidate selection by the fitting device 12 is preferably carried out based on the “function format”.
The transmission to the hearing aid 3 as well as the signal processing within the hearing aid 3 as well as the storing in the non-volatile memory of the hearing aid 3 is preferably carried out based on the “coefficient format” because of data size considerations and its suitability as a basis for signal processing.
In the candidate selection process, it may be determined that a particular compensator filter dataset C[ff] is a good filter candidate and should be applied as C[cf] in the hearing aid 3. Since it is not practical to calculate C[ct] directly from C[ff] a kind of backtracking is carried out. It is determined which CRAW%[cf] and which g compensator filter dataset C[ff] is based on. C[ct] is then calculated based on the equation C[ct]=CRAW[cf]*g.
The magnitude function |S|dB can be divided into three frequency ranges. In a first range below fmin there is the low frequency overshoot LOS. In a second range between fmin and fmax there is the actual occlusion reduction. In a third range above fmax is the high frequency overshoot HOS, which is typically at 1 to 3 kHz. Occlusion reduction in a particular frequency range is always accompanied by amplification below and above this range. This behavior is called waterbed effect. More formally it is called “Bode's integral theorem”. A large LOS may result in an unpleasant perception of footfall sounds. There is an area A1 between the f-axis and the LOS, an area A2 between the f-axis and negative section of the |S|dB-curve and an area A3 between the f-axis and the HOS. The sum of overshoot areas A1 and A3 is just as large as A2. The area A2 is equal to an absolute value |Sint| of the above defined Sint. The larger A2, the stronger the occlusion reduction. fmin and fmax can be defined as bordering the frequency range where |S|dB is below 0 dB. However, it is to be noted that |S|dB may also be smaller than 0 dB in small or negligible frequency ranges below and above the primary occlusion reduction frequency range. The range between fmin and above fmax may therefore be referred to as the “substantial frequency range where |S|dB is below 0 dB”.
For estimating the quality of a particular sensitivity S, as already indicated above, the various properties of S, and in particular of |S|, may be regarded, for example the shown Smin, Smax, SSmax, Δf, fmin, fmax, fS=Smin, A2, φmax and fφ=φmax and the not shown Sint and SRMS, in particular in relation to further values such as the shown Sthres, Starget, Sbound, ftarget, f1S, f2S, fφTarget, f1φ, f2φ, F0, F0min and F0max and the not shown OE, |OE|, fOE=OEmax, f1OE, f2OE, GERMS, VE, |VE| and fVE.
The parameter Smax may be treated in a special way. Smax should not exceed the upper bound Sbound because otherwise the system may become unstable which results typically in a whistling in the frequency range of the HOS and/or or in a rumbling in the frequency range of the LOS. Therefore the scaling factor g, i.e. g1 or g2, is preferably selected such that Smax is below or at the bound Sbound. The latter applies for the two curves shown in the Bode plot, i.e. g1 is equal to a maximum scaling factor g1max and g2 is equal to a maximum scaling factor g2max.
The parameter Smin may also be treated in a special way. Smin is a good indicator for the strength of the active occlusion control. A threshold Sthres may be used to assure a minimum strength. For rendering a compensator filter dataset “preferred” Smin<Sthres may have to apply. Further, a target value Starget may be defined for Smin to specify a target strength. Smin depends on g1. A target scaling factor gtarget can be defined as being the g1 for which Smin is equal or close to Starget. The scaling factor g1 of the curve shown in the Bode plot is smaller than gtarget. Accordingly Smin is several decibels larger than Starget. Scaling the compensator filter dataset candidate C1 with gtarget would result in a Smax above Sbound. The system would show substantial artifacts and would not be substantially robust against destabilization any more. Therefor the compensator filter dataset C1*gtarget should never be used in the actual hearing aid. However, it may be used for applying selection criteria. If this results in C1 being a preferred candidate, C1 would be scaled with gmax instead of gtarget before being employed and/or evaluated in the actual hearing aid.
The parameter SSmax indicates the maximum steepness of the magnitude curve |S|dB. The maximum steepness is typically in the frequency range towards the HOS. A large maximum steepness should be avoided because it may cause artifacts.
It is to be noted that within this document C is defined such that the filter must be configured with “−C”, i.e. with C having a negative sign, as shown in
It is further to be noted that within this document S is defined such that a smaller magnitudes mean less occlusion. However, S could be just as well defined such that larger magnitudes mean less occlusion. S according to the alternative definition is the multiplicative inverse of S according the primary definition. In case of such an alternative definition statements regarding S and formulas comprising S would have to be adapted accordingly. In particular the formula S=1/(1+P*C) would have to be changed to S=1+P*C. For magnitudes expressed in decibels the additive inverse would have to be used, i.e. “|S|dB” would have to be changed to “−|S|dB”. Claims containing such statements and/or formulas are to be interpreted such that they cover both definitions of S. Their substantial meaning is not changed by an alteration of the definition of S. The same applies in a similar manner for H, E, C, P, OE, OM and VE.
The fitting device may be configured for providing various graphical information about the fitting process and the fitting result, for example Bode plots of complex functions, graphs of spectral functions and bar or pie charts of continuous parameters or ratings. Diagrams may show for example characteristics of P, CRAW, C, S, OE, VE, OM, F0L, F0Spectrum and {R1, R2, R3 . . . }, in particular in relation to each other and/or for different compensator filter datasets in the same diagram. For example the Bode plot of different S shown in
It is estimated that only a small percentage hearing aid users may benefit from an active occlusion control, even if it is optimally fitted. Therefore, in a preferred embodiment of the invention, a benefit assessment is carried out at one or more stages of the method. The subjective benefit can be assessed, as already described above, by comparing one or more configurations having active occlusion control, such as configurations “A”, “B” and “C”, with a configuration “Ø” not having it. Besides of that or instead of that an automatic benefit assessment may be carried out to determine if a substantial benefit can be provided to the user 31 by the active occlusion control feature. If no substantial benefit can be provided the system outputs a corresponding message. The message can for example be an acoustic and/or visual message presented by the fitting device 12. One potential reason for insufficient benefit may be that the user has a relatively strong low frequency hearing loss such that he or she does not perceive occlusion sounds in the first place. Best candidates for occlusion control are individuals having mild hearing losses. Hence, the benefit assessment may comprise the step of analyzing the user's hearing loss or audiogram, in particular by checking if the hearing loss is less than 40 dB at a set of frequencies relevant for occlusion, in particular at 125 Hz, 250 Hz and/or 500 Hz. Further measures may be properties of the complex frequency-dependent plant transfer function P, of the objective frequency-dependent occlusion effect function OE or |OE|, of the frequency-dependent vent effect function VE or |VE| and/or of the fundamental frequency F0 or fundamental frequency range {F0min, F0max} of the user. The feature is useless if there is no substantial or no occlusion effect in the first place, for example if the vent 10 is sufficiently large and if there is no need to reduce its size. Once a compensator filter dataset C has been determined, it is possible to calculate and assess values indicative of the strength of the obtainable occlusion modification, in particular Smin, OMmin, A2, Sint and/or Δf as well as fmin, fmax and/or fS=Smin in relation to F0, F0min and/or F0max. The assessments are preferably carried out as soon as the necessary data becomes available, in particular directly after a corresponding acoustic measurement. Hearing loss data may be available before inserting the hearing aid for the first time, and in particular before selecting a hearing aid.
Although the inventions disclosed herein have been described in terms of the preferred embodiments above, numerous modifications and/or additions to the above-described preferred embodiments would be readily apparent to one skilled in the art. It is intended that the scope of the present inventions extend to all such modifications and/or additions and that the scope of the present inventions is limited solely by the claims set forth below.
Number | Date | Country | Kind |
---|---|---|---|
12159767 | Mar 2012 | EP | regional |