# Patent

## Title
MAC design for wireless hot-spot networks

## Abstract

In widely deployed wireless “hot-spot” networks, nodes frequently join or leave, inelastic/elastic and saturated/non-satuarted flows coexist. In such dynamic and diverse environments, it is challenging to maximize the channel utilization while providing satisfactory user experiences. In this invention, one proposes a novel contention-on-demand (CoD) MAC scheme to address this problem. The CoD scheme consists of a fixed-CW algorithm, a dynamic-CW algorithm, and an admission control unit. The fixed-CW algorithm allows elastic flows to access limited system bandwidth; the dynamic-CW algorithm enables inelastic flows to contend for channel on demand and quickly adapt to network change; and the admission control unit rejects overloaded traffic for providing good user experiences. One then performs an asymptotic analysis to develop a simple but efficient admission control rule. Finally, extensive simulations verify that the scheme is very effective and the theoretical result is very accurate.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of U.S. Provisional Patent Application No. 62/243,004, filed on Oct. 17, 2015, which is incorporated by reference herein in its entirety.


COPYRIGHT NOTICE
A portion of the disclosure of this patent document contains material, which is subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or the patent disclosure, as it appears in the Patent and Trademark Office patent file or records, but otherwise reserves all copyright rights whatsoever.
FIELD OF THE INVENTION
The present invention relates to a method for providing admission control to a wireless network system.
BACKGROUND
The following references are cited in the specification. Disclosures of these references are incorporated herein by reference in their entirety.
LIST OF REFERENCES


    
    
        [1]. ANSI/IEEE Std 802.11, Part 11: wireless LAN medium access control (MAC) and physical layer (PHY) specifications, 1999 Edition (R2007).
        [2]. A. Malika, J. Qadirb, B. Ahmada, K. A. Yauc, and U. Ullaha, “Qos in IEEE 802.11-based wireless networks: A contemporary review,” Journal of Network and Computer Applications, vol. 55, pp. 24-46, 2015.
        [3]. M. Burton, “Hotspot 2.0 and the next generation hotspot,” April 2013. Available: https://www.cwnp.com/hotspot-2-0-and-the-next-generationhotspot.
        [4]. K. Fitchard, “Hotspot 2.0 inches its way into public wi-fi networks,” March 2014. Available: https://gigaom.com/2014/03/07/hotspot-2-0-inches-its-way-into-public-wi-fi-networks/.
        [5]. https://en.wikipedia.org/wiki/Hotspot (Wi-Fi).
        [6]. Q. L. Zhao, D. H. K. Tsang, and T. Sakurai, “A simple critical-offeredload-based CAC scheme for IEEE 802.11 DCF networks,” IEEE/ACM Transactions on Networking, vol. 19, no. 5, pp. 1485-1498, October 2011.
        [7]. D. Chiu and R. Jain, “Analysis of the increase and decrease algorithms for congestion avoidance in computer networks,” Journal of Computer Networks and ISDN, vol. 17, pp. 1-14, 1989.
        [8]. V. Ramaiyan, A. Kumar, and E. Altman, “Fixed point analysis of single cell IEEE 802.11e WLANs: Uniqueness and multistability,” IEEE/ACM Transactions on Networking, vol. 16, no. 5, pp. 1080-1093, October 2008. The corresponding technical report is avalable: http://ece.iisc.ernet.in/anurag/papers/anurag/ramaiyan-eta105fixedpoint-general.pdf.gz.
        [9]. Q. L. Zhao, D. H. K. Tsang, and T. Sakurai, “A comprehensive and accurate non-saturated IEEE 802.11e EDCA model for an arbitrary buffer size,” IEEE Transactions on Mobile Computing, vol. 12, no. 12, pp. 2455-2469, 2013.
        [10]. M. Heusse, F. Rousseau, R. Guillier, and A. Duda, “Idle sense: An optimal access method for high throughput and fairness in rate diverse wireless lans,” Proceedings of ACM Sigcomm, 2005.
        [11]. S. H. Nguyen, H. L. Vu, and L. L. H. Andrew, “Performance analysis of ieee 802.11 wlans with saturated and unsaturated sources,” IEEE T. Vehicular Technology, vol. 61, no. 1, pp. 333-345, 2012.
        [12]. G. Bianchi, “Performance analysis of the IEEE 802.11 distributed coordination function,” IEEE Journal on Selected Areas in Communications, vol. 18, no. 3, pp. 535-547, March 2000.
        [13]. http://www.isi.edu/nsnam/ns/ns-build.html.
        [14]. D. Malone, K. Duffy, and D. Leith, “Modeling the 802.11 distributed coordination function in non-saturated heterogeneous conditions,” IEEE/ACM Transactions on Networking, vol. 15, no. 1, pp. 159-172, February 2007.
        [15]. R. M. Corlessa, G. H. Gonnet, D. E. G. Hare, D. J. Jeffrey, and D. E. Knuth, “On the lambert w function,” Adv. Comput. Math., vol. 5, pp. 329-359, 1996.
        [16]. M. Nassiri, M. Heusse, and A. Duda, “A novel access method for supporting absolute and proportional priorities in 802.11 WLANs,” INFOCOM, pp. 709-717, 2008.
        [17]. Z. Ma, Q. Zhao, and T. H. Luan, “Providing utility-optimal throughput guarantees in wireless lans,” IEEE Transactions Vehicular Technology. Accepted for publication.
        [18]. A. Kumar, E. Altman, D. Miorandi, and M. Goyal, “New insights from a fixed point analysis of single cell ieee 802.11 wlans,” IEEE/ACM Transactions on Networking, vol. 15, no. 3, pp. 588-601, March 2007.
    
    


IEEE 802.11 standard [1] has become the dominating solution for wireless “hot-spot” networks [2][3][4][5] at many public places, such as hotels, airports, restaurants, and malls. The salient features of “hot-spot” networks include (i) nodes frequently join and leave, and (ii) traffic is very diverse and heterogeneous. For example, in “hot-spot” networks, inelastic applications (such as video, audio) and elastic applications (such as reading Web pages, e-mail) coexist, and saturated (where a node always has packets to transmit) and non-saturated flows coexist.
However, current 802.11 networks have the following drawbacks: (i) the network capacity is not fully exploited since the protocol parameters such as contention window (CW) is statically configured, failing to adapt to the change in the node number and the packet size, (ii) the system performance (such as throughput and delay) deteriorates significantly if the offered load slightly exceeds a critical value [6], (iii) there is not an admission control unit that limits the admitted traffic loads. These drawbacks naturally lead to a very bad user experience in such dynamic and diverse “hot-spot” networks if traffic becomes congested.
In related works, [8][9], respectively, analyzed the saturated and non-saturation performance of 802.11 networks. [10] proposed a scheme to maximize the channel utilization of saturated 802.11 networks. [11] developed a performance model to analyze the 802.11 performance when non-saturated and saturated flows coexist. However, none of these works consider the objectives: designing a distributed scheme for maximizing the channel utilization and providing satisfactory QoS when inelastic and elastic flows coexist.
SUMMARY OF THE INVENTION
In this invention, one targets to devise a contention-on-demand (CoD) MAC scheme for fully utilizing the available system bandwidth and providing satisfactory user experiences, in wireless “hot-spot” networks. The CoD scheme consists of two distributed algorithms (i.e., a fixed-CW algorithm and a dynamic-CW algorithm) and an admission control unit. The fixed-CW algorithm is a variant of the 802.11 inherent binary-exponential-backoff (BEB) algorithm; this algorithm limits the use of the system bandwidth and is adopted by each elastic flow. The dynamic-CW algorithm is built on the famous additive-increase-multiplicative-decrease (AIMD) rule [7]; this algorithm will enable each inelastic flow to contend for neither more nor less than its required bandwidth, quickly adapt to the change of network states, as well as maximize the channel utilization. The admission control unit limits the admitted traffic load for providing satisfactory user experiences. One then performs an asymptotic analysis on the system throughput of the inelastic flows. On this basis, one develops a simple but efficient admission control rule. Finally, extensive simulations verify that the CoD scheme is very effective and the theoretical result is very accurate.
An aspect of the present invention is to provide methods for providing admission control to a wireless network system are set forth as preferred examples.
According to an embodiment of the present invention, a computer-implemented method for providing admission control to a wireless network system, the system having one or more high-priority (HP) nodes, one or more low-priority (LP) nodes, contending for access to an access point, the method comprises:

    
    
        categorizing the system into one LP access category (AC) and I HP ACs, wherein the LP AC has n0 LP nodes, and each of the LP nodes has a same packet size L0 and generates a random backoff count uniformly distributed in [0,CW0] for each of new transmission or retransmission, where CW0 is a pre-configured contention window (CW) size required by a fixed-CW algorithm, and wherein the HP AC i, where 1≦i≦I, has ni HP nodes, where n1+L+nI=n is a total number of HP nodes, and each of the HP nodes in the HP AC i has a same packet size Li and a same packet arrival rate λi, and generates a random backoff count uniformly distributed in [0,CWD] for each of new transmission or retransmission, where CWD is a contention window size dynamically set by a dynamic-CW algorithm;
        setting all of the LP and HP nodes having the packet size with L;
        determining an optimal asymptotic HP attempt rate as follows:
    
    






  
    β
    opt
    +
  
  =
  
    
      
        W
        0
      
      ⁡
      
        (
        
          
            
              C
              0
            
            ⁡
            
              (
              
                σ
                -
                
                  T
                  c
                
              
              )
            
          
          
            eT
            c
          
        
        )
      
    
    +
    1
  




    
    
        where w0(·) is a Lambert W(z) function, W(z)eW(z)=z, for any complex number z,
        σ is a length of a media access control MAC slot,
        C0(1−β0)n0,
        n0 is the total number of LP nodes,
        β0 is an average attempt rate per slot for each of the LP AC nodes,
        Tc is a mean time in slots for an unsuccessful transmission;
        determining an optimal asymptotic HP throughput Γ1(βopt+) by substituting the βopt+ for β+ in Γ1(β+) as follows:
    
    






  
    
      Γ
      1
    
    ⁡
    
      (
      
        β
        +
      
      )
    
  
  =
  
    
      
        e
        
          -
          
            β
            +
          
        
      
      ⁢
      
        β
        +
      
      ⁢
      
        C
        0
      
      ⁢
      L
    
    
      
        T
        c
      
      +
      
        
          
            e
            
              -
              
                β
                +
              
            
          
          ⁡
          
            [
            
              σ
              -
              
                T
                c
              
            
            ]
          
        
        ⁢
        
          C
          0
        
      
    
  




    
    
        where Γ1(β+) is an asymptotic total HP throughput, and
        β+ is a total asymptotic HP attempt rate;
        determining one or more numerical values of operating parameters including λi such that Σi=1IniλiL<Γ1(βopt+),
        where Σi=1IniλiL is a total traffic load of all of the HP ACs;
        determining the CWD based on a requirement of total delay of a flow.
    
    


Preferably, the total delay of the flow is related to the determined numerical values of λi.
Preferably, the total traffic load of all of the HP ACs is just slightly below the optimal asymptotic HP throughput.
According to another embodiment of the present invention, a computer-implemented method for providing admission control to a wireless network system, the system having one or more high-priority (HP) nodes, one or more low-priority (LP) nodes, contending for access to an access point, the method comprises:

    
    
        categorizing the system into one LP access category (AC) and I HP ACs, wherein the LP AC has n0 LP nodes, and each of the LP nodes has a same packet size L0 and generates a random backoff count uniformly distributed in [0, CW0] for each of new transmission or retransmission, where CW0 is a pre-configured contention window (CW) size required by a fixed-CW algorithm, and wherein the HP AC i, where 1≦i≦I, has ni HP nodes, where n1+L+nI=n is a total number of HP nodes, and each of the HP nodes in the HP AC i has a same packet size Li and a same packet arrival rate λi, and generates a random backoff count uniformly distributed in [0,CWD] for each of new transmission or retransmission, where CWD is a contention window size dynamically set by a dynamic-CW algorithm;
        determining an optimal asymptotic aggregate attempt rate as follows:
    
    






  
    β
    opt
    +
  
  =
  
    
      
        W
        0
      
      ⁡
      
        (
        
          
            
              
                C
                0
              
              ⁡
              
                (
                
                  σ
                  -
                  
                    T
                    c
                  
                
                )
              
            
            +
            
              
                C
                1
              
              ⁡
              
                (
                
                  
                    T
                    
                      s
                      0
                    
                  
                  -
                  
                    T
                    c
                  
                
                )
              
            
          
          
            eT
            c
          
        
        )
      
    
    +
    1
  




    
    
        where w0(·) is a Lambert W(z) function, W(z)eW(z)=z, for any complex number z,
        σ is a length of a media access control MAC slot,
        C0(1−β0)n0,
        C1n0β0(1−β0)n0−1,
        n0 is the total number of LP nodes,
        β0 is an average attempt rate per slot for each of the LP AC nodes,
        Tc is a mean time in slots for an unsuccessful transmission,
        Ts0 is a mean time in slots of a successful transmission for each of the AC 0 nodes,
        determining a maximum HP throughput Γ2(βopt+,r1,L,rI) by substituting the βopt+ for β+ in Γ2(β+,r1,L,rI) as follows:
    
    






  
    
      Γ
      2
    
    ⁡
    
      (
      
        
          β
          +
        
        ,
        
          r
          1
        
        ,
        L
        ,
        
          r
          l
        
      
      )
    
  
  =
  
    
      
        β
        +
      
      ⁢
      
        e
        
          -
          
            β
            +
          
        
      
      ⁢
      
        C
        0
      
      ⁢
      
        
          D
          0
        
        ⁡
        
          (
          
            
              r
              1
            
            ,
            L
            ,
            
              r
              l
            
          
          )
        
      
    
    
      
        
          
            
              
                e
                
                  -
                  
                    β
                    +
                  
                
              
              ⁢
              
                
                  C
                  0
                
                ⁡
                
                  [
                  
                    σ
                    +
                    
                      
                        β
                        +
                      
                      ⁢
                      
                        
                          D
                          1
                        
                        ⁡
                        
                          (
                          
                            
                              r
                              1
                            
                            ,
                            L
                            ,
                            
                              r
                              l
                            
                          
                          )
                        
                      
                    
                    +
                    
                      
                        C
                        1
                      
                      ⁢
                      
                        T
                        
                          s
                          0
                        
                      
                      ⁢
                      
                        /
                      
                      ⁢
                      
                        C
                        0
                      
                    
                  
                  ]
                
              
            
            +
          
        
      
      
        
          
            
              [
              
                1
                -
                
                  
                
                ⁢
                
                  
                    e
                    
                      -
                      
                        β
                        +
                      
                    
                  
                  ⁡
                  
                    (
                    
                      
                        C
                        0
                      
                      +
                      
                        C
                        1
                      
                      +
                      
                        
                          β
                          +
                        
                        ⁢
                        
                          C
                          0
                        
                      
                    
                    )
                  
                
              
              ]
            
            ⁢
            
              T
              c
            
          
        
      
    
  




    
    
        where Γ2(β+,r1,L,rI) is a maximum system throughput,
        β+ is a total asymptotic HP attempt rate,
    
    






  
    
      
        D
        0
      
      ⁡
      
        (
        
          
            r
            1
          
          ,
          L
          ,
          
            r
            I
          
        
        )
      
    
    =
    
      
        
          ∑
          
            i
            =
            1
          
          I
        
        ⁢
        
          r
          i
        
      
      
        
          ∑
          
            i
            =
            1
          
          I
        
        ⁢
        
          
            r
            i
          
          
            L
            i
          
        
      
    
  
  ,
  
    
      
        D
        1
      
      ⁡
      
        (
        
          
            r
            1
          
          ,
          L
          ,
          
            r
            I
          
        
        )
      
    
    =
    
      
        
          ∑
          
            i
            =
            1
          
          I
        
        ⁢
        
          
            
              r
              i
            
            
              L
              i
            
          
          ⁢
          
            T
            
              s
              i
            
          
        
      
      
        
          ∑
          
            i
            =
            1
          
          I
        
        ⁢
        
          
            r
            i
          
          
            L
            i
          
        
      
    
  
  ,




    
    
        ri is a ratio between HP AC i and AC 1 throughput,
        Tsi is a mean time in slots of a successful transmission for each of the AC i nodes,
        determining one or more numerical values of operating parameters including λi such that Σi=1IniλiLi<Γ2(βopt+,r1,L,rI)
        where Σi=1IniλiLi is a total traffic load of all of the HP ACs;
        determining the CWD based on a requirement of total delay of a flow.
    
    


Preferably, the total delay of the flow is related to the determined numerical values of λi.
Preferably, the total traffic load of all of the HP ACs is just slightly below the maximum HP throughput.



BRIEF DESCRIPTION OF THE DRAWINGS
Embodiments of the present invention are described in more detail hereinafter with reference to the drawings, in which:
FIG. 1 depicts admission region of AC 1 when n2=4, 10, 16, 22, 28 for homogeneous traffic;
FIG. 2 depicts total HP throughput for homogeneous traffic when n2=4, 10, 16, 22, 28 for homogeneous traffic;
FIG. 3 depicts admission region of AC 1 when n2=4, 10, 16 for heterogeneous traffic;
FIG. 4 depicts total HP throughput when n2=4, 10, 16 for heterogeneous traffic;
FIG. 5a depicts the average delay vs running time for AC 0; FIG. 5b depicts the average delay vs running time for AC 1; and FIG. 5c depicts the average delay vs running time for AC 2;
FIG. 6 depicts the average time between two consecutive successful transmission from all HP nodes; and
FIG. 7 depicts an overview of CoD hardware implementation for each node and the AP.



DETAILED DESCRIPTION
In the following description, methods for providing admission control to a wireless network system are set forth as preferred examples. It will be apparent to those skilled in the art that modifications, including additions and/or substitutions may be made without departing from the scope and spirit of the invention. Specific details may be omitted so as not to obscure the invention; however, the disclosure is written to enable one skilled in the art to practice the teachings herein without undue experimentation.
The present invention is devoted to designing a novel scheme for maximizing channel utilization and providing satisfactory QoS in wireless “hot-spot” networks. Then a theoretical model is developed to analyze the proposed scheme. Extensive simulations verify that the scheme is very effective and the theoretical result is very accurate.
Accordingly, Section A presents the proposed CoD scheme. Section B performs an asymptotic analysis for admission control. Section E validates the effectiveness of the proposed scheme and the accuracy of the theoretical result.
A. The Proposed CoD Scheme
In this section, one presents the proposed contention-on-demand (CoD) MAC scheme.
CoD is a parameterized quality-of-service (QoS)-oriented scheme. CoD employs the CSMA feature in IEEE 802.11 DCF, but replaces the BEB algorithm by one fixed-CW algorithm and one dynamic-CW algorithm. In the fixed-CW algorithm, the CW keeps unchanged for all flows that adopt this algorithm. The fixed-CW algorithm provides best-effort service for low-priority (LP) flows or elastic flows. In the dynamic-CW algorithm, the CW is dynamically adjusted by the AIMD algorithm, so as to adapt to QoS requirements of each flow that adopts this algorithm. The dynamic-CW algorithm provides QoS service for high-priority (HP) flows or inelastic flows. Once the CW is determined, each node chooses a backoff count that is uniformly distributed in the CW, and then performs backoff like 802.11 DCF.
In addition, CoD also employs an admission control unit to limit the total traffic load of HP flows. When the system resource is available, AIMD has the property of converging to equal values of the control variable [7]. The admission control is used to ensure that the system resource is available.
Let CWD denote current CW size. Below, one first presents the fixed-CW and dynamic-CW algorithms, and then present the admission control.
A.1. Fixed-CW Algorithm
The fixed-CW algorithm requires a pre-configured CW size denoted by CW0. In the algorithm, CWD=CW0=400.
A.2. Dynamic-CW Algorithm
The dynamic-CW algorithm adjusts the CW size according to the well-known AIMD algorithm. This algorithm requires five pre-configured parameters: CWinit (the initial CW size), CWmin (the minimum CW size), CWmax (the maximum CW size), fai (the additive-increase factor), and fmd (the multiplicative-decrease factor).










Algorithm1 Dynamic-CW (delayi, pktIntval)

















 1:  Init: CWinit=400; CWmin=10; CWmax=103; fai=10;



fmd=0.93; delayAVG=0; CWa=CWinit.



 2:  if Transmit the i-th packet then



 3:   CWD=CWinit if i=1; CWa otherwise.



 4:  end if



 5:  if Receive ACK of the i-th packet then



 6:   delayAVG+=(delayi−delayAVG)/i++.



 7:   if delayAVG<pktIntval then



 8:   CWa=CWa+fai.



 9:   else



 10:   CWa=CWa × fmd.



 11:  end if



 12:  CWa=max(min(CWa, CWmax), CWmin)



 13: end if










Algorithm 1 presents an example making the total delay converge to the packet arrival interval of a flow. In this way, one desires to guarantee the throughput stability (namely, the throughput of a flow is equal to its offered load). In this algorithm, one has two inputs: delay i  denoting the total delay of the i-th acknowledged packet of a flow, and pktIntval denoting the average packet arrival interval of the flow. When a flow is admitted to transmit its packets, CWD is set to CWinit for its first packet transmission and is set to CWa for the subsequent transmissions. CWa is set by the AIMD algorithm, in order to ensure that the total average delay, delayAVG, converges to pktIntval. To this end, one first calculates delayAVG (in line 6 of Algorithm 1), and then compare delayAVG with pktIntval. If delayAVG<pktIntval, one increases CWa (in line 8) and decrease CWa (in line 10) otherwise.
One now explains the calculation of delayAVG (in line 6). Let delayAVGn denote the average of first n total delays, namely




  
    delayAVG
    n
  
  =
  
    
      
        
          delay
          1
        
        +
        ...
        +
        
          delay
          n
        
      
      n
    
    .
  




Then, one has




  
    
      
        
          
            
              
                delayAVG
                
                  n
                  +
                  1
                
              
              =
              
              ⁢
              
                
                  
                    delay
                    1
                  
                  +
                  ...
                  +
                  
                    delay
                    
                      n
                      +
                      1
                    
                  
                
                
                  n
                  +
                  1
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    delay
                    
                      n
                      +
                      1
                    
                  
                  +
                  
                    n
                    ×
                    
                      delayAVG
                      n
                    
                  
                
                
                  n
                  +
                  1
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    delay
                    
                      n
                      +
                      1
                    
                  
                  +
                  
                    
                      (
                      
                        n
                        +
                        1
                      
                      )
                    
                    ×
                    
                      delayAVG
                      n
                    
                  
                  -
                  
                    delayAVG
                    n
                  
                
                
                  n
                  +
                  1
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  delayAVG
                  n
                
                +
                
                  
                    
                      
                        delay
                        
                          n
                          +
                          1
                        
                      
                      -
                      
                        delayAVG
                        n
                      
                    
                    
                      n
                      +
                      1
                    
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




Eq. (1) manifests that one only needs to store delayAVGn (instead of delay1, . . . , delayn), so that one can calculate delayAVGn+1 by delayn+1. The expression of delayAVG (in line 6) illustrates the computational form of (1) in programming.
In general, Algorithm 1 can make the total delay converge to the delay target of a flow. Therefore, one may replace the second input parameter in Algorithm 1, pktIntval, by delayTarget denoting the total delay requirement of a flow. When delayTarget is set to pktIntval, Algorithm 1 guarantees the throughput stability of a flow.
A.3. Admission Control
To ensure that the dynamic-CW algorithm can provide satisfactory QoS service for HP flows, one requires that the total traffic load, Σi=11niλiLi, should be below a upper bound Λ, namely,

Σi=1IniλiLi<Λ,

where one assumes that HP traffic is classified into I access categories (ACs), ni, λi, and Li, respectively, represent the node number, the packet arrival rate, and the packet size of AC i.
In the following sections, one seeks the optimal upper bound Λ.
B. Asymptotic HP Throughput
In this section, one performs an asymptotic analysis to derive the optimal upper bound Λ mentioned in Section A.3.
One assumes that the system, running the CoD schemes, consists of one LP AC and I HP ACs. All data packets from HP and LP nodes are transmitted to the AP, and the AP acts purely as the receiver of data packets.
The LP AC 0 has n0 nodes. Each LP node has the same packet size L0 and always generates a random backoff count uniformly distributed in [0,CW0] for each new transmission or retransmission, where CW0>1. one assumes that each LP node is in saturation operation (i.e., the node always has packets to transmit) because here one studies the maximum stable throughput that HP ACs can achieve, regardless of how the LP offered loads vary.
Each HP AC i , 1≦i≦I, has ni nodes, where n1+L+nI=n. Each HP i node has the same packet size Li and packet arrival rate λi, and always generates a random backoff count uniformly distributed in [0,CWD] for each new transmission or retransmission, where CWD is dynamically set by the dynamic-CW algorithm.
B.1 Exact HP Throughput
Let β0 be the average attempt rate per slot for each LP AC node on the condition that the buffer is not empty. From [12], one has β0=2/(CW0+1).
Let βi, i=1,L,I, be the average attempt rate per slot for each HP AC i node on the condition that all LP AC nodes will not transmit.
Let Ω be the mean time that elapses for one decrement of the backoff counter. The generic slot Ω takes different values depending on whether the slot is idle, interrupted by a successful transmission, or a collision. In terms of β0,β1,L,βI, one has




  
    
      
        
          
            Ω
            =
            
              
                
                  (
                  
                    1
                    -
                    
                      P
                      b
                    
                  
                  )
                
                ⁢
                σ
              
              +
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  I
                
                ⁢
                
                  
                
                ⁢
                
                  
                    P
                    
                      s
                      i
                    
                  
                  ⁢
                  
                    T
                    
                      s
                      i
                    
                  
                
              
              +
              
                
                  P
                  c
                
                ⁢
                
                  T
                  c
                
              
            
          
          ,
          
            

          
          ⁢
          where
        
        ⁢
        
          

        
        ⁢
        
          
            
              P
              b
            
            =
            
              1
              -
              
                
                  ∏
                  
                    j
                    =
                    0
                  
                  I
                
                ⁢
                
                  
                
                ⁢
                
                  
                    (
                    
                      1
                      -
                      
                        β
                        j
                      
                    
                    )
                  
                  
                    n
                    j
                  
                
              
            
          
          ,
          
            

          
          ⁢
          
            
              P
              
                s
                i
              
            
            =
            
              
                n
                i
              
              ⁢
              
                β
                i
              
              ⁢
              
                
                  
                    ∏
                    
                      j
                      =
                      0
                    
                    I
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      (
                      
                        1
                        -
                        
                          β
                          j
                        
                      
                      )
                    
                    
                      n
                      j
                    
                  
                
                
                  (
                  
                    1
                    -
                    
                      β
                      i
                    
                  
                  )
                
              
            
          
          ,
          
            

          
          ⁢
          
            
              P
              c
            
            =
            
              
                P
                b
              
              -
              
                
                  ∑
                  
                    j
                    =
                    0
                  
                  I
                
                ⁢
                
                  
                
                ⁢
                
                  
                    P
                    
                      s
                      j
                    
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




In (2), σ is the length of one MAC slot, i.e., σ=1 slot, Tsi is the mean time (in slots) of a successful transmission for each AC i node; Tc is the mean time (in slots) for an unsuccessful transmission; Pb is the probability of a busy slot; Psi is the probability of a successful packet transmission for each AC i node; and Pc is the probability of an unsuccessful packet transmission.
Exact HP Throughput: Let Γi(n,β1,L,βI) be the HP AC i throughput, which is defined to be the number of bits transmitted successfully by all HP AC i nodes in the time interval of Ω. Let Γ(n,β1,L,βI) be the total HP throughput. one has




  
    
      
        
          
            Γ
            i
          
          ⁡
          
            (
            
              n
              ,
              
                β
                1
              
              ,
              L
              ,
              
                β
                I
              
            
            )
          
        
        =
        
          
            
              
                
                  P
                  
                    s
                    i
                  
                
                ⁢
                
                  L
                  i
                
              
              Ω
            
            .
            
              

            
            ⁢
            
              Γ
              ⁡
              
                (
                
                  n
                  ,
                  
                    β
                    1
                  
                  ,
                  L
                  ,
                  
                    β
                    I
                  
                
                )
              
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              I
            
            ⁢
            
              
            
            ⁢
            
              
                
                  Γ
                  i
                
                ⁡
                
                  (
                  
                    n
                    ,
                    
                      β
                      1
                    
                    ,
                    L
                    ,
                    
                      β
                      I
                    
                  
                  )
                
              
              .
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




B.2. Asymptotic HP Throughput

One starts with the asymptotic assumption below.
Asymptotic assumption: To perform asymptotic analysis, one assumes that

n→∞ and βi+limn→∞niβi≦0   (4)

where βi+ is called the asymptotic attempt rate of HP AC i nodes.
Let us define

    
    
        C0(1−β0)n0,
        C1n0β0(1−β0)n0−1,
        β+Σi=1Iβi+
    
    


where β+ is called the total asymptotic HP attempt rate. Then, under assumption (4), applying Poisson approximation to (2), one has




  
    
      
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              
            
            ⁢
            
              P
              b
            
          
          =
          
            1
            -
            
              
                e
                
                  -
                  
                    β
                    +
                  
                
              
              ⁢
              
                C
                0
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              
            
            ⁢
            
              P
              
                s
                i
              
            
          
          =
          
            
              β
              i
              +
            
            ⁢
            
              e
              
                -
                
                  β
                  +
                
              
            
            ⁢
            
              C
              0
            
          
        
        ,
        
          
            for
            ⁢
            
              
            
            ⁢
            1
          
          <
          i
          <
          I
        
        ,
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                I
              
              ⁢
              
                
              
              ⁢
              
                P
                
                  s
                  i
                
              
            
          
          =
          
            
              β
              +
            
            ⁢
            
              e
              
                -
                
                  β
                  +
                
              
            
            ⁢
            
              C
              0
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              P
              
                s
                0
              
            
          
          =
          
            
              e
              
                -
                
                  β
                  +
                
              
            
            ⁢
            
              C
              1
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              P
              c
            
          
          =
          
            1
            -
            
              
                e
                
                  -
                  
                    β
                    +
                  
                
              
              ⁡
              
                (
                
                  
                    C
                    0
                  
                  +
                  
                    C
                    1
                  
                  +
                  
                    
                      β
                      +
                    
                    ⁢
                    
                      C
                      0
                    
                  
                
                )
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                I
              
              ⁢
              
                
                  P
                  
                    s
                    i
                  
                
                ⁢
                
                  T
                  
                    s
                    i
                  
                
              
            
          
          =
          
            
              e
              
                -
                
                  β
                  +
                
              
            
            ⁢
            
              C
              0
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                I
              
              ⁢
              
                
              
              ⁢
              
                
                  β
                  i
                  +
                
                ⁢
                
                  T
                  
                    s
                    i
                  
                
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            Ω
          
          =
          
            
              
                
                  
                    
                      e
                      
                        -
                        
                          β
                          +
                        
                      
                    
                    ⁢
                    
                      
                        C
                        0
                      
                      ⁡
                      
                        [
                        
                          σ
                          +
                          
                            
                              ∑
                              
                                i
                                =
                                1
                              
                              I
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              
                                β
                                i
                                +
                              
                              ⁢
                              
                                T
                                
                                  s
                                  i
                                
                              
                            
                          
                          +
                          
                            
                              C
                              1
                            
                            ⁢
                            
                              
                                T
                                
                                  s
                                  0
                                
                              
                              /
                              
                                C
                                0
                              
                            
                          
                        
                        ]
                      
                    
                  
                  +
                
              
            
            
              
                
                  
                    [
                    
                      1
                      -
                      
                        
                          e
                          
                            -
                            
                              β
                              +
                            
                          
                        
                        ⁡
                        
                          (
                          
                            
                              C
                              0
                            
                            +
                            
                              C
                              1
                            
                            +
                            
                              
                                β
                                +
                              
                              ⁢
                              
                                C
                                0
                              
                            
                          
                          )
                        
                      
                    
                    ]
                  
                  ⁢
                  
                    
                      T
                      c
                    
                    .
                  
                
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




Asymptotic HP throughput: Let




  
    Γ
    ⁡
    
      (
      
        
          β
          1
          +
        
        ,
        L
        ,
        
          β
          I
          +
        
      
      )
    
  
  ⁢
  
    
  
  ⁢
  
    
      lim
      
        n
        →
        ∞
      
    
    ⁢
    
      
    
    ⁢
    
      Γ
      ⁡
      
        (
        
          n
          ,
          
            β
            1
          
          ,
          L
          ,
          
            β
            I
          
        
        )
      
    
  




be the total asymptotic HP throughput. From (5) and (3), Γ(β1+,L,βI+), is given by





  
    
      
        Γ
        ⁡
        
          (
          
            
              β
              1
              +
            
            ,
            L
            ,
            
              β
              I
              +
            
          
          )
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        =
        
          
            
              e
              
                -
                
                  β
                  +
                
              
            
            ⁢
            
              C
              0
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  1
                
                I
              
              ⁢
              
                
              
              ⁢
              
                
                  β
                  i
                  +
                
                ⁢
                
                  L
                  i
                
              
            
          
          
            
              
                
                  
                    
                      e
                      
                        -
                        
                          β
                          +
                        
                      
                    
                    ⁢
                    
                      
                        C
                        0
                      
                      ⁡
                      
                        [
                        
                          σ
                          +
                          
                            
                              ∑
                              
                                i
                                =
                                1
                              
                              I
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              
                                β
                                i
                                +
                              
                              ⁢
                              
                                T
                                
                                  s
                                  i
                                
                              
                            
                          
                          +
                          
                            
                              C
                              1
                            
                            ⁢
                            
                              
                                T
                                
                                  s
                                  0
                                
                              
                              /
                              
                                C
                                0
                              
                            
                          
                        
                        ]
                      
                    
                  
                  +
                
              
            
            
              
                
                  
                    [
                    
                      1
                      -
                      
                        
                          e
                          
                            -
                            
                              β
                              +
                            
                          
                        
                        ⁡
                        
                          (
                          
                            
                              C
                              0
                            
                            +
                            
                              C
                              1
                            
                            +
                            
                              
                                β
                                +
                              
                              ⁢
                              
                                C
                                0
                              
                            
                          
                          )
                        
                      
                    
                    ]
                  
                  ⁢
                  
                    T
                    c
                  
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




where Li is the packet size for AC i node.
C. Admission Control for Homogeneous Traffic
In this section, one focuses on homogeneous traffic, where all nodes have the same packet size L, namely L0=L1= . . . =LI≡L. Below, one first optimizes the total asymptotic HP attempt rate β+, and then specify the admission control rule.
C.1. Optimal Asymptotic HP Attempt Rate
For homogeneous traffic, one can further assume Tc=Tsi. This assumption states that the ACK timeout matches one successful transmission time. It has been implemented in NS2 [13] and has widely been used in previous works such as [14]. Then, Ω in (2) and




  
    lim
    
      n
      →
      ∞
    
  
  ⁢
  
    
  
  ⁢
  Ω




in (5) reduce to





  
    
      
        
          Ω
          =
          
            
              
                (
                
                  1
                  -
                  
                    P
                    b
                  
                
                )
              
              ⁢
              σ
            
            +
            
              
                P
                b
              
              ⁢
              
                T
                c
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              lim
              
                n
                →
                ∞
              
            
            ⁢
            
              
            
            ⁢
            Ω
          
          =
          
            
              T
              c
            
            +
            
              
                
                  e
                  
                    -
                    
                      β
                      +
                    
                  
                
                ⁡
                
                  [
                  
                    σ
                    -
                    
                      T
                      c
                    
                  
                  ]
                
              
              ⁢
              
                
                  C
                  0
                
                .
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Therefore, Γ(β1+,L,βI+) in (6) can be expressed in terms of β+. Let Γ1(β+)Γ(β1+,L,βI+). one has




  
    
      
        
          
            Γ
            1
          
          ⁡
          
            (
            
              β
              +
            
            )
          
        
        =
        
          
            
              
                e
                
                  -
                  
                    β
                    +
                  
                
              
              ⁢
              
                β
                +
              
              ⁢
              
                C
                0
              
              ⁢
              L
            
            
              
                T
                c
              
              +
              
                
                  
                    e
                    
                      -
                      
                        β
                        +
                      
                    
                  
                  ⁡
                  
                    [
                    
                      σ
                      -
                      
                        T
                        c
                      
                    
                    ]
                  
                
                ⁢
                
                  C
                  0
                
              
            
          
          .
        
      
    
    
      
        (
        9
        )
      
    
  




Let βopt+ represent the optimal β+ that maximizes the asymptotic total HP throughput (9). Theorem 1 explicitly expresses βopt+.
Theorem 1. Under assumption (4), if Tc=Tsi, βopt+ is




  
    
      
        
          
            β
            opt
            +
          
          =
          
            
              
                W
                0
              
              ⁡
              
                (
                
                  
                    
                      C
                      0
                    
                    ⁡
                    
                      (
                      
                        σ
                        -
                        
                          T
                          c
                        
                      
                      )
                    
                  
                  
                    e
                    ⁢
                    
                      
                    
                    ⁢
                    
                      T
                      c
                    
                  
                
                )
              
            
            +
            1
          
        
        ,
      
    
    
      
        (
        10
        )
      
    
  




where W0(·) is one branch of the Lambert W(z) function [15], W(z)eW(z)=z, for any complex number z.
Proof To maximize Γ1(β+), one sets the first derivative of Γ1(β+) in (9) with respect to β+ to zero. This leads to

Tcβ+=C0(σ−Tc)e−β++Tc 

Tc(β+−1)=C0(σ−Tc)e−β+

(β+−1)e(β+−1)=C0(σ−Tc)e−1/Tc.

Then β+−1=W0(C0(σ−Tc)/e−1/Tc) or W−1(C0(σ−Tc)e−1/Tc). one has βopt+=W0(C0(σ−Tc)e−1/Tc)+1≧0, since only W0(C0(σ−Tc)e−1/Tc)>−1 for C0(σ−Tc)e−1/Tc∈(−1/e, 0).
When n0=0 and I=1 (i.e., the system has 1 HP AC and has not LP nodes), βopt+ reduces to the solution to (10) in [10].
C.2. Admission Control
To provide satisfactory QoS, the total traffic load of all HP ACs, Σi=1IniλiL, should be below the optimal asymptotic HP throughput Γ1(βopt+), namely,

βi=1IniλiL<Γ1(βopt+),   (11)

where Γ1(·) is given in (9).
From (9)-(10), one knows that for homogeneous traffic, the maximum total HP throughput keeps unchanged, regardless of how the node number and the packet arrival rate of HP nodes vary; therefore, one only needs to compute the total offered load for admission control, as shown in (11).
D. Admission Control for Heterogeneous Traffic
In this section, one focuses on heterogeneous traffic, where different ACs have different packet sizes. Below, one first optimizes the total asymptotic HP attempt rate β+, and then specify the admission control rule.
D.1. Optimal Asymptotic Aggregate Attempt Rate
Let βopt+ represent the optimal β+ that maximizes the asymptotic total HP throughput (6). Theorem 2 explicitly expresses βopt+.
Theorem 2. Under assumption (4), if




  
    
      
        
          
            
              [
              
                
                  
                    C
                    0
                  
                  ⁡
                  
                    (
                    
                      σ
                      -
                      
                        T
                        c
                      
                    
                    )
                  
                
                +
                
                  
                    C
                    1
                  
                  ⁡
                  
                    (
                    
                      
                        T
                        
                          s
                          0
                        
                      
                      -
                      
                        T
                        c
                      
                    
                    )
                  
                
              
              ]
            
            /
            
              T
              c
            
          
          ∈
          
            (
            
              
                -
                1
              
              ,
              0
            
            )
          
        
        ,
        
          

        
        ⁢
        
          
            
              β
              opt
              +
            
            ⁢
            
              
            
            ⁢
            is
            ⁢
            
              
            
            ⁢
            
              β
              opt
              +
            
          
          =
          
            
              
                W
                0
              
              ⁡
              
                (
                
                  
                    
                      
                        C
                        0
                      
                      ⁡
                      
                        (
                        
                          σ
                          -
                          
                            T
                            c
                          
                        
                        )
                      
                    
                    +
                    
                      
                        C
                        1
                      
                      ⁡
                      
                        (
                        
                          
                            T
                            
                              s
                              0
                            
                          
                          -
                          
                            T
                            c
                          
                        
                        )
                      
                    
                  
                  
                    e
                    ⁢
                    
                      
                    
                    ⁢
                    
                      T
                      c
                    
                  
                
                )
              
            
            +
            1.
          
        
      
    
    
      
        (
        12
        )
      
    
  




where W0(·) is one branch of the Lambert W(z) function [15], W(z)eW(z)=z, for any complex number z.
(i) Theorem 2 manifests that the optimal asymptotic HP attempt rate is independent of the HP traffic characteristics, because βopt+ in (12) only depends on the LP transmission time Ts0 and the common collision time Tc. (ii) Theorem 1 is a special case of Theorem 2, because (12) reduces to (10) if one sets Ts0=Tc.
D.2. Admission Control
In this section, for heterogeneous traffic, one develops an admission control rule. Before this, one first calculates the maximum system throughput.
From (6) and Theorem 2, the maximum system throughput depends on not only the optimal total HP attempt rate but also the HP packet size. To express the maximum system throughput in terms of βopt+, one adopts a key approximation, βi<<1, which is widely used in the related literatures such as [16][17]. The approximation holds true since βi represents the per-node attempt rate in a very short slot (e.g., 1 slot=20 μs in 802.11b) and therefore it is generally much less than 1.
Let




  
    r
    i
  
  ⁢
  
    =
    △
  
  ⁢
  
    
      lim
      
        n
        →
        ∞
      
    
    ⁢
    
      
        r
        i
      
      
        Γ
        1
      
    
  




be the ratio between HP AC i and AC 1 throughput. With the approximation βi<<1, and the assumption (4), one has





  
    
      
        
          r
          i
        
        =
        
        ⁢
        
          
            lim
            
              n
              →
              ∞
            
          
          ⁢
          
            
              Γ
              i
            
            
              Γ
              1
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            lim
            
              n
              →
              ∞
            
          
          ⁢
          
            
              
                P
                
                  s
                  i
                
              
              ⁢
              
                L
                i
              
            
            
              
                P
                
                  s
                  1
                
              
              ⁢
              
                L
                1
              
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            lim
            
              n
              →
              ∞
            
          
          ⁢
          
            [
            
              
                
                  
                    n
                    i
                  
                  ⁢
                  
                    β
                    i
                  
                
                
                  
                    n
                    1
                  
                  ⁢
                  
                    β
                    1
                  
                
              
              ⁢
              
                
                  1
                  -
                  
                    β
                    1
                  
                
                
                  1
                  -
                  
                    β
                    i
                  
                
              
              ⁢
              
                
                  L
                  i
                
                
                  L
                  1
                
              
            
            ]
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            
              β
              i
              +
            
            
              β
              1
              +
            
          
          ⁢
          
            
              
                L
                i
              
              
                L
                1
              
            
            .
          
        
      
    
  




Then, one has




  
    
      β
      i
      +
    
    =
    
      
        
          r
          i
        
        ⁢
        
          L
          1
        
        ⁢
        
          β
          1
          +
        
      
      
        L
        i
      
    
  
  ,
  
    

  
  ⁢
  
    
      β
      +
    
    =
    
      
        
          ∑
          
            i
            =
            1
          
          I
        
        ⁢
        
          
        
        ⁢
        
          β
          i
          +
        
      
      =
      
        
          ∑
          
            i
            =
            1
          
          I
        
        ⁢
        
          
        
        ⁢
        
          
            
              r
              i
            
            
              L
              i
            
          
          ⁢
          
            L
            1
          
          ⁢
          
            β
            1
            +
          
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        L
        1
      
      ⁢
      
        β
        1
        +
      
    
    =
    
      
        
          β
          +
        
        
          
            ∑
            
              i
              =
              1
            
            I
          
          ⁢
          
            
          
          ⁢
          
            
              r
              i
            
            
              L
              i
            
          
        
      
      .
    
  




Hence, βi+ can be expressed in terms of β+.




  
    
      
        
          β
          i
          +
        
        =
        
          
            
              
                r
                i
              
              ⁢
              
                L
                1
              
              ⁢
              
                β
                1
                +
              
            
            
              L
              i
            
          
          =
          
            
              
                r
                i
              
              
                L
                i
              
            
            ⁢
            
              
                
                  β
                  +
                
                
                  
                    Σ
                    
                      i
                      =
                      1
                    
                    I
                  
                  ⁢
                  
                    
                      r
                      i
                    
                    
                      L
                      i
                    
                  
                
              
              .
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




Let us define D0(r1,L,rI) and D1(r1,L,rI) as follows.




  
    
      
        
          
            
              D
              0
            
            ⁡
            
              (
              
                
                  r
                  1
                
                ,
                L
                ,
                
                  r
                  I
                
              
              )
            
          
          =
          
            
              
                Σ
                
                  i
                  =
                  1
                
                I
              
              ⁢
              
                r
                i
              
            
            
              
                Σ
                
                  i
                  =
                  1
                
                I
              
              ⁢
              
                
                  r
                  i
                
                
                  L
                  i
                
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            
              D
              1
            
            ⁡
            
              (
              
                
                  r
                  1
                
                ,
                L
                ,
                
                  r
                  I
                
              
              )
            
          
          =
          
            
              
                
                  Σ
                  
                    i
                    =
                    1
                  
                  I
                
                ⁢
                
                  
                    r
                    i
                  
                  
                    L
                    i
                  
                
                ⁢
                
                  T
                  
                    s
                    i
                  
                
              
              
                
                  Σ
                  
                    i
                    =
                    1
                  
                  I
                
                ⁢
                
                  
                    r
                    i
                  
                  
                    L
                    i
                  
                
              
            
            .
          
        
      
    
    
      
        (
        14
        )
      
    
  




Substituting (13) and (14) into (6), one has Γ(β1+,L,βI+)≈Γ2(β+,{ri}), where




  
    
      
        
          
            Γ
            2
          
          ⁡
          
            (
            
              
                β
                +
              
              ,
              
                r
                1
              
              ,
              L
              ,
              
                r
                I
              
            
            )
          
        
        =
        
          
            
              
                β
                +
              
              ⁢
              
                e
                
                  -
                  
                    β
                    +
                  
                
              
              ⁢
              
                C
                0
              
              ⁢
              
                
                  D
                  0
                
                ⁡
                
                  (
                  
                    
                      r
                      1
                    
                    ,
                    L
                    ,
                    
                      r
                      I
                    
                  
                  )
                
              
            
            
              
                
                  
                    
                      
                        e
                        
                          -
                          
                            β
                            +
                          
                        
                      
                      ⁢
                      
                        
                          C
                          0
                        
                        ⁡
                        
                          [
                          
                            σ
                            +
                            
                              
                                β
                                +
                              
                              ⁢
                              
                                
                                  D
                                  1
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      r
                                      1
                                    
                                    ,
                                    L
                                    ,
                                    
                                      r
                                      I
                                    
                                  
                                  )
                                
                              
                            
                            +
                            
                              
                                C
                                1
                              
                              ⁢
                              
                                
                                  T
                                  
                                    s
                                    0
                                  
                                
                                /
                                
                                  C
                                  0
                                
                              
                            
                          
                          ]
                        
                      
                    
                    +
                  
                
              
              
                
                  
                    
                      [
                      
                        1
                        -
                        
                          
                            e
                            
                              -
                              
                                β
                                +
                              
                            
                          
                          ⁢
                          
                            (
                            
                              
                                C
                                0
                              
                              +
                              
                                C
                                1
                              
                              +
                              
                                
                                  β
                                  +
                                
                                ⁢
                                
                                  C
                                  0
                                
                              
                            
                            )
                          
                        
                      
                      ]
                    
                    ⁢
                    
                      T
                      c
                    
                  
                
              
            
          
          .
        
      
    
    
      
        (
        15
        )
      
    
  




CAC rule: Let Λi=niλiLi denote the total traffic load of HP AC i. For admission control, one will ensure that the system throughput Γi of HP AC i is equal to its total offered load Λi, namely, Γi=Λi. Then, one can approximately calculate ri as follows.




  
    
      
        
          
            r
            i
          
          ≈
          
            
              Γ
              i
            
            
              Γ
              1
            
          
        
        =
        
          
            
              Λ
              i
            
            
              Λ
              1
            
          
          =
          
            
              
                
                  n
                  i
                
                ⁢
                
                  λ
                  i
                
                ⁢
                
                  L
                  i
                
              
              
                
                  n
                  1
                
                ⁢
                
                  λ
                  1
                
                ⁢
                
                  L
                  1
                
              
            
            .
          
        
      
    
    
      
        (
        16
        )
      
    
  




To provide satisfactory QoS, given Ci, Li, and Di, the total traffic load of all HP ACs, Σi=1IniλiL, should be below the maximum HP throughput Γ2(βopt+,r1,L,rI), namely,

Σi=1IniλiLi<Γ2(βopt+,r1,L,rI).   (17)

From the CAC rule, one has the following observations.
1 From (15), one knows that for heterogeneous traffic, the maximum total HP throughput varies as the ratio ri between offered loads varies; therefore, one needs to compute both the total offered load and the maximum total HP throughput for admission control, as shown in (17). This is a striking difference from the case of homogeneous traffic.

    
    
        The CAC rule (17) is applicable for homogeneous traffic, because (17) reduces to (11) if each AC has the same packet size and Tc=Tsi.
        The CAC rule (17) is also applicable for homogeneous IEEE 802.11 DCF networks. The reasons are: the CoD system becomes the homogeneous IEEE 802.11 DCF network, if n0=0 and I=1 (hence C0=1, C1=0), Li≡L, and the only AC adopts the binary-exponential-backoff algorithm. As a result, one has
    
    






  
    β
    opt
    +
  
  =
  
    ln
    ⁢
    
      b
      
        b
        -
        1
      
    
  




[18][6]; (15) reduces to (22) in [6]; and the CAC rule (17) becomes the CAC rule in DCF (i.e., Table II in [6]).

E. Model Verification








TABLE 1





Parameters for 802.11 b basic mode.



















CW0
32
Header
243 μs =
Mheader + Pheader +






RouteHeader


m/M
5/7
Ts
=
Header + Ltm + SIFS + δ +






ACK + δ + DIFS












σ
1
slot
Ts
=
Ts


δ
0
μs
Ltm
=
L bytes @ Rdata


SIFS
10
μs
ACK
304 μs =
24 bytes @ Rbasic + 14







bytes @ Rbasic


DIFS
50
μs
Mheader
 22 μs =
26 bytes @ Rdata + 4 bytes







@ Rdata


Rdata
11
Mbps
Pheader
192 μs =
24 bytes @ Rbasic


Rbasic
1
Mbps
RouteHeader
 29 μs =
40 bytes @ Rdata









In this section, one demonstrates the effectiveness of the proposed CAC scheme. one uses the 802.11 simulator in ns2 version 2.28 [13] as a validation tool, and set the protocol parameters to the default values for 802.11b, as listed in Table 1, where a slot is equal to 20 μs and δ denotes the propagation delay. In the simulation, one uses the Dumb Agent routing protocol. Each simulation value is an average over 5 simulation runs, where each run was for 100 seconds. A buffer size of 1000 packets is used in the simulation to mimic an infinite buffer.
In the experiment, one considers 4 ACs: AC 0, AC 1, AC 2, AC 3. one runs simulations for homogeneous traffic (where [L0, L1, L2, L3]=[500,500,500,500]) and heterogeneous traffic (where [L0, L1, L2, L3]=[500,80,400,800]). The other parameter settings are as follows: [n0,n1,n2,n3]=[10,*,*, 5] and [λ0,λ1,λ2,λ3]=[*,40,20,10] (unit: packets per second), where * represents the adjustable parameter value. In simulation, one considers two cases. In the first case, AC 0 nodes are in saturation operation while other nodes are in non-saturation operations, where one sets λ0=400 to mimic the saturation operation and label the simulation results with “sim_AC0_sat”. In the second case, all nodes are in non-saturation operation, where one sets λ0=40 for the non-saturation operation and label the simulation results with “sim_AC0_nsat”.
E.1. Homogeneous Tragic
In this section, one considers homogeneous traffic.
FIG. 1 plots the admission region of AC 1 when n2=4,10,16,22,28. In this experiment, one changes n2 and then find the maximum allowable n1. For the theoretical result, n1 is calculated by (11). For the simulation result, if accepting a new AC 1 node will cause the total throughput of AC 1 nodes to be less than their total offered loads, n1 is set to the current number of AC 1 nodes excluding the new AC 1 node. From this figure, one can see that n1 decreases as n2 increases. One has the following observations.
When AC 0 is in saturation operation, the theoretical results slightly overestimate the corresponding simulation results. Concretely speaking, each simulation result is less than 1 than the corresponding theoretical result for n2=4, 16, 22, 28, and the former is less than 2 than the latter for n2=10.
When AC 0 is in non-saturation operation, the theoretical results slightly underestimate the corresponding simulation results. Concretely speaking, each simulation result is more than 1 than the corresponding theoretical result for n2=4, 10, 16, 22, and the two results are equal for n2=28.
The above observations imply that the maximum allowable node number is sensitive to the traffic regime (i.e., saturation or non-saturation). One explanation is: one uses the basic mode to resolve collision for homogeneous traffic, so that collisions occur during packet transmission and therefore the overhead is large if the packet size is long; as a result, the overall overhead when the offered load approaches the capacity of the system but the system is still in non-saturation, might be obviously less than that when the system is in saturation operation.
FIG. 2 plots the total HP throughput when n2=4, 10, 16, 22, 28, where the theoretical result is calculated by (15). From this figure, one can see that the maximum total HP throughput remains unchanged, regardless of how n2 varies. For the simulation results, one has the following observations.
When AC 0 is in saturation operation, the maximum allowable traffic loads are slightly less than the corresponding theoretical result, because the maximum allowable n1 is less than the corresponding theoretical result (as shown in FIG. 1).
When AC 0 is in non-saturation operation, the maximum allowable traffic loads slightly fluctuate around the corresponding theoretical result, because the maximum allowable n1 is slightly larger than or equal to the corresponding theoretical result (as shown in FIG. 1).
E.2. Heterogeneous Traffic
In this section, one considers heterogeneous traffic.
FIG. 3 plots the admission region of AC 1 when n2=4, 10, 16. In this experiment, one changes n2 and then find the maximum allowable n1. For the theoretical result, n1 is calculated by (16). For the simulation result, if accepting a new AC 1 node will cause the total throughput of AC 1 nodes to be less than their total offered loads, n1 is set to the current number of AC 1 nodes excluding the new AC 1 node. From this figure, one can see that n1 decreases as n2 increases. One has the following observations.
Whether AC 0 is in saturation operation or in non-saturation operation, both simulation results are the same for each value of n2; and the theoretical results slightly underestimate the corresponding simulation results. It implies that the maximum allowable node number is insensitive to the traffic regime (i.e., saturation or non-saturation). This is a distinct difference from the case of homogeneous traffic. One explanation is: one uses the RTS/CTS mechanism to resolve collision for heterogeneous traffic, and hence the RTS/CTS overhead is fixed regardless of the packet size; as a result, the overall overhead when the offered load approaches the capacity of the system but the system is still in non-saturation, is almost the same as that when the system is in saturation operation.
FIG. 4 plots the total HP throughput when n2=4, 10, 16, where the theoretical result is calculated by (9). From this figure, one can see that for the theoretical result, the maximum total HP throughput increases as n2 increases. The reason is that AC 2 has a large packet size than AC 1 and therefore increasing the number of AC 2 will increase the throughput. This also implies that the system capacity varies with the heterogeneous traffic characteristics (i.e., the packet size, the node number, and the packet arrival rate). This is another distinct difference from the case of homogeneous traffic. In addition, one also observes that (i) when AC 0 is in saturation and non-saturation operations, both simulation values of the total HP throughput are almost the same for each value of n2, because both simulation values of the maximum allowable n1 are equal, as shown in FIG. 3; (ii) the theoretical results are slightly larger than the corresponding simulation results.
E.3. Convergence Speed
In this section, one demonstrates the convergence speed of the CoD scheme. One considers heterogeneous traffic, when n2=10, n1=7, and all other parameter settings are the same as those in FIG. 3. In the simulation, the targets of the total delay, respectively, are 25 ms, 50 ms, and 100 ms for AC 0, AC 1, AC 2. For each AC i (i=0; 1; 2), one is concerned with how fast that the CoD scheme makes the average total delay converge to its target value.
FIG. 5(a)-(c), respectively, plot the average total delay vs. the running time for AC 0, AC 1, and AC 2, where the average total delay is calculated by (1). From the three subfigures, one can see that the average delay of each AC first increases for a while, then decreases quickly, and finally converges to its target delay after 60 seconds. The reasons are explained as follows. At the beginning, each AC node has the same intimal CW value of 400, which are too large for the target delay. Then all AC nodes simultaneously start to decrease their CWs. When all CW sizes become small, they will cause many packet collisions, leading to large delays. After that, the CoD scheme makes all ACs adjust their respective CWs according to their respective delay targets. As a result, the average delay starts to decrease and finally converges to the target value.
According to an embodiment of the present invention, an FPGA-based implementation is specified for CoD. 802.11 MAC protocol consists of (a) CSMA/CA and (b) the binary exponential algorithm. The MAC design keeps (a) unchanged but replaces (b) by CoD. FIG. 7 illustrates the implementation framework for each node and the AP, which includes off-the-shelf components (i.e., 802.11 PHY, 802.11 LLC, and 802.11 MAC supporting CSMA/CA), and the newly proposed CoD module.
As illustrated below, one finishes the proof of Theorem 2.
Proof of Theorem 2: one first defines four types of virtual slots: (i) an idle MAC slot with probability 1−Pb, (ii) one successful transmission time of an LP node with probability Ps0, one collision time with probability Pc, and one successful transmission time of an HP node with probability PshΣi=1IPsi, where 1−Pb+Ps0+Pc+Psh=1.
Let To denote the interval between when one successful transmission from HP nodes ends and when the next successful transmission from HP nodes begins, as shown in FIG. 6. Clearly, To contains the first three types of virtual slots.
Let To denote the mean of To. To maximize the effective bandwidth occupied by HP nodes, one should minimize To.
Below, one first expresses To, and then find the optimal β+ that minimizes To.
Let Xo denote the number of the virtual slots during To. Since a virtual slot during To appears with probability 1−Psh, Xo follows a geometric distribution with parameter Psh and therefore its mean Xo=1/Psh−1.
Let Ωo be a random variable representing the length of a virtual slot during To. Ωo takes three types of values depending on the thee types of virtual slots during To. In terms of Pb, Ps0, and Pc, one defines Ωo as follows:




  
    
      
        
          Ω
          o
        
        =
        
          {
          
            
              
                σ
              
              
                
                  w
                  .
                  p
                  .
                
              
              
                
                  
                    
                      1
                      -
                      
                        P
                        b
                      
                    
                    
                      1
                      -
                      
                        P
                        b
                      
                      +
                      
                        P
                        
                          s
                          0
                        
                      
                      +
                      
                        P
                        c
                      
                    
                  
                  ,
                
              
            
            
              
                
                  T
                  
                    S
                    0
                  
                
              
              
                
                  w
                  .
                  p
                  .
                
              
              
                
                  
                    
                      P
                      
                        s
                        0
                      
                    
                    
                      1
                      -
                      
                        P
                        b
                      
                      +
                      
                        P
                        
                          s
                          0
                        
                      
                      +
                      
                        P
                        c
                      
                    
                  
                  ,
                
              
            
            
              
                
                  T
                  c
                
              
              
                
                  w
                  .
                  p
                  .
                
              
              
                
                  
                    
                      P
                      c
                    
                    
                      1
                      -
                      
                        P
                        b
                      
                      +
                      
                        P
                        
                          s
                          0
                        
                      
                      +
                      
                        P
                        c
                      
                    
                  
                  ,
                
              
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




where Ts0, Tc, Pb, Ps0, Pc are defined in (2). Then the mean Ωo can be easily calculated by (18).
Then To is equal to the mean time of a virtual slot, Ωo, times the mean number of virtual slots Xo. Noting that 1−Pb+Ps0+Pc+Psh=1, one has




  
    
      
        
          
            
              
                
                  T
                  _
                
                o
              
              =
              
              ⁢
              
                
                  
                    Ω
                    _
                  
                  o
                
                ·
                
                  
                    X
                    _
                  
                  o
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    
                      σ
                      ⁡
                      
                        (
                        
                          1
                          -
                          
                            P
                            b
                          
                        
                        )
                      
                    
                    +
                    
                      
                        T
                        
                          s
                          0
                        
                      
                      ⁢
                      
                        P
                        
                          s
                          0
                        
                      
                    
                    +
                    
                      
                        T
                        c
                      
                      ⁢
                      
                        P
                        c
                      
                    
                  
                  
                    1
                    -
                    
                      P
                      b
                    
                    +
                    
                      P
                      
                        s
                        0
                      
                    
                    +
                    
                      P
                      c
                    
                  
                
                ·
                
                  
                    1
                    -
                    
                      P
                      
                        s
                        h
                      
                    
                  
                  
                    P
                    
                      s
                      h
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    σ
                    ⁡
                    
                      (
                      
                        1
                        -
                        
                          P
                          b
                        
                      
                      )
                    
                  
                  +
                  
                    
                      T
                      
                        s
                        0
                      
                    
                    ⁢
                    
                      P
                      
                        s
                        0
                      
                    
                  
                  +
                  
                    
                      T
                      c
                    
                    ⁢
                    
                      P
                      c
                    
                  
                
                
                  P
                  
                    s
                    h
                  
                
              
            
          
        
      
    
    
      
        (
        19
        )
      
    
  




Define T0+(β+)limn→∞T0+. Let C2σ−Tc+(Ts0−Tc)C1/C0. Taking n→∞ for both sides in (19) and applying (5), one has




  
    
      
        
          
            
              T
              _
            
            o
            +
          
          ⁡
          
            (
            
              β
              +
            
            )
          
        
        =
        
        ⁢
        
          
            
              
                
                  
                    σ
                    ⁡
                    
                      [
                      
                        
                          e
                          
                            -
                            
                              β
                              +
                            
                          
                        
                        ⁢
                        
                          C
                          0
                        
                      
                      ]
                    
                  
                  +
                  
                    
                      T
                      
                        s
                        0
                      
                    
                    ⁡
                    
                      [
                      
                        
                          e
                          
                            -
                            
                              β
                              +
                            
                          
                        
                        ⁢
                        
                          C
                          1
                        
                      
                      ]
                    
                  
                  +
                
              
            
            
              
                
                  
                    T
                    c
                  
                  ⁡
                  
                    [
                    
                      1
                      -
                      
                        
                          e
                          
                            -
                            
                              β
                              +
                            
                          
                        
                        ⁡
                        
                          (
                          
                            
                              C
                              0
                            
                            +
                            
                              C
                              1
                            
                            +
                            
                              
                                β
                                +
                              
                              ⁢
                              
                                C
                                0
                              
                            
                          
                          )
                        
                      
                    
                    ]
                  
                
              
            
          
          
            
              β
              +
            
            ⁢
            
              e
              
                -
                
                  β
                  +
                
              
            
            ⁢
            
              C
              0
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            σ
            +
            
              
                T
                
                  s
                  0
                
              
              ⁢
              
                
                  C
                  1
                
                
                  C
                  0
                
              
            
            +
            
              
                T
                c
              
              ⁡
              
                [
                
                  
                    
                      e
                      
                        β
                        +
                      
                    
                    
                      C
                      0
                    
                  
                  -
                  
                    (
                    
                      1
                      +
                      
                        
                          C
                          1
                        
                        
                          C
                          0
                        
                      
                      +
                      
                        β
                        +
                      
                    
                    )
                  
                
                ]
              
            
          
          
            β
            +
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            
              [
              
                σ
                -
                
                  T
                  c
                
                +
                
                  
                    (
                    
                      
                        T
                        
                          s
                          0
                        
                      
                      -
                      
                        T
                        c
                      
                    
                    )
                  
                  ⁢
                  
                    
                      C
                      1
                    
                    
                      C
                      0
                    
                  
                
              
              ]
            
            +
            
              
                T
                c
              
              ⁡
              
                [
                
                  
                    
                      e
                      
                        β
                        +
                      
                    
                    
                      C
                      0
                    
                  
                  -
                  
                    β
                    +
                  
                
                ]
              
            
          
          
            β
            +
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            
              C
              2
            
            +
            
              
                T
                c
              
              ⁡
              
                [
                
                  
                    
                      e
                      
                        β
                        +
                      
                    
                    
                      C
                      0
                    
                  
                  -
                  
                    β
                    +
                  
                
                ]
              
            
          
          
            β
            +
          
        
      
    
  




To minimize To+(β+), one sets the first derivative of To+(β+) with respect to β+ to zero. This leads to



  
    
      
        
          T
          c
        
        ⁡
        
          (
          
            
              
                e
                
                  β
                  +
                
              
              
                C
                0
              
            
            -
            1
          
          )
        
      
      ⁢
      
        β
        +
      
    
    =
    
      
        C
        2
      
      +
      
        
          T
          c
        
        ⁡
        
          (
          
            
              
                e
                
                  β
                  +
                
              
              
                C
                0
              
            
            -
            
              β
              +
            
          
          )
        
      
    
  




  
    
      
        T
        
          
            
          
          ⁢
          c
        
      
      ⁢
      
        
          e
          
            β
            +
          
        
        
          C
          0
        
      
      ⁢
      
        β
        +
      
    
    =
    
      
        C
        2
      
      +
      
        
          T
          c
        
        ⁢
        
          
            e
            
              β
              +
            
          
          
            C
            0
          
        
      
    
  




  
    
      
        T
        c
      
      ⁢
      
        β
        +
      
    
    =
    
      
        
          C
          0
        
        ⁢
        
          C
          2
        
        ⁢
        
          e
          
            -
            
              β
              +
            
          
        
      
      +
      
        T
        c
      
    
  




  
    
      
        T
        c
      
      ⁡
      
        (
        
          
            β
            +
          
          -
          1
        
        )
      
    
    =
    
      
        
          C
          0
        
        ⁢
        
          C
          2
        
        ⁢
        
          
            e
            
              -
              
                β
                +
              
            
          
          ⁢
          
            

          
          (
          
            
              β
              +
            
            -
            1
          
          )
        
        ⁢
        
          e
          
            (
            
              
                β
                +
              
              -
              1
            
            )
          
        
      
      =
      
        
          C
          0
        
        ⁢
        
          C
          2
        
        ⁢
        
          
            e
            
              -
              1
            
          
          /
          
            
              T
              c
            
            .
          
        
      
    
  



Then β+−1=W0(C0C2e−1/Tc) or W−1(C0C2e−1/Tc). One has βopt+=W0(C0C2e−1/Tc)+1≧0, since only W0(C0C2e−1/Tc)>−1 for C0C2e−1/Tc∈(−1/e, 0).
As a result, one obtains (12) since




  
    
      
        C
        0
      
      ⁢
      
        C
        2
      
    
    
      e
      ⁢
      
        
      
      ⁢
      
        T
        c
      
    
  
  =
  
    
      
        
          
            C
            0
          
          ⁡
          
            (
            
              σ
              -
              
                T
                c
              
            
            )
          
        
        +
        
          
            C
            1
          
          ⁡
          
            (
            
              
                T
                
                  s
                  0
                
              
              -
              
                T
                c
              
            
            )
          
        
      
      
        e
        ⁢
        
          
        
        ⁢
        
          T
          c
        
      
    
    .
  




The embodiments disclosed herein may be implemented using general purpose or specialized computing devices, computer processors, or electronic circuitries including but not limited to digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA), and other programmable logic devices configured or programmed according to the teachings of the present disclosure. Computer instructions or software codes running in the general purpose or specialized computing devices, computer processors, or programmable logic devices can readily be prepared by practitioners skilled in the software or electronic art based on the teachings of the present disclosure.
In some embodiments, the present invention includes computer storage media having computer instructions or software codes stored therein which can be used to program computers or microprocessors to perform any of the processes of the present invention. The storage media can include, but is not limited to, floppy disks, optical discs, Blu-ray Disc, DVD, CD-ROMs, and magneto-optical disks, ROMs, RAMs, flash memory devices, or any type of media or devices suitable for storing instructions, codes, and/or data.
The present invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The present embodiment is therefore to be considered in all respects as illustrative and not restrictive. The scope of the invention is indicated by the appended claims rather than by the foregoing description, and all changes that come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.

## Claims

1. A computer-implemented method for providing admission control to a wireless network system, the system having one or more high-priority (HP) nodes, one or more low-priority (LP) nodes, contending for access to an access point, the method comprising:
categorizing the system into one LP access category (AC) and I HP ACs, wherein the LP AC has n0 LP nodes, and each of the LP nodes has a same packet size L0 and generates a first random backoff count uniformly distributed in [0,CW0] for each of new transmission or retransmission, where CW0 is a pre-configured contention window (CW) size required by a fixed-CW algorithm, and wherein the HP AC i, where 1≦i≦I, has ni HP nodes, where n1+L+nI=n is a total number of HP nodes, and each of the HP nodes in the HP AC i has a same packet size Li and a same packet arrival rate λi, and generates a second random backoff count uniformly distributed in [0,CWD] for each of new transmission or retransmission, where CWD is a contention window size dynamically set by a dynamic-CW algorithm;
setting all of the LP and HP nodes having the packet size with L;
determining an optimal asymptotic HP attempt rate as follows:





  
    β
    opt
    +
  
  =
  
    
      
        W
        0
      
      ⁡
      
        (
        
          
            
              C
              0
            
            ⁡
            
              (
              
                σ
                -
                
                  T
                  c
                
              
              )
            
          
          
            e
            ⁢
            
              
            
            ⁢
            
              T
              c
            
          
        
        )
      
    
    +
    1
  



where w0(·) is a Lambert W(z) function, W(z)eW(z)=z, for any complex number z,
σ is a length of a media access control MAC slot,
C0(1−β0)n0,
n0 is the total number of LP nodes,
β0 is an average attempt rate per slot for each of the LP AC nodes,
Tc is a mean time in slots for an unsuccessful transmission;
determining an optimal asymptotic HP throughput Γ1(βopt+) by substituting the βopt+ for β+ in Γ1(β+) as follows:





  
    
      Γ
      1
    
    ⁡
    
      (
      
        β
        +
      
      )
    
  
  =
  
    
      
        e
        
          -
          
            β
            +
          
        
      
      ⁢
      
        β
        +
      
      ⁢
      
        C
        0
      
      ⁢
      L
    
    
      
        T
        c
      
      +
      
        
          
            e
            
              -
              
                β
                +
              
            
          
          ⁡
          
            [
            
              σ
              -
              
                T
                c
              
            
            ]
          
        
        ⁢
        
          C
          0
        
      
    
  



where Γ1(β+) is an asymptotic total HP throughput, and
β+ is a total asymptotic HP attempt rate;
determining one or more numerical values of operating parameters including λi such that Σi=1IniλiL<Γ1(βopt+),
where Σi=1IniλiL is a total traffic load of all of the HP ACs;
determining the CWD based on a requirement of total delay of a flow.

2. The method of claim 1, wherein the total delay of the flow is related to the determined numerical values of λi.

3. The method of claim 1, wherein the fixed-CW algorithm is a variant of binary-exponential-backoff (BEB) algorithm.

4. The method of claim 1, wherein the CW0>1.

5. The method of claim 1, wherein the β0 is calculated as follows:


β0=2/(CW0+1).

6. The method of claim 1, wherein the admission control is for homogeneous traffic.

7. A computer-implemented method for providing admission control to a wireless network system, the system having one or more high-priority (HP) nodes, one or more low-priority (LP) nodes, contending for access to an access point, the method comprising:
categorizing the system into one LP access category (AC) and I HP ACs, wherein the LP AC has n0 LP nodes, and each of the LP nodes has a same packet size L0 and generates a first random backoff count uniformly distributed in [0, CW0] for each of new transmission or retransmission, where CW0 is a pre-configured contention window (CW) size required by a fixed-CW algorithm, and wherein the HP AC i, where 1≦i≦I, has ni HP nodes, where n1+L+nI=n is a total number of HP nodes, and each of the HP nodes in the HP AC i has a same packet size Li and a same packet arrival rate λi, and generates a second random backoff count uniformly distributed in [0,CWD] for each of new transmission or retransmission, where CWD is a contention window size dynamically set by a dynamic-CW algorithm;
determining an optimal asymptotic aggregate attempt rate as follows:





  
    β
    opt
    +
  
  =
  
    
      
        W
        0
      
      ⁡
      
        (
        
          
            
              
                C
                0
              
              ⁡
              
                (
                
                  σ
                  -
                  
                    T
                    c
                  
                
                )
              
            
            +
            
              
                C
                1
              
              ⁡
              
                (
                
                  
                    T
                    
                      s
                      0
                    
                  
                  -
                  
                    T
                    c
                  
                
                )
              
            
          
          
            e
            ⁢
            
              
            
            ⁢
            
              T
              c
            
          
        
        )
      
    
    +
    1
  



where w0(·) is a Lambert W(z) function, W(z)eW(z)=z, for any complex number z,
σ is a length of a media access control MAC slot,
C0(1−β0)n0,
C1n0β0(1−β0)n0−1,
n0 is the total number of LP nodes,
β0 is an average attempt rate per slot for each of the LP AC nodes,
Tc is a mean time in slots for an unsuccessful transmission,
Ts0 is a mean time in slots of a successful transmission for each of the AC 0 nodes,
determining a maximum HP throughput Γ2(βopt+,r1,L,rI) by substituting the βopt+ for β+ in Γ2(β+,r1,L,rI) as follows:





  
    
      Γ
      2
    
    ⁡
    
      (
      
        
          β
          +
        
        ,
        
          r
          1
        
        ,
        L
        ,
        
          r
          I
        
      
      )
    
  
  =
  
    
      
        β
        +
      
      ⁢
      
        e
        
          -
          
            β
            +
          
        
      
      ⁢
      
        C
        0
      
      ⁢
      
        
          D
          0
        
        ⁡
        
          (
          
            
              r
              1
            
            ,
            L
            ,
            
              r
              I
            
          
          )
        
      
    
    
      
        
          
            
              
                e
                
                  -
                  
                    β
                    +
                  
                
              
              ⁢
              
                
                  C
                  0
                
                ⁡
                
                  [
                  
                    σ
                    +
                    
                      
                        β
                        +
                      
                      ⁢
                      
                        
                          D
                          1
                        
                        ⁡
                        
                          (
                          
                            
                              r
                              1
                            
                            ,
                            L
                            ,
                            
                              r
                              I
                            
                          
                          )
                        
                      
                    
                    +
                    
                      
                        C
                        1
                      
                      ⁢
                      
                        
                          T
                          
                            s
                            0
                          
                        
                        /
                        
                          C
                          0
                        
                      
                    
                  
                  ]
                
              
            
            +
          
        
      
      
        
          
            
              [
              
                1
                -
                
                  
                    e
                    
                      -
                      
                        β
                        +
                      
                    
                  
                  ⁢
                  
                    (
                    
                      
                        C
                        0
                      
                      +
                      
                        C
                        1
                      
                      +
                      
                        
                          β
                          +
                        
                        ⁢
                        
                          C
                          0
                        
                      
                    
                    )
                  
                
              
              ]
            
            ⁢
            
              T
              c
            
          
        
      
    
  



where Γ2(β+,r1,L,rI) is a maximum system throughput,
β+ is a total asymptotic HP attempt rate,





  
    
      
        D
        0
      
      ⁡
      
        (
        
          
            r
            1
          
          ,
          L
          ,
          
            r
            I
          
        
        )
      
    
    =
    
      
        
          Σ
          
            i
            =
            1
          
          I
        
        ⁢
        
          r
          i
        
      
      
        
          Σ
          
            i
            =
            1
          
          I
        
        ⁢
        
          
            r
            i
          
          
            L
            i
          
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        D
        1
      
      ⁡
      
        (
        
          
            r
            1
          
          ,
          L
          ,
          
            r
            I
          
        
        )
      
    
    =
    
      
        
          Σ
          
            i
            =
            1
          
          I
        
        ⁢
        
          
            r
            i
          
          
            L
            i
          
        
        ⁢
        
          T
          
            s
            i
          
        
      
      
        
          Σ
          
            i
            =
            1
          
          I
        
        ⁢
        
          
            r
            i
          
          
            L
            i
          
        
      
    
  
  ,



ri is a ratio between HP AC i and AC 1 throughput,
Tsi is a mean time in slots of a successful transmission for each of the AC i nodes,
determining one or more numerical values of operating parameters including λ1 such that Σi=1IniλiLi<Γ1(βopt+,r1,L,rI)
where Σi=1IniλiLi is a total traffic load of all of the HP ACs;
determining the CWD based on a requirement of total delay of a flow.

8. The method of claim 7, wherein the total delay of the flow is related to the determined numerical values of λi.

9. The method of claim 7, wherein the fixed-CW algorithm is a variant of binary-exponential-backoff (BEB) algorithm.

10. The method of claim 7, wherein the CW0>1.

11. The method of claim 7, wherein the β0 is calculated as follows:


β0=2/(CW0+1).

12. The method of claim 7, wherein the admission control is for heterogeneous traffic.