# Patent

## Title
Multi-antenna channel estimation method based on polyphase decomposition

## Abstract

A multi-antenna channel estimation method based on polyphase decomposition includes: receiving frequency domain received signals transformed using discrete fourier transformation (DFT) in pilot symbols; performing phase correction on the frequency domain received signals; performing polyphase decomposition on the frequency domain received signals which are corrected using phase correction and acquiring polyphase signals; performing interpolation on the polyphase signals and acquiring the estimation values of the multi-antenna channel parameters with various linear combinations on each frequency; acquiring decorrelation array based on the pilot structure of the transmission antenna and decorrelating the estimation values of the multi-antenna channel parameters with various linear combination on each frequency using the decorrelation array and acquiring channel parameters of the pilot symbols on each frequency; acquiring channel parameters of data symbols based on the channel parameters of the pilot symbols. With the present invention, the inversion problem in multi-antenna channel estimation is avoided with lower complexity, and results of the multi-antenna channel estimation method become more accurate since the DFT and filtering are performed in the interpolation after polyphase decomposition.

## Background
FIELD OF THE INVENTION
The present invention relates to the communication field, particularly, it relates to a channel estimation method.
BACKGROUND OF THE INVENTION
To meet the increasing demand of future mobile communication systems for system capacity, spectrum utilization, and data transmission rate, etc., in the existing LTE (Long Term Evolution) standard, MIMO (Multiple Input Multiple Output) technique and OFDM (Orthogonal Frequency Division Multiplexing) technique are employed. Wherein, the MIMO technique can double or even triple the system capacity and spectrum utilization without increasing the bandwidth. On the other hand, the OFDM technique is a typical technique applied in wideband transmission, and it has many advantages, including high resistance to multi-path fading, insensitivity to narrow-band interference and noise, and flexible bandwidth expansion, etc.
In wireless communication systems, the coherent detection technique is currently employed to achieve better receiving performance. Since channel parameters must be obtained to accomplish coherent detection, the channel estimation must be performed at the receiving end. In order to estimate the channel parameters timely and accurately, a pilot symbol assisted channel estimation method is currently employed in the actual communication systems. In MIMO systems, the multiple channel parameters must be estimated for each receiving channel; In OFDM systems, the channel parameter to be estimated usually indicates the frequency response parameter of channel.
There are a variety of channel estimation methods for MIMO systems which employ code division pilot on the same time-frequency resource. In the existing channel estimation methods, the accurate estimation of channel parameters always involves the complex matrix operations, such as large-dimension matrix inversion; such complex matrix operations will bring huge load on the entire communication system and degrade the processing speed of the system. To solve the above problems, the present invention proposes a channel estimation method, which utilizes polyphase extraction and decorrelation methods to reduce the order of complexity in the implementation of the MIMO channel estimation method, and the discrete cosine transform and filtering methods are employed for interpolation after polyphase extraction so as to improve the accuracy of the MIMO channel estimation method.
SUMMARY OF THE INVENTION
Technical Problem
The present invention provides a MIMO channel estimation method based on polyphase decomposition to overcome the drawbacks of the prior art, such as poor performance and high complexity of the MIMO channel estimation. The method provided in the present invention utilizes polyphase decomposition, interpolation, and decorrelation, etc., and avoids inversion operations in the MIMO channel estimation at a lower order of complexity, and thereby improves the performance of the MIMO channel estimation.
Technical Scheme
In the embodiments of the present invention, a MIMO channel estimation method based on polyphase decomposition is provided, and the method comprises:

    
    
        a. receiving frequency domain receiving signals at the pilot symbols processed through discrete Fourier transformation;
        b. performing phase correction for the frequency domain receiving signals;
        c. performing polyphase decomposition for the frequency domain receiving signals treated by the phase correction to obtain polyphase signals;
        d. performing interpolation for the polyphase signals to obtain estimated values of different linear combinations of the MIMO channel parameters at the frequency points;
        e. obtaining the decorrelation matrix according to the pilot structure of the transmitting antenna, and utilizing the decorrelation matrix to perform decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points to obtain the channel parameters of the pilot symbols at the frequency points;
        f. obtaining the channel parameters at the data symbols according to the channel parameters at the pilot symbols.
    
    


The procedure of phase correction for the frequency domain receiving signals comprises:
Performing the phase correction for the frequency domain receiving signals with a known fundamental sequence of pilot symbols, wherein, the amplitude of pilot symbols is normalized, and the vector of receiving signals treated by phase correction is {tilde over (y)}= XHy,
wherein, X=Diag{ x}, x represents the vector of the fundamental sequence of pilot symbols, Diag{•} represents a diagonal matrix that takes the vector in the curly brackets as the diagonal elements, ( )H represents the conjugate transposition of the matrix, y=[y1, y2, . . . , yNb]T represents the frequency domain receiving signals, yk represents the element k of the signal vector y, i.e., it represents the receiving signal at the frequency point k, k=1, . . . , Nb, Nb represents the number of pilot symbols and data symbols in the frequency domain direction, and ( )T represents the transposition of the matrix.
The procedure of performing polyphase decomposition for the receiving signals processed through the phase correction to obtain polyphase signals comprises:
The polyphase decomposition is to decompose the vector {tilde over (y)} of receiving signals after phase correction into P groups of polyphase signals, starting from the frequency point P before the vector with an extraction interval P, wherein, the polyphase signal vector in group p is sp,
Wherein, sp=[sp,1, sp,2, . . . , sp,N]T, sp,i={tilde over (y)}i′, i=1, 2, . . . , N, i′=(i−1)× P+p, P=P/GCD(n2, . . . , nNT, P), ni represents the cyclic time shift address of the transmitting antenna pilot i, P represents the number of the cyclic time shifts that are supported by the system at an equal interval, NT represents the number of transmitting antennae, and GCD(•) represents the greatest common divisor of all elements in the brackets.
The procedure of the interpolation for the polyphase signals comprises:
Performing spline interpolation separately for each group of the polyphase signals, or performing linear interpolation separately for each group of the polyphase signals, or performing interpolation separately for each group of polyphase signals by means of discrete cosine transformation and filtering so as to obtain the estimated value of linear combination of the MIMO channel parameters at the frequency points corresponding to each polyphase signal.
The procedure of interpolation by discrete cosine transformation and filtering comprises:

    
    
        5a. performing discrete cosine transformation for each group of the polyphase signal vectors: and the discrete cosine transformation utilizes type II discrete cosine transformation with the following transformation matrix:
    
    





  
    
      
        
          [
          
            C
            N
            II
          
          ]
        
        
          k
          ,
          l
        
      
      =
      
        
          w
          k
        
        ⁢
        cos
        ⁢
        
          
            
              π
              ⁡
              
                (
                
                  k
                  -
                  1
                
                )
              
            
            ⁢
            
              (
              
                l
                -
                0.5
              
              )
            
          
          N
        
        ⁢
        k
      
    
    ,
    
      l
      =
      1
    
    ,
    2
    ,
    ...
    ⁢
    
      
    
    ,
    N
  




  
    
      w
      k
    
    =
    
      {
      
        
          
            
              1
              /
              
                N
              
            
          
          
            
              k
              =
              1
            
          
        
        
          
            
              
                2
              
              /
              
                N
              
            
          
          
            
              k
              ≠
              1
            
          
        
      
    
  



Wherein, CNII is type II discrete cosine transformation matrix in length N, N is the number of frequency points contained in each group of the polyphase signals, k, l, and w are intermediate variables;
After performing the discrete cosine transformation for each group of the polyphase signal, the obtained signal vector in the discrete cosine transformation domain is spD=CNIIsp;

    
    
        5b. performing single-point filtering for the signal vector after transformation: performing single-point filtering for the signal vector in the discrete cosine transformation domain to obtained the filtered signal vector sp=ΓpspD,
    
    


Wherein, Γp is the single-point filtering matrix corresponding to the polyphase signal p, and it is a Nb×Nb diagonal matrix; if the statistical information of the signal vector in the discrete cosine transformation domain is known, the diagonal element k in the single-point filtering matrix Γp can be denoted as:




  
    
      
        [
        
          Γ
          p
        
        ]
      
      
        k
        ,
        k
      
    
    =
    
      
        
          E
          ⁢
          
            {
            
              
                
                
                  s
                  
                    p
                    ,
                    k
                  
                  D
                
                
              
              2
            
            }
          
        
        -
        
          σ
          z
          2
        
      
      
        E
        ⁢
        
          {
          
            
              
              
                s
                
                  p
                  ,
                  k
                
                D
              
              
            
            2
          
          }
        
      
    
  
  ,
  
    k
    =
    1
  
  ,
  2
  ,
  ...
  ⁢
  
    
  
  ,
  
    N
    b
  




Wherein, sp,kD represents the element k of vector spD, E{ } represent the operation of expectation, | |2 represents the square of the absolute value, σz2 represents noise variance, and Nb represents the number of pilot symbols and data symbols in the frequency domain direction;
If the statistical information of the signal vector in the discrete cosine transformation domain is unknown, the single-point filtering matrix Γp is the identity matrix of Nb×Nb;

    
    
        5c. performing extended inverse discrete cosine transformation and mapping for the filtered signal vector to obtain the estimated values of linear combination of the MIMO channel parameters corresponding to each polyphase signal at the frequency points; wherein, the extended inverse discrete cosine transformation and mapping for the filtered signal vector comprises:
    
    


The vector of estimated values of the linear combination of the MIMO channel parameters corresponding to the group p of polyphase signals at the frequency points can be denoted as




  
    
      
        s
        ~
      
      p
    
    =
    
      
        Λ
        p
      
      ⁢
      
        C
        _
      
      ⁢
      
        
          s
          _
        
        p
        D
      
    
  
  ,




wherein, Λp is the Nb×(Nb+ P−1) mapping matrix corresponding to polyphase signal p, and C represents the extended inverse discrete cosine transformation matrix,




  
    
      Λ
      p
    
    =
    
      
        Diag
        ⁢
        
          
            
              {
              
                
                  
                    0
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    0
                  
                  
                    }
                    
                      p
                      -
                      
                        1
                        ↑
                      
                    
                  
                
                ,
                
                  
                    1
                    ,
                    1
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    1
                  
                  
                    }
                    
                      
                        N
                        b
                      
                      ↑
                    
                  
                
                ,
                
                  
                    0
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    0
                  
                  
                    }
                    
                      
                        P
                        _
                      
                      -
                      
                        p
                        ↑
                      
                    
                  
                
              
              }
            
            ⁢
            
              

            
            [
            
              C
              _
            
            ]
          
          
            k
            ,
            l
          
        
      
      =
      
        
          w
          l
        
        ⁢
        
          cos
          (
          
            
              π
              ⁡
              
                (
                
                  l
                  -
                  1
                
                )
              
            
            ⁢
            
              (
              
                
                  
                    k
                    -
                    
                      P
                      _
                    
                  
                  
                    N
                    b
                  
                
                +
                
                  1
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    N
                  
                
              
              )
            
          
          )
        
      
    
  




  
    
      w
      l
    
    =
    
      {
      
        
          
            
              
                
                  
                    1
                    /
                    
                      N
                    
                  
                
                
                  
                    l
                    =
                    1
                  
                
              
              
                
                  
                    
                      2
                    
                    /
                    
                      N
                    
                  
                
                
                  
                    
                      l
                      ≠
                      1
                    
                    ,
                  
                
              
            
            ⁢
            
              

            
            ⁢
            l
          
          =
          1
        
        ,
        2
        ,
        ...
        ⁢
        
          
        
        ,
        N
        ,
        
          k
          =
          1
        
        ,
        2
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            N
            b
          
          +
          
            P
            _
          
          -
          1
        
        ,
      
    
  



Wherein, w, l, and k are intermediate variables, the estimated values of different linear combinations of the MIMO channel parameters at the frequency points can be denoted in the form of a matrix: {tilde over (S)}=[{tilde over (s)}1, {tilde over (s)}2, . . . , {tilde over (s)} P].
The procedure of obtaining the decorrelation matrix on the basis of the pilot structure of the transmitting antenna comprises:
Calculating the decorrelation matrix U according to the periodicity of the cyclic shift vector in the pilot structure of the transmitting antenna:




  U
  =
  
    
      
        1
        
          
            P
            _
          
        
      
      ⁡
      
        [
        
          
            I
            
              
                N
                T
              
              ×
              
                N
                T
              
            
          
          ,
          
            O
            
              
                N
                T
              
              ×
              
                P
                _
              
            
          
        
        ]
      
    
    ·
    
      
        [
        
          
            w
            1
          
          ,
          
            w
            2
          
          ,
          ...
          ⁢
          
            
          
          ,
          
            w
            
              N
              T
            
          
        
        ]
      
      *
    
  




Wherein, INT×NT represents the identity matrix of NT×NT, ONT× P represents the all-zero matrix of NT× P, wi represents the cyclic shift vector of the pilot symbols of the transmitting antenna i, and [ ]* represents conjugate operation of the matrix.
The procedure of performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points with a decorrelation matrix to obtain the channel parameters of the pilot symbols at the frequency points comprises:
Performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points with a decorrelation matrix U to obtain the estimated values of frequency domain response parameter in the channels of the transmitting antennae, wherein, the estimated value of the frequency domain response parameter in the channel of transmitting antenna i is ĥi={tilde over (S)}ui, wherein, ui represents the column i of the decorrelation matrix U.
The procedure of obtaining the channel parameters at the data symbols according to the channel parameters at the pilot symbols comprises: obtaining the frequency domain response parameter of the channel at the data symbols by linear interpolation on the basis of the frequency domain response parameter of the channel at the pilot symbols.
Beneficial Effects:
In the MIMO channel estimation method provided in the embodiments of the present invention, the frequency domain receiving signals at the receiving end are polyphase mixed signals due to the cyclic shift characteristic of the pilot in the transmitting antenna, and polyphase decomposition and decorrelation measures are utilized to decrease the order of complexity in the implementation of the MIMO channel estimation method. In addition, the discrete cosine transformation and filtering measures are employed for the interpolation after polyphase decomposition, therefore, the MIMO channel estimation method is more accurate.



BRIEF DESCRIPTION OF THE DRAWINGS
To make clear of the technical schemes in the embodiments of the present invention, the drawings used in the description of the embodiments or the prior art hereunder will be introduced briefly. Apparently, the drawings described below only illustrate some embodiments of the present invention. Those skilled in the art can obtain the drawings of other embodiments on the basis of these drawings without creative work.
FIG. 1 is the flow diagram of a channel estimation method provided in an embodiment of the present invention;
FIG. 2 is the flow diagram of another channel estimation method provided in an embodiment of the present invention.



DETAILED DESCRIPTION OF THE EMBODIMENTS
To get a better understand of the technical schemes of the present invention, the technical schemes in the embodiments of the present invention hereunder will be described clearly and completely with the reference to the accompanying drawings. Apparently, the embodiments described below are only some embodiments of the present invention, instead of all embodiments of the present invention. Those skilled in the art can obtain other embodiments without creative work on the basis of the embodiments provided here however, all these embodiments shall be deemed as falling into the protected scope of the present invention.
FIG. 1 shows the flow diagram of a MIMO channel estimation method provided in an embodiment of the present invention, wherein, the method comprises:
101. Receiving frequency domain receiving signals at the pilot symbols processed through discrete Fourier transformation.
In the embodiment of the present invention, the receiving signals at the pilot symbol on the receiving antenna can be treated to remove the cyclic prefix before they are processed through discrete Fourier transformation to obtain frequency domain receiving signals, and the channel estimator can perform further treatment for the frequency domain signals after it receives the frequency domain signals.
102. Performing phase correction for the frequency domain receiving signals.
In the embodiment of the present invention, the amplitude of pilot symbols is normalized, and the receiving end knows the fundamental sequence of the pilot symbols. Therefore, the fundamental sequence of the pilot symbols can be utilized to correct the phase for the receiving signals.
103. Performing polyphase decomposition for the frequency domain receiving signals treated by the phase correction to obtain polyphase signals.
In the embodiment of the present invention, the polyphase decomposition refers to extract several groups of signal vectors from the vector of the frequency domain receiving signals, starting from different frequency points with an equal interval; each group of extracted signal vectors is referred to a group of polyphase signals, and the groups of polyphase signals constitute the polyphase signals.
Decompose the receiving signal vector {tilde over (y)} of after phase correction into P groups of polyphase signals, starting from the frequency point P before the vector with an extraction interval P, wherein, the polyphase signal vector in group p is sp,
Wherein, sp=[sp,1, sp,2, . . . , Sp,N]T, sp,i={tilde over (y)}i′, i=1, 2, . . . , N, i′=(i−1)× P+p, P=P/GCD(n2, . . . , nNT, P), ni represents the cyclic time shift address of the transmitting antenna pilot i, P represents the number of cyclic time shifts that are supported by the system at an equal interval, NT represents the number of transmitting antennae, and GCD(•) represents the greatest common divisor of all elements in the brackets.
104. Performing interpolation for the polyphase signals to obtain estimated values of different linear combinations of the MIMO channel parameters at the frequency points.
In the embodiment of the present invention, the interpolation is performed separately for each group of the polyphase signals after the polyphase decomposition, and performing the interpolation by means of piecewise linear interpolation or spline interpolation to obtain the estimated values of linear combination of the MIMO channel parameters corresponding to the polyphase signals at the frequency points. Specifically, the vector of estimated values of linear combination of the MIMO channel parameters corresponding to the group p of polyphase signals at the frequency points can be denoted as {tilde over (s)}p=[{tilde over (s)}p,1, {tilde over (s)}p,2, . . . , {tilde over (s)}p,Nb]T, wherein, {tilde over (s)}p,k represents the estimated value of linear combination of the MIMO channel parameters corresponding to group p of polyphase signals at frequency point k; the estimated values of all linear combinations of the MIMO channel parameters at the frequency points can be denoted in the form of a matrix {tilde over (S)}=[{tilde over (s)}1, {tilde over (s)}2, . . . , {tilde over (s)} P].
105. Obtaining the decorrelation matrix according to the pilot structure of the transmitting antenna, and utilizing the decorrelation matrix to perform decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points to obtain the channel parameters of the pilot symbols at the frequency points.
In the embodiment of the present invention, the decorrelation matrix U can be calculated as follows:




  U
  =
  
    
      
        1
        
          
            P
            _
          
        
      
      ⁡
      
        [
        
          
            I
            
              
                N
                T
              
              ×
              
                N
                T
              
            
          
          ,
          
            O
            
              
                N
                T
              
              ×
              
                P
                _
              
            
          
        
        ]
      
    
    ·
    
      
        [
        
          
            w
            1
          
          ,
          
            w
            2
          
          ,
          ...
          ⁢
          
            
          
          ,
          
            w
            
              N
              T
            
          
        
        ]
      
      *
    
  




Wherein, INT×NT represents the identity matrix of NT×NT, ONT× P represents the all-zero matrix of NT× P, [ ]* represents conjugate operation of the matrix, and wi represents the cyclic shift vector of the pilot symbols of the transmitting antenna i. Performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points with a decorrelation matrix U to obtain the estimated values of the frequency domain response parameter in the channels of the transmitting antennae, wherein, the estimated value of the frequency domain response parameter in the channel of the transmitting antenna i is ĥi={tilde over (S)}ui, wherein, ui represents the column i of the decorrelation matrix U.
106. Obtaining the channel parameters at the data symbols according to the channel parameters at the pilot symbols.
In the MIMO channel estimation method provided in the embodiment of the present invention, the receiving signals at the receiving end is a polyphase mixed signal due to the cyclic shift characteristic of the pilot in the transmitting antenna, and polyphase decomposition, interpolation, and decorrelation measures are utilized. Therefore, the MIMO channel estimation method is more accurate and the order of complexity in implementation is lower.
FIG. 2 shows the flow diagram of another MIMO channel estimation method provided in another embodiment of the present invention. The method can be applied in the scenario of SD-OFDMA (Space Division-Orthogonal Frequency Division Multiple Access) system. In this embodiment, NT transmitting antennae and NR receiving antennae are considered. In this case, the transmission channels are MIMO (Multiple Input Multiple Output) channels, so it should be noted that the present invention is not limited to the scenario of MIMO, though the embodiment is described in an example of MIMO channels. Specifically, the channel estimation method can comprise:
201. Constructing a two-dimensional pilot data structure that covers time domain and frequency domain for the transmitting end;
In this embodiment, a two-dimensional pilot data structure is constructed at the transmitting end first. Specifically, a pilot sequence and the arrangement of the pilots can be constructed. In this embodiment, a block-type pilot can be used, wherein, the number of time domain symbols per frame or time slot is Ns, the number of pilot symbols in the time domain direction is Np, the length of frequency domain subcarrier is Nc, the number of frequency domain subcarriers occupied by pilot symbols and data symbols is Nb, and usually Nb<Nc. The fundamental sequence of pilot symbols at the transmitting end can be generated from a CAZAC (Constant Amplitude Zero Auto Correlation) sequence, and the pilot sequences of different transmitting antennae can be obtained from different cyclic shifts of the fundamental sequence, in addition, the pilot sequences of different transmitting antennae occupy the same time-frequency resource. Suppose the fundamental sequence of pilot symbol np in a frame is xnp, then the pilot sequences of all transmitting antennae at the pilot symbol are in the following cyclic shift relation with the fundamental sequence xnp of pilot symbols:




  
    
      
        {
        
          
            
              
                
                  
                    x
                    
                      i
                      ,
                      
                        n
                        p
                      
                    
                  
                  =
                  
                    
                      W
                      i
                    
                    ⁢
                    
                      
                        x
                        _
                      
                      
                        n
                        p
                      
                    
                  
                
              
            
            
              
                
                  
                    W
                    i
                  
                  =
                  
                    Diag
                    ⁢
                    
                      {
                      
                        w
                        i
                      
                      }
                    
                  
                
              
            
            
              
                
                  
                    w
                    i
                  
                  =
                  
                    
                      [
                      
                        1
                        ,
                        
                          e
                          jα
                        
                        ,
                        
                          e
                          
                            j2α
                            i
                          
                        
                        ,
                        
                          e
                          
                            
                              j
                              ⁡
                              
                                (
                                
                                  
                                    N
                                    b
                                  
                                  -
                                  1
                                
                                )
                              
                            
                            ⁢
                            
                              α
                              i
                            
                          
                        
                      
                      ]
                    
                    T
                  
                
              
            
            
              
                
                  
                    α
                    i
                  
                  =
                  
                    2
                    ⁢
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        n
                        i
                      
                      /
                      P
                    
                  
                
              
            
          
          ,
          
            i
            =
            1
          
          ,
          2
          ,
          ...
          ⁢
          
            
          
          ,
          
            N
            T
          
        
      
    
    
      
        (
        1
        )
      
    
  




Wherein, xi,np represents the pilot sequence of transmitting antenna i at pilot symbol np in a frame, Diag{•} represents a diagonal matrix that takes the vector in the curly brackets as the diagonal elements, and ni represents the cyclic time shift address of transmitting antenna i, the value range as follows:




  
    
      
        
          n
          i
        
        =
        
          {
          
            
              
                0
              
              
                
                  i
                  =
                  1
                
              
            
            
              
                
                  1
                  ,
                  2
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  
                    P
                    -
                    1
                  
                
              
              
                
                  
                    i
                    =
                    2
                  
                  ,
                  
                    ...
                    ⁢
                    
                      
                    
                    ⁢
                    
                      N
                      T
                    
                  
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




In addition, if i≠j, then ni≠nj is always true. Under that condition, the pilot sequences of different transmitting antennae are different to each other; P represents the number of cyclic time shifts that supported by the system at an equal interval, αi represents the cyclic shift factor of transmitting antenna i, wi represents the cyclic shift vector of pilot of transmitting antenna i, Wi represents the cyclic shift matrix of pilot of transmitting antenna i.
202. Performing discrete Fourier transformation for the receiving signals to obtain the frequency domain receiving signals.
For NR receiving antennae, there are NR receiving channels, and for each receiving channel, the frequency domain response parameter of the channel to be estimated is NT×Nb. At the receiving end, the receiving signal vector of the pilot symbol np, which is obtained by removing the cyclic prefix of the signal received by receiving channel r and performing discrete Fourier transformation for the signal, can be denoted as:




  
    
      
        
          y
          
            r
            ,
            
              n
              p
            
          
        
        =
        
          
            
              Ω
              T
            
            ⁢
            
              F
              N
            
            ⁢
            
              Y
              
                r
                ,
                
                  n
                  p
                
              
            
          
          =
          
            
              
                ∑
                
                  i
                  =
                  1
                
                
                  N
                  T
                
              
              ⁢
              
                
              
              ⁢
              
                Diag
                ⁢
                
                  
                    {
                    
                      x
                      
                        i
                        ,
                        
                          n
                          p
                        
                      
                    
                    }
                  
                  ·
                  
                    h
                    
                      r
                      ,
                      i
                      ,
                      
                        n
                        p
                      
                    
                  
                
              
            
            +
            
              z
              
                r
                ,
                
                  n
                  p
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Wherein, ΩT represents the inverse mapping matrix of subcarriers, FN represents the normalized DFT matrix of N×N, N represents the total number of frequency domain subcarriers, Yr,np represents the receiving signal vector of receiving antenna r at pilot symbol np after the cyclic prefix is removed, and yr,np=[y(r,np,1), y(r, np,2), . . . , y(r,np,Nb)]T, y(r,np,k) represents the receiving signals of receiving antenna r at frequency point k at pilot symbol np, hr,i,np represents the column vector composed of the frequency domain response parameter of the channel from transmitting antenna i to receiving antenna r at pilot symbol np, zr,np represents complex Gaussian additive white noise, with the presumption that the mean value of distribution parameter of the additive Gaussian noise is 0, and the variance of the distribution parameter is σz2. In view that the channel estimation is performed separately from each pilot symbol in each receiving channel, the label np and r can be omitted, and thus the formula 3 can be denoted as:




  
    
      
        y
        =
        
          
            
              ∑
              
                i
                =
                1
              
              
                N
                T
              
            
            ⁢
            
              
            
            ⁢
            
              
                X
                i
              
              ⁢
              
                h
                i
              
            
          
          +
          z
        
      
    
    
      
        (
        4
        )
      
    
  




Wherein, Xi=Diag{xi}, the formula (4) represents the basic model of receiving signals in a receiving channel of a pilot band.
203. Performing phase correction for the frequency domain receiving signals at the pilot symbols to obtain the receiving signals processed through the phase correction.
Sending the constructed pilot data structure to the receiving end so that both the transmitting end and the receiving end know the pilot data structure. Since the amplitude of the pilot symbols is normalized, the fundamental sequence of the pilot symbols can be utilized at the receiving end to perform phase correction for the frequency domain receiving signals. Specifically, the receiving signals at the pilot symbols processed through the phase correction can be denoted as:




  
    
      
        
          y
          ~
        
        =
        
          
            
              
                X
                _
              
              H
            
            ⁢
            y
          
          =
          
            
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  
                    N
                    T
                  
                
                ⁢
                
                  
                    
                      X
                      _
                    
                    H
                  
                  ⁢
                  
                    X
                    i
                  
                  ⁢
                  
                    h
                    i
                  
                
              
              +
              
                
                  
                    X
                    _
                  
                  H
                
                ⁢
                z
              
            
            =
            
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  
                    N
                    T
                  
                
                ⁢
                
                  
                    W
                    i
                  
                  ⁢
                  
                    h
                    i
                  
                
              
              +
              
                z
                ~
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




Wherein, X=Diag{ x}, x represents the vector of fundamental sequence of the pilot symbols, ( )H represents the conjugate transposition of the matrix, i.e., {tilde over (z)}= XHz. It can be seen from formula (5), that the elements of receiving signal vector {tilde over (y)} at each frequency point processed through the phase correction can be represented by the sum of noise items and a linear combination of the frequency domain response parameters of the channels of NT transmitting antennae at the corresponding frequency point and, whereas each linear combination is determined jointly by all cyclic shift vectors.
204. Performing polyphase decomposition for the receiving signals processed through the phase correction to obtain polyphase signals.
The elements of the receiving signal vector at each frequency point can be represented by the sum of the noise items and a linear combination of the frequency domain response parameter of the channels of all transmitting antennae at the corresponding frequency point. And, the linear combination at each frequency point is determined by all cyclic shift vectors. According to the matrix properties of Wi, the linear combination above is not irregular, instead, the vector Wi of the diagonal elements in wi has periodicity, with a period as P/GCD(ni,P), wherein, GCD(•) represents the greatest common divisor of all elements in the brackets. Therefore, the linear combination varies periodically over a period P=P/GCD (n2, . . . , nNT,P) as the variety of the frequency point.
The polyphase decomposition is to decompose the vector {tilde over (y)} of receiving signals after phase correction into P groups of polyphase signals, starting from the frequency point P before the vector with an extraction interval P, wherein, the polyphase signal vector in group p is sp, wherein, sp=[sp,1, sp,2, . . . , sp,N]T, sp,i={tilde over (y)}i′, i=1, 2, . . . , N, i′=(i−1)× P+p, and N=Nb/ P represents the number of frequency points of the polyphase signals in group p, i and i′ are intermediate variables.
205. Performing interpolation for the polyphase signals after polyphase decomposition to obtain the estimated values of different linear combinations of the MIMO channel parameters at the frequency points.
Interpolation is performed separately for each group of the polyphase signals after polyphase decomposition to obtain the estimated values of linear combinations of the MIMO channel parameters corresponding to the polyphase signals at the frequency points. Specifically, the vector sp of group p of the polyphase signals after the polyphase decomposition only comprises the estimated values of linear combination of the MIMO channel parameters corresponding to the group p of the polyphase signals at frequency points with ID in the collection S, wherein the collection S is defined as:

S={k|k=(p−1)mod  P+1,1≦k≦Nb},

Interpolation is performed according to the vector sp of group p of the polyphase signals to obtain the vector {tilde over (s)}p of estimated values of linear combination of the MIMO channel parameters corresponding to the group p of the polyphase signals at the rest frequency points. Here, the interpolation can be performed with any method that can be contrived easily, such as linear interpolation, spline interpolation, etc.; in this embodiment, an interpolation method based on discrete cosine transformation and filtering is provided:
Performing Type II discrete cosine transformation for the vector sp of group p of the polyphase signals after polyphase decomposition to obtain a signal vector spD=CNIIsp in the discrete cosine transformation domain, wherein, CNII represents the matrix of Type II discrete cosine transformation, and the elements in the matrix are denoted as follows:




  
    
      
        
          
            
              
                [
                
                  C
                  N
                  II
                
                ]
              
              
                k
                ,
                l
              
            
            =
            
              
                w
                k
              
              ⁢
              cos
              ⁢
              
                
                  
                    π
                    ⁡
                    
                      (
                      
                        k
                        -
                        1
                      
                      )
                    
                  
                  ⁢
                  
                    (
                    
                      l
                      -
                      0.5
                    
                    )
                  
                
                N
              
              ⁢
              
                
              
              ⁢
              k
            
          
          ,
          
            l
            =
            1
          
          ,
          2
          ,
          ...
          ⁢
          
            
          
          ,
          N
        
        ⁢
        
          

        
        ⁢
        
          
            w
            k
          
          =
          
            {
            
              
                
                  
                    1
                    /
                    
                      N
                    
                  
                
                
                  
                    k
                    =
                    1
                  
                
              
              
                
                  
                    
                      2
                    
                    /
                    
                      N
                    
                  
                
                
                  
                    k
                    ≠
                    1
                  
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




Perform filtering for the vector after the transformation. Specifically, single-point filtering can be performed for the signal vector in the discrete cosine transformation domain:

 sp=ΓpspD  (7)

Wherein, Γp is the single-point filtering matrix corresponding to polyphase signal p, and it is a diagonal matrix of Nb×Nb. In this embodiment, different filtering matrixes can also be selected as appropriate according to the actual circumstance:
If the statistical information of the signal vector is known, the diagonal element k in Γp can be denoted as:




  
    
      
        
          
            
              [
              
                Γ
                p
              
              ]
            
            
              k
              ,
              k
            
          
          =
          
            
              
                E
                ⁢
                
                  {
                  
                    
                      
                      
                        s
                        
                          p
                          ,
                          k
                        
                        D
                      
                      
                    
                    2
                  
                  }
                
              
              -
              
                σ
                z
                2
              
            
            
              E
              ⁢
              
                {
                
                  
                    
                    
                      s
                      
                        p
                        ,
                        k
                      
                      D
                    
                    
                  
                  2
                
                }
              
            
          
        
        ,
        
          k
          =
          1
        
        ,
        2
        ,
        ...
        ⁢
        
          
        
        ,
        
          N
          b
        
      
    
    
      
        (
        8
        )
      
    
  




Wherein, sp,kD represents the element k of vector spD.
If the statistical information of the channel vector is unknown, then, Γp is the identity matrix of Nb×Nb.
If the statistical information is known, the noise in the signals after transformation can be removed by single-point filtering. Of course, the filtering method in the embodiments of the present invention is not limited to single-point filtering, which is to say, other types of filtering can also be used. For example, the filtering based on the minimum mean square error criterion can be used, however, such filtering is multi-point filtering corresponding to single point filtering, and the filtering matrix involves matrix inversion operation, therefore, that filtering method is more complex.
Performing extended inverse discrete cosine transformation and mapping for the filtered signal vector to obtain a vector {tilde over (s)}p=ΛpC spD after discrete cosine interpolation, wherein, C represents the matrix of extended inverse discrete cosine transformation, and the elements in the matrix can be denoted as follows:




  
    
      
        
          
            
              [
              
                C
                _
              
              ]
            
            
              k
              ,
              l
            
          
          =
          
            
              w
              l
            
            ⁢
            
              cos
              ⁡
              
                (
                
                  
                    π
                    ⁡
                    
                      (
                      
                        l
                        -
                        1
                      
                      )
                    
                  
                  ⁢
                  
                    (
                    
                      
                        
                          k
                          -
                          
                            P
                            _
                          
                        
                        
                          N
                          b
                        
                      
                      +
                      
                        1
                        
                          2
                          ⁢
                          
                            
                          
                          ⁢
                          N
                        
                      
                    
                    )
                  
                
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            w
            l
          
          =
          
            {
            
              
                
                  
                    
                      
                        
                          1
                          /
                          
                            N
                          
                        
                      
                      
                        
                          l
                          =
                          1
                        
                      
                    
                    
                      
                        
                          
                            2
                          
                          /
                          
                            N
                          
                        
                      
                      
                        
                          l
                          ≠
                          
                            1
                            ′
                          
                        
                      
                    
                  
                  ⁢
                  
                    

                  
                  ⁢
                  l
                
                =
                1
              
              ,
              2
              ,
              ...
              ⁢
              
                
              
              ,
              N
              ,
              
                k
                =
                1
              
              ,
              2
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  N
                  b
                
                +
                
                  P
                  _
                
                +
                1
              
              ,
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




Wherein, the matrix Λp is the Nb×(Nb+ P−1) mapping matrix corresponding to polyphase signal p, and the elements in the matrix can be denoted as follows:




  
    
      
        
          Λ
          p
        
        =
        
          Diag
          ⁢
          
            {
            
              
                
                  0
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  0
                
                
                  }
                  
                    p
                    -
                    
                      1
                      ↑
                    
                  
                
              
              ,
              
                
                  1
                  ,
                  1
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  1
                
                
                  }
                  
                    
                      N
                      b
                    
                    ↑
                  
                
              
              ,
              
                
                  0
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  0
                
                
                  }
                  
                    
                      P
                      _
                    
                    -
                    
                      p
                      ↑
                    
                  
                
              
            
            }
          
        
      
    
    
      
        (
        10
        )
      
    
  




With the interpolation method described above, the vector of estimated values of linear combination of the MIMO channel parameters corresponding to group p of the polyphase signals at the frequency points can be obtained, and the vector can be denoted as {tilde over (s)}p. Then the estimated values of different linear combinations of the MIMO channel parameters at the frequency points can be denoted as matrix {tilde over (S)}=[{tilde over (s)}1, {tilde over (s)}2, . . . , {tilde over (s)} p].
206. Obtaining the decorrelation matrix according to the pilot structure of the transmitting antenna. Performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points with a decorrelation matrix to obtain the channel parameters of the pilot symbols at the frequency points.
Each row of the elements in matrix {tilde over (S)} are estimated values of different linear combinations of frequency domain response parameter of the MIMO channel at the corresponding frequency point, and the estimated values of different linear combinations of frequency domain response parameter of the MIMO channel can be decorrelated by means of a decorrelation matrix to recover the frequency domain response parameter of the MIMO channel.
The decorrelation matrix U can be calculated as follows:




  
    
      
        U
        =
        
          
            
              1
              
                
                  P
                  _
                
              
            
            ⁡
            
              [
              
                
                  I
                  
                    
                      N
                      T
                    
                    ×
                    
                      N
                      T
                    
                  
                
                ,
                
                  O
                  
                    
                      N
                      T
                    
                    ×
                    
                      P
                      _
                    
                  
                
              
              ]
            
          
          ·
          
            
              [
              
                
                  w
                  1
                
                ,
                
                  w
                  2
                
                ,
                ...
                ⁢
                
                  
                
                ,
                
                  w
                  
                    N
                    T
                  
                
              
              ]
            
            *
          
        
      
    
    
      
        (
        11
        )
      
    
  




Wherein, INT×NT represents the identity matrix of NT×NT, ONT× P represents the all-zero matrix of NT× P, and [ ]* represents conjugate operation of the matrix. Performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters with a decorrelation matrix U to obtain the estimated values of the frequency domain response parameter in the channels of the transmitting antennae, wherein, the estimated value of the frequency domain response parameter in the channel of transmitting antenna i is:

ĥi={tilde over (S)}ui  (12)

Wherein, ui represents the column i in the decorrelation matrix U.
207. Obtaining the channel parameters at the data symbols by interpolation or fitting according to the channel parameters at the pilot symbols.
For example, in this procedure, the channel parameters at the data symbols can be obtained by linear interpolation:
Suppose two adjacent pilots are at pilot symbol n1 and n2 in a frame, the frequency domain response parameter of channel of transmitting antenna i at the pilot symbols is denoted as ĥi,n2 and ĥi,n2, respectively. Then the estimated value of frequency domain response parameter of channel of the transmitting antenna at data symbol n can be denoted as:




  
    
      
        
          
            h
            ^
          
          
            i
            ,
            n
          
        
        =
        
          
            
              
                
                  
                    h
                    ^
                  
                  
                    i
                    ,
                    
                      n
                      2
                    
                  
                
                -
                
                  
                    h
                    ^
                  
                  
                    i
                    ,
                    
                      n
                      1
                    
                  
                
              
              
                
                  n
                  2
                
                -
                
                  n
                  1
                
              
            
            ⁢
            
              (
              
                
                  n
                  1
                
                -
                n
              
              )
            
          
          +
          
            
              h
              ^
            
            
              i
              ,
              
                n
                1
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




In addition, there are other implementations for obtaining the frequency domain response parameter of the channel at data symbols according to the frequency domain response parameter of the channel at pilot symbols. For example, mean value substitution method can be used.
In the MIMO channel estimation method provided in the embodiments of the present invention, the receiving signals at the receiving end are polyphase mixed signals due to the cyclic shift characteristic of the pilot in the transmitting antenna. Polyphase decomposition and decorrelation measures are utilized to decrease the order of complexity in the implementation of the MIMO channel estimation method. In addition, the discrete cosine transformation and filtering measures are employed for the interpolation after polyphase decomposition, therefore, the MIMO channel estimation method is more accurate.
The embodiments of the present invention can be implemented by software, and the software program can be stored in a readable storage media, such as a hard disk, cache, or CD in a computer.
While the present invention has been illustrated and described with reference to some preferred embodiments, the present invention is not limited to these. Those skilled in the art should recognize that various variations and modifications can be made without departing from the spirit and scope of the present invention as defined by the accompanying claims.

## Claims

1. A MIMO channel estimation method based on polyphase decomposition, comprising:
a. Receiving frequency domain receiving signals at the pilot symbols processed through discrete Fourier transformation;
b. Performing phase correction for the frequency domain receiving signals;
c. Performing polyphase decomposition for the frequency domain receiving signals treated by the phase correction to obtain polyphase signals;
d. Performing interpolation for the polyphase signals to obtain estimated values of different linear combinations of the MIMO channel parameters at the frequency points;
e. Obtaining a decorrelation matrix according to the pilot structure of a transmitting antenna, and utilizing the decorrelation matrix to perform decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points to obtain the channel parameters of the pilot symbols at the frequency points;
f. Obtaining the channel parameters at the data symbols according to the channel parameters at the pilot symbols.

2. The MIMO channel estimation method based on polyphase decomposition according to claim 1, wherein, the procedure of phase correction for the frequency domain receiving signals comprises:
Performing the phase correction for the frequency domain receiving signals with a known fundamental sequence of the pilot symbols, wherein, the amplitude of pilot symbol is normalized, and the vector of the receiving signals treated by phase correction is {tilde over (y)}= XHy,
wherein, X=Diag{ x}, x represents the vector of fundamental sequence of the pilot symbols, Diag{•} represents a diagonal matrix that takes the vector in the curly brackets as the diagonal elements, ( )H represents the conjugate transposition of the matrix, y=[y1, y2, . . . , yNb]T represents the frequency domain receiving signals, yk represents the element k of the signal vector y, i.e., it represents the receiving signal at the frequency point k, and k=1, . . . , Nb, Nb represents the number of pilot symbols and data symbols in the frequency domain direction, and ( )T represents the transposition of the matrix.

3. The MIMO channel estimation method based on polyphase decomposition according to claim 1, wherein, the procedure of performing polyphase decomposition for the receiving signals processed through the phase correction to obtain polyphase signals comprises:
The polyphase decomposition is to decompose the vector {tilde over (y)} of receiving signals after phase correction into P groups of polyphase signals, starting from the frequency point P before the vector with an extraction interval P, wherein, the polyphase signal vector in group p is sp,
wherein, sp=[sp,1, sp,2, . . . , sp,N]T, sp,i={tilde over (y)}i′, i=1, 2, . . . , N, i′=(i−1)× P+p, P=P/GCD(n2, . . . , nNT, P), ni represents the cyclic time shift address of the transmitting antenna pilot i, P represents the number of cyclic time shifts that are supported by the system at an equal interval, NT represents the number of transmitting antennae, and GCD(•) represents the greatest common divisor of all elements in the brackets.

4. The MIMO channel estimation method based on polyphase decomposition according to claim 1, wherein, the procedure of interpolation for the polyphase signals comprises:
Performing spline interpolation separately for each group of the polyphase signals, or
performing linear interpolation separately for each group of the polyphase signals, or
performing interpolation separately for each group of polyphase signals by means of discrete cosine transformation and filtering to obtain the estimated value of linear combination of the MIMO channel parameters at the frequency points corresponding to each polyphase signal.

5. The MIMO channel estimation method based on polyphase decomposition according to claim 4, wherein, the procedure of interpolation by discrete cosine transformation and filtering comprises:
5a. performing discrete cosine transformation for each group of polyphase signal vectors: the discrete cosine transformation utilizes type II discrete cosine transformation with the following transformation matrix:




  
    
      
        
          [
          
            C
            N
            II
          
          ]
        
        
          k
          ,
          l
        
      
      =
      
        
          w
          k
        
        ⁢
        cos
        ⁢
        
          
            
              π
              ⁡
              
                (
                
                  k
                  -
                  1
                
                )
              
            
            ⁢
            
              (
              
                l
                -
                0.5
              
              )
            
          
          N
        
        ⁢
        
          
        
        ⁢
        k
      
    
    ,
    
      l
      =
      1
    
    ,
    2
    ,
    ...
    ⁢
    
      
    
    ,
    N
  




  
    
      w
      k
    
    =
    
      {
      
        
          
            
              1
              /
              
                N
              
            
          
          
            
              k
              =
              1
            
          
        
        
          
            
              
                2
              
              /
              
                N
              
            
          
          
            
              k
              ≠
              1
            
          
        
      
    
  


Wherein, CNII is type II discrete cosine transformation matrix in length N, N is the number of frequency points contained in each group of the polyphase signals, k, l, and w are intermediate variables;
After performing the discrete cosine transformation for each group of the polyphase signal, the obtained signal vector in the discrete cosine transformation domain is spD=CNIIsp;
5b. performing single-point filtering for the signal vector after transformation: performing single-point filtering for the signal vector in the discrete cosine transformation domain to obtain the filtered signal vector sp=ΓpspD,
Wherein, Γp is the single-point filtering matrix corresponding to the polyphase signal p, and it is a Nb×Nb diagonal matrix; if the statistical information of the signal vector in the discrete cosine transformation domain is known, the diagonal element k in the single-point filtering matrix Γp can be denoted as:





  
    
      
        [
        
          Γ
          p
        
        ]
      
      
        k
        ,
        k
      
    
    =
    
      
        
          E
          ⁢
          
            {
            
              
                
                
                  s
                  
                    p
                    ,
                    k
                  
                  D
                
                
              
              2
            
            }
          
        
        -
        
          σ
          z
          2
        
      
      
        E
        ⁢
        
          {
          
            
              
              
                s
                
                  p
                  ,
                  k
                
                D
              
              
            
            2
          
          }
        
      
    
  
  ,
  
    k
    =
    1
  
  ,
  2
  ,
  ...
  ⁢
  
    
  
  ,
  
    N
    b
  



Wherein sp,kD represents the element k of vector spD, E{ } represent the operation of expectation, | |2 represents the square of the absolute value, σz2 represents noise variance, and Nb represents the number of pilot symbols and data symbols in the frequency domain direction;
If the statistical information of the signal vector in the discrete cosine transformation domain is unknown, the single-point filtering matrix Γp is the identity matrix of Nb×Nb;
5c. performing extended inverse discrete cosine transformation and mapping for the filtered signal vector to obtain the estimated values of linear combination of the MIMO channel parameters corresponding to each polyphase signal at the frequency points;
wherein, the extended inverse discrete cosine transformation and mapping for the filtered signal vector comprise:
The vector of estimated values of the linear combination of the MIMO channel parameters corresponding to the group p of polyphase signals at the frequency points can be denoted as {tilde over (s)}p=ApC spD, wherein, Λp is the Nb×(Nb+ P−1) mapping matrix corresponding to polyphase signal p, and C represents the extended inverse discrete cosine transformation matrix,




  
    
      Λ
      p
    
    =
    
      
        Diag
        ⁢
        
          
            
              {
              
                
                  
                    0
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    0
                  
                  
                    }
                    
                      p
                      -
                      
                        1
                        ↑
                      
                    
                  
                
                ,
                
                  
                    1
                    ,
                    1
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    1
                  
                  
                    }
                    
                      
                        N
                        b
                      
                      ↑
                    
                  
                
                ,
                
                  
                    0
                    ,
                    ...
                    ⁢
                    
                      
                    
                    ,
                    0
                  
                  
                    }
                    
                      
                        P
                        _
                      
                      -
                      
                        p
                        ↑
                      
                    
                  
                
              
              }
            
            ⁢
            
              

            
            [
            
              C
              _
            
            ]
          
          
            k
            ,
            l
          
        
      
      =
      
        
          w
          l
        
        ⁢
        
          cos
          ⁡
          
            (
            
              
                π
                ⁡
                
                  (
                  
                    l
                    -
                    1
                  
                  )
                
              
              ⁢
              
                (
                
                  
                    
                      k
                      -
                      
                        P
                        _
                      
                    
                    
                      N
                      b
                    
                  
                  +
                  
                    1
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      N
                    
                  
                
                )
              
            
            )
          
        
      
    
  




  
    
      w
      l
    
    =
    
      {
      
        
          
            
              
                
                  
                    1
                    /
                    
                      N
                    
                  
                
                
                  
                    l
                    =
                    1
                  
                
              
              
                
                  
                    
                      2
                    
                    /
                    
                      N
                    
                  
                
                
                  
                    l
                    ≠
                    
                      1
                      ′
                    
                  
                
              
            
            ⁢
            
              

            
            ⁢
            l
          
          =
          1
        
        ,
        2
        ,
        ...
        ⁢
        
          
        
        ,
        N
        ,
        
          k
          =
          1
        
        ,
        2
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            N
            b
          
          +
          
            P
            _
          
          -
          1
        
        ,
      
    
  


Wherein, w, l, and k are intermediate variables, the estimated values of different linear combinations of the MIMO channel parameters at the frequency points can be denoted in the form of a matrix: {tilde over (S)}=[{tilde over (s)}1, {tilde over (s)}2, . . . , {tilde over (s)} p].

6. The MIMO channel estimation method based on polyphase decomposition according to claim 1, wherein, the procedure of obtaining a decorrelation matrix according to the pilot structure of the transmitting antenna comprises:
Calculating the decorrelation matrix U according to the periodicity of the cyclic shift vector in the pilot structure of the transmitting antenna:





  U
  =
  
    
      
        1
        
          
            P
            _
          
        
      
      ⁡
      
        [
        
          
            I
            
              
                N
                T
              
              ×
              
                N
                T
              
            
          
          ,
          
            O
            
              
                N
                T
              
              ×
              
                P
                _
              
            
          
        
        ]
      
    
    ·
    
      
        [
        
          
            w
            1
          
          ,
          
            w
            2
          
          ,
          ...
          ⁢
          
            
          
          ,
          
            w
            
              N
              T
            
          
        
        ]
      
      *
    
  



Wherein, INT×NT represents the identity matrix of NT×NT, ONT× p represents the all-zero matrix of NT× P, wi represents the cyclic shift vector of the pilot symbols of the transmitting antenna i, and [ ]* represents conjugate operation of the matrix.

7. The MIMO channel estimation method based on polyphase decomposition according to claim 1, wherein, the procedure of performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points with a decorrelation matrix to obtain the channel parameters of the pilot symbols at the frequency points comprises:
Performing decorrelation for the estimated values of different linear combinations of the MIMO channel parameters at the frequency points with a decorrelation matrix U to obtain the estimated values of frequency domain response parameter in the channels of the transmitting antennae, wherein, the estimated value of the frequency domain response parameter in the channel of transmitting antenna i is ĥi={tilde over (S)}ui, where, ui represents the column i of the decorrelation matrix U.

8. The MIMO channel estimation method based on polyphase decomposition according to claim 1, wherein, the procedure of obtaining the channel parameters at the data symbols according to the channel parameters at the pilot symbols comprises: obtaining the frequency domain response parameter of the channel at the data symbols by linear interpolation on the basis of the frequency domain response parameter of the channel at the pilot symbols.