# Patent

## Title
System and method for enhancing throughput in an additive gaussian noise channel with a predetermined rate set and unknown interference

## Abstract

A system for, and method of, determining channel throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates and a related system for, and method of, maximizing channel throughput. In one embodiment, the system includes: (1) a noise level calculator configured to calculate associated optimal noise levels for at least some rate subsets of the set of allowed transmission rates and (2) a code feasibility calculator associated with the noise level calculator and configured to determine whether feasible superposition codes exist for the at least some rate subsets given the associated optimal noise levels and thereby identify ones of the at least some rate subsets that are suitable for subsequent wireless transmissions.

## Background
TECHNICAL FIELD OF THE INVENTION
The present invention is directed, in general, to wireless networks and, more specifically, to a system and method for determining throughput for an additive white Gaussian noise (AWGN) channel with a predetermined rate set and unknown interference.
BACKGROUND OF THE INVENTION
Over the last several years, wireless networks have become indispensable tools for carrying voice and data communications for both business and personal use. Wireless networks partition a given geographical region into cells. Each cell is given a base station that serves as the anchor point to the network for all users in the cell. Users are assigned to the different cells based on, for example, channel conditions and network loading.
In current networks, coordination between base stations is not considered; users are only controlled by their assigned base station. However, due to the nature of wireless transmission, users in neighboring cells interfere with the transmission in a given cell. In general, this out-of-cell interference is unknown and uncontrollable. Unfortunately, the performance of the channels in the cell critically depends on the out-of-cell interference since the achievable transmission rate of each channel is related to its signal-to-interference-plus-noise ratio experienced at the receiver.
The network is then faced with the task of determining its maximum average rate of reliable transmission under the uncertainty on the interference level. Cover, in “Broadcast Channels,” IEEE Transactions on Information Theory, vol. IT-18, pp. 2-14, January 1972 (incorporated herein by reference), defined such compound channels and showed their equivalence to broadcast channels. Cover proposed and analyzed a solution, which he called “superposition coding,” for reliable transmission over such channels. In particular, Cover showed that superposition coding achieves larger throughput than time-sharing techniques.
Subsequently, Bergmans, in “A Simple Converse for Broadcast Channels with Additive White Gaussian Noise Maximum Likelihood Decoding,” IEEE Transactions on Information Theory, vol. IT-20, pp. 279-280, March 1974 (incorporated herein by reference), demonstrated the optimality of superposition coding. The main idea behind superposition coding, as it pertains to channel throughput, is to make sure that some data is correctly received even when the interference is large; additional information can be received when the interference is small.
Shamai, in “A Broadcast Strategy for the Gaussian Slowly Fading Channel,” in Proceedings of the IEEE International Symposium on Information Theory 1997, p. 150 (incorporated herein by reference), calculates the maximum throughput achievable under superposition coding, but shows that achieving the maximum throughput requires an infinite number of coding levels. Liu, et al., in “Optimal Rate Allocation for Superposition Coding in Quasi-static Fading Channels,” in Proceedings of IEEE International Symposium on Information Theory 2002, p. 111 (incorporated herein by reference), consider only a finite number of levels in the superposition coding and show that the performance achieved with two levels is fairly close to the maximum throughput. However, Liu, et al., assume that the transmitter does not have any information on the state of the channel. On the other hand, Whiting, et al., in “Optimal Encoding over Uncertain Channels with Decoding Delay Constraints,” in Proceedings of the IEEE International Symposium on Information Theory 2000, p. 430 (incorporated herein by reference), assume the transmitter to have delayed feedback knowledge about the channel state and a finite decoding delay. Whiting, et al., then provide the maximum achievable average rate using a broadcast approach.
All of the above-referenced prior art approaches have assumed that the transmission rates are of infinite granularity. However, in practical scenarios (due to standard specifications, practical implementations and the cost of building commercial transmitters and receivers), the set of allowed transmission rates is finite and predetermined, and the transmitters are constrained to only use some of these specified rates.
What is needed in the art is a technique for determining the throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates. What is further needed in the art is a system and a method for performing such technique.
SUMMARY OF THE INVENTION
To address the above-discussed deficiencies of the prior art, the present invention provides, in one aspect, a system for determining channel throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates. In one embodiment, the system includes: (1) a noise level calculator configured to calculate associated optimal noise levels for at least some rate subsets of the set of allowed transmission rates and (2) a code feasibility calculator associated with the noise level calculator and configured to determine whether feasible superposition codes exist for the at least some rate subsets given the associated optimal noise levels and thereby identify ones of the at least some rate subsets that are suitable for subsequent wireless transmissions.
In another aspect, the present invention provides a method of determining channel throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates. In one embodiment, the method includes: (1) calculating associated optimal noise levels for at least some rate subsets of the set of allowed transmission rates and (2) determining whether feasible superposition codes exist for the at least some rate subsets given the associated optimal noise levels and thereby identifying ones of the at least some rate subsets that are suitable for subsequent wireless transmissions.
It is assumed that the channel can reliably support at least one of the rates in the pre-determined rate set with an arbitrarily small probability of decoding error. Assuming that the rates need to be chosen from a finite set of predetermined rates leads to a combinatorial optimization problem and significantly contributes to the complexity of the problem. A sequence of propositions culminating in an application of the Kuhn-Tucker theorem to solve a constrained optimization problem will be set forth. This yields an algorithm that optimizes both the subset of transmission rates that ought to be used, and the power allocations to the various rates and the corresponding noise power levels.
The foregoing has outlined, rather broadly, preferred and alternative features of the present invention so that those skilled in the art may better understand the detailed description of the invention that follows. Additional features of the invention will be described hereinafter that form the subject of the claims of the invention. Those skilled in the art should appreciate that they can readily use the disclosed conception and specific embodiment as a basis for designing or modifying other structures for carrying out the same purposes of the present invention. Those skilled in the art should also realize that such equivalent constructions do not depart from the spirit and scope of the invention in its broadest form.



BRIEF DESCRIPTION OF THE DRAWINGS
For a more complete understanding of the present invention, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which:
FIG. 1 illustrates a block diagram of one embodiment of a system for determining the throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates constructed according to the principles of the present invention;
FIG. 2 illustrates a flow diagram of one embodiment of a method of determining the throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates carried out according to the principles of the present invention;
FIG. 3 illustrates a graph of signal-to-noise ratio (SNR) versus optimal throughput using superposition coding for different values of mean noise power 1/α, when the pdf of the noise power is exponential;
FIG. 4 illustrates a graph of SNR versus optimal throughput using superposition coding, both when the optimal rate subset for each SNR is unrestricted, and when it is forced to be of size 1 (single-rate coding), for α=0.7; and
FIG. 5 illustrates a graph of SNR versus the sum of the rates in the optimal rate subset for superposition and single-rate coding, for α=0.7.



DETAILED DESCRIPTION
A description of the channel model used herein will now be described, along with the assumptions made and a mathematical formulation of the problem of interest. A single user, point-to-point channel, modeled as a discrete-time, flat fading, additive white Gaussian noise channel with noise variance σ2 will be considered. The interference from other users' transmissions is assumed to be time-varying and modeled as a zero-mean, Gaussian random variable (r.v.), which is independent of the thermal noise of the channel. The interference power I, i.e., the variance of this Gaussian r.v., is itself an r.v. chosen according to a known probability density function which depends on the location and the transmission characteristics of the users. The joint effect of the interference and thermal noise will be considered as an equivalent noise process, the power of which is given as N=I+σ2, since the interference is assumed to be independent of the thermal noise. A new noise power N is then chosen according to a known probability density function (pdf) fN(•). This pdf is assumed to be monotone decreasing and have infinite support, taken to be [0, ∝), without loss of generality. It is implicitly assumed that interference cancellation between the different transmitting users is not employed at the receiver. In addition, it is assumed that the received power is upper-bounded by P. Perfect power control is employed to ensure that the received power remains constant and independent of channel condition and the level of interference. This requires that the transmitter has perfect and instantaneous knowledge of the channel gain between the transmitter and the receiver. Finally, it is assumed that a set of M possible rates exists from which the transmitter must choose its transmission rate(s).
The rate set is denoted by R. The elements of an arbitrary subset RK containing K of these M rates is denoted by RK≡{R1, R2, . . . , RK}. Reuse of a rate in a subset is not permitted, i.e., all elements of RK must be distinct.
Having set forth the applicable channel model, superposition coding will now be described. Given a rate subset RK of size K, K noise power levels exist, denoted (N1, . . . , NK). Without loss of generality, it is assumed that they are labeled in increasing order:

0≦N1≦N2≦ . . . ≦NK.  (1)

The codewords with rate Rk are allocated power Pk, k=1, . . . , K. A superposition decoding scheme operates as follows: a permutation π of the set {1, . . . , K}→{1, . . . , K} is chosen, and an attempt is made to decode the codeword with rate Rπ(K) at the highest noise level NK. If decoding is successful, the decoded codeword is subtracted from the received codeword. Then, an attempt is made to decode the codeword with rate Rπ(K−1) at the next highest noise level NK−1 and that decoded codeword subtracted. This proceeds until all codewords are successfully decoded. In other words, each rate Rj has a power allocation Pj, such that at noise power level Nk, the codewords of rate Rπ(k) can be decoded, assuming that the codewords with rates Rπ(k+1), . . . , Rπ(K) are already decoded and those with rates Rπ(1), . . . , Rπ(k−1) are not yet decoded:




  
    
      
        
          
            R
            
              π
              ⁡
              
                (
                k
                )
              
            
          
          =
          
            
              1
              2
            
            ⁢
            
              log
              (
              
                1
                +
                
                  
                    P
                    
                      π
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                  
                  
                    
                      
                        ∑
                        
                          i
                          =
                          1
                        
                        
                          k
                          -
                          1
                        
                      
                      ⁢
                      
                        P
                        
                          π
                          ⁡
                          
                            (
                            i
                            )
                          
                        
                      
                    
                    +
                    
                      N
                      k
                    
                  
                
              
              )
            
          
        
        ,
        
          k
          =
          1
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          K
          .
        
      
    
    
      
        (
        2
        )
      
    
  




Note that equation (2) implicitly assumes that Shannon's expression for the capacity of an AWGN channel sets the maximum transmission rate. The natural logarithm is used, so the rate is expressed as nats/s/Hz (“nats” being natural, or ln2, bits). Equation (2) associates the rate Rπ(k) with the noise power level Nk, k=1, . . . , K. Note that equation (2) is equivalent to





  
    
      
        
          
            P
            
              π
              ⁡
              
                (
                k
                )
              
            
          
          =
          
            
              γ
              
                π
                ⁡
                
                  (
                  k
                  )
                
              
            
            ⁡
            
              [
              
                
                  
                    ∑
                    
                      i
                      =
                      1
                    
                    
                      k
                      -
                      1
                    
                  
                  ⁢
                  
                    P
                    
                      π
                      ⁡
                      
                        (
                        i
                        )
                      
                    
                  
                
                +
                
                  N
                  k
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




γk=exp(2Rk)−1, k=1, . . . , K.  (4)

Thus, given π, R1, . . . , RK and N1, . . . , NK, equation (3) completely determines P1, . . . , PK. Further, it is clear from equation (3) that Pk≧0, k=1, . . . , K. Of course, the total power allocation cannot exceed the available total power P:




  
    
      
        
          P
          ≥
          
            
              P
              1
            
            +
            ...
            +
            
              
                P
                K
              
              ⁢
              
                
                  
                    ∑
                    
                      k
                      -
                      1
                    
                  
                  K
                
                ⁢
                
                  
                    q
                    
                      π
                      ⁡
                      
                        (
                        k
                        )
                      
                    
                  
                  ⁢
                  
                    N
                    k
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        5
        )
      
    
  




where equation (11) gives qπ(k), k=1, . . . , K.

Formally, for a given total received power P, rate set RK={R1, . . . , RK} and permutation π, if the noise power levels N=(N1, . . . , NK) satisfy equations (1) and (5), (RK, π, N, P) is called a “superposition coding scheme supporting the rate set RK” with decoding order π, noise power levels N, individual power allocations defined by equation (3), and total (received) power at most P.
The throughput using the superposition coding scheme (RK, π, N, P) is




  
    
      
        
          
            T
            ⁡
            
              (
              
                
                  R
                  K
                
                ,
                π
                ,
                N
              
              )
            
          
          =
          
            
              ∑
              
                k
                =
                1
              
              K
            
            ⁢
            
              
                R
                
                  π
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
              ⁢
              
                
                  ∫
                  0
                  
                    N
                    k
                  
                
                ⁢
                
                  
                    
                      f
                      N
                    
                    ⁡
                    
                      (
                      u
                      )
                    
                  
                  ⁢
                  
                    d
                    u
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        6
        )
      
    
  




where fN(•) is the pdf of the noise power. Given a rate set R of size M and total received power P, the goal is to find the superposition coding scheme (RK, π, N, P) that maximizes the throughput given in equation (6). In other words,





  
    
      
        
          max
          
            
              K
              ,
              
                R
                1
              
              ,
              ...
              ,
              
                R
                K
              
              ,
              π
              ,
              
                N
                1
              
              ,
              ...
              ,
              
                N
                K
              
            
            
              
                
                  1
                  ≤
                  K
                  ≤
                  M
                
                
                  
                    (
                    
                      
                        R
                        1
                      
                      ,
                      ...
                      ,
                      
                        R
                        K
                      
                    
                    )
                  
                  ⊆
                  R
                
              
              
                
                  0
                  ≤
                  
                    N
                    1
                  
                  ≤
                  ...
                  ≤
                  
                    N
                    K
                  
                
                
                  
                    
                      ∑
                      
                        k
                        =
                        1
                      
                      k
                    
                    ⁢
                    
                      
                        q
                        
                          π
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                      
                      ⁢
                      
                        N
                        k
                      
                    
                  
                  ≤
                  P
                
              
            
          
        
        ⁢
        
          
            ∑
            
              k
              =
              1
            
            K
          
          ⁢
          
            
              R
              
                π
                ⁡
                
                  (
                  k
                  )
                
              
            
            ⁢
            
              
                ∫
                0
                
                  N
                  k
                
              
              ⁢
              
                
                  
                    f
                    N
                  
                  ⁡
                  
                    (
                    u
                    )
                  
                
                ⁢
                
                  d
                  u
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




which is equivalent to





  
    
      
        
          
            max
            
              K
              ∈
              
                {
                
                  1
                  ,
                  
                    
                  
                  ⁢
                  ...
                  ⁢
                  
                    
                  
                  ,
                  M
                
                }
              
            
          
          ⁢
          
            
              max
              
                
                  {
                  
                    
                      R
                      1
                    
                    ,
                    
                      
                    
                    ⁢
                    ...
                    ⁢
                    
                      
                    
                    ,
                    
                      R
                      K
                    
                  
                  }
                
                ⊆
                R
              
            
            ⁢
            
              
                max
                
                  
                    π
                    :
                    
                      {
                      
                        1
                        ,
                        
                          
                        
                        ⁢
                        ...
                        ⁢
                        
                          
                        
                        ,
                        K
                      
                      }
                    
                  
                  ->
                  
                    {
                    
                      1
                      ,
                      
                        
                      
                      ⁢
                      ...
                      ⁢
                      
                        
                      
                      ,
                      K
                    
                    }
                  
                
              
              ⁢
              
                
                  
                    max
                    
                      
                        N
                        1
                      
                      ,
                      
                        
                      
                      ⁢
                      ...
                      ⁢
                      
                        
                      
                      ,
                      
                        N
                        K
                      
                    
                  
                  
                    
                      0
                      ≤
                      
                        N
                        1
                      
                      ≤
                      
                        
                      
                      ⁢
                      ...
                      ⁢
                      
                        
                      
                      ≤
                      
                        N
                        K
                      
                    
                    
                      
                        
                          ∑
                          
                            k
                            =
                            1
                          
                          K
                        
                        ⁢
                        
                          
                            q
                            
                              π
                              ⁡
                              
                                (
                                k
                                )
                              
                            
                          
                          ⁢
                          
                            N
                            K
                          
                        
                      
                      ≤
                      P
                    
                  
                
                ⁢
                
                  T
                  (
                  
                    
                  
                  ⁢
                  
                    
                      R
                      1
                    
                    ,
                    
                      
                    
                    ⁢
                    ...
                    ,
                    
                      R
                      K
                    
                    ,
                    π
                    ,
                    
                      N
                      1
                    
                    ,
                    ...
                    ,
                    
                      N
                      K
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        8
        )
      
    
  




where equation (6) gives T. The first three optimizations in equation (8) are computed by searching over all 2M−1 nonempty subsets of the given rate set R, and all K! permutations on each subset of size K, K=1, . . . , M. The innermost optimization for a given subset {R1, . . . , RK} and π is that of maximizing the throughput subject to equations (1) and (5):





  
    
      
        
          max
          
            
              
                N
                1
              
              ,
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ,
              
                
                  N
                  K
                
                :
              
            
            
              
                0
                ≤
                
                  N
                  1
                
                ≤
                
                  
                
                ⁢
                ...
                ⁢
                
                  
                
                ≤
                
                  N
                  K
                
              
              
                
                  
                    ∑
                    
                      k
                      =
                      1
                    
                    K
                  
                  ⁢
                  
                    
                      q
                      
                        π
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                    ⁢
                    
                      N
                      K
                    
                  
                
                ≤
                P
              
            
          
        
        ⁢
        
          
            ∑
            
              k
              =
              1
            
            K
          
          ⁢
          
            
              R
              
                π
                ⁡
                
                  (
                  k
                  )
                
              
            
            ⁢
            
              
                ∫
                0
                
                  N
                  k
                
              
              ⁢
              
                
                  
                    f
                    N
                  
                  ⁡
                  
                    (
                    u
                    )
                  
                
                ⁢
                
                  d
                  u
                
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




If the maximum throughput from the optimization problem of equation (9) is denoted T*(R1, . . . , RK, π), the following three propositions and theorem result:
Proposition 1: Let (N*1, . . . , N*K) achieve the maximum in equation (9), and suppose N1*=0. Then

T*(R1, . . . , RK, π)=T*(R′1, . . . , R′K−1, π′)

where Rπ(1) is dropped from the original rate subset RK to yield the reduced rate subset RK−1={R′1, . . . , R′K−1}, and π′ is the restriction of the original permutation π to RK−1.

Proposition 2: Suppose the noise power probability density function fN(•) is continuous and monotone decreasing. Then, for any rate subset {R1, . . . , RK} and permutation π, the throughput T(R1, . . . , RK, π, N1, . . . , NK), viewed as a function of N1, . . . , NK, is a concave continuously differentiable mapping from the open convex set RK++ into R.
Proposition 3: Given a rate subset {R1, . . . , RK} and permutation π: {1, . . . , K}→{1, . . . , K}, and noise power levels N1, . . . , NK:

    
    
        (a) The values of Pk obtained from equation (3) satisfy the following property:
    
    






  
    
      
        
          
            
              ∑
              
                i
                =
                1
              
              K
            
            ⁢
            
              P
              
                π
                ⁡
                
                  (
                  i
                  )
                
              
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              K
            
            ⁢
            
              
                q
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
              ⁢
              
                N
                i
              
            
          
        
        ,
      
    
    
      
        (
        10
        )
      
    
  
  
    
      
        
          
            q
            
              π
              ⁡
              
                (
                i
                )
              
            
          
          =
          
            
              γ
              
                π
                ⁡
                
                  (
                  i
                  )
                
              
            
            ⁢
            
              
                ∏
                
                  j
                  =
                  
                    i
                    +
                    1
                  
                
                K
              
              ⁢
              
                (
                
                  1
                  +
                  
                    γ
                    
                      π
                      ⁡
                      
                        (
                        j
                        )
                      
                    
                  
                
                )
              
            
          
        
        ,
        
          1
          ≤
          i
          ≤
          K
        
      
    
    
      
        (
        11
        )
      
    
  




    
    
        (b) For all m≧1,
    
    






  
    
      
        
          
            
              R
              
                π
                ⁡
                
                  (
                  i
                  )
                
              
            
            
              q
              
                π
                ⁡
                
                  (
                  i
                  )
                
              
            
          
          <
          
            
              
                ∑
                
                  j
                  =
                  0
                
                m
              
              ⁢
              
                R
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      j
                    
                    )
                  
                
              
            
            
              
                ∑
                
                  j
                  =
                  0
                
                m
              
              ⁢
              
                q
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      j
                    
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        12
        )
      
    
  




Theorem: Let the pdf of the noise power, fN(•), be monotone decreasing with infinite support. For a given rate subset R1, . . . , RK, permutation π, and maximum received power P, T(R1, . . . , RK, π, N) has a global maximum at (N*1, . . . , N*K) [over the set of all N* with 0<N1≦N2≦ . . . ≦NK] if and only if λ>0 exists such that




  
    
      
        
          
            
              
                f
                N
              
              ⁡
              
                (
                
                  N
                  i
                  *
                
                )
              
            
            ⁢
            
              
                R
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
              
                q
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
            
          
          =
          λ
        
        ,
        
          i
          =
          1
        
        ,
        ...
        ⁢
        
          
        
        ,
        K
        ,
      
    
    
      
        (
        13
        )
      
    
  
  
    
      
        
          
            ∑
            
              i
              =
              1
            
            K
          
          ⁢
          
            
              q
              
                π
                ⁡
                
                  (
                  i
                  )
                
              
            
            ⁢
            
              N
              i
              *
            
          
        
        =
        
          P
          .
        
      
    
    
      
        (
        14
        )
      
    
  




From equation (13), it follows that N*1>0 guarantees λ>0. Thus the only condition to be checked for a critical point (N*1, . . . , N*K) obtained by solving equation (13) is




  
    
      
        
          N
          1
          *
        
        =
        
          
            
              P
              -
              
                
                  ∑
                  
                    i
                    =
                    2
                  
                  K
                
                ⁢
                
                  
                    q
                    
                      π
                      ⁡
                      
                        (
                        i
                        )
                      
                    
                  
                  ⁢
                  
                    N
                    i
                    *
                  
                
              
            
            
              q
              
                π
                ⁡
                
                  (
                  i
                  )
                
              
            
          
          >
          0
        
      
    
    
      
        (
        15
        )
      
    
  




The theorem leads to a general algorithm for computing optimal noise levels for any decreasing noise pdf. However, the algorithm can be significantly simplified if the noise pdf is assumed to decrease in an exponential manner. Under the additional assumption that fN(•) is exponential, the optimal noise levels can be explicitly calculated in a recursive fashion for any permutation. This yields an algorithm to determine the maximum throughput and the throughput-maximizing set of discrete noise levels and the associated power allocation.
Let the noise power pdf be denoted by fN(n)=α exp(−αn), for n>0, where α>0 is a known parameter. Equation (13), evaluated for consecutive noise levels, can now be used to find a relation between N*i and N*1:




   
  
    
      
        
          
            
              
                
                  f
                  N
                
                ⁡
                
                  (
                  
                    N
                    i
                    *
                  
                  )
                
              
              ⁢
              
                
                  R
                  
                    π
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
                
                  q
                  
                    π
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
              
            
            =
            
              λ
              =
              
                
                  
                    
                      
                        f
                        N
                      
                      ⁡
                      
                        (
                        
                          N
                          
                            i
                            +
                            1
                          
                          *
                        
                        )
                      
                    
                    ⁢
                    
                      
                        R
                        
                          π
                          ⁡
                          
                            (
                            
                              i
                              +
                              1
                            
                            )
                          
                        
                      
                      
                        q
                        
                          π
                          ⁡
                          
                            (
                            
                              i
                              +
                              1
                            
                            )
                          
                        
                      
                    
                  
                  ⁢
                  
                    

                  
                  ⇔
                  
                    N
                    
                      i
                      +
                      1
                    
                    *
                  
                
                =
                
                  
                    N
                    i
                    *
                  
                  +
                  
                    x
                    
                      π
                      ⁡
                      
                        (
                        i
                        )
                      
                    
                  
                
              
            
          
          ,
          
            

          
          ⁢
          
            
              ⇔
              
                N
                i
                *
              
            
            =
            
              
                N
                1
                *
              
              +
              
                
                  ∑
                  
                    j
                    =
                    1
                  
                  
                    i
                    -
                    1
                  
                
                ⁢
                
                  x
                  
                    π
                    ⁡
                    
                      (
                      j
                      )
                    
                  
                
              
            
          
          ,
        
      
      
        
          (
          16
          )
        
      
    
  




in which the parameters xπ(i) are calculated according to the following equation:





  
    
      
        
          
            x
            
              π
              ⁡
              
                (
                i
                )
              
            
          
          =
          
            
              1
              α
            
            ⁢
            
              log
              ⁡
              
                (
                
                  
                    
                      R
                      
                        π
                        ⁡
                        
                          (
                          
                            i
                            +
                            1
                          
                          )
                        
                      
                    
                    
                      R
                      
                        π
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                    
                  
                  ⁢
                  
                    
                      
                        γ
                        
                          π
                          ⁡
                          
                            (
                            i
                            )
                          
                        
                      
                      ⁡
                      
                        (
                        
                          1
                          +
                          
                            γ
                            
                              π
                              ⁡
                              
                                (
                                
                                  i
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        )
                      
                    
                    
                      γ
                      
                        π
                        ⁡
                        
                          (
                          
                            i
                            +
                            1
                          
                          )
                        
                      
                    
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        17
        )
      
    
  




with N*1 determined from the power constraint:





   
  
    
      
        
          
            
              ∑
              
                j
                =
                1
              
              K
            
            ⁢
            
              
                q
                
                  π
                  ⁡
                  
                    (
                    j
                    )
                  
                
              
              ⁢
              
                N
                j
                *
              
            
          
          =
          
            
              P
              ⁢
              
                

              
              ⇔
              
                N
                1
                *
              
            
            =
            
              
                
                  P
                  -
                  
                    
                      ∑
                      
                        j
                        =
                        1
                      
                      
                        K
                        -
                        1
                      
                    
                    ⁢
                    
                      
                        (
                        
                          
                            ∑
                            
                              i
                              =
                              
                                j
                                +
                                1
                              
                            
                            K
                          
                          ⁢
                          
                            q
                            
                              π
                              ⁡
                              
                                (
                                i
                                )
                              
                            
                          
                        
                        )
                      
                      ⁢
                      
                        x
                        
                          π
                          ⁡
                          
                            (
                            i
                            )
                          
                        
                      
                    
                  
                
                
                  
                    ∑
                    
                      j
                      =
                      1
                    
                    K
                  
                  ⁢
                  
                    q
                    
                      π
                      ⁡
                      
                        (
                        j
                        )
                      
                    
                  
                
              
              .
            
          
        
      
      
        
          (
          18
          )
        
      
    
  




The only remaining feasibility condition of the above solution that remains to be checked is:




  
    
      
        
          
            N
            1
            *
          
          >
          0
        
        ⇔
        
          
            
              ∑
              
                j
                =
                1
              
              
                K
                -
                1
              
            
            ⁢
            
              
                (
                
                  
                    ∑
                    
                      i
                      =
                      
                        j
                        +
                        1
                      
                    
                    K
                  
                  ⁢
                  
                    q
                    
                      π
                      ⁡
                      
                        (
                        i
                        )
                      
                    
                  
                
                )
              
              ⁢
              
                x
                
                  π
                  ⁡
                  
                    (
                    j
                    )
                  
                
              
            
          
          <
          
            P
            .
          
        
      
    
    
      
        (
        19
        )
      
    
  




The above results may now be put together to result in the following algorithm for determining the throughput-optimal set of transmission rates and the corresponding set of discrete noise levels when the noise pdf is exponentially decreasing.

    
    
        Step 1: Let K=1 (this is the number of transmission rates chosen from the original rate set).
        Step 2: Choose a subset of K rates from the original set, denoted {R1, . . . , RK}.
        Step 3: For each permutation π, do the following:
        
            3.1: Determine the parameters qπ(i) and xπ(i) for the chosen rate subset, for i=1, . . . , K−1, from equations (11) and (17) respectively.
            3.2: If equation (19) is satisfied, then a feasible superposition code exists supporting this rate subset with decoding order specified by π, and the noise levels can be calculated from equations (17), (18) and (16), in that order. Finally the achieved throughput for the chosen rate subset R is given by:
        
        
    
    






  
    
      
        
          
            T
            *
          
          ⁡
          
            (
            
              
                R
                1
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                R
                K
              
              ,
              π
            
            )
          
        
        =
        
          
            ∑
            
              i
              =
              1
            
            K
          
          ⁢
          
            
              [
              
                1
                -
                
                  e
                  
                    
                      -
                      α
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      N
                      i
                      *
                    
                  
                
              
              ]
            
            ⁢
            
              
                R
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
              .
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




    
    
        
        
            3.3: If on the other hand equation (19) is not satisfied, no superposition code exists using decoding order π supporting the chosen rate subset for the given power constraint.
            3.4: Repeat Steps 3.1-3.3 for each permutation on the given rate subset {R1, . . . , RK} chosen in Step 2.
        
        
        Step 4: For a given K, choose all possible rate subsets of size K and repeat the operations in Steps 3.1-3.4.
        Step 5: Repeat the operations in Steps 2 through 4 for all possible sizes of the rate subset by varying K from 1 to M, where M is the number of rates in the original rate subset, from which the rate subset has to be chosen. At the end of this step, the throughput values for all possible feasible rate subsets and permutations should be in hand.
        Step 6: Choose the largest value of the achievable throughputs obtained in Steps 3.1-3.4 for the different rate subsets under consideration.
    
    


Note that the above procedure can be repeated for different received power constraints P to obtain the maximum achieved throughput under superposition coding as a function of the received power, or equivalently as a function of the transmit power (if the relation between transmit and receive power is known).
Further, the algorithm works for arbitrary monotone decreasing fN(•) if equations (17), (18) and (16) are all replaced by equations (13) and (14) to solve for the optimal noise power levels and if equation (19) is replaced by equation (15) in which the appropriate values for the noise power levels are used. Finally, in Step 3.2, the throughput T for a given rate subset and a given permutation π is calculated according to equation (6).
Having set forth an algorithm for determining the throughput-optimal set of transmission rates and the corresponding set of discrete noise levels, a system configured to operate in accordance with the principles of the algorithm may now be illustrated. Accordingly, FIG. 1 illustrates a block diagram of one embodiment of a system for determining the throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates constructed according to the principles of the present invention.
A database 110 contains the finite, predetermined set of allowed transmission rates from which the subsets are chosen for analysis. In the illustrated embodiment, a conventional or later-established industry transmission standard determines the set. Proprietary transmission rates may alternatively occupy the set.
A noise level calculator 120 is coupled to the database 110. The noise level calculator 120 is configured to calculate associated optimal noise levels for at least some rate subsets of the set of allowed transmission rates. In doing so, the illustrated embodiment of the noise level calculator 120 chooses subsets containing one or more transmission rates, but of fewer members than the entire set.
In the illustrated embodiment, the noise level calculator is configured to calculate associated optimal noise levels for all rate subsets of the set of allowed transmission rates as the algorithm set forth above does. Further, the illustrated embodiment of the noise level calculator employs equations that at least approximate equations (11) and (17) above to find the parameters qπ(i) and xπ(i) for each rate subset. N1 and Ni can also be calculated as desired. “At least approximate” means that the equations may not be exactly as those set forth above, but achieve substantially the same result. Of course, any suitable equation or equations for determining optimal noise levels fall within the broad scope of the present invention.
A code feasibility calculator 130 is associated with the noise level calculator 120. The code feasibility calculator 130 is configured to determine whether feasible superposition codes exist for the at least some rate subsets given the associated optimal noise levels. The rate subsets lacking feasible superposition codes are disregarded. The rate subsets having feasible superposition codes are suitable for subsequent wireless transmissions.
The illustrated embodiment of the code feasibility calculator 130 is further configured to determine throughputs associated with the ones of the at least some rate subsets. Further, the illustrated embodiment of the code feasibility calculator 130 employs an equation that at least approximates equation (20), respectively equation (6). Of course, any suitable equation or equations for determining throughputs fall within the broad scope of the present invention.
The illustrated embodiment of the code feasibility calculator 130 is further configured to select a one of the at least some rate subsets having a maximum associated throughput. It is assumed that maximizing throughput is the desired objective, though the broad scope of the present invention is not so limited.
Having selected the one (feasible) rate subset having the maximum associated throughput, the illustrated embodiment of the code feasibility calculator 130 is further configured to calculate power allocations for rates in the selected rate subset. The code feasibility calculator 130 can then provide the appropriate rate subset and power allocation for the various rates in the subset to a wireless network 140 for subsequent transmissions.
Turning now to FIG. 2, illustrated is a flow diagram of one embodiment of a method of determining the throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates carried out according to the principles of the present invention.
The method begins in a start step 210, wherein it is desired to determine channel throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates. The method proceeds to a step 220 in which associated optimal noise levels are calculated for at least some rate subsets of the set of allowed transmission rates. Next, in a step 230, the existence of feasible superposition codes for the at least some rate subsets given the associated optimal noise levels is determined. Rate subsets that are suitable for subsequent wireless transmissions are thereby identified.
Then, in a step 240, associated throughputs for the suitable rate subsets are determined. Next, in a step 250, the rate subset having the maximum associated throughput is selected. The method ends in an end step 260 wherein, in the illustrated embodiment, power levels for the rates in the selected rate subset are determined and given over to a wireless network for subsequent transmissions.
Some numerical results will now be presented to illustrate the performance of the superposition coding scheme when the pdf of noise power is exponential with parameter α, for various choices of α. The set of available rates (again expressed in nats/s/Hz) is {2, 1, 2⁄3, 1⁄2, 1⁄3, 1⁄4, 1⁄5}. The power P and the noise power levels N*1, . . . , N*K are normalized by dividing them by the mean noise power N=1/α. Thus P may be replaced in the optimization problem by the signal-to-noise ratio (SNR) P/ N, where N*1, . . . , N*K now refer to the noise power levels normalized by the mean.
The results are plotted in FIGS. 3, 4 and 5. In each of them, SNR is plotted versus the maximal throughput for this SNR. The maximal throughput is the throughput given by equation (20) for the best rate subset and the best permutation over this rate subset, with N*1 . . . , N*K obtained from equations (17), (18) and (16).
Turning to FIG. 3, shown are throughput versus SNR curves for different values of the mean noise power I/α, when the pdf of the noise power is exponential. For a given SNR, the maximal throughput increases as mean noise power decreases, i.e., as α increases. Turning to FIG. 4, α is set equal to 0.7 and the maximal throughput is plotted using a superposition code versus SNR, and this is compared to the maximal throughput when only one rate is permitted to be used out of the rate subset. The difference in throughput between the two curves represents the gain due to superposition coding.
For low SNRs, the two throughputs are close, but start to diverge as SNR increases. The reason for this widening gap is that as more power is available, the superposition code can take advantage of it to use rate subsets with finer granularity than if only use a single rate is used. In particular, the slope of the single-rate throughput curve changes discontinuously at an SNR of approximately 18 dB. At this point, the SNR is high enough that a higher rate is feasible as the single-rate code, so this curve then starts to climb up toward the superposition code curve.
In FIG. 5, the sum of the rates in the optimal rate subset is plotted for α=0.7, for both superposition and single-rate coding versus SNR. As available power rises, the sum of the rates in the optimal rate subset rises as well, except for the initial drop in the superposition code sum-rate at very low SNR. However, since the single-rate code can only use one rate, this optimal rate rises in bigger jumps than for the superposition code, which can use higher rates when power is available, with finer granularity than if only a single rate is used.
In general, the size of the optimal rate subset used in the superposition code increases with SNR. However, a pattern could not be readily identified in the optimal permutations over these optimal rate subsets. In fact, situations arose where the optimal rate subset was the same at two different SNRs, but the optimal permutations were different.
Although the present invention has been described in detail, those skilled in the art should understand that they can make various changes, substitutions and alterations herein without departing from the spirit and scope of the invention in its broadest form.

## Claims

1. A system for determining channel throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates, comprising:
a noise level calculator configured to calculate associated optimal noise levels for at least some rate subsets of said set of allowed transmission rates; and
a code feasibility calculator associated with said noise level calculator and configured to determine whether feasible superposition codes exist for said at least some rate subsets given said associated optimal noise levels and thereby identify ones of said at least some rate subsets that are suitable for subsequent wireless transmissions.

2. The system as recited in claim 1 wherein said noise level calculator is configured to calculate said associated optimal noise levels using an equation that at least approximates:




  
    
      q
      
        π
        ⁡
        
          (
          i
          )
        
      
    
    =
    
      
        γ
        
          π
          ⁡
          
            (
            i
            )
          
        
      
      ⁢
      
        
          ∏
          
            j
            =
            
              i
              +
              1
            
          
          K
        
        ⁢
        
          (
          
            1
            +
            
              γ
              
                π
                ⁡
                
                  (
                  j
                  )
                
              
            
          
          )
        
      
    
  
  ,
  
    
  
  ⁢
  
    1
    ≤
    i
    ≤
    K
  
  ,



where:
K is a number of rates in a subset Rk of said predetermined set of allowed transmission rates,
γπ(i)=exp(2Rπ(i))−1, for a permutation π of said subset, and
Rπ(i) is the ith transmission rate in said permutation.

3. The system as recited in claim 1 wherein said noise level calculator is configured to calculate said associated optimal noise levels using an equation that at least approximates:




  
    
      x
      
        π
        ⁡
        
          (
          i
          )
        
      
    
    =
    
      
        1
        α
      
      ⁢
      
        log
        ⁡
        
          (
          
            
              
                R
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
              
                R
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
            
            ⁢
            
              
                
                  γ
                  
                    π
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
                ⁡
                
                  (
                  
                    1
                    +
                    
                      γ
                      
                        π
                        ⁡
                        
                          (
                          
                            i
                            +
                            1
                          
                          )
                        
                      
                    
                  
                  )
                
              
              
                γ
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
            
          
          )
        
      
    
  
  ,



where:
α is an adjustable parameter,
γπ(i)=exp(2Rπ(i))−1, for a permutation π of a subset Rk of said predetermined set of allowed transmission rates, and
Rπ(i) is an ith transmission rate in said permutation.

4. The system as recited in claim 1 wherein said noise level calculator is configured to calculate relationships between noise levels using equations that at least approximate:




  
    
      N
      1
      *
    
    =
    
      
        
          
            P
            -
            
              
                ∑
                
                  j
                  =
                  1
                
                
                  K
                  -
                  1
                
              
              ⁢
              
                
                  (
                  
                    
                      ∑
                      
                        i
                        =
                        
                          j
                          +
                          1
                        
                      
                      K
                    
                    ⁢
                    
                      q
                      
                        π
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                    
                  
                  )
                
                ⁢
                
                  x
                  
                    π
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
              
            
          
          
            
              ∑
              
                j
                =
                1
              
              K
            
            ⁢
            
              q
              
                π
                ⁡
                
                  (
                  j
                  )
                
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        
          N
          i
          *
        
      
      =
      
        
          N
          1
          *
        
        +
        
          
            ∑
            
              j
              =
              1
            
            
              i
              -
              1
            
          
          ⁢
          
            x
            
              π
              ⁡
              
                (
                j
                )
              
            
          
        
      
    
  
  ,



where:
P is a total received power,
K is a number of rates in a subset Rk of said predetermined set of allowed transmission rates,






  
    
      q
      
        π
        ⁡
        
          (
          i
          )
        
      
    
    =
    
      
        γ
        
          π
          ⁡
          
            (
            i
            )
          
        
      
      ⁢
      
        
          ∏
          
            j
            =
            
              i
              +
              1
            
          
          K
        
        ⁢
        
          (
          
            1
            +
            
              γ
              
                π
                ⁡
                
                  (
                  j
                  )
                
              
            
          
          )
        
      
    
  
  ,
  
    1
    ≤
    i
    ≤
    K
  
  ,
  
    

  
  ⁢
  
    
      x
      
        π
        ⁡
        
          (
          j
          )
        
      
    
    =
    
      
        1
        α
      
      ⁢
      
        log
        ⁡
        
          (
          
            
              
                R
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
              
                R
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
            
            ⁢
            
              
                
                  γ
                  
                    π
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
                ⁡
                
                  (
                  
                    1
                    +
                    
                      γ
                      
                        π
                        ⁡
                        
                          (
                          
                            i
                            +
                            1
                          
                          )
                        
                      
                    
                  
                  )
                
              
              
                γ
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
            
          
          )
        
      
    
  
  ,




α is an adjustable parameter,
γπ(i)=exp(2Rπ(i))−1, for a permutation π of said subset, and
Rπ(i) is an lth transmission rate in said permutation.

5. The system as recited in claim 1 wherein said code feasibility calculator is further configured to determine throughputs associated with said ones of said at least some rate subsets.

6. The system as recited in claim 5 wherein said code feasibility calculator is configured to determine throughputs using an equation that at least approximates:




  
    
      
        T
        *
      
      ⁡
      
        (
        
          
            R
            1
          
          ,
          ...
          ⁢
          
            
          
          ,
          
            R
            K
          
          ,
          π
        
        )
      
    
    =
    
      
        ∑
        
          i
          =
          1
        
        K
      
      ⁢
      
        
          [
          
            1
            -
            
              e
              
                
                  -
                  α
                
                ⁢
                
                  
                
                ⁢
                
                  N
                  i
                  *
                
              
            
          
          ]
        
        ⁢
        
          R
          
            π
            ⁡
            
              (
              i
              )
            
          
        
      
    
  
  ,



where:
Rk is a kth transmission rate in a subset of K of said predetermined set of allowed transmission rates,
π is a permutation of said subset,
Rπ(i) is an ith transmission rate in said permutation,
α is an adjustable parameter, and
Ni* is a noise power level associated with said ith transmission rate, wherein a set of N1*, . . . ,Nk* is associated with a maximum of said channel throughput.

7. The system as recited in claim 5 wherein said code feasibility calculator is further configured to select a one of said at least some rate subsets having a maximum associated throughput.

8. The system as recited in claim 1 wherein said code feasibility calculator is further configured to calculate power allocations for rates in at least one of said at least some rate subsets.

9. The system as recited in claim 8 wherein said code feasibility calculator is further configured to calculate said power allocations using an equation that at least approximates:




  
    
      P
      
        π
        ⁢
        
          
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        γ
        
          π
          ⁢
          
            
          
          ⁢
          
            (
            k
            )
          
        
      
      ⁡
      
        [
        
          
            
              ∑
              
                i
                =
                1
              
              
                k
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              P
              
                π
                ⁢
                
                  
                
                ⁢
                
                  (
                  i
                  )
                
              
            
          
          +
          
            N
            k
          
        
        ]
      
    
  
  ,



where:
γπ(k)=exp(2Rπ(k))−1, for a permutation π of a subset Rk of said predetermined set of allowed transmission rates,
Pπ(i) is an allocated power associated with an ith transmission rate in said permutation, and
Nk is a noise power level associated with a kth transmission rate in said subset.

10. The system as recited in claim 1 wherein said at least some rate subsets include subsets of single rates.

11. The system as recited in claim 1 wherein said noise level calculator is configured to calculate associated optimal noise levels for all rate subsets of said set.

12. The system as recited in claim 1 wherein a transmission standard determines said set.

13. A method of determining channel throughput achievable under superposition coding given a finite, predetermined set of allowed transmission rates, comprising:
calculating associated optimal noise levels for at least some rate subsets of said set of allowed transmission rates; and
determining whether feasible superposition codes exist for said at least some rate subsets given said associated optimal noise levels and thereby identifying ones of said at least some rate subsets that are suitable for subsequent wireless transmissions.

14. The method as recited in claim 13 wherein said calculating comprises calculating said associated optimal noise levels using an equation that at least approximates:




  
    
      q
      
        π
        ⁢
        
          
        
        ⁢
        
          (
          i
          )
        
      
    
    =
    
      
        γ
        
          π
          ⁢
          
            
          
          ⁢
          
            (
            i
            )
          
        
      
      ⁢
      
        
          ∏
          
            j
            =
            
              i
              +
              1
            
          
          K
        
        ⁢
        
          
        
        ⁢
        
          (
          
            1
            +
            
              γ
              
                π
                ⁢
                
                  
                
                ⁢
                
                  (
                  j
                  )
                
              
            
          
          )
        
      
    
  
  ,
  
    
  
  ⁢
  
    1
    ≤
    i
    ≤
    K
  
  ,



where:
K is a number of rates in a subset Rk of said predetermined set of allowed transmission rates,
γπ(i)=exp(2Rπ(i))−1, for a permutation π of said subset, and
Rπ(i) is the ith transmission rate in said permutation π.

15. The method as recited in claim 13 wherein said calculating comprises calculating said associated optimal noise levels using an equation that at least approximates:




  
    
      x
      
        π
        ⁢
        
          
        
        ⁢
        
          (
          i
          )
        
      
    
    =
    
      
        1
        α
      
      ⁢
      
        log
        ⁡
        
          (
          
            
              
                R
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
              
                R
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    (
                    i
                    )
                  
                
              
            
            ⁢
            
              
                
                  γ
                  
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    
                      (
                      i
                      )
                    
                  
                
                ⁡
                
                  (
                  
                    1
                    +
                    
                      γ
                      
                        π
                        ⁢
                        
                          
                        
                        ⁢
                        
                          (
                          
                            i
                            +
                            1
                          
                          )
                        
                      
                    
                  
                  )
                
              
              
                γ
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
            
          
          )
        
      
    
  
  ,



where:
α is an adjustable parameter,
γπ(i)=exp(2Rπ(i))−1, for a permutation π of a subset Rk of said predetermined set of allowed transmission rates, and
Rπ(i) is an ith transmission rate in said permutation.

16. The method as recited in claim 13 wherein said calculating comprises calculating relationships between noise levels using equations that at least approximate:




  
    
      N
      1
      *
    
    =
    
      
        
          
            P
            -
            
              
                ∑
                
                  j
                  =
                  1
                
                
                  K
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                
                  (
                  
                    
                      ∑
                      
                        i
                        =
                        
                          j
                          +
                          1
                        
                      
                      K
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      q
                      
                        π
                        ⁢
                        
                          
                        
                        ⁢
                        
                          (
                          i
                          )
                        
                      
                    
                  
                  )
                
                ⁢
                
                  x
                  
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    
                      (
                      i
                      )
                    
                  
                
              
            
          
          
            
              ∑
              
                j
                =
                1
              
              K
            
            ⁢
            
              
            
            ⁢
            
              q
              
                π
                ⁢
                
                  
                
                ⁢
                
                  (
                  j
                  )
                
              
            
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        
          N
          i
          *
        
      
      =
      
        
          N
          1
          *
        
        +
        
          
            ∑
            
              j
              =
              1
            
            
              i
              -
              1
            
          
          ⁢
          
            
          
          ⁢
          
            x
            
              π
              ⁢
              
                
              
              ⁢
              
                (
                j
                )
              
            
          
        
      
    
  
  ,



where:
P is a total received power,
K is a number of rates in a subset Rk of said predetermined set of allowed transmission rates,






  
    
      q
      
        π
        ⁡
        
          (
          i
          )
        
      
    
    =
    
      
        γ
        
          π
          ⁡
          
            (
            i
            )
          
        
      
      ⁢
      
        
          ∏
          
            j
            =
            
              i
              +
              1
            
          
          K
        
        ⁢
        
          (
          
            1
            +
            
              γ
              
                π
                ⁡
                
                  (
                  j
                  )
                
              
            
          
          )
        
      
    
  
  ,
  
    1
    ≤
    i
    ≤
    K
  
  ,
  
    

  
  ⁢
  
    
      x
      
        π
        ⁡
        
          (
          j
          )
        
      
    
    =
    
      
        1
        α
      
      ⁢
      
        log
        ⁡
        
          (
          
            
              
                R
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
              
                R
                
                  π
                  ⁡
                  
                    (
                    i
                    )
                  
                
              
            
            ⁢
            
              
                
                  γ
                  
                    π
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
                ⁡
                
                  (
                  
                    1
                    +
                    
                      γ
                      
                        π
                        ⁡
                        
                          (
                          
                            i
                            +
                            1
                          
                          )
                        
                      
                    
                  
                  )
                
              
              
                γ
                
                  π
                  ⁡
                  
                    (
                    
                      i
                      +
                      1
                    
                    )
                  
                
              
            
          
          )
        
      
    
  
  ,




α is an adjustable parameter,
γπ(i)=exp(2Rπ(i))−1, for a permutation π of said subset, and
Rπ(i) is an ith transmission rate in said permutation.

17. The method as recited in claim 13 wherein said determining comprises determining throughputs associated with said ones of said at least some rate subsets.

18. The method as recited in claim 17 wherein said determining comprises determining throughputs using an equation that at least approximates:




  
    
      
        T
        *
      
      ⁡
      
        (
        
          
            R
            1
          
          ,
          ...
          ⁢
          
            
          
          ,
          
            R
            K
          
          ,
          π
        
        )
      
    
    =
    
      
        ∑
        
          i
          =
          1
        
        K
      
      ⁢
      
        
      
      ⁢
      
        
          [
          
            1
            -
            
              e
              
                
                  -
                  α
                
                ⁢
                
                  
                
                ⁢
                
                  N
                  i
                  *
                
              
            
          
          ]
        
        ⁢
        
          R
          
            π
            ⁢
            
              
            
            ⁢
            
              (
              i
              )
            
          
        
      
    
  
  ,



where:
Rk is a kth transmission rate in a subset of K of said predetermined set of allowed transmission rates,
π is a permutation of said subset,
Rπ(i) is an ith transmission rate in said permutation,
α is an adjustable parameter, and
Ni* is a noise power level associated with said ith transmission rate, wherein a set of N1*, . . . ,Nk* is associated with a maximum of said channel throughput.

19. The method as recited in claim 17 further comprising selecting a one of said at least some rate subsets having a maximum associated throughput.

20. The method as recited in claim 13 wherein said determining comprises calculating power allocations for rates in at least one of said at least some rate subsets.

21. The method as recited in claim 20 wherein said calculating comprises calculating said power allocations using an equation that at least approximates:




  
    
      P
      
        π
        ⁢
        
          
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        γ
        
          π
          ⁢
          
            
          
          ⁢
          
            (
            k
            )
          
        
      
      ⁡
      
        [
        
          
            
              ∑
              
                i
                =
                1
              
              
                k
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              P
              
                π
                ⁢
                
                  
                
                ⁢
                
                  (
                  i
                  )
                
              
            
          
          +
          
            N
            k
          
        
        ]
      
    
  
  ,



where:
γπ(k)=exp(2Rπ(k))−1, for a permutation π of a subset Rk of said predetermined set of allowed transmission rates,
Pπ(i) is an allocated power associated with an ith rate in said permutation, and
Nk is a noise power level associated with a kth transmission rate in said subset.

22. The method as recited in claim 13 wherein said at least some rate subsets include subsets of single rates.

23. The method as recited in claim 13 further comprising calculating associated optimal noise levels for all rate subsets of said set.

24. The method as recited in claim 13 wherein a transmission standard determines said set.