# Patent

## Title
Reciprocity calibration for multiple-input multiple-output systems

## Abstract

Systems and associated methods for reciprocity calibration of multiple-input multiple-output (MIMO) wireless communication are disclosed herein. In one embodiment, a method for reciprocity calibration of the MIMO system includes transmitting a pilot symbol by a transmitter (TX) of the reference antenna and receiving the pilot symbol by receivers (RXes) of antennas of a base station as ri,0 pilot symbols. (Index “i” denotes individual antenna “i” of the base station, and “0” denotes the reference antenna.) The method further includes transmitting the received pilot symbols by TXes of the antennas of the base station, receiving the pilot symbols transmitted by the antennas of the base station by the reference antenna as r0,i pilot symbols, and calculating non-reciprocity compensation factors as




  
    
      r
      
        i
        ,
        0
      
    
    
      r
      
        0
        ,
        i
      
    
  
  .

## Background
TECHNICAL FIELD
The disclosed embodiments are generally related to the field of multiple-input multiple-output (MIMO) wireless communication between base stations and mobile device users. More specifically, the disclosed embodiments are related to calibration of the MIMO systems.
BACKGROUND
As spectral bandwidth becomes an ever more valuable commodity for radio communication systems, techniques are needed to effectively use the available bandwidth. MIMO wireless technology is one such technique. MIMO uses multiple transmitters and multiple receivers to increase data throughput, e.g., to increase data transfer bandwidth, and/or to improve the signal-to-noise (SNR) ratio, e.g., to reduce the signal fading. With respect to increased data transfer bandwidth, the MIMO increases the capacity of a data channel by increasing the number of receive (RX) and transmit (TX) antennas to take advantage of signal multipathing. For example, multiple pairs of transmit/receive antennas can be deployed in the communication system to transmit multiple spatial streams (e.g., wireless signals or symbols). These data streams can be transmitted simultaneously and on a single frequency by using proper multiplexing (MUX) at the TX side of the communication system. The combination of multiple, independent data paths increases the throughput of a MIMO system. At the RX side of the MIMO system, the transferred data are recombined into the original data stream using suitable demultiplexing (deMUX). Typically, an increase of the number of RX and TX antenna pairs linearly increases the throughput of the MIMO system.
Furthermore, the MIMO systems can improve SNR ratio of the wireless data. For example, the same wireless signals can take multiple paths from the TX to the RX as the wireless signals reflect from objects that are in the path of, or adjacent to, the TX and RX. In absence of MIMO, these multiple paths would generally increase signal interference at the RX. However, with MIMO systems these additional paths can be utilized to improve the SNR of the wireless signals. For example, in a relatively simple implementation, a MIMO system may identify and select the strongest signal among the multiple available signals, while discarding other signals. In more advanced implementations, a MIMO system may accept multiple wireless signals as representing different versions of the same wireless signal having slight delays because of the reflection of the signal or because of differing distances traveled by the various paths from the TX antenna to the RX antenna. Generally, multiple versions of the same signal undergo statistically independent fading because the signals travel along different physical paths. When these different versions of the signal are summed up (or otherwise combined using proper adjustments for differing time delays or phase changes) into a resulting signal, the SNR improves due to statistical independence of the signal fading.
However, for practical MIMO systems, it is difficult to properly characterize different data channels (e.g., to calibrate gains and delays from one antenna to another). Furthermore, data channels may behave differently based on the direction of data transfer (e.g., whether the antenna sends the signal or receives it, or whether the wireless signal is transferred from antenna 1 to antenna 2 or from antenna 2 to antenna 1). Accordingly, a need remains for improved characterization of data channels in MIMO systems.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic view of a MIMO system in accordance with an embodiment of the present technology.
FIG. 2 is a graph diagram of upload (DL) antenna calibration in accordance with an embodiment of the present technology.
FIG. 3 is a graph diagram of upload (UL) antenna calibration in accordance with an embodiment of the present technology.



DETAILED DESCRIPTION
The following disclosure describes various embodiments of systems and associated methods for calibrating data channels in multiple-input multiple-output (MIMO) communication systems. In some embodiments, MIMO communication systems can be used for both upload (UL) and download (DL) of the wireless data (also referred to as symbols). With these systems, an individual antenna of a MIMO transceiver (TRX) at a base station can be connected to a receiver (RX) when the antenna receives data from user equipment (UE), e.g., a mobile phone, or to a transmitter (TX) when the antenna sends data to the UE. However, the RX and TX channels of the same TRX are not necessarily reciprocal, e.g., the channels may have different delays, phase distortions, amplitude amplification, etc. Therefore, in at least some embodiments, dedicated calibration symbols (also referred to as pilot symbols or calibration tones) can be transmitted from one of the antennas of the base station to other antennas to calibrate the MIMO communication system by estimating the behavior of the RX/TX propagation channels. In some embodiments, a reciprocity compensation matrix (CT) can be calculated by transmitting and receiving pilot symbols and estimating their amplitude and phase. The application of a reciprocity compensation matrix can improve reciprocity of the data channels, thereby, e.g., improving system throughput or signal to noise ratio (SNR).
In some embodiments, the pilot symbols can be transmitted by an auxiliary (AUX) antenna that is specifically dedicated to the calibration of the MIMO system. In some embodiments, the pilot symbols can be transmitted outside of the frequency bandwidth that is dedicated to normal data payload. With these embodiments, the pilot symbols may be transmitted either continuously or intermittently. In some other embodiments, the pilot symbols may be transmitted within the frequency bandwidth that is reserved for the data payload by, for example, intermittently pausing the transmission of the payload data and transmitting the pilot symbols. In at least some embodiments, the power (e.g., the power spectral density) of the pilot symbols may be comparable to the power of the symbols of the payload data to minimize phase noise and compression at the RXes.
FIG. 1 is a schematic view of a MIMO system 1000 in accordance with a disclosed embodiment. The illustrated MIMO system includes up to NB TRXes 105/106 at a base station 100B, and up to NM TRXes 155/156 at a user equipment (UE) side 100M. At the base station 100B, multiple TXes 105b-105c may transmit signals from their respective antennas 120a-120c to one or more antennas 170a-170b connected to RXes 155a-155b of the UE side 100M. In some embodiments, beamforming may be used to direct signals from the antennas 120a-120c over multiple paths 150 toward the UE side 100M. The signals received by the RXes 155a-155b of the UE side can be further processed and sent to user equipment as, e.g., audio signal in phone conversation. Conversely, TXes 156a-156b can send signals from the antennas 170a-170b at the UE side 100M to the antennas 120a-120c and their corresponding RXes 106a-106c of the base station 100B. The received signals at the base station may be audio signals, digital images, text signals, etc., that are retransmitted to, e.g., a different antenna (e.g., a different user) at the UE side 100M. In some embodiments, the antennas 170a-170b at the UE side 100M can be used for beamforming to, for example, improve signal transmission toward a selected antenna of the base station 100B. As an example, the TRX 105a/106a includes a TX 105a and an RX 106a that can be individually connected to an antenna 120a through a switch 110a. The illustrated TX 105a may send symbol S0 having amplitude fT0 and phase ejΦT0 to one or more antennas 170 at the UE side 100M of the MIMO system 1000. The symbol S0 may be received as, for example, symbol rB,NM−1P at an RX 155b at the UE side 100M. In another example, a TX 156a may send a pilot symbol SP,0 having amplitude gT0 and phase ejφR0 (or other, non-pilot symbol S0) to, for example, RX 106c that receives the sent symbol as a symbol rB,NB−1P having amplitude




  f
  
    R
    
      
        N
        B
      
      -
      1
    
  




and phase ejΦRNB−1. Many other combinations of transmitted and received symbols are possible with the MIMO system 1000. Two TRXes 160a, 160b are illustrated at the UE 100M, but the UE may include a larger number of the TRXes. In the illustrated embodiment, the TXes 105a and 105c of the base station 100B and the TX 156b of the UE 100M send symbols to their respective antennas, while the RX 106b of the base station 100B and the RX 155a of the UE 100M receive symbols, but other combinations of TX/RX are also possible, depending on particular data transfer or calibration needs. In at least some embodiments, the multiple antennas of the MIMO system 100 may increase data throughput (e.g., data bandwidth) and/or improve the SNR of the transmitted symbols. For example, the symbols S0 transmitted by an antenna 120a at the base station 100B can be received by antennas 170a-170b at the UE side, and may be recombined to improve either data throughput of the system or the SNR or both. In at least some embodiments, the antennas of the base station side 100B may be used for beamforming, i.e., for directing the transmit/receive symbols in a preferred direction.

For a MIMO system that operates in a time domain division (TDD) multiplexing mode, a combined radio propagation channel H (e.g., a transfer function of the channel) that operates between NB antennas at the base station 100B and NM antennas at the UE 100M can be assumed as reciprocal, i.e., it can be assumed that:

HUL=HDLT  (1)

where UL stands for data upload from the base station 100B to the UE 100M, and DL stands for data download from the UE 100M to the base station 100B. Therefore, the MIMO systems described by equation (1) rely on channel reciprocity by, for example, assuming that the TX 105b and RX 106b (or other TX/RX pairs of a transceiver) are reciprocal. However, in some practical applications, the reciprocity between the TX and RX that use the same antenna pair (e.g., are connected to a same antenna) is not present, therefore requiring calibration for improved performance of the overall MIMO system 1000. A non-reciprocal combined radio propagation channel H can be characterized as:

UL≠DLT  (2)

where UL and DL stand for observed (e.g., measured) values of data upload and download, respectively.

In some embodiments, the improved reciprocity may be more important for the base station 100B of the MIMO system 1000 than for the UE side 100M. Therefore, in at least some embodiments, the base station 100B may be calibrated as follows. The received signal at the base station 100B (e.g., the signal received at the RXes 106a-106c) can be represented in a matrix form as:

rBP=ULSP  (3.1)

or, when UL is expressed as a function of HUL, equation (3.1) can be rewritten as:

rBP=BRHULMTSP  (3.2)

where





  
    B
    R
  
  =
  
    diag
    ⁢
    
      {
      
        
          
            f
            
              R
              0
            
          
          ⁢
          
            e
            
              jφ
              
                R
                0
              
            
          
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            f
            
              R
              
                
                  N
                  B
                
                -
                1
              
            
          
          ⁢
          
            e
            
              jφ
              
                R
                
                  
                    N
                    B
                  
                  -
                  1
                
              
            
          
        
      
      }
    
  




(a diagonal matrix of amplitude/phase for the RXes of the base station 100B),





  
    M
    T
  
  =
  
    diag
    ⁢
    
      {
      
        
          
            g
            
              T
              0
            
          
          ⁢
          
            e
            
              jφ
              
                T
                0
              
            
          
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            g
            
              T
              
                
                  N
                  M
                
                -
                1
              
            
          
          ⁢
          
            e
            
              jφ
              
                T
                
                  
                    N
                    M
                  
                  -
                  1
                
              
            
          
        
      
      }
    
  




(a diagonal matrix of amplitude/phase for the TXes of the UE 100M), and SP is a NM×1 vector of transmitted symbols (e.g., pilots) during the upload from the base station to the UE.

Analogously to equations (3.1) and (3.2) above, the received signals at the UE 100M can be represented as:

rMP=DLSP  (4.1)

or as:

rMP=MRHULTBTSP  (4.2)

where





  
    M
    R
  
  =
  
    diag
    ⁢
    
      {
      
        
          
            g
            
              R
              0
            
          
          ⁢
          
            e
            
              jφ
              
                R
                0
              
            
          
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            g
            
              R
              
                
                  N
                  M
                
                -
                1
              
            
          
          ⁢
          
            e
            
              jφ
              
                R
                
                  
                    N
                    M
                  
                  -
                  1
                
              
            
          
        
      
      }
    
  




(a diagonal matrix of amplitude/phase for the RXes of the UE 100M), and





  
    B
    T
  
  =
  
    diag
    ⁢
    
      {
      
        
          
            f
            
              T
              0
            
          
          ⁢
          
            e
            
              jφ
              
                T
                0
              
            
          
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            f
            
              T
              
                
                  N
                  B
                
                -
                1
              
            
          
          ⁢
          
            e
            
              jφ
              
                T
                
                  
                    N
                    B
                  
                  -
                  1
                
              
            
          
        
      
      }
    
  




(a diagonal matrix of amplitude/phase for the TXes of the base station 100B), and SP is a NB×1 vector of transmitted symbols (e.g., pilots) during the upload from the base station to the UE.

As discussed with respect to equation (2), the estimated (e.g., measured) values of UL and DLT are not necessarily the same. However, since we know from equation (1) that the propagation channel can be modeled as HUL=HDLT, the equations (3.2) and (4.2) can be recombined as:

BR−1ULMT−1=BT−TDLTMR−T  (5)

and can be further reduced to:

DL=MRMT−TULTB−TRBT  (6)

Since by definition M=MRMT−T and B=BR−TBT, equation (6) can be further simplified to:

DL=MULTB  (7)

The diagonal matrices M (for the UEs) and B (for the base station) can be further rewritten as:





  
    
      
        M
        =
        
          diag
          ⁢
          
            {
            
              
                
                  
                    g
                    
                      R
                      0
                    
                  
                  
                    g
                    
                      T
                      0
                    
                  
                
                ⁢
                
                  e
                  
                    j
                    ⁡
                    
                      (
                      
                        
                          φ
                          
                            R
                            0
                          
                        
                        -
                        
                          φ
                          
                            T
                            0
                          
                        
                      
                      )
                    
                  
                
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  
                    g
                    
                      R
                      
                        R
                        
                          
                            N
                            M
                          
                          -
                          1
                        
                      
                    
                  
                  
                    g
                    
                      T
                      
                        R
                        
                          
                            N
                            M
                          
                          -
                          1
                        
                      
                    
                  
                
                ⁢
                
                  e
                  
                    j
                    ⁡
                    
                      (
                      
                        
                          φ
                          
                            R
                            
                              
                                N
                                M
                              
                              -
                              1
                            
                          
                        
                        -
                        
                          φ
                          
                            T
                            
                              
                                N
                                M
                              
                              -
                              1
                            
                          
                        
                      
                      )
                    
                  
                
              
            
            }
          
        
      
    
    
      
        (
        8.1
        )
      
    
  




which can be simplified to:





  
    
      
        M
        =
        
          diag
          ⁢
          
            {
            
              
                
                  G
                  0
                
                ⁢
                
                  e
                  
                    jΔφ
                    0
                  
                
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  G
                  
                    
                      N
                      M
                    
                    -
                    1
                  
                
                ⁢
                
                  e
                  
                    jΔφ
                    
                      
                        N
                        M
                      
                      -
                      1
                    
                  
                
              
            
            }
          
        
      
    
    
      
        (
        8.2
        )
      
    
  




Similarly, the matrix B can be rewritten and simplified as:





  
    
      
        B
        =
        
          diag
          ⁢
          
            {
            
              
                
                  F
                  0
                
                ⁢
                
                  e
                  
                    jΔφ
                    0
                  
                
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  F
                  
                    
                      N
                      B
                    
                    -
                    1
                  
                
                ⁢
                
                  e
                  
                    jΔφ
                    
                      
                        N
                        B
                      
                      -
                      1
                    
                  
                
              
            
            }
          
        
      
    
    
      
        (
        9
        )
      
    
  




where





  
    
      F
      i
    
    =
    
      
        f
        
          T
          i
        
      
      
        f
        
          R
          i
        
      
    
  
  ,








  
    
      e
      
        jΔφ
        i
      
    
    =
    
      e
      
        j
        ⁡
        
          (
          
            
              φ
              
                T
                i
              
            
            -
            
              φ
              
                R
                i
              
            
          
          )
        
      
    
  
  ,




i=0, . . . , NB−1. Equation (7) provides a relationship between the measured upload (UL) and download (DL) propagation channels.

In at least some embodiments, the orthogonality among the antennas at the UE side 100M is preserved despite lack of reciprocity of a set of single antennas (SSs). For example, differences in phase and gain between TX and RX paths can still exist without causing a loss in orthogonality provided that these differences in phase/gain are uniform across the antennas of the base station 100B. In at least some embodiments, these uniform differences across the antennas of the base station as well as the lack of reciprocity within an individual TRX can be accounted for by the channel calibration at the UE side 100M. Such a channel calibration is generally known as a relative channel calibration. Some embodiments of the relative channel calibration are described below with reference to FIGS. 2 and 3.
FIG. 2 is a graph diagram 2000 of antenna calibration in accordance with some embodiments. In the illustrated embodiment, the horizontal axis represents frequency and the vertical axis represents amplitude of data symbols. A total bandwidth of the illustrated DL data transmission is denoted as 210. Within the total bandwidth 210, a range of frequencies 220 is a frequency bandwidth of the data payload (e.g., the symbols of the DL data frame) for a sample antenna K of the base station (e.g., corresponding to antenna 120b in FIG. 1). Frequency bandwidths 230 and 240 may correspond to frequency guard bands for the antenna K. In some embodiments, in addition to transmitting the symbols within the frequency bandwidth 220, the antenna K may also transmit one or more calibration symbols (also called pilot symbols or calibration tones). An example of a pilot symbol is shown as a narrow band tone K inside frequency bandwidth 240. Other pilot symbols may be transmitted by other antennas, for example pilot symbols 1, 2 or N may be transmitted by their corresponding antennas 1, 2 or N. In some embodiments, a dedicated auxiliary (AUX) antenna can be used for the pilot symbols thereby preserving the bandwidth of the antennas of the base station for the data payload. In some embodiments, other antennas on the base station side may have comparable bandwidths. The pilot symbols (e.g., the pilot symbol K) can be received by other antennas of the base station (e.g., antennas 120a, 120c). In some embodiments, the power of the pilot symbol(s) may be comparable to the power of the symbols within the DL data frame to reduce phase noise and compression/desense of the RXes. The pilot symbols can be used to estimate the non-reciprocity compensation CT, as explained in more detail below.
FIG. 3 is a graph diagram 3000 of upload (UL) antenna calibration in accordance with an embodiment of the present technology. In the illustrated embodiment, the symbols of the UL data frame are transmitted by UE antennas (e.g., antennas 170 or 170b) within a frequency bandwidth 310. A frequency bandwidth 320 corresponds to the data frame (e.g., symbols in the data payload). In some embodiments, an auxiliary (AUX) antenna (e.g., an antenna dedicated to calibration) may generate symbols illustrated in the graph 3000. Frequency bandwidths 330 and 340 may be reserved for the pilot symbols. For example, the AUX antenna may transmit dedicated pilot symbols 1, 2, . . . , K, . . . N to the corresponding antennas 1-N at the base station. In other embodiments, the AUX antenna may transmit only one or several pilot symbols for one or several antennas of the base station. In some embodiments, an antenna of the UE side may play a role of the AUX antenna for calibration purposes. In some embodiments, the power of the pilot symbol(s) may be comparable to the power of the symbols within the DL data frame to keep different antennas of the base station balanced and/or to minimize the phase noise and compression at the RXes of the base station. In some embodiments, the symbols within the data frame frequency bandwidth (e.g., frequency bandwidths 220 in FIG. 2 and/or 320 in FIG. 3) may be used as the pilot symbols (in every data frame or intermittently in some data frames).
In at least some embodiments, a compensation can be employed to minimize or eliminate the non-reciprocity between the RX and TX path of the TRXes of, for example, base station 100B. For, example a beamforming W in the matrix form NB×NM can be applied to a set of single antennas (SSs). When the non-reciprocity compensation CT (having the matrix form NB×NB) is introduced, the received signal at the UE side 100M can be represented as:

rM=DLCTWS  (10)

The non-reciprocity compensation CT can be subjected to a constraint:

BCT=αejδI  (11)

where I is an identity matrix, α is an amplitude, and ejδ is a phase of the wireless signal. In at least some embodiments, the effective channel observed on pilot symbols on the UE side 100M can be treated as being indeed reciprocal, allowing the received signal rM from equation (10) to be expressed as:

rM=MULTBCTWS  (12.1)

When equation (11) is combined with equation (12.1), the following expression for rM is obtained:

rM=αejδMULTWS  (12.2)

At this point, it becomes possible to compensate for lack of reciprocity of the single antennas (SSs) based on auxiliary frequency non-selective complex gain, which can be based by channel estimation performed on DL symbol transfer. In some embodiments, the signal can be sent from the base station 100B to the UE side 100M using maximal ratio transmission (MRT) beamforming or zero-force (ZF) beamforming. For example, for the UE side impairment matrix WMRT using MRT beamforming, we can write:

WMRT=UL*  (13.1)

or

WMRT=(ULT)H  (13.2)

and

rM=αejδMULTWS  (13.3)

Similarly, for ZF beamforming, the following estimates are valid for the impairment matrix WZF and the received signal rM at the UE side:

WZF=(ULT)T((ULT)(ULT)T)−1  (14.1)

or

WZF=UL(ULTUL)−1  (14.2)

or

rM=ejΔMULTUL(ULTUL)−1S  (14.3)

Equation (14.3) can be reduced to:

rM=αejδMS  (14.4)

As shown in equations (13.3) and (14.3), in either the MRT or ZF beamforming (or other modes of wireless data transfer), the beamforming antennas may remain orthogonal even in absence of reciprocity within the TRXes of the individual antennas.
Without a loss of generality, a calibration procedure can be developed as follows. We can assume:




  
    
      
        
          δ
          =
          
            Δ
            ⁢
            
              
            
            ⁢
            
              φ
              0
            
          
        
        ,
        
          α
          =
          
            
              F
              0
            
            =
            
              
                f
                
                  T
                  0
                
              
              
                f
                
                  R
                  0
                
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




Therefore, the NB×NB non-reciprocity compensation matrix CT (for the base station 100B) can be written as:

CT=αejΦ0B−1  (16.1)

or

CT=diag{1,CT1, . . . , CTBB−1}  (16.2)

The individual components of the CT can be written as:





  
    
      
        
          C
          T
          i
        
        =
        
          
            
              f
              
                T
                0
              
            
            
              f
              
                R
                0
              
            
          
          ⁢
          
            
              f
              
                R
                i
              
            
            
              f
              
                T
                i
              
            
          
          ⁢
          
            e
            
              j
              ⁡
              
                (
                
                  
                    φ
                    
                      T
                      0
                    
                  
                  -
                  
                    φ
                    
                      R
                      0
                    
                  
                
                )
              
            
          
          ⁢
          
            e
            
              -
              
                j
                ⁡
                
                  (
                  
                    
                      φ
                      
                        T
                        i
                      
                    
                    -
                    
                      φ
                      
                        R
                        i
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




Therefore, the individual component CTi represents a ratio with a numerator being a signal sent from an antenna “0” of the base station (e.g., the antenna 105a) and received by antenna “i” of the base station (e.g., the antenna 120b) and a denominator being a signal sent from the antenna “i” and received by the antenna “0.” Some possible paths for the signals sent from the TXes to the RXes of the base station are shown as paths 14b in FIG. 1.

In at least some embodiments, the components CTi of the non-reciprocity compensation matrix CT, can be used for the calibration procedure that uses pilot symbols. As explained with reference to FIGS. 2 and 3, in some embodiments the calibration procedure may involve transmitting pilot symbols from one of the NB antennas (e.g., antenna 0) of the base station 100B and receiving the signal ri,0 on the remaining antennas (e.g., antennas 1 to NB−1) of the base station. Therefore, in this scenario, the antenna 0 may at least temporarily be the reference antenna that transmits signal:




  
    
      
        
          r
          
            0
            ,
            i
          
        
        =
        
          
            f
            
              R
              0
            
          
          ⁢
          
            e
            
              jφ
              
                R
                0
              
            
          
          ⁢
          
            H
            
              0
              ,
              i
            
          
          ⁢
          
            f
            
              T
              i
            
          
          ⁢
          
            e
            
              jφ
              ⁢
              
                
              
              ⁢
              
                T
                i
              
            
          
          ⁢
          
            S
            C
          
        
      
    
    
      
        (
        18
        )
      
    
  




Furthermore, the same set of pilot symbols that was received by the receiving antennas 1 to NB−1 may be transmitted the receiving antennas and received by the reference antenna 0 as:





  
    
      
        
          r
          
            0
            ,
            i
          
        
        =
        
          
            f
            
              R
              0
            
          
          ⁢
          
            e
            
              jφ
              
                R
                0
              
            
          
          ⁢
          
            H
            
              0
              ,
              i
            
          
          ⁢
          
            f
            
              T
              i
            
          
          ⁢
          
            e
            
              jφ
              ⁢
              
                
              
              ⁢
              
                T
                i
              
            
          
          ⁢
          
            S
            C
          
        
      
    
    
      
        (
        19
        )
      
    
  




In at least some embodiments, the transmission from different antennas can be orthogonalized, either in time domain or in frequency domain. Because of the reciprocity of the radio propagation channel (e.g., as in equation (1) above), the followings is true for the antenna 0:

H0,i=Hi,0  (20)

Therefore, for an individual TRX “i,” the non-reciprocity compensation can be expressed as:





  
    
      
        
          
            r
            
              i
              ,
              0
            
          
          
            r
            
              0
              ,
              i
            
          
        
        =
        
          
            
              f
              
                R
                i
              
            
            ⁢
            
              e
              
                jφ
                
                  R
                  i
                
              
            
            ⁢
            
              H
              
                i
                ,
                0
              
            
            ⁢
            
              f
              
                T
                0
              
            
            ⁢
            
              e
              
                jφ
                
                  T
                  0
                
              
            
            ⁢
            
              S
              C
            
          
          
            
              f
              
                R
                0
              
            
            ⁢
            
              e
              
                jφ
                
                  R
                  0
                
              
            
            ⁢
            
              H
              
                0
                ,
                i
              
            
            ⁢
            
              f
              
                T
                i
              
            
            ⁢
            
              e
              
                jφ
                
                  T
                  i
                
              
            
            ⁢
            
              S
              C
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




When Hi,0 and H0,i are eliminated from equation 21 based on equation 20, and when ri,0/r0,i are replaced with CTi, equation 21 can be rewritten as:





  
    
      
        
          C
          T
          i
        
        =
        
          
            
              f
              
                T
                0
              
            
            
              f
              
                R
                0
              
            
          
          ⁢
          
            
              f
              
                R
                i
              
            
            
              f
              
                T
                i
              
            
          
          ⁢
          
            e
            
              j
              ⁡
              
                (
                
                  
                    φ
                    
                      T
                      0
                    
                  
                  -
                  
                    φ
                    
                      R
                      0
                    
                  
                
                )
              
            
          
          ⁢
          
            e
            
              -
              
                j
                ⁡
                
                  (
                  
                    
                      φ
                      
                        T
                        i
                      
                    
                    -
                    
                      φ
                      
                        R
                        i
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        22
        )
      
    
  




The individual non-reciprocity compensation factor CTi can be calculated using pilot symbols for each individual antenna of, for example, the base station 100B by following the above steps. For example, the individual component CTi in equation (22) are ratios with the numerator being a pilot signal sent from an antenna “0” of the base station and received by antenna “i” of the base and a denominator being a signal sent from the antenna “i” and received by the antenna “0.” The individual non-reciprocity compensation factors CTi may be combined to the non-reciprocity compensation matrix CT, and can be applied by a computing module to, for example, equation 10 to eliminate or at least minimize non-reciprocity of the TRX.

In some embodiments, the pilot symbols described with reference to FIGS. 2 and 3 may be transmitted periodically by, for example, interrupting the flow of the payload data in order to perform a periodic recalibration of the MIMO system 1000. With these embodiments, multiple symbols can be used to accelerate the calibration process while the payload data are not transmitted. In other embodiments, the pilot symbols may be transmitted and received in parallel with transmitting and receiving the payload data, therefore resulting in the calibration that does not interrupt the normal flow of the payload data. In some embodiments, the transmission of the payload data may be reduced, but not interrupted, to optimize the balance between the pilot (calibration) symbol and the payload data.
From the foregoing, it will be appreciated that specific embodiments of the invention have been described herein for purposes of illustration, but that various modifications may be made without deviating from the scope of the various embodiments of the invention. For example, one or more of the UE antennas can be brought in relatively close proximity to the base station antennas to transmit pilot symbols for the base station antennas. Furthermore, while various advantages and features associated with certain embodiments of the disclosure have been described above in the context of those embodiments, other embodiments may also exhibit such advantages and/or features, and not all embodiments need necessarily exhibit such advantages and/or features to fall within the scope of the disclosure. Accordingly, the disclosure is not limited, except as by the appended claims.

## Claims

1. A method for reciprocity calibration of a multiple-input multiple-output (MIMO) system, the method comprising:
receiving, by a base station, a first set of pilot symbols by receivers (RXes) of the base station based on a pilot symbol transmitted from a transmitter (TX) of at least one reference antenna, wherein each of the first set of pilot symbols is received as an ri,0 pilot symbol by an individual antenna of the base station, where index “i” denotes individual antenna “i” of the base station, and index “0” denotes the at least one reference antenna;
transmitting, by the base station, the received first set of pilot symbols by TXes of the base station, wherein the transmitted received first set of pilot symbols are received by an RX of the at least one reference antenna as a second set of r0,i pilot symbols;
calculating non-reciprocity compensation factors based on the first set of pilot symbols and the second set of pilot symbols as





  
    
      r
      
        i
        ,
        0
      
    
    
      r
      
        0
        ,
        i
      
    
  
  ,




and
applying the non-reciprocity compensation factors to payload data.

2. The method of claim 1 wherein the at least one reference antenna is one of the antennas of the base station.

3. The method of claim 1 wherein the at least one reference antenna is an antenna at a user equipment (UE) side.

4. The method of claim 1 wherein the at least one reference antenna is a dedicated auxiliary (AUX) antenna.

5. The method of claim 1, wherein the non-reciprocity compensation factors are applied to payload data of the base station.

6. The method of claim 5 wherein the non-reciprocity compensation factors are applied as rM=HDLCTWS, wherein rM denotes received symbols by the antennas of the base station, HDL denotes an observed transfer function, CT denotes a matrix combination of the non-reciprocity compensation factors, and WS denotes a beamformed payload data.

7. The method of claim 1 wherein the pilot symbol transmitted by the TX of the at least one reference antenna is outside of a frequency bandwidth of a payload data.

8. The method of claim 1 wherein the pilot symbol transmitted by the TX of the at least one reference antenna has a power spectral density that is generally comparable to that of symbols of a payload data.

9. The method of claim 1 wherein the pilot symbol transmitted by the TX of the at least one reference antenna is a first pilot symbol, the method further comprising transmitting a second pilot symbol by the TX of the at least one reference antenna.

10. The method of claim 1, further comprising transmitting a payload data from the TXes of the base station to RXes of a user equipment (UE) side, wherein the pilot symbol and the data payload are transmitted within same time domain and within different frequency bandwidths.

11. The method of claim 10 wherein the pilot symbol and the payload data are within same frequency bandwidth.

12. The method of claim 11 wherein transmitting the payload data pauses while transmitting the pilot symbol.

13. The method of claim 1, further comprising transmitting a payload data from the TXes of the base station to RXes of a user equipment (UE) side, wherein transmitting the payload data pauses when the pilot symbol is transmitted.

14. A multiple-input multiple-output (MIMO) system for wireless signals, comprising:
a reference TX and at least one reference antenna, the reference TX and the at least one reference antenna configured to transmit a pilot symbol;
a base station comprising a plurality of base station antennas, wherein the plurality of base stations antennas and antennas of a user equipment (UE) side are configured to exchange the wireless signals, and wherein RXes of the plurality of base station antennas are configured to receive a first set of pilot symbols based on the pilot symbol transmitted from the reference TX and the at least one reference antenna, wherein each pilot symbol of the first set of pilot symbols is received as the ri,0 pilot symbol of the RX of individual antenna “i” of the base station, and wherein TXes of the plurality of base station antennas are configured to transmit the received first set of pilot symbols, and wherein an RX of the at least one reference antenna is configured to receive a second set of pilot symbols, wherein each pilot symbol of the second set of pilot symbols is received as the r0,i pilot symbol of the RX of the at least one reference antenna, wherein index “i” denotes individual antenna “i” of the base station, and index “0” denotes the at least one reference antenna;
wherein the (MIMO) system is configured to calculate non-reciprocity compensation factors based on the first set of pilot symbols and the second set of pilot symbols as





  
    
      r
      
        i
        ,
        0
      
    
    
      r
      
        0
        ,
        i
      
    
  
  ,




and
wherein the non-reciprocity compensation factors are applied to payload data.

15. The system of claim 14 wherein the at least one reference antenna is one of the antennas of the base station.

16. The system of claim 14 wherein the at least one reference antenna is an antenna of the user equipment (UE) side.

17. The system of claim 14 wherein the at least one reference antenna is a dedicated auxiliary (AUX) antenna.

18. The system of claim 14 wherein the computing module is further configured to apply the non-reciprocity compensation factors as rM=HDLCTWS, wherein rM denotes received symbols by the antennas of the base station, HDL denotes an observed transfer function, CT denotes a matrix combination of the non-reciprocity compensation factors, and WS denotes a beamformed payload data.

19. The system of claim 14 wherein the pilot symbol is one of a plurality of pilot symbols.

20. The system of claim 14 wherein the pilot symbol transmitted by the TX of the at least one reference antenna is outside of a frequency bandwidth of a payload data.

21. The system of claim 14 wherein the pilot symbol transmitted by the TX of the at least one reference antenna and a payload data are within same frequency bandwidth.

22. The system of claim 14 wherein the pilot symbol transmitted by the TX of the at least one reference antenna has a power spectral density that is generally comparable to that of symbols of a payload data.