# Patent

## Title
Method for characterising emitters by the association of parameters related to the same radio emitter

## Abstract

Method of using one or more transmitters and/or one or more parameters associated with a transmitter by using a reception station comprising a device suitable for measuring over time a set of K parameters dependent on the transmitters associated with vectors {circumflex over (η)}k representative of the transmitters for 1≦k≦K. The method includes a step of extracting the parameter or parameters consisting in grouping together by transmitter the parameters which are associated therewith by means of a technique of independent component analysis.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
The present Application is based on International Application No. PCT/EP2005/055541, filed on Oct. 25, 2005, which in turn corresponds to French Application No. 0411448, filed on Oct. 27, 2004, and priority is hereby claimed under 35 USC §119 based on these applications. Each of these applications are hereby incorporated by reference in their entirety into the present application.
FIELD OF THE INVENTION
The invention relates notably to a passive method of locating fixed or mobile transmitters on the ground. The objective is notably to determine the position of one or more transmitters on the ground on the basis of a mobile reception system. It also applies in respect of configurations where the transmitters are at altitude and the reception system on the ground.
BACKGROUND OF THE INVENTION
The associating of parameters also makes it possible to characterize a transmitter by its signal-to-noise ratio. It is also possible to characterize an FH (Frequency Hopping) by its plateau durations, its frequencies of appearance and its direction of arrival. A modulation can also be characterized by identifying for example its amplitude state and phase state.
FIG. 1 illustrates airborne location with a mobile reception system, the transmitter 1 is at the position (x0, y0, z0), the carrier 2 at the instant tk is at the position (xk, yk, zk) and sees the transmitter 1 at the angle of incidence (θ(tk, x0, y0, z0), Δ(tk, x0, y0, z0)). The angles θ(t, x0, y0, z0) and Δ(t, x0, y0, z0) evolve over time and depend on the position of the transmitter and also the trajectory of the reception system. The angles θ(t, x0, y0, z0) and Δ(t, x0, y0, z0) are charted with the aid of an array of N antennas that can be fixed under the carrier as shown by FIG. 2.
The antennas Ai of the array receive the sources with a phase and amplitude that depend on the angle of incidence of the sources and also the position of the antennas.
Antenna processing techniques generally utilize the spatial diversity of the sources (or transmitters): use of the spatial position of the antennas of the array to better utilize the differences in incidence and in distance of the sources. Antenna processing breaks down into two main areas of activity:
1—Spatial filtering, illustrated in FIG. 3, the objective of which is to extract either the modulated signals sm(t), or the symbols contained in the signal (Demodulation). This filtering consists in combining the signals received on the array of sensors so as to form an optimal reception antenna for one of the sources. Spatial filtering can be blind or cooperative.

It is cooperative when there is a priori knowledge about the transmitted signals (directions of arrival, symbol sequences, etc.) and it is blind in the contrary case. It is used for blind source separation, filtering matched to direction of arrival (beam formation) or to replicas, multi-sensor MODEM (demodulation), etc.

2—The objective of estimating the parameters of the transmitters is to determine various parameters such as: their Doppler frequencies, their bit rates, their modulation indices, their positions (xm,ym), their incidences (θm,Δm) and their direction vectors a(θm,Δm) (response of the array of sensors to a source with direction (θm,Δm)) etc.

For example, goniometry and blind identification procedures exist in this area:
The objective of goniometry is to determine the incidences θ(t, xm, ym, zm) in 1D or the pair of incidences (θ(t, xm, ym, zm),Δ(t, xm, ym, zm)) in 2D. For this purpose, goniometry algorithms use the observations arising from the antennas or sensors. When the waves from all the transmitters propagate in the same plane, it suffices to apply a 1D goniometry, in other cases, a 2D goniometry.
The objective of blind identification procedures (ICA) is notably to determine the direction vectors a(θm,Δm) of each of the transmitters.
The known location techniques according to the prior art generally use histogram techniques to group the parameters together. However, these techniques have the drawback of requiring a priori knowledge about the standard deviation of the parameters in order to fix the stepsize of the histogram.
SUMMARY OF THE INVENTION
The invention relies on a new approach consisting notably in judiciously associating the parameters related to one and the same transmitter.
The invention relates to a method for characterizing one or more transmitters and/or one or more parameters associated with a transmitter by using a reception station comprising a device adapted to measuring over time a set of K parameters dependent on the transmitters associated with vectors {circumflex over (η)}k representative of the transmitters for 1≦k≦K characterized in that it comprises at least one step of extracting the parameter or parameters consisting in grouping together by transmitter the parameters which are associated therewith by means of a technique of independent component analysis.
The method according to the invention has notably the following advantages:

    
    
        It does not require any parameter of settings and no a priori knowledge about the statistics of the parameters,
        it makes it possible to count the number of transmitters on the basis of the number of dominant values of a covariance matrix of the parameters ({circumflex over (R)}xx): use of the techniques for detecting the number of sources of the goniometry,
        it makes it possible to identify as many parameter vectors (ηm) as desired,
        it makes it possible to use the step of associating the parameters of different nature such as the incidence of the sources, their signal-to-noise ratio or else their direction vectors,
        it makes it possible to determine the mean incidence of each incident transmitter on the basis of the measured incidences,
        to perform the mean location of each incident transmitter on the basis of the measured direction vectors or the measured incidences,
        to extract the phase states of a modulation on the basis of the real and imaginary parts of the signal of a linear modulation,
        to separate the Frequency Hopping, FH, signals by measuring the plateau durations and incidences: an FH is characterized by a single incidence and a single plateau.
    
    


Still other objects and advantages of the present invention will become readily apparent to those skilled in the art from the following detailed description, wherein the preferred embodiments of the invention are shown and described, simply by way of illustration of the best mode contemplated of carrying out the invention. As will be realized, the invention is capable of other and different embodiments, and its several details are capable of modifications in various obvious aspects, all without departing from the invention. Accordingly, the drawings and description thereof are to be regarded as illustrative in nature, and not as restrictive.



BRIEF DESCRIPTION OF THE DRAWINGS
The present invention is illustrated by way of example, and not by limitation, in the figures of the accompanying drawings, wherein elements having the same reference numeral designations represent like elements throughout and wherein:
FIG. 1 an example of airborne location with a mobile reception system,
FIG. 2 an array of 5 antennas,
FIG. 3 a diagram of spatial filtering by beam formation in one direction,
FIGS. 4, 5 and 6 a numerical example of the use of the method according to the invention.



DETAILED DESCRIPTION OF THE INVENTION
The following example is given in conjunction with FIG. 1, comprising a transmitter 1 to be located by using an aircraft 2 equipped with devices allowing the measurement of parameters associated with the transmitters, and a processor suitable for executing the steps according to the invention.
In the presence of M transmitters, a location system measures, for example over time, a set of K parameters (representative of the transmitting sources) characterized by the vectors ηk for 1≦k≦K. The vectors ηk are for example composed of the azimuth θk and of the signal-to-noise ratio SNRk of one of the transmitters at the instant tk:ηk=[θk SNRk]T where (T) denotes the transpose of a vector.
This vector can also be composed of the direction vector a(θk) of one of the sources and of its signal-to-noise ratio: ηk=[a(θk)T SNRk]T.
More generally the kth measurement {circumflex over (η)}k is marred by an error and is associated with the mth transmitter in the following manner:

{circumflex over (η)}k•=ηm+ek for 1≦k≦K and 1≦m≦M  (1)

where ek is the noise vector associated with the kth measurement and ηm the parameter vector associated with the mth transmitter.

The invention consists notably in extracting the M vectors ηm associated with a transmitter m in the midst of the K measurements {circumflex over (η)}k.

Procedure for Identifying the Parameter Vectors of the Principal Sources

By virtue of the location system with which the aircraft is equipped or more generally a system performing parameter measurements, K measurements {circumflex over (η)}k are available. The prime objective of the invention is to identify the M vectors ηm associated with the M incident transmitters.
For this purpose, the method comprises a first step consisting in transforming the measured vectors {circumflex over (η)}k into vectors f({circumflex over (η)}k) of larger dimension.
For a goniometry system in azimuth θk where the vector {circumflex over (η)}k is equal to [SNRk]T (azimuth and signal-to-noise ratio), the transformation step consists in performing the following bijective transformation:




  
    
      
        
          f
          ⁡
          
            (
            
              
                
                  η
                  ^
                
                k
              
              =
              
                [
                
                  
                    
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                  
                  
                    
                      
                        SNR
                        k
                      
                    
                  
                
                ]
              
            
            )
          
        
        =
        
          [
          
            
              
                1
              
            
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      j
                      ⁢
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                    )
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      L
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                    )
                  
                
              
            
            
              
                
                  SNR
                  k
                
              
            
          
          ]
        
      
    
    
      
        (
        2
        )
      
    
  




For a goniometry system in azimuth θk and in elevation Δk where the vector representative of the set of parameters K measured for the M transmitters {circumflex over (η)}k•=[θkΔkSNRk]T the transformation step consists in performing the following bijective transformation:




  
    
      
        
          f
          ⁡
          
            (
            
              
                
                  η
                  ^
                
                k
              
              =
              
                [
                
                  
                    
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                  
                  
                    
                      
                        
                          Δ
                          ^
                        
                        k
                      
                    
                  
                  
                    
                      
                        SNR
                        k
                      
                    
                  
                
                ]
              
            
            )
          
        
        =
        
          [
          
            
              
                1
              
            
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      j
                      ⁢
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                    )
                  
                
              
            
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      j
                      ⁢
                      
                        
                          Δ
                          ^
                        
                        k
                      
                    
                    )
                  
                
              
            
            
              
                
                  
                    exp
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            θ
                            ^
                          
                          k
                        
                      
                      )
                    
                  
                  ⁢
                  
                    exp
                    ⁡
                    
                      (
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            Δ
                            ^
                          
                          k
                        
                      
                      )
                    
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  SNR
                  k
                
              
            
          
          ]
        
      
    
    
      
        (
        3
        )
      
    
  




For a goniometry system in the presence of FH signals whose incidences {circumflex over (θ)}k and plateau durations {circumflex over (T)}k have been measured, the vector {circumflex over (η)}k may be written: {circumflex over (η)}k•=•[{circumflex over (θ)}k{circumflex over (T)}k]T. The transformation step for {circumflex over (η)}k consists in performing the following bijective transformation:




  
    
      
        
          f
          ⁡
          
            (
            
              
                
                  η
                  ^
                
                k
              
              =
              
                [
                
                  
                    
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                  
                  
                    
                      
                        
                          T
                          ^
                        
                        k
                      
                    
                  
                
                ]
              
            
            )
          
        
        =
        
          [
          
            
              
                1
              
            
            
              
                
                  
                    θ
                    ^
                  
                  k
                
              
            
            
              
                
                  
                    T
                    ^
                  
                  k
                
              
            
          
          ]
        
      
    
    
      
        (
        4
        )
      
    
  




For a system seeking to extract the phase states of a transmitter from the signal x(kT) of the BPSK, the vector {circumflex over (η)}k may be written:




  
    
      
        
          
            η
            ^
          
          k
        
        =
        
          [
          
            
              
                
                                    ⁢
                  
                    (
                    
                      x
                      ⁡
                      
                        (
                        
                          k
                          ⁢
                          
                            
                          
                          ⁢
                          T
                        
                        )
                      
                    
                    )
                  
                
              
            
            
              
                
                  J
                  ⁡
                  
                    (
                    
                      x
                      ⁡
                      
                        (
                        
                          k
                          ⁢
                          
                            
                          
                          ⁢
                          T
                        
                        )
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        5
        )
      
    
  




Where (z) and I(z) denote the real and imaginary parts of the complex z and T is the symbol rate. In the case of a BPSK transmitter with 2 phase states, M=2 states are present, such that:





  
    
      
        
          η
          1
        
        =
        
          
            
              [
              
                
                  
                    
                      cos
                      ⁡
                      
                        (
                        φ
                        )
                      
                    
                  
                
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        φ
                        )
                      
                    
                  
                
              
              ]
            
            ⁢
            
              
            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            
              η
              2
            
          
          =
          
            [
            
              
                
                  
                    -
                    
                      cos
                      ⁡
                      
                        (
                        φ
                        )
                      
                    
                  
                
              
              
                
                  
                    -
                    
                      sin
                      ⁡
                      
                        (
                        φ
                        )
                      
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        6
        )
      
    
  




The determination of the vectors η1 and η2 will make it possible to deduce the phase rotation φ of the BPSK. In this case it is possible to construct the following vector f({circumflex over (η)}k):




  
    
      
        
          f
          ⁡
          
            (
            
              
                η
                ^
              
              k
            
            )
          
        
        =
        
          [
          
            
              
                1
              
            
            
              
                
                                    ⁢
                  
                    (
                    
                      x
                      ⁡
                      
                        (
                        
                          k
                          ⁢
                          
                            
                          
                          ⁢
                          T
                        
                        )
                      
                    
                    )
                  
                
              
            
            
              
                
                  J
                  ⁡
                  
                    (
                    
                      x
                      ⁡
                      
                        (
                        
                          k
                          ⁢
                          
                            
                          
                          ⁢
                          T
                        
                        )
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        7
        )
      
    
  




The length of the vector f({circumflex over (η)}k) determines the maximum number of identifiable transmitters. In an application of goniometric type that estimates the incidences θ or (θ, Δ), it is possible to say for example that this maximum number will not exceed the number of sensors of the array that enabled the goniometry.
On the basis of the K≧dim(•f({circumflex over (η)}k)) vectors f({circumflex over (η)}k) the method thereafter calculates the following covariance matrix:




  
    
      
        
          
            R
            ^
          
          xx
        
        =
        
          
            1
            K
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
                
                  [
                  
                    
                      f
                      ⁡
                      
                        (
                        
                          
                            η
                            ^
                          
                          k
                        
                        )
                      
                    
                    ⊗
                    
                      f
                      ⁡
                      
                        (
                        
                          
                            η
                            ^
                          
                          k
                        
                        )
                      
                    
                  
                  ]
                
                ⁡
                
                  [
                  
                    
                      f
                      ⁡
                      
                        (
                        
                          
                            η
                            ^
                          
                          k
                        
                        )
                      
                    
                    ⊗
                    
                      f
                      ⁡
                      
                        (
                        
                          
                            η
                            ^
                          
                          k
                        
                        )
                      
                    
                  
                  ]
                
              
              H
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where  denotes the Kronecker product such that uv=[u(1)vTu(2)vT . . . ] and (•)H the conjugate transpose.

This matrix {circumflex over (R)}xx may also be written:





  
    
      
        
          
            R
            ^
          
          xx
        
        =
        
          
            
              ∑
              
                m
                =
                1
              
              M
            
            ⁢
            
              r
              ⁢
              
                
              
              ⁢
              
                
                  
                    
                      ρ
                      m
                    
                    ⁡
                    
                      [
                      
                        
                          f
                          ⁡
                          
                            (
                            
                              η
                              m
                            
                            )
                          
                        
                        ⊗
                        
                          f
                          ⁡
                          
                            (
                            
                              η
                              m
                            
                            )
                          
                        
                      
                      ]
                    
                  
                  ⁡
                  
                    [
                    
                      
                        f
                        ⁡
                        
                          (
                          
                            η
                            m
                          
                          )
                        
                      
                      ⊗
                      
                        f
                        ⁡
                        
                          (
                          
                            η
                            m
                          
                          )
                        
                      
                    
                    ]
                  
                
                H
              
            
          
          +
          
            R
            bb
          
        
      
    
    
      
        (
        9
        )
      
    
  




where Rbb is a noise matrix and ρm is the number of vectors f({circumflex over (η)}k) associated with the vector f(ηm).

The matrix {circumflex over (R)}xx is the covariance matrix of the K observations f({circumflex over (η)}k)({circumflex over (η)}k) According to equation (9) it reduces to the covariance matrix of the signatures f(ηm)(ηm) of the M transmitters. Knowing that the signatures f(ηm)(ηm) are all different since the transmitters are associated with different parameters, the principal components of the matrix {circumflex over (R)}xx (eigenvectors associated with the M largest eigenvalues) define the same space as the M signatures f(ηm)(ηm) of the transmitters: {circumflex over (R)}xx is completely related to the vector space of the M signatures f(ηm)(ηm) of the transmitters.
The rank of the matrix {circumflex over (R)}xx is thus equal to the number of transmitters M. This rank can be determined from the eigenvalues of this matrix.
Thus in the presence of a vector f({circumflex over (η)}k) of dimension N, the matrix {circumflex over (R)}xx is of dimension N2×N2 and it is then possible to identify at most N2 transmitters.
The method thereafter comprises a step of identifying the transformed vectors f(ηm) on the basis of {circumflex over (R)}xx so as thereafter to deduce therefrom the parameter vectors ηm of each of the transmitters.
For this purpose the first operation consists in decomposing the matrix {circumflex over (R)}xx into eigenelements to obtain its eigenvalues.
On the basis of the eigenvalues of the matrix, it is possible to determine the number of sources M by applying, for example, the procedure described in reference [4] or any other “enumeration” procedure which makes it possible to count the number of principal components of the matrix {circumflex over (R)}xx. This number in the given example is related to the number M of transmitters.
On the basis of the M eigenelements associated with the largest eigenvalues λm it is possible to determine the square root of the matrix {circumflex over (R)}xx:

{circumflex over (R)}xx1/2=EsΛs1/2=BΩ1/2UH such that B=[f(η1)(η1) . . . f(ηM)(ηM)] and Ω=diag{ρ1 . . . ρM}  (10)

where diag{ . . . } is a diagonal matrix composed of the elements of { . . . }, Es and Λs=diag{(λ1 . . . λM} are composed respectively of the eigenvectors and eigenvalues of {circumflex over (R)}xx associated with the M largest eigenvalues: λ1≧ . . . ≧λM.

The columns of the matrix B are composed of the signatures f(ηm)(ηm) of each of the transmitters.

The matrix U is unitary (UHU=IM where IM is the identity matrix of dimension N×N). Knowing that the columns of the square root {circumflex over (R)}xx1/2 are in the same space as the columns f(ηm)(ηm) of the matrix B, the matrix U is a change of basis matrix. The matrix U is moreover unitary since its columns are mutually orthogonal vectors. Subsequently in the description the method will use this orthogonality property to identify the matrix U. For the identification of U the method will moreover use the redundant structure of B which is related to the Kronecker product .
The determination of U is done for example by utilizing the redundant structure of the matrix B, i.e.:




  
    
      
        
          B
          =
          
            [
            
              
                
                  
                    A
                    ⁢
                    
                      
                    
                    ⁢
                    
                      Φ
                      1
                    
                  
                
              
              
                
                  ⋮
                
              
              
                
                  
                    A
                    ⁢
                    
                      
                    
                    ⁢
                    
                      Φ
                      N
                    
                  
                
              
            
            ]
          
        
        ⁢
        
          
        
        ⁢
        
          

        
        ⁢
        with
        ⁢
        
          
            
          
          ⁢
          
            
          
        
        ⁢
        
          

        
        ⁢
        
          
        
        ⁢
        
          A
          =
          
            [
            
              
                f
                ⁡
                
                  (
                  
                    η
                    1
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                f
                ⁡
                
                  (
                  
                    η
                    M
                  
                  )
                
              
            
            ]
          
        
        ⁢
        
          
        
        ⁢
        
          

        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        
          

        
        ⁢
        
          
        
        ⁢
        
          
            Φ
            n
          
          =
          
            diag
            ⁢
            
              {
              
                
                  
                    f
                    n
                  
                  ⁡
                  
                    (
                    
                      η
                      1
                    
                    )
                  
                
                ⁢
                
                  
                
                ⁢
                ...
                ⁢
                
                  
                
                ⁢
                
                  
                    f
                    n
                  
                  ⁡
                  
                    (
                    
                      η
                      M
                    
                    )
                  
                
              
              }
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




where fn(ηm) is the nth component of the vector f(ηm) of dimension N×1.

Under these conditions the matrix {circumflex over (R)}xx1/2 is composed of N sub-blocks Γn such that:





  
    
      
        
          
            R
            ^
          
          xx
          
            1
            /
            2
          
        
        =
        
          
            [
            
              
                
                  
                    Γ
                    1
                  
                
              
              
                
                  ⋮
                
              
              
                
                  
                    Γ
                    N
                  
                
              
            
            ]
          
          =
          
            
              B
              ⁢
              
                
              
              ⁢
              
                U
                H
              
              ⁢
              with
              ⁢
              
                
              
              ⁢
              
                Γ
                n
              
            
            =
            
              A
              ⁢
              
                
              
              ⁢
              
                Φ
                n
              
              ⁢
              
                U
                H
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




The columns of each matrix Γn are in the same vector space as the sought-after signatures f(ηm) of each of the M transmitters. The matrices Γn differ by change of basis matrices which are equal to a diagonal matrix to within the sought-after unitary matrix U. These properties of the matrices Γn depend on the redundant structure of the matrix B.

Accordingly the following matrices Ψij:

Ψij=Γi#Γj=UΦi−1ΦjUH  (13)

all have the matrix U as matrix of eigenvectors (# denotes the pseudo inverse such that Γ#=(ΓHΓ)−1ΓH).

The method uses the unitary character of the matrix U to identify it: The columns of U are orthogonal vectors.

Under these conditions, to determine the unitary matrix U, the joint diagonalization of the JADE procedure described for example in reference [5] of the following matrices or of any other procedure known to the person skilled in the art is performed:

Ψij for (1≦i≦N and j<i)  (14)

Once the matrix U has been estimated, it is possible to deduce therefrom the matrix B to within an amplitude by performing according to relation (10):

{circumflex over (B)}=BΩ1/2={circumflex over (R)}xx1/2U=[f(η1)(η1)√{square root over (ρ1)} . . . f(ηM)(ηM)√{square root over (ρM)}  (15)

Knowing that the mth column of the matrix {circumflex over (B)} may be written bm=[bm1T . . . bmNT]T=f(ηm)(ηm)√{square root over (ρm)} it is transformed into the following matrix Bm:

Bm=[bm1 . . . bmN]=f(ηm)f(ηm)T√{square root over (ρm)}  (16)

where the bmi are vectors of dimensions N×1.
Knowing that the 1st component of f(ηm) is equal to 1, f(ηm) is deduced from Bm by taking the singular vector em of Bm associated with the largest singular value and by performing f(ηm)=em/em (1) where em (1) is the first component of the vector em. The method performs this normalization since the vectors f(ηk) are always constructed with a first component equal to 1. These steps of constructing f(ηk) and normalizing em make it possible to resolve the phase ambiguity of the singular vectors em.
The M principal vectors ηm are deduced from the transformed vectors f(ηm) since f(•) is bijective.
Having identified the M principal vectors ηm, the method compiles for example statistics on the components of each of the vectors. For an application of location type this step makes it possible notably to give in addition to the mean position of the transmitter, an error bracket for the estimation of the position.
For example for the vector ηm=.[θm SNRm]T (1≦m≦M), the method determines the statistics of the azimuth θm (bias and standard deviation) so as to give the value of the azimuth within a bracket.
The first step consists in determining the set Φm of vectors {circumflex over (η)}k associated with the transmitter having mean vector ηm.
Knowing that:




  
    
      
        
          
            
              
                (
                
                  
                    F
                    H
                  
                  ⁢
                  F
                
                )
              
              
                -
                1
              
            
            ⁢
            
              F
              H
            
            ⁢
            
              f
              ⁡
              
                (
                
                  η
                  m
                
                )
              
            
          
          =
          
            
              δ
              m
            
            =
            
              [
              
                
                  
                    
                      0
                      
                        m
                        -
                        1
                      
                    
                  
                
                
                  
                    1
                  
                
                
                  
                    
                      0
                      
                        M
                        -
                        m
                      
                    
                  
                
              
              ]
            
          
        
        ⁢
        
          

        
        ⁢
        with
        ⁢
        
          

        
        ⁢
        
          F
          =
          
            [
            
              
                f
                ⁡
                
                  (
                  
                    η
                    1
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                f
                ⁡
                
                  (
                  
                    η
                    M
                  
                  )
                
              
            
            ]
          
        
      
    
    
      
        (
        17
        )
      
    
  




where 0L is a zero vector of dimension L×1.

Thus all the components of δm are zero with the exception of the mth which equals 1. It should be noted that the filtering matrix (FHF)−1FH is a separator of the transmitters: By applying this filter to a signature f(ηm) of the mth transmitter only the component associated with this transmitter is nonzero.
Thus to determine the set Φm to which the vector {circumflex over (η)}k belongs, the method uses the property of equation (17) (where the objective of the matrix (FHF)−1FH is to separate the transmitters) by calculating the vector βk of equation (14) which ought to be close to δm when f({circumflex over (η)}k) is associated with the mth source.

(FHF)−1FHf({circumflex over (η)}k)=βk  (18)

The vector {circumflex over (η)}k belongs to the set Φm of the mth source if the mth component βk(m) of largest modulus satisfies: |βk(m)|>α.

The threshold α is chosen close to 1 (typically α=0.9).

The steps of the method of constructing the sets Φm in the presence of K vectors {circumflex over (η)}k comprise for example the following steps:

Step R.1 k=1 and initialization of the M sets Φm to Ø (empty set),

Step R.2 Calculation of the vector βk using equation (18).

Step R.3 Search for the component βk(imax) such that: |βk(imax).>βk(i) for i≠ imax.

Step R.4 If |βk(imax)|>α then Φimax={Φimax{circumflex over (η)}k},

Step R.5 k←k+1,

Step R.6 If k<K return to step R.2.

Once the M sets Φm={{circumflex over (η)}k close to ηm} have been determined, the method calculates a statistic of the components of the vector ηm for example as Mean Square Error or MSE.
The Mean Square Error (MSE) of the ith component of ηm may be written:




  
    
      
        
          
            
              
                MSE
                m
              
              ⁡
              
                (
                i
                )
              
            
            2
          
          =
          
            
              1
              
                card
                ⁡
                
                  (
                  
                    Φ
                    m
                  
                  )
                
              
            
            ⁢
            
              
                ∑
                
                  
                    
                      η
                      ^
                    
                    k
                  
                  ∈
                  
                    Φ
                    m
                  
                
              
              ⁢
              
                
                  
                    (
                    
                      
                        
                          
                            η
                            ^
                          
                          k
                        
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                      -
                      
                        
                          mean
                          m
                        
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                    
                    )
                  
                  2
                
                ⁢
                
                  
                
                ⁢
                with
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              mean
              m
            
            ⁡
            
              (
              i
              )
            
          
          =
          
            
              1
              
                card
                ⁡
                
                  (
                  
                    Φ
                    m
                  
                  )
                
              
            
            ⁢
            
              
                ∑
                
                  
                    
                      η
                      ^
                    
                    k
                  
                  ∈
                  
                    Φ
                    m
                  
                
              
              ⁢
              
                
                  
                    η
                    ^
                  
                  k
                
                ⁡
                
                  (
                  i
                  )
                
              
            
          
        
      
    
    
      
        (
        19
        )
      
    
  




where card (Φm) is the cardinal of the set Φm and meanm(i) is the mean value that has to be close to ηm(i).

In the example of FIGS. 4 and 5 the vector {circumflex over (η)}k=[θkTSNRk]T and the function f(•) satisfy:




  
    
      
        
          f
          ⁡
          
            (
            
              
                
                  η
                  ^
                
                k
              
              =
              
                [
                
                  
                    
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                  
                  
                    
                      
                        SNR
                        k
                      
                    
                  
                
                ]
              
            
            )
          
        
        =
        
          [
          
            
              
                1
              
            
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      j
                      ⁢
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                    )
                  
                
              
            
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      2
                      ⁢
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          θ
                          ^
                        
                        k
                      
                    
                    )
                  
                
              
            
            
              
                
                  SNR
                  k
                
              
            
          
          ]
        
      
    
    
      
        (
        20
        )
      
    
  




The following FIG. 4 shows the distribution of the goniometry plots in the (θk, SNRk) space:
and FIG. 5, the evolution of these plots over time:
The method gives for the M=2 sources:


    
    
        θ1=66.06° and SNR2=27.78 dB
        θ2=77.77° and SNR2=28.17 dB
    
    


In FIGS. 4 and 5, the mean values estimated by the method appear as a solid line.
In the example of FIG. 6 the vector {circumflex over (η)}k=[a(θk)T SNRk]T and the function f(•) satisfy:




  
    
      
        
          f
          ⁡
          
            (
            
              
                η
                ^
              
              k
            
            )
          
        
        =
        
          [
          
            
              
                1
              
            
            
              
                
                  a
                  ⁡
                  
                    (
                    
                      
                        θ
                        ^
                      
                      k
                    
                    )
                  
                
              
            
            
              
                
                  
                    SNR
                    k
                  
                  /
                  
                    (
                    
                      
                        max
                        ⁡
                        
                          (
                          
                            SNR
                            k
                          
                          )
                        
                      
                      -
                      
                        min
                        ⁡
                        
                          (
                          
                            SNR
                            k
                          
                          )
                        
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        21
        )
      
    
  




In FIG. 6 the coefficients ck=1Ha({circumflex over (θ)}k)| are represented dotted and the coefficients cm=|1H a(θm)| as a solid line, where the direction vector a(sm) has been deduced from the M vectors sm estimated by the method.
The vector 1 is composed of 1 s.
The method has detected M=3 categories of sources.
FIG. 6 shows that two of the categories are permanently present while the last is present in a much more sporadic manner.
These examples show that the method is applied in an independent manner to the type of parameters of the sources.
Without departing from the scope of the invention, the method can be applied in respect of directions of arrival θm, direction vectors a(θm) or else signal-to-noise ratios SNRm.
It will be readily seen by one of ordinary skill in the art that the present invention fulfils all of the objects set forth above. After reading the foregoing specification, one of ordinary skill in the art will be able to affect various changes, substitutions of equivalents and various aspects of the invention as broadly disclosed herein. It is therefore intended that the protection granted hereon be limited only by definition contained in the appended claims and equivalents thereof.
BIBLIOGRAPHY


    [1] L. ALBERA, A. FERREOL, P. CHEVALIER and P. COMON. GRETSI 2003, Paris, September 2003, “ICAR, un algorithme de ICA à convergence rapide, robuste au bruit” [ICAR, a fast convergence, noise-robust ICA algorithm].
    [2] L. ALBERA, A. FERREOL and P. CHEVALIER. ICA2003, Nara (Japan), April 2003, Sixth order blind identification of undetermined mixtures (SIRBI) of sources.
    [3] P. COMON, Signal Processing, Elsevier, April 1994, vol. 36”, No. 3, pp 287-314, Independent Component Analysis, a new concept ̃?.
    [4] O. MICHEL, P. LARZABAL and H. CLERGEOT Test de détection du nombre de sources corrélées pour les méthodes HR en traitement d'antenne [Test for detecting the number of correlated sources for HR procedures in antenna processing]. GRETSI 91 in Juans les Pins.
    [5] J. F. CARDOSO, A. SOULOUMIAC, IEE Proceedings-F, Vol. 140, No. 6, pp. 362-370, Dec. 1993. Blind beamforming for non-gaussian signals.

## Claims

1. A method of using one or more transmitters and/or one or more parameters associated with a transmitter by using a reception station having a device suitable for measuring over time a set of K parameters dependent on the transmitting associated with vectors {circumflex over (η)}k representative of the transmitters for 1<k<K comprising the steps of:
extracting the parameter or parameters consisting in grouping together by transmitter the parameters which are associated therewith by means of a technique of independent component analysis, wherein the step of associating the parameters for each transmitter M comprises at least the following steps:

Step No. 1 transforming the vectors {circumflex over (η)}k representative of the set of the K parameters for the M transmitters into vectors f({circumflex over (n)}k) by using a bijective function {tilde over (f)}( ) and where the 1st component of f({circumflex over (η)}k) being equal to 1,
Step No. 2 determining a covariance matrix {circumflex over (R)}xx the vectors f({circumflex over (n)}k) and decomposing it into eigenelements, so as to obtain its eigenvalues, where {circumflex over (R)}xx is a matrix related to the vector space of the M signatures f(ηm){circle around (×)}f(ηm) of the transmitters,
Step No. 3 calculating the rank M of the matrix {circumflex over (R)}xx its eigenvalues found in step 2,
Step No. 4 calculating a square root of {circumflex over (R)}xx its M dominant eigenelements: {circumflex over (R)}xx1/2=EsΛs1/2 where Es corresponds to the eigenvectors and Λs to the M largest eigenvalues,
Step No. 5 deducing from {circumflex over (R)}xx1/2 the matrices Γn according to




  
    
      R
      ^
    
    xx
    
      1
      /
      2
    
  
  =
  
    
      [
      
        
          
            
              Γ
              1
            
          
        
        
          
            ⋮
          
        
        
          
            
              Γ
              N
            
          
        
      
      ]
    
    =
    
      B
      ⁢
      
        
      
      ⁢
      
        U
        H
      
    
  




with Γn=AΦnUH and from the Γn calculating Ψij according to Ψij=Γi#Γj=UΨi−1ΨjUH for (1≦i≦N and j<i), where the columns of each matrix Γn define the same space vector as the signatures f(ηm) of the M transmitters,
Step No. 6 identifying unitary matrix U by joint diagonalization of the Ψij for (1≦i≦N and j<i),
Step No. 7 determining the vectors f(ηm) from the columns bm of the matrix {circumflex over (R)}xx1/2 U where U is the unitary matrix: Transformation of column bm into the matrix Bm according to Bm=[bm1 . . . bmN]=f(ηm)f(ηm)T√{square root over (ρm)} and extraction of this matrix from the singular vector em associated with the largest singular value so as to perform f(ηm)=em/em(1) with em(1) the component of the vector em, and
Step No. 8 applying the inverse transform of the function {tilde over (f)}( ) so as to deduce therefrom the vectors ηm.

2. The method as claimed in claim 1 comprising at least one step of evaluating the statistics of the components of the vectors ηm found.

3. The method as claimed in claim 2 comprising a step of evaluating the mean and the standard deviation of the incidences of each of the transmitters.

4. The method as claimed in claim 2 comprising at least the following steps:
determining the M sets Φm of vectors {circumflex over (η)}k associated with the principal vector ηm; and
determining the statistics of the components of the vector ηm from the set Φm by using a Mean Square Error procedure.