# Patent

## Title
Method for demodulating a signal

## Abstract

A method for demodulating a signal in a receiver including at least two antennas, each receiving a signal, the received signals corresponding to a single emitted signal including symbol frames in which certain symbols, referred to as driver symbols, are known to the receiver, the method including: noise whitening in order to form two combined signals, the noise components of which are separate; normalizing the noise components of the combined signals in order to form two signals, the noise components of which are separate and have equal average standards; and performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, the noise components of which are separate and which have equal average standards, wherein during the first noise-whitening step for forming the two combined signals, the latter are determined retroactively using the maximum criterion.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This is the U.S. National Stage of PCT/FR2014/000041, filed Feb. 20, 2014, which in turn claims priority to French Patent Application No. 1351578, filed Feb. 22, 2013, the entire contents of all applications are incorporated herein by reference in their entireties.
TECHNICAL FIELD
This invention relates to a method for demodulating a signal.
It also relates to a multi-antenna receiving equipment and a corresponding computer programme.
The field of the invention is that of the encoding of digital data, intended to be transmitted in particular in the presence of transmission noise, and of the decoding of said digital data after transmission.
The invention more particularly but not in a limited manner relates to the field of optimising the transmission of digital data, for example via a broadband radio network.
PRIOR ART
A multi-antenna receiver comprises at least two receiving antennas in order to receive replicas of the same emitted signal, introducing as such additional diversity having for effect to improve the quality of the reception. The additional diversity provided is either a space diversity if the antennas are sufficiently spaced apart, or a polarisation diversity if the antennas are collocated and polarised differently, or both a portion of one and of the other of these two types of diversity.
The invention relates more particularly to multi-carrier transmission systems, for example of the OFDM type (Orthogonal Frequency Division Multiplexing). This multi-carrier transmission provides a frequency diversity that is separate from the space diversity and/or polarisation diversity provided by the use of several receiving antennas.
This type of transmission is increasingly used. It has in particular been adopted in the LTE (“Long Term Evolution”), TEDS (“TETRA Enhanced Data Service”), DAB (“Digital Audio Broadcasting”) and DVB-T (“Digital Video Broadcasting-Terrestrial”) systems.
In OFDM transmission systems, the data is in general organised into frames. In each frame, certain symbols, distributed in the temps-frequency plane are inserted among the useful information intended for the receiver. These symbols, referred to as “driver symbols” are known to the emitter and to the receiver. They are generally used for the purposes of synchronisation and estimation of the propagation channel.
The noise and interferences that the signal is subjected to during its propagation between the emitter and the receiver degrade the reception of the signal. An insufficient processing of the noise component and of the interference then generates a high demodulation and decoding error rate.
Conventionally, the noise component and the interference are processed by comparing a signal resulting from the combination of the signals received by the antennas of the receiver with an estimate of the signal emitted. This estimate of the signal emitted can be known a priori to the receiver by using for example the driver symbols. The estimate of the signal emitted can also be obtained through a prior treatment of the signals received on the various antennas, for example by means of a technique known as signal demodulation with a maximum combination of the signal ratio (MRC, for “Maximum Ratio Combining”).
Patent application FR 11 61062 describes a method for reducing interference. Although the method for reducing interference described in this document is of higher performance than the conventional methods for reducing interference in the case of the signals received with noise from strong interference, the method may not provide performance in certain cases wherein the noise component of the signals received contains little interference.
The purpose of this invention is to overcome at lease this problem raised by prior art.
DISCLOSURE OF THE INVENTION
Such an objective is achieved with a method for demodulating a signal in a receiver comprising at least two antennas each receiving a signal transmitted through an associated radio propagation channel, the signals received corresponding to the same emitted signal comprising time and frequency distributed symbol frames in which certain symbols, referred to as driver symbols, are known to the receiver, said method comprising the steps of:

    
    
        noise whitening in order to form two combined signals the noise components of which are separate;
        normalising the noise components of the combined signals in order to form two signals the noise components of which are separate and have equal average standards;
        performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, the noise components of which are separate and which have equal average standards
    
    


characterised in that the first step of noise whitening for forming the two combined signals comprises sub-steps of:

    
    
        weighting of each one of the signals received with respective first weighting vectors associated with a respective antenna of the receiver, the signal associated with the first antenna being weighted by a vector w1 and the signal associated with the second antenna being weighted by a vector w2,
        combining weighted received signals in order to form a first combined signal (c1),
        weighting a reference signal comprising said driver symbols with another weighting vector (Wd),
        comparing the first combined signal and the weighted reference signal in order to form an error (ε), and
        determining weighting vectors (w1, w2 and wd) using the maximum a posteriori criterion by maximising the probability of the occurring of said weighting vectors conditionally with the error obtained,
        weighting of each signal received with second weighting vectors, the signal received on the first antenna being weighted by the conjugate of the complex vector w2 and the signal received on the first antenna being weighted by the opposite of the conjugate of the complex vector w1,
        combining received signals weighted by the second weighting vectors in order to form a second combined signal (c2).
    
    


Driver symbol means any symbol known to the receiver, in particular both the driver symbols such as designated in the standards of known OFDM systems such as LTE and TEDS, and/or synchronisation symbols and/or symbols already demodulated and decided by the receiver in a prior processing step.
Thanks to the use of the maximum a posteriori approach (this can be referred to as MAP, for Maximum A Posteriori), this invention makes it possible to optimally reduce the level of the interferences in a multi-antenna receiver and as such improve the Signal to Interference plus Noise Ratio (SINR).
The method described hereinabove makes it possible to obtain two signals c1 and c2 of which the noise components are orthogonal between them. The signal c1 has a level of interference that is optimally reduced and the signal c2 of which the noise is orthogonal on each symbol therefore has a maximum level of interference. No information is lost when this combination is considered.
Thanks to the use of the signal-to-noise ratio maximum combination approach (Maximum Ratio Combining or MRC) in a combination of signals such has just been described, i.e. of two signals with one having a noise component of which the level of interference is optimally reduced, the other a noise component comprising a maximum level of interference, and to the normalising of the noise components of these two signals, the method proposed responds to the problem raised by prior art.
As such, the method performs in all of the cases, including when the noise component of the signals received contains little interference.
In a preferred embodiment, the respective coefficients (α1 and α2) of the channel seen by the useful signal associated with a respective antenna of the receiver and used in the step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, r′ and r′′, obtained at the end of the noise-whitening step, can be determined using the maximum a posteriori criterion by maximising the probability of the occurring of the channel (α1 and α2) conditionally with the error present on the one hand in the signal c1 and on the other hand in the signal c2, by taking into account the propagation channel seen in the signal c1 and in the signal c2.
Preferably, the step of determining the whitening step can comprise the sub-steps of:

    
    
        calculating a covariance matrix GB of the propagation channel;
        calculating a covariance matrix GD of the product, symbol by symbol, of two propagation channels; and
        determining weighting vectors using covariance matrices GB and GD.
    
    


Indeed, the weighting vectors are homogeneous with the propagation channel. The covariance matrices GB and GD make it possible to represent the time and frequency constraints of this channel.
In addition, the step of determining the whitening step can further comprise a sub-step of decomposing the covariance matrix GB into eigenvectors according to the relationship GB=CBΛBCBH wherein CB is a matrix of eigenvectors of the matrix GB and ΛB is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB and a sub-step of decomposing the covariance matrix GD into eigenvectors according to the relationship GD=CDΛDCDH wherein CD is a matrix of eigenvectors of the matrix GD and Λ is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD.
These decompositions make it possible to simplify the calculations of the weighting vectors.
Advantageously, the step of determining the whitening step can further comprise a sub-step of extracting a matrix C′B from the matrix CB and a matrix Λ′B from the matrix ΛB, the matrix Λ′B comprising a determined number nB′ of eigenvalues of ΛB and the matrix C′B′ comprising the eigenvectors associated with these nB′ eigenvalues, and a sub-step of extracting a matrix C′D from the matrix CD and a matrix Λ′D from the matrix ΛD, the matrix Λ′D comprising a determined number nD′ of eigenvalues of ΛD and the matrix C′D comprising the eigenvectors associated with these nD′ eigenvalues, the weighting vectors then being determined using matrices of eigenvectors C′B and C′D and of eigenvalues Λ′B and Λ′D.
Preferably, the nB′ eigenvalues retained are the largest eigenvalues of the matrix ΛB and the n′D eigenvalues retained are the largest eigenvalues of the matrix ΛD.
This makes it possible to reduce the complexity of the calculation of the weighting vectors while still minimising performance losses thanks to the fact that the eigenvectors retained in the matrices C′B and C′D represent a major portion of the energy, represented by the sum of the eigenvalues retained in the matrices Λ′B and Λ′D.
Preferably, n′B and n′D are less than or equal to the number of driver symbols in each frame.
In addition, the covariance matrix GB can be the Kronecker product of a covariance matrix GB,T in the time domain and of a covariance matrix GB,F in the frequency domain and the covariance matrix GD is the Kronecker product of a covariance matrix GD,T in the time domain and of a covariance matrix GD,F in the frequency domain and the step of determining (28) can comprise the sub-steps of:

    
    
        decomposing the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,TH and GB,F=CB,FΛB,FCB,FH, wherein:
        CB,T is a matrix of eigenvectors of the matrix GB,T;
        ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
        CB,F is a matrix of eigenvectors of the matrix GB,F and
        ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
        extracting of a matrix from the matrix C′B,T and a matrix ΛB,T from the matrix ΛB,T the matrix ΛB,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues;
        extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
        decomposing of the covariance matrices GD,T and GD,F into eigenvectors according to the relationships GD,T=CD,TΛD,TH and GD,F=CD,FΛD,FH, wherein:
        CD,T is a matrix of eigenvectors of the matrix GD,T;
        ΛD,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,T;
        CD,F is a matrix of eigenvectors of the matrix GD,F; and
        ΛD,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,F;
        extracting of a matrix C′D,T from the matrix CD,T and a matrix Λ′D,T from the matrix ΛD,T, the matrix Λ′D,T comprising a determined number nDt′ of the largest eigenvalues of ΛD,T and the matrix C′D,T comprising the eigenvectors associated with these nDt′ largest eigenvalues; and
        extracting of a matrix C′D,F from the matrix CD,F and a matrix ΛD,F from the matrix ΛD,F, the matrix Λ′D,F comprising a determined number nDf′ of the largest eigenvalues of ΛD,F and the matrix C′D,F comprising the eigenvectors associated with these nDf′ largest eigenvalues,
    
    


the weighting vectors then being determined using matrices of eigenvectors C′B,T, C′B,F, C′D,F and of eigenvalues Λ′B,T, ΛB,F, Λ′D,T, ΛD,F.
The Kronecker product of two matrices A and B, the matrix A having for components (aij), where i is an integer between 1 and m and j is an integer between 1 and n, is the matrix product noted as A{circle around (×)}B and defined by the following expression:




  
    A
    ⊗
    B
  
  =
  
    
      (
      
        
          
            
              
                a
                11
              
              ⁢
              B
            
          
          
            ...
          
          
            ...
          
          
            
              
                a
                
                  1
                  ⁢
                  n
                
              
              ⁢
              B
            
          
        
        
          
            ⋮
          
          
            ...
          
          
            ...
          
          
            ⋮
          
        
        
          
            
              
                a
                
                  m
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              ⁢
              B
            
          
          
            ...
          
          
            ...
          
          
            
              
                a
                mn
              
              ⁢
              B
            
          
        
      
      )
    
    .
  




This use of the separability of the covariance matrix between the time domain and the frequency domain makes it possible to further reduce the complexity of the calculation of the weighting vectors while still guaranteeing performance very close to that obtained by using the entire covariance matrix.
In addition, the step of determining the signal demodulating step can comprise the sub-steps of:

    
    
        calculating a covariance matrix GB of the propagation channel;
        determining weighting vectors using the covariance matrix GB.
    
    


Advantageously, the step of determining the signal demodulating step can further comprise a sub-step of decomposing the covariance matrix GB into eigenvectors according to the relationship GB=CBΛBCBH wherein CB is a matrix of eigenvectors of the matrix GB and ΛB is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB.
Preferably, the step of determining the signal demodulating step can further comprise a sub-step of extracting a matrix C′B from the matrix CB and a matrix Λ′B from the matrix ΛB, the matrix Λ′B comprising a determined number nB′ of eigenvalues of ΛB and the matrix C′B comprising the eigenvectors associated with these nB′ eigenvalues, the weighting vectors then being determined using the matrix of eigenvectors C′B and of eigenvalues Λ′B.
In addition, the covariance matrix GB can be the Kronecker product of a covariance matrix GB,T in the time domain and of a covariance matrix GB,F in the frequency domain and the step of determining the signal demodulating step comprises the sub-steps of:

    
    
        decomposing of the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,TH and GB,F=CB,FΛB,FCB,FH, wherein:
        CB,T is a matrix of eigenvectors of the matrix GB,T;
        ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
        CB,F is a matrix of eigenvectors of the matrix GB,F; and
        ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
        extracting of a matrix C′B,T from the matrix CB,T and a matrix Λ′B,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues;
        extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
    
    


the weighting vectors then being determined using matrices of eigenvectors C′B,T, C′B,F and of eigenvalues Λ′B,T, Λ′B,F.
Preferably, the signal emitted is a multi-carrier signal, in particular an OFDM signal.
The invention can also apply to single-carrier systems.
According to another aspect of the invention, a receiving equipment is proposed comprising at least two antennas each receiving a signal transmitted through an associated radio propagation channel, the signals received corresponding to the same emitted signal comprising time and frequency distributed symbol frames in which certain symbols, referred to as driver symbols, are known to the receiver, said method comprising means for:

    
    
        noise whitening in order to form two combined signals the noise components of which are separate;
        normalising the noise components of the combined signals in order to form two signals the noise components of which are separate and have equal average standards;
        step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, the noise components of which are separate and which have equal average standards;
    
    


characterised in that the means of noise whitening for forming the two combined signals comprises means for:

    
    
        weighting of each one of the signals received with respective first weighting vectors associated with a respective antenna of the receiver, the signal associated with the first antenna being weighted by a vector w1 and the signal associated with the second antenna being weighted by a vector w2,
        combining weighted received signals in order to form a first combined signal (c1).
        weighting a reference signal comprising said driver symbols with another weighting vector (Wd),
        comparing the first combined signal and the weighted reference signal in order to form an error (ε), and
        determining weighting vectors (w1, W2 and Wd) using the maximum a posteriori criterion by maximising the probability of the occurring of said weighting vectors conditionally with the error obtained,
        weighting of each signal received with second weighting vectors, the signal received on the first antenna being weighted by the conjugate of the complex vector w2 and the signal received on the first antenna being weighted by the opposite of the conjugate of the complex vector w1, combining received signals weighted by the second weighting vectors in order to form a second combined signal (c2).
    
    


Advantageously, the receiving equipment according to the invention can furthermore comprise means for determining the respective coefficients (α1 and α2) of the channel seen by the useful signal associated with a respective antenna of the receiver and used in the step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, r′ and r′′, obtained at the end of the noise-whitening step, can be determined using the maximum a posteriori criterion by maximising the probability of the occurring of the channel (α1 and α2) conditionally with the error present on the one hand in the signal c1 and on the other hand in the signal c2, by taking into account the propagation channel seen in the signal c1 and in the signal c2.
According to another aspect of the invention, a computer programme is proposed comprising instructions for the implementation of the method according to the invention when the programme is executed by at least one processor.



DESCRIPTION OF THE FIGURES AND EMBODIMENTS
Other advantages and particularities of the invention shall appear when reading the detailed description of implementations and embodiments that are in no way restricted, and the following annexed drawings:
FIG. 1 is a diagram of a receiving equipment according to the invention;
FIG. 2 is a diagram showing the operating principle of the method of demodulation;
FIG. 3 is a flow chart showing a first part of the operation of the method of demodulation according to the invention;
FIG. 4 shows a useful signal in a projection base; and
FIG. 5 is a flow chart showing a second part of the operation of the method of demodulation according to the invention.



These embodiments being in no way restricted, alternatives of the invention can in particular be considered that comprise only a selection of the characteristics described in what follows separated from the other characteristics described (even if this selection is isolated within a sentence comprising these other characteristics), if this selection of characteristics is sufficient to confer a technical advantage or for differentiating the invention in relation to prior art. This selection comprises at least one more preferably functional characteristic without structural details, or with only one portion of the structural details if this portion only is sufficient to confer a technical advantage or to differentiate the invention in relation to prior art.
FIG. 1 shows a receiving equipment 2 of an OFDM transmission system, such as for example a base station or a mobile terminal.
The receiver 2 comprises two antennas 4, 6 for the reception of two signals corresponding to the same OFDM signal transmitted from an emitter through two radio propagation channels associated with the antennas 4, 6. The channels associated with the antennas 4, 6 are supposed to respond to identical physical constraints.
The signal OFDM emitted by the emitter is organised into time and frequency distributed frames of symbols among which certain symbols, referred to as driver symbols, are known to the receiver 2 and are stored in a memory 8 of said receiver 2. Each frame comprises as such n symbols with nf sub-carriers and nt temps-symbols, n being equal to the product of nf and of nt.
The receiver 2 comprises an input module 10 comprising amplification, baseband formatting, sampling and guard interval suppression stages.
The receiver 2 further comprises means for the time-frequency conversion 12 of the signals received by the antennas 4, 6 and processed by the input module 10 in order to switch them from the time domain to the frequency domain. These means for time-frequency conversion 12 implement a Fast Fourier Transform (FFT).
The receiver 2 further comprises an interference reduction module 14 making it possible to reduce the level of interference in a useful signal resulting from the combination of the signals received by the two antennas 4, 6. These interferences can be caused, by way of examples, by the presence of scramblers emitting on the same radio channel as the emitter.
The interference reduction module 14 is able to provide data symbols wherein the contribution of the interference is minimised.
A demodulator 16 of the receiver 2 makes it possible to demodulate these data symbols into demodulated bits according to the same modulation technique as that used in the emitter. The demodulator 16 is furthermore arranged in order to determine a similarity of each demodulated bit. The similarity of a bit has a negative or positive soft value, in comparison with a hard value such as the binary value “1” or “0”, in order to indicate that the demodulator 16 delivers actual floating values that each have a sign which imposes a later decision, by a decoder 18 of the receiver 2, on the state of the corresponding bit, i.e. a decision on the “hard” value “0” or “1”,
The decoder 18 as such makes it possible to decode the demodulated bits supplied by the demodulator 16 according to the previously determined similarities. The decoder 18 implements a recoding that corresponds to the encoding used when the signal is emitted, for example a convolutive decoding that corrects the errors using the Viterbi algorithm.
FIG. 2 shows the principle of the method for reducing interference used by the interference reduction module 14.
By noting as s1 and s2 the signals received respectively on the antennas 4 and 6, the principle used by this invention consists in applying weightings w1 and w2 respectively to the signals s1 and s2, then in combining the two weighted signals, for example by adding them together, in order to form a signal c1 from which is subtracted an estimate of the emitted signal d weighted by a weighting wd. The resulting difference ε shows a residual error.
The method of the invention uses advantageously the maximum a posteriori approach in order to calculate the most probable weightings knowing this error.
The steps of this method are described in reference to the flow chart of FIG. 3.
During a first step 20, the interference reduction module 14 weights the signals s1 and s2 with first respective weighting vectors (w1; w2) associated with a respective antenna (4, 6) of the receiver (2). The first weighting vectors (w1; w2) define a minimum direction, noted as DIRmin, for the interference as is shown in FIG. 4. The content of the first weighting vectors (w1; w2) remains to be determined.
In the step 22, the interference reduction module 14 combines, here by adding them together, the weighted signals s1 and s2 in order to form a first combined signal (c1).
As such, we have: c1=w1s1+w2s2.
In the step 24, the interference reduction module 14 weights a reference signal d, comprising the driver symbols, with another weighting vector wd of which the content is to be determined.
The vectors w1, w2 and wd are column vectors each containing as many lines as the signal received in a frame contains symbols, i.e. n lines.
In the step 26, the interference reduction module 14 determines an error ε that corresponds to the difference between the first combined signal c1 and the weighted reference signal.
As such by noting as w the vector consisting of a concatenation of vectors w1, w2 and wd and R the matrix comprising the signals received and the reference signal defines according to the expression:




  R
  =
  
    [
    
      
        
          
            
              [
              
                
                  
                    
                      s
                      
                        1
                        ,
                        0
                        ,
                        0
                      
                    
                  
                  
                    0
                  
                  
                    ...
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      s
                      
                        1
                        ,
                        ...
                      
                    
                  
                  
                    0
                  
                  
                    ...
                  
                
                
                  
                    ...
                  
                  
                    0
                  
                  
                    
                      s
                      
                        1
                        ,
                        ...
                      
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    ...
                  
                  
                    0
                  
                  
                    
                      s
                      
                        1
                        ,
                        
                          n
                          f
                        
                        ,
                        
                          n
                          i
                        
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      s
                      
                        2
                        ,
                        0
                        ,
                        0
                      
                    
                  
                  
                    0
                  
                  
                    ...
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      s
                      
                        2
                        ,
                        ...
                      
                    
                  
                  
                    0
                  
                  
                    ...
                  
                
                
                  
                    ...
                  
                  
                    0
                  
                  
                    
                      s
                      
                        2
                        ,
                        ...
                      
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    ...
                  
                  
                    0
                  
                  
                    
                      s
                      
                        2
                        ,
                        
                          n
                          f
                        
                        ,
                        
                          n
                          i
                        
                      
                    
                  
                
              
              ]
            
          
        
      
      
        
          
            [
            
              
                
                  
                    -
                    
                      d
                      
                        0
                        ,
                        0
                      
                    
                  
                
                
                  0
                
                
                  ...
                
                
                  0
                
              
              
                
                  0
                
                
                  
                    -
                    d
                  
                
                
                  0
                
                
                  ...
                
              
              
                
                  ...
                
                
                  0
                
                
                  
                    -
                    d
                  
                
                
                  0
                
              
              
                
                  0
                
                
                  ...
                
                
                  0
                
                
                  
                    -
                    
                      d
                      
                        
                          n
                          f
                        
                        ,
                        
                          n
                          i
                        
                      
                    
                  
                
              
            
            ]
          
        
      
    
    ]
  




the error ε is a vector expressed by ε=R·w.
In the step 28, the interference reduction module 14 calculates the vector w1, w2 and wd by using the MAP approach.
This approach consists in maximising the probability of the occurring of the weighting vector w. This probability if conditional to the observation of the error ε.
By using the Bayes formula, this probability is equal to the probability that the error ε is observed conditionally to the probability that the weighting is equal to the vector w, this probability being multiplied by the probability that the weighting vector w is carried out. This is expressed by the following relationship:

P(W)=f(w/ε)∝f(ε/w)·f(w),

wherein P and f indicate probabilities and the symbol ∝ indicates the “equivalent” relationship.
However, the weighting vectors w1, w2 and wd are linked to the propagation channel. This channel is conditioned by constraints relating to its maximum time spread, due to reflections on distant obstacles, and to its maximum frequency spread, due to the speed of the receiving equipment and to the carrier frequency, i.e. the Doppler spread.
On the one hand, in the frequency domain, the frequency spread of the channel, called the Doppler spread, due to the reflections on near obstacles is limited. This spread is between −FD and +FD, wherein FD is the maximum Doppler frequency given by the relationship




  
    
      F
      D
    
    =
    
      
        v
        c
      
      ·
      
        F
        p
      
    
  
  ,




wherein v is the speed of the receiver 2, c is the speed of light, and Fp is the carrier frequency. The components of the frequency spectrum of the propagation channel according to the frequency axis are therefore between these limits −FD and +FD.

On the other hand, the time spread of the channel, due to the reflections on the distant obstacles, is limited. This time spread depends on the frequency band used and of the environment. By way of examples, at a carrier frequency of 400 MHz, in an urban environment the time spread is about 5 μs while in a mountainous environment, this spread is about 15 μs. The components of the time response of the channel are therefore between fixed limits for given environmental conditions.
The limits of the frequency spectrum and of the time response of the channel are known to the receiver 2 and are stored in the memory 8.
These characteristics of the propagation channel are represented by a matrix referred to as the global covariance matrix G.
Due to the constraints expressed in the global covariance matrix G, the probability of the occurrence of the weighting vector w is given by the expression:




  
    f
    ⁡
    
      (
      w
      )
    
  
  =
  
    β
    ×
    
      e
      
        
          
            -
            1
          
          /
          2
        
        ⁢
        
          w
          H
        
        ⁢
        
          G
          
            -
            1
          
        
        ⁢
        w
      
    
  




wherein β is a constant and the notation XH indicates that it is a conjugate and transposed matrix X.
Furthermore, under the hypothesis that the noise component of the channel id of the Gaussian type, the conditional probability of observing the error ε is given by the expression:




  
    f
    ⁡
    
      (
      
        ɛ
        /
        w
      
      )
    
  
  =
  
    α
    ×
    
      e
      
        
          -
          
            
              
              
                R
                ·
                w
              
              
            
            2
          
        
        /
        
          (
          
            2
            ⁢
            
              σ
              2
            
          
          )
        
      
    
  




wherein α is a constant and σ2 represents the variance of the noise component in the signal corresponding to the signals received on the various weighted and combined antennas.
As such the probability of the occurring, expressing the MAP criterion, of the weighting vector w is given by the following expression [1]:




  
    
      
        
          P
          ⁡
          
            (
            w
            )
          
        
        =
        
          α
          ×
          
            e
            
              -
              
                
                  
                    
                    
                      R
                      ·
                      w
                    
                    
                  
                  2
                
                
                  (
                  
                    2
                    ·
                    
                      σ
                      2
                    
                  
                  )
                
              
            
          
          ×
          β
          ×
          
            
              e
              
                
                  
                    
                      -
                      1
                    
                    /
                    2
                  
                  ·
                  
                    w
                    H
                  
                
                ⁢
                
                  G
                  
                    -
                    1
                  
                
                ⁢
                w
              
            
            .
          
        
      
    
    
      
        [
        1
        ]
      
    
  




Furthermore, the variance, or power, σ2 is equal to:

σ2=σ12∥w1∥2+σ22∥w2∥2,

wherein σ12 is the variance of the noise component on the first antenna 4 and σ22 is the variance of the noise component on the second antenna 6.
Assuming that the power of the noise component is the same on each one of the two antennas, i.e. that σ12=σ22 the variance of the total noise component is then equal to

σ2=σ12(∥w1∥2+∥w2∥2)(.

By taking the opposite of the logarithm of the expression [1], the following expression [2] is obtained:




  
    
      
        
          
            L
            ⁡
            
              (
              
                P
                ⁡
                
                  (
                  w
                  )
                
              
              )
            
          
          =
          
            
              
                
                  
                  
                    R
                    ·
                    w
                  
                  
                
                2
              
              
                
                  σ
                  1
                  2
                
                ⁡
                
                  (
                  
                    
                      
                        
                        
                          w
                          1
                        
                        
                      
                      2
                    
                    +
                    
                      
                        
                        
                          w
                          2
                        
                        
                      
                      2
                    
                  
                  )
                
              
            
            +
            
              
                w
                H
              
              ⁢
              
                G
                
                  -
                  1
                
              
              ⁢
              w
            
            +
            cst
          
        
        ,
      
    
    
      
        [
        2
        ]
      
    
  




wherein cst designates a constant value.
The interference reduction module 14 attempts to minimise this logarithm L(P(w)).
Moreover, the covariance matrix G is a diagonal matrix by blocks constituted by the concatenation of covariance matrices corresponding to each one of the weighting vectors w1, w2 and wd.
The covariance matrix of w1 is the same as that of w2 given that the two weighting vectors w1 and w2 are both homogeneous with a propagation channel, on the case of two receiving antennas. This covariance matrix is noted as GB and shows the time and frequency constraints relative to the propagation channel.
The weighting vector wd is homogeneous with the product, symbol by symbol, of two propagation channels, on the case of two receiving antennas. The corresponding covariance matrix GD shows the time and frequency constraints relative to such a product.
The covariance matrix G can therefore be written as:




  G
  =
  
    
      (
      
        
          
            
              G
              B
            
          
          
            0
          
          
            0
          
        
        
          
            0
          
          
            
              G
              B
            
          
          
            0
          
        
        
          
            0
          
          
            0
          
          
            
              G
              D
            
          
        
      
      )
    
    .
  




Furthermore, the following constraint is imposed:

k2=∥w1∥2+∥w2∥2=cst2.

This constraint does not limit the solutions to the minimisation problem of the expression [2]. Indeed, regardless of the value of k2, the solution w does not change given that changing this constant would only multiply the solution by a scalar constant.
In what follows, k2 is set to 2. This makes it possible to consider that the weightings w1 and w2 are each of unitary power. Then the matrix GB is the normalised covariance matrix, i.e. obtained with a unitary average power of the channel, and the matrix GD is the covariance matrix corresponding to the product, symbol by symbol, of two channels of unitary power. Consequently, σ1 shows the inverse′ of the signal-to-noise ratio observed on any of the antennas.
It is possible to set a target value for the signal-to-noise ratio, with this value not changing regardless of the effective value of the puissance of the noise component and of the useful signal power. Other solutions are possible as, for example, estimating the signal-to-noise ratio as the information received is received and demodulated, for example on the driver symbols which are known to the emitter and to the receiver. In this case it is possible to dynamically adapt the signal-to-noise ratio in the process of estimating the channel.
The expression [2] can then be written as:




  
    L
    ⁡
    
      (
      
        P
        ⁡
        
          (
          w
          )
        
      
      )
    
  
  =
  
    
      
        
          
          
            R
            ·
            w
          
          
        
        2
      
      
        2
        ·
        
          σ
          1
          2
        
      
    
    +
    
      
        w
        H
      
      ⁢
      
        G
        
          -
          1
        
      
      ⁢
      w
    
    +
    
      cst
      .
    
  




The minimisation problem of the step 28 is as such an optimisation problem with constraint, which is resolved with the Lagrange multipliers according to the following relationship:




  
    
      ∇
      w
    
    ⁢
    
      (
      
        L
        ⁡
        
          (
          
            P
            ⁡
            
              (
              w
              )
            
          
          )
        
      
      )
    
  
  =
  
    
      
        (
        
          
            
              
                R
                H
              
              ⁢
              R
            
            
              2
              ·
              
                σ
                1
                2
              
            
          
          +
          
            G
            
              -
              1
            
          
        
        )
      
      ⁢
      w
    
    -
    
      μ
      ⁢
      
        
      
      ⁢
      
        N
        H
      
      ⁢
      Nw
    
  




wherein ∇w designates the gradient in relation to weighting vector w and the matrix N is the identity matrix for the two signals received at the antennas 4, 6 and the zero matrix for the reference signal.
In other terms, the matrix N can be written as:




  N
  =
  
    
      (
      
        
          
            1
          
          
            0
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            0
          
        
        
          
            0
          
          
            ⋱
          
          
            ⋱
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            ⋮
          
        
        
          
            ⋮
          
          
            ⋱
          
          
            1
          
          
            ⋱
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            ⋮
          
        
        
          
            ⋮
          
          
            
              
            
          
          
            ⋱
          
          
            1
          
          
            ⋱
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            ⋮
          
        
        
          
            ⋮
          
          
            
              
            
          
          
            
              
            
          
          
            ⋱
          
          
            ⋱
          
          
            ⋱
          
          
            
              
            
          
          
            
              
            
          
          
            ⋮
          
        
        
          
            ⋮
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            ⋱
          
          
            1
          
          
            ⋱
          
          
            
              
            
          
          
            ⋮
          
        
        
          
            ⋮
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            ⋱
          
          
            0
          
          
            ⋱
          
          
            ⋮
          
        
        
          
            ⋮
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            
              
            
          
          
            ⋱
          
          
            ⋱
          
          
            0
          
        
        
          
            0
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
          
            0
          
          
            0
          
        
      
      )
    
    .
  




The minimisation of the gradient then gives the following expression [3]:




  
    
      
        
          
            
              ∇
              w
            
            ⁢
            
              (
              
                L
                ⁡
                
                  (
                  
                    P
                    ⁡
                    
                      (
                      w
                      )
                    
                  
                  )
                
              
              )
            
          
          =
          
            
              0
              ⇒
              
                
                  (
                  
                    
                      
                        
                          R
                          H
                        
                        ⁢
                        R
                      
                      
                        2
                        ·
                        
                          σ
                          1
                          2
                        
                      
                    
                    +
                    H
                  
                  )
                
                ⁢
                w
              
            
            =
            
              μ
              ⁢
              
                
              
              ⁢
              
                N
                H
              
              ⁢
              Nw
            
          
        
        ,
      
    
    
      
        [
        3
        ]
      
    
  




wherein the matrix H is equal to the inverse of the covariance matrix G.
The weighting vector w solution of the minimisation problem is as such a generalised eigenvector both of the term




  (
  
    
      
        
          R
          H
        
        ⁢
        R
      
      
        2
        ·
        
          σ
          1
          2
        
      
    
    +
    H
  
  )




and of the term NHN and μ is the associated generalised eigenvalue.

Furthermore, the solution must satisfy the relationship L(P(w))=μ+cst.
As such, the weighting vector w minimising L(P(w)) is the generalised eigenvector corresponding to the smallest generalised eigenvalue μ.
According to a preferred embodiment, in order to reduce the complexity of the calculations, the covariance matrices GB and GD are decomposed in the step 28 into eigenvectors and eigenvalues in the following way:




  {
  
    
      
        
          
            G
            =
            
              (
              
                
                  
                    
                      G
                      B
                    
                  
                  
                    0
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      G
                      B
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    0
                  
                  
                    
                      G
                      D
                    
                  
                
              
              )
            
          
        
      
      
        
          
            {
            
              
                
                  
                    
                      
                        
                          G
                          B
                        
                        =
                        
                          C
                          B
                        
                      
                      ,
                      
                        Λ
                        B
                      
                      ,
                      
                        C
                        B
                        H
                      
                    
                  
                
                
                  
                    
                      
                        
                          G
                          D
                        
                        =
                        
                          C
                          D
                        
                      
                      ,
                      
                        Λ
                        D
                      
                      ,
                      
                        C
                        D
                        H
                      
                    
                  
                
              
              ⇒
              
                {
                
                  
                    
                      
                        
                          
                            G
                            B
                            
                              -
                              1
                            
                          
                          =
                          
                            C
                            B
                          
                        
                        ,
                        
                          Λ
                          B
                          
                            -
                            1
                          
                        
                        ,
                        
                          C
                          B
                          H
                        
                      
                    
                  
                  
                    
                      
                        
                          
                            G
                            D
                            
                              -
                              1
                            
                          
                          =
                          
                            C
                            D
                          
                        
                        ,
                        
                          Λ
                          D
                          
                            -
                            1
                          
                        
                        ,
                        
                          C
                          D
                          H
                        
                      
                    
                  
                
              
            
          
        
      
    
    ,
  




where the matrices CB and CD are matrices of eigenvectors of the matrices GB and GD respectively and the matrices ΛB and ΛD are the diagonal matrices of corresponding eigenvalues respectively.
The covariance matrix G can as such be written G=CΛCH wherein C is a matrix of eigenvectors of the matrix G and Λ is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix G.
As such, the inverse matrix of the covariance matrix G is equal to:




  {
  
    
      
        
          
            
              G
              
                -
                1
              
            
            =
            
              (
              
                
                  
                    
                      
                        C
                        B
                      
                      ,
                      
                        Λ
                        B
                        
                          -
                          1
                        
                      
                      ,
                      
                        C
                        B
                        H
                      
                    
                  
                  
                    0
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      
                        C
                        B
                      
                      ,
                      
                        Λ
                        B
                        
                          -
                          1
                        
                      
                      ,
                      
                        C
                        B
                        H
                      
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    0
                  
                  
                    
                      
                        C
                        D
                      
                      ,
                      
                        Λ
                        D
                        
                          -
                          1
                        
                      
                      ,
                      
                        C
                        D
                        H
                      
                    
                  
                
              
              )
            
          
        
      
      
        
          soit
        
      
      
        
          
            
              
                G
                
                  -
                  1
                
              
              =
              C
            
            ,
            
              Λ
              
                -
                1
              
            
            ,
            
              C
              H
            
          
        
      
    
    .
  




According to an embodiment, the complexity of the calculation is further reduced by retained only certain eingenvalues of the covariance matrix G as well as the corresponding eigenvalues. In this case, a matrix C′ is extracted from the matrix C and a matrix Λ′ is extracted from the matrix Λ, the matrix Λ′ comprising a determined number n′ of eigenvalues of Λ and the matrix C′ comprising the eigenvectors associated with these n′ eigenvalues. The weighting vector w is then determined using matrices of eigenvectors C′ and of eigenvalues Λ′.
In a preferred embodiment, the covariance matrix G is a diagonal matrix par bloc, constituted of the matrices GB, GB and GD, and for which each one of the matrices GB and GD is the Kronecker product of a covariance matrix GB,T, respectively GD,T, in the time domain and of a covariance matrix GB,F, respectively GD,F, in the frequency domain and the step of determining comprises the sub-steps of:

    
    
        decomposing of the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,TH and GB,F=CB,FΛB,FCB,FH, wherein:
        CB,T is a matrix of eigenvectors of the matrix GB,T;
        ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
        CB,F is a matrix of eigenvectors of the matrix GB,F; and
        ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
        extracting of a matrix C′B,T from the matrix CB,T and a matrix Λ′VB,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues; and
        extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
        decomposing of covariance matrices GD,T and GD,F into eigenvectors according to the relationships GD,T=CD,TΛD,TCD,TH and GD,F=CD,FΛD,FCD,FH, wherein:
        CD,T is a matrix of eigenvectors of the matrix GD,T;
        ΛD,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,T;
        CD,F is a matrix of eigenvectors of the matrix GD,F; and
        ΛD,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,F;
        extracting of a matrix C′D,T from the matrix CD,T and a matrix Λ′D,T from the matrix ΛD,T, the matrix Λ′D,T comprising a determined number nDt′ of the largest eigenvalues of ΛD,T and the matrix C′D,T comprising the eigenvectors associated with these nDt′ largest eigenvalues; and
        extracting of a matrix C′D,F from the matrix CD,F and a matrix Λ′D,F from the matrix ΛD,F, the matrix Λ′D,F comprising a determined number nDf′ of the largest eigenvalues of ΛD,F and the matrix C′D,F comprising the eigenvectors associated with these nDf′ largest eigenvalues.
    
    


The weighting vector w is then determined using matrices of eigenvectors C′B,T, C′B,F, C′D,T, C′D,F and of eigenvalues Λ′B,T, Λ′B,FΛ′D,T, Λ′D,F.
Returning to the expression [3] and by using the decomposing into eigenvectors and into eingenvalues of the covariance matrix, the following is obtained:




  
    
      (
      
        
          
            
              R
              H
            
            ⁢
            R
          
          
            σ
            2
          
        
        +
        H
      
      )
    
    ⁢
    w
  
  =
  
    μ
    ⁢
    
      
    
    ⁢
    
      N
      H
    
    ⁢
    Nw
    ⁢
    
      {
      
        
          
            
              
                ⇔
                
                  
                    (
                    
                      
                        
                          
                            
                              R
                              H
                            
                            ⁢
                            R
                          
                          
                            ⁣
                            
                              2
                              ,
                              
                                σ
                                2
                              
                            
                          
                        
                        +
                        C
                      
                      ,
                      
                        Λ
                        
                          -
                          1
                        
                      
                      ,
                      
                        C
                        H
                      
                    
                    )
                  
                  ⁢
                  w
                
              
              =
              
                μ
                ⁢
                
                  
                
                ⁢
                
                  N
                  H
                
                ⁢
                Nw
              
            
          
        
        
          
            
              
                ⇔
                
                  
                    (
                    
                      
                        
                          
                            C
                            H
                          
                          ⁢
                          
                            R
                            H
                          
                          ⁢
                          RC
                        
                        
                          2
                          ,
                          
                            σ
                            2
                          
                        
                      
                      +
                      
                        Λ
                        
                          -
                          1
                        
                      
                    
                    )
                  
                  ⁢
                  y
                
              
              =
              
                μ
                ⁢
                
                  
                
                ⁢
                
                  C
                  H
                
                ⁢
                
                  N
                  H
                
                ⁢
                NCy
              
            
          
        
      
    
  




wherein y=CH w wherein y is the concatenation of vectors y1, y2, and yd:




  y
  =
  
    
      [
      
        
          
            
              y
              1
            
          
        
        
          
            
              y
              2
            
          
        
        
          
            
              y
              d
            
          
        
      
      ]
    
    =
    
      
        
          C
          H
        
        ⁢
        w
      
      =
      
        
          [
          
            
              
                
                  C
                  B
                  H
                
              
              
                0
              
              
                0
              
            
            
              
                0
              
              
                
                  C
                  B
                  H
                
              
              
                0
              
            
            
              
                0
              
              
                0
              
              
                
                  C
                  D
                  H
                
              
            
          
          ]
        
        ·
        
          [
          
            
              
                
                  w
                  1
                
              
            
            
              
                
                  w
                  2
                
              
            
            
              
                
                  w
                  d
                
              
            
          
          ]
        
      
    
  




The problem can be further simplified by separating the equations in y1 and y2 from that in yd. This can be done by multiplying the equation on the left by M=I−N, respectively the equation on the left by N.



  
    
      
        (
        
          
            
              1
              
                2
                ⁢
                
                  
                
                ⁢
                
                  σ
                  2
                
              
            
            ⁢
            
              MC
              H
            
            ⁢
            
              R
              H
            
            ⁢
            RC
          
          +
          
            M
            ⁢
            
              
            
            ⁢
            
              Λ
              
                -
                1
              
            
          
        
        )
      
      ⁢
      
        y
        _
      
    
    =
    
      
        μ
        ⁢
        
          
        
        ⁢
        
          MC
          H
        
        ⁢
        NC
        ⁢
        
          y
          _
        
      
      =
      0
    
  




  
    
      
        (
        
          
            
              1
              
                2
                ⁢
                
                  
                
                ⁢
                
                  σ
                  2
                
              
            
            ⁢
            
              NC
              H
            
            ⁢
            
              R
              H
            
            ⁢
            RC
          
          +
          
            N
            ⁢
            
              
            
            ⁢
            
              Λ
              
                -
                1
              
            
          
        
        )
      
      ⁢
      
        y
        _
      
    
    =
    
      
        μ
        ⁢
        
          
        
        ⁢
        
          NC
          H
        
        ⁢
        NC
        ⁢
        
          y
          _
        
      
      =
      
        μ
        ⁢
        
          
        
        ⁢
        
          C
          H
        
        ⁢
        NC
        ⁢
        
          y
          _
        
      
    
  



The solution for the first equation in yd gives:



  
    
      
    
    ⁢
    
      
        
          
            (
            
              
                
                  C
                  D
                  H
                
                ⁢
                
                  d
                  H
                
                ⁢
                d
                ⁢
                
                  
                
                ⁢
                
                  C
                  D
                
              
              +
              
                2
                ⁢
                
                  
                
                ⁢
                
                  σ
                  2
                
                ⁢
                
                  Λ
                  D
                  
                    -
                    1
                  
                
              
            
            )
          
          ⁢
          
            y
            d
          
        
        =
        
          
            [
            
              
                
                  
                    
                      C
                      D
                      H
                    
                    ⁢
                    
                      d
                      H
                    
                    ⁢
                    
                      s
                      1
                    
                    ⁢
                    
                      C
                      B
                    
                  
                
                
                  
                    
                      C
                      D
                      H
                    
                    ⁢
                    
                      d
                      H
                    
                    ⁢
                    
                      s
                      2
                    
                    ⁢
                    
                      C
                      B
                    
                  
                
              
            
            ]
          
          ⁢
          
            
              y
              12
            
            _
          
        
      
      ,
      
        

      
      ⁢
      
        
      
      ⁢
      
        
          
            y
            12
          
          _
        
        ≡
        
          [
          
            
              
                
                  y
                  1
                
              
            
            
              
                
                  y
                  2
                
              
            
          
          ]
        
      
    
  




  
    Where
    ⁢
    
      
    
    ⁢
    
      (
      
        
          
            
              1
              
                2
                ⁢
                
                  
                
                ⁢
                
                  σ
                  2
                
              
            
            ⁡
            
              [
              
                
                  
                    
                      
                        C
                        B
                        H
                      
                      ⁢
                      
                        s
                        1
                        H
                      
                      ⁢
                      
                        s
                        1
                      
                      ⁢
                      
                        C
                        B
                      
                    
                  
                  
                    
                      
                        C
                        B
                        H
                      
                      ⁢
                      
                        s
                        1
                        H
                      
                      ⁢
                      
                        s
                        2
                      
                      ⁢
                      
                        C
                        B
                      
                    
                  
                
                
                  
                    
                      
                        C
                        B
                        H
                      
                      ⁢
                      
                        s
                        2
                        H
                      
                      ⁢
                      
                        s
                        1
                      
                      ⁢
                      
                        C
                        B
                      
                    
                  
                  
                    
                      
                        C
                        B
                        H
                      
                      ⁢
                      
                        s
                        2
                        H
                      
                      ⁢
                      
                        s
                        2
                      
                      ⁢
                      
                        C
                        B
                      
                    
                  
                
              
              ]
            
          
          -
          
            
              
                1
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  
                    σ
                    2
                  
                
              
              ⁡
              
                [
                
                  
                    
                      
                        
                          C
                          B
                          H
                        
                        ⁢
                        
                          s
                          1
                          H
                        
                        ⁢
                        d
                        ⁢
                        
                          
                        
                        ⁢
                        
                          C
                          D
                        
                      
                    
                  
                  
                    
                      
                        
                          C
                          B
                          H
                        
                        ⁢
                        
                          s
                          2
                          H
                        
                        ⁢
                        d
                        ⁢
                        
                          
                        
                        ⁢
                        
                          C
                          D
                        
                      
                    
                  
                
                ]
              
            
            ⁢
            
              
                (
                
                  
                    
                      C
                      D
                      H
                    
                    ⁢
                    
                      d
                      H
                    
                    ⁢
                    d
                    ⁢
                    
                      
                    
                    ⁢
                    
                      C
                      D
                    
                  
                  +
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      σ
                      2
                    
                    ⁢
                    
                      Λ
                      D
                      
                        -
                        1
                      
                    
                  
                
                )
              
              
                -
                1
              
            
            ⁢
            
               
              
                
                  [
                  
                    
                      
                        
                          
                            C
                            D
                            H
                          
                          ⁢
                          
                            d
                            H
                          
                          ⁢
                          
                            s
                            1
                          
                          ⁢
                          
                            C
                            B
                          
                        
                      
                      
                        
                          
                            C
                            D
                            H
                          
                          ⁢
                          
                            d
                            H
                          
                          ⁢
                          
                            s
                            2
                          
                          ⁢
                          
                            C
                            B
                          
                        
                      
                    
                  
                  ]
                
                +
                
                  
                    [
                    
                      
                        
                          
                            Λ
                            B
                          
                        
                        
                          0
                        
                      
                      
                        
                          0
                        
                        
                          
                            Λ
                            B
                          
                        
                      
                    
                    ]
                  
                  
                    -
                    1
                  
                
              
              )
            
            ⁢
            
              
                y
                12
              
              _
            
          
        
        =
        
          μ
          ⁢
          
            
          
          ⁢
          
            
              y
              12
            
            _
          
        
      
    
  



After having determined the vector y using the preceding expression, the interference reduction module 14 calculates the weighting vector w according to the relationship:

w=C·y. 

In the step 30, the interference reduction module 14 weights the signals received s1 and s2 respectively by weighting vectors w*2 and −w*1. As such the vector (w1, w2), referred to as the first weighting vector, formed by the first weighting vectors (w1, w2) is orthogonal (i.e. the hermitian product is zero) on each symbol to the vector (w*2, −w*1), referred to as the second weighting vector, formed by the second weighting vectors (w*2, −w*1). The second weighting vectors (w*2, −w*1) define a maximum direction, noted as DIRmax, for the interference as is shown in FIG. 4.
The notation X* indicates that it is a conjugate matrix X.
In the step 32, the interference reduction module 14 combines, here by adding them together, the weighted signals c1 and c2 by the second weighting vectors (w*2, −w*1) in order to form a combined signal c2.
As such, we have: c2=w*2s1−w*1s2.
Two signals c1 and c2 are thus obtained of which the noise components are not correlated as they are orthogonal.
Then, a step of normalising (not shown) of the noise components of the combined signals (c1; c2) in order to form two signals (r′; r′′) the noise components of which are separate and have equal average standards is carried out.
The following steps of the method are described in reference to the flow chart of FIG. 5.
On note (α1, α2) the respective coefficients of the channel seen by the useful signal associated with a respective antenna (4, 6) of the receiver (2).
The vectors α1, α2 are calculated by using the MAP approach.
During a first step 40, the interference reduction module 14 weights the signals received s1 and s2 respectively by weighting vectors w1 and w2 determined previously.
In the step 42, the interference reduction module 14 combines, here by adding them together, the weighted signals s1 and s2 in order to form a combined signal c1.
As such, the following signal is obtained: c1=w1s1+w2s2.
In the step 44, the interference reduction module 14 weights a reference signal d, comprising the driver symbols, with another weighting vector wd.
It is observed that wd can be obtained by combining a weighting of the first weighting vector (w1, w2) determined during the step of determining (28) of the step of whitening E1 by the channel vector (α1, α2). As such, ideally: wd=w1α1+w2α2.
The vectors α1, α2 are column vectors each containing as many lines as the signal received in a frame contains symbols, i.e. n lines.
In the step 46, the interference reduction module 14 determines an error ε that corresponds to the difference between the combined signal obtained c1=w1s1+w2s2 and the weighted reference signal d of wd=w1α1+w2α2.
As such by noting as α the vector consisting of a concatenation of vectors w1, w2 and wd=w1α1+w2α2 and the matrix comprising the signals received as defined hereinabove, the error ε still exists which is the vector expressed by ε=Rα.
In the step 48, the interference reduction module 14 calculates the vector α1, α2 by using the MAP approach.
This approach consists in maximising the probability of the occurring of the weighting vector α. This probability is conditional to the observation of the error ε.
By using the Bayes formula, this probability is equal to the probability that the error ε is observed conditionally to the probability that the weighting is equal to the vector α, this probability being multiplied by the probability that the weighting vector α is carried out. This is expressed by the following relationship:

P(α)=f(α/ε)∝f(ε/α)·f(α),

wherein P and f indicate probabilities and the symbol ∝ indicates the “equivalent to” relationship.
As previously, the limits of the frequency spectrum and of the time response of the channel are known to the receiver 2 and are stored in the memory 8.
These characteristics of the propagation channel are represented by a matrix referred to as the global covariance matrix G2.




  {
  
    
      
        
          
            
              G
              2
            
            =
            
              (
              
                
                  
                    
                      G
                      B
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      G
                      B
                    
                  
                
              
              )
            
          
        
      
      
        
          
            
              
                G
                B
              
              =
              
                C
                B
              
            
            ,
            
              Λ
              B
            
            ,
            
              
                
                  C
                  B
                  H
                
                ⇒
                
                  G
                  B
                  
                    -
                    1
                  
                
              
              =
              
                C
                B
              
            
            ,
            
              Λ
              B
              
                -
                1
              
            
            ,
            
              C
              B
              H
            
          
        
      
    
     
  




Due to the constraints expressed in the global covariance matrix G2, the probability of the occurrence of the weighting vector α is given by the expression:




  
    f
    ⁡
    
      (
      α
      )
    
  
  =
  
    
      β
      ′
    
    ×
    
      e
      
        
          
            
              -
              1
            
            /
            2
          
          ·
          
            α
            H
          
        
        ⁢
        
          G
          2
          
            -
            1
          
        
        ⁢
        α
      
    
  




wherein β′ is a constant.
Furthermore, under the hypothesis that the noise component of the channel id of the Gaussian type, the conditional probability of observing the error ε is given by the expression:




  
    f
    ⁡
    
      (
      
        α
        /
        w
      
      )
    
  
  =
  
    
      α
      ′
    
    ×
    
      
        e
        
          -
          
            
              
              
                R
                ·
                α
              
              
            
            2
          
        
      
      /
      
        (
        
          2
          ,
          
            σ
            
              -
              2
            
          
        
        )
      
    
  




wherein α′ is a constant and α′2 represents the variance of the noise component in the signal corresponding to the signals received on the various weighted and combined antennas.
As such the probability of the occurring, expressing the MAP criterion, of the weighting vector a is given by the following expression [1]:




  
    
      
        
          P
          ⁡
          
            (
            α
            )
          
        
        =
        
          
            α
            ′
          
          ⁢
          
            e
            
              -
              
                
                  
                    
                    
                      R
                      ·
                      α
                    
                    
                  
                  2
                
                
                  (
                  
                    2
                    ,
                    
                      σ
                      2
                    
                  
                  )
                
              
            
          
          ×
          
            β
            ′
          
          ×
          
            e
            
              
                
                  -
                  1
                
                /
                2
              
              ,
              
                
                  α
                  H
                
                ⁢
                
                  G
                  2
                  
                    -
                    1
                  
                
                ⁢
                α
              
            
          
        
      
    
    
      
        [
        4
        ]
      
    
  




By taking the opposite of the logarithm of the expression [4] the following expression [5] is obtained:




  
    
      
        
          L
          ⁡
          
            (
            
              P
              ⁡
              
                (
                α
                )
              
            
            )
          
        
        =
        
          
            
              
                
                
                  R
                  ,
                  α
                
                
              
              2
            
            
              σ
              2
            
          
          +
          
            
              α
              H
            
            ⁢
            
              G
              2
              
                -
                1
              
            
            ⁢
            α
          
          +
          
            cst
            ′
          
        
      
    
    
      
        [
        5
        ]
      
    
  




wherein cst′ designates a constant value.
The interference reduction module 14 attempts to minimise this logarithm L(P(α)), which reverts to minimising:




  
    
      
        
        
          
            
              [
              
                
                  
                    
                      w
                      1
                    
                  
                  
                    
                      w
                      2
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      s
                      1
                    
                  
                
                
                  
                    
                      s
                      2
                    
                  
                
              
              ]
            
          
          ⁢
          
            d
            ·
            
              [
              
                
                  
                    
                      w
                      1
                    
                  
                  
                    
                      w
                      2
                    
                  
                
              
              ]
            
            ·
            
              [
              
                
                  
                    
                      α
                      1
                    
                  
                
                
                  
                    
                      α
                      2
                    
                  
                
              
              ]
            
          
        
        
      
      2
    
    
      
        
        
          b
          min
        
        
      
      2
    
  
  +
  
    
      
        
        
          
            
              [
              
                
                  
                    
                      w
                      2
                      *
                    
                  
                  
                    
                      -
                      
                        w
                        1
                        *
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      
                        s
                        1
                      
                      →
                    
                  
                
                
                  
                    
                      
                        s
                        2
                      
                      →
                    
                  
                
              
              ]
            
          
          -
          
            d
            ·
            
              
                [
                
                  
                    
                      
                        w
                        2
                        *
                      
                    
                    
                      
                        -
                        
                          w
                          1
                          *
                        
                      
                    
                  
                
                ]
              
              ⁡
              
                [
                
                  
                    
                      
                        α
                        1
                      
                    
                  
                  
                    
                      
                        α
                        2
                      
                    
                  
                
                ]
              
            
          
        
        
      
      2
    
    
      
        
        
          b
          max
        
        
      
      2
    
  
  +
  
     
    
      
        
          [
          
            
              
                
                  α
                  1
                  *
                
              
              
                
                  α
                  2
                  *
                
              
            
          
          ]
        
        ⁡
        
          [
          
            
              
                
                  G
                  B
                  
                    -
                    1
                  
                
              
              
                0
              
            
            
              
                0
              
              
                
                  G
                  B
                  
                    -
                    1
                  
                
              
            
          
          ]
        
      
      ·
      
        [
        
          
            
              
                α
                1
              
            
          
          
            
              
                α
                2
              
            
          
        
        ]
      
    
  




Avec bmin the noise component plus interference along the axis DIRmin and bmax the noise component plus interference along the axis DIRmax.
The resolution of this equation makes it possible to find the coefficients α1, α2.
In order to minimise the preceding expression, its gradient will simply be cancelled.
The following relationship is thus obtained which will make it possible to find α1 and α2:




  
    
      
        -
        
          (
          
            
              
                1
                
                  
                    
                    
                      b
                      min
                    
                    
                  
                  2
                
              
              ·
              
                [
                
                  
                    
                      
                        w
                        1
                        *
                      
                    
                  
                  
                    
                      
                        w
                        2
                        *
                      
                    
                  
                
                ]
              
              ·
              
                d
                H
              
              ·
              
                [
                
                  
                    
                      
                        w
                        1
                      
                    
                    
                      
                        w
                        2
                      
                    
                  
                
                ]
              
            
            +
            
              
                1
                
                  
                    
                    
                      b
                      max
                    
                    
                  
                  2
                
              
              ·
              
                [
                
                  
                    
                      
                        w
                        2
                      
                    
                  
                  
                    
                      
                        -
                        
                          w
                          1
                        
                      
                    
                  
                
                ]
              
              ·
              
                d
                H
              
              ·
              
                [
                
                  
                    
                      
                        w
                        2
                        *
                      
                    
                    
                      
                        -
                        
                          w
                          1
                          *
                        
                      
                    
                  
                
                ]
              
            
          
          )
        
      
      ⁢
      
        (
        
          
            [
            
              
                
                  
                    s
                    1
                  
                
              
              
                
                  
                    s
                    2
                  
                
              
            
            ]
          
          -
          
            d
            ·
            
              [
              
                
                  
                    
                      α
                      1
                    
                  
                
                
                  
                    
                      α
                      2
                    
                  
                
              
              ]
            
          
        
        )
      
    
    +
    
      
        [
        
          
            
              
                G
                B
                
                  -
                  1
                
              
            
            
              0
            
          
          
            
              0
            
            
              
                G
                B
                
                  -
                  1
                
              
            
          
        
        ]
      
      ·
      
        [
        
          
            
              
                α
                1
              
            
          
          
            
              
                α
                2
              
            
          
        
        ]
      
    
  
  =
  0




For this, in the same way as hereinabove, the following is posited:




  z
  =
  
    
      [
      
        
          
            
              z
              1
            
          
        
        
          
            
              z
              2
            
          
        
      
      ]
    
    =
    
      
        
          C
          2
          H
        
        ·
        
          [
          
            
              
                
                  α
                  1
                
              
            
            
              
                
                  α
                  2
                
              
            
          
          ]
        
      
      =
      
        
          [
          
            
              
                
                  C
                  B
                  H
                
              
              
                0
              
            
            
              
                0
              
              
                
                  C
                  B
                  H
                
              
            
          
          ]
        
        ·
        
          [
          
            
              
                
                  α
                  1
                
              
            
            
              
                
                  α
                  2
                
              
            
          
          ]
        
      
    
  




As such, the equation becomes:



  
    
      
        (
        
          
            
              1
              
                
                  
                  
                    b
                    min
                  
                  
                
                2
              
            
            ·
            
              [
              
                
                  
                    
                      w
                      1
                      *
                    
                  
                
                
                  
                    
                      w
                      2
                      *
                    
                  
                
              
              ]
            
            ·
            
              d
              H
            
            ·
            
              [
              
                
                  
                    
                      w
                      1
                    
                  
                  
                    
                      w
                      2
                    
                  
                
              
              ]
            
          
          +
          
            
              1
              
                
                  
                  
                    b
                    max
                  
                  
                
                2
              
            
            ·
            
              [
              
                
                  
                    
                      w
                      2
                    
                  
                
                
                  
                    
                      -
                      
                        w
                        1
                      
                    
                  
                
              
              ]
            
            ·
            
              d
              H
            
            ·
            
              [
              
                
                  
                    
                      w
                      2
                      *
                    
                  
                  
                    
                      -
                      
                        w
                        1
                        *
                      
                    
                  
                
              
              ]
            
          
        
        )
      
      ⁢
      
        (
        
          
            [
            
              
                
                  
                    s
                    1
                  
                
              
              
                
                  
                    s
                    2
                  
                
              
            
            ]
          
          -
          
            d
            ·
            
              [
              
                
                  
                    
                      C
                      B
                    
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    
                      C
                      B
                    
                  
                
              
              ]
            
            ·
            
              [
              
                
                  
                    
                      z
                      1
                    
                  
                
                
                  
                    
                      z
                      2
                    
                  
                
              
              ]
            
          
        
        )
      
    
    =
    
      
        [
        
          
            
              
                
                  C
                  B
                
                ,
                
                  Λ
                  B
                  
                    -
                    1
                  
                
              
            
            
              0
            
          
          
            
              0
            
            
              
                
                  C
                  B
                
                ,
                
                  Λ
                  B
                  
                    -
                    1
                  
                
              
            
          
        
        ]
      
      ·
      
        [
        
          
            
              
                z
                1
              
            
          
          
            
              
                z
                2
              
            
          
        
        ]
      
    
  




  
    
      
    
    ⁢
    
      
        with
        ⁢
        
          
            :
          
          ⁢
          
            

          
          ⁢
          
            
          
          [
          
            
              
                
                  α
                  1
                
              
            
            
              
                
                  α
                  2
                
              
            
          
          ]
        
      
      =
      
        
          [
          
            
              
                
                  C
                  B
                
              
              
                0
              
            
            
              
                0
              
              
                
                  C
                  B
                
              
            
          
          ]
        
        ·
        
          [
          
            
              
                
                  z
                  1
                
              
            
            
              
                
                  z
                  2
                
              
            
          
          ]
        
      
    
  



This therefore makes it possible to obtain α1 and α2.
All of the preceding remarks on the time/frequency separability of the covariance matrices GB, and therefore CB, as well as the reducing of these matrices with C′ can as such be applied.
Only the matrices GB and CB appear in the calculations. The matrices GD and CD no longer appear in the calculations.
As previously, the matrices GB,T , GB,F, then G′B,T G′B,F, are introduced again. All of the methods for resolving the system presented in the first resolution of Maximum A Posteriori can be applied in the same way (without the matrix GD).
According to an embodiment, the complexity of the calculation is further reduced by retained only certain eingenvalues of the covariance matrix G2 as well as the corresponding eigenvalues. In this case, a matrix C′ is extracted from the matrix C2 and a matrix Λ′ is extracted from the matrix Λ2, the matrix Λ′ comprising a determined number n′ of eigenvalues of Λ2 and the matrix C′ comprising the eigenvectors associated with these n′ eigenvalues. The weighting vector α is then determined using matrices of eigenvectors C′ and of eigenvalues Λ′.
In a preferred embodiment, the covariance matrix G2 is a diagonal matrix par bloc, constituted of the matrices GB, GB, and for which each one of the matrices GB is the Kronecker product of a covariance matrix GB,T, in the time domain and of a covariance matrix GB,F, in the frequency domain and the step of determining comprises the sub-steps of:

    
    
        decomposing of the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,TH and GD,F=CB,FΛB,FCB,FH, wherein:
        CB,T, is a matrix of eigenvectors of the matrix GB,T;
        ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
        CB,F is a matrix of eigenvectors of the matrix GB,F; and
        ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
        extracting of a matrix C′B,T from the matrix CB,T and a matrix Λ′B,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues; and
        extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
    
    


The weighting vector α is then determined using matrices of eigenvectors C′B,T, C′B,F, and of eigenvalues ΛB,T and ΛB,F.
The coefficients α1, α2 are then used in the step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals (r′, r′′).
Of course, the invention is not limited to the examples that have just been described and many arrangements can be made to these examples without leaving the scope of the invention.
In particular, the approach was described for two antennas but it can be generalised to a greater number of antennas.

## Claims

1. A method for demodulating a signal in a receiver comprising at least two antennas each receiving a signal transmitted through an associated radio propagation channel, the signals received corresponding to a same emitted signal comprising time and frequency distributed symbol frames in which driver symbols are known to the receiver, said method comprising:
noise whitening to form two combined signals that have noise components which are separate;
normalising the noise components of the two combined signals in order to form two signals (r′, r′′) that have noise components which are separate and have equal average standards;
performing signal demodulation with a maximum combination of a signal-to-noise ratio on the two signals (r′, r′′), the noise components of which are separate and which have equal average standards,
wherein the noise whitening to form the two combined signals comprises:
weighting each of the signals received with respective first weighting vectors associated with a respective antenna of the receiver, the signal associated with the first antenna being weighted by a weighting vector w1 and the signal associated with the second antenna being weighted by a weighting vector w2,
combining the weighted received signals in order to form a first combined signal (c1),
weighting a reference signal comprising said driver symbols with another weighting vector,
comparing the first combined signal and the weighted reference signal in order to form an error,
determining the weighting vectors (w1; w2) using a maximum a posteriori criterion by maximising a probability of occurring of said weighting vectors (w1; w2) conditionally with the error obtained,
weighting each signal received with: second weighting vectors (w*2, −w*1), the signal received on the first antenna being weighted by a conjugate of a complex vector w2 and the signal received on the second antenna being weighted by an opposite of the conjugate of a complex vector w1,
combining signals received weighted by the second weighting vectors in order to form a second combined signal (c2),
wherein respective coefficients of the propagation channel corresponding to a useful signal associated with the respective antenna of the receiver and used in performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals (r′, r′′), obtained at an end of the noise-whitening step, are determined using the maximum a posteriori criterion by maximising the probability of the occurring of the channel conditionally with the error present on the one hand in the first combined signal (c1) and on the other hand in the second combined signal (c2), by taking into account the propagation channel seen in the first combined signal (c1) and in the second combined signal (c2).

2. The method according to claim 1, wherein the determining comprises:
calculating a covariance matrix GB of the propagation channel;
calculating a covariance matrix GD of a product symbol by symbol of two propagation channels; and
determining weighting vectors using covariance matrices GB and GD.

3. The method according to claim 2, wherein the determining further comprises decomposing the covariance matrix GB into eigenvectors according to a relationship GB=CBΛBCBH wherein CB is a matrix of eigenvectors of the matrix GB and ΛB is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB and decomposing the covariance matrix GD into eigenvectors according to a relationship GD=CDΛDCDH wherein CD is a matrix of eigenvectors of the matrix GD and ΛD is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD.

4. The method according to claim 3, wherein the determining further comprises extracting a matrix C′B from the matrix CB and a matrix Λ′B from the matrix ΛB, the matrix Λ′B comprising a determined number nB′ of eigenvalues of ΛB and the matrix C′B comprising the eigenvectors associated with these nB′ eigenvalues, and extracting a matrix C′D from the matrix CD and a matrix Λ′D from the matrix ΛD, the matrix Λ′D comprising a determined number nD′ of eigenvalues of ΛD and the matrix C′D comprising the eigenvectors associated with the nD′ eigenvalues, the weighting vectors then being determined using matrices of eigenvectors C′B and C′D and of eigenvalues Λ′B and Λ′D.

5. The method according to claim 3, wherein the covariance matrix GB is a Kronecker product of a covariance matrix GB,T in a time domain and of a covariance matrix GB,F in a frequency domain and the covariance matrix GD is a Kronecker product of a covariance matrix GD,T in the time domain and of a covariance matrix GD,F in the frequency domain and the determining comprises:
decomposing the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,TH and GB,F=CB,FΛB,FCB,FH, wherein:
CB,T is a matrix of eigenvectors of the matrix GB,T;
ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
CB,F is a matrix of eigenvectors of the matrix GB,F and
ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
extracting a matrix C′B,T from the matrix CB,T and a matrix Λ′B,F from the matrix ΛB,T the matrix Λ′B,T comprising a determined number nBt′ of largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with said nBt′ largest eigenvalues;
extracting a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with said nBf′ largest eigenvalues,
decomposing the covariance matrices GD,T and GD,F into eigenvectors according to the relationships GD,T=CD,TΛD,TCD,TH and GD,F=CD,FΛD,FCD,FH, wherein:
CD,T is a matrix of eigenvectors of the matrix GD,T;
ΛD,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,T;
CD,F is a matrix of eigenvectors of the matrix GD,F; and
ΛD,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,F;
extracting a matrix C′D,T from the matrix CD,T and a matrix Λ′D,T from the matrix ΛD,T, the matrix Λ′D,T comprising a determined number nDt′ of largest eigenvalues of ΛD,T and the matrix C′D,T comprising the eigenvectors associated with said nDt′ largest eigenvalues; and
extracting a matrix C′D,F from the matrix CD,F and a matrix Λ′D,F from the matrix ΛD,F, the matrix Λ′D,F comprising a determined number nDf′ of the largest eigenvalues of ΛD,F and the matrix C′D,F comprising the eigenvectors associated with said nDf′ largest eigenvalues,
the weighting vectors then being determined using matrices of eigenvectors C′B,T, C′B,F, C′D,T, C′D,F and of eigenvalues Λ′B,T, Λ′B,F, Λ′D,T, Λ′D,F.

6. The method according to claim 1, wherein determining the signal demodulating comprises:
calculating a covariance matrix GB of the propagation channel;
determining weighting vectors using the covariance matrix GB.

7. The method according to claim 6, wherein the determining the signal demodulating further comprises decomposing the covariance matrix GB into eigenvectors according to a relationship GB=CBΛBCBH wherein CB is a matrix of eigenvectors of the matrix GB and ΛB is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB.

8. The method according to claim 7, wherein the determining the signal demodulating further comprises extracting a matrix C′B from the matrix CB and a matrix Λ′B from the matrix ΛB, the matrix Λ′B comprising a determined number nB′ of eigenvalues of ΛB and the matrix C′B comprising the eigenvectors associated with said nB′ eigenvalues, the weighting vectors then being determined using the matrix of eigenvectors C′B and of eigenvalues Λ′B.

9. The method according to claim 7, wherein the covariance matrix GB is a Kronecker product of a covariance matrix GB,T in a time domain and of a covariance matrix GB,F in a frequency domain and the determining the signal demodulating comprises:
decomposing covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,TH and GB,F=CB,FΛB,FCB,FH, wherein:
CB,T is a matrix of eigenvectors of the matrix GB,T;
ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
CB,F is a matrix of eigenvectors of the matrix GB,F; and
ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
extracting a matrix C′B,T from the matrix CB,T and a matrix Λ′B,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with said nBt′ largest eigenvalues;
extracting a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these said nBf′ largest eigenvalues,
the weighting vectors then being determined using matrices of eigenvectors C′B,T, C′B,F and of eigenvalues Λ′B,T, Λ′B,F.

10. The method according to claim 1, wherein the signal emitted is a multi-carrier signal.

11. The method according to claim 10, wherein the signal emitted is an OFDM signal.

12. A receiving equipment comprising at least two antennas each receiving a signal transmitted through an associated radio propagation channel, the signals received corresponding to a same emitted signal comprising time and frequency distributed symbol frames in which driver symbols are known to the receiving equipment, said method comprising means for:
noise whitening to form two combined signals that have noise components of which are separate;
normalising the noise components of the combined signals in order to form two signals (r′, r′′) that have noise components which are separate and have equal average standards;
performing signal demodulation with a maximum combination of the signal-to-noise ratio on the two signals (r′, r′′), the noise components of which are separate and which have equal average standards,
wherein the means of noise whitening to form the two combined signals comprises means for:
weighting each of the signals received with respective first weighting vectors (w1; w2) associated with a respective antenna of the receiving equipment, the signal associated with the first antenna being weighted by a weighting vector w1 and the signal associated with the second antenna being weighted by a weighting vector w2,
combining the weighted received signals in order to form a first combined signal (c1),
weighting a reference signal comprising said driver symbols with another weighting vector,
comparing the first combined signal and the weighted reference signal in order to form an error,
determining the weighting vectors (w1; w2) using a maximum a posteriori criterion by maximising a probability of occurring of said weighting vectors (w1; w2) conditionally with the error obtained,
weighting each signal received with second respective weighting vectors (w*2, −w*1), the signal received on the first antenna being weighted by a conjugate of a complex vector w2 and the signal received on the second antenna being weighted by an opposite of the conjugate of a complex vector w1,
combining the signals received weighted by the second weighting vectors in order to form a second combined signal,
wherein the receiving equipment further comprises means for determining respective coefficients of the propagation channel corresponding to a useful signal associated with the respective antenna of the receiving equipment and used in performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals (r′, r′′), obtained at an end of the noise-whitening step, are determined using the maximum a posteriori criterion by maximising the probability of the occurring of the channel conditionally with the error present on the one hand in the first combined signal c1 and on the other hand in the second combined signal c2, by taking into account the propagation channel seen in the first combined signal c1 and in the second combined signal c2.

13. A non-transitory computer readable medium comprising instructions for the implementation of the method according to claim 1 when the instructions are executed by at least one processor.