US20070030982A1 - Interference suppression techniques - Google Patents
Interference suppression techniques Download PDFInfo
- Publication number
- US20070030982A1 US20070030982A1 US11/545,256 US54525606A US2007030982A1 US 20070030982 A1 US20070030982 A1 US 20070030982A1 US 54525606 A US54525606 A US 54525606A US 2007030982 A1 US2007030982 A1 US 2007030982A1
- Authority
- US
- United States
- Prior art keywords
- acoustic
- output signal
- sensor signals
- correlation
- components
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 107
- 230000001629 suppression Effects 0.000 title abstract description 4
- 230000006870 function Effects 0.000 claims description 40
- 230000005284 excitation Effects 0.000 claims description 35
- 238000012545 processing Methods 0.000 claims description 32
- 239000011159 matrix material Substances 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 12
- 238000012544 monitoring process Methods 0.000 claims description 11
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 12
- 230000002452 interceptive effect Effects 0.000 abstract description 10
- 238000000605 extraction Methods 0.000 abstract description 2
- 239000000872 buffer Substances 0.000 description 15
- 230000004807 localization Effects 0.000 description 14
- 238000012360 testing method Methods 0.000 description 13
- 238000005070 sampling Methods 0.000 description 9
- 238000013459 approach Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 239000007943 implant Substances 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 210000003128 head Anatomy 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000001143 conditioned effect Effects 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 210000000959 ear middle Anatomy 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000408659 Darpa Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 208000016354 hearing loss disease Diseases 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 239000012723 sample buffer Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R25/00—Deaf-aid sets, i.e. electro-acoustic or electro-mechanical hearing aids; Electric tinnitus maskers providing an auditory perception
- H04R25/40—Arrangements for obtaining a desired directivity characteristic
- H04R25/407—Circuits for combining signals of a plurality of transducers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
- G10L2021/02161—Number of inputs available containing the signal or the noise to be suppressed
- G10L2021/02165—Two microphones, one receiving mainly the noise signal and the other one mainly the speech signal
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2201/00—Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
- H04R2201/40—Details of arrangements for obtaining desired directional characteristic by combining a number of identical transducers covered by H04R1/40 but not provided for in any of its subgroups
- H04R2201/403—Linear arrays of transducers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2225/00—Details of deaf aids covered by H04R25/00, not provided for in any of its subgroups
- H04R2225/43—Signal processing in hearing aids to enhance the speech intelligibility
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2430/00—Signal processing covered by H04R, not provided for in its groups
- H04R2430/20—Processing of the output signals of the acoustic transducers of an array for obtaining a desired directivity characteristic
Definitions
- the present invention is directed to the processing of acoustic signals, and more particularly, but not exclusively, relates to techniques to extract an acoustic signal from a selected source while suppressing interference from other sources using two or more microphones.
- the difficulty of extracting a desired signal in the presence of interfering signals is a long-standing problem confronted by acoustic engineers.
- This problem impacts the design and construction of many kinds of devices such as systems for voice recognition and intelligence gathering.
- Especially troublesome is the separation of desired sound from unwanted sound with hearing, aid devices.
- hearing aid devices do not permit selective amplification of a desired sound when contaminated by noise from a nearby source. This problem is even more severe when the desired sound is a speech signal and the nearby noise is also a speech signal produced by other talkers.
- “noise” refers not only to random or nondeterministic signals, but also to undesired signals and signals interfering with the perception of a desired signal.
- One form of the present invention includes a unique signal processing technique using two or more microphones.
- Other forms include unique devices and methods for processing acoustic signals.
- FIG. 1 is a diagrammatic view of a signal processing system.
- FIG. 2 is a diagram further depicting selected aspects of the system of FIG. 1 .
- FIG. 3 is a flow chart of a routine for operating the system of FIG. 1 .
- FIGS. 4 and 5 depict other embodiments of the present invention corresponding to hearing aid and computer voice recognition applications of the system of FIG. 1 , respectively.
- FIG. 6 is a diagrammatic view of an experimental setup of the system of FIG. 1 .
- FIG. 7 is a graph of magnitude versus time of a target speech signal and two interfering speech signals.
- FIG. 8 is a graph of magnitude versus time of a composite of the speech signals of FIG. 7 before processing, an extracted signal corresponding to the target speech signal of FIG. 7 , and a duplicate of the target speech signal of FIG. 7 for comparison.
- FIG. 9 is a graph providing line plots for regularization factor (M) values of 1.001, 1.005, 1.01, and 1.03 in terms of beamwidth versus frequency.
- FIG. 10 is a flowchart of a procedure that can be performed with the system of FIG. 1 either with or without the routine of FIG. 3 .
- FIGS. 11 and 12 are graphs illustrating the efficacy of the procedure of FIG. 10 .
- FIG. 1 illustrates an acoustic signal processing system 10 of one embodiment of the present invention.
- System 10 is configured to extract a desired acoustic excitation from acoustic source 12 in the presence of interference or noise from other sources, such as acoustic sources 14 , 16 .
- System 10 includes acoustic sensor array 20 .
- sensor array 20 includes a pair of acoustic sensors 22 , 24 within the reception range of sources 12 , 14 , 16 .
- Acoustic sensors 22 , 24 are arranged to detect acoustic excitation from sources 12 , 14 , 16 .
- Sensors 22 , 24 are separated by distance D as illustrated by the like labeled line segment along lateral axis T.
- Lateral axis T is perpendicular to azimuthal axis AZ.
- Midpoint M represents the halfway point along distance D from sensor 22 to sensor 24 .
- Axis AZ intersects midpoint M and acoustic source 12 .
- Axis AZ is designated as a point of reference (zero degrees) for sources 12 , 14 , 16 in the azimuthal plane and for sensors 22 , 24 .
- sources 14 , 16 define azimuthal angles 14 a , 16 a relative to axis AZ of about +22° and ⁇ 65°, respectively.
- acoustic source 12 is at 0° relative to axis AZ.
- the “on axis” alignment of acoustic source 12 with axis AZ selects it as a desired or target source of acoustic excitation to be monitored with system 10 .
- the “off-axis” sources 14 , 16 are treated as noise and suppressed by system 10 , which is explained in more detail hereinafter.
- sensors 22 , 24 can be moved to change the position of axis AZ.
- the designated monitoring direction can be adjusted by changing a direction indicator incorporated in the routine of FIG. 3 as more fully described below. For these operating modes, it should be understood that neither sensor 22 nor 24 needs to be moved to change the designated monitoring direction, and the designated monitoring direction need not be coincident with axis AZ.
- sensors 22 , 24 are omnidirectional dynamic microphones.
- a different type of microphone such as cardioid or hypercardioid variety could be utilized, or such different sensor type can be utilized as would occur to one skilled in the art.
- more or fewer acoustic sources at different azimuths may be present; where the illustrated number and arrangement of sources 12 , 14 , 16 is provided as merely one of many examples. In one such example, a room with several groups of individuals engaged in simultaneous conversation may provide a number of the sources.
- Sensors 22 , 24 are operatively coupled to processing subsystem 30 to process signals received therefrom.
- sensors 22 , 24 are designated as belonging to left channel L and right channel R, respectively.
- the analog time domain signals provided by sensors 22 , 24 to processing subsystem 30 are designated x L (t) and x R (t) for the respective channels L and R.
- Processing subsystem 30 is operable to provide an output signal that suppresses interference from sources 14 , 16 in favor of acoustic excitation detected from the selected acoustic source 12 positioned along axis AZ. This output signal is provided to output device 90 for presentation to a user in the form of an audible or visual signal which can be further processed.
- Processing subsystem 30 includes signal conditioner/filters 32 a and 32 b to filter and condition input signals x L (t) and x R (t) from sensors 22 , 24 ; where t represents time. After signal conditioner/filter 32 a and 32 b , the conditioned signals are input to corresponding Analog-to-Digital (A/D) converters 34 a , 34 b to provide discrete signals x L (Z) and x R (z), for channels L and R, respectively; where z indexes discrete sampling events. The sampling rate f S is selected to provide desired fidelity for a frequency range of interest. Processing subsystem 30 also includes digital circuitry 40 comprising processor 42 and memory 50 . Discrete signals x L (z) and x R (z) are stored in sample buffer 52 of memory 50 in a First-In-First-Out (FIFO) fashion.
- FIFO First-In-First-Out
- Processor 42 can be a software or firmware programmable device, a state logic machine, or a combination of both programmable and dedicated hardware. Furthermore, processor 42 can be comprised of one or more components and can include one or more Central Processing Units (CPUs). In one embodiment, processor 42 is in the form of a digitally programmable, highly integrated semiconductor chip particularly suited for signal processing. In other embodiments, processor 42 may be of a general purpose type or other arrangement as would occur to those skilled in the art.
- CPUs Central Processing Units
- memory 50 can be variously configured as would occur to those skilled in the art.
- Memory 50 can include one or more types of solid-state electronic memory, magnetic memory, or optical memory of the volatile and/or nonvolatile variety.
- memory can be integral with one or more other components of processing subsystem 30 and/or comprised of one or more distinct components.
- Processing subsystem 30 can include any oscillators, control clocks, interfaces, signal conditioners, additional filters, limiters, converters, power supplies, communication ports, or other types of components as would occur to those skilled in the art to implement the present invention.
- subsystem 30 is provided in the form of a single microelectronic device.
- routine 140 is illustrated.
- Digital circuitry 40 is configured to perform routine 140 .
- Processor 42 executes logic to perform at least some the operations of routine 140 .
- this logic can be in the form of software programming instructions, hardware, firmware, or a combination of these.
- the logic can be partially or completely stored on memory 50 and/or provided with one or more other components or devices.
- processing subsystem 30 in the form of signals that are carried by a transmission medium such as a computer network or other wired and/or wireless communication network.
- routine 140 begins with initiation of the A/D sampling and storage of the resulting discrete input samples x L (z) and x R (z) in buffer 52 as previously described. Sampling is performed in parallel with other stages of routine 140 as will become apparent from the following description. Routine 140 proceeds from stage 142 to conditional 144 . Conditional 144 tests whether routine 140 is to continue. If not, routine 140 halts. Otherwise, routine 140 continues with stage 146 . Conditional 144 can correspond to an operator switch, control signal, or power control associated with system 10 (not shown).
- a fast discrete fourier transform (FFT) algorithm is executed on a sequence of samples x L (Z) and x R (z) and stored in buffer 54 for each channel L and R to provide corresponding frequency domain signals X L (k) and X R (k); where k is an index to the discrete frequencies of the FFTs (alternatively referred to as “frequency bins” herein).
- the set of samples x L (z) and x R (z) upon which an FFT is performed can be described in terms of a time duration of the sample data. Typically, for a given sampling rate f S , each FFT is based on more than 100 samples.
- FFT calculations include application of a windowing technique to the sample data.
- a windowing technique utilizes a Hamming window.
- data windowing can be absent or a different type utilized, the FFT can be based on a different sampling approach, and/or a different transform can be employed as would occur to those skilled in the art.
- the resulting spectra X L (k) and X R (k) are stored in FFT buffer 54 of memory 50 . These spectra are generally complex-valued.
- reception of acoustic excitation emanating from a desired direction can be improved by weighting and summing the input signals in a manner arranged to minimize the variance (or equivalently, the energy) of the resulting output signal while under the constraint that signals from the desired direction are output with a predetermined gain.
- Minimizing the variance in this case is equivalent to minimizing the output power of off-axis sources, as related by the optimization goal of relationship (2) that follows: Min W ⁇ E ⁇ ⁇ ⁇ Y ⁇ ( k ) ⁇ 2 ⁇ ( 2 ) where Y(k) is the output signal described in connection with relationship (1).
- e is a two element vector which corresponds to the desired direction.
- the elements of vector e can be selected to monitor along a desired direction that is not coincident with axis AZ.
- vector e becomes complex-valued to represent the appropriate time/phase delays between sensors 22 , 24 that correspond to acoustic excitation off axis AZ.
- vector e operates as the direction indicator previously described.
- alternative embodiments can be arranged to select a desired acoustic excitation source by establishing a different geometric relationship relative to axis AZ.
- the direction for monitoring a desired source can be disposed at a nonzero azimuthal angle relative to axis AZ.
- Procedure 520 described in connection with the flowchart of FIG. 10 hereinafter provides an example of a localization/tracking routine that can be used in conjunction with routine 140 to steer vector e.
- e is the vector associated with the desired reception direction
- R(k) is the correlation matrix for the k th frequency
- W(k) is the optimal weight vector for the k th frequency
- the superscript “ ⁇ 1” denotes the matrix inverse.
- the correlation matrix R(k) can be estimated from spectral data obtained via a number “F” of fast discrete Fourier transforms (FFTs) calculated over a relevant time interval.
- the correlation matrix for the k th frequency, R(k) is expressed by the following relationship (5):
- X ll (k), X lr (k), X rl (k), and X rr (k) represent the weighted sums for purposes of compact expression. It should be appreciated that the elements of the R(k) matrix are nonlinear, and therefore Y(k) is a nonlinear function of the inputs.
- stage 148 spectra X l (k) and X r (k) previously stored in buffer 54 are read from memory 50 in a First-In-First-Out (FIFO) sequence. Routine 140 then proceeds to stage 150 . In stage 150 , multiplier weights W L (k), W R (k) are applied to X l (k) and X r (k), respectively, in accordance with the relationship (1) for each frequency k to provide the output spectra Y(k). Routine 140 continues with stage 152 which performs an Inverse Fast Fourier Transform (IFFT) to change the Y(k) FFT determined in stage 150 into a discrete time domain form designated y(z).
- IFFT Inverse Fast Fourier Transform
- a Digital-to-Analog (D/A) conversion is performed with D/A converter 84 ( FIG. 2 ) to provide an analog output signal y(t).
- D/A converter 84 FIG. 2
- correspondence between Y(k) FFTs and output sample y(z) can vary. In one embodiment, there is one Y(k) FFT output for every y(z), providing a one-to-one correspondence. In another embodiment, there may be one Y(k) FFT for every 16 output samples y(z) desired, in which case the extra samples can be obtained from available Y(k) FFTs. In still other embodiments, a different correspondence may be established.
- signal y(t) is input to signal conditioner/filter 86 .
- Conditioner/filter 86 provides the conditioned signal to output device 90 .
- output device 90 includes an amplifier 92 and audio output device 94 .
- Device 94 may be a loudspeaker, hearing aid receiver output, or other device as would occur to those skilled in the art.
- system 10 processes a binaural input to produce an monaural output. In some embodiments, this output could be further processed to provide multiple outputs. In one hearing aid application example, two outputs are provided that deliver generally the same sound to each ear of a user. In another hearing aid application, the sound provided to each ear selectively differs in terms of intensity and/or timing to account for differences in the orientation of the sound source to each sensor 22 , 24 , improving sound perception.
- conditional 156 tests whether a desired time interval has passed since the last calculation of vector W(k). If this time period has not lapsed, then control flows to stage 158 to shift buffers 52 , 54 to process the next group of signals. From stage 158 , processing loop 160 closes, returning to conditional 144 . Provided conditional 144 remains true, stage 146 is repeated for the next group of samples of x L (z) and x R (Z) to determine the next pair of X L (k) and X R (k) FFTs for storage in buffer 54 .
- stages 148 , 150 , 152 , 154 are repeated to process previously stored X l (k) and X r (k) FFTs to determine the next Y(k) FFT and correspondingly generate a continuous y(t).
- buffers 52 , 54 are periodically shifted in stage 158 with each repetition of loop 160 until either routine 140 halts as tested by conditional 144 or the time period of conditional 156 has lapsed.
- routine 140 proceeds from the affirmative branch of conditional 156 to calculate the correlation matrix R(k) in accordance with relationship (5) in stage 162 . From this new correlation matrix R(k), an updated vector W(k) is determined in accordance with relationship (4) in stage 164 . From stage 164 , update loop 170 continues with stage 158 previously described, and processing loop 160 is re-entered until routine 140 halts per conditional 144 or the time for another recalculation of vector W(k) arrives.
- the time period tested in conditional 156 may be measured in terms of the number of times loop 160 is repeated, the number of FFTs or samples generated between updates, and the like. Alternatively, the period between updates can be dynamically adjusted based on feedback from an operator or monitoring device (not shown).
- routine 140 When routine 140 initially starts, earlier stored data is not generally available. Accordingly, appropriate seed values may be stored in buffers 52 , 54 in support of initial processing. In other embodiments, a greater number of acoustic sensors can be included in array 20 and routine 140 can be adjusted accordingly.
- vector e may be varied with frequency to change the desired monitoring direction or look-direction and correspondingly steer the array.
- H ⁇ ( W ) 1 2 ⁇ W ⁇ ( k ) H ⁇ R ⁇ ( k ) ⁇ W ⁇ ( k ) + ⁇ ⁇ ( e H ⁇ W ⁇ ( k ) - 1 ) ( 12 ) where the factor of one half (1 ⁇ 2) is introduced to simplify later math.
- relationship (5) may be expressed more compactly by absorbing the weighted sums into the terms X ll , X lr , X rl and X rr , and then renaming them as components of the correlation matrix R(k) per relationship (18):
- Its inverse may be expressed in relationship (19) as:
- R ⁇ ( k ) - 1 [ R 22 - R 12 - R 21 R 11 ] * 1 det ⁇ ( R ⁇ ( k ) ) ( 19 ) where det( ) is the determinant operator.
- a modified approach can be utilized in applications where gain differences between sensors of array 20 are negligible.
- an additional constraint is utilized.
- the desired weights satisfy relationship (25) as follows:
- the variance minimization goal and unity gain constraint for this alternative approach correspond to the following relationships (26) and (27), respectively:
- e H ⁇ [ 1 2 + Im ⁇ [ w 1 ] 1 2 + Im ⁇ [ w 2 ] ] 1 ( 27 )
- the weights determined in accordance with relationship (29) can be used in place of those determined with relationships (22), (23), and (24); where R 11 , R 12 , R 21 , R 22 , are the same as those described in connection with relationship (18). Under appropriate conditions, this substitution typically provides comparable results with more efficient computation.
- relationship (29) it is generally desirable for the target speech or other acoustic signal to originate from the on-axis direction and for the sensors to be matched to one another or to otherwise compensate for inter-sensor differences in gain.
- localization information about sources of interest in each frequency band can be utilized to steer sensor array 20 in conjunction with the relationship (29) approach. This information can be provided in accordance with procedure 520 more fully described hereinafter in connection with the flowchart of FIG. 10 .
- regularization factor M typically is slightly greater than 1.00 to limit the magnitude of the weights in the event that the correlation matrix R(k) is, or is close to being, singular, and therefore noninvertable. This occurs, for example, when time-domain input signals are exactly the same for F consecutive FFT calculations. It has been found that this form of regularization also can improve the perceived sound quality by reducing or eliminating processing artifacts common to time-domain beamformers.
- regularization factor M is a constant. In other embodiments, regularization factor M can be used to adjust or otherwise control the array beamwidth, or the angular range at which a sound of a particular frequency can impinge on the array relative to axis AZ and be processed by routine 140 without significant attenuation.
- Beamwidth —3 dB defines a beamwidth that attenuates the signal of interest by a relative amount less than or equal to three decibels (dB). It should be understood that a different attenuation threshold can be selected to define beamwidth in other embodiments of the present invention.
- FIG. 9 provides a graph of four lines of different patterns to represent constant values 1.001, 1.005, 1.01, and 1.03, of regularization factor M, respectively, in terms of beamwidth versus frequency.
- routine 140 Per relationship (30), as frequency increases, beamwidth decreases; and as regularization factor M increases, the beamwidth increases. Accordingly, in one alternative embodiment of routine 140 , regularization factor M is increased as a function of frequency to provide a more uniform beamwidth across a desired range of frequencies. In another embodiment of routine 140 , M is alternatively or additionally varied as a function of time. For example, if little interference is present in the input signals in certain frequency bands, the regularization factor M can be increased in those bands. It has been found that beamwidth increases in frequency bands with low or no inference commonly provide a better subjective sound quality by limiting the magnitude of the weights used in relationships ( 22 ), (23), and/or (29).
- this improvement can be complemented by decreasing regularization factor M for frequency bands that contain interference above a selected threshold. It has been found that such decreases commonly provide more accurate filtering, and better cancellation of interference.
- regularization factor M varies in accordance with an adaptive function based on frequency-band-specific interference.
- regularization factor M varies in accordance with one or more other relationships as would occur to those skilled in the art.
- system 210 includes eyeglasses G and acoustic sensors 22 and 24 .
- Acoustic sensors 22 and 24 are fixed to eyeglasses G in this embodiment and spaced apart from one another, and are operatively coupled to processor 30 .
- Processor 30 is operatively coupled to output device 190 .
- Output device 190 is in the form of a hearing aid earphone and is positioned in ear E of the user to provide a corresponding audio signal.
- processor 30 is configured to perform routine 140 or its variants with the output signal y(t) being provided to output device 190 instead of output device 90 of FIG. 2 .
- an additional output device 190 can be coupled to processor 30 to provide sound to another ear (not shown).
- This arrangement defines axis AZ to be perpendicular to the view plane of FIG. 4 as designated by the like labeled cross-hairs located generally midway between sensors 22 and 24 .
- the user wearing eyeglasses G can selectively receive an acoustic signal by aligning the corresponding source with a designated direction, such as axis AZ.
- a designated direction such as axis AZ.
- sources from other directions are attenuated.
- the wearer may select a different signal by realigning axis AZ with another desired sound source and correspondingly suppress a different set of off-axis sources.
- system 210 can be configured to operate with a reception direction that is not coincident with axis AZ.
- Processor 30 and output device 190 may be separate units (as depicted) or included in a common unit worn in the ear.
- the coupling between processor 30 and output device 190 may be an electrical cable or a wireless transmission.
- sensors 22 , 24 and processor 30 are remotely located relative to each other and are configured to broadcast to one or more output devices 190 situated in the ear E via a radio frequency transmission.
- sensors 22 , 24 are sized and shaped to fit in the ear of a listener, and the processor algorithms are adjusted to account for shadowing caused by the head, torso, and pinnae.
- This adjustment may be provided by deriving a Head-Related-Transfer-Function (HRTF) specific to the listener or from a population average using techniques known to those skilled in the art. This function is then used to provide appropriate weightings of the output signals that compensate for shadowing.
- HRTF Head-Related-Transfer-Function
- a cochlear implant is typically disposed in a middle ear passage of a user and is configured to provide electrical stimulation signals along the middle ear in a standard manner.
- the implant can include some or all of processing subsystem 30 to operate in accordance with the teachings of the present invention.
- one or more external modules include some or all of subsystem 30 .
- a sensor array associated with a hearing aid system based on a cochlear implant is worn externally, being arranged to communicate with the implant through wires, cables, and/or by using a wireless technique.
- FIG. 5 shows a voice input device 310 employing the present invention as a front end speech enhancement device for a voice recognition routine for personal computer C; where like reference numerals refer to like features.
- Device 310 includes acoustic sensors 22 , 24 spaced apart from each other in a predetermined relationship. Sensors 22 , 24 are operatively coupled to processor 330 within computer C.
- Processor . 330 provides an output signal for internal use or responsive reply via speakers 394 a , 394 b and/or visual display 396 ; and is arranged to process vocal inputs from sensors 22 , 24 in accordance with routine 140 or its variants.
- a user of computer C aligns with a predetermined axis to deliver voice inputs to device 310 .
- device 310 changes its monitoring direction based on feedback from an operator and/or automatically selects a monitoring direction based on the location of the most intense sound source over a selected period of time.
- the source localization/tracking ability provided by procedure 520 as illustrated in the flowchart of FIG. 10 can be utilized.
- the directionally selective speech processing features of the present invention are utilized to enhance performance of a hands-free telephone, audio surveillance device, or other audio system.
- the directional orientation of a sensor array relative to the target acoustic source changes. Without accounting for such changes, attenuation of the target signal can result. This situation can arise, for example, when a binaural hearing aid wearer turns his or her head so that he or she is not aligned properly with the target source, and the hearing aid does not otherwise account for this misalignment. It has been found that attenuation due to misalignment can be reduced by localizing and/or tracking one or more acoustic sources of interests.
- the flowchart of FIG. 10 illustrates procedure 520 to track and/or localize a desired acoustic source relative to a reference.
- Procedure 520 can be utilized for a hearing aid or in other applications such as a voice input device, a hands-free telephone, audio surveillance equipment, and the like—either in conjunction with or independent of previously described embodiments.
- Procedure 520 is described as follows in terms of an implementation with system 10 of FIG. 1 .
- processing system 30 can include logic to execute one or more stages and/or conditionals of procedure 520 as appropriate.
- a different arrangement can be used to implement procedure 520 as would occur to one skilled in the art.
- Procedure 520 starts with AID conversion in stage 522 in a manner like that described for stage 142 of routine 140 . From stage 522 , procedure 520 continues with stage 524 to transform the digital data obtained from stage 522 , such that “G” number of FFTs are provided each with “N” number of FFT frequency bins. Stages 522 and 524 can be executed in an ongoing fashion, buffering the results periodically for later access by other operations of procedure 520 in a parallel, pipelined, sequence-specific, or different manner as would occur to one skilled in the art.
- procedure 520 continues by entering frequency bin processing loop 530 and FFT processing loop 540 .
- loop 530 is nested within loop 540 .
- Loops 530 and 540 begin with stage 532 .
- routine 520 determines the difference in phase between channels L and R for the current frequency bin k of the FFT g, converts the phase difference to a difference in distance, and determines the ratio x(g,k) of this distance difference to the sensor spacing D in accordance with relationship (35).
- Ratio x(g,k) is used to find the signal angle of arrival ⁇ x , rounded to the nearest degree, in accordance with relationship (34).
- loop 540 closes, returning to stage 532 to process the new g and k combination. If conditional test 542 is affirmative, then all N bins for each of the G number of FFTs have been processed, and loops 530 and 540 are exited.
- the elements of array P( ⁇ ) provide a measure of the likelihood that an acoustic source corresponds to a given direction (azimuth in this case). By examining P( ⁇ ), an estimate of the spatial distribution of acoustic sources at a given moment in time is obtained. From loops 530 , 540 , procedure 520 continues with stage 550 .
- the PEAKS operation of relationship (36) can use a number of peak-finding algorithms to locate maxima of the data, including optionally smoothing the data and other operations.
- procedure 520 continues with stage 552 in which one or more peaks are selected.
- the peak closest to the on-axis direction typically corresponds to the desired source.
- Procedure 520 continues from stage 554 to conditional 560 . Conditional 560 tests whether procedure 520 is to continue or not. If the conditional 560 test is true, procedure 520 loops back to stage 522 . If the conditional 560 test is false, procedure 520 halts.
- the peak closest to axis AZ is selected, and utilized to steer array 20 by adjusting steering vector e.
- vector e is modified for each frequency bin k so that it corresponds to the closest peak direction ⁇ tar .
- routine 140 the modified steering vector e of relationship (38) can be substituted into relationship (4) of routine 140 to extract a signal originating from direction ⁇ tar .
- procedure 520 can be integrated with routine 140 to perform localization with the same FFT data.
- the AID conversion of stage 142 can be used to provide digital data for subsequent processing by both routine 140 and procedure 520 .
- some or all of the FFTs obtained for routine 140 can be used to provide the G FFTs for procedure 520 .
- beamwidth modifications can be combined with procedure 520 in various applications either with or without routine 140 .
- the indexed execution of loops 530 and 540 can be at least partially performed in parallel with or without routine 140 .
- one or more transformation techniques are utilized in addition to or as an alternative to fourier transforms in one or more forms of the invention previously described.
- wavelet transform which mathematically breaks up the time-domain waveform into many simple waveforms, which may vary widely in shape.
- wavelet basis functions are similarly shaped signals with logarithmically spaced frequencies. As frequency rises, the basis functions become shorter in time duration with the inverse of frequency.
- wavelet transforms represent the processed signal with several different components that retain amplitude and phase information. Accordingly, routine 140 and/or routine 520 can be adapted to use such alternative or additional transformation techniques.
- any signal transform components that provide amplitude and/or phase information about different parts of an input signal and have a corresponding inverse transformation can be applied in addition to or in place of FFTs.
- Routine 140 and the variations previously described generally adapt more quickly to signal changes than conventional time-domain iterative-adaptive schemes.
- the F number of FFTs associated with correlation matrix R(k) may provide a more desirable result if it is not constant for all signals (alternatively designated the correlation length F).
- the correlation length F may be designated the correlation length F.
- a smaller correlation length F is best for rapidly changing input signals, while a larger correlation length F is best for slowly changing input signals.
- a varying correlation length F can be implemented in a number of ways.
- filter weights are determined using different parts of the frequency-domain data stored in the correlation buffers. For buffer storage in the order of the time they are obtained (First-In, First-Out (FIFO) storage), the first half of the correlation buffer contains data obtained from the first half of the subject time interval and the second half of the buffer contains data from the second half of this time interval.
- FIFO First-In, First-Out
- filter coefficients can be obtained using both R 1 (k) and R 2 (k). If the weights differ significantly for some frequency band k between R 1 (k) and R 2 (k), a significant change in signal statistics may be indicated. This change can be quantified by examining the change in one weight through determining the magnitude and phase change of the weight and then using these quantities in a function to select the appropriate correlation length F.
- ⁇ M ( k )
- F (k ) The correlation length F for some frequency bin k is now denoted as F(k).
- ⁇ A(k) and ⁇ M(k) increase, indicating a change in the data, the output of the function decreases.
- F(k) is limited between c min (k) and c min (k), so that the correlation length can vary only within a predetermined range. It should also be understood that F(k) may take different forms, such as a nonlinear function or a function of other measures of the input signals.
- the adaptive correlation length process described in connection with relationships (39)-(44) can be incorporated into the correlation matrix stage 162 and weight determination stage 164 for use in a hearing aid, such as that described in connection with FIG. 4 , or other applications like surveillance equipment, voice recognition systems, and hands-free telephones, just to name a few.
- Logic of processing subsystem 30 can be adjusted as appropriate to provide for this incorporation.
- the adaptive correlation length process can be utilized with the relationship (29) approach to weight computation, the dynamic beamwidth regularization factor variation described in connection with relationship (30) and FIG. 9 , the localization/tracking procedure 520 , alternative transformation embodiments, and/or such different embodiments or variations of routine 140 as would occur to one skilled in the art.
- the application of adaptive correlation length can be operator selected and/or automatically applied based on one or more measured parameters as would occur to those skilled in the art.
- One further embodiment includes: detecting acoustic excitation with a number of acoustic sensors that provide a number of sensor signals; establishing a set of frequency components for each of the sensor signals; and determining an output signal representative of the acoustic excitation from a designated direction. This determination includes weighting the set of frequency components for each of the sensor signals to reduce variance of the output signal and provide a predefined gain of the acoustic excitation from the designated direction.
- a hearing aid in another embodiment, includes a number of acoustic sensors in the presence of multiple acoustic sources that provide a corresponding number of sensor signals. A selected one of the acoustic sources is monitored. An output signal representative of the selected one of the acoustic sources is generated. This output signal is a weighted combination of the sensor signals that is calculated to minimize variance of the output signal.
- a still further embodiment includes: operating a voice input device including a number of acoustic sensors that provide a corresponding number of sensor signals; determining a set of frequency components for each of the sensor signals; and generating an output signal representative of acoustic excitation from a designated direction.
- This output signal is a weighted combination of the set of frequency components for each of the sensor signals calculated to minimize variance of the output signal.
- a further embodiment includes an acoustic sensor array operable to detect acoustic excitation that includes two or more acoustic sensors each operable to provide a respective one of a number of sensor signals. Also included is a processor to determine a set of frequency components for each of the sensor signals and generate an output signal representative of the acoustic excitation from a designated direction. This output signal is calculated from a weighted combination of the set of frequency components for each of the sensor signals to reduce variance of the output signal subject to a gain constraint for the acoustic excitation from the designated direction.
- a further embodiment includes: detecting acoustic excitation with a number of acoustic sensors that provide a corresponding number of signals; establishing a number of signal transform components for each of these signals; and determining an output signal representative of acoustic excitation from a designated direction.
- the signal transform components can be of the frequency domain type.
- a determination of the output signal can include weighting the components to reduce variance of the output signal and provide a predefined gain of the acoustic excitation from the designated direction.
- a hearing aid is operated that includes a number of acoustic sensors. These sensors provide a corresponding number of sensor signals. A direction is selected to monitor for acoustic excitation with the hearing aid. A set of signal transform components for each of the sensor signals is determined and a number of weight values are calculated as a function of a correlation of these components, an adjustment factor, and the selected direction. The signal transform components are weighted with the weight values to provide an output signal representative of the acoustic excitation emanating from the direction.
- the adjustment factor can be directed to correlation length or a beamwidth control parameter just to name a few examples.
- a hearing aid is operated that includes a number of acoustic sensors to provide a corresponding number of sensor signals.
- a set of signal transform components are provided for each of the sensor signals and a number of weight values are calculated as a function of a correlation of the transform components for each of a number of different frequencies. This calculation includes applying a first beamwidth control value for a first one of the frequencies and a second beamwidth control value for a second one of the frequencies that is different than the first value.
- the signal transform components are weighted with the weight values to provide an output signal.
- acoustic sensors of the hearing aid provide corresponding signals that are represented by a plurality of signal transform components.
- a first set of weight values are calculated as a function of a first correlation of a first number of these components that correspond to a first correlation length.
- a second set of weight values are calculated as a function of a second correlation of a second number of these components that correspond to a second correlation length different than the first correlation length.
- An output signal is generated as a function of the first and second weight values.
- acoustic excitation is detected with a number of sensors that provide a corresponding number of sensor signals.
- a set of signal transform components is determined for each of these signals.
- At least one acoustic source is localized as a function of the transform components.
- the location of one or more acoustic sources can be tracked relative to a reference.
- an output signal can be provided as a function of the location of the acoustic source determined by localization and/or tracking, and a correlation of the transform components.
- FIG. 6 illustrates the experimental set-up for testing the present invention.
- the algorithm has been tested with real recorded speech signals, played through loudspeakers at different spatial locations relative to the receiving microphones in an anechoic chamber.
- a pair of microphones 422 , 424 (Sennheiser MKE 2-60) with an inter-microphone distance D of 15 cm, were situated in a listening room to serve as sensors 22 , 24 .
- Various loudspeakers were placed at a distance of about 3 feet from the midpoint M of the microphones 422 , 424 corresponding to different azimuths.
- One loudspeaker was situated in front of the microphones that intersected axis AZ to broadcast a target speech signal (corresponding to source 12 of FIG. 2 ).
- Several loudspeakers were used to broadcast words or sentences that interfere with the listening of target speech from different azimuths.
- Microphones 422 , 424 were each operatively coupled to a Mic-to-Line preamp 432 (Shure FP-11).
- the output of each preamp 432 was provided to a dual channel volume control 434 provided in the form of an audio preamplifier (Adcom GTP-5511).
- the output of volume control 434 was fed into A/D converters of a Digital Signal Processor (DSP) development board 440 provided by Texas Instruments (model number TI-C6201 DSP Evaluation Module (EVM)).
- DSP Digital Signal Processor
- Development board 440 includes a fixed-point DSP chip (model number TMS320C62) running at a clock speed of 133 MHz with a peak throughput of 1064 MIPS (millions of instructions per second).
- This DSP executed software configured to implement routine 140 in real-time.
- the sampling frequency for these experiments was about 8 kHz with 16-bit A/D and D/A conversion.
- the FFT length was 256 samples, with an FFT calculated every 16 samples.
- the computation leading to the characterization and extraction of the desired signal was found to introduce a delay in a range of about 10-20 milliseconds between the input and output.
- FIGS. 7 and 8 each depict traces of three acoustic signals of approximately the same energy.
- the target signal trace is shown between two interfering signals traces broadcast from azimuths 22° and ⁇ 65°, respectively. These azimuths are depicted in FIG. 1 .
- the target sound is a prerecorded voice from a female (second trace), and is emitted by the loudspeaker located near 0°.
- One interfering sound is provided by a female talker (top trace of FIG. 7 ) and the other interfering sound is provided by a male talker (bottom trace of FIG. 7 ).
- the phrase repeated by the corresponding talker is reproduced above the respective trace.
- Routine 140 as embodied in board 440 , processed this contaminated signal with high fidelity and extracted the target signal by markedly suppressing the interfering sounds. Accordingly, intelligibility of the target signal was restored as illustrated by the second trace. The intelligibility was significantly improved and the extracted signal resembled the original target signal reproduced for comparative purposes as the bottom trace of FIG. 8 .
- FIGS. 11 and 12 are computer generated image graphs of simulated results for procedure 520 . These graphs plot localization results of azimuth in degrees versus time in seconds. The localization results are plotted as shading, where the darker the shading, the stronger the localization result at that angle and time. Such simulations are accepted by those skilled in the art to indicate efficacy of this type of procedure.
- FIG. 11 illustrates the localization results when the target acoustic source is generally stationary with a direction of about 10° off-axis.
- the actual direction of the target is indicated by a solid black line.
- FIG. 12 illustrates the localization results for a target with a direction that is changing sinusoidally between +10° and ⁇ 10°, as might be the case for a hearing aid wearer shaking his or her head.
- the actual location of the source is again indicated by a solid black line.
- the localization technique of procedure 520 accurately indicates the location of the target source in both cases because the darker shading matches closely to the actual location lines. Because the target source is not always producing a signal free of interference overlap, localization results may be strong only at certain times. In FIG. 12 , these stronger intervals can be noted at about 0.2, 0.7, 0.9, 1.25, 1.7, and 2.0 seconds. It should be understood that the target location can be readily estimated between such times.
Landscapes
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Neurosurgery (AREA)
- Circuit For Audible Band Transducer (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Amplifiers (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
- Transition And Organic Metals Composition Catalysts For Addition Polymerization (AREA)
Abstract
System (10) is disclosed including an acoustic sensor array (20) coupled to processor (42). System (10) processes inputs from array (20) to extract a desired acoustic signal through the suppression of interfering signals. The extraction/suppression is performed by modifying the array (20) inputs in the frequency domain with weights selected to minimize variance of the resulting output signal while maintaining unity gain of signals received in the direction of the desired acoustic signal. System (10) may be utilized in hearing aids, voice input devices, surveillance devices, and other applications.
Description
- The present application is a continuation-in-part of U.S. patent application Ser. No. No. 09/568,430 filed on May 10, 2000, and is related to: U.S. patent application Ser. No. No. 09/193,058 filed on 16 Nov. 1998, which is a continuation-in-part of U.S. patent application Ser. No. 08/666,757 filed Jun. 19, 1996 (now U.S. Pat. No. 6,222,927 B1);
- U.S. patent application Ser. No. 09/568,435 filed on May 10, 2000; and U.S. patent application Ser. No. 09/805,233 filed on Mar. 13, 2001, which is a continuation of International Patent Application Number PCT/US99/26965, all of which are hereby incorporated by reference.
- The U.S. Government has a paid-up license in this invention and the right in limited circumstances to require the patent owner to license others on reasonable terms as provided for by DARPA Contract No. ARMY SUNY240-6762A and National Institutes of Health Contract No. R21DC04840.
- The present invention is directed to the processing of acoustic signals, and more particularly, but not exclusively, relates to techniques to extract an acoustic signal from a selected source while suppressing interference from other sources using two or more microphones.
- The difficulty of extracting a desired signal in the presence of interfering signals is a long-standing problem confronted by acoustic engineers. This problem impacts the design and construction of many kinds of devices such as systems for voice recognition and intelligence gathering. Especially troublesome is the separation of desired sound from unwanted sound with hearing, aid devices. Generally, hearing aid devices do not permit selective amplification of a desired sound when contaminated by noise from a nearby source. This problem is even more severe when the desired sound is a speech signal and the nearby noise is also a speech signal produced by other talkers. As used herein, “noise” refers not only to random or nondeterministic signals, but also to undesired signals and signals interfering with the perception of a desired signal.
- One form of the present invention includes a unique signal processing technique using two or more microphones. Other forms include unique devices and methods for processing acoustic signals.
- Further embodiments, objects, features, aspects, benefits, forms, and advantages of the present invention shall become apparent from the detailed drawings and descriptions provided herein.
-
FIG. 1 is a diagrammatic view of a signal processing system. -
FIG. 2 is a diagram further depicting selected aspects of the system ofFIG. 1 . -
FIG. 3 is a flow chart of a routine for operating the system ofFIG. 1 . -
FIGS. 4 and 5 depict other embodiments of the present invention corresponding to hearing aid and computer voice recognition applications of the system ofFIG. 1 , respectively. -
FIG. 6 is a diagrammatic view of an experimental setup of the system ofFIG. 1 . -
FIG. 7 is a graph of magnitude versus time of a target speech signal and two interfering speech signals. -
FIG. 8 is a graph of magnitude versus time of a composite of the speech signals ofFIG. 7 before processing, an extracted signal corresponding to the target speech signal ofFIG. 7 , and a duplicate of the target speech signal ofFIG. 7 for comparison. -
FIG. 9 is a graph providing line plots for regularization factor (M) values of 1.001, 1.005, 1.01, and 1.03 in terms of beamwidth versus frequency. -
FIG. 10 is a flowchart of a procedure that can be performed with the system ofFIG. 1 either with or without the routine ofFIG. 3 . -
FIGS. 11 and 12 are graphs illustrating the efficacy of the procedure ofFIG. 10 . - While the present invention can take many different forms, for the purpose of promoting an understanding of the principles of the invention, reference will now be made to the embodiments illustrated in the drawings and specific language will be used to describe the same. It will nevertheless be understood that no limitation of the scope of the invention is thereby intended. Any alterations and further modifications of the described embodiments, and any further applications of the principles of the invention as described herein are contemplated as would normally occur to one skilled in the art to which the invention relates.
-
FIG. 1 illustrates an acousticsignal processing system 10 of one embodiment of the present invention.System 10 is configured to extract a desired acoustic excitation fromacoustic source 12 in the presence of interference or noise from other sources, such asacoustic sources System 10 includesacoustic sensor array 20. For the example illustrated,sensor array 20 includes a pair ofacoustic sensors sources Acoustic sensors sources -
Sensors sensor 22 tosensor 24. Axis AZ intersects midpoint M andacoustic source 12. Axis AZ is designated as a point of reference (zero degrees) forsources sensors sources azimuthal angles acoustic source 12 is at 0° relative to axis AZ. In one mode of operation ofsystem 10, the “on axis” alignment ofacoustic source 12 with axis AZ selects it as a desired or target source of acoustic excitation to be monitored withsystem 10. In contrast, the “off-axis”sources system 10, which is explained in more detail hereinafter. To adjust the direction being monitored,sensors FIG. 3 as more fully described below. For these operating modes, it should be understood that neithersensor 22 nor 24 needs to be moved to change the designated monitoring direction, and the designated monitoring direction need not be coincident with axis AZ. - In one embodiment,
sensors sources -
Sensors subsystem 30 to process signals received therefrom. For the convenience of description,sensors sensors processing subsystem 30 are designated xL(t) and xR(t) for the respective channels L andR. Processing subsystem 30 is operable to provide an output signal that suppresses interference fromsources acoustic source 12 positioned along axis AZ. This output signal is provided tooutput device 90 for presentation to a user in the form of an audible or visual signal which can be further processed. - Referring additionally to
FIG. 2 , a diagram is provided that depicts other details ofsystem 10.Processing subsystem 30 includes signal conditioner/filters 32 a and 32 b to filter and condition input signals xL(t) and xR(t) fromsensors filter converters Processing subsystem 30 also includesdigital circuitry 40 comprisingprocessor 42 andmemory 50. Discrete signals xL(z) and xR(z) are stored insample buffer 52 ofmemory 50 in a First-In-First-Out (FIFO) fashion. -
Processor 42 can be a software or firmware programmable device, a state logic machine, or a combination of both programmable and dedicated hardware. Furthermore,processor 42 can be comprised of one or more components and can include one or more Central Processing Units (CPUs). In one embodiment,processor 42 is in the form of a digitally programmable, highly integrated semiconductor chip particularly suited for signal processing. In other embodiments,processor 42 may be of a general purpose type or other arrangement as would occur to those skilled in the art. - Likewise,
memory 50 can be variously configured as would occur to those skilled in the art.Memory 50 can include one or more types of solid-state electronic memory, magnetic memory, or optical memory of the volatile and/or nonvolatile variety. Furthermore, memory can be integral with one or more other components ofprocessing subsystem 30 and/or comprised of one or more distinct components. -
Processing subsystem 30 can include any oscillators, control clocks, interfaces, signal conditioners, additional filters, limiters, converters, power supplies, communication ports, or other types of components as would occur to those skilled in the art to implement the present invention. In one embodiment,subsystem 30 is provided in the form of a single microelectronic device. - Referring also to the flow chart of
FIG. 3 , routine 140 is illustrated.Digital circuitry 40 is configured to perform routine 140.Processor 42 executes logic to perform at least some the operations ofroutine 140. By way of nonlimiting example, this logic can be in the form of software programming instructions, hardware, firmware, or a combination of these. The logic can be partially or completely stored onmemory 50 and/or provided with one or more other components or devices. By way of nonlimiting example, such logic can be provided toprocessing subsystem 30 in the form of signals that are carried by a transmission medium such as a computer network or other wired and/or wireless communication network. - In
stage 142, routine 140 begins with initiation of the A/D sampling and storage of the resulting discrete input samples xL(z) and xR(z) inbuffer 52 as previously described. Sampling is performed in parallel with other stages of routine 140 as will become apparent from the following description.Routine 140 proceeds fromstage 142 to conditional 144. Conditional 144 tests whether routine 140 is to continue. If not, routine 140 halts. Otherwise, routine 140 continues withstage 146. Conditional 144 can correspond to an operator switch, control signal, or power control associated with system 10 (not shown). - In
stage 146, a fast discrete fourier transform (FFT) algorithm is executed on a sequence of samples xL(Z) and xR(z) and stored inbuffer 54 for each channel L and R to provide corresponding frequency domain signals XL(k) and XR(k); where k is an index to the discrete frequencies of the FFTs (alternatively referred to as “frequency bins” herein). The set of samples xL(z) and xR(z) upon which an FFT is performed can be described in terms of a time duration of the sample data. Typically, for a given sampling rate fS, each FFT is based on more than 100 samples. Furthermore, forstage 146, FFT calculations include application of a windowing technique to the sample data. One embodiment utilizes a Hamming window. In other embodiments, data windowing can be absent or a different type utilized, the FFT can be based on a different sampling approach, and/or a different transform can be employed as would occur to those skilled in the art. After the transformation, the resulting spectra XL(k) and XR(k) are stored inFFT buffer 54 ofmemory 50. These spectra are generally complex-valued. - It has been found that reception of acoustic excitation emanating from a desired direction can be improved by weighting and summing the input signals in a manner arranged to minimize the variance (or equivalently, the energy) of the resulting output signal while under the constraint that signals from the desired direction are output with a predetermined gain. The following relationship (1) expresses this linear combination of the frequency domain input signals:
Y(k)=W L*(k)X L(k)+W R*(k)X R(k)=W H(k)X(k); (1)
where:
Y(k) is the output signal in frequency domain form, WL(k) and WR(k) are complex valued multipliers (weights) for each frequency k corresponding to channels L and R, the superscript “*” denotes the complex conjugate operation, and the superscript “H” denotes taking the Hermitian of a vector. For this approach, it is desired to determine an “optimal” set of weights WL(k) and WR(k) to minimize variance of Y(k). Minimizing the variance generally causes cancellation of sources not aligned with the desired direction. For the mode of operation where the desired direction is along axis AZ, frequency components which do not originate from directly ahead of the array are attenuated because they are not consistent in phase across the left and right channels L, R, and therefore have a larger variance than a source directly ahead. Minimizing the variance in this case is equivalent to minimizing the output power of off-axis sources, as related by the optimization goal of relationship (2) that follows:
where Y(k) is the output signal described in connection with relationship (1). In one form, the constraint requires that “on axis” acoustic signals from sources along the axis AZ be passed with unity gain as provided in relationship (3) that follows:
e H W(k)=1 (3)
Here e is a two element vector which corresponds to the desired direction. When this direction is coincident with axis AZ,sensors source 12 of the illustrated embodiment, the vector e is real-valued with equal weighted elements—for instance eH=[0.5 0.5]. In contrast, if the selected acoustic source is not on axis AZ, thensensors - In an additional or alternative mode of operation, the elements of vector e can be selected to monitor along a desired direction that is not coincident with axis AZ. For such operating modes, vector e becomes complex-valued to represent the appropriate time/phase delays between
sensors sensor Procedure 520 described in connection with the flowchart ofFIG. 10 hereinafter provides an example of a localization/tracking routine that can be used in conjunction with routine 140 to steer vector e. - For inputs XL(k) and XR(k) that generally correspond to stationary random processes (which is typical of speech signals over small periods of time), the following weight vector W(k) relationship (4) can be determined from relationships (2) and (3):
where e is the vector associated with the desired reception direction, R(k) is the correlation matrix for the kth frequency, W(k) is the optimal weight vector for the kth frequency and the superscript “−1” denotes the matrix inverse. The derivation of this relationship is explained in connection with a general model of the present invention applicable to embodiments with more than twosensors array 20. - The correlation matrix R(k) can be estimated from spectral data obtained via a number “F” of fast discrete Fourier transforms (FFTs) calculated over a relevant time interval. For the two channel L, R embodiment, the correlation matrix for the kth frequency, R(k), is expressed by the following relationship (5):
where Xl is the FFT in the frequency buffer for the left channel L and Xr is the FFT in the frequency buffer for right channel R obtained from previously stored FFTs that were calculated from an earlier execution ofstage 146; “n” is an index to the number “F” of FFTs used for the calculation; and “M” is a regularization parameter. The terms Xll(k), Xlr(k), Xrl(k), and Xrr(k) represent the weighted sums for purposes of compact expression. It should be appreciated that the elements of the R(k) matrix are nonlinear, and therefore Y(k) is a nonlinear function of the inputs. - Accordingly, in
stage 148 spectra Xl(k) and Xr(k) previously stored inbuffer 54 are read frommemory 50 in a First-In-First-Out (FIFO) sequence.Routine 140 then proceeds to stage 150. Instage 150, multiplier weights WL(k), WR(k) are applied to Xl(k) and Xr(k), respectively, in accordance with the relationship (1) for each frequency k to provide the output spectra Y(k).Routine 140 continues withstage 152 which performs an Inverse Fast Fourier Transform (IFFT) to change the Y(k) FFT determined instage 150 into a discrete time domain form designated y(z). Next, instage 154, a Digital-to-Analog (D/A) conversion is performed with D/A converter 84 (FIG. 2 ) to provide an analog output signal y(t). It should be understood that correspondence between Y(k) FFTs and output sample y(z) can vary. In one embodiment, there is one Y(k) FFT output for every y(z), providing a one-to-one correspondence. In another embodiment, there may be one Y(k) FFT for every 16 output samples y(z) desired, in which case the extra samples can be obtained from available Y(k) FFTs. In still other embodiments, a different correspondence may be established. - After conversion to the continuous time domain form, signal y(t) is input to signal conditioner/
filter 86. Conditioner/filter 86 provides the conditioned signal tooutput device 90. As illustrated inFIG. 2 ,output device 90 includes anamplifier 92 andaudio output device 94.Device 94 may be a loudspeaker, hearing aid receiver output, or other device as would occur to those skilled in the art. It should be appreciated thatsystem 10 processes a binaural input to produce an monaural output. In some embodiments, this output could be further processed to provide multiple outputs. In one hearing aid application example, two outputs are provided that deliver generally the same sound to each ear of a user. In another hearing aid application, the sound provided to each ear selectively differs in terms of intensity and/or timing to account for differences in the orientation of the sound source to eachsensor - After
stage 154, routine 140 continues with conditional 156. In many applications it may not be desirable to recalculate the elements of weight vector W(k) for every Y(k). Accordingly, conditional 156 tests whether a desired time interval has passed since the last calculation of vector W(k). If this time period has not lapsed, then control flows to stage 158 to shiftbuffers stage 158,processing loop 160 closes, returning to conditional 144. Provided conditional 144 remains true,stage 146 is repeated for the next group of samples of xL(z) and xR(Z) to determine the next pair of XL(k) and XR(k) FFTs for storage inbuffer 54. Also, with each execution ofprocessing loop 160, stages 148, 150, 152, 154 are repeated to process previously stored Xl(k) and Xr(k) FFTs to determine the next Y(k) FFT and correspondingly generate a continuous y(t). In this manner buffers 52, 54 are periodically shifted instage 158 with each repetition ofloop 160 until either routine 140 halts as tested by conditional 144 or the time period of conditional 156 has lapsed. - If the test of conditional 156 is true, then routine 140 proceeds from the affirmative branch of conditional 156 to calculate the correlation matrix R(k) in accordance with relationship (5) in
stage 162. From this new correlation matrix R(k), an updated vector W(k) is determined in accordance with relationship (4) instage 164. Fromstage 164,update loop 170 continues withstage 158 previously described, andprocessing loop 160 is re-entered untilroutine 140 halts per conditional 144 or the time for another recalculation of vector W(k) arrives. Notably, the time period tested in conditional 156 may be measured in terms of the number oftimes loop 160 is repeated, the number of FFTs or samples generated between updates, and the like. Alternatively, the period between updates can be dynamically adjusted based on feedback from an operator or monitoring device (not shown). - When routine 140 initially starts, earlier stored data is not generally available. Accordingly, appropriate seed values may be stored in
buffers array 20 and routine 140 can be adjusted accordingly. For this more general form, the output can be expressed by relationship (6) as follows:
Y(k)=W H(k)X(k) (6)
where the X(k) is a vector with an entry for each of “C” number of input channels and the weight vector W(k) is of like dimension. Equation (6) is the same at equation (1) but the dimension of each vector is C instead of 2. The output power can be expressed by relationship (7) as follows:
E[Y(k)2 ]=E[W(k)HX(k)XH(k)W(k)]=W(k)H R(k)W(k) (7)
where the correlation matrix R(k) is square with “C×C” dimensions. The vector e is the steering vector describing the weights and delays associated with a desired monitoring direction and is of the form provided by relationships (8) and (9) that follow:
where C is the number of array elements, c is the speed of sound in meters per second, and θ is the desired “look direction.” Thus, vector e may be varied with frequency to change the desired monitoring direction or look-direction and correspondingly steer the array. With the same constraint regarding vector e as described by relationship (3), the problem can be summarized by relationship (10) as follows:
This problem can be solved using the method of Lagrange multipliers generally characterized by relationship (11) as follows:
where the cost function is the output power, and the constraint is as listed above for vector e. A general vector solution begins with the Lagrange multiplier function H(W) of relationship (12):
where the factor of one half (½) is introduced to simplify later math. Taking the gradient of H(W) with respect to W(k), and setting this result equal to zero, relationship (13) results as follows:
∇w H(W)=R(k)W(k)+eλ=0 (13)
Also, relationship (14) follows:
W(k)=−R(k)−1 eλ (14)
Using this result in the constraint equation relationships (15) and (16) that follow:
e H └−R(k)−1 eλ┘=1 (15)
λ=−[e H R(k)−1 e] −1 (16)
and using relationship (14), the optimal weights are as set forth in relationship (17):
W opt =R(k)−4 e[e H R(k)−1 e] −1 (17)
Because the bracketed term is a scalar, relationship (4) has this term in the denominator, and thus is equivalent. - Returning to the two variable case for the sake of clarity, relationship (5) may be expressed more compactly by absorbing the weighted sums into the terms Xll, Xlr, Xrl and Xrr, and then renaming them as components of the correlation matrix R(k) per relationship (18):
Its inverse may be expressed in relationship (19) as:
where det( ) is the determinant operator. If the desired monitoring direction is perpendicular to the sensor array, e=[0.5 0.5]T, the numerator of relationship (4) may then be expressed by relationship (20) as:
Using the previous result, the denominator is expressed by relationship (21) as:
Canceling out the common factor of the determinant, the simplified relationship (22) is completed as:
It can also be expressed in terms of averages of the sums of correlations between the two channels in relationship (23) as:
where wl(k) and wr(k) are the desired weights for the left and right channels, 5 respectively, for the kth frequency, and the components of the correlation matrix are now expressed by relationships (24) as:
just as in relationship (5). Thus, after computing the averaged sums (which may be kept as running averages), computational load can be reduced for this two channel embodiment. - In a further variation of routine 140, a modified approach can be utilized in applications where gain differences between sensors of
array 20 are negligible. For this approach, an additional constraint is utilized. For a two-sensor arrangement with a fixed on-axis steering direction and negligible inter-sensor gain differences, the desired weights satisfy relationship (25) as follows:
The variance minimization goal and unity gain constraint for this alternative approach correspond to the following relationships (26) and (27), respectively:
By inspection, when eH=[1 1], relationship (27) reduces to relationship (28) as follows:
Im[w1]=−Im[w2] (28)
Solving for desired weights subject to the constraint in relationship (27) and using relationship (28) results in the following relationship (29): - The weights determined in accordance with relationship (29) can be used in place of those determined with relationships (22), (23), and (24); where R11, R12, R21, R22, are the same as those described in connection with relationship (18). Under appropriate conditions, this substitution typically provides comparable results with more efficient computation. When relationship (29) is utilized, it is generally desirable for the target speech or other acoustic signal to originate from the on-axis direction and for the sensors to be matched to one another or to otherwise compensate for inter-sensor differences in gain. Alternatively, localization information about sources of interest in each frequency band can be utilized to steer
sensor array 20 in conjunction with the relationship (29) approach. This information can be provided in accordance withprocedure 520 more fully described hereinafter in connection with the flowchart ofFIG. 10 . - Referring to relationship (5), regularization factor M typically is slightly greater than 1.00 to limit the magnitude of the weights in the event that the correlation matrix R(k) is, or is close to being, singular, and therefore noninvertable. This occurs, for example, when time-domain input signals are exactly the same for F consecutive FFT calculations. It has been found that this form of regularization also can improve the perceived sound quality by reducing or eliminating processing artifacts common to time-domain beamformers.
- In one embodiment, regularization factor M is a constant. In other embodiments, regularization factor M can be used to adjust or otherwise control the array beamwidth, or the angular range at which a sound of a particular frequency can impinge on the array relative to axis AZ and be processed by routine 140 without significant attenuation. This beamwidth is typically larger at lower frequencies than higher frequencies, and can be expressed by the following relationship (30):
r=1−M, where M is the regularization factor, as in relationship (5), c represents the speed of sound in meters per second (m/s),f represents frequency in Hertz (Hz), D is the distance between microphones in meters (m). For relationship (30), Beamwidth—3 dB defines a beamwidth that attenuates the signal of interest by a relative amount less than or equal to three decibels (dB). It should be understood that a different attenuation threshold can be selected to define beamwidth in other embodiments of the present invention.FIG. 9 provides a graph of four lines of different patterns to represent constant values 1.001, 1.005, 1.01, and 1.03, of regularization factor M, respectively, in terms of beamwidth versus frequency. - Per relationship (30), as frequency increases, beamwidth decreases; and as regularization factor M increases, the beamwidth increases. Accordingly, in one alternative embodiment of routine 140, regularization factor M is increased as a function of frequency to provide a more uniform beamwidth across a desired range of frequencies. In another embodiment of routine 140, M is alternatively or additionally varied as a function of time. For example, if little interference is present in the input signals in certain frequency bands, the regularization factor M can be increased in those bands. It has been found that beamwidth increases in frequency bands with low or no inference commonly provide a better subjective sound quality by limiting the magnitude of the weights used in relationships (22), (23), and/or (29). In a further variation, this improvement can be complemented by decreasing regularization factor M for frequency bands that contain interference above a selected threshold. It has been found that such decreases commonly provide more accurate filtering, and better cancellation of interference. In still another embodiment, regularization factor M varies in accordance with an adaptive function based on frequency-band-specific interference. In yet further embodiments, regularization factor M varies in accordance with one or more other relationships as would occur to those skilled in the art.
- Referring to
FIG. 4 , one application of the various embodiments of the present invention is depicted as hearingaid system 210; where like reference numerals refer to like features. In one embodiment,system 210 includes eyeglasses G andacoustic sensors Acoustic sensors processor 30.Processor 30 is operatively coupled tooutput device 190.Output device 190 is in the form of a hearing aid earphone and is positioned in ear E of the user to provide a corresponding audio signal. Forsystem 210,processor 30 is configured to perform routine 140 or its variants with the output signal y(t) being provided tooutput device 190 instead ofoutput device 90 ofFIG. 2 . As previously discussed, anadditional output device 190 can be coupled toprocessor 30 to provide sound to another ear (not shown). This arrangement defines axis AZ to be perpendicular to the view plane ofFIG. 4 as designated by the like labeled cross-hairs located generally midway betweensensors - In operation, the user wearing eyeglasses G can selectively receive an acoustic signal by aligning the corresponding source with a designated direction, such as axis AZ. As a result, sources from other directions are attenuated. Moreover, the wearer may select a different signal by realigning axis AZ with another desired sound source and correspondingly suppress a different set of off-axis sources. Alternatively or additionally,
system 210 can be configured to operate with a reception direction that is not coincident with axis AZ. -
Processor 30 andoutput device 190 may be separate units (as depicted) or included in a common unit worn in the ear. The coupling betweenprocessor 30 andoutput device 190 may be an electrical cable or a wireless transmission. In one alternative embodiment,sensors processor 30 are remotely located relative to each other and are configured to broadcast to one ormore output devices 190 situated in the ear E via a radio frequency transmission. - In a further hearing aid embodiment,
sensors - Another hearing aid system embodiment is based on a cochlear implant. A cochlear implant is typically disposed in a middle ear passage of a user and is configured to provide electrical stimulation signals along the middle ear in a standard manner. The implant can include some or all of
processing subsystem 30 to operate in accordance with the teachings of the present invention. Alternatively or additionally, one or more external modules include some or all ofsubsystem 30. Typically a sensor array associated with a hearing aid system based on a cochlear implant is worn externally, being arranged to communicate with the implant through wires, cables, and/or by using a wireless technique. - Besides various forms of hearing aids, the present invention can be applied in other configurations. For instance,
FIG. 5 shows avoice input device 310 employing the present invention as a front end speech enhancement device for a voice recognition routine for personal computer C; where like reference numerals refer to like features.Device 310 includesacoustic sensors Sensors processor 330 within computer C. Processor .330 provides an output signal for internal use or responsive reply viaspeakers visual display 396; and is arranged to process vocal inputs fromsensors device 310. In another mode of operation,device 310 changes its monitoring direction based on feedback from an operator and/or automatically selects a monitoring direction based on the location of the most intense sound source over a selected period of time. Alternatively or additionally, the source localization/tracking ability provided byprocedure 520 as illustrated in the flowchart ofFIG. 10 can be utilized. In still another voice input application, the directionally selective speech processing features of the present invention are utilized to enhance performance of a hands-free telephone, audio surveillance device, or other audio system. - Under certain circumstances, the directional orientation of a sensor array relative to the target acoustic source changes. Without accounting for such changes, attenuation of the target signal can result. This situation can arise, for example, when a binaural hearing aid wearer turns his or her head so that he or she is not aligned properly with the target source, and the hearing aid does not otherwise account for this misalignment. It has been found that attenuation due to misalignment can be reduced by localizing and/or tracking one or more acoustic sources of interests. The flowchart of
FIG. 10 illustratesprocedure 520 to track and/or localize a desired acoustic source relative to a reference.Procedure 520 can be utilized for a hearing aid or in other applications such as a voice input device, a hands-free telephone, audio surveillance equipment, and the like—either in conjunction with or independent of previously described embodiments.Procedure 520 is described as follows in terms of an implementation withsystem 10 ofFIG. 1 . For this embodiment,processing system 30 can include logic to execute one or more stages and/or conditionals ofprocedure 520 as appropriate. In other embodiments, a different arrangement can be used to implementprocedure 520 as would occur to one skilled in the art. -
Procedure 520 starts with AID conversion instage 522 in a manner like that described forstage 142 ofroutine 140. Fromstage 522,procedure 520 continues withstage 524 to transform the digital data obtained fromstage 522, such that “G” number of FFTs are provided each with “N” number of FFT frequency bins.Stages procedure 520 in a parallel, pipelined, sequence-specific, or different manner as would occur to one skilled in the art. With the FFTs fromstage 524, an array of localization results, P(γ), can be described in terms of relationships (31)-(35) as follows:
where the operator “INT” returns the integer part of its operand, L(g,k) and R(g,k) are the frequency-domain data from channels L and R, respectively, for the kth FFT frequency bin of the gth FFT, Mthr(k) is a threshold value for the frequency-domain data in FFT frequency bin k, the operator “ROUND” returns the nearest integer degree of its operand, c is the speed of sound in meters per second, fs is the sampling rate in Hertz, and D is the distance (in meters) between the two sensors ofarray 20. For these relationships, array P(γ) is defined with 181 azimuth location elements, which correspond to directions −90° to +90° in 1° increments. In other embodiments, a different resolution and/or location indication technique can be used. - From
stage 524,procedure 520 continues withindex initialization stage 526 in which index g to the G number of FFTs and index k to the N frequency bins of each FFT are set to one and zero, (g=1, k=0), respectively. Fromstage 526,procedure 520 continues by entering frequencybin processing loop 530 andFFT processing loop 540. For this example,loop 530 is nested withinloop 540.Loops stage 532. - For an off-axis acoustic source, the corresponding signal travels different distances to reach each of the
sensors array 20. Generally, these different distances cause a phase difference between channels L and R at some frequency. Instage 532, routine 520 determines the difference in phase between channels L and R for the current frequency bin k of the FFT g, converts the phase difference to a difference in distance, and determines the ratio x(g,k) of this distance difference to the sensor spacing D in accordance with relationship (35). Ratio x(g,k) is used to find the signal angle of arrival θx, rounded to the nearest degree, in accordance with relationship (34). - Conditional 534 is next encountered to test whether the signal energy level in channels L and R have more energy than a threshold level Mthr and the value of x(g,k) was one for which a valid angle of arrival could be calculated. If both conditions are met, then in stage 535 a value of one is added to the corresponding element of P(γ), where y=θx.
Procedure 520 proceeds fromstage 535 to conditional 536. If neither condition of conditional 534 is met, then P(γ) is not modified, andprocedure 520 bypassesstage 535, continuing with conditional 536. - Conditional 536 tests if all the frequency bins have been processed, that is whether index k equals N, the total number of bins. If not (conditional 536 test is negative),
procedure 520 continues withstage 537 in which index k is incremented by one (k=k+1). Fromstage 537,loop 530 closes, returning tostage 532 to process the new g and k combination. If the conditional 536 test is affirmative, conditional 542 is next encountered, which tests if all FFTs have been processed, that is whether index g equals G number of FFTs. If not (conditional 542 is negative),procedure 520 continues withstage 544 to increment g by one (g=g+1) and to reset k to zero (k=0). Fromstage 544,loop 540 closes, returning tostage 532 to process the new g and k combination. Ifconditional test 542 is affirmative, then all N bins for each of the G number of FFTs have been processed, andloops - With the conclusion of processing by
loops loops procedure 520 continues withstage 550. - In
stage 550, the elements of array P(γ) having the greatest relative values, or “peaks,” are identified in accordance with relationship (36) as follows:
p(l)=PEAKS(P(γ), γlim, Pthr) (36)
where p(l) is direction of the lth peak in the function P(γ) for values of γ between ±γlim (a typical value for γlim is 10°, but this may vary significantly) and for which the peak values are above the threshold value Pthr. The PEAKS operation of relationship (36) can use a number of peak-finding algorithms to locate maxima of the data, including optionally smoothing the data and other operations. - From
stage 550,procedure 520 continues withstage 552 in which one or more peaks are selected. When tracking a source that was initially on-axis, the peak closest to the on-axis direction typically corresponds to the desired source. The selection of this closest peak can be performed in accordance with relationship (37) as follows:
where θtar is the direction angle of the chosen peak. Regardless of the selection criteria,procedure 520 proceeds to stage 554 to apply the selected peak or peaks.Procedure 520 continues fromstage 554 to conditional 560. Conditional 560 tests whetherprocedure 520 is to continue or not. If the conditional 560 test is true,procedure 520 loops back tostage 522. If the conditional 560 test is false,procedure 520 halts. - In an application relating to routine 140, the peak closest to axis AZ is selected, and utilized to steer
array 20 by adjusting steering vector e. In this application, vector e is modified for each frequency bin k so that it corresponds to the closest peak direction θtar. For a steering direction of θtar, the vector e can be represented by the following relationship (38), which is a simplified version of relationships (8) and (9):
where k is the FFT frequency bin number, D is the distance in meters betweensensors routine 140, the modified steering vector e of relationship (38) can be substituted into relationship (4) ofroutine 140 to extract a signal originating from direction θtar. Likewise,procedure 520 can be integrated with routine 140 to perform localization with the same FFT data. In other words, the AID conversion ofstage 142 can be used to provide digital data for subsequent processing by both routine 140 andprocedure 520. Alternatively or additionally, some or all of the FFTs obtained for routine 140 can be used to provide the G FFTs forprocedure 520. Moreover, beamwidth modifications can be combined withprocedure 520 in various applications either with or withoutroutine 140. In still other embodiments, the indexed execution ofloops routine 140. - In a further embodiment, one or more transformation techniques are utilized in addition to or as an alternative to fourier transforms in one or more forms of the invention previously described. One example is the wavelet transform, which mathematically breaks up the time-domain waveform into many simple waveforms, which may vary widely in shape. Typically wavelet basis functions are similarly shaped signals with logarithmically spaced frequencies. As frequency rises, the basis functions become shorter in time duration with the inverse of frequency. Like fourier transforms, wavelet transforms represent the processed signal with several different components that retain amplitude and phase information. Accordingly, routine 140 and/or routine 520 can be adapted to use such alternative or additional transformation techniques. In general, any signal transform components that provide amplitude and/or phase information about different parts of an input signal and have a corresponding inverse transformation can be applied in addition to or in place of FFTs.
-
Routine 140 and the variations previously described generally adapt more quickly to signal changes than conventional time-domain iterative-adaptive schemes. In certain applications where the input signal changes rapidly over a small interval of time, it may be desired to be more responsive to such changes. For these applications, the F number of FFTs associated with correlation matrix R(k) may provide a more desirable result if it is not constant for all signals (alternatively designated the correlation length F). Generally, a smaller correlation length F is best for rapidly changing input signals, while a larger correlation length F is best for slowly changing input signals. - A varying correlation length F can be implemented in a number of ways. In one example, filter weights are determined using different parts of the frequency-domain data stored in the correlation buffers. For buffer storage in the order of the time they are obtained (First-In, First-Out (FIFO) storage), the first half of the correlation buffer contains data obtained from the first half of the subject time interval and the second half of the buffer contains data from the second half of this time interval. Accordingly, the correlation matrices R1(k) and R2(k) can be determined for each buffer half according to relationships (39) and (40) as follows:
R(k) can be obtained by summing correlation matrices R1(k) and R2(k). - Using relationship (4) of
routine 140, filter coefficients (weights) can be obtained using both R1(k) and R2(k). If the weights differ significantly for some frequency band k between R1(k) and R2(k), a significant change in signal statistics may be indicated. This change can be quantified by examining the change in one weight through determining the magnitude and phase change of the weight and then using these quantities in a function to select the appropriate correlation length F. The magnitude difference is defined according to relationship (41) as follows:
ΔM(k)=||w 1.1(k)|−|w1.2(k)|| (41)
where w1.1(k) and w1.2(k) are the weights calculated for the left channel using R1(k) and R2(k), respectively. The angle difference is defined according to relationship (42) as follows:
where the factor of ±2π is introduced to provide the actual phase difference in the case of a ±2π jump in the phase of one of the angles. - The correlation length F for some frequency bin k is now denoted as F(k). An example function is given by the following relationship (43):
F(k)=max(b(k)·ΔA(k)+d(k)·ΔM(k)+c max(k), C min(k)) (43)
where Cmin(k) represents the minimum correlation length, cmax(k) represents the maximum correlation length and b(k) and d(k) are negative constants, all for the kth frequency band. Thus, as ΔA(k) and ΔM(k) increase, indicating a change in the data, the output of the function decreases. With proper choice of b(k) and d(k), F(k) is limited between cmin(k) and cmin(k), so that the correlation length can vary only within a predetermined range. It should also be understood that F(k) may take different forms, such as a nonlinear function or a function of other measures of the input signals. - Values for function F(k) are obtained for each frequency bin k. It is possible that a small number of correlation lengths may be used, so in each frequency bin k the correlation length that is closest to F1(k) is used to form R(k). This closest value is found using relationship (44) as follows:
where imin, is the index for the minimized function F(k) and c(i) is the set of possible correlation length values ranging from cmin to cmax. - The adaptive correlation length process described in connection with relationships (39)-(44) can be incorporated into the
correlation matrix stage 162 andweight determination stage 164 for use in a hearing aid, such as that described in connection withFIG. 4 , or other applications like surveillance equipment, voice recognition systems, and hands-free telephones, just to name a few. Logic ofprocessing subsystem 30 can be adjusted as appropriate to provide for this incorporation. Optionally, the adaptive correlation length process can be utilized with the relationship (29) approach to weight computation, the dynamic beamwidth regularization factor variation described in connection with relationship (30) andFIG. 9 , the localization/tracking procedure 520, alternative transformation embodiments, and/or such different embodiments or variations of routine 140 as would occur to one skilled in the art. The application of adaptive correlation length can be operator selected and/or automatically applied based on one or more measured parameters as would occur to those skilled in the art. - Many other further embodiments of the present invention are envisioned. One further embodiment includes: detecting acoustic excitation with a number of acoustic sensors that provide a number of sensor signals; establishing a set of frequency components for each of the sensor signals; and determining an output signal representative of the acoustic excitation from a designated direction. This determination includes weighting the set of frequency components for each of the sensor signals to reduce variance of the output signal and provide a predefined gain of the acoustic excitation from the designated direction.
- In another embodiment, a hearing aid includes a number of acoustic sensors in the presence of multiple acoustic sources that provide a corresponding number of sensor signals. A selected one of the acoustic sources is monitored. An output signal representative of the selected one of the acoustic sources is generated. This output signal is a weighted combination of the sensor signals that is calculated to minimize variance of the output signal.
- A still further embodiment includes: operating a voice input device including a number of acoustic sensors that provide a corresponding number of sensor signals; determining a set of frequency components for each of the sensor signals; and generating an output signal representative of acoustic excitation from a designated direction. This output signal is a weighted combination of the set of frequency components for each of the sensor signals calculated to minimize variance of the output signal.
- Yet a further embodiment includes an acoustic sensor array operable to detect acoustic excitation that includes two or more acoustic sensors each operable to provide a respective one of a number of sensor signals. Also included is a processor to determine a set of frequency components for each of the sensor signals and generate an output signal representative of the acoustic excitation from a designated direction. This output signal is calculated from a weighted combination of the set of frequency components for each of the sensor signals to reduce variance of the output signal subject to a gain constraint for the acoustic excitation from the designated direction.
- A further embodiment includes: detecting acoustic excitation with a number of acoustic sensors that provide a corresponding number of signals; establishing a number of signal transform components for each of these signals; and determining an output signal representative of acoustic excitation from a designated direction. The signal transform components can be of the frequency domain type. Alternatively or additionally, a determination of the output signal can include weighting the components to reduce variance of the output signal and provide a predefined gain of the acoustic excitation from the designated direction.
- In yet another embodiment, a hearing aid is operated that includes a number of acoustic sensors. These sensors provide a corresponding number of sensor signals. A direction is selected to monitor for acoustic excitation with the hearing aid. A set of signal transform components for each of the sensor signals is determined and a number of weight values are calculated as a function of a correlation of these components, an adjustment factor, and the selected direction. The signal transform components are weighted with the weight values to provide an output signal representative of the acoustic excitation emanating from the direction. The adjustment factor can be directed to correlation length or a beamwidth control parameter just to name a few examples.
- For a further embodiment, a hearing aid is operated that includes a number of acoustic sensors to provide a corresponding number of sensor signals. A set of signal transform components are provided for each of the sensor signals and a number of weight values are calculated as a function of a correlation of the transform components for each of a number of different frequencies. This calculation includes applying a first beamwidth control value for a first one of the frequencies and a second beamwidth control value for a second one of the frequencies that is different than the first value. The signal transform components are weighted with the weight values to provide an output signal.
- For another embodiment, acoustic sensors of the hearing aid provide corresponding signals that are represented by a plurality of signal transform components. A first set of weight values are calculated as a function of a first correlation of a first number of these components that correspond to a first correlation length. A second set of weight values are calculated as a function of a second correlation of a second number of these components that correspond to a second correlation length different than the first correlation length. An output signal is generated as a function of the first and second weight values.
- In another embodiment, acoustic excitation is detected with a number of sensors that provide a corresponding number of sensor signals. A set of signal transform components is determined for each of these signals. At least one acoustic source is localized as a function of the transform components. In one form of this embodiment, the location of one or more acoustic sources can be tracked relative to a reference. Alternatively or additionally, an output signal can be provided as a function of the location of the acoustic source determined by localization and/or tracking, and a correlation of the transform components.
- It is contemplated that various signal flow operators, converters, functional blocks, generators, units, stages, processes, and techniques may be altered, rearranged, substituted, deleted, duplicated, combined or added as would occur to those skilled in the art without departing from the spirit of the present inventions. It should be understood that the operations of any routine, procedure, or variant thereof can be executed in parallel, in a pipeline manner, in a specific sequence, as a combination of these appropriate to the interdependence of such operations on one another, or as would otherwise occur to those skilled in the art. By way of nonlimiting example, A/D conversion, D/A conversion, FFT generation, and FFT inversion can typically be performed as other operations are being executed. These other operations could be directed to processing of previously stored AID or signal transform components, such as
stages - The following experimental results provide nonlimiting examples, and should not be construed to restrict the scope of the present invention.
-
FIG. 6 illustrates the experimental set-up for testing the present invention. The algorithm has been tested with real recorded speech signals, played through loudspeakers at different spatial locations relative to the receiving microphones in an anechoic chamber. A pair ofmicrophones 422, 424 (Sennheiser MKE 2-60) with an inter-microphone distance D of 15 cm, were situated in a listening room to serve assensors microphones FIG. 2 ). Several loudspeakers were used to broadcast words or sentences that interfere with the listening of target speech from different azimuths. -
Microphones preamp 432 was provided to a dualchannel volume control 434 provided in the form of an audio preamplifier (Adcom GTP-5511). The output ofvolume control 434 was fed into A/D converters of a Digital Signal Processor (DSP)development board 440 provided by Texas Instruments (model number TI-C6201 DSP Evaluation Module (EVM)).Development board 440 includes a fixed-point DSP chip (model number TMS320C62) running at a clock speed of 133 MHz with a peak throughput of 1064 MIPS (millions of instructions per second). This DSP executed software configured to implement routine 140 in real-time. The sampling frequency for these experiments was about 8 kHz with 16-bit A/D and D/A conversion. The FFT length was 256 samples, with an FFT calculated every 16 samples. The computation leading to the characterization and extraction of the desired signal was found to introduce a delay in a range of about 10-20 milliseconds between the input and output. -
FIGS. 7 and 8 each depict traces of three acoustic signals of approximately the same energy. InFIG. 7 , the target signal trace is shown between two interfering signals traces broadcast fromazimuths 22° and −65°, respectively. These azimuths are depicted inFIG. 1 . The target sound is a prerecorded voice from a female (second trace), and is emitted by the loudspeaker located near 0°. One interfering sound is provided by a female talker (top trace ofFIG. 7 ) and the other interfering sound is provided by a male talker (bottom trace ofFIG. 7 ). The phrase repeated by the corresponding talker is reproduced above the respective trace. - Referring to
FIG. 8 , as revealed by the top trace, when the target speech sound is emitted in the presence of two interfering sources, its waveform (and power spectrum) is contaminated. This contaminated sound was difficult to understand for most listeners, especially those with hearing impairment.Routine 140, as embodied inboard 440, processed this contaminated signal with high fidelity and extracted the target signal by markedly suppressing the interfering sounds. Accordingly, intelligibility of the target signal was restored as illustrated by the second trace. The intelligibility was significantly improved and the extracted signal resembled the original target signal reproduced for comparative purposes as the bottom trace ofFIG. 8 . - These experiments demonstrate marked suppression of interfering sounds. The use of the regularization parameter (valued at approximately 1.03) effectively limited the magnitude of the calculated weights and results in an output with much less audible distortion when the target source is slightly off-axis, as would occur when the hearing aid wearer's head is slightly misaligned to the target talker. Miniaturization of this technology to a size suitable for hearing aids and other applications can be provided using techniques known to those skilled in the art.
-
FIGS. 11 and 12 are computer generated image graphs of simulated results forprocedure 520. These graphs plot localization results of azimuth in degrees versus time in seconds. The localization results are plotted as shading, where the darker the shading, the stronger the localization result at that angle and time. Such simulations are accepted by those skilled in the art to indicate efficacy of this type of procedure. -
FIG. 11 illustrates the localization results when the target acoustic source is generally stationary with a direction of about 10° off-axis. The actual direction of the target is indicated by a solid black line.FIG. 12 illustrates the localization results for a target with a direction that is changing sinusoidally between +10° and −10°, as might be the case for a hearing aid wearer shaking his or her head. The actual location of the source is again indicated by a solid black line. The localization technique ofprocedure 520 accurately indicates the location of the target source in both cases because the darker shading matches closely to the actual location lines. Because the target source is not always producing a signal free of interference overlap, localization results may be strong only at certain times. InFIG. 12 , these stronger intervals can be noted at about 0.2, 0.7, 0.9, 1.25, 1.7, and 2.0 seconds. It should be understood that the target location can be readily estimated between such times. - Experiments described herein are simply for the purpose of demonstrating operation of one form of a processing system of the present invention. The equipment, the speech materials, the talker configurations, and/or the parameters can be varied as would occur to those skilled in the art.
- Any theory, mechanism of operation, proof, or finding stated herein is meant to further enhance understanding of the present invention and is not intended to make the present invention in any way dependent upon such theory, mechanism of operation, proof, or finding. While the invention has been illustrated and described in detail in the drawings and foregoing description, the same is to be considered as illustrative and not restrictive in character, it being understood that only the selected embodiments have been shown and described and that all changes, modifications and equivalents that come within the spirit of the invention as defined herein or by the following claims are desired to be protected.
Claims (58)
1. A method, comprising:
detecting acoustic excitation with a number of acoustic sensors, the acoustic sensors providing a corresponding number of sensor signals;
establishing a number of frequency domain components for each of the sensor signals; and
determining an output signal representative of the acoustic excitation from a designated direction, said determining including weighting the components for each of the sensor signals to reduce variance of the output signal and provide a predefined gain of the acoustic excitation from the designated direction.
2. The method of claim 1 , wherein said determining includes minimizing the variance of the output signal and the predefined gain is approximately unity.
3. The method of claim 1 , further comprising changing the designated direction without moving any of the acoustic sensors and repeating said establishing and said determining after said changing.
4. The method of claim 1 , further comprising changing from the designated direction by moving one or more of the acoustic sensors and repeating said establishing and said determining after said changing.
5. The method of claim 1 , wherein said components correspond to fourier transforms and said weighting includes calculating a number of weights to minimize the variance of the output signal subject to a constraint that the predefined gain be generally maintained at unity, the weights being determined as a function of a frequency domain correlation matrix and a vector corresponding to the designated direction.
6. The method of claim 5 , further comprising recalculating the weights from time to time and repeating said establishing and said determining on an established basis.
7. The method of claim 1 , further comprising calculating said weights subject to a constraint of an insubstantial level of gain difference between the acoustic sensors.
8. The method of claim 1 , further comprising adjusting a correlation factor to control beamwidth as a function of frequency.
9. The method of claim 1 , further comprising calculating a number of correlation matrices and adaptively changing correlation length for one or more of the correlation matrices relative to at least one other of the correlation matrices.
10. The method of claim 1 , further comprising tracking location of at least one acoustic signal source as a function of a phase difference between the acoustic sensors.
11. The method of any of claims 1-10, further comprising providing a hearing aid with the acoustic sensors and a processor operable to perform said establishing and said determining.
12. The method of any of claims 1-10, wherein a voice input device includes the acoustic sensors and a processor operable to perform said establishing and said determining.
13. A method, comprising:
operating a hearing aid including a number of acoustic sensors in the presence of multiple acoustic sources, the acoustic sensors providing a corresponding number of sensor signals;
monitoring a selected one of the acoustic sources;
determining a set of frequency components for each of the sensor signals; and
generating an output signal representative of the selected one of the acoustic sources, the output signal being a weighted combination of the set of frequency components for each of the sensor signals calculated to minimize variance of the output signal.
14. The method of claim 13 , further comprising processing the output signal to provide at least one acoustic output to a user of the hearing aid.
15. A method, comprising:
operating a voice input device including a number of acoustic sensors, the acoustic sensors providing a corresponding number of sensor signals;
determining a set of frequency components for each of the sensor signals; and
generating an output signal representative of acoustic excitation from a designated direction, the output signal being a weighted combination of the set of frequency components for each of the sensor signals calculated to minimize variance of the output signal.
16. The method of claim 15 , wherein the voice input device is included in a voice recognition system for a computer.
17. The method of any of claims 13-16, wherein said generating includes calculating a number of weights as a function of a frequency domain correlation matrix and a vector corresponding to the designated direction.
18. The method of claim 17 , further comprising recalculating the weights from time to time.
19. The method of claim 17 , further comprising determining the weighted combination of the sensor signals as a function of a gain constraint associated with the designated direction.
20. The method of claim 17 , further comprising adjusting a correlation factor to control beamwidth as a function of frequency.
21. The method of claim 17 , further comprising adaptively changing correlation length.
22. A method, comprising:
operating a hearing aid including a number of acoustic sensors, the acoustic sensors providing a corresponding number of sensor signals;
selecting a direction to monitor for acoustic excitation with the hearing aid;
determining a set of signal transform components for each of the sensor signals;
calculating a number of weight values as a function of a correlation of the signal transform components, an adjustment factor, and the direction; and
weighting the signal transform components with the weight values to provide an output signal representative of the acoustic excitation emanating from the direction.
23. The method of claim 22 , wherein the transform components correspond to different frequencies and the adjustment factor has a first value for a first one of the frequencies and second value different than the first value for a second one of the frequencies to control beamwidth.
24. The method of claim 22 , wherein the adjustment factor corresponds to correlation length and further comprising determining a number of different correlations with correlation length adaptively changed in accordance with different values for the adjustment factor.
25. The method of claim 22 , further comprising:
determining a level of interference; and
adjusting the beamwidth of the hearing aid in response to the level of interference with the adjustment factor.
26. The method of claim 22 , further comprising:
determining a rate of change of at least one frequency of at least one of the sensor signals with respect to time; and
adjusting the correlation length in response to the rate of change with the adjustment factor.
27. A method, comprising:
operating a hearing aid including a number of acoustic sensors, the acoustic sensors providing a corresponding number of sensor signals;
providing a set of signal transform components for each of the sensor signals;
calculating a number of weight values as a function of a correlation of the transform components for each of a number different frequencies, said calculating including applying a first beamwidth control value for a first one of the frequencies and a second beamwidth control value for a second one of the frequencies different than the first beamwidth control value; and
weighting the signal transform components with the weight values to provide an output signal.
28. The method of claim 27 , further comprising selecting the first beamwidth value and the second beamwidth value to provide a generally constant beamwidth of the hearing aid over a predefined frequency range.
29. The method of claim 27 , wherein the first beamwidth value and the second beamwidth value differ in accordance with a difference in an amount of interference at the first one of the frequencies relative to the second one of the frequencies.
30. A method, comprising:
operating a hearing aid including a number of acoustic sensors, the acoustic sensors providing a corresponding number of sensor signals;
providing a first plurality of signal transform components for the sensor signals;
calculating a first set of weight values as a function of a first correlation of the first signal transform components corresponding to a first correlation length;
providing a second plurality of signal transform components for the sensor signals;
calculating a second set of weight values as a function of a second correlation of the second signal transform components corresponding to a second correlation length different that the first correlation length; and
generating an output signal as a function of the first weight values and the second weight values.
31. The method of claim 30 , wherein the first correlation length and the second correlation length differ in accordance with a difference in rate of change of at least one frequency of at least one of the sensor signals with respect to time.
32. The method of any of claims 22-31, wherein the number of sensors is two and the hearing aid has a single, monaural output.
33. The method of any of claims 22-31, wherein said calculating is performed to minimize output variance.
34. The method of any of claims 22-31, further comprising localizing a selected acoustic source relative to a reference as a function of the transform components.
35. The method of any of claims 22-31, wherein the transform components are of a fourier type.
36. A hearing aid system operable to perform the method of any of claims 22-31.
37. A method comprising:
detecting acoustic excitation with a number of acoustic sensors, the acoustic sensors providing a corresponding number of sensor signals;
establishing a set of signal transform components for each of the sensor signals;
tracking location of a source of the acoustic excitation relative to a reference as a function of the transform components; and
providing an output signal as a function of the location and a correlation of the transform components.
38. The method of claim 37 , wherein the number of sensors is two and said tracking includes determining a phase difference between the sensor signals.
39. The method of claim 37 , wherein the reference is a designated axis and the location is provided in the form of an azimuthal direction.
40. The method of claim 37 , wherein said tracking includes generating an array with a number of elements each corresponding to a different azimuth and detecting one or more peak values among the elements of the array.
41. The method of claim 37 , further comprising adjusting a beamwidth factor relative to frequency.
42. The method of claim 37 , further comprising calculating a number of different correlation matrices and adaptively changing correlation length of one or more of the matrices relative to at least one other of the matrices.
43. The method of claim 37 , further comprising steering a direction-indicating vector corresponding to the location.
44. The method of claim 37 , wherein said providing include generating the output signal by weighting the transform components to reduce variance of the output signal and provide a predefined gain.
45. A device operable to perform the method of any of claims 37-44.
46. A hearing aid system operable to perform the method of any of claims 37-44.
47. An apparatus, comprising:
an acoustic sensor array operable to detect acoustic excitation, said acoustic sensor array including two or more acoustic sensors each operable to provide a respective one of a number of sensor signals; and
a processor operable to determine a set of frequency components for each of said sensor signals and generate an output signal representative of the acoustic excitation from a designated direction, said output signal being calculated from a weighted combination of said set of frequency components for each of said sensor signals to reduce variance of said output signal subject to a gain constraint for the acoustic excitation from said designated direction.
48. The apparatus of claim 47 , wherein said processor is operable to calculate said weighted combination to generally minimize said variance of said output signal and generally maintain said gain at unity.
49. The apparatus of claim 47 , wherein said processor is operable to determine a number of signal weights as a function of a frequency domain correlation matrix and a vector corresponding to said designated direction.
50. An apparatus, comprising:
a first acoustic sensor operable to provide a first sensor signal;
a second acoustic sensor operable to provide a second sensor signal;
a processor operable to generate an output signal representative of acoustic excitation detected with said first acoustic sensor and said second acoustic sensor from a designated direction, said processor including:
means for transforming said first sensor signal to a first number of frequency domain transform components and said second sensor signal to a second number of frequency domain transform components,
means for weighting said first transform components to provide a corresponding number of first weighted components and said second transform components to provide a corresponding number of second weighted components as a function of variance of said output signal and a gain constraint for the acoustic excitation from said designated direction,
means for combining each of said first weighted components with a corresponding one of said second weighted components to provide a frequency domain form of said output signal; and
means for providing a time domain form of said output signal from said frequency domain form.
51. The apparatus of any of claims 47-50, wherein said processor includes means for steering said designated direction.
52. The apparatus of any of claims 47-50, further comprising at least one acoustic output device responsive to said output signal.
53. The apparatus of any of claims 47-50, wherein the apparatus is arranged as a hearing aid.
54. The apparatus of any of claims 47-50, wherein the apparatus is arranged as a voice input device.
55. The apparatus of any of claims 47-50, wherein said processor is operable to localize an acoustic excitation source relative to a reference.
56. The apparatus of any of claims 47-50, wherein said processor is operable to track location of an acoustic excitation source relative to an azimuthal plane.
57. The apparatus of any of claims 47-50, wherein said processor is operable to adjust a beamwidth control parameter with frequency.
58. The apparatus of any of claims 47-50, wherein said processor is operable to calculate a number of different correlation matrices and adaptively adjust correlation length of one or more of the matrices relative to at least one other of the matrices.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/545,256 US20070030982A1 (en) | 2000-05-10 | 2006-10-10 | Interference suppression techniques |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US56843000A | 2000-05-10 | 2000-05-10 | |
US10/290,137 US7613309B2 (en) | 2000-05-10 | 2002-11-07 | Interference suppression techniques |
US11/545,256 US20070030982A1 (en) | 2000-05-10 | 2006-10-10 | Interference suppression techniques |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/290,137 Continuation US7613309B2 (en) | 2000-05-10 | 2002-11-07 | Interference suppression techniques |
Publications (1)
Publication Number | Publication Date |
---|---|
US20070030982A1 true US20070030982A1 (en) | 2007-02-08 |
Family
ID=24271254
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/290,137 Expired - Fee Related US7613309B2 (en) | 2000-05-10 | 2002-11-07 | Interference suppression techniques |
US11/545,256 Abandoned US20070030982A1 (en) | 2000-05-10 | 2006-10-10 | Interference suppression techniques |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/290,137 Expired - Fee Related US7613309B2 (en) | 2000-05-10 | 2002-11-07 | Interference suppression techniques |
Country Status (9)
Country | Link |
---|---|
US (2) | US7613309B2 (en) |
EP (1) | EP1312239B1 (en) |
JP (1) | JP2003533152A (en) |
CN (1) | CN1440628A (en) |
AU (1) | AU2001261344A1 (en) |
CA (2) | CA2685434A1 (en) |
DE (1) | DE60125553T2 (en) |
DK (1) | DK1312239T3 (en) |
WO (1) | WO2001087011A2 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008106680A2 (en) * | 2007-03-01 | 2008-09-04 | Jerry Mahabub | Audio spatialization and environment simulation |
US20090106021A1 (en) * | 2007-10-18 | 2009-04-23 | Motorola, Inc. | Robust two microphone noise suppression system |
US20100027821A1 (en) * | 2006-10-10 | 2010-02-04 | Eghart Fischer | Hearing aid and method for operating a hearing aid |
US20100246831A1 (en) * | 2008-10-20 | 2010-09-30 | Jerry Mahabub | Audio spatialization and environment simulation |
US20130030801A1 (en) * | 2011-07-29 | 2013-01-31 | Qnx Software Systems Limited | Off-Axis Audio Suppressions in An Automobile Cabin |
US8638946B1 (en) | 2004-03-16 | 2014-01-28 | Genaudio, Inc. | Method and apparatus for creating spatialized sound |
EP2833651A1 (en) * | 2013-08-01 | 2015-02-04 | Siemens Medical Instruments Pte. Ltd. | Method for tracking of a sound source |
Families Citing this family (53)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7720229B2 (en) * | 2002-11-08 | 2010-05-18 | University Of Maryland | Method for measurement of head related transfer functions |
US7945064B2 (en) | 2003-04-09 | 2011-05-17 | Board Of Trustees Of The University Of Illinois | Intrabody communication with ultrasound |
US7076072B2 (en) * | 2003-04-09 | 2006-07-11 | Board Of Trustees For The University Of Illinois | Systems and methods for interference-suppression with directional sensing patterns |
EP1524879B1 (en) | 2003-06-30 | 2014-05-07 | Nuance Communications, Inc. | Handsfree system for use in a vehicle |
GB0321722D0 (en) * | 2003-09-16 | 2003-10-15 | Mitel Networks Corp | A method for optimal microphone array design under uniform acoustic coupling constraints |
US7283639B2 (en) * | 2004-03-10 | 2007-10-16 | Starkey Laboratories, Inc. | Hearing instrument with data transmission interference blocking |
WO2005109951A1 (en) * | 2004-05-05 | 2005-11-17 | Deka Products Limited Partnership | Angular discrimination of acoustical or radio signals |
US8139787B2 (en) * | 2005-09-09 | 2012-03-20 | Simon Haykin | Method and device for binaural signal enhancement |
US8345890B2 (en) | 2006-01-05 | 2013-01-01 | Audience, Inc. | System and method for utilizing inter-microphone level differences for speech enhancement |
US8194880B2 (en) * | 2006-01-30 | 2012-06-05 | Audience, Inc. | System and method for utilizing omni-directional microphones for speech enhancement |
US8744844B2 (en) | 2007-07-06 | 2014-06-03 | Audience, Inc. | System and method for adaptive intelligent noise suppression |
US9185487B2 (en) | 2006-01-30 | 2015-11-10 | Audience, Inc. | System and method for providing noise suppression utilizing null processing noise subtraction |
US8204252B1 (en) | 2006-10-10 | 2012-06-19 | Audience, Inc. | System and method for providing close microphone adaptive array processing |
DE102006018634B4 (en) * | 2006-04-21 | 2017-12-07 | Sivantos Gmbh | Hearing aid with source separation and corresponding method |
US8204253B1 (en) | 2008-06-30 | 2012-06-19 | Audience, Inc. | Self calibration of audio device |
US8150065B2 (en) | 2006-05-25 | 2012-04-03 | Audience, Inc. | System and method for processing an audio signal |
US8934641B2 (en) | 2006-05-25 | 2015-01-13 | Audience, Inc. | Systems and methods for reconstructing decomposed audio signals |
US8849231B1 (en) | 2007-08-08 | 2014-09-30 | Audience, Inc. | System and method for adaptive power control |
US8949120B1 (en) | 2006-05-25 | 2015-02-03 | Audience, Inc. | Adaptive noise cancelation |
ATE430975T1 (en) * | 2006-07-10 | 2009-05-15 | Harman Becker Automotive Sys | REDUCING BACKGROUND NOISE IN HANDS-FREE SYSTEMS |
JP5070873B2 (en) * | 2006-08-09 | 2012-11-14 | 富士通株式会社 | Sound source direction estimating apparatus, sound source direction estimating method, and computer program |
DE102006047982A1 (en) * | 2006-10-10 | 2008-04-24 | Siemens Audiologische Technik Gmbh | Method for operating a hearing aid, and hearing aid |
AU2007306432B2 (en) * | 2006-10-10 | 2012-03-29 | Sivantos Gmbh | Method for operating a hearing aid, and hearing aid |
DE102006047983A1 (en) | 2006-10-10 | 2008-04-24 | Siemens Audiologische Technik Gmbh | Processing an input signal in a hearing aid |
EP1912472A1 (en) * | 2006-10-10 | 2008-04-16 | Siemens Audiologische Technik GmbH | Method for operating a hearing aid and hearing aid |
JP4854533B2 (en) * | 2007-01-30 | 2012-01-18 | 富士通株式会社 | Acoustic judgment method, acoustic judgment device, and computer program |
US8259926B1 (en) | 2007-02-23 | 2012-09-04 | Audience, Inc. | System and method for 2-channel and 3-channel acoustic echo cancellation |
US8189766B1 (en) | 2007-07-26 | 2012-05-29 | Audience, Inc. | System and method for blind subband acoustic echo cancellation postfiltering |
US8352274B2 (en) | 2007-09-11 | 2013-01-08 | Panasonic Corporation | Sound determination device, sound detection device, and sound determination method for determining frequency signals of a to-be-extracted sound included in a mixed sound |
GB0720473D0 (en) * | 2007-10-19 | 2007-11-28 | Univ Surrey | Accoustic source separation |
US8180064B1 (en) | 2007-12-21 | 2012-05-15 | Audience, Inc. | System and method for providing voice equalization |
US8143620B1 (en) | 2007-12-21 | 2012-03-27 | Audience, Inc. | System and method for adaptive classification of audio sources |
US8194882B2 (en) | 2008-02-29 | 2012-06-05 | Audience, Inc. | System and method for providing single microphone noise suppression fallback |
US8355511B2 (en) | 2008-03-18 | 2013-01-15 | Audience, Inc. | System and method for envelope-based acoustic echo cancellation |
WO2009151578A2 (en) | 2008-06-09 | 2009-12-17 | The Board Of Trustees Of The University Of Illinois | Method and apparatus for blind signal recovery in noisy, reverberant environments |
US8521530B1 (en) | 2008-06-30 | 2013-08-27 | Audience, Inc. | System and method for enhancing a monaural audio signal |
US8774423B1 (en) | 2008-06-30 | 2014-07-08 | Audience, Inc. | System and method for controlling adaptivity of signal modification using a phantom coefficient |
TWI475896B (en) * | 2008-09-25 | 2015-03-01 | Dolby Lab Licensing Corp | Binaural filters for monophonic compatibility and loudspeaker compatibility |
DK2211579T3 (en) * | 2009-01-21 | 2012-10-08 | Oticon As | Transmission power control in a low power wireless communication system |
US9838784B2 (en) * | 2009-12-02 | 2017-12-05 | Knowles Electronics, Llc | Directional audio capture |
US9008329B1 (en) | 2010-01-26 | 2015-04-14 | Audience, Inc. | Noise reduction using multi-feature cluster tracker |
US8798290B1 (en) | 2010-04-21 | 2014-08-05 | Audience, Inc. | Systems and methods for adaptive signal equalization |
US9640194B1 (en) | 2012-10-04 | 2017-05-02 | Knowles Electronics, Llc | Noise suppression for speech processing based on machine-learning mask estimation |
US9078057B2 (en) | 2012-11-01 | 2015-07-07 | Csr Technology Inc. | Adaptive microphone beamforming |
US20140270219A1 (en) * | 2013-03-15 | 2014-09-18 | CSR Technology, Inc. | Method, apparatus, and manufacture for beamforming with fixed weights and adaptive selection or resynthesis |
US9536540B2 (en) | 2013-07-19 | 2017-01-03 | Knowles Electronics, Llc | Speech signal separation and synthesis based on auditory scene analysis and speech modeling |
EP2928210A1 (en) * | 2014-04-03 | 2015-10-07 | Oticon A/s | A binaural hearing assistance system comprising binaural noise reduction |
DE112015003945T5 (en) | 2014-08-28 | 2017-05-11 | Knowles Electronics, Llc | Multi-source noise reduction |
US9875081B2 (en) * | 2015-09-21 | 2018-01-23 | Amazon Technologies, Inc. | Device selection for providing a response |
DE102017206788B3 (en) * | 2017-04-21 | 2018-08-02 | Sivantos Pte. Ltd. | Method for operating a hearing aid |
US10482904B1 (en) | 2017-08-15 | 2019-11-19 | Amazon Technologies, Inc. | Context driven device arbitration |
CN110070709B (en) * | 2019-05-29 | 2023-10-27 | 杭州聚声科技有限公司 | Pedestrian crossing directional voice prompt system and method thereof |
CN115751737B (en) * | 2023-01-09 | 2023-04-25 | 南通源动太阳能科技有限公司 | Dish type heat collection heater for solar thermal power generation system and design method |
Citations (99)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4025721A (en) * | 1976-05-04 | 1977-05-24 | Biocommunications Research Corporation | Method of and means for adaptively filtering near-stationary noise from speech |
US4207441A (en) * | 1977-03-16 | 1980-06-10 | Bertin & Cie | Auditory prosthesis equipment |
US4334740A (en) * | 1978-09-12 | 1982-06-15 | Polaroid Corporation | Receiving system having pre-selected directional response |
US4354064A (en) * | 1980-02-19 | 1982-10-12 | Scott Instruments Company | Vibratory aid for presbycusis |
US4536887A (en) * | 1982-10-18 | 1985-08-20 | Nippon Telegraph & Telephone Public Corporation | Microphone-array apparatus and method for extracting desired signal |
US4611598A (en) * | 1984-05-30 | 1986-09-16 | Hortmann Gmbh | Multi-frequency transmission system for implanted hearing aids |
US4653606A (en) * | 1985-03-22 | 1987-03-31 | American Telephone And Telegraph Company | Electroacoustic device with broad frequency range directional response |
US4703506A (en) * | 1985-07-23 | 1987-10-27 | Victor Company Of Japan, Ltd. | Directional microphone apparatus |
US4742548A (en) * | 1984-12-20 | 1988-05-03 | American Telephone And Telegraph Company | Unidirectional second order gradient microphone |
US4752961A (en) * | 1985-09-23 | 1988-06-21 | Northern Telecom Limited | Microphone arrangement |
US4773095A (en) * | 1985-10-16 | 1988-09-20 | Siemens Aktiengesellschaft | Hearing aid with locating microphones |
US4845755A (en) * | 1984-08-28 | 1989-07-04 | Siemens Aktiengesellschaft | Remote control hearing aid |
US4858612A (en) * | 1983-12-19 | 1989-08-22 | Stocklin Philip L | Hearing device |
US4918737A (en) * | 1987-07-07 | 1990-04-17 | Siemens Aktiengesellschaft | Hearing aid with wireless remote control |
US4982434A (en) * | 1989-05-30 | 1991-01-01 | Center For Innovative Technology | Supersonic bone conduction hearing aid and method |
US4987897A (en) * | 1989-09-18 | 1991-01-29 | Medtronic, Inc. | Body bus medical device communication system |
US4988981A (en) * | 1987-03-17 | 1991-01-29 | Vpl Research, Inc. | Computer data entry and manipulation apparatus and method |
US5012520A (en) * | 1988-05-06 | 1991-04-30 | Siemens Aktiengesellschaft | Hearing aid with wireless remote control |
US5029216A (en) * | 1989-06-09 | 1991-07-02 | The United States Of America As Represented By The Administrator Of The National Aeronautics & Space Administration | Visual aid for the hearing impaired |
US5040156A (en) * | 1989-06-29 | 1991-08-13 | Battelle-Institut E.V. | Acoustic sensor device with noise suppression |
US5047994A (en) * | 1989-05-30 | 1991-09-10 | Center For Innovative Technology | Supersonic bone conduction hearing aid and method |
US5113859A (en) * | 1988-09-19 | 1992-05-19 | Medtronic, Inc. | Acoustic body bus medical device communication system |
US5245556A (en) * | 1992-09-15 | 1993-09-14 | Universal Data Systems, Inc. | Adaptive equalizer method and apparatus |
US5259032A (en) * | 1990-11-07 | 1993-11-02 | Resound Corporation | contact transducer assembly for hearing devices |
US5285499A (en) * | 1993-04-27 | 1994-02-08 | Signal Science, Inc. | Ultrasonic frequency expansion processor |
US5289544A (en) * | 1991-12-31 | 1994-02-22 | Audiological Engineering Corporation | Method and apparatus for reducing background noise in communication systems and for enhancing binaural hearing systems for the hearing impaired |
US5321332A (en) * | 1992-11-12 | 1994-06-14 | The Whitaker Corporation | Wideband ultrasonic transducer |
US5325436A (en) * | 1993-06-30 | 1994-06-28 | House Ear Institute | Method of signal processing for maintaining directional hearing with hearing aids |
US5383915A (en) * | 1991-04-10 | 1995-01-24 | Angeion Corporation | Wireless programmer/repeater system for an implanted medical device |
US5400409A (en) * | 1992-12-23 | 1995-03-21 | Daimler-Benz Ag | Noise-reduction method for noise-affected voice channels |
US5417113A (en) * | 1993-08-18 | 1995-05-23 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Leak detection utilizing analog binaural (VLSI) techniques |
US5430690A (en) * | 1992-03-20 | 1995-07-04 | Abel; Jonathan S. | Method and apparatus for processing signals to extract narrow bandwidth features |
US5454838A (en) * | 1992-07-27 | 1995-10-03 | Sorin Biomedica S.P.A. | Method and a device for monitoring heart function |
US5463694A (en) * | 1993-11-01 | 1995-10-31 | Motorola | Gradient directional microphone system and method therefor |
US5485515A (en) * | 1993-12-29 | 1996-01-16 | At&T Corp. | Background noise compensation in a telephone network |
US5495534A (en) * | 1990-01-19 | 1996-02-27 | Sony Corporation | Audio signal reproducing apparatus |
US5507781A (en) * | 1991-05-23 | 1996-04-16 | Angeion Corporation | Implantable defibrillator system with capacitor switching circuitry |
US5511128A (en) * | 1994-01-21 | 1996-04-23 | Lindemann; Eric | Dynamic intensity beamforming system for noise reduction in a binaural hearing aid |
US5550923A (en) * | 1994-09-02 | 1996-08-27 | Minnesota Mining And Manufacturing Company | Directional ear device with adaptive bandwidth and gain control |
US5574824A (en) * | 1994-04-11 | 1996-11-12 | The United States Of America As Represented By The Secretary Of The Air Force | Analysis/synthesis-based microphone array speech enhancer with variable signal distortion |
US5627799A (en) * | 1994-09-01 | 1997-05-06 | Nec Corporation | Beamformer using coefficient restrained adaptive filters for detecting interference signals |
US5651071A (en) * | 1993-09-17 | 1997-07-22 | Audiologic, Inc. | Noise reduction system for binaural hearing aid |
US5663727A (en) * | 1995-06-23 | 1997-09-02 | Hearing Innovations Incorporated | Frequency response analyzer and shaping apparatus and digital hearing enhancement apparatus and method utilizing the same |
US5706352A (en) * | 1993-04-07 | 1998-01-06 | K/S Himpp | Adaptive gain and filtering circuit for a sound reproduction system |
US5715319A (en) * | 1996-05-30 | 1998-02-03 | Picturetel Corporation | Method and apparatus for steerable and endfire superdirective microphone arrays with reduced analog-to-digital converter and computational requirements |
US5721783A (en) * | 1995-06-07 | 1998-02-24 | Anderson; James C. | Hearing aid with wireless remote processor |
US5734976A (en) * | 1994-03-07 | 1998-03-31 | Phonak Communications Ag | Micro-receiver for receiving a high frequency frequency-modulated or phase-modulated signal |
US5737430A (en) * | 1993-07-22 | 1998-04-07 | Cardinal Sound Labs, Inc. | Directional hearing aid |
US5755748A (en) * | 1996-07-24 | 1998-05-26 | Dew Engineering & Development Limited | Transcutaneous energy transfer device |
US5757932A (en) * | 1993-09-17 | 1998-05-26 | Audiologic, Inc. | Digital hearing aid system |
US5768392A (en) * | 1996-04-16 | 1998-06-16 | Aura Systems Inc. | Blind adaptive filtering of unknown signals in unknown noise in quasi-closed loop system |
US5793875A (en) * | 1996-04-22 | 1998-08-11 | Cardinal Sound Labs, Inc. | Directional hearing system |
US5825898A (en) * | 1996-06-27 | 1998-10-20 | Lamar Signal Processing Ltd. | System and method for adaptive interference cancelling |
US5831936A (en) * | 1995-02-21 | 1998-11-03 | State Of Israel/Ministry Of Defense Armament Development Authority - Rafael | System and method of noise detection |
US5833603A (en) * | 1996-03-13 | 1998-11-10 | Lipomatrix, Inc. | Implantable biosensing transponder |
US5878147A (en) * | 1996-12-31 | 1999-03-02 | Etymotic Research, Inc. | Directional microphone assembly |
US5889870A (en) * | 1996-07-17 | 1999-03-30 | American Technology Corporation | Acoustic heterodyne device and method |
US5991419A (en) * | 1997-04-29 | 1999-11-23 | Beltone Electronics Corporation | Bilateral signal processing prosthesis |
US6010532A (en) * | 1996-11-25 | 2000-01-04 | St. Croix Medical, Inc. | Dual path implantable hearing assistance device |
US6023514A (en) * | 1997-12-22 | 2000-02-08 | Strandberg; Malcolm W. P. | System and method for factoring a merged wave field into independent components |
US6068589A (en) * | 1996-02-15 | 2000-05-30 | Neukermans; Armand P. | Biocompatible fully implantable hearing aid transducers |
US6094150A (en) * | 1997-09-10 | 2000-07-25 | Mitsubishi Heavy Industries, Ltd. | System and method of measuring noise of mobile body using a plurality microphones |
US6104822A (en) * | 1995-10-10 | 2000-08-15 | Audiologic, Inc. | Digital signal processing hearing aid |
US6118882A (en) * | 1995-01-25 | 2000-09-12 | Haynes; Philip Ashley | Communication method |
US6137889A (en) * | 1998-05-27 | 2000-10-24 | Insonus Medical, Inc. | Direct tympanic membrane excitation via vibrationally conductive assembly |
US6141591A (en) * | 1996-03-06 | 2000-10-31 | Advanced Bionics Corporation | Magnetless implantable stimulator and external transmitter and implant tools for aligning same |
US6154552A (en) * | 1997-05-15 | 2000-11-28 | Planning Systems Inc. | Hybrid adaptive beamformer |
US6173062B1 (en) * | 1994-03-16 | 2001-01-09 | Hearing Innovations Incorporated | Frequency transpositional hearing aid with digital and single sideband modulation |
US6182018B1 (en) * | 1998-08-25 | 2001-01-30 | Ford Global Technologies, Inc. | Method and apparatus for identifying sound in a composite sound signal |
US6192134B1 (en) * | 1997-11-20 | 2001-02-20 | Conexant Systems, Inc. | System and method for a monolithic directional microphone array |
US6198693B1 (en) * | 1998-04-13 | 2001-03-06 | Andrea Electronics Corporation | System and method for finding the direction of a wave source using an array of sensors |
US6217508B1 (en) * | 1998-08-14 | 2001-04-17 | Symphonix Devices, Inc. | Ultrasonic hearing system |
US6223018B1 (en) * | 1996-12-12 | 2001-04-24 | Nippon Telegraph And Telephone Corporation | Intra-body information transfer device |
US6222927B1 (en) * | 1996-06-19 | 2001-04-24 | The University Of Illinois | Binaural signal processing system and method |
US6229900B1 (en) * | 1997-07-18 | 2001-05-08 | Beltone Netherlands B.V. | Hearing aid including a programmable processor |
US6243471B1 (en) * | 1995-03-07 | 2001-06-05 | Brown University Research Foundation | Methods and apparatus for source location estimation from microphone-array time-delay estimates |
US6261224B1 (en) * | 1996-08-07 | 2001-07-17 | St. Croix Medical, Inc. | Piezoelectric film transducer for cochlear prosthetic |
US6272229B1 (en) * | 1999-08-03 | 2001-08-07 | Topholm & Westermann Aps | Hearing aid with adaptive matching of microphones |
US6275596B1 (en) * | 1997-01-10 | 2001-08-14 | Gn Resound Corporation | Open ear canal hearing aid system |
US6283915B1 (en) * | 1997-03-12 | 2001-09-04 | Sarnoff Corporation | Disposable in-the-ear monitoring instrument and method of manufacture |
US20010031996A1 (en) * | 2000-04-13 | 2001-10-18 | Hans Leysieffer | At least partially implantable system for rehabilitation of a hearing disorder |
US6307945B1 (en) * | 1990-12-21 | 2001-10-23 | Sense-Sonic Limited | Radio-based hearing aid system |
US6342035B1 (en) * | 1999-02-05 | 2002-01-29 | St. Croix Medical, Inc. | Hearing assistance device sensing otovibratory or otoacoustic emissions evoked by middle ear vibrations |
US20020012438A1 (en) * | 2000-06-30 | 2002-01-31 | Hans Leysieffer | System for rehabilitation of a hearing disorder |
US20020019668A1 (en) * | 2000-08-11 | 2002-02-14 | Friedemann Stockert | At least partially implantable system for rehabilitation of a hearing disorder |
US20020029070A1 (en) * | 2000-04-13 | 2002-03-07 | Hans Leysieffer | At least partially implantable system for rehabilitation a hearing disorder |
US6380896B1 (en) * | 2000-10-30 | 2002-04-30 | Siemens Information And Communication Mobile, Llc | Circular polarization antenna for wireless communication system |
US6385323B1 (en) * | 1998-05-15 | 2002-05-07 | Siemens Audiologische Technik Gmbh | Hearing aid with automatic microphone balancing and method for operating a hearing aid with automatic microphone balancing |
US6389142B1 (en) * | 1996-12-11 | 2002-05-14 | Micro Ear Technology | In-the-ear hearing aid with directional microphone system |
US20020057817A1 (en) * | 2000-10-10 | 2002-05-16 | Resistance Technology, Inc. | Hearing aid |
US6390971B1 (en) * | 1999-02-05 | 2002-05-21 | St. Croix Medical, Inc. | Method and apparatus for a programmable implantable hearing aid |
US6397186B1 (en) * | 1999-12-22 | 2002-05-28 | Ambush Interactive, Inc. | Hands-free, voice-operated remote control transmitter |
US6421448B1 (en) * | 1999-04-26 | 2002-07-16 | Siemens Audiologische Technik Gmbh | Hearing aid with a directional microphone characteristic and method for producing same |
US6424721B1 (en) * | 1998-03-09 | 2002-07-23 | Siemens Audiologische Technik Gmbh | Hearing aid with a directional microphone system as well as method for the operation thereof |
US20020110255A1 (en) * | 2000-10-05 | 2002-08-15 | Killion Mead C. | Directional microphone assembly |
US20020131613A1 (en) * | 2001-03-13 | 2002-09-19 | Andreas Jakob | Method for establishing a binaural communication link and binaural hearing devices |
US20020141595A1 (en) * | 2001-02-23 | 2002-10-03 | Jouppi Norman P. | System and method for audio telepresence |
US6549586B2 (en) * | 1999-04-12 | 2003-04-15 | Telefonaktiebolaget L M Ericsson | System and method for dual microphone signal noise reduction using spectral subtraction |
US6717991B1 (en) * | 1998-05-27 | 2004-04-06 | Telefonaktiebolaget Lm Ericsson (Publ) | System and method for dual microphone signal noise reduction using spectral subtraction |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1105565A (en) | 1978-09-12 | 1981-07-21 | Kaufman (John G.) Hospital Products Ltd. | Electrosurgical electrode |
DE2924539C2 (en) | 1979-06-19 | 1983-01-13 | Fa. Carl Freudenberg, 6940 Weinheim | Polyolefin filament spunbond and process for its manufacture |
JPS5939198A (en) | 1982-08-27 | 1984-03-03 | Victor Co Of Japan Ltd | Microphone device |
AT379929B (en) | 1984-07-18 | 1986-03-10 | Viennatone Gmbh | HOERGERAET |
US5524056A (en) | 1993-04-13 | 1996-06-04 | Etymotic Research, Inc. | Hearing aid having plural microphones and a microphone switching system |
US5479522A (en) | 1993-09-17 | 1995-12-26 | Audiologic, Inc. | Binaural hearing aid |
US5473701A (en) | 1993-11-05 | 1995-12-05 | At&T Corp. | Adaptive microphone array |
US6002776A (en) | 1995-09-18 | 1999-12-14 | Interval Research Corporation | Directional acoustic signal processor and method therefor |
US5694474A (en) | 1995-09-18 | 1997-12-02 | Interval Research Corporation | Adaptive filter for signal processing and method therefor |
US6161046A (en) | 1996-04-09 | 2000-12-12 | Maniglia; Anthony J. | Totally implantable cochlear implant for improvement of partial and total sensorineural hearing loss |
US6317703B1 (en) | 1996-11-12 | 2001-11-13 | International Business Machines Corporation | Separation of a mixture of acoustic sources into its components |
US6178248B1 (en) | 1997-04-14 | 2001-01-23 | Andrea Electronics Corporation | Dual-processing interference cancelling system and method |
EP0802699A3 (en) * | 1997-07-16 | 1998-02-25 | Phonak Ag | Method for electronically enlarging the distance between two acoustical/electrical transducers and hearing aid apparatus |
FR2768290B1 (en) | 1997-09-10 | 1999-10-15 | France Telecom | ANTENNA FORMED OF A PLURALITY OF ACOUSTIC SENSORS |
US20010051776A1 (en) * | 1998-10-14 | 2001-12-13 | Lenhardt Martin L. | Tinnitus masker/suppressor |
US6167312A (en) | 1999-04-30 | 2000-12-26 | Medtronic, Inc. | Telemetry system for implantable medical devices |
ATE417483T1 (en) * | 2000-02-02 | 2008-12-15 | Bernafon Ag | CIRCUIT AND METHOD FOR ADAPTIVE NOISE CANCELLATION |
DE10018360C2 (en) * | 2000-04-13 | 2002-10-10 | Cochlear Ltd | At least partially implantable system for the rehabilitation of a hearing impairment |
-
2001
- 2001-05-10 AU AU2001261344A patent/AU2001261344A1/en not_active Abandoned
- 2001-05-10 EP EP01935234A patent/EP1312239B1/en not_active Expired - Lifetime
- 2001-05-10 DK DK01935234T patent/DK1312239T3/en active
- 2001-05-10 CN CN01812199A patent/CN1440628A/en active Pending
- 2001-05-10 CA CA2685434A patent/CA2685434A1/en not_active Abandoned
- 2001-05-10 CA CA002407855A patent/CA2407855C/en not_active Expired - Fee Related
- 2001-05-10 DE DE60125553T patent/DE60125553T2/en not_active Expired - Lifetime
- 2001-05-10 WO PCT/US2001/015047 patent/WO2001087011A2/en active IP Right Grant
- 2001-05-10 JP JP2001583102A patent/JP2003533152A/en active Pending
-
2002
- 2002-11-07 US US10/290,137 patent/US7613309B2/en not_active Expired - Fee Related
-
2006
- 2006-10-10 US US11/545,256 patent/US20070030982A1/en not_active Abandoned
Patent Citations (100)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4025721A (en) * | 1976-05-04 | 1977-05-24 | Biocommunications Research Corporation | Method of and means for adaptively filtering near-stationary noise from speech |
US4207441A (en) * | 1977-03-16 | 1980-06-10 | Bertin & Cie | Auditory prosthesis equipment |
US4334740A (en) * | 1978-09-12 | 1982-06-15 | Polaroid Corporation | Receiving system having pre-selected directional response |
US4354064A (en) * | 1980-02-19 | 1982-10-12 | Scott Instruments Company | Vibratory aid for presbycusis |
US4536887A (en) * | 1982-10-18 | 1985-08-20 | Nippon Telegraph & Telephone Public Corporation | Microphone-array apparatus and method for extracting desired signal |
US4858612A (en) * | 1983-12-19 | 1989-08-22 | Stocklin Philip L | Hearing device |
US4611598A (en) * | 1984-05-30 | 1986-09-16 | Hortmann Gmbh | Multi-frequency transmission system for implanted hearing aids |
US4845755A (en) * | 1984-08-28 | 1989-07-04 | Siemens Aktiengesellschaft | Remote control hearing aid |
US4742548A (en) * | 1984-12-20 | 1988-05-03 | American Telephone And Telegraph Company | Unidirectional second order gradient microphone |
US4653606A (en) * | 1985-03-22 | 1987-03-31 | American Telephone And Telegraph Company | Electroacoustic device with broad frequency range directional response |
US4703506A (en) * | 1985-07-23 | 1987-10-27 | Victor Company Of Japan, Ltd. | Directional microphone apparatus |
US4752961A (en) * | 1985-09-23 | 1988-06-21 | Northern Telecom Limited | Microphone arrangement |
US4773095A (en) * | 1985-10-16 | 1988-09-20 | Siemens Aktiengesellschaft | Hearing aid with locating microphones |
US4988981A (en) * | 1987-03-17 | 1991-01-29 | Vpl Research, Inc. | Computer data entry and manipulation apparatus and method |
US4988981B1 (en) * | 1987-03-17 | 1999-05-18 | Vpl Newco Inc | Computer data entry and manipulation apparatus and method |
US4918737A (en) * | 1987-07-07 | 1990-04-17 | Siemens Aktiengesellschaft | Hearing aid with wireless remote control |
US5012520A (en) * | 1988-05-06 | 1991-04-30 | Siemens Aktiengesellschaft | Hearing aid with wireless remote control |
US5113859A (en) * | 1988-09-19 | 1992-05-19 | Medtronic, Inc. | Acoustic body bus medical device communication system |
US5047994A (en) * | 1989-05-30 | 1991-09-10 | Center For Innovative Technology | Supersonic bone conduction hearing aid and method |
US4982434A (en) * | 1989-05-30 | 1991-01-01 | Center For Innovative Technology | Supersonic bone conduction hearing aid and method |
US5029216A (en) * | 1989-06-09 | 1991-07-02 | The United States Of America As Represented By The Administrator Of The National Aeronautics & Space Administration | Visual aid for the hearing impaired |
US5040156A (en) * | 1989-06-29 | 1991-08-13 | Battelle-Institut E.V. | Acoustic sensor device with noise suppression |
US4987897A (en) * | 1989-09-18 | 1991-01-29 | Medtronic, Inc. | Body bus medical device communication system |
US5495534A (en) * | 1990-01-19 | 1996-02-27 | Sony Corporation | Audio signal reproducing apparatus |
US5259032A (en) * | 1990-11-07 | 1993-11-02 | Resound Corporation | contact transducer assembly for hearing devices |
US6307945B1 (en) * | 1990-12-21 | 2001-10-23 | Sense-Sonic Limited | Radio-based hearing aid system |
US5383915A (en) * | 1991-04-10 | 1995-01-24 | Angeion Corporation | Wireless programmer/repeater system for an implanted medical device |
US5507781A (en) * | 1991-05-23 | 1996-04-16 | Angeion Corporation | Implantable defibrillator system with capacitor switching circuitry |
US5289544A (en) * | 1991-12-31 | 1994-02-22 | Audiological Engineering Corporation | Method and apparatus for reducing background noise in communication systems and for enhancing binaural hearing systems for the hearing impaired |
US5430690A (en) * | 1992-03-20 | 1995-07-04 | Abel; Jonathan S. | Method and apparatus for processing signals to extract narrow bandwidth features |
US5454838A (en) * | 1992-07-27 | 1995-10-03 | Sorin Biomedica S.P.A. | Method and a device for monitoring heart function |
US5245556A (en) * | 1992-09-15 | 1993-09-14 | Universal Data Systems, Inc. | Adaptive equalizer method and apparatus |
US5321332A (en) * | 1992-11-12 | 1994-06-14 | The Whitaker Corporation | Wideband ultrasonic transducer |
US5400409A (en) * | 1992-12-23 | 1995-03-21 | Daimler-Benz Ag | Noise-reduction method for noise-affected voice channels |
US5706352A (en) * | 1993-04-07 | 1998-01-06 | K/S Himpp | Adaptive gain and filtering circuit for a sound reproduction system |
US5285499A (en) * | 1993-04-27 | 1994-02-08 | Signal Science, Inc. | Ultrasonic frequency expansion processor |
US5325436A (en) * | 1993-06-30 | 1994-06-28 | House Ear Institute | Method of signal processing for maintaining directional hearing with hearing aids |
US5737430A (en) * | 1993-07-22 | 1998-04-07 | Cardinal Sound Labs, Inc. | Directional hearing aid |
US5417113A (en) * | 1993-08-18 | 1995-05-23 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Leak detection utilizing analog binaural (VLSI) techniques |
US5757932A (en) * | 1993-09-17 | 1998-05-26 | Audiologic, Inc. | Digital hearing aid system |
US5651071A (en) * | 1993-09-17 | 1997-07-22 | Audiologic, Inc. | Noise reduction system for binaural hearing aid |
US5463694A (en) * | 1993-11-01 | 1995-10-31 | Motorola | Gradient directional microphone system and method therefor |
US5485515A (en) * | 1993-12-29 | 1996-01-16 | At&T Corp. | Background noise compensation in a telephone network |
US5511128A (en) * | 1994-01-21 | 1996-04-23 | Lindemann; Eric | Dynamic intensity beamforming system for noise reduction in a binaural hearing aid |
US5734976A (en) * | 1994-03-07 | 1998-03-31 | Phonak Communications Ag | Micro-receiver for receiving a high frequency frequency-modulated or phase-modulated signal |
US6173062B1 (en) * | 1994-03-16 | 2001-01-09 | Hearing Innovations Incorporated | Frequency transpositional hearing aid with digital and single sideband modulation |
US5574824A (en) * | 1994-04-11 | 1996-11-12 | The United States Of America As Represented By The Secretary Of The Air Force | Analysis/synthesis-based microphone array speech enhancer with variable signal distortion |
US5627799A (en) * | 1994-09-01 | 1997-05-06 | Nec Corporation | Beamformer using coefficient restrained adaptive filters for detecting interference signals |
US5550923A (en) * | 1994-09-02 | 1996-08-27 | Minnesota Mining And Manufacturing Company | Directional ear device with adaptive bandwidth and gain control |
US6118882A (en) * | 1995-01-25 | 2000-09-12 | Haynes; Philip Ashley | Communication method |
US5831936A (en) * | 1995-02-21 | 1998-11-03 | State Of Israel/Ministry Of Defense Armament Development Authority - Rafael | System and method of noise detection |
US6243471B1 (en) * | 1995-03-07 | 2001-06-05 | Brown University Research Foundation | Methods and apparatus for source location estimation from microphone-array time-delay estimates |
US5721783A (en) * | 1995-06-07 | 1998-02-24 | Anderson; James C. | Hearing aid with wireless remote processor |
US5663727A (en) * | 1995-06-23 | 1997-09-02 | Hearing Innovations Incorporated | Frequency response analyzer and shaping apparatus and digital hearing enhancement apparatus and method utilizing the same |
US6104822A (en) * | 1995-10-10 | 2000-08-15 | Audiologic, Inc. | Digital signal processing hearing aid |
US6068589A (en) * | 1996-02-15 | 2000-05-30 | Neukermans; Armand P. | Biocompatible fully implantable hearing aid transducers |
US6141591A (en) * | 1996-03-06 | 2000-10-31 | Advanced Bionics Corporation | Magnetless implantable stimulator and external transmitter and implant tools for aligning same |
US5833603A (en) * | 1996-03-13 | 1998-11-10 | Lipomatrix, Inc. | Implantable biosensing transponder |
US5768392A (en) * | 1996-04-16 | 1998-06-16 | Aura Systems Inc. | Blind adaptive filtering of unknown signals in unknown noise in quasi-closed loop system |
US5793875A (en) * | 1996-04-22 | 1998-08-11 | Cardinal Sound Labs, Inc. | Directional hearing system |
US5715319A (en) * | 1996-05-30 | 1998-02-03 | Picturetel Corporation | Method and apparatus for steerable and endfire superdirective microphone arrays with reduced analog-to-digital converter and computational requirements |
US6222927B1 (en) * | 1996-06-19 | 2001-04-24 | The University Of Illinois | Binaural signal processing system and method |
US5825898A (en) * | 1996-06-27 | 1998-10-20 | Lamar Signal Processing Ltd. | System and method for adaptive interference cancelling |
US5889870A (en) * | 1996-07-17 | 1999-03-30 | American Technology Corporation | Acoustic heterodyne device and method |
US5755748A (en) * | 1996-07-24 | 1998-05-26 | Dew Engineering & Development Limited | Transcutaneous energy transfer device |
US6261224B1 (en) * | 1996-08-07 | 2001-07-17 | St. Croix Medical, Inc. | Piezoelectric film transducer for cochlear prosthetic |
US6010532A (en) * | 1996-11-25 | 2000-01-04 | St. Croix Medical, Inc. | Dual path implantable hearing assistance device |
US6389142B1 (en) * | 1996-12-11 | 2002-05-14 | Micro Ear Technology | In-the-ear hearing aid with directional microphone system |
US6223018B1 (en) * | 1996-12-12 | 2001-04-24 | Nippon Telegraph And Telephone Corporation | Intra-body information transfer device |
US5878147A (en) * | 1996-12-31 | 1999-03-02 | Etymotic Research, Inc. | Directional microphone assembly |
US6275596B1 (en) * | 1997-01-10 | 2001-08-14 | Gn Resound Corporation | Open ear canal hearing aid system |
US6283915B1 (en) * | 1997-03-12 | 2001-09-04 | Sarnoff Corporation | Disposable in-the-ear monitoring instrument and method of manufacture |
US5991419A (en) * | 1997-04-29 | 1999-11-23 | Beltone Electronics Corporation | Bilateral signal processing prosthesis |
US6154552A (en) * | 1997-05-15 | 2000-11-28 | Planning Systems Inc. | Hybrid adaptive beamformer |
US6229900B1 (en) * | 1997-07-18 | 2001-05-08 | Beltone Netherlands B.V. | Hearing aid including a programmable processor |
US6094150A (en) * | 1997-09-10 | 2000-07-25 | Mitsubishi Heavy Industries, Ltd. | System and method of measuring noise of mobile body using a plurality microphones |
US6192134B1 (en) * | 1997-11-20 | 2001-02-20 | Conexant Systems, Inc. | System and method for a monolithic directional microphone array |
US6023514A (en) * | 1997-12-22 | 2000-02-08 | Strandberg; Malcolm W. P. | System and method for factoring a merged wave field into independent components |
US6424721B1 (en) * | 1998-03-09 | 2002-07-23 | Siemens Audiologische Technik Gmbh | Hearing aid with a directional microphone system as well as method for the operation thereof |
US6198693B1 (en) * | 1998-04-13 | 2001-03-06 | Andrea Electronics Corporation | System and method for finding the direction of a wave source using an array of sensors |
US6385323B1 (en) * | 1998-05-15 | 2002-05-07 | Siemens Audiologische Technik Gmbh | Hearing aid with automatic microphone balancing and method for operating a hearing aid with automatic microphone balancing |
US6137889A (en) * | 1998-05-27 | 2000-10-24 | Insonus Medical, Inc. | Direct tympanic membrane excitation via vibrationally conductive assembly |
US6717991B1 (en) * | 1998-05-27 | 2004-04-06 | Telefonaktiebolaget Lm Ericsson (Publ) | System and method for dual microphone signal noise reduction using spectral subtraction |
US6217508B1 (en) * | 1998-08-14 | 2001-04-17 | Symphonix Devices, Inc. | Ultrasonic hearing system |
US6182018B1 (en) * | 1998-08-25 | 2001-01-30 | Ford Global Technologies, Inc. | Method and apparatus for identifying sound in a composite sound signal |
US6342035B1 (en) * | 1999-02-05 | 2002-01-29 | St. Croix Medical, Inc. | Hearing assistance device sensing otovibratory or otoacoustic emissions evoked by middle ear vibrations |
US6390971B1 (en) * | 1999-02-05 | 2002-05-21 | St. Croix Medical, Inc. | Method and apparatus for a programmable implantable hearing aid |
US6549586B2 (en) * | 1999-04-12 | 2003-04-15 | Telefonaktiebolaget L M Ericsson | System and method for dual microphone signal noise reduction using spectral subtraction |
US6421448B1 (en) * | 1999-04-26 | 2002-07-16 | Siemens Audiologische Technik Gmbh | Hearing aid with a directional microphone characteristic and method for producing same |
US6272229B1 (en) * | 1999-08-03 | 2001-08-07 | Topholm & Westermann Aps | Hearing aid with adaptive matching of microphones |
US6397186B1 (en) * | 1999-12-22 | 2002-05-28 | Ambush Interactive, Inc. | Hands-free, voice-operated remote control transmitter |
US20020029070A1 (en) * | 2000-04-13 | 2002-03-07 | Hans Leysieffer | At least partially implantable system for rehabilitation a hearing disorder |
US20010031996A1 (en) * | 2000-04-13 | 2001-10-18 | Hans Leysieffer | At least partially implantable system for rehabilitation of a hearing disorder |
US20020012438A1 (en) * | 2000-06-30 | 2002-01-31 | Hans Leysieffer | System for rehabilitation of a hearing disorder |
US20020019668A1 (en) * | 2000-08-11 | 2002-02-14 | Friedemann Stockert | At least partially implantable system for rehabilitation of a hearing disorder |
US20020110255A1 (en) * | 2000-10-05 | 2002-08-15 | Killion Mead C. | Directional microphone assembly |
US20020057817A1 (en) * | 2000-10-10 | 2002-05-16 | Resistance Technology, Inc. | Hearing aid |
US6380896B1 (en) * | 2000-10-30 | 2002-04-30 | Siemens Information And Communication Mobile, Llc | Circular polarization antenna for wireless communication system |
US20020141595A1 (en) * | 2001-02-23 | 2002-10-03 | Jouppi Norman P. | System and method for audio telepresence |
US20020131613A1 (en) * | 2001-03-13 | 2002-09-19 | Andreas Jakob | Method for establishing a binaural communication link and binaural hearing devices |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8638946B1 (en) | 2004-03-16 | 2014-01-28 | Genaudio, Inc. | Method and apparatus for creating spatialized sound |
US20100027821A1 (en) * | 2006-10-10 | 2010-02-04 | Eghart Fischer | Hearing aid and method for operating a hearing aid |
US8325957B2 (en) | 2006-10-10 | 2012-12-04 | Siemens Audiologische Technik Gmbh | Hearing aid and method for operating a hearing aid |
WO2008106680A3 (en) * | 2007-03-01 | 2008-10-16 | Jerry Mahabub | Audio spatialization and environment simulation |
US9197977B2 (en) | 2007-03-01 | 2015-11-24 | Genaudio, Inc. | Audio spatialization and environment simulation |
US9271080B2 (en) | 2007-03-01 | 2016-02-23 | Genaudio, Inc. | Audio spatialization and environment simulation |
WO2008106680A2 (en) * | 2007-03-01 | 2008-09-04 | Jerry Mahabub | Audio spatialization and environment simulation |
US20090046864A1 (en) * | 2007-03-01 | 2009-02-19 | Genaudio, Inc. | Audio spatialization and environment simulation |
EP2207168A3 (en) * | 2007-10-18 | 2010-10-20 | Motorola, Inc. | Robust two microphone noise suppression system |
EP2183853A4 (en) * | 2007-10-18 | 2010-11-03 | Motorola Inc | Robust two microphone noise suppression system |
US8046219B2 (en) | 2007-10-18 | 2011-10-25 | Motorola Mobility, Inc. | Robust two microphone noise suppression system |
US20090106021A1 (en) * | 2007-10-18 | 2009-04-23 | Motorola, Inc. | Robust two microphone noise suppression system |
EP2183853A1 (en) * | 2007-10-18 | 2010-05-12 | Motorola, Inc. | Robust two microphone noise suppression system |
WO2009051959A1 (en) | 2007-10-18 | 2009-04-23 | Motorola, Inc. | Robust two microphone noise suppression system |
KR101171494B1 (en) * | 2007-10-18 | 2012-08-07 | 모토로라 모빌리티, 인크. | Robust two microphone noise suppression system |
US20100246831A1 (en) * | 2008-10-20 | 2010-09-30 | Jerry Mahabub | Audio spatialization and environment simulation |
US8520873B2 (en) | 2008-10-20 | 2013-08-27 | Jerry Mahabub | Audio spatialization and environment simulation |
US8818800B2 (en) * | 2011-07-29 | 2014-08-26 | 2236008 Ontario Inc. | Off-axis audio suppressions in an automobile cabin |
US20140348333A1 (en) * | 2011-07-29 | 2014-11-27 | 2236008 Ontario Inc. | Off-axis audio suppressions in an automobile cabin |
US20130030801A1 (en) * | 2011-07-29 | 2013-01-31 | Qnx Software Systems Limited | Off-Axis Audio Suppressions in An Automobile Cabin |
US9437181B2 (en) * | 2011-07-29 | 2016-09-06 | 2236008 Ontario Inc. | Off-axis audio suppression in an automobile cabin |
EP2833651A1 (en) * | 2013-08-01 | 2015-02-04 | Siemens Medical Instruments Pte. Ltd. | Method for tracking of a sound source |
Also Published As
Publication number | Publication date |
---|---|
DE60125553D1 (en) | 2007-02-08 |
CA2407855A1 (en) | 2001-11-15 |
DK1312239T3 (en) | 2007-04-30 |
US7613309B2 (en) | 2009-11-03 |
WO2001087011A3 (en) | 2003-03-20 |
JP2003533152A (en) | 2003-11-05 |
AU2001261344A1 (en) | 2001-11-20 |
CA2407855C (en) | 2010-02-02 |
CN1440628A (en) | 2003-09-03 |
WO2001087011A2 (en) | 2001-11-15 |
CA2685434A1 (en) | 2001-11-15 |
DE60125553T2 (en) | 2007-10-04 |
EP1312239B1 (en) | 2006-12-27 |
EP1312239A2 (en) | 2003-05-21 |
US20030138116A1 (en) | 2003-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7613309B2 (en) | Interference suppression techniques | |
US7076072B2 (en) | Systems and methods for interference-suppression with directional sensing patterns | |
Lotter et al. | Dual-channel speech enhancement by superdirective beamforming | |
US9113247B2 (en) | Device and method for direction dependent spatial noise reduction | |
US8565446B1 (en) | Estimating direction of arrival from plural microphones | |
US8180067B2 (en) | System for selectively extracting components of an audio input signal | |
US8098844B2 (en) | Dual-microphone spatial noise suppression | |
Lockwood et al. | Performance of time-and frequency-domain binaural beamformers based on recorded signals from real rooms | |
US20140003635A1 (en) | Audio signal processing device calibration | |
US8213623B2 (en) | Method to generate an output audio signal from two or more input audio signals | |
US20080201138A1 (en) | Headset for Separation of Speech Signals in a Noisy Environment | |
US8615392B1 (en) | Systems and methods for producing an acoustic field having a target spatial pattern | |
JP2013543987A (en) | System, method, apparatus and computer readable medium for far-field multi-source tracking and separation | |
JP2001516196A (en) | Electronic beam forming method of acoustic signal and acoustic sensor device | |
Neo et al. | Robust microphone arrays using subband adaptive filters | |
US20130253923A1 (en) | Multichannel enhancement system for preserving spatial cues | |
EP1065909A2 (en) | Noise canceling microphone array | |
EP4161105A1 (en) | Spatial audio filtering within spatial audio capture | |
Farmani et al. | Sound source localization for hearing aid applications using wireless microphones | |
US11470429B2 (en) | Method of operating an ear level audio system and an ear level audio system | |
Chatlani et al. | Spatial noise reduction in binaural hearing aids | |
Zhang et al. | A frequency domain approach for speech enhancement with directionality using compact microphone array. | |
Nordholm et al. | Hands‐free mobile telephony by means of an adaptive microphone array | |
Zhang et al. | A compact-microphone-array-based speech enhancement algorithm using auditory subbands and probability constrained postfilter | |
Ristimäki | Hajautettu mikrofoniryhmäjärjestelmä kahdensuuntaisessa äänikommunikaatiossa |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |