Riding on the Primary: A New Spectrum Sharing Paradigm for WirelessPowered IoT Devices
Abstract
In this paper, a new spectrum sharing model referred to as riding on the primary (ROP) is proposed for wirelesspowered IoT devices with ambient backscatter communication capabilities. The key idea of ROP is that the secondary transmitter harvests energy from the primary signal, then modulates its information bits to the primary signal, and reflects the modulated signal to the secondary receiver without violating the primary system’s interference requirement. Compared with the conventional spectrum sharing model, the secondary system in the proposed ROP not only utilizes the spectrum of the primary system but also takes advantage of the primary signal to harvest energy and to carry its information. In this paper, we investigate the performance of such a spectrum sharing system under fading channels. To be specific, we maximize the ergodic capacity of the secondary system by jointly optimizing the transmit power of the primary signal and the reflection coefficient of the secondary ambient backscatter. Different (ideal/practical) energy consumption models, different (peak/average) transmit power constraints, different types (fixed/dynamically adjustable) reflection coefficient, different primary system’s interference requirements (rate/outage) are considered. Optimal power allocation and reflection coefficient are obtained for each scenario.
Cognitive radio, spectrum sharing, backscatter communication, IoT, wirelesspowered network, energy harvesting, ergodic capacity, power allocation, optimization.
1 Introduction
Internet of Things (IoT) is a key application scenario of the fifth generation (5G) mobile communication systems. It covers a wide range of use cases, such as smart home, smart wearables, smart farming, smart manufacturing, smart utilities, and smart city, which enable new business opportunities and new operational considerations for 5G. With the diverse use cases anticipated in IoT, the types of IoT devices are expected to diversified, and the characteristics and demands of different IoT devices are expected to vary a lot. Some of the devices, such as sensors and RadioFrequency Identification (RFID) tags, are expected to be simple, small, low power, low throughput field devices. For this kind of IoT devices, a key requirement from the industry [2] is that the power consumption should be very low and the battery life should be as long as ten years for extreme use cases. In these situations, energy harvesting, with potential to provide a perpetual power supply, becomes an attractive approach to prolong these devices’ battery lifetime. Classic sources for energy harvesting include solar and wind. Recently, ambient radio signal [3][6] is receiving much research attention as a new viable source for energy harvesting, supported by the advantage that the wireless signals can carry both energy and information.
The backscatter communications technology used in RFID systems is a realworld application of energy harvesting from RF signals. In a typical backscatter system [7, 8], the reader transmits a RF sinusoidal signal to a passive tag. The passive tag harvests RF energy from the signal to power its circuit, modulates its information bits onto the received sinusoidal signal by intentionally changing its amplitude and/or phase which is realized by changing its antenna impedance, and reflects the modulated signal back to the reader. In [9], ambient backscatter, which is able to harvest energy from and transmits information over the ambient RF signals (e.g. TV signals), was proposed. In [10], WiFi backscatter that uses the existing WiFi infrastructure to provide internet connectivity for RFpowered devices was proposed. In [11] and [12], maximumlikelihood detection was studied for an ambient backscatter system in which the tag adopts differential modulation. In [13], the modulator and the decoder design for backscattering over the ambient orthogonal frequency division multiplexing (OFDM) signals was studied. It was shown in recent work [14] and [15] that power harvested from ambient RF signals is sufficient to support the daily communication of a batteryless sensor through dedicated decoder and signal design. In [16], it was further shown a network architecture with combined wireless power transfer and backscattering communication technology is able to a large and dense IoT network. However, these aforementioned works mainly focus on the hardware and decoder design but lacks fundamental system analysis from theoretical aspects. Besides, the mutual influence (such as the interference) between the backscatter system and the primary system are not considered.
In this paper, we introduce the ambient backscatter communication technology to the cognitive radio (CR) system, and propose a new spectrum sharing model based on that. The proposed spectrum sharing model applies to CR systems with conventional primary communication systems and ambientbackscatterbased secondary systems. The key idea is that the secondary transmitter harvests energy from the primary signal, then modulates its information bits to the primary signal, and reflects the modulated signal to the secondary receiver without violating the primary system’s interference requirement [17]. The main contribution of this paper is summarized as follows:

We propose a new spectrum sharing model as Riding on the Primary (ROP). The differences between the proposed ROP and the existing technologies are as follows: (i) Compared with conventional backscatter communication systems, the reader (secondary receiver) in our system does not need to generate and transmit a RF sinusoidal signal, which can reduce its power consumption and prolongs its battery life. (ii) Compared with ambient backscatter communication systems, the interference from the wirelesspowered tag (secondary transmitter) to the primary system is taken into consideration when designing the system. (iii) Compared with conventional spectrum sharing systems, the secondary system not only utilizes the spectrum of the primary system but also takes advantage of the primary system’s signal transmission to carry its information.

We investigate the performance of the proposed ROP system under fading channels. To be specific, we maximize the ergodic capacity of the secondary system by jointly optimizing the transmit power of the primary signal and the reflection coefficient of the secondary ambient backscatter. Different (ideal/practical) energy consumption models, different (peak/average) transmit power constraints, different types (fixed/dynamically adjustable) reflection coefficient, different primary system’s interference requirements (rate/outage) are considered. Optimal power allocation and optimal reflection coefficient are obtained for each case. It is worth pointing out that certain degree of cooperation is needed between the primary and secondary system in order to do the joint optimization. For example, the primary system needs the knowledge of the secondary system (e.g., channel information) to optimize the transmit power, whereas the secondary system needs the knowledge of the primary system to optimize the reflection coefficient.

We show by numerical examples that dynamically adjustable reflection coefficient will result in a better system performance than a fixed reflection coefficient. The system performance under the average transmit power constraint is in general better than that under the peak transmit power constraint. Applying the primary transmission outage probability constraint to protect the primary system usually leads to a higher secondary ergodic capacity than adopting the primary transmission rate constraint. All these findings can serve as the guidance for designing highperformance practical ROP system.
The rest of this paper is organized as follows. Section II introduces the system model of the proposed ROP spectrum sharing system. Section III presents the basic problem formulation. Section IV investigates the capacity maximization problem under the practical energy consumption model. Section V studies the capacity maximization problem for fixed reflection coefficient and average transmit power constraint. Section VI considers the capacity maximization problem under the primary transmission outage constraint. Then, numerical results are given in Section VII to verify the proposed studies. Finally, Section VIII concludes the paper.
2 System Model
2.1 Channel Model
In this paper, we consider a spectrum sharing communication system consists of a primary communication pair and a secondary transmission pair. The primary communication pair is a conventional communication system consisting of a RF source (e.g. Base Stations, TV towers, WiFi APs) and a receiver (e.g., cell phones, TV receivers). The secondary communication pair is an ambient backscatter communication system which consists of a wirelesspowered passive tag and a batterypowered reader. For ease of explanation, we denote the RF source and the receiver of the primary system as the primary transmitter (PT) and the primary receiver (PR), respectively. The wirelesspowered tag and the reader of the backscatter system are denoted as secondary transmitter (ST) and secondary receiver (SR), respectively. In this paper, we consider the block fading channel model [17], where the channel coefficients remain the same for each block but may change from one block to another. As shown in Fig. 1, the channel power gains for the fading block , from the PT to the PR, from the ST to the SR are denoted by and , respectively. The channel power gains for cross channels for fading block , i.e., from the PT to the ST, from the PT to the SR, and from the ST to the PR, are denoted by , and , respectively.
2.2 Transmission Model
Transmitted signal at the PT. Let denote the transmitted signal of the PT at th symbol of the th block where , and denote the transmit power for fading block . Then, the transmitted signal of the PT for the th symbol of block is given by
(1) 
Transmitted signal at the ST. In fading block , the signal received at the ST from the PT is . Note that the noise at the ST (Tag) is neglected as [9, 10, 11, 12, 14, 18] since the ontag integrated circuit only includes passive components. The power of the received signal at the ST from the PT is . Part energy of the received signal is absorbed by the ST to power its circuit operation. The remaining part of the received signal is modified and backscattered to the reader. For convenience, we refer to this splitting factor as the reflection coefficient, and denote it by where . Then, the energy of the transmitted signal of the ST can be denoted as . Let where denote the ST’s own signal, then the transmitted signal of the ST for the th symbol of block is given by
(2) 
Note that we assume there is no signal processing delay of the backscatter circuit, i.e., there is no time delay between the transmitted signal and the received signal of the ST. This assumption is widely used in backscatter communication research literatures [9][13].
Received signal at the PR. Let denote the received signal at the PR for the th symbol of block , then we have
(3) 
where denotes the Gaussian receiving noise at the PR with zero mean and variance . Then, the instantaneous received signaltointerferenceplusnoise ratio (SINR) at the PR for block denoted by is given by
(4) 
Received signal at the SR. Let denote the received signal at the SR for the th symbol of block , then we have
(5) 
where denotes the Gaussian receiving noise at the SR with zero mean and variance . In this paper, we assume that SR decodes the received signal by performing successive interference cancellation (SIC). SIC is a wellknown physical layer technique [19, 20]. Briefly, SIC is the ability of a receiver to receive two or more signals concurrently. The SIC receiver decodes the stronger signal first, subtracts it from the combined signal, and extracts the weaker one from the residue. For the system setup considered in this paper, the secondary system is an ambient backscatter system. The strength of the signal received from the ST is in general much lower than that received from the PT (e.g., TV/WiFi signals). Thus, the SIC procedure at the SR is decoding the primary signal first and then subtracting it from the received signal before decoding its own signal. Thus, the instantaneous received SNR at the SR for the block denoted by is given by
(6) 
Note that we assume that the SR performs SIC decoding while the PR does not. This is due the fact that SIC needs to decode the stronger signal first. As mentioned above, the strength of the signal received from the ST is in general much lower than that received from the PT. Thus, SIC is not applicable at the PR.
3 Ergodic Capacity Maximization
Under the system model given in Section 2, the ergodic capacity of the secondary system can be written as
(7) 
where denotes the statistic expectation, and it is taken over the joint fading states of the fading block . In this paper, our objective is to maximize the ergodic capacity of the secondary system while guaranteing the performance of the primary system. In the following, we introduce the constraints that need to be considered when optimizing this network.
PT’s transmit power constraint. Let denote the maximum transmit power of the PT, then the peak transmit power constraint can be written as
(8) 
The reflection coefficient constraint. Since the tag is a passive device, thus the energy harvested and reflected from the tag must be equal to the energy received from the primary signal. Thus, the reflection coefficient must satisfy the following constraint
(9) 
PR’s rate constraint. To guarantee the quality of service (QoS) of the primary system, we assume that there is a minimum rate requirement, which can be written as
(10) 
where is the minimum rate of the primary system.
Tag’s (ST’s) circuit operation power requirement. As aforementioned, the ST harvests energy from the primary signal to power its circuit operation. Let be the minimum power that the ST needs to support its circuit operation, then the following constraint must be satisfied in order for the ST to work, i.e.,
(11) 
where is the energy harvesting efficiency coefficient.
In this paper, the objective is to optimize the performance of such a spectrum sharing system by jointly optimizing the transmit power of the PT and the reflection coefficient of the ST. The problems can be formulated as
(12)  
s.t.  (13) 
For notation convenience, the fading block number is dropped from now on.
P1 is not a convex optimization problem. A problem is a convex optimization problem if its objective function is either convex or concave, and its feasible set is a convex set. For P1, its objective function is neither convex nor concave. This can be verified by looking at its Hessian matrix (A function is convex if its Hessian matrix is positive semidefinite). According to the composition rule given in [21], for P1, we can determine the convexity of the objective function by investigating the convexity of . The Hessian matrix of is , which is not positive semidefinite. Thus, P1 is not a convex optimization problem, and the conventional convex optimization techniques can not be applied to solve P1. To solve P1, we first introduce the following Lemma.
Lemma 1. The largest that makes P1 feasible, denoted by , is given by
(14) 
Proof.
The constraint (10) can be rewritten as
(15) 
It is observed that the right hand side of (15) is an increasing function of . Besides, due to the fact that , the largest denoted by satisfies (15) is
(16) 
Proof.
First, it is observed from P1 that all the constraints are instantaneous constraints. Thus, maximizing the ergodic capacity is equivalent to maximizing the instantaneous transmission rate, i.e., . For any given in the feasible region, the instantaneous rate is a monotonically increasing function with respect to , and it attains the maximum value when . It is observed that for any given feasible , the instantaneous rate is a monotonically increasing function with respect to . Thus, should be chosen as the largest that makes P1 feasible, i.e., given by (14), and it is shown in the proof of Lemma 1 that is obtained when . Thus, it is clear that P1 is maximized when and . ∎
4 Practical Energy Consumption Model
In this section, we consider a more practical energy consumption model of the ST’s backscatter circuit, which is
(21) 
where denotes the static energy consumption when the circuit is on, and denotes the dynamic energy consumption which is a function of its transmission rate. In practice, the dynamic energy consumption is in general proportional to the transmission rate. This is due to the following fact. The backscatter transmitter maps its bit sequence to RF waveforms by adjusting the load impedance of the antenna. The backscatter can control the rate at which it will generate the modulation symbols by controlling the switching frequencies on the SPDT (single pole double throw) switches [14]. In general, to achieve a higher data rate, the backscatter needs a higher frequency operation on the switches, which will cost more energy. Thus, we model the dynamic energy consumption by , where is a constant conversion parameter that relates the transmission rate with the energy consumption. In this section, under this energy consumption model, we reinvestigate the optimization problem for this spectrum sharing system.
(22)  
s.t.  (23) 
The constraint (21) is a hyperfunction with respect to and , which makes the problem difficult to solve. Thus, to solve P2, we first present the following two propositions.
Proposition 1. Let denote the largest that satisfies the constraint (21) for a given , then can be obtained by solving the following equation:
(24) 
Proof.
It is observed that the left hand side of (21) (i.e., ) is a monotonically decreasing function with respect to , while the right hand side of (21) () is a monotonically increasing function with respect to . Thus, it is easy to observe that the largest is the intersection point of two curves, which is the solution of (24). Proposition 1 is thus proved. ∎
Proposition 2. Let and be the solution of (24) when and , respectively. Then, we have
(25) 
Proof.
Since and be the solution of (24) when and , respectively. We have
(26) 
(27) 
Then, in the following, we prove Proposition 2 by contradiction. Assume when . Then, under this presumption, it is clear that the left hand side of (28) is strictly positive, while the right hand side of (28) is zero or negative. This contradicts with the fact that the left hand side of (28) should be equal to the right hand side of (28). Thus, our presumption does not hold. Thus, it follows that Proposition 2 is thus proved. ∎
However, unlike P1, we cannot further prove . Thus, the approach used to solve P1 can not be applied here. Thus, to solve P2, we first consider the following problem, which is
(29)  
s.t.  (30) 
Let and be the solution of (24) when and , respectively. Then, from Proposition 2, it follows Since the objective function is an increasing function with respect to , it is clear that the objective function attains its maximum value at . Thus, the optimal solution of P2a can be obtained as
(31)  
(32) 
Now, we return to P2. It is clear that the constraint (10) can be rewritten as It is observed that its right hand side is an increasing function of . Besides, due to the fact that , the largest satisfying denoted by is
(33) 
Based on these results, we are now able to solve P2, and the solution is summarized in the following theorem.
Theorem 2. The optimal solution of P2 is given by
(34)  
(35) 
Proof.
To prove Theorem 2, we consider the following two cases.
Case 1: . In this case, the optimal solution of P2 is the same as that of P3a, which is and .
Case 2: . In this case, the optimal must satisfy the condition . Now, we look at the following equation,
(36) 
The above equation comes from the fact that is the solution of (24) when . Since , thus if we replace with , the left hand side of (36) will increase, while the right hand side of (36) will decrease. Thus, it follows
(37) 
which indicates and is a feasible solution of P2. Since is the largest feasible and is the largest feasible , thus the optimal solution is and .
Summarizing the above results, Theorem 2 follows. ∎
5 Fixed Reflection Coefficient and Average Transmit Power Constraint
Fixed Reflection Coefficient. In practice, for the purpose of circuit design simplicity, the reflection coefficient of the tag is designed to be fixed, i.e., the reflection cannot be dynamically changed in each fading block. To capture this fact, we introduce the following constraint
(38) 
Average transmit power constraint. In practice, there is always a longterm power budget of the PT, and an average power constraint usually applies, which can be written as
(39) 
where the statistic expectation is taken over the joint fading states of the fading block.
In the following subsections, we reinvestigate P1 and P2 under the above two practical constraints.
5.1 Ideal Energy Consumption Model
(40)  
s.t.  (41) 
The objective function of P3 is the same as P1, and thus it is easy to verify that P3 is not a convex optimization problem. Therefore, it cannot be solved directly using convex optimization techniques. Thus, for solving P3, we first consider P3 under a given . For given , P3 can be rewritten as
(42)  
s.t.  (43)  
(44)  
(45) 
Note the fading block number is dropped from now on for notation convenience.
P3a is a convex optimization problem since its objective function is concave and all the constraints are affine. To solve this problem, we first at the feasibility of the problem. Note that the constraint (44) is infeasible if , i.e., no matter how is chosen, (44) cannot be satisfied for such fading block. Thus, to save power, the optimal power allocation for such fading block is . When feasible, it can be shown that constraints (44) and (45) can be rewritten as and , respectively. Thus, (44) and (45) can be replaced by the following constraint
(46) 
where
(47) 
Theorem 3. The optimal solution of P3a is given by
(48) 
where can be obtained by solving , and is given by (47).
Proof.
By introducing the dual variable associated with the average transmit power constraint, the partial Lagrangian of P3a problem is expressed as
(49) 
where is nonnegative Lagrange dual variable associated with constraints (43).
Let denote the set of . The dual function is then expressed as The Lagrange dual problem is then defined as . The duality gap is zero for the convex optimization problem addressed here, and thus solving its dual problem is equivalent to solving the original problem. The duality gap is zero if and only if strong duality holds. If the primal problem is convex and satisfies the Slater¡¯s condition [21], then strong duality holds. The Slater’s condition reduces to feasibility when the constraints are all linear equalities and inequalities [21]. P3a is a convex problem and all its constraints are linear inequalities. It also can be verified that P3a is feasible under our constraints. Thus, the duality gap is zero for P3a. Therefore, according to the KarushKuhnTucker (KKT) conditions, the optimal solutions needs to satisfy the following equations:
(50)  
(51) 
For a fixed , by dual decomposition, the dual function can be decomposed into a series of similar subdualfunctions each for one fading state. For a particular fading state, the problem can be shown equivalent to
(52)  
s.t.  (53) 
It is easy to observe that the optimal solution of this subproblem is if . Thus, in the following, we consider the optimal solution of this subproblem under the condition that .
The Lagrangian of this subproblem is
(54) 
where and are nonnegative Lagrange dual variables associated with the constraints and , respectively. Since the problem is convex, KKT conditions are sufficient to obtain its optimal solution. Thus, in the following, we investigate its KKT conditions:
(55)  
(56)  
(57)  
(58)  
(59) 
Now, we derive the optimal solution by solving these KKT conditions. To solve these KKT conditions, we consider the following two cases:
Thus, combining the results obtained in Case 1 to 4, the optimal solution for the subproblem can be summarized as (48). Now, we have to find the optimal . As aforementioned, has to satisfy (87) and (89). It is observed that if , the optimal solution for the subproblem is . This definitely violates the constraint . Thus, it follows that . Then, according to the constraint , it follows that the optimal must satisfy . Theorem 2 is thus proved. ∎
With the optimal solution of P3a given in Theorem 3, the optimal solution of P3 can be obtained by performing a onedimension search for over the space .
5.2 Practical Energy Consumption Model
(63)  
s.t.  (64) 
The objective function of P4 is the same as P1, and thus it is easy to verify that P4 is not a convex optimization problem. Therefore, it cannot be solved directly using convex optimization techniques. Thus, before solving P4, we first present the following proposition.
Proposition 3. For a given , the constraint (21) can be rewritten as
(65) 
where is the positive solution of .
Proof.
For the convenience of exposition, we introduce the following two functions and . It is easy to observe that both and are monotonically increasing functions with respect to , and the increasing rate of is larger than that of . It is also observed that . Thus, and must have a positive crossing point . For any larger than , we have . Proposition 3 is thus proved. ∎
It can be shown that constraints (10) can be rewritten as . Thus, (10) and (21) can be replaced by the following constraint
(66) 
where
(67) 
Using the same approach as P3, for a given , the optimal power allocation of P4 can be summarized as follows.
Theorem 4. For a given , the optimal power allocation of P4 is given by
(68) 
where can be obtained by solving , and is given by (67).
With the optimal solution of P4 under given obtained in Theorem 4, the optimal solution of P4 can be obtained by performing a onedimension search for over the the space .
6 Primary Transmission Outage Constraint
In previous sections, we use the PR’s rate constraint to guarantee the quality of service (QoS) of the primary system, i.e., there is a minimum rate requirement that the primary transmission has to fulfil. However, this constraint is too strict. In practice, certain ratio of transmission outage is usually acceptable. Thus, in this section, we introduce the primary transmission outage constraint [22], which is mathematically defined as
(69) 
where denotes the probability of the event, and denotes the maximum outage probability that is acceptable by the primary system.
6.1 Peak Transmit Power Constraint
We first consider the peak transmit power constraint case. Using the primary transmission outage constraint to replace the PR’s rate constraint, and considering the fixed reflection coefficient, the problem can be formulated as
(70)  
s.t.  (71) 
The objective function of P5 is the same as P1, and thus it is easy to verify that P5x is not a convex optimization problem. Therefore, it cannot be solved directly using convex optimization techniques. We first consider P5 under a given , i.e., where is a constant. For expression convenience, we denote the problem of P5 under a given as P5a.
Theorem 5. The optimal power allocation of P5a, i.e., P5 under a given , is given by
(72) 
Proof.
For a given , the constraint (11) can be rewritten as . For a given fading block , it is easy to observe that P5a is infeasible when . This means that there does not exist a feasible that can satisfy the ST’s circuit operation power requirement, which indicates the secondary system does not work in this fading block. Thus, the contribution of this fading block to the secondary transmission’s ergodic capacity is zero. The interference from the secondary system to the primary system is also zero. Thus, channel inversion power allocation or the maximum transmit power should be adopted for the primary system to minimize the possibility of outage for this fading block. Therefore, for simplicity, the optimal power allocation for such a fading block can be obtained as
(73) 
In the following, we consider the case that . To solve the problem, we introduce the following indicator function
(74) 
For notation convenience, the fading block number is dropped from now on. Then, it can be shown that the constraint given in (69) is equivalent to
(75) 
By introducing the dual variable associated with (75), the partial Lagrangian of P5 under given can be expressed as
(76) 
where is nonnegative Lagrange dual variable associated with the constraint .
Let denote the set of . The dual function is then expressed as According to the KarushKuhnTucker (KKT) conditions, the optimal solution needs to satisfy the following equations:
(77)  
(78) 
For a fixed , by dual decomposition, the dual function can be decomposed into a series of similar subdualfunctions each for one fading state. For a particular fading state, the problem can be shown equivalent to
(79)  
s.t.  (80)  
(81) 
Since is an increasing function, the optimal solution of P5b is if .
In the following, we consider the optimal solution of this subproblem under the condition that . Fisrt, it can be observed that can be rewritten as
Case 1: .
(82) 
Case 2: .
(83) 
Let us first consider Case 1. In this case, the objective function of P5b becomes , which is an increasing function with respect to . Thus, the optimal solution of P5b is .
Now we consider Case 2. It is observed from (83) that is a step function with respect to p for this case. The critical point is . Thus, the objective function of P5b has a form as illustrated in Fig. 2. It can be observed from Fig. 2 that the function is in general an increasing function of . Thus, the optimal solution of P5b is .
Combining all these results, the optimal solution for P5a can be obtained as (72). Theorem 5 is thus proved. ∎
With the optimal solution of P5a (i.e, P5 under given ) obtained in Theorem 5, the optimal solution of P5 can be obtained by performing a onedimension search for over the space .
6.2 Average Transmit Power Constraint
Now, we investigate the problem under the average transmit power constraint, which is
(84)  
s.t.  (85) 
We first consider P6 under a given , i.e., where is a constant. For expression convenience, we denote the problem of P6 under a given as P6a.
Same as solving P5, it can be shown that (69) can be rewritten as (75). By introducing the dual variables associated with (75) and (39), the partial Lagrangian of P6a under given can be expressed as
(86) 
where and are nonnegative Lagrange dual variables associated with the constraints and , respectively.
Let denote the set of . The dual function is then expressed as According to the KarushKuhnTucker (KKT) conditions, the optimal solution needs to satisfy the following equations:
(87)  
(88)  
(89)  
(90) 
For a fixed , by dual decomposition, the dual function can be decomposed into a series of similar subdualfunctions each for one fading state. For a particular fading state, the problem can be shown equivalent to
(91)  
s.t.  (92) 
P6a can be solved by iteratively solving P6b for all fading states with fixed and , and updating these dual variables via subgradient based method, e.g., the ellipsoid method [21], for which the details are omitted here for brevity.
In the following, we derive the solution to P6b. To solve the problem, we first introduce the following function
(93) 
It is easy to show that is a concave function with respect to , and attains its maximum value when is equal to
(94) 
Note that the objective function in P6b now becomes . Now we look at . It has been shown in the previous subsection that can be rewritten as (82) and (83).
Thus, when , we have . Then, the objective function of P6b becomes , which attains its maximum value at the same as . Thus, the optimal solution of P6b for this case is .
When , as shown in (83), is a step function with respect to , and its critical point is . In the following, we consider P6b for this case. For the convenience of discussion, we denote the critical point of as , i.e.,
(95) 
For the convenience of discussion, we rewrite the constraint (92) as , where
(96) 
Let denote the optimal solution of P6b. The following discussions are then made on :