# Patent

## Title
Apparatus and method for estimating the direction of arrival of a source signal using a microphone array

## Abstract

An apparatus for estimating the direction of arrival (“DOA”) of a source signal includes microphones for transducing the source signal, amplifiers in signal communication with a corresponding one of the microphones, analog-to-digital converters in signal communication with a corresponding one of the amplifiers, and a digital processor in signal communication with the analog-to-digital converters for estimating the direction of arrival of the source signal; where a corresponding method for estimating the DOA of the source signal at one of several microphones includes estimating the measured signal spectral covariance matrix (“Rx”), and estimating the angle of the DOA responsive to the measured signal spectral covariance matrix.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of U.S. Provisional Application Ser. No. 60/386,760, filed Jun. 5, 2002 and entitled “Apparatus and Method for Estimating the Direction Of Arrival of a Source Signal using a Microphone Array”, which is incorporated herein by reference in its entirety.


BACKGROUND
The present disclosure relates to the estimation of the direction of arrival of a source signal using a microphone array. Applications may include Adaptive Signal Processing schemes for Hearing Aids, car kits, mobile communications, voice controlled devices, and the like.
There are various estimators known for estimating the Direction Of Arrival (“DOA”) of a given source. With respect to a fixed system of coordinates, the DOA is defined as the set of angles that define the location of the source, up to the intrinsic symmetry of the array. Thus, for the 3-microphone array shown in FIG. 1, the only angle that can be estimated is theta, the angle between the incidence line and the microphone line. With such an array, any source located on a hyperboloid symmetric with respect to the microphone axis will correspond to the same delay between microphone signals. Such procedures are commonly known to those of ordinary skill in the pertinent art.
SUMMARY
These and other drawbacks and disadvantages of the prior art are addressed by an apparatus and method for estimating the direction of arrival of a source signal using a microphone array.
An apparatus for estimating the direction of arrival (“DOA”) of a source signal includes microphones for transducing the source signal, amplifiers in signal communication with a corresponding one of the microphones, analog-to-digital converters in signal communication with a corresponding one of the amplifiers, and a digital processor in signal communication with the analog-to-digital converters for estimating the direction of arrival of the source signal.
A corresponding method for estimating the DOA of the source signal at one of several microphones includes estimating the measured signal spectral covariance matrix (“Rx”), and estimating the angle of the DOA responsive to the measured signal spectral covariance matrix. These and other aspects, features and advantages of the present disclosure will become apparent from the following description of exemplary embodiments, which is to be read in connection with the accompanying drawings.
These and other aspects, features and advantages of the present disclosure will become apparent from the following description of exemplary embodiments, which is to be read in connection with the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
The present disclosure teaches an apparatus and method for estimating the direction of arrival of a source signal using a microphone array in accordance with the following exemplary figures, in which:
FIG. 1 shows a schematic geometric diagram of a 3-Microphone Array and the Direction Of Arrival (“DOA”);
FIG. 2 shows a schematic block diagram of a DOA Estimator according to an illustrative embodiment of the present disclosure; and
FIG. 3 shows a schematic hardware diagram of a DOA Estimator according to an illustrative embodiment of the present disclosure.



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
The present disclosure presents an array signal processing apparatus and method for estimating the direction of arrival (“DOA”) of a source signal in a noisy environment using a multi-sensor array. Embodiments of this disclosure may be used as a component block in general adaptive algorithms or schemes.
Possible applications include: Adaptive Signal Processing schemes for Hearing Aids, car kits, mobile communication, voice controlled devices, and the like. The present disclosure provides a method and apparatus to estimate the angle theta defining the DOA of a source. The presently disclosed approach makes novel use of statistical assumptions and optimization.
As shown in FIG. 1, a 3-Microphone Array with an incident DOA is indicated generally by the reference numeral 100. The exemplary array includes a first microphone 102, a second microphone 104 disposed a fixed distance d from the first microphone, and a third microphone 106 disposed a fixed distance d from the second microphone along a line intersecting the first and second microphones. A signal source 108 is disposed at an angle theta relative to the line of the microphones. The angle theta represents the DOA.
The mixing model is assumed as:

x1(t)=s(t)+n1(t)  (1)

. . .   (2)

xD(t)=s(t−(D−1)δτ)+nD(t)  (3)

where D is the number of microphones, δ=cos (θ) is a direct measure of the DOA, τ=dfx/c with d distance between adjacent microphones, fx is the sampling frequency, and c is the sound speed; x1(t), . . . , xD(t) are the signals as measured by the microphones, s(t) are the source signals, and n1(t), . . . , nD(t) are the noise signals.

The distance d between microphones is assumed to be very small (smaller than c/fx, so that τ≦=1 and therefore the noise signals are not independent. The following assumption is made regarding their statistics:
Assumption 1: The joint random variables (N 1 (w), . . . , N D (w)), which are the Fourier transforms of ( n 1 . . . , n D), are Gaussian distributed with zero mean and covariance matrix:




  
    
      
        
          
            R
            n
          
          ⁡
          
            (
            w
            )
          
        
        =
        
          
            
              ρ
              n
            
            ⁡
            
              (
              w
              )
            
          
          ⁡
          
            [
            
              
                
                  1
                
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                        
                        )
                      
                    
                    
                      w
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          2
                          ⁢
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                        
                        )
                      
                    
                    
                      2
                      ⁢
                      w
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                
              
              
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                        
                        )
                      
                    
                    
                      w
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                
                
                  1
                
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                        
                        )
                      
                    
                    
                      w
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                
              
              
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          2
                          ⁢
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                        
                        )
                      
                    
                    
                      2
                      ⁢
                      w
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                        
                        )
                      
                    
                    
                      w
                      ⁢
                      
                        
                      
                      ⁢
                      τ
                    
                  
                
                
                  1
                
              
            
            ]
          
        
      
    
    
      
        (
        4
        )
      
    
  




It is also assumed that:
Assumption 2: The source signal Fourier transform S (w) is Gaussian distributed with zero mean and spectral power Rs(w).
The presently disclosed method optimally estimates (in the sense of maximum likelihood) the parameters θ, Rx and ρn. The block diagram of this solution is depicted in FIG. 2.
Turning to FIG. 2, a DOA Estimator according to an illustrative embodiment of the present disclosure is indicated generally by the reference numeral 200. The estimator 200 includes an Rx block 210 for receiving an output vector x from the microphones. The block 210 is coupled in signal communication with a DOA block 212 for providing the angle theta.
Turning now to FIG. 3, a DOA Estimator according to an illustrative embodiment of the present disclosure is indicated generally by the reference numeral 300. The DOA estimator 300 includes an exemplary array of microphones comprising a first microphone 302, a second microphone 304 and a third microphone 306. The first microphone 302 is coupled in signal communication to a first amplifier 322, which, in turn, is coupled in signal communication to a first analog-to-digital converter 332. The output of the analog-to-digital converter 332 is coupled in signal communication to a digital processor 340.
Similarly, the second microphone 304 is coupled in signal communication to a second amplifier 324, which, in turn, is coupled in signal communication to a second analog-to-digital converter 334. The output of the analog-to-digital converter 334 is also coupled in signal communication to the digital processor 340. Likewise, the third microphone 306 is coupled in signal communication to a third amplifier 326, which, in turn, is coupled in signal communication to a third analog-to-digital converter 336. The output of the analog-to-digital converter 336 is also coupled in signal communication to the digital processor 340.
Thus, a physical embodiment of the scheme is presented in FIG. 3. The component blocks of FIG. 2 are now described, and a physical embodiment is presented as in FIG. 3.
Components of the DOA Estimator are as follows: The DOA estimator comprises two blocks: the first block 210 of FIG. 2 provides for estimation of Rx, and the second block 212 of FIG. 2 provides for estimation of the DOA. The microphones 302, 304 and 306 of FIG. 3 convert acoustic pressure into electrical signals, which are next amplified by the Amplifiers 322, 324 and 326, respectively, and converted into the digital domain by the A/D Converters 332, 334 and 336, respectively. The digitized measured signals are denoted by x1, . . . , xD, where D is the number of microphones. In the following exemplary embodiment, the description of the algorithms is for the case D=3, but this is by no means a restriction to the applicability of the principle. The rational for the formulas (eqns. 7–17) is also presented.
Referring to FIG. 2, the measured signal spectral covariance matrix Rx is estimated as follows:

    
    
        1. Initialize Rx(w)=0, for M values of w,
    
    






  
    w
    =
    0
  
  ,
  
    
      2
      ⁢
      π
    
    M
  
  ,
  ...
  ⁢
  
    
  
  ,
  
    
      
        2
        ⁢
        
          π
          ⁡
          
            (
            
              M
              -
              1
            
            )
          
        
      
      M
    
    ;
  




    
    
        2. Using the last M samples, x1(t), x1(t−1), . . . , x1(t−M+1), x2(t), x2(t−1), . . . , x2(t−M+1),x3(t), x3(t−1), . . . , x3(t−M+1), compute the windowed Fourier transform:
    
    






  
    
      
        
          
            
              X
              k
            
            ⁡
            
              (
              w
              )
            
          
          =
          
            
              ∑
              
                l
                =
                0
              
              
                M
                -
                1
              
            
            ⁢
            
              
                w
                ⁡
                
                  (
                  l
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    l
                  
                  )
                
              
              ⁢
              
                e
                
                  
                    -
                    i
                  
                  ⁢
                  
                    
                  
                  ⁢
                  wl
                
              
            
          
        
        ,
        
          k
          =
          1
        
        ,
        2
        ,
        3
        ,
        
          w
          =
          0
        
        ,
        
          
            2
            ⁢
            π
          
          M
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            2
            ⁢
            
              π
              ⁡
              
                (
                
                  M
                  -
                  1
                
                )
              
            
          
          M
        
      
    
    
      
        (
        5
        )
      
    
  




where w(0), . . . , w(M−1) is the window, and w is the frequency.

    
    
        3. Update the matrices Rx (w)
    
    






  
    
      
        
          
            R
            x
            updated
          
          ⁡
          
            (
            w
            )
          
        
        =
        
          
            
              (
              
                1
                -
                α
              
              )
            
            ⁢
            
              
                R
                x
                previous
              
              ⁡
              
                (
                w
                )
              
            
          
          +
          
            
              α
              ⁡
              
                [
                
                  
                    
                      
                        
                          X
                          1
                        
                        ⁡
                        
                          (
                          w
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          X
                          2
                        
                        ⁡
                        
                          (
                          w
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          X
                          3
                        
                        ⁡
                        
                          (
                          w
                          )
                        
                      
                    
                  
                
                ]
              
            
            ⁡
            
              [
              
                
                  
                    
                  
                  
                    
                      X
                      1
                    
                    ⁡
                    
                      (
                      w
                      )
                    
                  
                
                ⁢
                
                  
                    
                  
                  
                    
                      X
                      2
                    
                    ⁡
                    
                      (
                      w
                      )
                    
                  
                
                ⁢
                
                  
                    
                  
                  
                    
                      X
                      3
                    
                    ⁡
                    
                      (
                      w
                      )
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




where α is a learning rate.

Referring again to FIG. 3, the exemplary embodiment of a DOA (θ) Estimator is described with respect to the following algorithm:

    
    
        1. For each frequency w in
    
    






  
    S
    =
    
      {
      
        0
        ,
        
          
            2
            ⁢
            π
          
          M
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            2
            ⁢
            
              π
              ⁡
              
                (
                
                  M
                  -
                  1
                
                )
              
            
          
          M
        
      
      }
    
  
  ,




repeate the following steps:

    
    
        (a) Set the 3×3 matrix
    
    






  
    
      
        
          M
          ⁡
          
            (
            w
            )
          
        
        ⁡
        
          [
          
            
              
                1
              
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        w
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                      )
                    
                  
                  
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
              
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        2
                        ⁢
                        w
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                      )
                    
                  
                  
                    2
                    ⁢
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
              
            
            
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        w
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                      )
                    
                  
                  
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
              
              
                1
              
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        w
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                      )
                    
                  
                  
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
              
            
            
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        2
                        ⁢
                        w
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                      )
                    
                  
                  
                    2
                    ⁢
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
              
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        w
                        ⁢
                        
                          
                        
                        ⁢
                        τ
                      
                      )
                    
                  
                  
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                  
                
              
              
                1
              
            
          
          ]
        
      
    
    
      
        (
        7
        )
      
    
  




    
    
        where τ=dfx/c,d is the distance between two adjacent microphones, fs is the sampling frequency, and c is the speed of sound;
        (b) Solve the third order polynomial equation:

det(λM−Rx)=0  (8)

        where λ is the unknown.
        (c) Order the solution λ3≧λ2λ1 (There are three solutions because equation (8) is of order three, and, in turn, this is due to the number of microphones D.) Set

ρn=λ1  (9)

        (d) Solve the third order polynomial equation

det(μI−Rx+ρnM)=0  (10)
        
            where μ is the unknown, and I is the 3×3 identity matrix.
        
        
        (e) Order the solutions μ3≧μ2≧μ1 (note μ1=0); set Rx(w)=μ3(w).
        (f) Find a nonzero solution z=[z1 z2 z3]T of the linear system

μ3I−Rx−ρnM)z=0  (11)

        (g) Define
    
    






  
    
      
        
          
            w
            1
          
          ⁡
          
            (
            w
            )
          
        
        =
        
          
            z
            2
          
          
            z
            1
          
        
      
    
    
      
        (
        12
        )
      
    
  
  
    
      
        
          w
          2
        
        (
        
          
            
              2
              ⁢
              w
            
            +
            
              
                
                  2
                  ⁢
                  π
                
                M
              
              ⁢
              
                
              
              ⁢
              mod
              ⁢
              
                
              
              ⁢
              2
              ⁢
              π
            
          
          =
          
            
              z
              3
            
            
              z
              1
            
          
        
      
    
    
      
        (
        13
        )
      
    
  
  
    
      
        
          
            w
            2
          
          ⁡
          
            (
            
              
                2
                ⁢
                w
              
              +
              
                
                  
                    2
                    ⁢
                    π
                  
                  M
                
                ⁢
                
                  
                
                ⁢
                mod
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                π
              
            
            )
          
        
        =
        0
      
    
    
      
        (
        14
        )
      
    
  




    
    
        for wεS.
    
    






  
    
      
        
          J
          ⁡
          
            (
            δ
            )
          
        
        =
        
          real
          ⁢
          
            
          
          ⁢
          
            (
            
              
                ∑
                
                  
                    w
                    ∈
                    S
                  
                  ⁢
                  
                    
                  
                
                
                  
                
              
              ⁢
              
                
                  e
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    
                      δτ
                      /
                      N
                    
                  
                
                ⁡
                
                  (
                  
                    
                      
                        w
                        1
                      
                      ⁡
                      
                        (
                        w
                        )
                      
                    
                    +
                    
                      
                        w
                        2
                      
                      ⁡
                      
                        (
                        w
                        )
                      
                    
                  
                  )
                
              
            
            )
          
        
      
    
    
      
        (
        15
        )
      
    
  




    
    
        for δ=−N,−N+1, . . . N.
        3. Find the maximum of J for δ over {−N,−N+1, . . . , N}, say δmax,

δmax=arg maxδε{−N,−N+1, . . . ,N}J(δ)  (16)

        4. Compute the DOA θ through
    
    






  
    
      
        θ
        =
        
          arccos
          ⁡
          
            (
            
              
                δ
                max
              
              N
            
            )
          
        
      
    
    
      
        (
        17
        )
      
    
  




In operation, the exemplary embodiment of FIG. 3 is implemented with the following exemplary apparatus, including an array of 3 (in this exemplary case) closely arranged microphones 302, 304 and 306 where the distance between two adjacent microphones is around 1 cm; amplifiers 322, 324 and 326 and A/D converters 332, 334 and 336, one of each for each microphone; and a digital processor 340.
The computations presented in the previous section may be implemented very efficiently so that the total computational load is fairly small. The following features are also pertinent:

    
    
        1. The third order algebraic equations (8) and (11) have exact solutions in closed form.
        2. the equation (11) is of the form μ3+αμ2+bμ=0; thus one of the solutions is μ1=0, and the computation of the other two roots involves just a couple of operations.
    
    


Exemplary values for the parameters are summarized as follows:

    
    
        1. D=3[microphones]
        2. M=128[samples]
        3. α=0.01
        4. N=100
        5. d=1[cm]
        6. fs=20000[Hz]
        7. w is the Hamming window of size M;
    
    


It shall be understood that the above-mentioned values are merely exemplary. Thus, those of ordinary skill in the pertinent art may contemplate various other values without diminishing the generality of the algorithm or of the apparatus.
In this section we present the rationale for the expressions (7–17). The likelihood of the parameters θ, Rs(·), ρn(w) given the model (1–3) is




  
    
      
        
          L
          ⁡
          
            (
            
              θ
              ,
              
                
                  R
                  s
                
                ⁡
                
                  (
                  ·
                  )
                
              
              ,
              
                
                  ρ
                  n
                
                ⁡
                
                  (
                  ·
                  )
                
              
            
            )
          
        
        =
        
          
            ∏
            
              w
              ∈
              S
            
            
              
            
          
          ⁢
          
            
          
          ⁢
          
            
              1
              
                
                  π
                  
                    3
                    /
                    2
                  
                
                ⁢
                
                  det
                  ⁡
                  
                    (
                    
                      
                        
                          R
                          s
                        
                        ⁢
                        
                          v
                          *
                        
                        ⁢
                        v
                      
                      +
                      
                        
                          ρ
                          n
                        
                        ⁢
                        M
                      
                    
                    )
                  
                
              
            
            ⁢
            
              exp
              ⁡
              
                (
                
                  -
                  
                    Trace
                    ⁡
                    
                      (
                      
                        
                          
                            R
                            x
                          
                          ⁡
                          
                            (
                            
                              
                                
                                  R
                                  s
                                
                                ⁢
                                
                                  vv
                                  *
                                
                              
                              +
                              
                                
                                  ρ
                                  n
                                
                                ⁢
                                M
                              
                            
                            )
                          
                        
                        
                          -
                          1
                        
                      
                      )
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




where we have omitted the explicit dependency on w of Rsρn, M, v, Rx for the sake of clarity, and





  
    
      
        v
        =
        
          [
          
            1
            ⁢
            
              
            
            ⁢
            
              e
              
                
                  -
                  i
                
                ⁢
                
                  
                
                ⁢
                w
                ⁢
                
                  
                
                ⁢
                τ
                ⁢
                
                  
                
                ⁢
                
                  cos
                  ⁡
                  
                    (
                    θ
                    )
                  
                
                ⁢
                
                  
                
                ⁢
                
                  e
                  
                    
                      -
                      2
                    
                    ⁢
                    
                      
                    
                    ⁢
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    w
                    ⁢
                    
                      
                    
                    ⁢
                    τ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      cos
                      ⁡
                      
                        (
                        θ
                        )
                      
                    
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        19
        )
      
    
  




({circumflex over (θ)}, {circumflex over (R)}S(·), {circumflex over (ρ)}n(·))=argmaxθ,RS(·)L(θ,Rs(·),ρn(·))  (20)

For sufficiently many samples, the approximate solution to the above estimator can be found by using the Covariance Matching Estimation Technique (“COMET”), together with the Extended Invariance Principle as known in the art. This amounts to solving at every frequency w the problem

argminz,Rs,ρnTrace(Rx−Rxzz*−ρnM)2)  (21)

where z=[1 z2 z3] is the unknown vector, constraine to be 1 on the first component. This problem is solved, under the additional positivity constraint Rszz*+ρnM≦Rx in 8–11. Once z has been obtained at every frequency, the solution for θ is obtained through the optimization





  
    
      
        
          θ
          ^
        
        =
        
          
            arg
            ⁢
            
              
            
            ⁢
            
              
                min
                θ
              
              ⁢
              
                
                  ∑
                  
                    w
                    ∈
                    S
                  
                
                ⁢
                
                  
                    
                    
                      
                        
                          z
                          2
                        
                        ⁡
                        
                          (
                          w
                          )
                        
                      
                      -
                      
                        e
                        
                          
                            -
                            i
                          
                          ⁢
                          
                            
                          
                          ⁢
                          w
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            cos
                            ⁡
                            
                              (
                              θ
                              )
                            
                          
                        
                      
                    
                    
                  
                  2
                
              
            
          
          +
          
            
            
              
                
                  z
                  3
                
                ⁡
                
                  (
                  w
                  )
                
              
              -
              
                e
                
                  
                    -
                    2
                  
                  ⁢
                  
                    
                  
                  ⁢
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  w
                  ⁢
                  
                    
                  
                  ⁢
                  τ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    cos
                    ⁡
                    
                      (
                      θ
                      )
                    
                  
                
              
            
            
          
        
      
    
    
      
        (
        22
        )
      
    
  




These and other features and advantages of the present disclosure may be readily ascertained by one of ordinary skill in the pertinent art based on the teachings herein. It is to be understood that the teachings of the present disclosure may be implemented in various forms of hardware, software, firmware, special purpose processors, or combinations thereof.
Most preferably, the teachings of the present disclosure are implemented as a combination of hardware and software. Moreover, the software is preferably implemented as an application program tangibly embodied on a program storage unit. The application program may be uploaded to, and executed by, a machine comprising any suitable architecture. Preferably, the machine is implemented on a computer platform having hardware such as one or more central processing units (“CPU”), a random access memory (“RAM”), and input/output (“I/O”) interfaces. The computer platform may also include an operating system and microinstruction code. The various processes and functions described herein may be either part of the microinstruction code or part of the application program, or any combination thereof, which may be executed by a CPU. In addition, various other peripheral units may be connected to the computer platform such as an additional data storage unit and a printing unit.
It is to be further understood that, because some of the constituent system components and methods depicted in the accompanying drawings are preferably implemented in software, the actual connections between the system components or the process function blocks may differ depending upon the manner in which the present disclosure is programmed. Given the teachings herein, one of ordinary skill in the pertinent art will be able to contemplate these and similar implementations or configurations of the present disclosure.
Although the illustrative embodiments have been described herein with reference to the accompanying drawings, it is to be understood that the present disclosure is not limited to those precise embodiments, and that various changes and modifications may be effected therein by one of ordinary skill in the pertinent art without departing from the scope or spirit of the present disclosure. All such changes and modifications are intended to be included within the scope of the present disclosure as set forth in the appended claims.

## Claims

1. An apparatus for estimating the direction of arrival (“DOA”) of a source signal, the apparatus comprising:
a plurality of microphones for transducing the source signal;
a plurality of amplifiers, each amplifier in signal communication with a corresponding one of the plurality of microphones;
a plurality of analog-to-digital converters, each converter in signal communication with a corresponding one of the plurality of amplifiers; and
a digital processor in signal communication with the plurality of analog-to-digital converters for estimating the direction of arrival of the source signal and optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein the digital processor comprises a DOA portion for estimating the angle of the direction of arrival.

2. An apparatus as defined in claim 1 wherein the source signal is the dominant source signal of a plurality of source signals.

3. An apparatus as defined in claim 1 wherein the digital processor comprises an optimization portion for optimizing the DOA solution in the statistical sense of maximum likelihood.

4. An apparatus as defined in claim 1 wherein the digital processor comprises an Rx portion for estimating the measured signal spectral covariance matrix (“Rx”).

5. An apparatus as defined in claim 1 wherein the digital processor differentiates between any source located on a hyperboloid symmetric with respect to the microphone axis that corresponds to the same delay between microphone signals.

6. A method for estimating the direction of arrival (“DOA”) of a source signal at one of a plurality of microphones, the method comprising:
estimating the measured signal spectral covariance matrix (“Rx”);
estimating the angle of the direction of arrival (“DOA”) responsive to the measured signal spectral covariance matrix; and
optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein estimating the measured signal spectral covariance matrix (“Rx”) comprises:
receiving signal samples;
initializing the measured signal spectral covariance matrix;
computing a windowed Fourier transform responsive to the signal samples; and
updating the measured signal spectral covariance matrix responsive to the windowed Fourier transform.

7. A method as defined in claim 6, further comprising differentiating between any source located on a hyperboloid symmetric with respect to the microphone axis that corresponds to the same delay between microphone signals.

8. A method for estimating the direction of arrival (“DOA”) of a source signal at one of a plurality of microphones, the method comprising:
estimating the measured signal spectral covariance matrix (“Rx”);
estimating the angle of the direction of arrival (“DOA”) responsive to the measured signal spectral covariance matrix; and
optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein estimating the angle of the direction of arrival (“DOA”) comprises:
selecting a frequency;
forming a matrix responsive to each of the selected frequency and the distance between two adjacent microphones;
solving a polynomial equation responsive to each of the formed matrix and the measured signal spectral covariance matrix to obtain a plurality of solutions, the plurality equaling the number of microphones;
ordering the plurality of solutions;
finding a non-zero solution; and
computing the DOA in response to the ordered and non-zero solutions.

9. A method as defined in claim 8 wherein the estimated angle of the direction of arrival (“DOA”) is optimized in the statistical sense of maximum likelihood.

10. A method as defined in claim 8 wherein the estimated angle of the direction of arrival (“DOA”) corresponds to the source signal that is the dominant source signal of a plurality of source signals.

11. An apparatus for estimating the direction of arrival (“DOA”) of a sonic source signal, the apparatus comprising:
covariance means for estimating the measured signal spectral covariance matrix (“Rx”);
direction of arrival means for estimating the angle of the direction of arrival (“DOA”) responsive to the measured signal spectral covariance matrix; and
optimization means for optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein the covariance means for estimating the measured signal spectral covariance matrix (“Rx”) comprises:
microphone means for receiving signal samples;
initialization means for initializing the measured signal spectral covariance matrix;
Fourier means for computing a windowed Fourier transform responsive to the signal samples; and
update means for updating the measured signal spectral covariance matrix responsive to the windowed Fourier transform.

12. An apparatus as defined in claim 11, further comprising differentiation means for differentiating between any source located on a hyperboloid symmetric with respect to the microphone axis that corresponds to the same delay between microphone signals.

13. An apparatus for estimating the direction of arrival (“DOA”) of a sonic source signal, the apparatus comprising:
covariance means for estimating the measured signal spectral covariance matrix (“Rx”);
direction of arrival means for estimating the angle of the direction of arrival (“DOA”) responsive to the measured signal spectral covariance matrix; and
optimization means for optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein the direction of arrival means for estimating the angle of the direction of arrival (“DOA”) comprises:
selection means for selecting a frequency;
matrix means for forming a matrix responsive to each of the selected frequency and the distance between two adjacent microphones;
solution means for solving a polynomial equation responsive to each of the formed matrix and the measured signal spectral covariance matrix to obtain a plurality of solutions, the plurality equaling the number of microphones;
ordering means for ordering the plurality of solutions;
non-zero means for finding a non-zero solution; and
computation means for computing the DOA in response to the ordered and non-zero solutions.

14. An apparatus as defined in claim 13, further comprising dominant source means for estimating the angle of the direction of arrival (“DOA”) corresponding to the source signal that is the dominant source signal of a plurality of source signals.

15. An apparatus as defined in claim 13, further comprising optimization means for optimizing the estimated angle of the direction of arrival (“DOA”) in the statistical sense of maximum likelihood.

16. A program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform program steps for estimating the direction of arrival (“DOA”) of a source signal at one of a plurality of microphones, the program steps comprising:
estimating the measured signal spectral covariance matrix (“Rx”);
estimating the angle of the direction of arrival (“DOA”) responsive to the measured signal spectral covariance matrix; and
optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein the program steps for estimating the measured signal spectral covariance matrix (“Rx”) comprise:
receiving signal samples;
initializing the measured signal spectral covariance matrix;
computing a windowed Fourier transform responsive to the signal samples; and
updating the measured signal spectral covariance matrix responsive to the windowed Fourier transform.

17. A program storage device as defined in claim 16, the program steps further comprising differentiating between any source located on a hyperboloid symmetric with respect to the microphone axis that corresponds to the same delay between microphone signals.

18. A program storage device readable by machine, tangibly embodying a program of instructions executable by the machine to perform program steps for estimating the direction of arrival (“DOA”) of a source signal at one of a plurality of microphones, the program steps comprising:
estimating the measured signal spectral covariance matrix (“Rx”);
estimating the angle of the direction of arrival (“DOA”) responsive to the measured signal spectral covariance matrix; and
optimizing the DOA solution in the statistical sense of maximum likelihood,
wherein the program steps for estimating the angle of the direction of arrival (“DOA”) comprise:
selecting a frequency;
forming a matrix responsive to each of the selected frequency and the distance between two adjacent microphones;
solving a polynomial equation responsive to each of the formed matrix and the measured signal spectral covariance matrix to obtain a plurality of solutions, the plurality equaling the number of microphones;
ordering the plurality of solutions;
finding a non-zero solution; and
computing the DOA in response to the ordered and non-zero solutions.

19. A program storage device as defined in claim 18 wherein the estimated angle of the direction of arrival (“DOA”) corresponds to the source signal that is the dominant source signal of a plurality of source signals.

20. A program storage device as defined in claim 18 wherein the estimated angle of the direction of arrival (“DOA”) is optimized in the statistical sense of maximum likelihood.