AU2009308442A1 - Acoustic Voice Activity Detection (AVAD) for electronic systems - Google Patents
Acoustic Voice Activity Detection (AVAD) for electronic systems Download PDFInfo
- Publication number
- AU2009308442A1 AU2009308442A1 AU2009308442A AU2009308442A AU2009308442A1 AU 2009308442 A1 AU2009308442 A1 AU 2009308442A1 AU 2009308442 A AU2009308442 A AU 2009308442A AU 2009308442 A AU2009308442 A AU 2009308442A AU 2009308442 A1 AU2009308442 A1 AU 2009308442A1
- Authority
- AU
- Australia
- Prior art keywords
- microphone
- signal
- speech
- physical
- filter
- 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
- 238000001514 detection method Methods 0.000 title claims description 49
- 230000000694 effects Effects 0.000 title claims description 40
- 230000004044 response Effects 0.000 claims description 144
- 238000000034 method Methods 0.000 claims description 122
- 230000003044 adaptive effect Effects 0.000 claims description 52
- 238000012545 processing Methods 0.000 claims description 46
- 230000008569 process Effects 0.000 claims description 10
- 230000003068 static effect Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 description 46
- 238000010586 diagram Methods 0.000 description 24
- 230000001629 suppression Effects 0.000 description 21
- 230000000875 corresponding effect Effects 0.000 description 16
- 230000006870 function Effects 0.000 description 14
- 230000005534 acoustic noise Effects 0.000 description 13
- 230000008859 change Effects 0.000 description 11
- 238000010276 construction Methods 0.000 description 10
- 238000012549 training Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 7
- 230000033001 locomotion Effects 0.000 description 7
- 238000012546 transfer Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 230000006978 adaptation Effects 0.000 description 5
- 238000003491 array Methods 0.000 description 5
- 230000007423 decrease Effects 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 5
- 238000009472 formulation Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 230000002829 reductive effect Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 102100026436 Regulator of MON1-CCZ1 complex Human genes 0.000 description 3
- 230000001413 cellular effect Effects 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- 230000001934 delay Effects 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 101710180672 Regulator of MON1-CCZ1 complex Proteins 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 229920000547 conjugated polymer Polymers 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910044991 metal oxide Inorganic materials 0.000 description 2
- 150000004706 metal oxides Chemical class 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 102000008482 12E7 Antigen Human genes 0.000 description 1
- 108010020567 12E7 Antigen Proteins 0.000 description 1
- 101100494265 Caenorhabditis elegans best-15 gene Proteins 0.000 description 1
- 101100172892 Caenorhabditis elegans sec-8 gene Proteins 0.000 description 1
- 102000003712 Complement factor B Human genes 0.000 description 1
- 108090000056 Complement factor B Proteins 0.000 description 1
- 101000893549 Homo sapiens Growth/differentiation factor 15 Proteins 0.000 description 1
- 101000692878 Homo sapiens Regulator of MON1-CCZ1 complex Proteins 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001010 compromised effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011982 device technology Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- ZMRUPTIKESYGQW-UHFFFAOYSA-N propranolol hydrochloride Chemical compound [H+].[Cl-].C1=CC=C2C(OCC(O)CNC(C)C)=CC=CC2=C1 ZMRUPTIKESYGQW-UHFFFAOYSA-N 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 210000001260 vocal cord Anatomy 0.000 description 1
- 230000001755 vocal effect Effects 0.000 description 1
Classifications
-
- 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
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/78—Detection of presence or absence of voice signals
-
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Circuit For Audible Band Transducer (AREA)
Description
WO 2010/048635 PCT/US2009/062129 ACOUSTIC VOICE ACTIVITY DETECTION (AVAD) FOR ELECTRONIC SYSTEMS 5 RELATED APPLICATIONS This application claims the benefit of United States (US) Patent 10 Application number 61/108,426, filed October 24, 2008. This application is a continuation in part of US Patent Application number 11/805,987, filed May 25, 2007. This application is a continuation in part of US Patent Application number 12/139,333, filed June 13, 2008. 15 TECHNICAL FIELD The disclosure herein relates generally to noise suppression. In particular, this disclosure relates to noise suppression systems, devices, and methods for use in acoustic applications. 20 BACKGROUND The ability to correctly identify voiced and unvoiced speech is critical to many speech applications including speech recognition, speaker verification, noise suppression, and many others. In a typical acoustic 25 application, speech from a human speaker is captured and transmitted to a receiver in a different location. In the speaker's environment there may exist one or more noise sources that pollute the speech signal, the signal of interest, with unwanted acoustic noise. This makes it difficult or impossible for the receiver, whether human or machine, to understand the user's 30 speech. Typical methods for classifying voiced and unvoiced speech have relied mainly on the acoustic content of single microphone data, which is 1 WO 2010/048635 PCT/US2009/062129 plagued by problems with noise and the corresponding uncertainties in signal content. This is especially problematic with the proliferation of portable communication devices like mobile telephones. There are methods known in the art for suppressing the noise present in the speech 5 signals, but these normally require a robust method of determining when speech is being produced. Non-acoustic methods have been employed successfully in commercial products such as the Jawbone headset produced by Aliphcom, Inc., San Francisco, California (Aliph), but an acoustic-only solution is desired in some cases (e.g., for reduced cost, as a supplement 10 to the non-acoustic sensor, etc.). INCORPORATION BY REFERENCE Each patent, patent application, and/or publication mentioned in this specification is herein incorporated by reference in its entirety to the same 15 extent as if each individual patent, patent application, and/or publication was specifically and individually indicated to be incorporated by reference. BRIEF DESCRIPTION OF THE DRAWINGS Figure 1 is a configuration of a two-microphone array with speech 20 source S, under an embodiment. Figure 2 is a block diagram of V 2 construction using a fixed p(z), under an embodiment. Figure 3 is a block diagram of V 2 construction using an adaptive p(z), under an embodiment. 25 Figure 4 is a block diagram of V 1 construction, under an embodiment. Figure 5 is a flow diagram of acoustic voice activity detection, under an embodiment. Figure 6 shows experimental results of the algorithm using a fixed 30 beta when only noise is present, under an embodiment. Figure 7 shows experimental results of the algorithm using a fixed beta when only speech is present, under an embodiment. 2 WO 2010/048635 PCT/US2009/062129 Figure 8 shows experimental results of the algorithm using a fixed beta when speech and noise is present, under an embodiment. Figure 9 shows experimental results of the algorithm using an adaptive beta when only noise is present, under an embodiment. 5 Figure 10 shows experimental results of the algorithm using an adaptive beta when only speech is present, under an embodiment. Figure 11 shows experimental results of the algorithm using an adaptive beta when speech and noise is present, under an embodiment. Figure 12 is a block diagram of a NAVSAD system, under an embodiment. 10 Figure 13 is a block diagram of a PSAD system, under an embodiment, Figure 14 is a block diagram of a denoising subsystem, referred to herein as the Pathfinder system, under an embodiment, Figure 15 is a flow diagram of a detection algorithm for use in detecting voiced and unvoiced speech, under an embodiment. 15 Figures 16A, 16B, and 17 show data plots for an example in which a subject twice speaks the phrase "pop pan", under an embodiment, Figure 16A plots the received GEMS signal for this utterance along with the mean correlation between the GEMS signal and the Mic 1 signal and the threshold T1 used for voiced speech detection, under an embodiment. 20 Figure 16B plots the received GEMS signal for this utterance along with the standard deviation of the GEMS signal and the threshold T2 used for voiced speech detection, under an embodiment. Figure 17 plots voiced speech detected from the acoustic or audio signal, along with the GEMS signal and the acoustic noise; no unvoiced speech is detected in 25 this example because of the heavy background babble noise, under an embodiment. Figure 18 is a microphone array for use under an embodiment of the PSAD system. Figure 19 is a plot of AM versus d 1 for several Ad values, under an embodiment, 30 Figure 20 shows a plot of the gain parameter as the sum of the absolute values of H 1 (z) and the acoustic data or audio from microphone 1, under an embodiment 3 WO 2010/048635 PCT/US2009/062129 Figure 21 is an alternative plot of acoustic data presented in Figure 20, under an embodiment., Figure 22 is a two-microphone adaptive noise suppression system, under an embodiment. 5 Figure 23 is a generalized two-microphone array (DOMA) including an array and speech source S configuration, under an embodiment. Figure 24 is a system for generating or producing a first order gradient microphone V using two omnidirectional elements 01 and 02, under an embodiment. Figure 25 is a block diagram for a DOMA including two physical microphones 10 configured to form two virtual microphones V 1 and V 2 , under an embodiment, Figure 26 is a block diagram for a DOMA including two physical microphones configured to form N virtual microphones V 1 through VN, where N is any number greater than one, under an embodiment. Figure 27 is an example of a headset or head-worn device that includes the 15 DOMA, as described herein, under an embodiment. Figure 28 is a flow diagram for denoising acoustic signals using the DOMA, under an embodiment. Figure 29 is a flow diagram for forming the DOMA, under an embodiment. Figure 30 is a plot of linear response of virtual microphone V 2 with P = 0.8 to a 20 1 kHz speech source at a distance of 0.1 m, under an embodiment. Figure 31 is a plot of linear response of virtual microphone V 2 with P = 0.8 to a 1 kHz noise source at a distance of 1.0 m, under an embodiment. Figure 32 is a plot of linear response of virtual microphone V 1 with P = 0 8 to a 1 kHz speech source at a distance of 0.1 m, under an embodiment 25 Figure 33 is a plot of linear response of virtual microphone V 1 with P = 0.8 to a 1 kHz noise source at a distance of 1.0 m, under an embodiment. Figure 34 is a plot of linear response of virtual microphone V 1 with P = 0.8 to a speech source at a distance of 0 1 m for frequencies of 100, 500, 1000, 2000, 3000, and 4000 Hz, under an embodiment, 30 Figure 35 is a plot showing comparison of frequency responses for speech for the array of an embodiment and for a conventional cardioid microphone, under an embodiment. 4 WO 2010/048635 PCT/US2009/062129 Figure 36 is a plot showing speech response for V 1 (top, dashed) and V 2 (bottom, solid) versus B with d, assumed to be 0.1 m, under an embodiment, under an embodiment. Figure 37 is a plot showing a ratio of V 1
/V
2 speech responses shown in Figure 5 31 versus B, under an embodiment Figure 38 is a plot of B versus actual d, assuming that ds = 10 cm and theta = 0, under an embodiment. Figure 39 is a plot of B versus theta with ds = 10 cm and assuming ds = 10 cm, under an embodiment. 10 Figure 40 is a plot of amplitude (top) and phase (bottom) response of N(s) with B = 1 and D = -7.2 pLsec, under an embodiment, Figure 41 is a plot of amplitude (top) and phase (bottom) response of N(s) with B = 1.2 and D = -7.2 pLsec, under an embodiment. Figure 42 is a plot of amplitude (top) and phase (bottom) response of the effect 15 on the speech cancellation in V 2 due to a mistake in the location of the speech source with q1 = 0 degrees and q2 = 30 degrees, under an embodiment. Figure 43 is a plot of amplitude (top) and phase (bottom) response of the effect on the speech cancellation in V 2 due to a mistake in the location of the speech source with q1 = 0 degrees and q2 = 45 degrees, under an embodiment. 20 Figure 44 shows experimental results for a 2do = 19 mm array using a linear p of 0.83 and B1 = B2 = 1 on a Bruel and Kjaer Head and Torso Simulator (HATS) in very loud (-85 dBA) music/speech noise environment. DETAILED DESCRIPTION 25 Acoustic Voice Activity Detection (AVAD) methods and systems are described herein. The AVAD methods and systems, which include algorithms or programs, use microphones to generate virtual directional microphones which have very similar noise responses and very dissimilar speech responses. The ratio of the energies of the virtual microphones is 30 then calculated over a given window size and the ratio can then be used with a variety of methods to generate a VAD signal. The virtual microphones can be constructed using either a fixed or an adaptive filter. The adaptive filter generally results in a more accurate and noise-robust 5 WO 2010/048635 PCT/US2009/062129 VAD signal but requires training. In addition, restrictions can be placed on the filter to ensure that it is training only on speech and not on environmental noise. In the following description, numerous specific details are introduced 5 to provide a thorough understanding of, and enabling description for, embodiments. One skilled in the relevant art, however, will recognize that these embodiments can be practiced without one or more of the specific details, or with other components, systems, etc. In other instances, well-known structures or operations are not shown, or are not described in 10 detail, to avoid obscuring aspects of the disclosed embodiments. Figure 1 is a configuration of a two-microphone array of the AVAD with speech source S, under an embodiment. The AVAD of an embodiment uses two physical microphones (01 and 02) to form two virtual microphones (V 1 and V 2 ). The virtual microphones of an embodiment are 15 directional microphones, but the embodiment is not so limited. The physical microphones of an embodiment include omnidirectional microphones, but the embodiments described herein are not limited to omnidirectional microphones. The virtual microphone (VM) V 2 is configured in such a way that it has minimal response to the speech of the user, while 20 Vi is configured so that it does respond to the user's speech but has a very similar noise magnitude response to V 2 , as described in detail herein. The PSAD VAD methods can then be used to determine when speech is taking place. A further refinement is the use of an adaptive filter to further minimize the speech response of V 2 , thereby increasing the speech energy 25 ratio used in PSAD and resulting in better overall performance of the AVAD. The PSAD algorithm as described herein calculates the ratio of the energies of two directional microphones M 1 and M 2 : 30 6 WO 2010/048635 PCT/US2009/062129 where the "z" indicates the discrete frequency domain and "i" ranges from the beginning of the window of interest to the end, but the same relationship holds in the time domain. The summation can occur over a window of any length; 200 samples at a sampling rate of 8 kHz has been 5 used to good effect. Microphone M 1 is assumed to have a greater speech response than microphone M 2 . The ratio R depends on the relative strength of the acoustic signal of interest as detected by the microphones. For matched omnidirectional microphones (i.e. they have the same response to acoustic signals for all spatial orientations and frequencies), 10 the size of R can be calculated for speech and noise by approximating the propagation of speech and noise waves as spherically symmetric sources. For these the energy of the propagating wave decreases as 16: 15 The distance di is the distance from the acoustic source to M 1 , d 2 is the distance from the acoustic source to M 2 , and d = d 2 -di (see Figure 1). It is assumed that 01 is closer to the speech source (the user's mouth) so that d is always positive. If the microphones and the user's mouth are all 20 on a line, then d = 2do, the distance between the microphones. For matched omnidirectional microphones, the magnitude of R, depends only on the relative distance between the microphones and the acoustic source. For noise sources, the distances are typically a meter or more, and for speech sources, the distances are on the order of 10 cm, but the distances 25 are not so limited. Therefore for a 2-cm array typical values of R are: .~ 12 cm~ d 1 10 cm d, 102 cm R, =-~ t=1.02 di 100 cm 7 WO 2010/048635 PCT/US2009/062129 where the "S" subscript denotes the ratio for speech sources and "N" the ratio for noise sources. There is not a significant amount of separation between noise and speech sources in this case, and therefore it would be 5 difficult to implement a robust solution using simple omnidirectional microphones. A better implementation is to use directional microphones where the second microphone has minimal speech response. As described herein, such microphones can be constructed using omnidirectional microphones 10 01 and 02: where a(z) is a calibration filter used to compensate 02's response so that 15 it is the same as 01, p(z) is a filter that describes the relationship between 01 and calibrated 02 for speech, and y is a fixed delay that depends on the size of the array. There is no loss of generality in defining a(z) as above, as either microphone may be compensated to match the other. For this configuration V 1 and V 2 have very similar noise response magnitudes and 20 very dissimilar speech response magnitudes if where again d = 2do and c is the speed of sound in air, which is temperature dependent and approximately 25 1 c = 31.3 .1 -___ 273-15 sec 8 WO 2010/048635 PCT/US2009/062129 where T is the temperature of the air in Celsius. The filter p(z) can be calculated using wave theory to be - [21 5 where again dk is the distance from the user's mouth to Ok. Figure 2 is a block diagram of V 2 construction using a fixed p(z), under an embodiment. This fixed (or static) p works sufficiently well if the calibration filter a(z) is accurate and di and d 2 are accurate for the user. This fixed-p algorithm, 10 however, neglects important effects such as reflection, diffraction, poor array orientation (i.e. the microphones and the mouth of the user are not all on a line), and the possibility of different di and d 2 values for different users. The filter P(z) can also be determined experimentally using an 15 adaptive filter. Figure 3 is a block diagram of V 2 construction using an adaptive p(z), under an embodiment, where: 20 The adaptive process varies ) to minimize the output of V 2 when only speech is being received by 01 and 02. A small amount of noise may be tolerated with little ill effect, but it is preferred that only speech is being received when the coefficients of e.() are calculated. Any adaptive process may be used; a normalized least-mean squares (NLMS) algorithm was 25 used in the examples below. The V 1 can be constructed using the current value for j() or the fixed filter .fz) can be used for simplicity. Figure 4 is a block diagram of
V
1 construction, under an embodiment. 9 WO 2010/048635 PCT/US2009/062129 Now the ratio R is where double bar indicates norm and again any size window may be used. If j -) has been accurately calculated, the ratio for speech should be 5 relatively high (e.g., greater than approximately 2) and the ratio for noise should be relatively low (e.g., less than approximately 1.1). The ratio calculated will depend on both the relative energies of the speech and noise as well as the orientation of the noise and the reverberance of the environment. In practice, either the adapted filter (z) or the static filter 10 b(z) may be used for V 1 (z) with little effect on R - but it is important to use the adapted filter Iz) in V 2 (z) for best performance. Many techniques known to those skilled in the art (e.g., smoothing, etc.) can be used to make R more amenable to use in generating a VAD and the embodiments herein are not so limited. 15 The ratio R can be calculated for the entire frequency band of interest, or can be calculated in frequency subbands. One effective subband discovered was 250 Hz to 1250 Hz, another was 200 Hz to 3000 Hz, but many others are possible and useful. Once generated, the vector of the ratio R versus time (or the matrix 20 of R versus time if multiple subbands are used) can be used with any detection system (such as one that uses fixed and/or adaptive thresholds) to determine when speech is occurring. While many detection systems and methods are known to exist by those skilled in the art and may be used, the method described herein for generating an R so that the speech is 25 easily discernable is novel. It is important to note that the R does not depend on the type of noise or its orientation or frequency content; R simply depends on the V 1 and V 2 spatial response similarity for noise and spatial response dissimilarity for speech. In this way it is very robust and can operate smoothly in a variety of noisy acoustic environments. 10 WO 2010/048635 PCT/US2009/062129 Figure 5 is a flow diagram of acoustic voice activity detection 500, under an embodiment. The detection comprises forming a first virtual microphone by combining a first signal of a first physical microphone and a second signal of a second physical microphone 502. The detection 5 comprises forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone 504. The detection comprises forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal, and summing the first intermediate signal and the second signal 506. The 10 detection comprises generating an energy ratio of energies of the first virtual microphone and the second virtual microphone 508. The detection comprises detecting acoustic voice activity of a speaker when the energy ratio is greater than a threshold value 510. The accuracy of the adaptation to the p(z) of the system is a factor in 15 determining the effectiveness of the AVAD. A more accurate adaptation to the actual p(z) of the system leads to lower energy of the speech response in V 2 , and a higher ratio R. The noise (far-field) magnitude response is largely unchanged by the adaptation process, so the ratio R will be near unity for accurately adapted beta. For purposes of accuracy, the system 20 can be trained on speech alone, or the noise should be low enough in energy so as not to affect or to have a minimal affect the training. To make the training as accurate as possible, the coefficients of the filter p(z) of an embodiment are generally updated under the following conditions, but the embodiment is not so limited: speech is being produced 25 (requires a relatively high SNR or other method of detection such as an Aliph Skin Surface Microphone (SSM) as described in United States Patent Application number 10/769,302, filed January 30, 2004, which is incorporated by reference herein in its entirety); no wind is detected (wind can be detected using many different methods known in the art, such as 30 examining the microphones for uncorrelated low-frequency noise); and the current value of R is much larger than a smoothed history of R values (this ensures that training occurs only when strong speech is present). These 11 WO 2010/048635 PCT/US2009/062129 procedures are flexible and others may be used without significantly affecting the performance of the system. These restrictions can make the system relatively more robust. Even with these precautions, it is possible that the system 5 accidentally trains on noise (e.g., there may be a higher likelihood of this without use of a non-acoustic VAD device such as the SSM used in the Jawbone headset produced by Aliph, San Francisco, California). Thus, an embodiment includes a further failsafe system to preclude accidental training from significantly disrupting the system. The adaptive P is limited 10 to certain values expected for speech. For example, values for di for an ear-mounted headset will normally fall between 9 and 14 centimeters, so using an array length of 2do =2.0 cm and Equation 2 above, ad d CI which means that 15 0,82 02 8. The magnitude of the P filter can therefore be limited to between approximately 0.82 and 0.88 to preclude problems if noise is present during training. Looser limits can be used to compensate for inaccurate 20 calibrations (the response of omnidirectional microphones is usually calibrated to one another so that their frequency response LIs I sam to the same acoustic source - if the calibration is not completely accurate the virtual microphones may not form properly). Similarly, the phase of the P filter can be limited to be what is 25 expected from a speech source within +- 30 degrees from the axis of the array. As described herein, and with reference to Figure 1, d: - di (seconds) C 12 WO 2010/048635 PCT/US2009/062129 d, 2 d'.7 arCC S C d1 = d - 2dsde'~ cocsti -d where ds is the distance from the midpoint of the array to the speech 5 source. Varying ds from 10 to 15 cm and allowing 0 to vary between 0 and +- 30 degrees, the maximum difference in y results from the difference of y at 0 degrees (58.8 psec) and y at +-30 degrees for ds = 10 cm (50.8 [tsec). This means that the maximum expected phase difference is 58.8 50.8 = 8.0 [tsec, or 0.064 samples at an 8 kHz sampling rate. Since 10 = 2,ft- = fi3Ox10 -) ra d the maximum phase difference realized at 4 kHz is only 0.2 rad or about 11.4 degrees, a small amount, but not a negligible one. Therefore the P filter should almost linear phase, but some allowance made for differences 15 in position and angle. In practice a slightly larger amount was used (0.071 samples at 8 kHz) in order to compensate for poor calibration and diffraction effects, and this worked well. The limit on the phase in the example below was implemented as the ratio of the central tap energy to the combined energy of the other taps: ph ase limit ratio 20 where P is the current estimate. This limits the phase by restricting the effects of the non-center taps. Other ways of limiting the phase of the beta filter are known to those skilled in the art and the algorithm presented here is not so limited. 25 Embodiments are presented herein that use both a fixed p(z) and an adaptive p(z), as described in detail above. In both cases, R was calculated using frequencies between 250 and 3000 Hz using a window size of 200 samples at 8 kHz. The results for V 1 (top plot), V 2 (middle plot), R (bottom plot, 13 WO 2010/048635 PCT/US2009/062129 solid line, windowed using a 200 sample rectangular window at 8 kHz) and the VAD (bottom plot, dashed line) are shown in Figures 6-11. Figures 6-11 demonstrate the use of a fixed beta filter p(z) in conditions of only noise (street and bus noise, approximately 70 dB SPL at the ear), only 5 speech (normalized to 94 dB SPL at the mouth reference point (MRP)), and mixed noise and speech, respectively. A Bruel & Kjaer Head and Torso Simulator (HATS) was used for the tests and the omnidirectional microphones mounted on HATS' ear with the midline of the array approximately 11 cm from the MRP. The fixed beta filter used was 10 pF(Z)=0.
8 2 , where the "F" subscript indicates a fixed filter. The VAD was calculated using a fixed threshold of 1.5. Figure 6 shows experimental results of the algorithm using a fixed beta when only noise is present, under an embodiment. The top plot is V 1 , the middle plot is V 2 , and the bottom plot is R (solid line) and the VAD 15 result (dashed line) versus time. Examining Figure 6, the response of both V 1 and V 2 are very similar, and the ratio R is very near unity for the entire sample. The VAD response has occasional false positives denoted by spikes in the R plot (windows that are identified by the algorithm as containing speech when they do not), but these are easily removed using 20 standard pulse removal algorithms and/or smoothing of the R results. Figure 7 shows experimental results of the algorithm using a fixed beta when only speech is present, under an embodiment. The top plot is
V
1 , the middle plot is V 2 , and the bottom plot is R (solid line) and the VAD result (dashed line) versus time. The R ratio is between approximately 2 25 and approximately 7 on average, and the speech is easily discernable using the fixed threshold. These results show that the response of the two virtual microphones to speech are very different, and indeed that ratio R varies from 2-7 during speech. There are very few false positives and very few false negatives (windows that contain speech but are not identified as 30 speech windows). The speech is easily and accurately detected. Figure 8 shows experimental results of the algorithm using a fixed beta when speech and noise is present, under an embodiment. The top 14 WO 2010/048635 PCT/US2009/062129 plot is V 1 , the middle plot is V 2 , and the bottom plot is R (solid line) and the VAD result (dashed line) versus time. The R ratio is lower than when no noise is present, but the VAD remains accurate with only a few false positives. There are more false negatives than with no noise, but the 5 speech remains easily detectable using standard thresholding algorithms. Even in a moderately loud noise environment (Figure 8) the R ratio remains significantly above unity, and the VAD once again returns few false positives. More false negatives are observed, but these may be reduced using standard methods such as smoothing of R and allowing the VAD to 10 continue reporting voiced windows for a few windows after R is under the threshold. Results using the adaptive beta filter are shown in Figures 9-11. The adaptive filter used was a five-tap NLMS FIR filter using the frequency band from 100 Hz to 3500 Hz. A fixed filter of z- 0 4 3 is used to filter 01 so 15 that 01 and 02 are aligned for speech before the adaptive filter is calculated. The adaptive filter was constrained using the methods above using a low p limit of 0.73, a high p limit of 0.98, and a phase limit ratio of 0.98. Again a fixed threshold was used to generate the VAD result from the ratio R, but in this case a threshold value of 2.5 was used since the R 20 values using the adaptive beta filter are normally greater than when the fixed filter is used. This allows for a reduction of false positives without significantly increasing false negatives. Figure 9 shows experimental results of the algorithm using an adaptive beta when only noise is present, under an embodiment. The top 25 plot is V 1 , the middle plot is V 2 , and the bottom plot is R (solid line) and the VAD result (dashed line) versus time, with the y-axis expanded to 0 50. Again, V 1 and V 2 are very close in energy and the R ratio is near unity. Only a single false positive was generated. Figure 10 shows experimental results of the algorithm using an 30 adaptive beta when only speech is present, under an embodiment. The top plot is V 1 , the middle plot is V 2 , and the bottom plot is (solid line) and the VAD result (dashed line) versus time, expanded to 0-50. The V 2 15 WO 2010/048635 PCT/US2009/062129 response is greatly reduced using the adaptive beta, and the R ratio has increased from the range of approximately 2-7 to the range of approximately 5-30 on average, making the speech even simpler to detect using standard thresholding algorithms. There are almost no false 5 positives or false negatives. Therefore, the response of V 2 to speech is minimal, R is very high, and all of the speech is easily detected with almost no false positives. Figure 11 shows experimental results of the algorithm using an adaptive beta when speech and noise is present, under an embodiment. The top plot is V 1 , 10 the middle plot is V 2 , and the bottom plot is R (solid line) and the VAD result (dashed line) versus time, with the y-axis expanded to 0-50. The R ratio is again lower than when no noise is present, but this R with significant noise present results in a VAD signal that is about the same as the case using the fixed beta with no noise present. This shows that use of the adaptive beta allows the 15 system to perform well in higher noise environments than the fixed beta. Therefore, with mixed noise and speech, there are again very few false positives and fewer false negatives than in the results of Figure 8, demonstrating that the adaptive filter can outperform the fixed filter in the same noise environment. In practice, the adaptive filter has proven to be significantly more sensitive to speech 20 and less sensitive to noise. Detecting voiced and unvoiced speech using both acoustic and nonacoustic sensors Systems and methods for discriminating voiced and unvoiced speech from bacgund-I noise aren noA recrribed including a Pathfindeir Se A i election MAII~'' - A -. ''- __ - I (III~kV. I. Izi - 1A -1. Speech rtivifvt~ rin 25 (PSAD) system, referenced above, and a Non-Acoustic Sensor Voiced Speech Activity Detection (NAVSAD) system. The noise removal and reduction methods provided herein, while allowing for the separation and classification of unvoiced and voiced human speech from background noise, address the shortcomings of typical systems known in the art by cleaning acoustic signals of interest without distortion. 30 Figure 12 is a block diagram of a NAVSAD system 1200, under an embodiment. The NAVSAD system couples microphones 1210 and sensors 1220 to at least one processor 1230. The sensors1220 of an embodiment include voicing activity detectors or non-acoustic sensors, The processor 1230 controls subsystems including a detection subsystem 1250, referred to herein as a detection algorithm, and 16 WO 2010/048635 PCT/US2009/062129 a denoising subsystem 1240, Operation of the denoising subsystem 1240 is described in detail in the Related Applications. The NAVSAD system works extremely well in any background acoustic noise environment Figure 13 is a block diagram of a PSAD system 1300, under an embodiment 5 The PSAD system couples microphones 1210 to at least one processor 1230 The processor 1230 includes a detection subsystem 1250, referred to herein as a detection algorithm, and a denoising subsystem 1240 The PSAD system is highly sensitive in low acoustic noise environments and relatively insensitive in high acoustic noise environments. The PSAD can operate independently or as a backup to the 10 NAVSAD, detecting voiced speech if the NAVSAD fails. Note that the detection subsystems 1250 and denoising subsystems 1240 of both the NAVSAD and PSAD systems of an embodiment are algorithms controlled by the processor 1230, but are not so limited. Alternative embodiments of the NAVSAD and PSAD systems can include detection subsystems 1250 and/or denoising 15 subsystems 1240 that comprise additional hardware, firmware, software, and/or combinations of hardware, firmware, and software. Furthermore, functions of the detection subsystems 1250 and denoising subsystems 1240 may be distributed across numerous components of the NAVSAD and PSAD systems. Figure 14 is a block diagram of a denoising subsystem 1400, referred to herein 20 as the Pathfinder system, under an embodiment. The Pathfinder system is briefly described below, and is described in detail in the Related Applications. Two microphones Mic 1 and Mic 2 are used in the Pathfinder system, and Mic 1 is considered the "signal" microphone, With reference to Figure 12, the Pathfinder system 1400 is equivalent to the NAVSAD system 1200 when the voicing activity 25 detector (VAD) 1420 is a non-acoustic voicing sensor 1220 and the noise removal subsystem 1440 includes the detection subsystem 1250 and the denoising subsystem 1240 With reference to Figure 13, the Pathfinder system 1400 is equivalent to the PSAD system 1300 in the absence of the VAD 1420, and when the noise removal subsystem 1440 includes the detection subsystem 1250 and the denoising subsystem 30 1240. The NAVSAD and PSAD systems support a two-level commercial approach in which (i) a relatively less expensive PSAD system supports an acoustic approach that functions in most low- to medium-noise environments, and (ii) a NAVSAD system adds a non-acoustic sensor to enable detection of voiced speech in any environment. 17 WO 2010/048635 PCT/US2009/062129 Unvoiced speech is normally not detected using the sensor, as it normally does not sufficiently vibrate human tissue, However, in high noise situations detecting the unvoiced speech is not as important, as it is normally very low in energy and easily washed out by the noise, Therefore in high noise environments the unvoiced speech 5 is unlikely to affect the voiced speech denoising. Unvoiced speech information is most important in the presence of little to no noise and, therefore, the unvoiced detection should be highly sensitive in low noise situations, and insensitive in high noise situations. This is not easily accomplished, and comparable acoustic unvoiced detectors known in the art are incapable of operating under these environmental 10 constraints. The NAVSAD and PSAD systems include an array algorithm for speech detection that uses the difference in frequency content between two microphones to calculate a relationship between the signals of the two microphones. This is in contrast to conventional arrays that attempt to use the time/phase difference of each 15 microphone to remove the noise outside of an "area of sensitivity". The methods described herein provide a significant advantage, as they do not require a specific orientation of the array with respect to the signal. Further, the systems described herein are sensitive to noise of every type and every orientation, unlike conventional arrays that depend on specific noise 20 orientations,, Consequently, the frequency-based arrays presented herein are unique as they depend only on the relative orientation of the two microphones themselves with no dependence on the orientation of the noise and signal with respect to the microphones. This results in a robust signal processing system with respect to the type of noise, microphones, and orientation between the noise/signal source and the 25 microphones. The systems described herein use the information derived from the Pathfinder noise suppression system and/or a non-acoustic sensor described in the Related Applications to determine the voicing state of an input signal, as described in detail below The voicing state includes silent, voiced, and unvoiced states. The NAVSAD 30 system, for example, includes a non-acoustic sensor to detect the vibration of human tissue associated with speech. The non-acoustic sensor of an embodiment is a General Electromagnetic Movement Sensor (GEMS) as described briefly below and in detail in the Related Applications, but is not so limited Alternative embodiments, 18 WO 2010/048635 PCT/US2009/062129 however, may use any sensor that is able to detect human tissue motion associated with speech and is unaffected by environmental acoustic noise The GEMS is a radio frequency device (2.4 GHz) that allows the detection of moving human tissue dielectric interfaces. The GEMS includes an RF interferometer 5 that uses homodyne mixing to detect small phase shifts associated with target motion. In essence, the sensor sends out weak electromagnetic waves (less than 1 milliwatt) that reflect off of whatever is around the sensor. The reflected waves are mixed with the original transmitted waves and the results analyzed for any change in position of the targets. Anything that moves near the sensor will cause a change in phase of the 10 reflected wave that will be amplified and displayed as a change in voltage output from the sensor, A similar sensor is described by Gregory C. Burnett (1999) in "The physiological basis of glottal electromagnetic micropower sensors (GEMS) and their use in defining an excitation function for the human vocal tract"; Ph.D Thesis, University of California at Davis, 15 Figure 15 is a flow diagram of a detection algorithm 1250 for use in detecting voiced and unvoiced speech, under an embodiment. With reference to Figures 12 and 13, both the NAVSAD and PSAD systems of an embodiment include the detection algorithm 1250 as the detection subsystem 1250, This detection algorithm 1250 operates in real-time and, in an embodiment, operates on 20 millisecond 20 windows and steps 10 milliseconds at a time, but is not so limited. The voice activity determination is recorded for the first 10 milliseconds, and the second 10 milliseconds functions as a "look-ahead" buffer. While an embodiment uses the 20/10 windows, alternative embodiments may use numerous other combinations of window values. (onsideratinn Was given to a number of multi-dimensional factors in developing 25 the detection algorithm1250. The biggest consideration was to maintaining the effectiveness of the Pathfinder denoising technique, described in detail in the Related Applications and reviewed herein, Pathfinder performance can be compromised if the adaptive filter training is conducted on speech rather than on noise. It is therefore important not to exclude any significant amount of speech from the VAD to keep such 30 disturbances to a minimum, Consideration was also given to the accuracy of the characterization between voiced and unvoiced speech signals, and distinguishing each of these speech signals from noise signals, This type of characterization can be useful in such applications as speech recognition and speaker verification 19 WO 2010/048635 PCT/US2009/062129 Furthermore, the systems using the detection algorithm of an embodiment function in environments containing varying amounts of background acoustic noise. If the non-acoustic sensor is available, this external noise is not a problem for voiced speech. However, for unvoiced speech (and voiced if the non-acoustic sensor is not 5 available or has malfunctioned) reliance is placed on acoustic data alone to separate noise from unvoiced speech. An advantage inheres in the use of two microphones in an embodiment of the Pathfinder noise suppression system, and the spatial relationship between the microphones is exploited to assist in the detection of unvoiced speech. However, there may occasionally be noise levels high enough that 10 the speech will be nearly undetectable and the acoustic-only method will fail In these situations, the non-acoustic sensor (or hereafter just the sensor) will be required to ensure good performance. In the two-microphone system, the speech source should be relatively louder in one designated microphone when compared to the other microphone, Tests have 15 shown that this requirement is easily met with conventional microphones when the microphones are placed on the head, as any noise should result in an H 1 with a gain near unity.. Regarding the NAVSAD system, and with reference to Figure 12 and Figure 14, the NAVSAD relies on two parameters to detect voiced speech. These two 20 parameters include the energy of the sensor in the window of interest, determined in an embodiment by the standard deviation (SD), and optionally the cross-correlation (XCORR) between the acoustic signal from microphone I and the sensor data. The energy of the sensor can be determined in any one of a number of ways, and the SD is just one convenient way to determine the energy. 25 For the sensor, the SD is akin to the energy of the signal, which normally corresponds quite accurately to the voicing state, but may be susceptible to movement noise (relative motion of the sensor with respect to the human user) and/or electromagnetic noise. To further differentiate sensor noise from tissue motion, the XCORR can be used. The XCORR is only calculated to 15 delays, which 30 corresponds to just under 2 milliseconds at 8000 Hz. The XCORR can also be useful when the sensor signal is distorted or modulated in some fashion. For example, there are sensor locations (such as the jaw or back of the neck) where speech production can be detected but where the signal may have incorrect or distorted time-based information That is, they may not have 20 WO 2010/048635 PCT/US2009/062129 well defined features in time that will match with the acoustic waveform However, XCORR is more susceptible to errors from acoustic noise, and in high (<0 dB SNR) environments is almost useless, Therefore it should not be the sole source of voicing information 5 The sensor detects human tissue motion associated with the closure of the vocal folds, so the acoustic signal produced by the closure of the folds is highly correlated with the closures, Therefore, sensor data that correlates highly with the acoustic signal is declared as speech, and sensor data that does not correlate well is termed noise. The acoustic data is expected to lag behind the sensor data by about 10 0.1 to 0.8 milliseconds (or about 1-7 samples) as a result of the delay time due to the relatively slower speed of sound (around 330 m/s) However, an embodiment uses a 15-sample correlation, as the acoustic wave shape varies significantly depending on the sound produced, and a larger correlation width is needed to ensure detection. The SD and XCORR signals are related, but are sufficiently different so that the 15 voiced speech detection is more reliable. For simplicity, though, either parameter may be used. The values for the SD and XCORR are compared to empirical thresholds, and if both are above their threshold, voiced speech is declared. Example data is presented and described below. Figures 16A, 16B, and 17 show data plots for an example in which a subject 20 twice speaks the phrase "pop pan", under an embodiment. Figure 16A plots the received GEMS signal 1602 for this utterance along with the mean correlation 1604 between the GEMS signal and the Mic 1 signal and the threshold T1 used for voiced speech detection. Figure 16B plots the received GEMS signal 1602 for this utterance along With the standard deviation 1606 of the GEMS Sinal andl the threshold T2 used 25 for voiced speech detection. Figure 17 plots voiced speech 1702 detected from the acoustic or audio signal 1708, along with the GEMS signal 1704 and the acoustic noise 1706; no unvoiced speech is detected in this example because of the heavy background babble noise 1706 The thresholds have been set so that there are virtually no false negatives, and only occasional false positives A voiced speech 30 activity detection accuracy of greater than 99% has been attained under any acoustic background noise conditions,, The NAVSAD can determine when voiced speech is occurring with high degrees of accuracy due to the non-acoustic sensor data. However, the sensor offers little assistance in separating unvoiced speech from noise, as unvoiced speech 21 WO 2010/048635 PCT/US2009/062129 normally causes no detectable signal in most non-acoustic sensors If there is a detectable signal, the NAVSAD can be used, although use of the SD method is dictated as unvoiced speech is normally poorly correlated, In the absence of a detectable signal use is made of the system and methods of the Pathfinder noise 5 removal algorithm in determining when unvoiced speech is occurring. A brief review of the Pathfinder algorithm is described below, while a detailed description is provided in the Related Applications. With reference to Figure 14, the acoustic information coming into Microphone 1 is denoted by m 1 (n), the information coming into Microphone 2 is similarly labeled 10 m 2 (n), and the GEMS sensor is assumed available to determine voiced speech areas In the z (digital frequency) domain, these signals are represented as Mj(z) and M 2 (z). Then M, (z)= S(z)+ N 2 (z)
M
2 (z) = N(z)+ S 2 (z) with 15
N
2 (z)= N(z)H, (z) S2(z)= S(z)H2(z) so that M, (z) = S(z)+ N(z)H, (z)
M
2 (z)= N(z) + S(z)H 2 (z) (1) This is the general case for all two microphone systems. There is always going to be some leakage of noise into Mic 1, and some leakage of signal into Mic 2 Equation 1 20 has four unknowns and only two relationships and cannot be solved explicitly. However, there is another way to solve for some of the unknowns in Equation 1. Examine the case where the signal is not being generated - that is, where the GEMS signal indicates voicing is not occurring In this case, s(n) = S(z) = 0, and Equation I reduces to 25 M,, (z) = N(z)H, (z)
M
2 n (z) = N(z) where the n subscript on the M variables indicate that only noise is being received This leads to M,, (z) = M 2 n (z)H, (z)
M
1 (z) (2) HV, (Z) = ,,"Z
M
9 (z) 22 WO 2010/048635 PCT/US2009/062129
H
1 (z) can be calculated using any of the available system identification algorithms and the microphone outputs when only noise is being received The calculation can be done adaptively, so that if the noise changes significantly H 1 (z) can be recalculated quickly. 5 With a solution for one of the unknowns in Equation 1, solutions can be found for another, H 2 (z), by using the amplitude of the GEMS or similar device along with the amplitude of the two microphones. When the GEMS indicates voicing, but the recent (less than 1 second) history of the microphones indicate low levels of noise, assume that n(s) = N(z) - 0. Then Equation 1 reduces to 10 M' (z) = S(z) M,, (z) = S(z)H 2 (z) which in turn leads to M'- (Z) = M (Z)z H2(Z)= M" (z) Mv'f(z) which is the inverse of the H 1 (z) calculation, but note that different inputs are being used. 15 After calculating H1(z) and H 2 (z) above, they are used to remove the noise from the signal Rewrite Equation 1 as S(z)= M (z)- N(z)H, (z) N(z) = M2 (z) - S(z)H2 (z) S(z)= M (z)-[v1 2 (z)-S(z)H 2 (z)]H,(z) S(z)[1 - H 2 (z)H, (z)]= M, (z)- M 2 (z)H, (z) and solve for S(z) as: M,) /1(z) - M2(z)H, (Z) S(z)= 1 - H 2 (z)H (z) 20 In practice H 2 (z) is usually quite small, so that H 2 (z)H,(z) <<1, and S(z) - M,(Z)- M (z)HI(z), obviating the need for the H 2 (z) calculation. With reference to Figure 13 and Figure 14, the PSAD system is described, As sound waves propagate, they normally lose energy as they travel due to diffraction 25 and dispersion. Assuming the sound waves originate from a point source and radiate 23 WO 2010/048635 PCT/US2009/062129 isotropically, their amplitude will decrease as a function of 1/r, where r is the distance from the originating point. This function of 1/r proportional to amplitude is the worst case, if confined to a smaller area the reduction will be less. However it is an adequate model for the configurations of interest, specifically the propagation of noise 5 and speech to microphones located somewhere on the user's head, Figure 18 is a microphone array for use under an embodiment of the PSAD system. Placing the microphones Mic 1 and Mic 2 in a linear array with the mouth on the array midline, the difference in signal strength in Mic 1 and Mic 2 (assuming the microphones have identical frequency responses) will be proportional to both di and 10 Ad. Assuming a 1/r (or in this case 1/d) relationship, it is seen that Mic1 d +±Ad AM =AH 1 (z)c , Mic2l d] where AM is the difference in gain between Mic 1 and Mic 2 and therefore H 1 (z), as above in Equation 2. The variable d 1 is the distance from Mic 1 to the speech or noise source 15 Figure 19 is a plot 1900 of AM versus d, for several Ad values, under an embodiment, It is clear that as Ad becomes larger and the noise source is closer, AM becomes larger. The variable Ad will change depending on the orientation to the speech/noise source, from the maximum value on the array midline to zero perpendicular to the array midline From the plot 1900 it is clear that for small Ad and 20 for distances over approximately 30 centimeters (cm), AM is close to unity, Since most noise sources are farther away than 30 cm and are unlikely to be on the midline on the array, it is probable that when calculating H 1 (z) as above in Equation 2, AM (or equivalently the gain of H 1 (z)) will be close to unity, Conversely, for noise sources that are close (within a few centimeters), there could be a substantial difference in 25 gain depending on which microphone is closer to the noise If the "noise" is the user speaking, and Mic 1 is closer to the mouth than Mic 2, the gain increases Since environmental noise normally originates much farther away from the user's head than speech, noise will be found during the time when the gain of
H
1 (z) is near unity or some fixed value, and speech can be found after a sharp rise in 30 gain The speech can be unvoiced or voiced, as long as it is of sufficient volume compared to the surrounding noise. The gain will stay somewhat high during the speech portions, then descend quickly after speech ceases. The rapid increase and 24 WO 2010/048635 PCT/US2009/062129 decrease in the gain of H1(z) should be sufficient to allow the detection of speech under almost any circumstances The gain in this example is calculated by the sum of the absolute value of the filter coefficients. This sum is not equivalent to the gain, but the two are related in that a rise in the sum of the absolute value reflects a rise in the 5 gain. As an example of this behavior, Figure 20 shows a plot 2000 of the gain parameter 2002 as the sum of the absolute values of H 1 (z) and the acoustic data 2004 or audio from microphone 1 The speech signal was an utterance of the phrase "pop pan", repeated twice, The evaluated bandwidth included the frequency range from 10 2500 Hz to 3500 Hz, although 1500Hz to 2500 Hz was additionally used in practice Note the rapid increase in the gain when the unvoiced speech is first encountered, then the rapid return to normal when the speech ends. The large changes in gain that result from transitions between noise and speech can be detected by any standard signal processing techniques. The standard deviation of the last few gain calculations 15 is used, with thresholds being defined by a running average of the standard deviations and the standard deviation noise floor. The later changes in gain for the voiced speech are suppressed in this plot 2000 for clarity Figure 21 is an alternative plot 2100 of acoustic data presented in Figure 20. The data used to form plot 2000 is presented again in this plot 2100, along with audio 20 data 2104 and GEMS data 2106 without noise to make the unvoiced speech apparent. The voiced signal 2102 has three possible values. 0 for noise, 1 for unvoiced, and 2 for voiced. Denoising is only accomplished when V = 0. It is clear that the unvoiced speech is captured very well, aside from two single dropouts in the unvo iced detection near the end of each "pop". However, these single-window 25 dropouts are not common and do not significantly affect the denoising algorithm. They can easily be removed using standard smoothing techniques. What is not clear from this plot 2100 is that the PSAD system functions as an automatic backup to the NAVSAD. This is because the voiced speech (since it has the same spatial relationship to the mics as the unvoiced) will be detected as unvoiced 30 if the sensor or NAVSAD system fail for any reason. The voiced speech will be misclassified as unvoiced, but the denoising will still not take place, preserving the quality of the speech signal. However, this automatic backup of the NAVSAD system functions best in an environment with low noise (approximately 10+ dB SNR), as high amounts (10 dB of 25 WO 2010/048635 PCT/US2009/062129 SNR or less) of acoustic noise can quickly overwhelm any acoustic-only unvoiced detector, including the PSAD. This is evident in the difference in the voiced signal data 1702 and 2102 shown in plots 1700 and 1200 of Figures 17 and 21, respectively, where the same utterance is spoken, but the data of plot 1700 shows no 5 unvoiced speech because the unvoiced speech is undetectable, This is the desired behavior when performing denoising, since if the unvoiced speech is not detectable then it will not significantly affect the denoising process. Using the Pathfinder system to detect unvoiced speech ensures detection of any unvoiced speech loud enough to distort the denoising 10 Regarding hardware considerations, and with reference to Figure 18, the configuration of the microphones can have an effect on the change in gain associated with speech and the thresholds needed to detect speech. In general, each configuration will require testing to determine the proper thresholds, but tests with two very different microphone configurations showed the same thresholds and other 15 parameters to work well. The first microphone set had the signal microphone near the mouth and the noise microphone several centimeters away at the ear, while the second configuration placed the noise and signal microphones back-to-back within a few centimeters of the mouth. The results presented herein were derived using the first microphone configuration, but the results using the other set are virtually identical, 20 so the detection algorithm is relatively robust with respect to microphone placement A number of configurations are possible using the NAVSAD and PSAD systems to detect voiced and unvoiced speech, One configuration uses the NAVSAD system (non-acoustic only) to detect voiced speech along with the PSAD system to detect unvoiced speech; the PSAD also functions as a backup to the NAVSAD system 25 for detecting voiced speech An alternative configuration uses the NAVSAD system (non-acoustic correlated with acoustic) to detect voiced speech along with the PSAD system to detect unvoiced speech; the PSAD also functions as a backup to the NAVSAD system for detecting voiced speech Another alternative configuration uses the PSAD system to detect both voiced and unvoiced speech. 30 While the systems described above have been described with reference to separating voiced and unvoiced speech from background acoustic noise, there are no reasons more complex classifications can not be made. For more in-depth characterization of speech, the system can bandpass the information from Mic 1 and Mic 2 so that it is possible to see which bands in the Mic 1 data are more heavily 26 WO 2010/048635 PCT/US2009/062129 composed of noise and which are more weighted with speech, Using this knowledge, it is possible to group the utterances by their spectral characteristics similar to conventional acoustic methods; this method would work better in noisy environments. As an example, the "k" in "kick" has significant frequency content form 500 Hz 5 to 4000 Hz, but a "sh" in "she" only contains significant energy from 1700-4000 Hz Voiced speech could be classified in a similar manner, For instance, an /i/ ("ee") has significant energy around 300 Hz and 2500 Hz, and an /al ("ah") has energy at around 900 Hz and 1200 Hz. This ability to discriminate unvoiced and voiced speech in the presence of noise is, thus, very useful, 10 Dual omnidirectional microphone array (DOMA) A dual omnidirectional microphone array (DOMA) that provides improved noise suppression is now described. Compared to conventional arrays and algorithms, which seek to reduce noise by nulling out noise sources, the array of 1.5 an embodiment is used to form two distinct virtual directional microphones, as described in detail above. The two virtual microphones are configured to have very similar noise responses and very dissimilar speech responses. The only null formed by the DOMA is one used to remove the speech of the user from V 2 . The two virtual microphones of an embodiment can be paired with an adaptive filter 20 algorithm and/or VAD algorithm, as described in detail above, to significantly reduce the noise without distorting the speech, significantly improving the SNR of the desired speech over conventional noise suppression systems. The embodiments described herein are stable in operation, flexible with respect to virtual microphone pattern choice, and have proven to be robust with respect to 25 speech source-to-array distance and orientation as well as temperature and calibration techniques. Unless otherwise specified, the following terms used in describing the DOMA of an embodiment have the corresponding meanings in addition to any meaning or understanding they may convey to one skilled in the art. 30 The term "bleedthrough" means the undesired presence of noise during speech. The term "denoising" means removing unwanted noise from Micd, and also refers to the amount of reduction of noise energy in a signal in decibels (dB). The term "devoicing" means removing/distorting the desired speech from 35 Mic. 27 WO 2010/048635 PCT/US2009/062129 The term "directional microphone (DM)" means a physical directional microphone that is vented on both sides of the sensing diaphragm. The term "Mic1 (Ml)" means a general designation for an adaptive noise suppression system microphone that usually contains more speech than noise. 5 The term "Mic2 (M2)" means a general designation for an adaptive noise suppression system microphone that usually contains more noise than speech. The term "noise" means unwanted environmental acoustic noise. The term "null" means a zero or minima in the spatial response of a physical or virtual directional microphone. 10 The term "01" means a first physical omnidirectional microphone used to form a microphone array. The term "02" means a second physical omnidirectional microphone used to form a microphone array. The term "speech" means desired speech of the user. 15 The term "Skin Surface Microphone (SSM)" is a microphone used in an earpiece (e.g., the Jawbone earpiece available from Aliph of San Francisco, California) to detect speech vibrations on the user's skin. The term "Vi" means the virtual directional "speech" microphone, which has no nulls. 20 The term "V 2 " means the virtual directional "noise" microphone, which has a null for the user's speech. The term "Voice Activity Detection (VAD) signal" means a signal indicating when user speech is detected. The term "virtual microphones (VM)" or "virtual directional microphones" 25 means a microphone constructed using two or more omnidirectional microphones and associated signal processing. Figure 22 is a two-microphone adaptive noise suppression system 2200, under an embodiment. The two-microphone system 2200 including the combination of physical microphones MIC 1 and MIC 2 along with the processing 30 or circuitry components to which the microphones couple (described in detail below, but not shown in this figure) is referred to herein as the dual omnidirectional microphone array (DOMA) 2210, but the embodiment is not so limited. Referring to Figure 22, in analyzing the single noise source 2201 and the direct path to the microphones, the total acoustic information coming into MIC 1 35 (2202, which can be an physical or virtual microphone) is denoted by m 1 (n). The 28 WO 2010/048635 PCT/US2009/062129 total acoustic information coming into MIC 2 (2203, which can also be an physical or virtual microphone) is similarly labeled m 2 (n). In the z (digital frequency) domain, these are represented as M 1 (z) and M 2 (z). Then,
M
1 (z) = S(z)+N 2 (z) 5 M2 (z)=N(z)+S 2 (z) with
N
2 (z)=N(z)H ,(z)
S
2 (z)=S(z)H 2 (z), so that 10 M ,(z)=S(z)+N(z)H (z)
M
2 (z)=N(z)+S(z)H 2 (z). This is the general case for all two microphone systems. Equation 1 has four unknowns and only two known relationships and therefore cannot be solved 15 explicitly. However, there is another way to solve for some of the unknowns in Equation 1. The analysis starts with an examination of the case where the speech is not being generated, that is, where a signal from the VAD subsystem 2204 (optional) equals zero. In this case, s(n) = S(z) = 0, and Equation 1 reduces to 20 MI. (z)= N(z)H, (z) M 2N (z) =N(z) , where the N subscript on the M variables indicate that only noise is being received. This leads to 25 MIN(Z)=M 2 N(z)Hj(z) H ,(z) = MIN. Eq. 2
M
2 N (Z) The function H1(z) can be calculated using any of the available system identification algorithms and the microphone outputs when the system is certain 30 that only noise is being received. The calculation can be done adaptively, so that the system can react to changes in the noise. A solution is now available for H 1 (z), one of the unknowns in Equation 1. The final unknown, H 2 (z), can be determined by using the instances where speech is being produced and the VAD equals one. When this is occurring, but the recent 29 WO 2010/048635 PCT/US2009/062129 (perhaps less than 1 second) history of the microphones indicate low levels of noise, it can be assumed that n(s) = N(z) ~ 0. Then Equation 1 reduces to MIs(z)=S(z) M2s (z) =S(z)H2(z), 5 which in turn leads to M2s (z)=Mis (z)H2(z)
M
2 s (z)
H
2 (z)= , Mis(z) which is the inverse of the H 1 (z) calculation. However, it is noted that different 10 inputs are being used (now only the speech is occurring whereas before only the noise was occurring). While calculating H 2 (z), the values calculated for H1(z) are held constant (and vice versa) and it is assumed that the noise level is not high enough to cause errors in the H 2 (z) calculation. After calculating H 1 (z) and H 2 (z), they are used to remove the noise from 15 the signal. If Equation 1 is rewritten as S(z) = MI (z) - N(z)H , (z) N(z)=M 2 (z)--S(z)H 2 (z) S(z)=M 1 (z)-[M2(z)-S(z)H 2 (z)]H (z) 20 S(z)[1- H 2(z)Hi(z)]= M i(z)- M2 (z)Hli(z), then N(z) may be substituted as shown to solve for S(z) as z M 1 (z)-M 2 (z)H 1 (z) r I 1- H, (z)H2 (z) 25 If the transfer functions H 1 (z) and H 2 (z) can be described with sufficient accuracy, then the noise can be completely removed and the original signal recovered. This remains true without respect to the amplitude or spectral characteristics of the noise. If there is very little or no leakage from the speech source into M 2 , then H 2 (z)~~ 0 and Equation 3 reduces to 30 S(z)~M 1 (z)-M 2 (z)H,(z), Eq. 4 30 WO 2010/048635 PCT/US2009/062129 Equation 4 is much simpler to implement and is very stable, assuming
H
1 (z) is stable. However, if significant speech energy is in M 2 (z), devoicing can occur. In order to construct a well-performing system and use Equation 4, consideration is given to the following conditions: 5 R1. Availability of a perfect (or at least very good) VAD in noisy conditions R2. Sufficiently accurate Hi(z) R3. Very small (ideally zero) H 2 (z). R4. During speech production, Hi(z) cannot change substantially. 10 R5. During noise, H 2 (z) cannot change substantially. Condition R1 is easy to satisfy if the SNR of the desired speech to the unwanted noise is high enough. "Enough" means different things depending on the method of VAD generation. If a VAD vibration sensor is used, as in Burnett 15 7,256,048, accurate VAD in very low SNRs (-10 dB or less) is possible. Acoustic only methods using information from 01 and 02 can also return accurate VADs, but are limited to SNRs of -3 dB or greater for adequate performance. Condition R5 is normally simple to satisfy because for most applications the microphones will not change position with respect to the user's mouth very often 20 or rapidly. In those applications where it may happen (such as hands-free conferencing systems) it can be satisfied by configuring Mic2 so that H 2 (z)~ 0 . Satisfying conditions R2, R3, and R4 are more difficult but are possible given the right combination of V 1 and V 2 . Methods are examined below that have proven to be effective in satisfying the above, resulting in excellent noise 25 suppression performance and minimal speech removal and distortion in an embodiment. The DOMA, in various embodiments, can be used with the Pathfinder system as the adaptive filter system or noise removal. The Pathfinder system, available from AliphCom, San Francisco, CA, is described in detail in other patents 30 and patent applications referenced herein. Alternatively, any adaptive filter or noise removal algorithm can be used with the DOMA in one or more various alternative embodiments or configurations. When the DOMA is used with the Pathfinder system, the Pathfinder system generally provides adaptive noise cancellation by combining the two microphone 35 signals (e.g., Micd, Mic2) by filtering and summing in the time domain. The 31 WO 2010/048635 PCT/US2009/062129 adaptive filter generally uses the signal received from a first microphone of the DOMA to remove noise from the speech received from at least one other microphone of the DOMA, which relies on a slowly varying linear transfer function between the two microphones for sources of noise. Following processing of the 5 two channels of the DOMA, an output signal is generated in which the noise content is attenuated with respect to the speech content, as described in detail below. Figure 23 is a generalized two-microphone array (DOMA) including an array 2301/2302 and speech source S configuration, under an embodiment. 10 Figure 24 is a system 2400 for generating or producing a first order gradient microphone V using two omnidirectional elements 01 and 02, under an embodiment. The array of an embodiment includes two physical microphones 2301 and 2302 (e.g., omnidirectional microphones) placed a distance 2do apart and a speech source 2300 is located a distance d away at an angle of 0. This 15 array is axially symmetric (at least in free space), so no other angle is needed. The output from each microphone 2301 and 2302 can be delayed (zi and z 2 ), multiplied by a gain (A 1 and A 2 ), and then summed with the other as demonstrated in Figure 24. The output of the array is or forms at least one virtual microphone, as described in detail below. This operation can be over any 20 frequency range desired. By varying the magnitude and sign of the delays and gains, a wide variety of virtual microphones (VMs), also referred to herein as virtual directional microphones, can be realized. There are other methods known to those skilled in the art for constructing VMs but this is a common one and will be used in the enablement below. 25 As an example, Figure 25 is a block diagram for a DOMA 2500 including two physical microphones configured to form two virtual microphones V 1 and V 2 , under an embodiment. The DOMA includes two first order gradient microphones
V
1 and V 2 formed using the outputs of two microphones or elements 01 and 02 (2301 and 2302), under an embodiment. The DOMA of an embodiment includes 30 two physical microphones 2301 and 2302 that are omnidirectional microphones, as described above with reference to Figures 23 and 24. The output from each microphone is coupled to a processing component 2502, or circuitry, and the processing component outputs signals representing or corresponding to the virtual microphones V 1 and V 2 . 32 WO 2010/048635 PCT/US2009/062129 In this example system 2500, the output of physical microphone 2301 is coupled to processing component 2502 that includes a first processing path that includes application of a first delay z 11 and a first gain All and a second processing path that includes application of a second delay z 12 and a second gain A 12 . The 5 output of physical microphone 2302 is coupled to a third processing path of the processing component 2502 that includes application of a third delay z 21 and a third gain A 21 and a fourth processing path that includes application of a fourth delay z 22 and a fourth gain A 2 2 . The output of the first and third processing paths is summed to form virtual microphone V 1 , and the output of the second and fourth 10 processing paths is summed to form virtual microphone V 2 , As described in detail below, varying the magnitude and sign of the delays and gains of the processing paths leads to a wide variety of virtual microphones (VMs), also referred to herein as virtual directional microphones, can be realized. While the processing component 2502 described in this example includes four 15 processing paths generating two virtual microphones or microphone signals, the embodiment is not so limited. For example, Figure 26 is a block diagram for a DOMA 2600 including two physical microphones configured to form N virtual microphones V 1 through VN, where N is any number greater than one, under an embodiment. Thus, the DOMA can include a processing component 2602 having 20 any number of processing paths as appropriate to form a number N of virtual microphones. The DOMA of an embodiment can be coupled or connected to one or more remote devices. In a system configuration, the DOMA outputs signals to the remote devices. The remote devices include, but are not limited to, at least one 25 of cellular telephones, satellite telephones, portable telephones, wireline telephones, Internet telephones, wireless transceivers, wireless communication radios, personal digital assistants (PDAs), personal computers (PCs), headset devices, head-worn devices, and earpieces. Furthermore, the DOMA of an embodiment can be a component or 30 subsystem integrated with a host device. In this system configuration, the DOMA outputs signals to components or subsystems of the host device. The host device includes, but is not limited to, at least one of cellular telephones, satellite telephones, portable telephones, wireline telephones, Internet telephones, wireless transceivers, wireless communication radios, personal digital assistants 33 WO 2010/048635 PCT/US2009/062129 (PDAs), personal computers (PCs), headset devices, head-worn devices, and earpieces. As an example, Figure 27 is an example of a headset or head-worn device 2700 that includes the DOMA, as described herein, under an embodiment. The 5 headset 2700 of an embodiment includes a housing having two areas or receptacles (not shown) that receive and hold two microphones (e.g., 01 and 02). The headset 2700 is generally a device that can be worn by a speaker 2702, for example, a headset or earpiece that positions or holds the microphones in the vicinity of the speaker's mouth. The headset 2700 of an embodiment places a 10 first physical microphone (e.g., physical microphone 01) in a vicinity of a speaker's lips. A second physical microphone (e.g., physical microphone 02) is placed a distance behind the first physical microphone. The distance of an embodiment is in a range of a few centimeters behind the first physical microphone or as described herein (e.g., described with reference to Figures 22 15 26). The DOMA is symmetric and is used in the same configuration or manner as a single close-talk microphone, but is not so limited. Figure 28 is a flow diagram for denoising 2800 acoustic signals using the DOMA, under an embodiment. The denoising 2800 begins by receiving 2802 acoustic signals at a first physical microphone and a second physical microphone. 20 In response to the acoustic signals, a first microphone signal is output from the first physical microphone and a second microphone signal is output from the second physical microphone 2804. A first virtual microphone is formed 2806 by generating a first combination of the first microphone signal and the second microphone signal. A second virtual microphone is formed 2808 by generating a 25 second combination of the first microphone signal and the second microphone signal, and the second combination is different from the first combination. The first virtual microphone and the second virtual microphone are distinct virtual directional microphones with substantially similar responses to noise and substantially dissimilar responses to speech. The denoising 2800 generates 2810 30 output signals by combining signals from the first virtual microphone and the second virtual microphone, and the output signals include less acoustic noise than the acoustic signals. Figure 29 is a flow diagram for forming 2900 the DOMA, under an embodiment. Formation 2900 of the DOMA includes forming 2902 a physical 35 microphone array including a first physical microphone and a second physical 34 WO 2010/048635 PCT/US2009/062129 microphone. The first physical microphone outputs a first microphone signal and the second physical microphone outputs a second microphone signal. A virtual microphone array is formed 2904 comprising a first virtual microphone and a second virtual microphone. The first virtual microphone comprises a first 5 combination of the first microphone signal and the second microphone signal. The second virtual microphone comprises a second combination of the first microphone signal and the second microphone signal, and the second combination is different from the first combination. The virtual microphone array including a single null oriented in a direction toward a source of speech of a human speaker. 10 The construction of VMs for the adaptive noise suppression system of an embodiment includes substantially similar noise response in V 1 and V 2 . Substantially similar noise response as used herein means that H 1 (z) is simple to model and will not change much during speech, satisfying conditions R2 and R4 described above and allowing strong denoising and minimized bleedthrough. 15 The construction of VMs for the adaptive noise suppression system of an embodiment includes relatively small speech response for V 2 . The relatively small speech response for V 2 means that H- (z) e 0, which will satisfy conditions R3 and R5 described above. The construction of VMs for the adaptive noise suppression system of an 20 embodiment further includes sufficient speech response for Vi so that the cleaned speech will have significantly higher SNR than the original speech captured by 0. The description that follows assumes that the responses of the omnidirectional microphones 01 and 02 to an identical acoustic source have been normalized so that they have exactly the same response (amplitude and phase) to 25 that source. This can be accomplished using standard microphone array methods (such as frequency-based calibration) well known to those versed in the art. Referring to the condition that construction of VMs for the adaptive noise suppression system of an embodiment includes relatively small speech response for V 2 , it is seen that for discrete systems V 2 (z) can be represented as: 30
V
2 (z)= 0 2 (z)- zPO, (z) where d d2 35 WO 2010/048635 PCT/US2009/062129 y = d2 -djfi. (samples) dI = d -2dd 0 cos(9)+d d2 = d + 2ddo cos(0)+ d The distances di and d 2 are the distance from 01 and 02 to the speech source (see 5 Figure 23), respectively, and y is their difference divided by c, the speed of sound, and multiplied by the sampling frequency fs. Thus y is in samples, but need not be an integer. For non-integer y, fractional-delay filters (well known to those versed in the art) may be used. It is important to note that the P above is not the conventional P used to 10 denote the mixing of VMs in adaptive beamforming; it is a physical variable of the system that depends on the intra-microphone distance do (which is fixed) and the distance d, and angle 0, which can vary. As shown below, for properly calibrated microphones, it is not necessary for the system to be programmed with the exact P of the array. Errors of approximately 10-15% in the actual p (i.e. the p used by 15 the algorithm is not the P of the physical array) have been used with very little degradation in quality. The algorithmic value of P may be calculated and set for a particular user or may be calculated adaptively during speech production when little or no noise is present. However, adaptation during use is not required for nominal performance. 20 Figure 30 is a plot of linear response of virtual microphone V 2 with p 0.8 to a 1 kHz speech source at a distance of 0.1 m, under an embodiment. The null in the linear response of virtual microphone V 2 to speech is located at 0 degrees, where the speech is typically expected to be located. Figure 31 is a plot of linear response of virtual microphone V 2 with P = 0.8 to a 1 kHz noise source at a 25 distance of 1.0 m, under an embodiment. The linear response of V 2 to noise is devoid of or includes no null, meaning all noise sources are detected. The above formulation for V 2 (z) has a null at the speech location and will therefore exhibit minimal response to the speech. This is shown in Figure 30 for an array with do = 10.7 mm and a speech source on the axis of the array (0 = 0) 30 at 10 cm (p 0.8). Note that the speech null at zero degrees is not present for noise in the far field for the same microphone, as shown in Figure 31 with a noise source distance of approximately 1 meter. This insures that noise in front 36 WO 2010/048635 PCT/US2009/062129 of the user will be detected so that it can be removed. This differs from conventional systems that can have difficulty removing noise in the direction of the mouth of the user. The V 1 (z) can be formulated using the general form for V 1 (z): 5 Vl(z)=AOl (z).z-dA -UaO2(z) zdB Since
V
2 (z)= 0 2 (z)- z- pOI(z) 10 and, since for noise in the forward direction 02N(Z)= OIN(z).zy then
V
2 N (z)= 0 1N (z). - Z OIN (Z) V2N (Z) = (1 - P)(O IN (Z) - Z7 15 If this is then set equal to V 1 (z) above, the result is VIN (Z AOIN (Z)*Z - Z 1^ - B IN (Z) - Z - z Zd. =( -P)(O IN (Z)'Z thus we may set 20 dA =Y dB= 0 aA 1 GLB = to get 25 VI(z)=O(z)-z- -P0 2 (z) The definitions for V 1 and V 2 above mean that for noise H 1 (z) is: H (z)= V 1 (z) _ -pO 2 (z)+O 1 (z)- z-'
V
2 (z) 0 2 (z)- z~7pOI(z) 30 37 WO 2010/048635 PCT/US2009/062129 which, if the amplitude noise responses are about the same, has the form of an allpass filter. This has the advantage of being easily and accurately modeled, especially in magnitude response, satisfying R2. 5 This formulation assures that the noise response will be as similar as possible and that the speech response will be proportional to (1-p2). Since P is the ratio of the distances from 01 and 02 to the speech source, it is affected by the size of the array and the distance from the array to the speech source. Figure 32 is a plot of linear response of virtual microphone V 1 with p 0.8 10 to a 1 kHz speech source at a distance of 0.1 m, under an embodiment. The linear response of virtual microphone Vi to speech is devoid of or includes no null and the response for speech is greater than that shown in Figure 25. Figure 33 is a plot of linear response of virtual microphone V 1 with 3 0.8 to a 1 kHz noise source at a distance of 1.0 m, under an embodiment. The linear 15 response of virtual microphone V 1 to noise is devoid of or includes no null and the response is very similar to V 2 shown in Figure 26. Figure 34 is a plot of linear response of virtual microphone V 1 with p 0.8 to a speech source at a distance of 0.1 m for frequencies of 100, 500, 1000, 2000, 3000, and 4000 Hz, under an embodiment. Figure 35 is a plot showing 20 comparison of frequency responses for speech for the array of an embodiment and for a conventional cardioid microphone. The response of Vi to speech is shown in Figure 32, and the response to noise in Figure 33. Note the difference in speech response compared to V 2 shown in Figure 30 and the similarity of noise response shown in Figure 31. 25 Also note that the orientation of the speech response for Vi shown in Figure 32 is completely opposite the orientation of conventional systems, where the main lobe of response is normally oriented toward the speech source. The orientation of an embodiment, in which the main lobe of the speech response of V1 is oriented away from the speech source, means that the speech sensitivity of V 1 is lower 30 than a normal directional microphone but is flat for all frequencies within approximately +-30 degrees of the axis of the array, as shown in Figure 34. This flatness of response for speech means that no shaping postfilter is needed to restore omnidirectional frequency response. This does come at a price - as shown in Figure 35, which shows the speech response of V 1 with 3 = 0.8 and the speech 35 response of a cardioid microphone. The speech response of V 1 is approximately 0 38 WO 2010/048635 PCT/US2009/062129 to "13 dB less than a normal directional microphone between approximately 500 and 7500 Hz and approximately 0 to 10+ dB greater than a directional microphone below approximately 500 Hz and above 7500 Hz for a sampling frequency of approximately 16000 Hz. However, the superior noise suppression 5 made possible using this system more than compensates for the initially poorer SNR. It should be noted that Figures 30-33 assume the speech is located at approximately 0 degrees and approximately 10 cm, P = 0.8, and the noise at all angles is located approximately 1.0 meter away from the midpoint of the array. 10 Generally, the noise distance is not required to be 1 m or more, but the denoising is the best for those distances. For distances less than approximately 1 m, denoising will not be as effective due to the greater dissimilarity in the noise responses of V 1 and V 2 . This has not proven to be an impediment in practical use - in fact, it can be seen as a feature. Any "noise" source that is -10 cm away 15 from the earpiece is likely to be desired to be captured and transmitted. The speech null of V 2 means that the VAD signal is no longer a critical component. The VAD's purpose was to ensure that the system would not train on speech and then subsequently remove it, resulting in speech distortion. If, however, V 2 contains no speech, the adaptive system cannot train on the speech 20 and cannot remove it. As a result, the system can denoise all the time without fear of devoicing, and the resulting clean audio can then be used to generate a VAD signal for use in subsequent single-channel noise suppression algorithms such as spectral subtraction. In addition, constraints on the absolute value of
H
1 (z) (i.e. restricting it to absolute values less than two) can keep the system 25 from fully training on speech even if it is detected. In reality, though, speech can be present due to a mis-located V 2 null and/or echoes or other phenomena, and a VAD sensor or other acoustic-only VAD is recommended to minimize speech distortion. Depending on the application, P and y may be fixed in the noise suppression 30 algorithm or they can be estimated when the algorithm indicates that speech production is taking place in the presence of little or no noise. In either case, there may be an error in the estimate of the actual p and y of the system. The following description examines these errors and their effect on the performance of the system. As above, "good performance" of the system indicates that there is 35 sufficient denoising and minimal devoicing. 39 WO 2010/048635 PCT/US2009/062129 The effect of an incorrect P and y on the response of V 1 and V 2 can be seen by examining the definitions above: V, (z)= 0 1 (z). z-T - PT02(z)
V
2 (z)= O 2 (z)- z PTOI(z) 5 where PT and yT denote the theoretical estimates of p and y used in the noise suppression algorithm. In reality, the speech response of 02 is 0 2S (z)= PROIs (z) ZR 10 where PR and yR denote the real p and y of the physical system. The differences between the theoretical and actual values of p and y can be due to mis-location of the speech source (it is not where it is assumed to be) and/or a change in air temperature (which changes the speed of sound). Inserting the actual response 15 of 02 for speech into the above equations for V 1 and V 2 yields Vis(z)= 0is(z)[z- -P zT-PRZ
V
2 s(z)=OIs(z)p[Rz-YR _PTZ-7T If the difference in phase is represented by 20 YR -YT YD And the difference in amplitude as PR = BP 1 then Vis (z)= Os (z)z~YT - B z-YD V2s(Z)=p3TOs (Z)Z-YT [BzD I] Eq. 5 25 The speech cancellation in V 2 (which directly affects the degree of devoicing) and the speech response of Vi will be dependent on both B and D. An examination of the case where D = 0 follows. Figure 36 is a plot showing speech response for V 1 (top, dashed) and V 2 (bottom, solid) versus B with d, assumed to 30 be 0.1 m, under an embodiment. This plot shows the spatial null in V 2 to be 40 WO 2010/048635 PCT/US2009/062129 relatively broad. Figure 37 is a plot showing a ratio of V 1
/V
2 speech responses shown in Figure 31 versus B, under an embodiment. The ratio of V 1
/V
2 is above 10 dB for all 0.8 < B < 1.1, and this means that the physical P of the system need not be exactly modeled for good performance. Figure 38 is a plot of B versus 5 actual d assuming that d = 10 cm and theta = 0, under an embodiment. Figure 39 is a plot of B versus theta with d, = 10 cm and assuming ds = 10 cm, under an embodiment. In Figure 36, the speech response for V 1 (upper, dashed) and V 2 (lower, solid) compared to 01 is shown versus B when ds is thought to be approximately 10 10 cm and 9 = 0. When B = 1, the speech is absent from V 2 . In Figure 37, the ratio of the speech responses in Figure 31 is shown. When 0.8 < B < 1.1, the
V
1
/V
2 ratio is above approximately 10 dB - enough for good performance. Clearly, if D = 0, B can vary significantly without adversely affecting the performance of the system. Again, this assumes that calibration of the 15 microphones so that both their amplitude and phase response is the same for an identical source has been performed. The B factor can be non-unity for a variety of reasons. Either the distance to the speech source or the relative orientation of the array axis and the speech source or both can be different than expected. If both distance and angle 20 mismatches are included for B, then PR diR -2dSR dOCOS(0R)±d d +2d 0 COS(O T)+ d P1 d S R +2dSRdo COS(OR )T d Id - 2dSTdO CoS(OT )+d where again the T subscripts indicate the theorized values and R the actual 25 values. In Figure 38, the factor B is plotted with respect to the actual d, with the assumption that ds = 10 cm and 0 = 0. So, if the speech source in on-axis of the array, the actual distance can vary from approximately 5 cm to 18 cm without significantly affecting performance - a significant amount. Similarly, Figure 39 shows what happens if the speech source is located at a distance of approximately 30 10 cm but not on the axis of the array. In this case, the angle can vary up to approximately +-55 degrees and still result in a B less than 1.1, assuring good performance. This is a significant amount of allowable angular deviation. If there is both angular and distance errors, the equation above may be used to determine 41 WO 2010/048635 PCT/US2009/062129 if the deviations will result in adequate performance. Of course, if the value for Pr is allowed to update during speech, essentially tracking the speech source, then B can be kept near unity for almost all configurations. An examination follows of the case where B is unity but D is nonzero. This 5 can happen if the speech source is not where it is thought to be or if the speed of sound is different from what it is believed to be. From Equation 5 above, it can be sees that the factor that weakens the speech null in V 2 for speech is N(z)= BZ-YD _I or in the continuous s domain 10 N(s)= Be-Ds Since y is the time difference between arrival of speech at V 1 compared to V 2 , it can be errors in estimation of the angular location of the speech source with respect to the axis of the array and/or by temperature changes. Examining the 15 temperature sensitivity, the speed of sound varies with temperature as e = 331.3 +(0.606T) rn/s where T is degrees Celsius. As the temperature decreases, the speed of sound 20 also decreases. Setting 20 C as a design temperature and a maximum expected temperature range to -40 C to +60 C (-40 F to 140 F). The design speed of sound at 20 C is 343 m/s and the slowest speed of sound will be 307 m/s at -40 C with the fastest speed of sound 362 m/s at 60 C. Set the array length (2do) to be 21 mm. For speech sources on the axis of the array, the difference in travel time 25 for the largest change in the speed of sound is VtMAX - d 1 0,021 m = -K1 7.2 x10-6 seC c1 c2 3431m/s 307 m/s) or approximately 7 microseconds. The response for N(s) given B = 1 and D = 30 7.2 psec is shown in Figure 40. Figure 40 is a plot of amplitude (top) and phase (bottom) response of N(s) with B = 1 and D = -7.2 ksec, under an embodiment. The resulting phase difference clearly affects high frequencies more than low. The amplitude response is less than approximately -10 dB for all frequencies less 42 WO 2010/048635 PCT/US2009/062129 than 7 kHz and is only about -9 dB at 8 kHz. Therefore, assuming B = 1, this system would likely perform well at frequencies up to approximately 8 kHz. This means that a properly compensated system would work well even up to 8 kHz in an exceptionally wide (e.g., -40 C to 80 C) temperature range. Note that the 5 phase mismatch due to the delay estimation error causes N(s) to be much larger at high frequencies compared to low. If B is not unity, the robustness of the system is reduced since the effect from non-unity B is cumulative with that of non-zero D. Figure 41 shows the amplitude and phase response for B = 1.2 and D = 7.2 Vsec. Figure 41 is a plot 10 of amplitude (top) and phase (bottom) response of N(s) with B = 1.2 and D = 7.2 [sec, under an embodiment. Non-unity B affects the entire frequency range. Now N(s) is below approximately -10 dB only for frequencies less than approximately 5 kHz and the response at low frequencies is much larger. Such a system would still perform well below 5 kHz and would only suffer from slightly 15 elevated devoicing for frequencies above 5 kHz. For ultimate performance, a temperature sensor may be integrated into the system to allow the algorithm to adjust YT as the temperature varies. Another way in which D can be non-zero is when the speech source is riot where it is believed to be - specifically, the angle from the axis of the array to the 20 speech source is incorrect. The distance to the source may be incorrect as well, but that introduces an error in B, not D. Referring to Figure 23, it can be seen that for two speech sources (each with their own ds and 0) that the time difference between the arrival of the speech at 01 and the arrival at 02 is 25 At I(d 2 -di 1 -d 2 2 +d 21 ) C where de=d~si -2dsido cos(Oi)+ di d 1 1 = ds +2dsIdo cos(0 1 )+d0 d1 = d 2 -2ds 2 de cos(0 2 )+ d d22 = ds 2 +2ds 2 d 0 cos(0 2 )+di The V 2 speech cancellation response for 01 = 0 degrees and 02 = 30 degrees 30 and assuming that B = 1 is shown in Figure 42. Figure 42 is a plot of amplitude 43 WO 2010/048635 PCT/US2009/062129 (top) and phase (bottom) response of the effect on the speech cancellation in V 2 due to a mistake in the location of the speech source with q1 = 0 degrees and q2 = 30 degrees, under an embodiment. Note that the cancellation is still below -10 dB for frequencies below 6 kHz. The cancellation is still below approximately -10 5 dB for frequencies below approximately 6 kHz, so an error of this type will not significantly affect the performance of the system. However, if 02 is increased to approximately 45 degrees, as shown in Figure 43, the cancellation is below approximately -10 dB only for frequencies below approximately 2.8 kHz. Figure 43 is a plot of amplitude (top) and phase (bottom) response of the effect on the 10 speech cancellation in V 2 due to a mistake in the location of the speech source with q1 = 0 degrees and q2 = 45 degrees, under an embodiment. Now the cancellation is below -10 dB only for frequencies below about 2.8 kHz and a reduction in performance is expected. The poor V 2 speech cancellation above approximately 4 kHz may result in significant devoicing for those frequencies. 15 The description above has assumed that the microphones 01 and 02 were calibrated so that their response to a source located the same distance away was identical for both amplitude and phase. This is not always feasible, so a more practical calibration procedure is presented below. It is not as accurate, but is much simpler to implement. Begin by defining a filter a(z) such that: 20 z=c31 (z0,ria) where the "C" subscript indicates the use of a known calibration source. The simplest one to use is the speech of the user. Then 25 0 1 s ) =1 D C)0 2 1-) The microphone definitions are now: VI(z)= 0(z) z-Y - (z)cz)O2(z)
V
2 (z= (z)O 2 (z) - Z '(z) 1 O (z) 30 The P of the system should be fixed and as close to the real value as possible. In practice, the system is not sensitive to changes in P and errors of approximately +-5% are easily tolerated. During times when the user is 44 WO 2010/048635 PCT/US2009/062129 producing speech but there is little or no noise, the system can train a(z) to remove as much speech as possible. This is accomplished by: 1. Construct an adaptive system as shown in Figure 22 with pOis(z)z- in the 5 "MIC1" position, O 2 s(z) in the "MIC2" position, and a(z) in the H 1 (z) position. 2. During speech, adapt a(z) to minimize the residual of the system. 3. Construct V 1 (z) and V 2 (z) as above. 10 A simple adaptive filter can be used for a(z) so that only the relationship between the microphones is well modeled. The system of an embodiment trains only when speech is being produced by the user. A sensor like the SSM is invaluable in determining when speech is being produced in the absence of noise. If the speech source is fixed in position and will not vary significantly during use 15 (such as when the array is on an earpiece), the adaptation should be infrequent and slow to update in order to minimize any errors introduced by noise present during training. The above formulation works very well because the noise (far-field) responses of V 1 and V 2 are very similar while the speech (near-field) responses 20 are very different. However, the formulations for V 1 and V 2 can be varied and still result in good performance of the system as a whole. If the definitions for V 1 and
V
2 are taken from above and new variables B1 and B2 are inserted, the result is: V (z)=O(z z' -B 1 f 02(z)
V
2 (z)=0 2 (z)- z 4B 2 f 1 0 1 (z) 25 where B1 and B2 are both positive numbers or zero. If B1 and B2 are set equal to unity, the optimal system results as described above. If B1 is allowed to vary from unity, the response of V 1 is affected. An examination of the case where B2 is left at 1 and B1 is decreased follows. As B1 drops to approximately zero, V 1 30 becomes less and less directional, until it becomes a simple omnidirectional microphone when B1 = 0. Since B2 = 1, a speech null remains in V 2 , so very different speech responses remain for V 1 and V 2 . However, the noise responses are much less similar, so denoising will not be as effective. Practically, though, 45 WO 2010/048635 PCT/US2009/062129 the system still performs well. B1 can also be increased from unity and once again the system will still denoise well, just not as well as with B1 = 1. If B2 is allowed to vary, the speech null in V 2 is affected. As long as the speech null is still sufficiently deep, the system will still perform well. Practically 5 values down to approximately B2 = 0.6 have shown sufficient performance, but it is recommended to set B2 close to unity for optimal performance. Similarly, variables E and A may be introduced so that: 1- + I ')O 10 This formulation also allows the virtual microphone responses to be varied but retains the all-pass characteristic of H 1 (z). In conclusion, the system is flexible enough to operate well at a variety of 81 values, but B2 values should be close to unity to limit devoicing for best 15 performance. Experimental results for a 2do = 19 mm array using a linear P of 0.83 and B1 = B2 = 1 on a Bruel and Kjaer Head and Torso Simulator (HATS) in very loud (-85 dBA) music/speech noise environment are shown in Figure 44. The alternate microphone calibration technique discussed above was used to calibrate 20 the microphones. The noise has been reduced by about 25 dB and the speech hardly affected, with no noticeable distortion. Clearly the technique significantly increases the SNR of the original speech, far outperforming conventional noise suppression techniques. Embodiments described herein include a method comprising: forming 25 a first virtual microphone by combining a first signal of a first physical microphone and a second signal of a second physical microphone; forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone; forming a second virtual microphone by applying the filter to the first signal to generate a first 30 intermediate signal, and summing the first intermediate signal and the second signal; generating an energy ratio of energies of the first virtual microphone and the second virtual microphone; and detecting acoustic 46 WO 2010/048635 PCT/US2009/062129 voice activity of a speaker when the energy ratio is greater than a threshold value. The first virtual microphone and the second virtual microphone of an embodiment are distinct virtual directional microphones. 5 The first virtual microphone and the second virtual microphone of an embodiment have approximately similar responses to noise. The first virtual microphone and the second virtual microphone of an embodiment have approximately dissimilar responses to speech. The method of an embodiment comprises applying a calibration to at 10 least one of the first signal and the second signal. The calibration of an embodiment compensates a second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. The method of an embodiment comprises applying a delay to the first 15 intermediate signal. The delay of an embodiment is proportional to a time difference between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone. The forming of the first virtual microphone of an embodiment 20 comprises applying the filter to the second signal. The forming of the first virtual microphone of an embodiment comprises applying the calibration to the second signal. The forming of the first virtual microphone of an embodiment comprises applying the delay to the first signal. 25 The forming of the first virtual microphone by the combining of an embodiment comprises subtracting the second signal from the first signal. The filter of an embodiment is an adaptive filter. The method of an embodiment comprises adapting the filter to minimize a second virtual microphone output when only speech is being 30 received by the first physical microphone and the second physical microphone. 47 WO 2010/048635 PCT/US2009/062129 The adapting of an embodiment comprises applying a least-mean squares process. The method of an embodiment comprises generating coefficients of the filter during a period when only speech is being received by the first 5 physical microphone and the second physical microphone. The forming of the filter of an embodiment comprises generating a first quantity by applying a calibration to the second signal. The forming of the filter of an embodiment comprises generating a second quantity by applying the delay to the first signal. The forming of the filter of an 10 embodiment comprises forming the filter as a ratio of the first quantity to the second quantity. The generating of the energy ratio of an embodiment comprises generating the energy ratio for a frequency band. The generating of the energy ratio of an embodiment comprises 15 generating the energy ratio for a frequency subband. The frequency subband of an embodiment includes frequencies higher than approximately 200 Hertz (Hz). The frequency subband of an embodiment includes frequencies in a range from approximately 250 Hz to 1250 Hz. 20 The frequency subband of an embodiment includes frequencies in a range from approximately 200 Hz to 3000 Hz. The filter of an embodiment is a static filter. The forming of the filter of an embodiment comprises determining a first distance as distance between the first physical microphone and a 25 mouth of the speaker. The forming of the filter of an embodiment comprises determining a second distance as distance between the second physical microphone and the mouth. The forming of the filter of an embodiment comprises forming a ratio of the first distance to the second distance. 30 The method of an embodiment comprises generating a vector of the energy ratio versus time. 48 WO 2010/048635 PCT/US2009/062129 The first and second physical microphones of an embodiment are omnidirectional microphones. The method of an embodiment comprises positioning the first physical microphone and the second physical microphone along an axis and 5 separating the first physical microphone and the second physical microphone by a first distance. A midpoint of the axis of an embodiment is a second distance from a mouth of the speaker, wherein the mouth is located in a direction defined by an angle relative to the midpoint. 10 Embodiments described herein include a method comprising: forming a first virtual microphone; forming a filter by generating a first quantity by applying a calibration to a second signal of a second physical microphone, generating a second quantity by applying the delay to a first signal of a first physical microphone, and forming the filter as a ratio of the first 15 quantity to the second quantity; forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal, and summing the first intermediate signal and the second signal; and generating a ratio of energies of the first virtual microphone and the second virtual microphone and detecting acoustic voice activity using the 20 ratio. The first virtual microphone and the second virtual microphone of an embodiment have approximately similar responses to noise and approximately dissimilar responses to speech. The method of an embodiment comprises applying a calibration to at 25 least one of the first signal and the second signal, wherein the calibration compensates a second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. The method of an embodiment comprises applying a delay to the first 30 intermediate signal, wherein the delay is proportional to a time difference between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone. 49 WO 2010/048635 PCT/US2009/062129 The forming of the first virtual microphone of an embodiment comprises applying the filter to the second signal. The forming of the first virtual microphone of an embodiment comprises applying the calibration to the second signal. 5 The forming of the first virtual microphone of an embodiment comprises applying the delay to the first signal. The forming of the first virtual microphone by the combining of an embodiment comprises subtracting the second signal from the first signal. The filter of an embodiment is an adaptive filter. 10 The method of an embodiment comprises adapting the filter to minimize a second virtual microphone output when only speech is being received by the first physical microphone and the second physical microphone. The adapting of an embodiment comprises applying a least-mean 15 squares process. The method of an embodiment comprises generating coefficients of the filter during a period when only speech is being received by the first physical microphone and the second physical microphone. The generating of the ratio of an embodiment comprises generating 20 the ratio for a frequency band. The generating of the ratio of an embodiment comprises generating the ratio for a frequency subband. The method of an embodiment comprises generating a vector of the ratio versus time. 25 Embodiments described herein include a method comprising: forming a first virtual microphone by generating a first combination of a first signal and a second signal, wherein the first signal is received from a first physical microphone and the second signal is received from a second physical microphone; forming a filter by generating a first quantity by 30 applying a calibration to at least one of the first signal and the second signal, generating a second quantity by applying a delay to the first signal, and forming the filter as a ratio of the first quantity to the second quantity; 50 WO 2010/048635 PCT/US2009/062129 and forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal and summing the first intermediate signal and the second signal; and determining a presence of acoustic voice activity of a speaker when an energy ratio of energies of the 5 first virtual microphone and the second virtual microphone is greater than a threshold value. Embodiments described herein include an acoustic voice activity detection system comprising: a first virtual microphone comprising a first combination of a first signal and a second signal, wherein the first signal is 10 received from a first physical microphone and the second signal is received from a second physical microphone; a filter, wherein the filter is formed by generating a first quantity by applying a calibration to at least one of the first signal and the second signal, generating a second quantity by applying a delay to the first signal, and forming the filter as a ratio of the first 15 quantity to the second quantity; and a second virtual microphone formed by applying the filter to the first signal to generate a first intermediate signal and summing the first intermediate signal and the second signal, wherein acoustic voice activity of a speaker is determined to be present when an energy ratio of energies of the first virtual microphone and the 20 second virtual microphone is greater than a threshold value. The first virtual microphone and the second virtual microphone of an embodiment have approximately similar responses to noise and approximately dissimilar responses to speech. A calibration is applied to the second signal of an embodiment, 25 wherein the calibration compensates a second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. The delay is applied to the first intermediate signal of an embodiment, wherein the delay is proportional to a time difference 30 between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone. 51 WO 2010/048635 PCT/US2009/062129 The first virtual microphone of an embodiment is formed by applying the filter to the second signal. The first virtual microphone of an embodiment is formed by applying the calibration to the second signal. 5 The first virtual microphone of an embodiment is formed by applying the delay to the first signal. The first virtual microphone of an embodiment is formed by subtracting the second signal from the first signal. The filter of an embodiment is an adaptive filter. 10 The filter of an embodiment is adapted to minimize a second virtual microphone output when only speech is being received by the first physical microphone and the second physical microphone. Coefficients of the filter of an embodiment are generated during a period when only speech is being received by the first physical microphone 15 and the second physical microphone. The energy ratio of an embodiment comprises an energy ratio for a frequency band. The energy ratio of an embodiment comprises an energy ratio for a frequency subband. 20 Embodiments described herein include a device comprising: a first physical microphone generating a first signal; a second physical microphone generating a second signal; and a processing component coupled to the first physical microphone and the second physical microphone, the processing component forming a first virtual microphone, 25 the processing component forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone, the processing component forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal, and summing the first intermediate signal and the 30 second signal, the processing component detecting acoustic voice activity of a speaker when an energy ratio of energies of the first virtual 52 WO 2010/048635 PCT/US2009/062129 microphone and the second virtual microphone is greater than a threshold value. The device of an embodiment comprises applying a calibration to at least one of the first signal and the second signal. 5 The calibration of an embodiment compensates a second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. The device of an embodiment comprises applying a delay to the first intermediate signal. 10 The delay of an embodiment is proportional to a time difference between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone. The forming of the first virtual microphone of an embodiment comprises applying the filter to the second signal. 15 The forming of the first virtual microphone of an embodiment comprises applying the calibration to the second signal. The forming of the first virtual microphone of an embodiment comprises applying the delay to the first signal. The forming of the first virtual microphone by the combining of an 20 embodiment comprises subtracting the second signal from the first signal. The filter of an embodiment is an adaptive filter. The device of an embodiment comprises adapting the filter to minimize a second virtual microphone output when only speech is being received by the first physical microphone and the second physical 25 microphone. The adapting of an embodiment comprises applying a least-mean squares process. The device of an embodiment comprises generating coefficients of the filter during a period when only speech is being received by the first 30 physical microphone and the second physical microphone. The forming of the filter of an embodiment comprises generating a first quantity by applying a calibration to the second signal. The forming of 53 WO 2010/048635 PCT/US2009/062129 the filter of an embodiment comprises generating a second quantity by applying the delay to the first signal. The forming of the filter of an embodiment comprises forming the filter as a ratio of the first quantity to the second quantity. 5 The generating of the energy ratio of an embodiment comprises generating the energy ratio for a frequency band. The generating of the energy ratio of an embodiment comprises generating the energy ratio for a frequency subband. The frequency subband of an embodiment includes frequencies 10 higher than approximately 200 Hertz (Hz). The frequency subband of an embodiment includes frequencies in a range from approximately 250 Hz to 1250 Hz. The frequency subband of an embodiment includes frequencies in a range from approximately 200 Hz to 3000 Hz. 15 The filter of an embodiment is a static filter. The forming of the filter of an embodiment comprises determining a first distance as distance between the first physical microphone and a mouth of the speaker. The forming of the filter of an embodiment comprises determining a second distance as distance between the second 20 physical microphone and the mouth. The forming of the filter of an embodiment comprises forming a ratio of the first distance to the second distance. The device of an embodiment comprises generating a vector of the energy ratio versus time. 25 The first virtual microphone and the second virtual microphone of an embodiment are distinct virtual directional microphones. The first virtual microphone and the second virtual microphone of an embodiment have approximately similar responses to noise. The first virtual microphone and the second virtual microphone of an 30 embodiment have approximately dissimilar responses to speech. The first and second physical microphones of an embodiment are omnidirectional microphones. 54 WO 2010/048635 PCT/US2009/062129 The device of an embodiment comprises positioning the first physical microphone and the second physical microphone along an axis and separating the first physical microphone and the second physical microphone by a first distance. 5 A midpoint of the axis of an embodiment is a second distance from a mouth of the speaker, wherein the mouth is located in a direction defined by an angle relative to the midpoint. Embodiments described herein include a device comprising: a headset including at least one loudspeaker, wherein the headset attaches 10 to a region of a human head; a microphone array connected to the headset, the microphone array including a first physical microphone outputting a first signal and a second physical microphone outputting a second signal; and a processing component coupled to the first physical microphone and the second physical microphone, the processing 15 component forming a first virtual microphone, the processing component forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone, the processing component forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal, and summing the 20 first intermediate signal and the second signal, the processing component detecting acoustic voice activity of a speaker when an energy ratio of energies of the first virtual microphone and the second virtual microphone is greater than a threshold value. The AVAD can be a component of a single system, multiple systems, and/or 25 geographically separate systems. The AVAD can also be a subcomponent or subsystem of a single system, multiple systems, and/or geographically separate systems. The AVAD can be coupled to one or more other components (not shown) of a host system or a system coupled to the host system. One or more components of the AVAD and/or a corresponding system or 30 application to which the AVAD is coupled or connected includes and/or runs under and/or in association with a processing system. The processing system includes any collection of processor-based devices or computing devices operating together, or components of processing systems or devices, as is known in the art. 55 WO 2010/048635 PCT/US2009/062129 For example, the processing system can include one or more of a portable computer, portable communication device operating in a communication network, and/or a network server. The portable computer can be any of a number and/or combination of devices selected from among personal computers, cellular 5 telephones, personal digital assistants, portable computing devices, and portable communication devices, but is not so limited. The processing system can include components within a larger computer system. Aspects of the AVAD and corresponding systems and methods described herein may be implemented as functionality programmed into any of a variety of 10 circuitry, including programmable logic devices (PLDs), such as field programmable gate arrays (FPGAs), programmable array logic (PAL) devices, electrically programmable logic and memory devices and standard cell-based devices, as well as application specific integrated circuits (ASICs). Some other possibilities for implementing aspects of the AVAD and corresponding systems and 15 methods include: microcontrollers with memory (such as electronically erasable programmable read only memory (EEPROM)), embedded microprocessors, firmware, software, etc. Furthermore, aspects of the AVAD and corresponding systems and methods may be embodied in microprocessors having software based circuit emulation, discrete logic (sequential and combinatorial), custom 20 devices, fuzzy (neural) logic, quantum devices, and hybrids of any of the above device types. Of course the underlying device technologies may be provided in a variety of component types, e.g., metal-oxide semiconductor field-effect transistor (MOSFET) technologies like complementary metal-oxide semiconductor (CMOS), bipolar technologies like emitter-coupled logic (ECL), polymer 25 technologies (e.g., silicon-conjugated polymer and metal-conjugated polymer metal structures), mixed analog and digital, etc. It should be noted that any system, method, and/or other components disclosed herein may be described using computer aided design tools and expressed (or represented), as data and/or instructions embodied in various 30 computer-readable media, in terms of their behavioral, register transfer, logic component, transistor, layout geometries, and/or other characteristics. Computer-readable media in which such formatted data and/or instructions may be embodied include, but are not limited to, non-volatile storage media in various forms (e.g., optical, magnetic or semiconductor storage media) and carrier waves 35 that may be used to transfer such formatted data and/or instructions through 56 WO 2010/048635 PCT/US2009/062129 wireless, optical, or wired signaling media or any combination thereof. Examples of transfers of such formatted data and/or instructions by carrier waves include, but are not limited to, transfers (uploads, downloads, e-mail, etc.) over the Internet and/or other computer networks via one or more data transfer protocols 5 (e.g., HTTP, FTP, SMTP, etc.). When received within a computer system via one or more computer-readable media, such data and/or instruction-based expressions of the above described components may be processed by a processing entity (e.g., one or more processors) within the computer system in conjunction with execution of one or more other computer programs. 10 Unless the context clearly requires otherwise, throughout the description and the claims, the words "comprise," "comprising," and the like are to be construed in an inclusive sense as opposed to an exclusive or exhaustive sense; that is to say, in a sense of "including, but not limited to." Words using the singular or plural number also include the plural or singular number respectively. 15 Additionally, the words "herein," hereunderr," "above," "below," and words of similar import, when used in this application, refer to this application as a whole and not to any particular portions of this application. When the word "or" is used in reference to a list of two or more items, that word covers all of the following interpretations of the word: any of the items in the list, all of the items in the list 20 and any combination of the items in the list. The above description of embodiments of the AVAD and corresponding systems and methods is not intended to be exhaustive or to limit the systems and methods to the precise forms disclosed. While specific embodiments of, and examples for, the AVAD and corresponding systems and methods are described 25 herein for illustrative purposes, various equivalent modifications are possible within the scope of the systems and methods, as those skilled in the relevant art will recognize. The teachings of the AVAD and corresponding systems and methods provided herein can be applied to other systems and methods, not only for the systems and methods described above. 30 The elements and acts of the various embodiments described above can be combined to provide further embodiments. These and other changes can be made to the AVAD and corresponding systems and methods in light of the above detailed description. In general, in the following claims, the terms used should not be construed 35 to limit the AVAD and corresponding systems and methods to the specific 57 WO 2010/048635 PCT/US2009/062129 embodiments disclosed in the specification and the claims, but should be construed to include all systems that operate under the claims. Accordingly, the AVAD and corresponding systems and methods is not limited by the disclosure, but instead the scope is to be determined entirely by the claims. 5 While certain aspects of the AVAD and corresponding systems and methods are presented below in certain claim forms, the inventors contemplate the various aspects of the AVAD and corresponding systems and methods in any number of claim forms. Accordingly, the inventors reserve the right to add additional claims after filing the application to pursue such additional claim forms for other aspects 10 of the AVAD and corresponding systems and methods. 58
Claims (86)
1. A method comprising: forming a first virtual microphone by combining a first signal of a first 5 physical microphone and a second signal of a second physical microphone; forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone; forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal, and summing the first 10 intermediate signal and the second signal; generating an energy ratio of energies of the first virtual microphone and the second virtual microphone; and detecting acoustic voice activity of a speaker when the energy ratio is greater than a threshold value. 15
2. The method of claim 1, wherein the first virtual microphone and the second virtual microphone are distinct virtual directional microphones.
3. The method of claim 2, wherein the first virtual microphone and the 20 second virtual microphone have approximately similar responses to noise.
4. The method of claim 3, wherein the first virtual microphone and the second virtual microphone have approximately dissimilar responses to speech. 25
5. The method of claim 1, comprising applying a calibration to at least one of the first signal and the second signal.
6. The method of claim 5, wherein the calibration compensates a 30 second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. 59 WO 2010/048635 PCT/US2009/062129
7. The method of claim 5, comprising applying a delay to the first intermediate signal.
8. The method of claim 7, wherein the delay is proportional to a time 5 difference between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone.
9. The method of claim 8, wherein the forming of the first virtual microphone comprises applying the filter to the second signal. 10
10. The method of claim 9, wherein the forming of the first virtual microphone comprises applying the calibration to the second signal.
11. The method of claim 10, wherein the forming of the first virtual 15 microphone comprises applying the delay to the first signal.
12. The method of claim 11, wherein the forming of the first virtual microphone by the combining comprises subtracting the second signal from the first signal. 20
13. The method of claim 12, wherein the filter is an adaptive filter.
14. The method of claim 13, comprising adapting the filter to minimize a second virtual microphone output when only speech is being received by 25 the first physical microphone and the second physical microphone.
15. The method of claim 13, wherein the adapting comprises applying a least-mean squares process. 30
16. The method of claim 13, comprising generating coefficients of the filter during a period when only speech is being received by the first physical microphone and the second physical microphone. 60 WO 2010/048635 PCT/US2009/062129
17. The method of claim 13, wherein the forming of the filter comprises: generating a first quantity by applying a calibration to the second signal; 5 generating a second quantity by applying the delay to the first signal; forming the filter as a ratio of the first quantity to the second quantity. 10
18. The method of claim 17, wherein the generating of the energy ratio comprises generating the energy ratio for a frequency band.
19. The method of claim 17, wherein the generating of the energy ratio comprises generating the energy ratio for a frequency subband. 15
20. The method of claim 19, wherein the frequency subband includes frequencies higher than approximately 200 Hertz (Hz).
21. The method of claim 19, wherein the frequency subband includes 20 frequencies in a range from approximately 250 Hz to 1250 Hz.
22. The method of claim 19, wherein the frequency subband includes frequencies in a range from approximately 200 Hz to 3000 Hz. 25
23. The method of claim 12, wherein the filter is a static filter.
24. The method of claim 23, wherein the forming of the filter comprises: determining a first distance as distance between the first physical microphone and a mouth of the speaker; 30 determining a second distance as distance between the second physical microphone and the mouth; and forming a ratio of the first distance to the second distance. 61 WO 2010/048635 PCT/US2009/062129
25. The method of claim 1, comprising generating a vector of the energy ratio versus time. 5
26. The method of claim 1, wherein the first and second physical microphones are omnidirectional microphones.
27. The method of claim 1, comprising positioning the first physical microphone and the second physical microphone along an axis and 10 separating the first physical microphone and the second physical microphone by a first distance.
28. The method of claim 27, wherein a midpoint of the axis is a second distance from a mouth of the speaker, wherein the mouth is located in a 15 direction defined by an angle relative to the midpoint.
29. A method comprising: forming a first virtual microphone; forming a filter by generating a first quantity by applying a 20 calibration to a second signal of a second physical microphone, generating a second quantity by applying the delay to a first signal of a first physical microphone, and forming the filter as a ratio of the first quantity to the second quantity; forming a second virtual microphone by applying the filter to the first 25 signal to generate a first intermediate signal, and summing the first intermediate signal and the second signal; and generating a ratio of energies of the first virtual microphone and the second virtual microphone and detecting acoustic voice activity using the ratio. 30 62 WO 2010/048635 PCT/US2009/062129
30. The method of claim 29, wherein the first virtual microphone and the second virtual microphone have approximately similar responses to noise and approximately dissimilar responses to speech. 5
31. The method of claim 29, comprising applying a calibration to at least one of the first signal and the second signal, wherein the calibration compensates a second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. 10
32. The method of claim 29, comprising applying a delay to the first intermediate signal, wherein the delay is proportional to a time difference between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone. 15
33. The method of claim 29, wherein the forming of the first virtual microphone comprises applying the filter to the second signal.
34. The method of claim 33, wherein the forming of the first virtual 20 microphone comprises applying the calibration to the second signal.
35. The method of claim 34, wherein the forming of the first virtual microphone comprises applying the delay to the first signal. 25
36. The method of claim 35, wherein the forming of the first virtual microphone by the combining comprises subtracting the second signal from the first signal.
37. The method of claim 29, wherein the filter is an adaptive filter. 30 63 WO 2010/048635 PCT/US2009/062129
38. The method of claim 29, comprising adapting the filter to minimize a second virtual microphone output when only speech is being received by the first physical microphone and the second physical microphone. 5
39. The method of claim 37, wherein the adapting comprises applying a least-mean squares process.
40. The method of claim 37, comprising generating coefficients of the filter during a period when only speech is being received by the first 10 physical microphone and the second physical microphone.
41. The method of claim 29, wherein the generating of the ratio comprises generating the ratio for a frequency band. 15
42. The method of claim 29, wherein the generating of the ratio comprises generating the ratio for a frequency subband.
43. The method of claim 29, comprising generating a vector of the ratio versus time. 20
44. A method comprising: forming a first virtual microphone by generating a first combination of a first signal and a second signal, wherein the first signal is received from a first physical microphone and the second signal is received from a 25 second physical microphone; forming a filter by generating a first quantity by applying a calibration to at least one of the first signal and the second signal, generating a second quantity by applying a delay to the first signal, and forming the filter as a ratio of the first quantity to the second quantity; and 30 forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal and summing the first intermediate signal and the second signal; and 64 WO 2010/048635 PCT/US2009/062129 determining a presence of acoustic voice activity of a speaker when an energy ratio of energies of the first virtual microphone and the second virtual microphone is greater than a threshold value. 5
45. An acoustic voice activity detection system comprising: a first virtual microphone comprising a first combination of a first signal and a second signal, wherein the first signal is received from a first physical microphone and the second signal is received from a second physical microphone; 10 a filter, wherein the filter is formed by generating a first quantity by applying a calibration to at least one of the first signal and the second signal, generating a second quantity by applying a delay to the first signal, and forming the filter as a ratio of the first quantity to the second quantity; and 15 a second virtual microphone formed by applying the filter to the first signal to generate a first intermediate signal and summing the first intermediate signal and the second signal, wherein acoustic voice activity of a speaker is determined to be present when an energy ratio of energies of the first virtual microphone and the second virtual microphone is greater 20 than a threshold value.
46. The system of claim 45, wherein the first virtual microphone and the second virtual microphone have approximately similar responses to noise and approximately dissimilar responses to speech. 25
47. The system of claim 45, wherein a calibration is applied to the second signal, wherein the calibration compensates a second response of the second physical microphone so that the second response is equivalent to a first response of the first physical microphone. 30
48. The system of claim 45, wherein the delay is applied to the first intermediate signal, wherein the delay is proportional to a time difference 65 WO 2010/048635 PCT/US2009/062129 between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone.
49. The system of claim 45, wherein the first virtual microphone is 5 formed by applying the filter to the second signal.
50. The system of claim 49, wherein the first virtual microphone is formed by applying the calibration to the second signal. 10
51. The system of claim 50, wherein the first virtual microphone is formed by applying the delay to the first signal.
52. The system of claim 51, wherein the first virtual microphone is formed by subtracting the second signal from the first signal. 15
53. The system of claim 45, wherein the filter is an adaptive filter.
54. The system of claim 45, wherein the filter is adapted to minimize a second virtual microphone output when only speech is being received by 20 the first physical microphone and the second physical microphone.
55. The system of claim 45, wherein coefficients of the filter are generated during a period when only speech is being received by the first physical microphone and the second physical microphone. 25
56. The system of claim 45, wherein the energy ratio comprises an energy ratio for a frequency band.
57. The system of claim 45, wherein the energy ratio comprises an 30 energy ratio for a frequency subband.
58. A device comprising: 66 WO 2010/048635 PCT/US2009/062129 a first physical microphone generating a first signal; a second physical microphone generating a second signal; and a processing component coupled to the first physical microphone and the second physical microphone, the processing component forming a first 5 virtual microphone, the processing component forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone, the processing component forming a second virtual microphone by applying the filter to the first signal to generate a first intermediate signal, and summing the first intermediate 10 signal and the second signal, the processing component detecting acoustic voice activity of a speaker when an energy ratio of energies of the first virtual microphone and the second virtual microphone is greater than a threshold value. 15
59. The device of claim 58, comprising applying a calibration to at least one of the first signal and the second signal.
60. The device of claim 59, wherein the calibration compensates a second response of the second physical microphone so that the second 20 response is equivalent to a first response of the first physical microphone.
61. The device of claim 59, comprising applying a delay to the first intermediate signal. 25
62. The device of claim 61, wherein the delay is proportional to a time difference between arrival of the speech at the second physical microphone and arrival of the speech at the first physical microphone.
63. The device of claim 62, wherein the forming of the first virtual 30 microphone comprises applying the filter to the second signal. 67 WO 2010/048635 PCT/US2009/062129
64. The device of claim 63, wherein the forming of the first virtual microphone comprises applying the calibration to the second signal.
65. The device of claim 64, wherein the forming of the first virtual 5 microphone comprises applying the delay to the first signal.
66. The device of claim 65, wherein the forming of the first virtual microphone by the combining comprises subtracting the second signal from the first signal. 10
67. The device of claim 66, wherein the filter is an adaptive filter.
68. The device of claim 67, comprising adapting the filter to minimize a second virtual microphone output when only speech is being received by 15 the first physical microphone and the second physical microphone.
69. The device of claim 67, wherein the adapting comprises applying a least-mean squares process. 20
70. The device of claim 67, comprising generating coefficients of the filter during a period when only speech is being received by the first physical microphone and the second physical microphone.
71. The device of claim 67, wherein the forming of the filter comprises: 25 generating a first quantity by applying a calibration to the second signal; generating a second quantity by applying the delay to the first signal; forming the filter as a ratio of the first quantity to the second 30 quantity. 68 WO 2010/048635 PCT/US2009/062129
72. The device of claim 71, wherein the generating of the energy ratio comprises generating the energy ratio for a frequency band.
73. The device of claim 71, wherein the generating of the energy ratio 5 comprises generating the energy ratio for a frequency subband.
74. The device of claim 73, wherein the frequency subband includes frequencies higher than approximately 200 Hertz (Hz). 10
75. The device of claim 73, wherein the frequency subband includes frequencies in a range from approximately 250 Hz to 1250 Hz.
76. The device of claim 73, wherein the frequency subband includes frequencies in a range from approximately 200 Hz to 3000 Hz. 15
77. The device of claim 66, wherein the filter is a static filter.
78. The device of claim 77, wherein the forming of the filter comprises: determining a first distance as distance between the first physical 20 microphone and a mouth of the speaker; determining a second distance as distance between the second physical microphone and the mouth; and forming a ratio of the first distance to the second distance. 25
79. The device of claim 58, comprising generating a vector of the energy ratio versus time.
80. The device of claim 58, wherein the first virtual microphone and the second virtual microphone are distinct virtual directional microphones. 30
81. The device of claim 80, wherein the first virtual microphone and the second virtual microphone have approximately similar responses to noise. 69 WO 2010/048635 PCT/US2009/062129
82. The device of claim 81, wherein the first virtual microphone and the second virtual microphone have approximately dissimilar responses to speech. 5
83. The device of claim 58, wherein the first and second physical microphones are omnidirectional microphones.
84. The device of claim 58, comprising positioning the first physical 10 microphone and the second physical microphone along an axis and separating the first physical microphone and the second physical microphone by a first distance.
85. The device of claim 84, wherein a midpoint of the axis is a second 15 distance from a mouth of the speaker, wherein the mouth is located in a direction defined by an angle relative to the midpoint.
86. A device comprising: a headset including at least one loudspeaker, wherein the headset 20 attaches to a region of a human head; a microphone array connected to the headset, the microphone array including a first physical microphone outputting a first signal and a second physical microphone outputting a second signal; and a processing component coupled to the first physical microphone and 25 the second physical microphone, the processing component forming a first virtual microphone, the processing component forming a filter that describes a relationship for speech between the first physical microphone and the second physical microphone, the processing component forming a second virtual microphone by applying the filter to the first signal to 30 generate a first intermediate signal, and summing the first intermediate signal and the second signal, the processing component detecting acoustic voice activity of a speaker when an energy ratio of energies of the first 70 WO 2010/048635 PCT/US2009/062129 virtual microphone and the second virtual microphone is greater than a threshold value. 71
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2016202314A AU2016202314A1 (en) | 2008-10-24 | 2016-04-13 | Acoustic Voice Activity Detection (AVAD) for electronic systems |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10842608P | 2008-10-24 | 2008-10-24 | |
US61/108,426 | 2008-10-24 | ||
PCT/US2009/062129 WO2010048635A1 (en) | 2008-10-24 | 2009-10-26 | Acoustic voice activity detection (avad) for electronic systems |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
AU2016202314A Division AU2016202314A1 (en) | 2008-10-24 | 2016-04-13 | Acoustic Voice Activity Detection (AVAD) for electronic systems |
Publications (1)
Publication Number | Publication Date |
---|---|
AU2009308442A1 true AU2009308442A1 (en) | 2010-04-29 |
Family
ID=42119724
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
AU2009308442A Abandoned AU2009308442A1 (en) | 2008-10-24 | 2009-10-26 | Acoustic Voice Activity Detection (AVAD) for electronic systems |
AU2016202314A Abandoned AU2016202314A1 (en) | 2008-10-24 | 2016-04-13 | Acoustic Voice Activity Detection (AVAD) for electronic systems |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
AU2016202314A Abandoned AU2016202314A1 (en) | 2008-10-24 | 2016-04-13 | Acoustic Voice Activity Detection (AVAD) for electronic systems |
Country Status (5)
Country | Link |
---|---|
EP (1) | EP2353302A4 (en) |
CN (1) | CN102282865A (en) |
AU (2) | AU2009308442A1 (en) |
CA (1) | CA2741652A1 (en) |
WO (1) | WO2010048635A1 (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8019091B2 (en) | 2000-07-19 | 2011-09-13 | Aliphcom, Inc. | Voice activity detector (VAD) -based multiple-microphone acoustic noise suppression |
US9066186B2 (en) | 2003-01-30 | 2015-06-23 | Aliphcom | Light-based detection for acoustic applications |
US9099094B2 (en) | 2003-03-27 | 2015-08-04 | Aliphcom | Microphone array with rear venting |
WO2011140110A1 (en) * | 2010-05-03 | 2011-11-10 | Aliphcom, Inc. | Wind suppression/replacement component for use with electronic systems |
GB2499781A (en) * | 2012-02-16 | 2013-09-04 | Ian Vince Mcloughlin | Acoustic information used to determine a user's mouth state which leads to operation of a voice activity detector |
TWI474317B (en) | 2012-07-06 | 2015-02-21 | Realtek Semiconductor Corp | Signal processing apparatus and signal processing method |
CN106454673B (en) * | 2016-09-05 | 2019-01-22 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | Microphone array output signal adaptive calibration method based on RLS algorithm |
EP3291226B1 (en) | 2016-09-05 | 2020-11-04 | Unify Patente GmbH & Co. KG | A method of treating speech data, a device for handling telephone calls and a hearing device |
CN107889002B (en) * | 2017-10-30 | 2019-08-27 | 恒玄科技(上海)有限公司 | Neck ring bluetooth headset, the noise reduction system of neck ring bluetooth headset and noise-reduction method |
DE102019102414B4 (en) * | 2019-01-31 | 2022-01-20 | Harmann Becker Automotive Systems Gmbh | Method and system for detecting fricatives in speech signals |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7123727B2 (en) * | 2001-07-18 | 2006-10-17 | Agere Systems Inc. | Adaptive close-talking differential microphone array |
US8098844B2 (en) * | 2002-02-05 | 2012-01-17 | Mh Acoustics, Llc | Dual-microphone spatial noise suppression |
KR101434071B1 (en) * | 2002-03-27 | 2014-08-26 | 앨리프컴 | Microphone and voice activity detection (vad) configurations for use with communication systems |
DE10223544C1 (en) * | 2002-05-27 | 2003-07-24 | Siemens Audiologische Technik | Amplifier device for hearing aid with microphone and pick-up coil inputs, has amplifier provided with separate filters for acoustic and inductive feedback compensation |
ATE339757T1 (en) * | 2003-06-17 | 2006-10-15 | Sony Ericsson Mobile Comm Ab | METHOD AND DEVICE FOR VOICE ACTIVITY DETECTION |
US7916849B2 (en) * | 2004-06-02 | 2011-03-29 | Clearone Communications, Inc. | Systems and methods for managing the gating of microphones in a multi-pod conference system |
US7706549B2 (en) * | 2006-09-14 | 2010-04-27 | Fortemedia, Inc. | Broadside small array microphone beamforming apparatus |
WO2008157421A1 (en) * | 2007-06-13 | 2008-12-24 | Aliphcom, Inc. | Dual omnidirectional microphone array |
-
2009
- 2009-10-26 CN CN2009801515125A patent/CN102282865A/en active Pending
- 2009-10-26 AU AU2009308442A patent/AU2009308442A1/en not_active Abandoned
- 2009-10-26 WO PCT/US2009/062129 patent/WO2010048635A1/en active Application Filing
- 2009-10-26 CA CA2741652A patent/CA2741652A1/en not_active Abandoned
- 2009-10-26 EP EP09822855.4A patent/EP2353302A4/en not_active Withdrawn
-
2016
- 2016-04-13 AU AU2016202314A patent/AU2016202314A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
EP2353302A1 (en) | 2011-08-10 |
WO2010048635A1 (en) | 2010-04-29 |
CA2741652A1 (en) | 2010-04-29 |
AU2016202314A1 (en) | 2016-05-05 |
EP2353302A4 (en) | 2016-06-08 |
CN102282865A (en) | 2011-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8321213B2 (en) | Acoustic voice activity detection (AVAD) for electronic systems | |
US8326611B2 (en) | Acoustic voice activity detection (AVAD) for electronic systems | |
US11818534B2 (en) | Forming virtual microphone arrays using dual omnidirectional microphone array (DOMA) | |
US9263062B2 (en) | Vibration sensor and acoustic voice activity detection systems (VADS) for use with electronic systems | |
US20140126743A1 (en) | Acoustic voice activity detection (avad) for electronic systems | |
US8488803B2 (en) | Wind suppression/replacement component for use with electronic systems | |
US8452023B2 (en) | Wind suppression/replacement component for use with electronic systems | |
US10230346B2 (en) | Acoustic voice activity detection | |
AU2016202314A1 (en) | Acoustic Voice Activity Detection (AVAD) for electronic systems | |
US11627413B2 (en) | Acoustic voice activity detection (AVAD) for electronic systems | |
CA2798282A1 (en) | Wind suppression/replacement component for use with electronic systems | |
US20230379621A1 (en) | Acoustic voice activity detection (avad) for electronic systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MK5 | Application lapsed section 142(2)(e) - patent request and compl. specification not accepted |