Methods and arrangements for signal reception and parameter estimation

• US 4,965,732 A
• Filed: 11/02/1987
• Issued: 10/23/1990
• Est. Priority Date: 11/06/1985
• Status: Expired due to Term

• Current Assignee: Board of Trustees of the Leland Stanford Junior University (Stanford Management Co.)
• Original Assignee: Board of Trustees of the Leland Stanford Junior University (Stanford Management Co.)
• Inventors: Paulraj, Arogyaswami J., Kailath, Thomas, Roy, Richard H. III
• Primary Examiner(s): Smith, Jerry
• Assistant Examiner(s): Baker, Stephen M.
• Application Number:
  US07/116,094
• Time in Patent Office: 1,086 Days
• Field of Search: 364/807, 364/724.16, 364/456, 364/460, 364/516, 364/517, 342/147, 342/153, 342/195, 367/124, 367/129
• US Class Current: 342/147
• CPC Class Codes: G01S 3/14 Systems for determining direction or deviation from predetermined direction aerial arrangements for changing or varying the orientation or the shape of the directional pattern H01Q3/00; combinations of different interacting aerial units for giving a desired directional characteristic H01Q21/29; aerials or aerial systems providing at least two radiation patterns H01Q25/00
• CPC Class Codes: G01S 3/74 Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals systems in which the directions of different signals are determined sequentially and displayed simultaneously G01S3/04, G01S3/14
• CPC Class Codes: G01S 3/8006 Multi-channel systems specially adapted for direction-finding, i.e. having a single aerial system capable of giving simultaneous indications of the directions of different signals
• CPC Class Codes: G01S 3/802 Systems for determining direction or deviation from predetermined direction
• CPC Class Codes: G06F 2218/22 Source localisation; Inverse modelling

Abstract

The invention described herein relates generally to the field of signal processing for signal reception and parameter estimation. The invention has many applications such as frequency estimation and filtering, and array data processing, etc. For convenience, only applications of this invention to sensor array processing are described herein. The array processing problem addressed is that of signal parameter and waveform estimation utilizing data collected by an array of sensors. Unique to this invention is that the sensor array geometry and individual sensor characteristics need not be known. Also, the invention provides substantial advantages in computations and storage over prior methods. However, the sensors must occur in pairs such that the paired elements are identical except for a displacement which is the same for all pairs. These element pairs define two subarrays which are identical except for a fixed known displacement. The signals must also have a particular structure which in direction-of-arrival estimation applications manifests itself in the requirement that the wavefronts impinging on the sensor array be planar. Once the number of signals and their parameters are estimated, the array configurations can be determined and the signals individually extracted. The invention is applicable in the context of array data processing to a number of areas including cellular mobile communications, space antennas, sonobuoys, towed arrays of acoustic sensors, and structural analysis.

Claims

1. A method of detecting multiple signal sources and estimating parameters thereof comprising the following steps:

(a) providing an array of at least one group of a plurality of signal sensor pairs, the sensors in each pair being identical and the displacement between sensors of each pair in a group being equal, thereby defining two subarrays (X and Y),

(b) obtaining signal measurements with the sensor array so configured,

(c) processing said signal measurements from said two subarrays (X and Y) to identify the number of sources and estimate parameters thereof, including identifying eigenvalues from which source number and parameter estimates are based,

(d) solving the signal copy problem and determining array response (direction) vectors using the generalized eigenvectors, and

(e) estimating the array geometry from the said array response vectors.

2. The method as defined in claim 1 and further including a variation to improve numerical characteristics using generalized singular value decompositions of data matrices by:

(a) forming data matrices X and Y from the data from the subarrays,

(b) computing the generalized singular vectors of the matrix pair (X, Y) yielding X=U_x Σ

_x V* and Y=U_x Σ

_x V*,

(c) then computing the eigenvalues of Σ

_x^-1 U*U_y Σ

_y and locating those which lie on or near the unit circle, the number of which corresponding to the number of sources and the locations of which corresponding to the parameter estimates.

3. The method as defined by claim 1 wherein said step of identifying eigenvalues utilizes a total least-squares algorithm.

4. The method as defined by claim 3 wherein said step of identifying eigenvalues includes

obtaining an estimate of R_ZZ, denoted R_ZZ, from the measurements available,

computing the generalized eigen-decomposition

5. A methodd of locating signal sources and estimating source parameters comprising the following steps:

(a) providing an array of at least one group of a plurality of signal sensor pairs, the sensors in each pair being identical and the displacement between sensors of each pair in a group being equal, thereby defining two subarrays (X and Y),

(b) obtaining signal measurements with the sensor array so configured, and

(c) processing said signal measurements from said two subarrays (X and Y) to identify the number of sources and estimate parameters thereof, including a generalized singular value decomposition of data matrices comprising

forming the matrix Z from the available measurements,

computing the generalized singular value decomposition (GSVD) of ##EQU26## obtaining the signal subspace estimate S_Z =span E_Z where ##EQU27## computing the singular value decomposition (SVD) of ##EQU28## calculating the eigenvalues of

space="preserve" listing-type="equation">Φ

.sub.k =λ

.sub.k (-V.sub.XY V.sub.YY.sup.-1),

and

estimating the signal parameters θ

_k =f^-1 (Φ

_k).

6. For use in locating signal sources and estimating source parameters, apparatus for measuring signals from said sources comprising

an array of at least one group of a plurality of signal sensor pairs for generating signals, the sensors in each pair being identical and the displacement between sensors of each pair in a group being equal, thereby defining two subarrays (X and Y), and

signal processing means for processing said signals from said two subarrays (X and Y) to identify the number of sources and estimate parameters thereof, wherein said signal processing means

obtains an estimate of R_ZZ, denoted R_ZZ, from the measurements available,

computes the generalized eigen-decomposition

View Dependent Claims (7)

7. Apparatus as defined by claim 6 wherein said signal processing means

forms the matrix Z from the available measurements,

computes the GSVD of ##EQU32## obtains the signal subspace estimate S_Z =R{E_Z }, ##EQU33## computes the SVD of {[E_X |E_Y ]Σ

.sub.ψ

^1/2 }=UΣ

V*, ##EQU34## calculates the eigenvalues of

space="preserve" listing-type="equation">Φ

.sub.k =λ

.sub.k (-V.sub.XY V.sub.YY.sup.-1),

and

estimates the signal parameters θ

_k =f^-1 (Φ

_k).