Abstract:
In prior art, selection of precoding granularity is affected by two inter-restricting factors, i.e., precoding accuracy and channel estimation at a mobile terminal. To solve this problem, the present invention provides a method of maintaining coherency of a precoding channel in a communication network and an associated apparatus. During precoding, this method takes into account both channel coherency and system capacity. A base station adjusts phase and/or amplitude of a precoding matrix corresponding to each precoded unit to maintain coherency of associated information of the overall precoding channel. The associated information of the precoding channel includes, for example, CSI or eigenvalue matrix of the precoding channel. Afterwards, a mobile terminal performs channel estimation based on reference signals of multiple precoded units, thereby eliminating the limitation in prior art that a mobile terminal can perform channel estimation only within one or more resource block limited by a precoding granularity.

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention relates to a wireless communication system, and more particularly, to a method of processing a precoding matrix in a base station and associated apparatus. 
       BACKGROUND OF THE INVENTION 
       [0002]    Non-codebook based precoding is a promising technology for a time division duplex (TDD) long term evolution-advanced (LTE-A) system, due to inherent channel reciprocity, namely symmetry between an uplink frequency and a downlink frequency. The industry widely accepts the assumption of uplink-downlink reciprocity and channel estimation is effectively performed based on this assumption. 
         [0003]    In the case of non-codebook based precoding, the precoding matrix is obtained at the transmitting end. The transmitting end utilizes predicated channel status information (CSI) to calculate the precoding matrix. Common methods of calculating a precoding matrix include singular value decomposition (SVD), uniform channel decomposition (UCD) and the QR algorithm. 
         [0004]      FIG. 1  is a structural diagram illustrating a transmitter and a receiver in a SVD-based multiple input multiple output (MIMO) system. Given that the transmitter of a base station  1  has N antennas and the receiver of a mobile terminal  2  has M antennas, the dimension of the effective uplink CSI, i.e., the spatial channel matrix, is M×N, and the spatial channel matrix may be denoted as H M×N . H M×N  is processed according to SVD shown in formula (1): 
         [0000]      H=UDV H    (1)
 
         [0000]    where U and V are the left singular vector matrix and the right singular vector matrix of H, respectively. Both U and V are a unitary matrix, i.e., it follows UU H =I=VV H , where I is an identity matrix and (•) H  denotes hermitian operation resulting in a transposed complex conjugate. It follows that U ∈C N×N , i.e., the dimension of U is N×N, and that V ∈C M×M , i.e., the dimension of V is M×M. The rank r of the CSI matrix H satisfies r≦min(M,N). The diagonal matrix D can be represented as 
         [0000]    
       
         
           
             
               D 
               = 
               
                 
                   [ 
                   
                     
                       
                         
                           D 
                           r 
                         
                       
                       
                         0 
                       
                     
                     
                       
                         0 
                       
                       
                         0 
                       
                     
                   
                   ] 
                 
                 
                   N 
                   × 
                   M 
                 
               
             
             , 
           
         
       
     
         [0000]    where D r =diag(λ 1 ,λ 2 , . . . λ r ) and λ i  is singular value of H in descending order, i.e., λ 1 &gt;λ 2  &gt; . . . λ r . 
         [0005]    The right singular vector matrix V obtained after SVD is a linear precoding matrix. Each column of V is called an eigenvector of H H H, which is related to the eigenmode of the communication channel. If rank self-adaption is required, column vectors corresponding to greater singular values are selected from the right singular vector matrix V to compose the precoding matrix. 
         [0006]    Non-codebook based precoding requires a dedicated pilot, which means data symbols and pilot symbols are precoded together. Thus the receiving end can obtain the effective channel after precoding by only performing channel estimation, thereby facilitating data modulation. 
         [0007]    Since accurate CSI can be obtained due to reciprocity between the uplink and the downlink in a TDD system, non-codebook based precoding can provide additional precoding gain. Generally and theoretically, the smaller a precoding granularity is, the higher the corresponding precoding gain becomes. Precoding granularity is defined as a unit to be precoded, such as one or more resource block (RB).  FIG. 2  illustrates precoding performance corresponding to different precoding granularities in a single-layer beamforming case. As shown, given equal signal to noise ratio (SNR), the smaller the precoding granularity is utilized, the greater the system throughput becomes. A precoding granularity of 10 means that 10 RBs employ the same precoding matrix. However, there is difference between channel responses corresponding to the 10 RBs. Therefore, the greater a precoding granularity is set, the less accurately a precoding matrix weighting all the RBs in the precoding granularity matches the actual channel status of each RB in the precoded unit. Consequently, taking into account matching of a precoding matrix and a channel, it is desirable to have a smaller precoding granularity to obtain a greater precoding gain. 
         [0008]    However, in a practical system, performance gain of precoding is affected by channel estimation error (3GPP R1-092794). Since a smaller precoding granularity utilizes a lower reference signal (RS) power, accuracy of channel estimation is reduced. Therefore, an appropriately selected precoding granularity inevitably affects the system capacity. Furthermore, selection of precoding granularity is also an important issue for a multiple-user multiple input multiple output (MU-MIMO) system or a coordinative multiple point (CoMP) system sensitive to different multi-path delays between different user equipments or cells. 
         [0009]    The above conclusion relies on the fact that channel estimation can be performed only within a precoding granularity. This is because each precoding granularity corresponds to a different precoding matrix and a different precoding matrix will corrupt channel coherency among multiple precoded units. Therefore, taking precoding accuracy into account, it is desirable to have a smaller precoding granularity. On the other hand, taking channel estimation into account, it is desirable to have a greater precoding granularity. Consequently, these two factors inter-restrict. 
         [0010]    In a prior art solution, a base station dynamically monitors the channel to obtain its real-time status and then selects a corresponding precoding granularity according to information such as channel coherency, signal to interference and noise ratio (SINR), etc. Then, the base station transmits the selected precoding granularity to a mobile terminal. The mobile terminal performs channel estimation according to this indication within the resource blocks limited by the precoding granularity. The terminal needs to be notified of such indication real-time. Therefore, a lot of time-frequency resources are occupied. 
       SUMMARY OF THE INVENTION 
       [0011]    The present invention analyzes channel coherency among multiple precoded units. The so called coherency is statistical characteristics of a channel, i.e., frequency selective characteristics and time variant characteristics of the channel. Furthermore, the present invention provides a method of coherent time frequency precoding (CTFP). This method takes into account both channel coherency and system capacity. A base station (eNB) adjusts phase and/or amplitude of a precoding matrix corresponding to each precoded unit to maintain coherency of associated information of the overall precoding channel. The associated information of the precoding channel includes, for example. CSI or eigenvalue matrix of the precoding channel. Afterwards, a mobile terminal performs channel estimation based on reference signals of multiple precoded units, thereby eliminating the limitation in prior art that a mobile terminal can perform channel estimation only within one or more resource block limited by a precoding granularity. In other words, the base station can employ a precoding granularity as small as possible without affecting channel estimation at the side of the mobile terminal. 
         [0012]    According to a first aspect of the present invention, there is provided a method of transmitting a pilot and/or data weighted by a precoding matrix in a base station of a wireless communication system, wherein the base station obtains channel status information. The method includes: performing matrix decomposition to the channel status information to obtain an initial precoding matrix, wherein the initial precoding matrix is not unique; performing linear transformation to the initial precoding matrix such that information associated with an corrected precoding channel under the condition of weighting of the linearly transformed precoding matrix maintains coherency; transmitting to a mobile terminal the pilot and/or data weighted by the linearly transformed precoding matrix. 
         [0013]    According to a second aspect of the present invention, there is provided a processing apparatus for transmitting a pilot and/or data weighted by a precoding matrix in a base station of a wireless communication system, wherein the base station obtains channel status information. The apparatus includes: an initial precoding matrix obtaining means for performing matrix decomposition to the channel status information to obtain an initial precoding matrix, wherein the initial precoding matrix is not unique; an correcting means for performing linear transformation to the initial precoding matrix such that information associated with an corrected precoding channel under the condition of weighting of the linearly transformed precoding matrix maintains coherency; a transmitting means for transmitting to a mobile terminal the pilot and/or data weighted by the linearly transformed precoding matrix. 
         [0014]    The solution of the present invention enables a base station to perform precoding within a precoded unit of smaller precoding granularity, thereby increasing precoding gain. Moreover, it enables a mobile terminal to perform channel estimation based on reference signals of multiple precoded units and thus eliminates the limitation in prior art that a mobile terminal can perform channel estimation only within one or more time frequency resource block limited by a precoding granularity, thereby improving channel estimation performance of the mobile terminal. In addition, the base station doesn&#39;t need to provide the mobile terminal with indication about the precoding granularity, thereby saving corresponding signaling overhead. In a MU-MIMO and CoMP system, this solution enables selection of optimum precoding granularity between a base station and a user without considering channel characteristics of a coordinate base station or terminal. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0015]    Other features, objectives and advantages of the present invention will become more apparent by studying the following detailed description of non-limiting embodiments, with reference to the accompanying drawings, wherein below: 
           [0016]      FIG. 1  is a structural diagram illustrating a transmitter and a receiver in a SVD-based MIMO system. 
           [0017]      FIG. 2  illustrates the precoding performance corresponding to different precoding granularities in a single-layer beamforming case. 
           [0018]      FIG. 3  illustrates absolute value of amplitude of the channel after the precoding of an initial precoding matrix. 
           [0019]      FIG. 4  illustrates phase of the channel after the precoding of an initial precoding matrix. 
           [0020]      FIG. 5  is a flowchart illustrating a method according to an embodiment of the present invention. 
           [0021]      FIG. 6  illustrates a channel curve with respect to that of  FIG. 4 , after phase rotation 
           [0022]      FIGS. 7 and 8  illustrate a channel curve of a time domain response obtained by domain transforming to the channel of  FIG. 6 . 
           [0023]      FIG. 9  is a flowchart illustrating a method according to another embodiment of the present invention. 
           [0024]      FIG. 10  is a block diagram illustrating an apparatus according to an embodiment of the present invention. 
           [0025]      FIG. 11  is a block diagram illustrating an apparatus according to another embodiment of the present invention, 
       
    
    
       [0026]    Identical or similar reference signs in the drawings represent identical or similar components. 
       DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
       [0027]    The following formula can be obtained by performing singular value decomposition to a MIMO channel matrix: 
         [0000]      H i =U i D i V i   H    (2)
 
         [0000]    where i is the index of a precoded unit, D i  is an eigenvalue matrix, i.e., singular value matrix, and V i   H  is the hermitian transformation of V i . As is well-known, the weighting matrices U i  and V i  of SVD are not unique. For example, formula (1) remains equal if the first columns or the identical columns of both U i  and V i   H  are rotated by π±2kπ. 
         [0028]    Taking a case of two transmitting antennas (TX) and two receiving antennas (RX) for example,  FIGS. 3 and 4  are a frequency domain expression of an initial precoding channel (each sub-carrier) after the precoding of an initial precoding matrix, for a typical spatial channel model (SCM).  FIG. 3  illustrates absolute value of amplitude of the initial precoding channel after the precoding of an initial precoding matrix, and  FIG. 4  illustrates the phase of the initial precoding channel after the precoding of an initial precoding matrix. 
         [0029]    In  FIGS. 3 and 4 , h 11  represents the impulse response of a channel from an antenna TX 1  to an antenna RX 1 , h 12  represents the impulse response of a channel from the antenna TX 1  to an antenna RX 2 , h 21  represents the impulse response of a channel from an antenna TX 2  to the antenna RX 1 , and h 22  represents the impulse response of a channel from the antenna TX 2  to the antenna RX 2 . The channel matrix is 
         [0000]    
       
         
           
             
               H 
               i 
             
             = 
             
               
                 [ 
                 
                   
                     
                       
                         h 
                         11 
                       
                     
                     
                       
                         h 
                         12 
                       
                     
                   
                   
                     
                       
                         h 
                         21 
                       
                     
                     
                       
                         h 
                         22 
                       
                     
                   
                 
                 ] 
               
               . 
             
           
         
       
     
         [0000]    In this embodiment, the matrix V i  obtained after performing SVD to the channel matrix based on formula (1) serves as the initial precoding matrix of a channel corresponding to each precoded unit. Therefore, the initial precoding channel is H ip0 =H i V i =U i D i V i   H V i =U i D i , where i=0, 1, 2, 3 . . . 120 representing the index of a precoding granularity, p0 represents an initial precoding channel. 
         [0030]    In  FIGS. 3 and 4 , a DFT-based channel estimation algorithm is utilized to estimate channel parameters. As easily seen, absolute value of amplitude of the initial precoding channel is coherent/smooth, but there are some jump points or jump segment on the phase curve. Moreover, one column (e.g., an eigenvector) of the initial preceding channel has the same change segment. 
         [0031]    Note that in  FIGS. 3 and 4 , for purposes of illustration, each sub-carrier utilizes a precoding sub-channel corresponding to a different SVD. In other words, the precoding granularity shown in  FIGS. 3 and 4  is one sub-carrier. For precoding of each resource block, possible jump points occur between the edges of two neighboring resource blocks. When a precoded unit includes multiple resource blocks, possible jump points occur between the edges of two neighboring precoded units. 
         [0032]    Hereinafter is described how to compensate channel coherency corrupted due to use of a precoding matrix, according to different embodiments of the present invention. 
         [0033]    A. Keep a Precoding Channel Coherent 
         [0034]    Solution 1: Phase Rotation 
         [0035]    First, the following description is provided taking a case where the precoded unit includes one sub-carrier for example. A person skilled in the art would understand that a precoded unit may include multiple sub-carriers, such as one or more resource block. In an OFDM system, each resource block includes 12 sub-carriers. 
         [0036]    Firstly, channel status information is obtained. For a TDD system, the base station  1  may estimate the downlink channel according to the received uplink reference signal transmitted by the mobile terminal  2 , to obtain a channel matrix. And in a FDD system, the mobile terminal  2  measures the downlink channel and feeds back the measured downlink channel to the base station  1 , and thus the base station  1  can obtain the downlink channel matrix. 
         [0037]    Then, the base station  1  performs matrix decomposition to the downlink channel matrix, to obtain an initial precoding matrix. In a preferable embodiment, SVD is utilized to obtain an initial precoding matrix V. In a varied embodiment, QR decomposition may be utilized to obtain an initial precoding matrix Q, where Q represents an orthogonal matrix and R represents an upper triangular matrix. Note that the solution of these matrix decompositions is not unique. Therefore, the obtained corresponding precoding matrix is also not unique. 
         [0038]    Then, a phase rotation matrix is used, to keep coherent the channel precoded by the initial precoding matrices corresponding to multiple different precoded units. Therefore, the new precoding matrix is represented by: 
         [0000]      F i =V i G i    (3)
 
         [0000]    where G i  is a diagonal matrix utilized to adjust the phase of the initial precoding channel U i D i  such that channel coherence between multiple precoded units is recovered. 
         [0039]      FIG. 5  is a flowchart illustrating a method according to an embodiment of the present invention, 
         [0040]    In step S 50 , the base station  1  detects phase of an initial precoding channel at each sub-carrier and obtains phase of the initial precoding channel at the sub-carrier neighboring each sub-carrier. For example, the base station  1  hopes to investigate phase of one sub-carrier, which is called a target sub-carrier hereinafter. 
         [0041]    Then in step S 51 , the base station  1  compares the difference between the phase of the initial precoding channel at the target sub-carrier and the phase of an corrected (initial) precoding channel at a neighboring sub-carrier with (−2π, −π, 0, π, 2π). Optionally, the neighboring sub-carrier is the immediately preceding sub-carrier of the target sub-carrier. If the difference is closest to 0 or ±2π, it follows that the phase of the initial precoding channel at the target sub-carrier has not jumped with respect to that of the corrected (initial) precoding channel at the immediately preceding sub-carrier; otherwise, if the difference is closest to ±π, it follows that the phase of the initial precoding channel at the target sub-carrier has jumped with respect to that of the corrected (initial) precoding channel at the immediately preceding sub-carrier. For example, the differences between on one hand h 11  and h 21  in the initial precoding channel taking i as 1 and on the other hand h 11  and h 21  in the initial precoding channel taking i as 0 are closest to π, while the differences between on one hand h 12  and h 22  in the initial precoding channel taking i as 1 and on the other hand h 12  and h 22  in the initial precoding channel taking i as 0 are closest to 2π, then a diagonal matrix G 1  is constructed, 
         [0000]    
       
         
           
             
               G 
               1 
             
             = 
             
               
                 [ 
                 
                   
                     
                       
                         - 
                         1 
                       
                     
                     
                       0 
                     
                   
                   
                     
                       0 
                     
                     
                       1 
                     
                   
                 
                 ] 
               
               . 
             
           
         
       
     
         [0000]    Therefore, the corrected precoding channel H 1p1  weighted by a precoding matrix transformed by the diagonal matrix G 1  makes H 1p1  and H 0p1  continuous, where subscript p1 represents a corrected precoding channel. Then, phase of the initial preceding channel at the second sub-carrier is compared with phase of the corrected initial preceding channel at the first sub-carrier, Generally, H ip1 =H i V i G i =U i D i V i   H V i G i =U i D i G i , G i  makes continuous the phase of the corrected preceding matrix of the i-th preceded unit and that of the corrected preceding matrix of the i-1-th preceded unit. 
         [0042]    When the preceded unit includes multiple sub-carriers, preceding channels corresponding to sub-carriers with the same index within different preceded units are compared, because each sub-carrier in the same preceded unit utilizes the same preceding matrix. For example, in step S 50 ′, the base station  1  detects phase of the initial preceding channel at the first sub-carrier in each preceded unit. Then, according to the phase of the initial preceding channel at the first sub-carrier in the i-th preceded unit and the phase of the corrected (initial) preceding channel at the first sub-carrier in the i-1-th preceded unit, if discontinuous phase is found, in step S 51 ′, phase of each sub-carrier in the i-th preceded unit is adjusted accordingly. 
         [0043]      FIG. 6  illustrates a channel curve with respect to that of  FIG. 4 , after phase rotation. 
         [0044]    Due to complexity of mobile environment, a signal from a transmitter to a receiver usually includes multiple signal components resulting from reflection, diffraction, etc. And different signal components arrive at the receiver with different strength, time, and direction, which vary greatly in different environments. Due to different time-of-arrival of the different multipath components, the received signal is spread in the time domain. The basic feature of multipath propagation is that each path of signals arriving at the receiver has different attenuation factors and delays. The time domain spread of the received signal is called delay spread, which directly reflects the frequency selectivity of a channel (different frequency spectrums of a signal carry different power). Delay spread is defined as the maximum delay among multiple paths. Therefore, the time domain characteristics of the preceding channel with rotated phase, i.e., the corrected preceding channel, is further observed. 
         [0045]    Statistical characteristics of a preceding channel with rotated phase and a non-precoded channel are compared. Taking into account channel coherency in frequency domain, the present disclosure utilizes precoding of each resource block with rotated phase to approximate a more practical scenario. The precoding channel with rotated phase is transformed to the time domain (corresponding to 120-point IFFT allocated for 10 resource blocks). The observed statistical characteristics are shown in  FIGS. 7 and 8 . 
         [0046]    As shown in  FIGS. 7 and 8 , the corrected precoding channel after precoding has very similar statistical characteristics to that of a raw channel without precoding. Therefore, a user equipment can perform channel estimation at the allocated resources to determine parameters of the precoding channel. 
         [0047]    In solution 1, when selecting G i , division is mainly performed according to different sub-carriers in the frequency domain. In the above example, impulse responses of initial precoding channels of one or more sub-carrier in a precoding granularity and one or more sub-carrier in a neighboring precoding granularity are compared to perform correction. Besides frequency domain, time domain coherency may also be considered in phase rotation operation, if necessary. In other words, coherency between one or more symbol (time slot, or subframe) in a precoding granularity and one or more symbol (time slot, or subframe) in a neighboring precoding granularity is considered. Operations performed to keep time domain coherency between different precoding granularities are similar to those performed to keep frequency domain coherency. Obviously, correction for keeping time domain coherency may be performed after the correction for keeping frequency domain coherency, or may be performed separately. 
         [0048]    Solution 2: Amplitude Smoothing 
         [0049]    Firstly, channel status information is obtained. For a TDD system, the base station  1  may estimate the downlink channel according to the received uplink reference signal transmitted by the mobile terminal  2  to obtain a channel matrix. And in a FDD system, the mobile terminal  2  measures the downlink channel and feeds back the measured downlink channel to the base station  1 , and thus the base station  1  can obtain the downlink channel matrix. 
         [0050]    Then, the base station  1  performs matrix decomposition to the downlink channel matrix to obtain an initial precoding matrix. In a preferable embodiment, SVD is utilized to obtain an initial precoding matrix V. In a varied embodiment, QR decomposition may be utilized to obtain an initial precoding matrix Q, where Q represents an orthogonal matrix and R represents an upper triangular matrix. Note that the solution of these matrix decompositions is not unique. Therefore, the obtained corresponding precoding matrix is also not unique. 
         [0051]    Besides phase rotation, the matrix G i  may also he utilized to smooth amplitude of the initial precoding channel. Hereinafter is provided detailed description of an FFT-based smoothing solution with reference to  FIG. 9 . Phase smoothing may be performed in an overall initial precoding channel consisting of initial precoding channels corresponding to multiple precoding granularities. The overall initial precoding channel is denoted as H int precoded . 
         [0052]    Firstly, in step S 90 , the base station  1  performs inverse Fourier transformation (IFFT) to the channel H int precoded  and thus obtains a time domain channel impulse response h. 
         [0053]    Then, in step S 91 , the base station  1  truncates h, keeps a certain length (e.g., a length of acyclic prefix), and sets the points that are truncated off to zeros, thereby obtaining h clip . 
         [0054]    For example, according to an expected maximum multipath delay, the base station  1  may keep the points in h before a time point corresponding to the expected maximum multipath delay and set the remaining points to zeros. This corresponds to reducing maximum multipath delay. And the shorter the maximum multipath delay is, the smoother the amplitude of the frequency domain channel becomes. 
         [0055]    Then, in step S 92 , the base station  1  performs Fourier transformation (FFT) to h clip  to recover a smooth channel H smooth . 
         [0056]    Then, in step S 93 , a correcting matrix G is obtained according to H smooth , which can be expressed by 
         [0000]        G=H   smooth ·/( UD )   (4)
 
         [0000]    where ·/ denotes scalar quantity division (element-wise division), meaning each element in a matrix is abstracted, the value of each element in multiple precoding unit is divided. 
         [0057]    In solution  2 , when selecting G, division is mainly performed according to different sub-carriers in the frequency domain. Besides frequency domain, time domain coherency may also he considered in amplitude smoothing operation if necessary. In other words, coherency between one or more symbol (time slot, or subframe) in a precoding granularity and one or more symbol (time slot, or subframe) in a neighboring precoding granularity is considered. Operations performed to keep time domain coherency between different precoding granularities are similar to those performed to keep frequency domain coherency. Obviously, correction for keeping t me domain coherency may be performed after correction for keeping frequency domain coherency, or be performed separately. 
         [0058]    Solution 3: Combination of Phase Rotation and Amplitude Smoothing 
         [0059]    In a varied embodiment, solutions 1 and 2 may be combined to correct the initial precoding matrix. 
         [0060]    For example, after the operation of phase rotation, the operation of amplitude smoothing may be further performed to the matrix after phase rotation to further improve coherency of the precoding matrix. 
         [0061]    Scenario A is completely transparent to a terminal, and no modification needs to be made to the mobile terminal  2 . 
         [0062]    Since the corrected precoding channel under the condition of the weighting of the corrected precoding matrix satisfies coherency, in the three solutions, the mobile terminal  2  can perform unified channel estimation across different precoding granularities. 
         [0063]    B. Keep the Eigenvalue Matrix of the Precoding Channel Coherent 
         [0064]    The matrix G i  can be utilized to smooth the eigenvalue matrix D i  of the precoding channel, instead of U i D i , such that D i G i  becomes a flat fading diagonal matrix. 
         [0065]    Specifically, linear transformation may be performed to the initial precoding matrix UD, for example, by left multiplying the inverse matrix U −1  of U. Then, subsequent operations are performed with reference to one of solutions 1-3 in scenario A. Note that when solution 2 is utilized to perform amplitude correction, formula (4) should be modified to formula (5): G=H smooth ·/D. 
         [0066]    Since amplitude is adjusted, this solution for smoothing the eigenvalue matrix is actually a power allocation algorithm. 
         [0067]    Since the corrected precoding channel under the condition of the weighting of the corrected precoding matrix satisfies coherency, the mobile terminal  2  can perform unified channel estimation across different precoding granularities. 
         [0068]    Hereinafter is provided a description of the present invention in the apparatus perspective.  FIG. 10  is a block diagram illustrating an apparatus according to an embodiment of the present invention. The processing apparatus  10  in  FIG. 10  is located in the base station  1 . The processing apparatus  10  includes an initial precoding matrix obtaining means  100 , a correcting means  101 , and a transmitting means  102 . Wherein, the correcting means  101  includes a first correcting matrix obtaining means  1010 , a rotation means  1011 , a second correcting matrix obtaining means  1012 , an amplitude adjusting means  1013 , and an eigenvalue matrix obtaining means  1014 . 
         [0069]    The following description is provided taking a case where the precoded unit includes one sub-carrier for example. A person skilled in the art would understand that a precoded unit may include multiple sub-carriers, such as one or more resource block. In an OFDM system, each resource block includes 12 sub-carriers. 
         [0070]    Firstly, the initial precoding matrix obtaining means  100  obtains channel status information. For a TDD system, the initial precoding matrix obtaining means  100  may estimate the downlink channel according to the received uplink reference signal transmitted by the mobile terminal  2  to obtain a channel matrix. And in a FDD system, the mobile terminal  2  measures the downlink channel and feeds back the measured downlink channel to the initial precoding matrix obtaining means  100 , and thus the base station  1  can obtain the downlink channel matrix. 
         [0071]    Then, the initial precoding matrix obtaining means  100  performs matrix decomposition to the downlink channel matrix to obtain an initial precoding matrix. In a preferable embodiment, SVD is utilized to obtain the initial precoding matrix V. In a varied embodiment, QR decomposition may be utilized to obtain an initial precoding matrix Q, where Q represents an orthogonal matrix and R represents an upper triangular matrix. Note that the solution of these matrix decompositions is not unique. Therefore, the obtained corresponding precoding matrix is also not unique. 
         [0072]    Then, a correcting matrix is utilized to keep coherence of the channels precoded by the initial preceding matrix corresponding to multiple different precoded units. Therefore, the new precoding matrix is represented by: 
         [0000]      F i =V i G i    (3)
 
         [0000]    where G i  may be a diagonal matrix utilized to adjust the phase of the initial preceding channel U i D i  such that channel coherence between multiple precoded units is recovered. 
         [0073]    Solution I: Phase Rotation 
         [0074]    Firstly, the first correcting matrix obtaining means  1010  in the correcting means  101  detects the phase of an initial precoding channel at each sub-carrier and obtains the phase of an initial precoding channel at a sub-carrier neighboring each sub-carrier. For example, the base station  1  hopes to investigate the phase of one sub-carrier, which is called as a target sub-carrier hereinafter. 
         [0075]    Then, the difference between the phase of the initial precoding channel at the target sub-carrier and that of an corrected (initial) precoding channel at a neighboring sub-carrier is compared with (−2π, −π, 0, π, 2π). Optionally, the neighboring sub-carrier is the immediately preceding sub-carrier of the target sub-carrier. If the difference is closest to 0 or ±2π, it follows that the phase of the initial precoding channel at the target sub-carrier has not jumped with respect to that of the corrected (initial) precoding channel at the immediately preceding sub-carrier; otherwise, if the difference is closest to ±π, it follows that the phase of the initial precoding channel at the target sub-carrier has jumped with respect to that of the corrected (initial) precoding channel at the immediately preceding sub-carrier. For example, the differences between on one hand h 11  and h 21  in the initial precoding channel taking i as 1 and on the other hand h 11  and h 21  in the initial precoding channel taking i as 0 are closest to π, while the differences between on one hand h 12  and h 22  in the initial precoding channel taking i as 1 and on the other hand h 12  and h 22  in the initial precoding channel taking i as 0 are closest to 2π, then a diagonal matrix G 1  is constructed, 
         [0000]    
       
         
           
             
               G 
               1 
             
             = 
             
               
                 [ 
                 
                   
                     
                       
                         - 
                         1 
                       
                     
                     
                       0 
                     
                   
                   
                     
                       0 
                     
                     
                       1 
                     
                   
                 
                 ] 
               
               . 
             
           
         
       
     
         [0000]    Therefore, the corrected precoding channel H 1p1  weighted by the precoding matrix transformed by the phase rotation means  1011  via the diagonal matrix G 1  makes H 1p1  and H 0p1  continuous, where subscript p1 represents a corrected precoding channel. Then, phase of the initial precoding channel at the second sub-carrier is compared with phase of the corrected initial precoding channel at the first sub-carrier. Generally, H ip1 =H i V i G i =U i D i V i   H V i G i =U i D i G i , G i  makes continuous the phase of the corrected precoding matrix of the i-th precoded unit and that of he corrected precoding matrix of the i-1-th precoded unit. 
         [0076]    Then, the transmitting means  102  transmits to the mobile terminal  2  a pilot and/or data weighted by the precoding matrix corrected by the first correcting means. 
         [0077]    When the precoded unit includes multiple sub-carriers, precoding channels corresponding to sub-carriers with the same index from different precoded units are compared because each sub-carrier in the same precoded unit utilizes the same precoding matrix. For example, the first correcting matrix obtaining means  1010  detects phase of the initial precoding channel at the first sub-carrier in each precoded unit. Then, according to the phase of the initial precoding channel at the first sub-carrier in the i-th precoded unit and the phase of the corrected (initial) precoding channel at the first sub-carrier in the i-1-th precoded unit, if discontinuous phase is found, the rotation means  1011  adjusts the phase of each sub-carrier in the i-th precoded unit accordingly. 
         [0078]      FIG. 6  illustrates a channel curve with respect to that of  FIG. 4 , after phase rotation. 
         [0079]    Due to complexity of mobile environment, a signal from a transmitter to a receiver usually includes multiple signal components resulting from reflection, diffraction, etc. And different signal components arrive at the receiver with different strength, time, and direction, which vary greatly in different environments. Due to different time-of-arrival of different multipath components, the received signal is spread in the time domain. The basic feature of multipath propagation is that each path of signals arriving at the receiver has different attenuation factors and delays. The time domain spread of the received signal is called as delay spread, which directly reflects frequency selectivity of a channel (different frequency spectrum of a signal carry different power). Delay spread is defined as the maximum delay among multiple paths. Therefore, the time domain characteristics of a precoding channel with rotated phase, i.e., a corrected precoding channel, is further observed. 
         [0080]    Statistical characteristics of a precoding channel with rotated phase and a non-precoded channel are compared. Taking into account frequency domain channel coherency, the present disclosure utilizes precoding of each resource block with rotated phase to approximate a more practical scenario. The precoding channel with rotated phase is transformed to the time domain (120-point IFFT allocated for 10 resource blocks). The observed statistical characteristics are shown in  FIGS. 7 and 8 . 
         [0081]    As shown in  FIGS. 7 and 8 , the corrected precoding channel after precoding has very similar statistical characteristics to that of a raw channel without precoding. Therefore, a user equipment can perform channel estimation at the allocated resources to determine parameters of the precoding channel. 
         [0082]    In solution I, when selecting G i , division is mainly performed according to different sub-carriers in the frequency domain. In the above example, impulse responses of initial precoding channels of one or more sub-carrier in a precoding granularity and one or more sub-carrier in a neighboring precoding granularity are compared to perform correction. Besides frequency domain, time domain coherency may also be considered in phase rotation operation if necessary. In other words, coherency between one or more symbol (time slot, or subframe) in a precoding granularity and one or more symbol (time slot, or subframe) in a neighboring precoding granularity is considered. Operations performed to keep time domain coherency between different precoding granularities are similar to those performed to keep frequency domain coherency. Obviously, correction for keeping time domain coherency may be performed after correction for keeping frequency domain coherency, or be performed separately. 
         [0083]    Solution II: Amplitude Smoothing 
         [0084]    Firstly, channel status information is obtained. For a TDD system, a base station  1  may estimate the downlink channel according to the received uplink reference signal transmitted by a mobile terminal  2  to obtain a channel matrix. And in a FDD system, the mobile terminal  2  measures the downlink channel and feeds back the measured downlink channel to the base station  1 , and thus the base station  1  can obtain the downlink channel matrix. 
         [0085]    Then, the base station performs matrix decomposition to the downlink channel matrix o obtain an initial precoding matrix. In a preferable embodiment, SVD is utilized to obtain an initial precoding matrix V. In a modified embodiment, QR decomposition may be utilized to obtain an initial precoding matrix Q, where Q represents an orthogonal matrix and R represents an upper triangular matrix. Note that the solution of these matrix decompositions s not unique. Therefore, the obtained corresponding precoding matrix is also not unique. 
         [0086]    Besides phase rotation, the matrix G i  may also be utilized to smooth amplitude of an initial precoding channel. Hereinafter is provided detailed description of an FFT-based smoothing solution with reference to  FIG. 9 . Phase smoothing may be performed in an overall initial precoding channel consisting of initial precoding channels corresponding to multiple precoding granularties. The overall initial precoding channel is denoted as H int precoded . 
         [0087]    Firstly, an inverse Fourier transforming means (not shown) in the second correcting matrix obtaining means  1012  performs inverse Fourier transforming (IFFT) to the channel H int precoded  and thus obtains a time domain channel impulse response h. 
         [0088]    Then, a truncating means (not shown) in the second correcting matrix obtaining means  1012  truncates h, keeps a certain length (e.g., a length of a cyclic prefix), and sets the points that are truncated off to zeros, thereby obtaining h clip . 
         [0089]    For example, according to an expected maximum multipath delay, the truncating means may keep the points in h before a time point corresponding to the expected maximum multipath delay and set the remaining points to zeros. This corresponds to reducing maximum multipath delay. And the shorter the maximum multipath delay is, the smoother the amplitude of the frequency domain channel becomes. 
         [0090]    Then, an Fourier transforming means (not shown) in the second correcting matrix obtaining means  1012  performs Fourier transforming (FFT) to h clip  to recover a smooth channel H smooth . 
         [0091]    Then, the second correcting matrix obtaining means  1012  obtains a correcting matrix G according to H smooth , which can be expressed by 
         [0000]      G=H smooth ·/( UD )   (4)
 
         [0000]    where ·/ denotes scalar quantity division (element-wise division), meaning each element in a matrix is abstracted, the value of each element in multiple precoding unit is divided. 
         [0092]    Then, the amplitude smoothing means  1013  smoothes the amplitude of the initial precoding channel according to the precoding matrix corrected by the second correcting matrix G obtained by the second correcting means. 
         [0093]    Then, the transmitting means  101  weights the data and pilot with the corrected precoding matrix corrected by the second correcting matrix and transmits them to the mobile terminal  2 . 
         [0094]    In solution when selecting G, division is mainly performed according to different sub-carriers in the frequency domain. Besides frequency domain, time domain coherency may also be considered in amplitude smoothing operation if necessary. In other words, coherency between one or more symbol (time slot, or subframe) in a precoding granularity and one or more symbol (time slot, or subframe) in a neighboring precoding granularity is considered. Operations performed to keep time domain coherency between different precoding granularities are similar to those performed to keep frequency domain coherency. Obviously, correction for keeping time domain coherency may be performed after correction for keeping frequency, domain coherency, or be performed separately. 
         [0095]    Solution III: Combination of Phase Rotation and Amplitude Smoothing 
         [0096]    In a varied embodiment, solutions I and II may be combined to correct the initial precoding matrix. 
         [0097]    For example, after the operation of phase rotation, the operation of amplitude smoothing may be further performed to the matrix after phase rotation to further improve coherency of the precoding matrix. 
         [0098]    The above three solutions are completely transparent to a terminal and no modification needs to be made to the mobile terminal  2 . 
         [0099]    Since the corrected precoding channel under the condition of weighting of the corrected precoding matrix satisfies coherency, in the three solutions, the mobile terminal  2  can perform unified channel estimation across different precoding granularities. 
         [0100]    C. Keep Coherent the Eigenvalue Matrix of the Precoding Channel 
         [0101]    The matrix G i  can be utilized to smooth the eigenvalue matrix D i  of the precoding channel, instead of U i D i , such that D i G i  becomes a flat fading diagonal matrix. 
         [0102]    Specifically, referring to  FIG. 11 , the eigenvalue matrix obtaining means  1014  performs linear transformation to the initial precoding matrix UD, for example, by left multiplying the inverse matrix U −1  of U. Then, subsequent operations are performed with reference to one of solutions 1-3 in scenario A. Note that when solution 2 is utilized to perform amplitude correction, formula (4) should be modified to formula (5): G=H smooth ·/D. 
         [0103]    Since amplitude is adjusted, this solution for smoothing the eigenvalue matrix actually a power allocation algorithm. 
         [0104]    Since the corrected precoding channel under the condition of weighting of the corrected precoding matrix satisfies coherency, the mobile terminal  2  can perform unified channel estimation across different precoding granularities. 
         [0105]    Although the present invention is illustrated and described in the drawings and the above description, such illustration and description is explanatory and exemplary instead of limiting. Therefore, the present invention is not limited to the above embodiments. 
         [0106]    Other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims. In the claims, the word “comprising” does not exclude other elements or steps, and the indefinite article “a” or “an” does not exclude a plurality. In practical applications of the present invention, a single unit may fulfill the functions of several items recited in the claims. Any reference signs in the claims should not be construed as limiting the scope.