# Patent

## Title
Method and apparatus for beamforming with coupled antennas

## Abstract

A method and apparatus provide beamforming with coupled antennas. A first channel matrix can be transformed into a transformed second channel matrix. The first channel matrix can be a channel from a transmitter array of antennas of the transmitting device to at least one receiver antenna of a receiving device. A precoding matrix can be determined that maximizes the capacity for the transformed second channel matrix subject to a power constraint of the precoding matrix. The channel capacity maximizing precoding matrix for the transformed second channel matrix can be converted into an optimal precoding matrix for the first channel matrix. A signal can be received for transmission. The optimal precoding matrix can be applied to the signal to generate a precoded signal for transmission over a physical channel. The precoded signal can then be transmitted.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is related to an application entitled “Method and Apparatus for Equal Energy Codebooks for Antenna Arrays with Mutual Coupling,” U.S. application Ser. No. 14/855,693, filed on Sep. 16, 2015, Motorola Mobility docket number MM01601, and commonly assigned to the assignee of the present application, which is hereby incorporated by reference in its entirety.
BACKGROUND
1. Field
The present disclosure is directed to a method and apparatus for beamforming with coupled antennas.
2. Introduction
Presently, wireless communication devices communicate with other communication devices using wireless signals. Many wireless communication devices have multiple antennas that can transmit more focused signals to a receiving device using antenna beamforming. For a transmitter with M antennas transmitting to a receiver with a single antenna, an optimal beamformer maximizes the signal-to-interference ratio at the receiver subject to a constraint on transmit power. The solution for maximizing the signal-to-interference ratio is known as the eigenbeamformer due to the fact that the optimal transmit weights correspond to the eigenvector of a matrix having the largest eigenvalue. While this method is known, the impact of antenna coupling at the transmitter has typically not been considered in the formulation. However, the effect of the antenna correlation must be considered in the formulation of the eigenbeamformer as it affects the computation of the radiated power, and thus the constraint on radiated power.
The effect of mutual coupling has been considered to some extent. U.S. Pat. No. 7,457,590, (issued Nov. 25, 2008), which is hereby incorporated by reference, addresses how to compute the eigenbeamformer as a function of the channel received angular power profile ƒ(θ) in the case that the antenna pattern pk(θ) resulting from each antenna weighting of the form




  
    u
    k
  
  =
  
    
      
        [
        
          
            
              0
              ⁢
              ...0
            
            
              }
              
                k
                -
                1
              
            
          
          ⁢
          10
          ⁢
          ...0
        
        ]
      
      T
    
    .
  




In “Mutual Coupling in MIMO Wireless Systems: A Rigorous Network Theory Analysis” (Wallace, IEEE, p. 1317-1325, Jul. 4, 2004), which is hereby incorporated by reference, the eigenbeamformer is not explicitly calculated. However, the impact of mutual coupling on the calculation of the radiated power is addressed in computing the MIMO channel capacity in the case of mutual coupling at the transmitter, though only for the case in which the antennas are driven by ideal current sources. Unfortunately, the technique does not address coupled antennas driven by actual sources.
Thus, there is a need for a method and apparatus for improved beamforming with coupled antennas.



BRIEF DESCRIPTION OF THE DRAWINGS
In order to describe the manner in which advantages and features of the disclosure can be obtained, a description of the disclosure is rendered by reference to specific embodiments thereof which are illustrated in the appended drawings. These drawings depict only example embodiments of the disclosure and are not therefore to be considered to be limiting of its scope.
FIG. 1 is an example block diagram of a system according to a possible embodiment;
FIG. 2 is an example illustration of a two-port model for a two-element array according to a possible embodiment;
FIG. 3 is an example illustration of a Thevenin source model according to a possible embodiment;
FIG. 4 is an example illustration of a Norton source model according to a possible embodiment;
FIG. 5 is an example flowchart illustrating the operation of a transmitting device according to a possible embodiment; and
FIG. 6 is an example block diagram of an apparatus according to a possible embodiment.



DETAILED DESCRIPTION
Embodiments provide a method and apparatus for beamforming with coupled antennas. According to a possible embodiment, a first channel matrix H can be transformed into a transformed second channel matrix P−TSourceH. The first channel matrix H can be a channel from a transmitter array of antennas of the transmitting device to at least one receiver antenna of a receiving device. A precoding matrix W can be determined that maximizes the capacity for the transformed second channel matrix P−TSourceH subject to a power constraint of the precoding matrix W. The channel capacity maximizing precoding matrix W for the transformed second channel matrix can be converted into an optimal precoding matrix V for the first channel matrix. A signal can be received for transmission. The optimal precoding matrix V can be applied to the signal to generate a precoded signal for transmission over a physical channel. The precoded signal can then be transmitted.
FIG. 1 is an example block diagram of a system 100 according to a possible embodiment. The system 100 can include a transmitting device 110 and a receiving device 120. The transmitting device 110 can be a User Equipment (UE), a base station, an access point, or any other device that can transmit wireless signals. Similarly, the receiving device 120 can be a UE, a base station, an access point, or any other device that can receive wireless signals. A UE can be a wireless terminal, a portable wireless communication device, a smartphone, a cellular telephone, a flip phone, a personal digital assistant, a device having a subscriber identity module, a personal computer, a selective call receiver, a tablet computer, a laptop computer, or any other device that is capable of sending and receiving wireless communication signals.
The transmitting device 110 can include a precoding matrix determination controller 112, a codebook 114, and an antenna array 116. The precoding matrix determination controller 112 can be one element or can be distributed between different elements. For example, the precoding matrix determination controller 112 can be part of a processor, can be part of a transceiver, can be part of a precoder, can be part of other elements in a transmitting device, and/or can be distributed between combinations of elements in a transmitting device and/or over cloud computing. The receiving device 120 can include at least one antenna 122. For example, in some embodiments the receiving device 120 can have one antenna and in other embodiments the receiving device 120 can have an array of antennas.
In operation, a first channel matrix H can be transformed into a transformed second channel matrix P−TSourceH. The first channel matrix H can be a channel from a transmitter array of antennas 116 of the transmitting device 110 to at least one receiver antenna 122 of the receiving device 120. A precoding matrix W can be determined that maximizes the capacity for the transformed second channel matrix P−TSourceH subject to a power constraint of the precoding matrix W. The channel capacity maximizing precoding matrix W for the transformed second channel matrix can be converted into an optimal precoding matrix V for the first channel matrix. A signal can be received for transmission. The optimal precoding matrix V can be applied to the signal to generate a precoded signal for transmission over a physical channel. The precoded signal can then be transmitted.
For example, embodiments can provide for optimal beamforming with coupled antennas of the antenna array 116 and linear source models. A circuit model can be used to model the mutual coupling of the antennas in order to compute an eigenbeamformer for an arbitrary channel matrix H. As the eigenbeamformer can be a function of the source used to drive the antenna array, the eigenbeamformer can be computed for both Thevenin and Norton source models.
FIG. 2 is an example illustration of a two-port model 200 for a two-element array according to a possible embodiment. The two ports can correspond to two antennas in an antenna array. An M-port circuit can be used to model the vector voltage-current relationship for the M-ports of the M-element antenna array, which can be given by

V=ZI, 

where Z can be the M×M impedance matrix for the array. In the two port model 200, i1 and v1 denote the current and voltage for the first antenna, while i2 and v2 denote the current and voltage for the second antenna.

FIG. 3 is an example illustration of a Thevenin source model 300 according to a possible embodiment. FIG. 4 is an example illustration of a Norton source model 400 according to a possible embodiment. The two linear source models 300 and 400 can be considered for driving an antenna array. The Thevenin source model 300 can include ideal voltage sources vS in combination with series impedances ZS_Thev, while the Norton source model 400 can include an ideal current source iS in combination with a parallel shunt impedance ZS_Nor.
The Norton source can yield two-port currents i and voltages v which can be equal to that for the Thevenin source so long as

ZS_Thev=ZS_Nor and vS=ZS_TheviS 

where





  
    
      
        
          v
          S
        
        =
        
          [
          
            
              
                
                  v
                  
                    S
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
            
            
              
                
                  v
                  
                    S
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
              
            
          
          ]
        
      
    
    
      
        
          Z
          
            S
            ⁢
            _
            ⁢
            Thev
          
        
        =
        
          [
          
            
              
                
                  Z
                  
                    
                      S
                      ⁢
                      _
                      ⁢
                      Thev
                    
                    ,
                    1
                  
                
              
              
                0
              
            
            
              
                0
              
              
                
                  Z
                  
                    
                      S
                      ⁢
                      _
                      ⁢
                      Thev
                    
                    ,
                    2
                  
                
              
            
          
          ]
        
      
    
  




and





  
    
      
        
          
            i
            S
          
          =
          
            [
            
              
                
                  
                    i
                    
                      S
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                  
                
              
              
                
                  
                    i
                    
                      S
                      ⁢
                      
                        
                      
                      ⁢
                      2
                    
                  
                
              
            
            ]
          
        
      
      
        
          
            Z
            
              S
              ⁢
              _
              ⁢
              Nor
            
          
          =
          
            [
            
              
                
                  
                    Z
                    
                      
                        S
                        ⁢
                        _
                        ⁢
                        Nor
                      
                      ,
                      1
                    
                  
                
                
                  0
                
              
              
                
                  0
                
                
                  
                    Z
                    
                      
                        S
                        ⁢
                        _Nor
                      
                      ,
                      2
                    
                  
                
              
            
            ]
          
        
      
    
  
  .




For an M-element array, the peak radiated power (average power is one-half of peak) can be equal to the power delivered to the M-port device and can be given by

Re(vHi)=Re(iHZHi),

where Z can be the impedance matrix and i can be the vector of input currents. For the Thevenin source 300 with source voltage vS and source impedance ZS_Thev, the current vector at the input to the two-port device can be given by

i=(ZS_Thev+Z)−1vS.

Thus, the radiated power for the Thevenin source can be given by




  
    
      
        
          
            P
            
              rad
              ⁢
              _
              ⁢
              Thev
            
          
          ⁡
          
            (
            
              
                v
                S
              
              ,
              
                Z
                
                  S
                  ⁢
                  _
                  ⁢
                  Thev
                
              
              ,
              Z
            
            )
          
        
        =
        
        ⁢
        
          Re
          (
          
            
              (
              
                
                  
                    (
                    
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      +
                      Z
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  v
                  S
                
              
              )
            
            H
          
        
      
    
  
  
    
      
        
        ⁢
        
          Z
          ⁡
          
            (
            
              
                
                  (
                  
                    
                      Z
                      S_Thev
                    
                    +
                    Z
                  
                  )
                
                
                  -
                  1
                
              
              ⁢
              
                v
                S
              
            
            )
          
        
        )
      
    
  
  
    
      
        =
        
        ⁢
        
          Re
          (
          
            
              
                v
                S
                H
              
              ⁡
              
                (
                
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Thev
                    
                  
                  +
                  Z
                
                )
              
            
            
              -
              H
            
          
        
      
    
  
  
    
      
        
          
          ⁢
          
            
              
                Z
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Thev
                      
                    
                    +
                    Z
                  
                  )
                
              
              
                -
                1
              
            
            ⁢
            
              v
              S
            
          
          )
        
        .
      
    
  




This expression can be further simplified as




  
    
      
        
          
            P
            
              rad
              ⁢
              _
              ⁢
              Thev
            
          
          ⁡
          
            (
            
              
                v
                S
              
              ,
              
                Z
                
                  S
                  ⁢
                  _
                  ⁢
                  Thev
                
              
              ,
              Z
            
            )
          
        
        =
        
        ⁢
        
          
            1
            2
          
          ⁢
          
            (
            
              
                
                  
                    
                      
                        
                          
                            v
                            S
                            H
                          
                          ⁡
                          
                            (
                            
                              
                                Z
                                
                                  S
                                  ⁢
                                  _
                                  ⁢
                                  Thev
                                
                              
                              +
                              Z
                            
                            )
                          
                        
                        
                          -
                          H
                        
                      
                      ⁢
                      
                        
                          Z
                          ⁡
                          
                            (
                            
                              
                                Z
                                
                                  S
                                  ⁢
                                  _
                                  ⁢
                                  Thev
                                
                              
                              +
                              Z
                            
                            )
                          
                        
                        
                          -
                          1
                        
                      
                      ⁢
                      
                        v
                        S
                      
                    
                    +
                  
                
              
              
                
                  
                    
                      (
                      
                        
                          
                            
                              v
                              S
                              H
                            
                            ⁡
                            
                              (
                              
                                
                                  Z
                                  
                                    S
                                    ⁢
                                    _
                                    ⁢
                                    Thev
                                  
                                
                                +
                                Z
                              
                              )
                            
                          
                          
                            -
                            H
                          
                        
                        ⁢
                        
                          
                            Z
                            ⁡
                            
                              (
                              
                                
                                  Z
                                  
                                    S
                                    ⁢
                                    _
                                    ⁢
                                    Thev
                                  
                                
                                +
                                Z
                              
                              )
                            
                          
                          
                            -
                            1
                          
                        
                        ⁢
                        
                          v
                          S
                        
                      
                      )
                    
                    H
                  
                
              
            
            )
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            1
            2
          
          ⁢
          
            (
            
              
                
                  
                    v
                    S
                    H
                  
                  ⁡
                  
                    (
                    
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      +
                      Z
                    
                    )
                  
                
                
                  -
                  H
                
              
              ⁢
              
                (
                
                  Z
                  +
                  
                    Z
                    H
                  
                
                )
              
              ⁢
              
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Thev
                      
                    
                    +
                    Z
                  
                  )
                
                
                  -
                  1
                
              
              ⁢
              
                v
                S
              
            
            )
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            
              v
              S
              H
            
            ⁡
            
              (
              
                
                  
                    (
                    
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      +
                      Z
                    
                    )
                  
                  
                    -
                    H
                  
                
                ⁢
                
                  Re
                  ⁡
                  
                    (
                    Z
                    )
                  
                
                ⁢
                
                  
                    (
                    
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      +
                      Z
                    
                    )
                  
                  
                    -
                    1
                  
                
              
              )
            
          
          ⁢
          
            v
            S
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            v
            S
            H
          
          ⁢
          
            Q
            Thev
          
          ⁢
          
            v
            S
          
        
      
    
  




where

QThev=(ZS_Thev+Z)−HRe(Z)(ZS_Thev+Z)−1,

and where for passive linear networks,

Z=ZT

and thus,





  
    
      Z
      +
      
        Z
        H
      
    
    2
  
  =
  
    
      
        Z
        +
        
          
            (
            
              Z
              T
            
            )
          
          *
        
      
      2
    
    =
    
      
        
          Z
          +
          
            Z
            *
          
        
        2
      
      =
      
        
          Re
          ⁡
          
            (
            Z
            )
          
        
        .
      
    
  




For the Norton source model 400 with source currents iS and source impedance ZS_Nor, the antenna currents can be given by

i=Z−1(ZS_Nor−1+Z−1)−1iS,

where Z can be the impedance matrix for the array. Assuming that all of the power delivered to the array is radiated (i.e., no ohmic or other losses), the radiated power for the Norton source can be given by





  
    
      
        
          
            P
            
              
                rad
                ⁢
                _
                ⁢
                Nor
              
              ⁢
              
                _
                ⁢
                circuit
              
            
          
          ⁡
          
            (
            
              
                i
                S
              
              ,
              
                Z
                
                  S
                  ⁢
                  _
                  ⁢
                  Nor
                
              
              ,
              Z
            
            )
          
        
        =
        
        ⁢
        
          Re
          (
          
            
              
                
                  i
                  S
                  H
                
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Nor
                      
                      
                        -
                        1
                      
                    
                    +
                    
                      Z
                      
                        -
                        1
                      
                    
                  
                  )
                
              
              
                -
                H
              
            
            ⁢
            
              Z
              
                -
                H
              
            
            ⁢
            
              ZZ
              
                -
                1
              
            
          
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                
                  Z
                  
                    S
                    ⁢
                    _
                    ⁢
                    Nor
                  
                  
                    -
                    1
                  
                
                +
                
                  Z
                  
                    -
                    1
                  
                
              
              )
            
            
              -
              1
            
          
          ⁢
          
            i
            S
          
        
        )
      
    
  
  
    
      
        =
        
        ⁢
        
          Re
          (
          
            
              i
              S
              H
            
            ⁢
            
              
                
                  Z
                  
                    S
                    ⁢
                    _
                    ⁢
                    Nor
                  
                  H
                
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Nor
                      
                    
                    +
                    Z
                  
                  )
                
              
              
                -
                H
              
            
          
        
      
    
  
  
    
      
        
          
          ⁢
          
            
              
                Z
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Nor
                      
                    
                    +
                    Z
                  
                  )
                
              
              
                -
                1
              
            
            ⁢
            
              Z
              
                S
                ⁢
                _
                ⁢
                Nor
              
            
            ⁢
            
              i
              S
            
          
          )
        
        .
      
    
  




If the Norton source is equivalent to the Thevenin source in the previous section, then

ZS_NoriS=vS and ZS_Nor=ZS_Thev 

and the radiated power can be given by





  
    
      
        
          
            P
            
              
                rad
                ⁢
                _
                ⁢
                Nor
              
              ⁢
              
                _
                ⁢
                circuit
              
            
          
          ⁡
          
            (
            
              
                i
                S
              
              ,
              
                Z
                
                  S
                  ⁢
                  _
                  ⁢
                  Nor
                
              
              ,
              Z
            
            )
          
        
        =
        
        ⁢
        
          Re
          (
          
            
              i
              S
              H
            
            ⁢
            
              
                
                  Z
                  
                    S
                    ⁢
                    _
                    ⁢
                    Nor
                  
                  H
                
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Nor
                      
                    
                    +
                    Z
                  
                  )
                
              
              
                -
                H
              
            
          
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              Z
              ⁡
              
                (
                
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Nor
                    
                  
                  +
                  Z
                
                )
              
            
            
              -
              1
            
          
          ⁢
          
            Z
            
              S
              ⁢
              _
              ⁢
              Nor
            
          
          ⁢
          
            i
            S
          
        
        )
      
    
  
  
    
      
        =
        
        ⁢
        
          Re
          (
          
            
              
                v
                S
                H
              
              ⁡
              
                (
                
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Thev
                    
                  
                  +
                  Z
                
                )
              
            
            
              -
              H
            
          
        
      
    
  
  
    
      
        
          
          ⁢
          
            
              
                Z
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Thev
                      
                    
                    +
                    Z
                  
                  )
                
              
              
                -
                1
              
            
            ⁢
            
              v
              S
            
          
          )
        
        ,
      
    
  




which can be the same as for the Thevenin source model 300.

As in the case of the Thevenin source model 300, the expression for radiated power can be simplified as




  
    
      
        
          
            P
            
              
                rad
                ⁢
                _
                ⁢
                Nor
              
              ⁢
              
                _
                ⁢
                circuit
              
            
          
          ⁡
          
            (
            
              
                i
                S
              
              ,
              
                Z
                
                  S
                  ⁢
                  _
                  ⁢
                  Nor
                
              
              ,
              Z
            
            )
          
        
        =
        
        ⁢
        
          
            1
            2
          
          ⁢
          
            (
            
              
                
                  
                    
                      i
                      S
                      H
                    
                    ⁢
                    
                      
                        
                          Z
                          
                            S
                            ⁢
                            _
                            ⁢
                            Nor
                          
                          H
                        
                        ⁡
                        
                          (
                          
                            
                              Z
                              
                                S
                                ⁢
                                _
                                ⁢
                                Nor
                              
                            
                            +
                            Z
                          
                          )
                        
                      
                      
                        -
                        H
                      
                    
                    ⁢
                    
                      
                        Z
                        ⁡
                        
                          (
                          
                            
                              Z
                              
                                S
                                ⁢
                                _
                                ⁢
                                Nor
                              
                            
                            +
                            Z
                          
                          )
                        
                      
                      
                        -
                        1
                      
                    
                  
                
              
              
                
                  
                    
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Nor
                        
                      
                      ⁢
                      
                        i
                        S
                      
                    
                    +
                    
                      (
                      
                        
                          i
                          S
                          H
                        
                        ⁢
                        
                          
                            
                              Z
                              
                                S
                                ⁢
                                _
                                ⁢
                                Nor
                              
                              H
                            
                            ⁡
                            
                              (
                              
                                
                                  Z
                                  
                                    S
                                    ⁢
                                    _
                                    ⁢
                                    Nor
                                  
                                
                                +
                                Z
                              
                              )
                            
                          
                          
                            -
                            H
                          
                        
                      
                    
                  
                
              
              
                
                  
                    
                      
                        
                          
                            Z
                            ⁡
                            
                              (
                              
                                
                                  Z
                                  
                                    S
                                    ⁢
                                    _
                                    ⁢
                                    Nor
                                  
                                
                                +
                                Z
                              
                              )
                            
                          
                          
                            -
                            1
                          
                        
                        ⁢
                        
                          Z
                          
                            S
                            ⁢
                            _
                            ⁢
                            Nor
                          
                        
                        ⁢
                        
                          i
                          S
                        
                      
                      )
                    
                    H
                  
                
              
            
            )
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            1
            2
          
          ⁢
          
            (
            
              
                i
                S
                H
              
              ⁢
              
                
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Nor
                    
                    H
                  
                  ⁡
                  
                    (
                    
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Nor
                        
                      
                      +
                      Z
                    
                    )
                  
                
                
                  -
                  H
                
              
              ⁢
              
                (
                
                  Z
                  +
                  
                    Z
                    H
                  
                
                )
              
            
          
        
      
    
  
  
    
      
        
        ⁢
        
          
            
              (
              
                
                  Z
                  
                    S
                    ⁢
                    _
                    ⁢
                    Nor
                  
                
                +
                Z
              
              )
            
            
              -
              1
            
          
          ⁢
          
            Z
            
              S
              ⁢
              _
              ⁢
              Nor
            
          
          ⁢
          
            i
            S
          
        
        )
      
    
  
  
    
      
        =
        
        ⁢
        
          
            i
            S
            H
          
          (
          
            
              
                
                  Z
                  
                    S
                    ⁢
                    _
                    ⁢
                    Nor
                  
                  H
                
                ⁡
                
                  (
                  
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Nor
                      
                    
                    +
                    Z
                  
                  )
                
              
              
                -
                H
              
            
            ⁢
            
              Re
              ⁡
              
                (
                Z
                )
              
            
          
        
      
    
  
  
    
      
        
          
          ⁢
          
            
              
                (
                
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Nor
                    
                  
                  +
                  Z
                
                )
              
              
                -
                1
              
            
            ⁢
            
              Z
              
                S
                ⁢
                _
                ⁢
                Nor
              
            
          
          )
        
        ⁢
        
          i
          S
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            i
            S
            H
          
          ⁢
          
            Q
            Nor
          
          ⁢
          
            i
            S
          
        
      
    
  




where

QNor=ZS—NorH(ZS_Nor+Z)−HRe(Z)(ZS_Nor+Z)−1ZS_Nor.

For a transmitter with a Thevenin source and a receiver with a single receive antenna, the M×1 vector hThev(ZS_Thev) given by

hThev(ZS_Thev)=[hThev,0(ZS_Thev) . . . hThev,k(ZS_Thev) . . . hThev,M−1(ZS_Thev)]T 

can denote the channel observed at the receiver from each transmit antenna from a Thevenin source with series impedance ZS_Thev. More precisely, let hThev,k(ZS_Thev) denote the channel observed at the receiver when the voltage source vector





  
    v
    k
  
  =
  
    
      [
      
        
          
            0
            ⁢
            ...0
          
          
            }
            
              k
              -
              1
            
          
        
        ⁢
        10
        ⁢
        ...0
      
      ]
    
    T
  




is applied at the transmitter. The fact that the channel observed at the receiver can depend on the source impedance ZS_Thev can be due to the fact that the antenna pattern that results when the precoder vk is applied to the array is given by

vkT(ZS_Thev+Z)−1p(θ,φ),

where

p(θ,φ)=[p1(θ,φ),p2(θ,φ), . . . , pM(θ,φ)]T,

and pk(θ) can be the antenna pattern for the k-th antenna element when all of the other elements are removed from the array. For an arbitrary voltage source vector v, the signal y observed at the receiver can be given by

y=vThThev(ZS_Thev)+n, 

where n can denote a zero-mean complex Gaussian random variable with variance σ2.

The signal-to-noise ratio at the receiver can then be given by




  
    
      
        
        
          E
          ⁡
          
            (
            y
            )
          
        
        
      
      2
    
    
      σ
      2
    
  
  =
  
    
      
        
          
          
            
              v
              T
            
            ⁢
            
              
                h
                Thev
              
              ⁡
              
                (
                
                  Z
                  
                    S
                    ⁢
                    _
                    ⁢
                    Thev
                  
                
                )
              
            
          
          
        
        2
      
      
        σ
        2
      
    
    =
    
      
        
          
            
              v
              T
            
            ⁡
            
              (
              
                
                  
                    h
                    Thev
                  
                  ⁡
                  
                    (
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Thev
                      
                    
                    )
                  
                
                ⁢
                
                  
                    h
                    Thev
                    H
                  
                  ⁡
                  
                    (
                    
                      Z
                      
                        S
                        ⁢
                        _
                        ⁢
                        Thev
                      
                    
                    )
                  
                
              
              )
            
          
          ⁢
          
            v
            *
          
        
        
          σ
          2
        
      
      =
      
        
          
            
              
                v
                H
              
              ⁡
              
                (
                
                  
                    
                      h
                      Thev
                      *
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      )
                    
                  
                  ⁢
                  
                    
                      h
                      Thev
                      T
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      )
                    
                  
                
                )
              
            
            ⁢
            v
          
          
            σ
            2
          
        
        .
      
    
  




Thus, the optimal beamformer vopt can maximize

vH(h*Thev(ZS_Thev)hThevT(ZS_Thev))v 

subject to the constraint that

vHQThevv=1.

Since the matrix QThev is Hermitian and positive definite, this matrix can be factored as

QThev=PThevHPThev 

where we refer to PThevH and PThev as the left and right factors of the product.

If we define

w=PThevv, 

it then follows that

v=PThev−1w. 

With this notation, the objective of the beamformer can now be to maximize

wH(PThev−1)H(h*Thev(ZS_Thev)hThevT(ZS_Thev))PThev−1w 

subject to the constraint that

wHw=1.

The solution to this problem is known and can be obtained by setting w

w=eThev_max,

where eThev_max can be the eigenvector corresponding to the largest eigenvalue, λThev_max, of the matrix (PThev−1)H(h*Thev(ZS_Thev)hThevT(ZS_Thev))PThev−1. The optimal beamformer vopt can then be given by

vopt=PThev−1eThev_max.

The signal-to-noise ratio of the optimal beamformer vopt can be given by




  
    
      SNR
      Thev
    
    ⁡
    
      (
      
        v
        opt
      
      )
    
  
  =
  
    
      
        λ
        
          Thev
          ⁢
          _
          ⁢
          max
        
      
      
        σ
        2
      
    
    .
  




For the gain of the optimal beamformer relative to that of any other beamformer v for which the radiated power is normalized to unity, the required normalization of v can be given by




  
    v
    
      
        (
        
          
            v
            H
          
          ⁢
          
            Q
            Thev
          
          ⁢
          v
        
        )
      
      
        1
        /
        2
      
    
  
  .




The received signal-to-noise ratio for this beamformer can be given by




  
    
      SNR
      Thev
    
    ⁡
    
      (
      
        v
        
          
            (
            
              
                v
                H
              
              ⁢
              
                Q
                Thev
              
              ⁢
              v
            
            )
          
          
            1
            /
            2
          
        
      
      )
    
  
  =
  
    
      
        
          
            v
            H
          
          ⁡
          
            (
            
              
                
                  h
                  Thev
                  *
                
                ⁡
                
                  (
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Thev
                    
                  
                  )
                
              
              ⁢
              
                
                  h
                  Thev
                  T
                
                ⁡
                
                  (
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Thev
                    
                  
                  )
                
              
            
            )
          
        
        ⁢
        v
      
      
        
          (
          
            
              v
              H
            
            ⁢
            
              Q
              Thev
            
            ⁢
            v
          
          )
        
        ⁢
        
          σ
          2
        
      
    
    .
  




Thus, the gain of the optimal beamformer vopt to an arbitrary beamformer v can be given by




  
    
      
        SNR
        Thev
      
      ⁡
      
        (
        
          v
          opt
        
        )
      
    
    
      
        SNR
        Thev
      
      ⁡
      
        (
        
          
            (
            
              
                v
                H
              
              ⁢
              
                Q
                Thev
              
              ⁢
              v
            
            )
          
          
            
              -
              1
            
            /
            2
          
        
        )
      
    
  
  =
  
    
      
        
          λ
          
            Thev
            ⁢
            _
            ⁢
            max
          
        
        
          σ
          2
        
      
      ⁢
      
        
          (
          
            
              
                
                  v
                  H
                
                ⁡
                
                  (
                  
                    
                      
                        h
                        Thev
                        *
                      
                      ⁡
                      
                        (
                        
                          Z
                          
                            S
                            ⁢
                            _
                            ⁢
                            Thev
                          
                        
                        )
                      
                    
                    ⁢
                    
                      
                        h
                        Thev
                        T
                      
                      ⁡
                      
                        (
                        
                          Z
                          
                            S
                            ⁢
                            _
                            ⁢
                            Thev
                          
                        
                        )
                      
                    
                  
                  )
                
              
              ⁢
              v
            
            
              
                (
                
                  
                    v
                    H
                  
                  ⁢
                  
                    Q
                    Thev
                  
                  ⁢
                  v
                
                )
              
              ⁢
              
                σ
                2
              
            
          
          )
        
        
          -
          1
        
      
    
    =
    
      
        λ
        
          Thev
          ⁢
          _
          ⁢
          max
        
      
      ⁢
      
        
          
            
              v
              H
            
            ⁢
            
              Q
              Thev
            
            ⁢
            v
          
          
            
              
                v
                H
              
              ⁡
              
                (
                
                  
                    
                      h
                      Thev
                      *
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      )
                    
                  
                  ⁢
                  
                    
                      h
                      Thev
                      T
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Thev
                        
                      
                      )
                    
                  
                
                )
              
            
            ⁢
            v
          
        
        .
      
    
  




For a transmitter with a Norton source and a receiver with a single receive antenna, let the M×1 vector hNor(ZS_Nor) given by

hNor(ZS_Nor)=[hNor,0(ZS_Nor) . . . hNor,k(ZS_Nor) . . . hNor,M−1(ZS_Nor)]T 

denote the channel observed at the receiver from each transmit antenna from a Norton source with shunt impedance ZS_Nor. More precisely, let hNor,k(ZS_Nor) denote the channel observed at the receiver when the current source vector





  
    i
    k
  
  =
  
    
      [
      
        
          
            0
            ⁢
            ...0
          
          
            }
            
              k
              -
              1
            
          
        
        ⁢
        10
        ⁢
        ...0
      
      ]
    
    T
  




is applied at the transmitter. The fact that the channel observed at the receiver depends on the source impedance ZS_Nor can be due to the fact that the antenna pattern that results when the precoder ik is applied to the array is given by

iST(ZS_Nor−1+Z−1)−1Z−1p(θ,φ).

where

p(θ,φ)=[p1(θ,φ),p2(θ,φ), . . . , pM(θ,φ)]T.

and pk(θ) can be the antenna pattern for the k-th antenna element when all of the other elements are removed from the array. Thus, for an arbitrary current source vector i, the signal y observed at the receiver can be given by

y=iThNor(ZS_Nor)+n, 

where n denotes a zero-mean complex Gaussian random variable with variance σ2. The signal-to-noise ratio at the receiver can be then given by





  
    
      
        
          
            
              
              
                E
                ⁡
                
                  (
                  y
                  )
                
              
              
            
            2
          
          
            σ
            2
          
        
        =
        
          
            
              
                
                
                  
                    i
                    T
                  
                  ⁢
                  
                    
                      h
                      Nor
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Nor
                        
                      
                      )
                    
                  
                
                
              
              2
            
            
              σ
              2
            
          
          =
          
            
              
                
                  
                    i
                    T
                  
                  ⁡
                  
                    (
                    
                      
                        
                          h
                          Nor
                        
                        ⁡
                        
                          (
                          
                            Z
                            
                              S
                              ⁢
                              _
                              ⁢
                              Nor
                            
                          
                          )
                        
                      
                      ⁢
                      
                        
                          h
                          Nor
                          H
                        
                        ⁡
                        
                          (
                          
                            Z
                            
                              S
                              ⁢
                              _
                              ⁢
                              Nor
                            
                          
                          )
                        
                      
                    
                    )
                  
                
                ⁢
                
                  i
                  *
                
              
              
                σ
                2
              
            
            =
            
              
                
                  
                    
                      i
                      H
                    
                    ⁡
                    
                      (
                      
                        
                          
                            h
                            Nor
                            *
                          
                          ⁡
                          
                            (
                            
                              Z
                              
                                S
                                ⁢
                                _
                                ⁢
                                Nor
                              
                            
                            )
                          
                        
                        ⁢
                        
                          
                            h
                            Nor
                            T
                          
                          ⁡
                          
                            (
                            
                              Z
                              
                                S
                                ⁢
                                _
                                ⁢
                                Nor
                              
                            
                            )
                          
                        
                      
                      )
                    
                  
                  ⁢
                  i
                
                
                  σ
                  2
                
              
              .
            
          
        
      
    
    
      
        
      
    
  




Thus, the optimal beamformer iopt can maximize

iH(h*Nor(ZS_Nor)hNorT(ZS_Nor))i 

subject to the constraint that

iHQNori=1

Since the matrix QNor is Hermitian and positive definite, this matrix can be factored as

QNor=PNorHPNor 

where we refer to PNorH and PNor as the left and right factors of the product.

If we define

w=PNori, 

it then follows that

i=PNor−1w. 

With this notation, the objective of the beamformer can now be to maximize

wH(PNor−1)H(h*Nor(ZS_Nor)hNorT(ZS_Nor))PNor−1w 

subject to the constraint that

wHw=1.

The solution to this problem is known and can be obtained by

w=eNor_max,

where eNor_max can be the eigenvector corresponding to the largest eigenvalue, λNor_max, of the matrix (PNor−1)H(h*Nor(ZS_Nor)hNorT(ZS_Nor))PNor−1. The optimal beamformer iopt can then be given by

iopt=PNor−1eNor_max.

The signal-to-noise ratio of the optimal beamformer iopt can be given by




  
    SNR
    Nor
  
  =
  
    
      (
      
        i
        opt
      
      )
    
    =
    
      
        
          λ
          
            Nor
            ⁢
            _
            ⁢
            max
          
        
        
          σ
          2
        
      
      .
    
  




When the gain of the optimal beamformer relative to that of any other beamformer i for which the radiated power is normalized to unity, the required normalization of i can be given by




  
    i
    
      
        (
        
          
            i
            H
          
          ⁢
          
            Q
            Nor
          
          ⁢
          i
        
        )
      
      
        1
        /
        2
      
    
  
  .




The received signal-to-noise ratio for this beamformer can be given by




  
    
      SNR
      Nor
    
    ⁡
    
      (
      
        
          (
          
            
              v
              H
            
            ⁢
            
              Q
              Nor
            
            ⁢
            i
          
          )
        
        
          
            -
            1
          
          /
          2
        
      
      )
    
  
  =
  
    
      
        
          
            v
            H
          
          ⁡
          
            (
            
              
                
                  h
                  Nor
                  *
                
                ⁡
                
                  (
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Nor
                    
                  
                  )
                
              
              ⁢
              
                
                  h
                  Nor
                  T
                
                ⁡
                
                  (
                  
                    Z
                    
                      S
                      ⁢
                      _
                      ⁢
                      Nor
                    
                  
                  )
                
              
            
            )
          
        
        ⁢
        v
      
      
        
          (
          
            
              i
              H
            
            ⁢
            
              Q
              Nor
            
            ⁢
            i
          
          )
        
        ⁢
        
          σ
          2
        
      
    
    .
  




Thus, the gain of the optimal beamformer iopt relative to an arbitrary beamformer i can be given by




  
    
      
        SNR
        Nor
      
      ⁡
      
        (
        
          i
          opt
        
        )
      
    
    
      
        SNR
        Nor
      
      ⁡
      
        (
        
          
            (
            
              
                i
                H
              
              ⁢
              
                Q
                Thev
              
              ⁢
              i
            
            )
          
          
            
              -
              1
            
            /
            2
          
        
        )
      
    
  
  =
  
    
      
        
          λ
          
            Nor
            ⁢
            _
            ⁢
            max
          
        
        
          σ
          2
        
      
      ⁢
      
        
          (
          
            
              
                
                  i
                  H
                
                ⁡
                
                  (
                  
                    
                      
                        h
                        Nor
                        *
                      
                      ⁡
                      
                        (
                        
                          Z
                          
                            S
                            ⁢
                            _
                            ⁢
                            Nor
                          
                        
                        )
                      
                    
                    ⁢
                    
                      
                        h
                        Nor
                        T
                      
                      ⁡
                      
                        (
                        
                          Z
                          
                            S
                            ⁢
                            _
                            ⁢
                            Nor
                          
                        
                        )
                      
                    
                  
                  )
                
              
              ⁢
              i
            
            
              
                (
                
                  
                    i
                    H
                  
                  ⁢
                  
                    Q
                    Nor
                  
                  ⁢
                  i
                
                )
              
              ⁢
              
                σ
                2
              
            
          
          )
        
        
          -
          1
        
      
    
    =
    
      
        λ
        
          Nor
          ⁢
          _
          ⁢
          max
        
      
      ⁢
      
        
          
            
              i
              H
            
            ⁢
            
              Q
              Nor
            
            ⁢
            i
          
          
            
              
                i
                H
              
              ⁡
              
                (
                
                  
                    
                      h
                      Nor
                      *
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Nor
                        
                      
                      )
                    
                  
                  ⁢
                  
                    
                      h
                      Nor
                      T
                    
                    ⁡
                    
                      (
                      
                        Z
                        
                          S
                          ⁢
                          _
                          ⁢
                          Nor
                        
                      
                      )
                    
                  
                
                )
              
            
            ⁢
            i
          
        
        .
      
    
  




For application to Multiple Input Multiple Output (MIMO) capacity optimization, as noted previously, the problem of mutual coupling for MIMO wireless systems has been considered previously by evaluating the capacity of the MIMO channel when the transmitting antennas are coupled and are driven by ideal current sources (i.e., a Norton source with infinite shunt impedance). Embodiments can consider the capacity of the MIMO channel with mutual coupling and general linear source models.
For a Thevenin source model, the MIMO channel can have M transmit antennas and N receive antennas. The M×N channel matrix can be denoted by H, where Hi,j can denote the channel observed at the j-th receive antenna when the precoder vi is applied at the transmitter, where




  
    v
    i
  
  =
  
    
      
        [
        
          
            
              0
              ⁢
              ...0
            
            
              }
              
                i
                -
                1
              
            
          
          ⁢
          10
          ⁢
          ...0
        
        ]
      
      T
    
    .
  




The capacity of this channel subject to a power constraint is known in the case that the transmitted power corresponding to the precoding vector v is given by




  
    
      
        
        v
        
      
      2
      2
    
    =
    
      
        
          v
          H
        
        ⁢
        v
      
      =
      
        
          ∑
          
            i
            =
            1
          
          M
        
        ⁢
        
          
            
            
              v
              i
            
            
          
          2
        
      
    
  
  ,




and also in the case that the transmitted power corresponding to the multi-layer precoding matrix V with L layers can be given by





  
    
      
      V
      
    
    2
    2
  
  =
  
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        
          V
          
            ·
            
              ,
              l
            
          
          H
        
        ⁢
        
          V
          
            ·
            
              ,
              l
            
          
        
      
    
    =
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        
          ∑
          
            l
            =
            1
          
          M
        
        ⁢
        
          
            
              
              
                V
                
                  i
                  ,
                  l
                
              
              
            
            2
          
          .
        
      
    
  




However, this condition is not satisfied in general if the transmitter antennas are coupled.
As discussed above, let Z denote the M×M impedance matrix for the transmit array, and let ZS_Thev denote the M×M diagonal source impedance. From above, the transmitted power for the precoding vector v can be given by




  
    
      
        v
        H
      
      ⁢
      
        Q
        Thev
      
      ⁢
      v
    
    =
    
      
        ∑
        
          i
          =
          1
        
        M
      
      ⁢
      
        
          ∑
          
            j
            =
            1
          
          M
        
        ⁢
        
          
            v
            i
            *
          
          ⁢
          
            Q
            
              
                Thev
                ⁢
                _
                ⁢
                i
              
              ,
              j
            
          
          ⁢
          
            v
            j
          
        
      
    
  
  ,




where

QThev=(ZS_Thev+Z)−HRe(Z)(ZS_Thev+Z)−1.

Thus, for a multi-layer precoding matrix, the transmitted power can be given by




  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        V
        
          ·
          
            ,
            l
          
        
        H
      
      ⁢
      
        Q
        Thev
      
      ⁢
      
        V
        
          ·
          
            ,
            l
          
        
      
    
  
  =
  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        ∑
        
          i
          =
          1
        
        M
      
      ⁢
      
        
          ∑
          
            j
            =
            1
          
          M
        
        ⁢
        
          
            V
            
              i
              ,
              l
            
            *
          
          ⁢
          
            Q
            
              
                Thev
                ⁢
                _
                ⁢
                i
              
              ,
              j
            
          
          ⁢
          
            V
            
              j
              ,
              l
            
          
        
      
    
  




This problem can be solved by transforming it to a similar problem for which the capacity maximization solution is known. As in previous sections, the matrix QThev can be factored as

QThev=PThevHPThev 

and the vector w can be defined such that

w=PThevv. 

If the precoder v is applied at the transmitter, the channel observed at the receiver can be given by

vTH=wTPThev−TH=wT(PThev−TH).

Thus, the precoding vector v that maximizes the capacity of the channel H subject the power constraint

vHQThevv≦P 

can be given by

v=PThev−1w, 

where w maximizes the capacity of the channel PThev−TH subject to the constraint

vHQThevv=wHw≦P. 

Similarly, for the multi-layer precoding matrix V, a second multi-layer precoding matrix W can be defined such that

W=PThevV. 

If the multi-layer precoding matrix V is applied at the transmitter, the vector channel observed at the receiver can be given by

VTH=WTPThev−TH=WT(PThev−TH).

Thus, the multi-layer precoding matrix V which maximizes the capacity of the channel H subject the power constraint





  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        V
        
          ·
          
            ,
            l
          
        
        H
      
      ⁢
      
        Q
        Thev
      
      ⁢
      
        V
        
          ·
          
            ,
            l
          
        
      
    
  
  ≤
  P




can be given by

V=PThev−1W, 

where the multi-layer precoding matrix W maximizes the capacity of the channel PThev−TH subject to the constraint





  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        W
        
          ·
          
            ,
            l
          
        
        H
      
      ⁢
      
        W
        
          ·
          
            ,
            l
          
        
      
    
  
  ≤
  
    P
    .
  




Maximizing capacity subject to a power constraint for the case above that the transmitted power corresponding to a multi-layer precoder can be given by




  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        W
        
          ·
          
            ,
            l
          
        
        H
      
      ⁢
      
        W
        
          ·
          
            ,
            l
          
        
      
    
  
  ,




the precoder selection and power allocation which maximizes capacity is known. In particular, for a given precoder W and M×1 input signal vector d, the N×1 received signal vector y can be given by

y=(PThev−TH)TW d+n, 

where the N×1 vector n can denote a zero-mean complex Gaussian random vector with covariance matrix given by

E(nnH)=σ2IN×N,

where IN×N can be the identity matrix of dimension N.

The channel matrix PThev−TH can be expressed in terms of its singular value decomposition as

PThev−TH=UΣXH,

where U can be a unitary M×M matrix, Σ can be an M×N rectangular diagonal matrix with non-negative real numbers on the diagonal, and XH can be a unitary N×N matrix where the superscript H denotes the conjugate transpose of the matrix. The columns of U can be the left singular vectors of PThev−TH, the columns of X can be the right singular vectors of PThev−TH, and the diagonal elements of Σ can be the singular values of PThev−TH. Let {α1,α2, . . . ,αK} denote the singular values of PThev−TH in order from top left to bottom right, where K can be defined as the minimum of the number of antenna elements at the transmitter and the receiver, so that

K=min(M,N).

Similarly, let {u1,u2, . . . ,uM} denote the columns of U in order from left to right, which can be the left singular vectors of PThev−TH, and let {x1,x2, . . . ,xN} denote the columns of X in order from left to right, which can be the right singular vectors of PThev−TH.
Now, the real positive value β can be chosen such that




  
    P
    =
    
      
        ∑
        
          i
          =
          1
        
        K
      
      ⁢
      
        max
        ⁡
        
          (
          
            
              β
              -
              
                
                  σ
                  2
                
                
                  
                    
                    
                      α
                      i
                    
                    
                  
                  2
                
              
            
            ,
            0
          
          )
        
      
    
  
  ,




where P can be the constraint on transmit power. Now, define J as the set of indices {j1,j2, . . . ,jL}, for which





  
    
      β
      -
      
        
          σ
          2
        
        
          
            
            
              α
              
                j
                i
              
            
            
          
          2
        
      
    
    >
    0
  
  ,




where 0<L≦K. Then the optimal precoder W for the channel PThev−TH can have rank L, where the L columns of W are the L left singular vectors of PThev−TH belonging to the set J, so that

W=[uj1 uj2 . . . ujL].

For this precoder, the input signal vector d can have dimension L×1, so that

d=[d1 d2 . . . dL]T 

and the power allocated to the symbol di can be given by





  β
  -
  
    
      σ
      2
    
    
      
        
        
          α
          
            j
            i
          
        
        
      
      2
    
  




for 1≦i≦L. The resulting channel capacity C for the indicated precoder and power allocation can be then given by





  
    
      
        C
        =
        
          
            ∑
            
              i
              =
              1
            
            L
          
          ⁢
          
            
              log
              2
            
            ⁡
            
              (
              
                1
                +
                
                  
                    
                      
                        
                        
                          α
                          
                            j
                            i
                          
                        
                        
                      
                      2
                    
                    
                      σ
                      2
                    
                  
                  ⁢
                  
                    (
                    
                      β
                      -
                      
                        
                          σ
                          2
                        
                        
                          
                            
                            
                              α
                              
                                j
                                i
                              
                            
                            
                          
                          2
                        
                      
                    
                    )
                  
                
              
              )
            
          
        
      
    
  
  
    
      
        =
        
          
            ∑
            
              i
              =
              1
            
            L
          
          ⁢
          
            
              
                log
                2
              
              ⁡
              
                (
                
                  
                    β
                    
                      σ
                      2
                    
                  
                  ⁢
                  
                    
                      
                      
                        α
                        
                          j
                          i
                        
                      
                      
                    
                    2
                  
                
                )
              
            
            .
          
        
      
    
  




To summarize, the precoder that maximizes the capacity of the channel H subject to a power constraint P can be given by

V=PThev−1W, 

where W can maximize the capacity of the channel PThev−TH subject to the same power constraint P. Furthermore, the L columns of W can be the L left singular vectors of PThev−TH belonging to the set J, so that

W=[uj1 uj2 . . . ujL],

where J can be the set of indices {j1, j2, . . . , jL}, for which





  
    
      β
      -
      
        
          σ
          2
        
        
          
            
            
              α
              
                j
                i
              
            
            
          
          2
        
      
    
    >
    0
  
  ,




and the real positive value β can be chosen such that





  P
  =
  
    
      ∑
      
        i
        =
        1
      
      K
    
    ⁢
    
      
        max
        ⁡
        
          (
          
            
              β
              -
              
                
                  σ
                  2
                
                
                  
                    
                    
                      α
                      i
                    
                    
                  
                  2
                
              
            
            ,
            0
          
          )
        
      
      .
    
  




Finally, the input signal vector d can have dimension L×1, so that

d=[d1 d2 . . . dL]T 

and the power allocated to the symbol di can be given by





  β
  -
  
    
      σ
      2
    
    
      
        
        
          α
          
            j
            i
          
        
        
      
      2
    
  




for 1≦i≦L, and the resulting channel capacity C for the indicated precoder and power allocation can then be given by





  
    
      
        C
        =
        
          
            ∑
            
              i
              =
              1
            
            L
          
          ⁢
          
            
              log
              2
            
            ⁡
            
              (
              
                1
                +
                
                  
                    
                      
                        
                        
                          α
                          
                            j
                            i
                          
                        
                        
                      
                      2
                    
                    
                      σ
                      2
                    
                  
                  ⁢
                  
                    (
                    
                      β
                      -
                      
                        
                          σ
                          2
                        
                        
                          
                            
                            
                              α
                              
                                j
                                i
                              
                            
                            
                          
                          2
                        
                      
                    
                    )
                  
                
              
              )
            
          
        
      
    
  
  
    
      
        =
        
          
            ∑
            
              i
              =
              1
            
            L
          
          ⁢
          
            
              
                log
                2
              
              ⁡
              
                (
                
                  
                    β
                    
                      σ
                      2
                    
                  
                  ⁢
                  
                    
                      
                      
                        α
                        
                          j
                          i
                        
                      
                      
                    
                    2
                  
                
                )
              
            
            .
          
        
      
    
  




From the discussion above, the solution for the Norton source model can be straightforward. If the precoder i is applied at the transmitter, the channel observed at the receiver can be given by iTH. For the Norton source with diagonal source impedance ZS_Nor, we have

QNor=ZS_NorH(ZS_Nor+Z)−HRe(Z)(ZS_Nor+Z)−1ZS_Nor 

and

QNor=PNorHPNor.

If we define

w=PNori 

then the channel observed at the receiver can be expressed as

iTH=wTPNor−TH=wT(PNor−TH).

Thus, the precoding vector i which maximizes the capacity of the channel H subject the power constraint

iHQNori≦P 

can be given by

i=PNor−1w, 

where w can maximize the capacity of the channel PNor−TH subject to the constraint

iHQNori=wHw≦P. 

Similarly, for the multi-layer precoding matrix I, a second multi-layer precoding matrix W can be defined such that

W=PNorI. 

If the multi-layer precoding matrix V is applied at the transmitter, the vector channel observed at the receiver can be given by

ITH=WTPNor−TH=WT(PNor−TH).

Thus, the multi-layer precoding matrix V that maximizes the capacity of the channel H subject the power constraint




  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        I
        
          ·
          
            ,
            l
          
        
        H
      
      ⁢
      
        Q
        Nor
      
      ⁢
      
        I
        
          ·
          
            ,
            l
          
        
      
    
  
  ≤
  P




can be given by

I=PNor−1W, 

where the multi-layer precoding matrix W can maximize the capacity of the channel PNor−TH subject to the constraint





  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      
        W
        
          ·
          
            ,
            l
          
        
        H
      
      ⁢
      
        W
        
          ·
          
            ,
            l
          
        
      
    
  
  ≤
  
    P
    .
  




FIG. 5 is an example flowchart 500 illustrating the operation of a transmitting device, such as the device 110 or the device 120, according to a possible embodiment. For example, the method of the flowchart 500 can be performed in a user equipment, in a base station, or in any other device that uses precoders and has a transmitter. At 510, the flowchart 500 can begin. At 520, channel measurements of a channel between a transmitting device and a receiving device can be ascertained. For example, the channel measurements can be made at the transmitting device or the channel measurements can be made at the receiving device and signaled back to the transmitting device.
At 530, a first channel matrix H can be transformed into a transformed second channel matrix P−TSourceH, where the term “−T” can indicate an inverse transpose. The first channel matrix H can be a channel from a transmitter array of antennas of the transmitting device to at least one receiver antenna of the receiving device. The transmitter array of antennas can be mutually coupled in that voltage or current applied to one antenna element can induce a voltage or current on another antenna element in the transmitter array of antennas. The first channel matrix H can be based on the channel measurements of the channel between the transmitting device and the receiving device. For example, the first channel matrix H can be based on channel reciprocity, which can be based on the transmitting device taking measurements of reference symbols transmitted by a target receiving device. The first channel matrix H can also be based on channel measurements taken at the target receiving device and signaled back to the transmitting device.
A transformation matrix P−TSource for the transformation of the second channel P−TSourceH can be an inverse transpose of a right factor of a factorization of a Hermitian and non-negative definite matrix. The transformed second channel matrix P−TSourceH can be a product of the transformation matrix P−TSource and the first channel matrix H. The Hermitian and non-negative definite matrix can be a function of a source model of the transmitter, a source impedance of a transmitter of the transmitting device, and an impedance matrix of the transmitter array antennas. PSource can be based on

QThev=PSourceHPSource,

where

QThev=(ZS_Thev+Z)−HRe(Z)(ZS_Thev+Z)−1,

and

where ZS_Thev can be a diagonal matrix of a transmitter source impedances of the transmitting device and Z can be an impedance matrix of the transmitter array antennas. PSource can also be based on

QNor=PSourceHPSource,

where

QNor=ZS_NorH(ZS_Nor=Z)−HRe(Z)(ZS_Nor+Z)−1ZS—Nor,

and where ZS_Nor can be a diagonal matrix of a transmitter source impedances of the transmitting device and Z can be an impedance matrix of the transmitter array antennas.

At 540, a precoding matrix W can be determined that maximizes the capacity for the transformed second channel matrix P−TSourceH subject to a power constraint of the precoding matrix W. The precoding matrix can be a one dimensional matrix, such as a vector for a single receive antenna, or can be a multidimensional matrix, such as a two -dimensional matrix for multiple receive antennas. Similarly the channel matrix can be a one dimensional matrix or a multidimensional matrix. Each column vector of the precoding matrix W can be a left singular vector of the transformed second channel matrix P−TSourceH. A power constraint P on the precoding matrix W can be expressed as




  
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        
          W
          
            ·
            
              ,
              l
            
          
          H
        
        ⁢
        
          W
          
            ·
            
              ,
              l
            
          
        
      
    
    ≤
    P
  
  ,




where P can denote the power constraint which can be a real positive number and L can be the number of columns in W and can also be the number of transmission layers. P−TSource can be P−TThev, P−TNorton, and/or can be based on other models of an antenna array.

At 550, the channel capacity maximizing precoding matrix W for the transformed second channel matrix can be converted into an optimal precoding matrix V for the first channel matrix. An optimal beamformer vopt in the optimal precoding matrix V can be given by vopt=PThev−1eThev_max, where eThev_max can be the eigenvector corresponding to the largest eigenvalue, λThev_max, of the matrix (PThev−1)H(h*Thev(ZS_Thev)hThevT(ZS_Thev))PThev−1. An optimal beamformer iopt in the optimal precoding matrix V can also be given by iopt=PNor−1eNor_max, where eNor_max can be the eigenvector corresponding to the largest eigenvalue, λNor_max, of the matrix (PNor−1)H(h*Nor(ZS_Nor)hNorT(ZS_Nor))PNor−1.
At 560, a signal can be received for transmission. The signal can be a vector signal. At 570, the optimal precoding matrix V can be applied to the signal to generate a precoded signal for transmission over a physical channel. Columns of the optimal precoding matrix V can include a plurality of precoding vectors. Applying can include multiplying the vector signal by the optimal precoding matrix V to generate a precoded signal vector for transmission over a physical channel. At 580, the precoded signal can be transmitted. At 560, the flowchart 500 can end.
It should be understood that, notwithstanding the particular steps as shown in the figures, a variety of additional or different steps can be performed depending upon the embodiment, and one or more of the particular steps can be rearranged, repeated or eliminated entirely depending upon the embodiment. Also, some of the steps performed can be repeated on an ongoing or continuous basis simultaneously while other steps are performed. Furthermore, different steps can be performed by different elements or in a single element of the disclosed embodiments.
FIG. 6 is an example block diagram of an apparatus 600, such as the device 110 or the device 120, according to a possible embodiment. The apparatus 600 can include a housing 610, a controller 620 within the housing 610, audio input and output circuitry 630 coupled to the controller 620, a display 640 coupled to the controller 620, a transceiver 650 coupled to the controller 620, a plurality of antennas 655 and 657, such as an array of antennas, coupled to the transceiver 650, a user interface 660 coupled to the controller 620, a memory 670 coupled to the controller 620, and a network interface 680 coupled to the controller 620. The apparatus 600 can also include additional elements or less elements depending on the device in which it is implemented. The apparatus 600 can perform the methods described in all the embodiments.
The display 640 can be a viewfinder, a liquid crystal display (LCD), a light emitting diode (LED) display, a plasma display, a projection display, a touch screen, or any other device that displays information. The transceiver 650 can include a transmitter and/or a receiver. The audio input and output circuitry 630 can include a microphone, a speaker, a transducer, or any other audio input and output circuitry. The user interface 660 can include a keypad, a keyboard, buttons, a touch pad, a joystick, a touch screen display, another additional display, or any other device useful for providing an interface between a user and an electronic device. The network interface 680 can be a Universal Serial Bus (USB) port, an Ethernet port, an infrared transmitter/receiver, an IEEE 1396 port, a WLAN transceiver, or any other interface that can connect an apparatus to a network, device, or computer and that can transmit and receive data communication signals. The memory 670 can include a random access memory, a read only memory, an optical memory, a flash memory, a removable memory, a hard drive, a cache, or any other memory that can be coupled to a wireless communication device. The plurality of antennas 655 and 657 can be considered a transmitter array of antennas when the transceiver 650 is transmitting signals. The transmitter array of antennas 655 and 657 can include two or more antennas. The transmitter array of antennas 655 and 657 can be mutually coupled in that one of a voltage and current applied to one antenna element induces a voltage or current on another antenna element in the transmitter array of antennas 655 and 657.
The apparatus 600 or the controller 620 may implement any operating system, such as Microsoft Windows®, UNIX®, or LINUX®, Android®, or any other operating system. Apparatus operation software may be written in any programming language, such as C, C++, Java or Visual Basic, for example. Apparatus software may also run on an application framework, such as, for example, a Java® framework, a .NET® framework, or any other application framework. The software and/or the operating system may be stored in the memory 670 or elsewhere on the apparatus 600. The apparatus 600 or the controller 620 may also use hardware to implement disclosed operations. For example, the controller 620 may be any programmable processor. Disclosed embodiments may also be implemented on a general-purpose or a special purpose computer, a programmed microprocessor or microprocessor, peripheral integrated circuit elements, an application-specific integrated circuit or other integrated circuits, hardware/electronic logic circuits, such as a discrete element circuit, a programmable logic device, such as a programmable logic array, field programmable gate-array, or the like. In general, the controller 620 may be any controller or processor device or devices capable of operating a communication device and implementing the disclosed embodiments.
In operation, the controller 620 can be configured to ascertain channel measurements of a channel between the apparatus 600 and a receiving device. For example, the channel measurements can be ascertained based on channel reciprocity, which can be based on the apparatus 600 taking measurements of reference symbols transmitted by a target receiving device. The channel measurements can also be ascertained based on channel measurements taken at a target receiving device and signaled back to the apparatus 600.
The controller 620 can transform a first channel matrix H into a transformed second channel matrix P−TsourceH. The first channel matrix H can be a channel from the transmitter array of antennas 655 and 657 to at least one receiver antenna of a receiving device. For example, the first channel matrix H can be based on the channel measurements of the channel between the transceiver 650 and the receiving device. A transformation matrix P−TSource for the transformation of the second channel matrix P−TSourceH can be an inverse transpose of a right factor of a factorization of a Hermitian and non-negative definite matrix. The transformed second channel matrix P−TSourceH can be a product of the transformation matrix P−TSource and the first channel matrix H. The Hermitian and non-negative definite matrix can be a function of a source model of the transmitter of the transceiver 650, a source impedance of the transmitter, and an impedance matrix of the transmitter array antennas 655 and 657. PSource can be based on

QThev=PSourceHPSource,

where

QThev=(ZS_Thev+Z)−HRe(Z)(ZS_Thev+Z)−1,

and where ZS_Thev can be a diagonal matrix of a transmitter source impedances of the transmitting device and Z can be an impedance matrix of the transmitter array antennas. PSource can also be based on

QNor=PSourceHPSource,

where

QNor=ZS_NorH(ZS_Nor+Z)−HRe(Z)(ZS_Nor+Z)−1ZS_Nor,

and where ZS_Nor can be a diagonal matrix of a transmitter source impedances of the transmitting device and Z can be an impedance matrix of the transmitter array antennas.

The controller 620 can determine a precoding matrix W that maximizes the capacity for the transformed second channel matrix P−TSourceH subject to a power constraint of the precoding matrix W. Each column vector of the precoding matrix W can be a left singular vector of the transformed second channel matrix P−TSourceH. The controller 620 can convert the channel capacity maximizing precoding matrix W for the transformed second channel matrix into an optimal precoding matrix V for the first channel matrix. Columns of the optimal precoding matrix V can include a plurality of precoding vectors.
The controller 620 can receive a signal for transmission. The signal can be a vector signal. The controller 620 can apply the optimal precoding matrix V to the signal to generate a precoded signal for transmission over a physical channel. The controller 620 can apply the optimal precoding matrix V by multiplying the vector signal by the optimal precoding matrix V to generate a precoded signal vector for transmission over a physical channel. The transceiver 650 can transmit the precoded signal over the physical channel via the transmitter array of antennas.
Some example embodiments above describe a two-port model for a two-element antenna array, and more generally, an M-port model for an M-element antenna array. In these example embodiments, impedance parameters (Z matrix) model the relationship of the voltages and the currents for this two-port, or M-port, antenna array as

v=Zi. 

There are also other equivalent parameters that can be used to model the antenna array. For example, other sets of parameters can include admittance parameters (Y), hybrid parameters (H), inverse hybrid parameters (G), ABCD parameters (ABCD), scattering parameters (S), scattering transfer parameters (T), and other parameters useful for modeling an antenna array. All of these models are equivalent, even if they look slightly different. For example, the admittance parameters (Y) can have the voltage -current relationship of

i=Yv 

so that

v=Y−1i 

and thus

Z=Y−1.

The last expression can give the relationship between the impedance parameters and the admittance parameters. Thus, if the antenna array is represented in terms of its admittance parameters, the precoder transformation may look slightly different, but is still equivalent to the transformation with the impedance parameters. For example, the transformation can be exactly the same with the exact same or similar mapping of a precoder to a transformed precoder, except that Z can be replaced everywhere by Y−1, and these can be exactly equal. Similarly, all the other parameters cases above can be converted to Z parameters and are thus equivalent.
The method of this disclosure can be implemented on a programmed processor. However, the controllers, flowcharts, and modules may also be implemented on a general purpose or special purpose computer, a programmed microprocessor or microcontroller and peripheral integrated circuit elements, an integrated circuit, a hardware electronic or logic circuit such as a discrete element circuit, a programmable logic device, or the like. In general, any device on which resides a finite state machine capable of implementing the flowcharts shown in the figures may be used to implement the processor functions of this disclosure.
While this disclosure has been described with specific embodiments thereof, it is evident that many alternatives, modifications, and variations will be apparent to those skilled in the art. For example, various components of the embodiments may be interchanged, added, or substituted in the other embodiments. Also, all of the elements of each figure are not necessary for operation of the disclosed embodiments. For example, one of ordinary skill in the art of the disclosed embodiments would be enabled to make and use the teachings of the disclosure by simply employing the elements of the independent claims. Accordingly, embodiments of the disclosure as set forth herein are intended to be illustrative, not limiting. Various changes may be made without departing from the spirit and scope of the disclosure.
In this document, relational terms such as “first,” “second,” and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. The phrase “at least one of” followed by a list is defined to mean one, some, or all, but not necessarily all of, the elements in the list. The terms “comprises,” “comprising,” or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. An element proceeded by “a,” “an,” or the like does not, without more constraints, preclude the existence of additional identical elements in the process, method, article, or apparatus that comprises the element. Also, the term “another” is defined as at least a second or more. The terms “including,” “having,” and the like, as used herein, are defined as “comprising.” Furthermore, the background section is written as the inventor's own understanding of the context of some embodiments at the time of filing and includes the inventor's own recognition of any problems with existing technologies and/or problems experienced in the inventor's own work.

## Claims

1. A method in a transmitting device, the method comprising:
transforming a first channel matrix H into a transformed second channel matrix P−TSourceH, where the first channel matrix H is a channel from a transmitter array of antennas of the transmitting device to at least one receiver antenna of a receiving device, and where a transformation matrix P−TSource is a function of a transmitter source impedance and an impedance matrix of the transmitter array of antennas;
determining a precoding matrix W that maximizes the capacity for the transformed second channel matrix P−TSourceH subject to a power constraint of the precoding matrix W;
converting the channel capacity maximizing precoding matrix W for the transformed second channel matrix into an optimal precoding matrix V for the first channel matrix;
receiving a signal for transmission;
applying the optimal precoding matrix V to the signal to generate a precoded signal for transmission over a physical channel; and
transmitting the precoded signal.

2. The method according to claim 1, wherein each column vector of the precoding matrix W is a left singular vector of the transformed second channel matrix P−TSourceH.

3. The method according to claim 1, wherein the transformation matrix P−TSource for the transformed second channel matrix P−TSourceH is an inverse transpose of a right factor of a factorization of a Hermitian and non-negative definite matrix.

4. The method according to claim 3, wherein the transformed second channel matrix is a product of the transformation matrix P−TSource and the first channel matrix H.

5. The method according to claim 3, wherein the Hermitian and non -negative definite matrix is a function of a source model of the transmitter, the transmitter source impedance of the transmitting device, and the impedance matrix of the transmitter array of antennas.

6. The method according to claim 3, wherein PSource is based on


QThev=PSourceHPSource,

where PSourceH is a left factor of a factorization of a Hermitian matrix QThev and PSource is a right factor of a factorization of the Hermitian matrix QThev,
where


QThev=(ZS_Thev+Z)−HRe(Z)(ZS_Thev+Z)−1, and


where ZS_Thev is a diagonal matrix of the transmitter source impedance of the transmitting device and Z is the impedance matrix of the transmitter array of antennas.

7. The method according to claim 3, wherein PSource is based on


QNorPSourceHPSource,

where PSourceH is a left factor of a factorization of a Hermitian matrix QNor and PSource is a right factor of a factorization of the Hermitian matrix QNor,
where


QNor=ZS_NorH(ZS_Nor+Z)−HRe(Z)(ZS_Nor+Z)−1ZS_Nor, and


where ZS_Nor is a diagonal matrix of a the transmitter source impedance of the transmitting device and Z is the impedance matrix of the transmitter array of antennas.

8. The method according to claim 1,
wherein the signal comprises a vector signal,
wherein columns of the optimal precoding matrix V include a plurality of precoding vectors, and
wherein applying comprises multiplying the vector signal by the optimal precoding matrix V to generate a precoded signal vector for transmission over a physical channel.

9. The method according to claim 1, wherein the transmitter array of antennas are mutually coupled in that voltage or current applied to one antenna element induces a voltage or current on another antenna element in the transmitter array of antennas.

10. The method according to claim 1, further comprising ascertaining channel measurements of a channel between the transmitting device and the receiving device,
wherein the first channel matrix H is based on the channel measurements of the channel between the transmitting device and the receiving device.

11. An apparatus comprising:
a transmitter array of antennas;
a transceiver coupled to the transmitter array of antennas; and
a controller coupled to the transceiver, the controller configured to
transform a first channel matrix H into a transformed second channel matrix P−TSourceH, where the first channel matrix H is a channel from the transmitter array of antennas to at least one receiver antenna of a receiving device, and where a transformation matrix P−TSource is a function of a transmitter source impedance and an impedance matrix of the transmitter array of antennas,
determine a precoding matrix W that maximizes the capacity for the transformed second channel matrix P−TSourceH subject to a power constraint of the precoding matrix W,
convert the channel capacity maximizing precoding matrix W for the transformed second channel matrix into an optimal precoding matrix V for the first channel matrix,
receive a signal for transmission, and
apply the optimal precoding matrix V to the signal to generate a precoded signal for transmission over a physical channel,

wherein the transceiver is configured to transmit the precoded signal over the physical channel via the transmitter array of antennas.

12. The apparatus according to claim 11, wherein each column vector of the precoding matrix W is a left singular vector of the transformed second channel matrix P−TSourceH.

13. The apparatus according to claim 11, wherein the transformation matrix P−TSource for the transformed second channel matrix P−TSourceH is an inverse transpose of a right factor of a factorization of a Hermitian and non -negative definite matrix.

14. The apparatus according to claim 13, wherein the transformed second channel matrix P−TSourceH is a product of the transformation matrix P−TSource and the first channel matrix H.

15. The apparatus according to claim 13, wherein the Hermitian and non -negative definite matrix is a function of a source model of a transmitter of the transceiver, the transmitter source impedance, and the impedance matrix of the transmitter array of antennas.

16. The apparatus according to claim 13, wherein PSource is based on


QThev=PSourceHPSource,

where PSourceH is a left factor of a factorization of a Hermitian matrix QThev and PSource is a right factor of a factorization of the Hermitian matrix QThev,
where


QThev=(ZS_Thev+Z)−HRe(Z)(ZS_Thev+Z)−1, and


where ZS_Thev is a diagonal matrix of the transmitter source impedance of a transmitter of the transceiver and Z is the impedance matrix of the transmitter array of antennas.

17. The apparatus according to claim 13, wherein PSource is based on


QNor=PSourceHPSource,

where PSourceH is a left factor of a factorization of a Hermitian matrix QNor and PSource is a right factor of a factorization of the Hermitian matrix QNor ,
where


QNor=ZS_NorH(ZS_Nor+Z)−HRe(Z)(ZS_Nor+Z)−1ZS—Nor, and


where ZS_Nor is a diagonal matrix of the transmitter source impedance of transmitter of the transceiver and Z is the impedance matrix of the transmitter array of antennas.

18. The apparatus according to claim 11,
wherein the signal comprises a vector signal,
wherein columns of the optimal precoding matrix V include a plurality of precoding vectors, and
wherein the controller is configured to apply the optimal precoding matrix V by multiplying the vector signal by the optimal precoding matrix V to generate a precoded signal vector for transmission over a physical channel.

19. The apparatus according to claim 11, wherein the transmitter array of antennas are mutually coupled in that one of a voltage and current applied to one antenna element induces a voltage or current on another antenna element in the transmitter array of antennas.

20. The method according to claim 11,
wherein the controller is configured to ascertain channel measurements of a channel between the transceiver and the receiving device, and
wherein the first channel matrix H is based on the channel measurements of the channel between the apparatus and the receiving device.

21. A method in a transmitting device, the method comprising:
ascertaining channel measurements of a channel between the transmitting device and a receiving device;
transforming a first channel matrix into a transformed second channel matrix, where the first channel matrix is a channel from a transmitter array of antennas of the transmitting device to at least one receiver antenna of the receiving device, where the first channel matrix is based on the channel measurements of the channel between the transmitting device and the receiving device, where the transmitter array of antennas are mutually coupled in that voltage or current applied to one antenna element induces a voltage or current on another antenna element in the transmitter array of antennas, and where a transformation matrix for the transforming is a function of a transmitter source impedance and an impedance matrix of the transmitter array of antennas;
determining a first precoding matrix that maximizes the capacity for the transformed second channel matrix subject to a power constraint of the first precoding matrix;
converting the channel capacity maximizing first precoding matrix for the transformed second channel matrix into a second precoding matrix for the first channel matrix;
receiving a signal for transmission;
applying the second precoding matrix to the signal to generate a precoded signal for transmission over a physical channel; and
transmitting the precoded signal.