Weighted Sum Rate Maximization in OFDM Based Cooperative Cognitive Radios : A Joint Optimization Approach

Improving utilization of the radio spectrum is the main goal of Cognitive Radio Networks (CRN). Recent studies made use of cooperative relay technology in cognitive networks, to increase transmission diversity gain. In this paper we consider an OFDM based cooperative cognitive network with a pair of Source-Destination nodes as the primary user (PU), and a pair of Source-Destination nodes—which is assisted with a relay—as the secondary (cognitive) user (SU). Both primary and secondary users share a same spectrum. In a two hop transmission, the source transmits in the first hop, and the half-duplex relay decodes the message, re-encodes and forwards it to the destination in the second hop on a different subcarrier. The cognitive network obeys an underlay paradigm where the SU is allowed to transmit simultaneously with PU, while its power is limited such that the interference caused for PU does not exceed a defined temperature. Under this constraint, a joint subcarrier pairing and power allocation is proposed for SU to maximize its weighted sum rate. The problem is transformed to a convex optimization problem and solved in the dual domain. Then an algorithm to achieve feasible solutions is used based on the optimization results. Through extensive simulations, we compare the spectrum utilization of the proposed approach with the existing ones, and show that interestingly the proposed method improves the weighted sum rate of SU.


Introduction
In the era of communications, daily increasing demands for wireless services and its high penetration in routine human activities is overloading the precious radio spectrum.Beside this bottleneck, and according to spectrum efficiency measurements recently revealed by FCC [1], the inefficient usage of available spectrum in fixed assignments scenarios causes the effort for more utilizing the scarce spectrum to be the focus of many researches.Emerging Cognitive Radio Networks (CRN)-coined by Joseph Mitola in his Doctoral Dissertation [2]-is becoming an efficient solution to this problem in recent years.Diverse aspects of CRNs have been studied in many articles since then.A detailed exposition of Cognitive Radios (CR) from the signal processing point of view has been presented in [3].Development of the IEEE 802.22 as the first Cognitive Radio Wireless Regional Area Network Standard is a worldwide effort to embark use of this novel technique [4].As it is categorized in [5], there are three well known paradigms for unlicensed use of spectrum in CRNs: Interweave, Overlay, and Under-lay.
Among different system models defined in recent studies, some made use of cooperative relay technology in cognitive networks to increase transmission diversity gain.In [6] a number of cooperative diversity cognitive models are discussed by combining a cognitive radio with the cooperative technique.A brief overview of applications of cooperative cognitive networks is presented in [7].Also cooperative transmission between secondary nodes to improve spatial diversity and improve spectrum diversity with a focus on the latter one is discussed, and a new MAC protocol is proposed [7].
A cognitive cooperative relaying scheme in which cognitive users assist PUs by relaying their information is proposed in [8]; where a Decode-and-Forward (DF) strategy is used, and the outage probability of PU is investigated.In [9] a cooperative cognitive network with two pairs of Source-Destination as primary and cognitive users is considered, where a cognitive relay assists both primary and cognitive users.There, a half-duplex Amplify-and-Forward (AF) relaying scheme is proposed, and the power allocation optimization problem for cognitive user and cognitive relay is considered [9].A joint subcarrier and power allocation algorithm for cooperative multiuser orthogonal frequency division multiplexing (OFDM) CR is presented in [10]; where a DF cooperative strategy is used, and a three-node relay network is simplified to a tow-node network with the equivalent channel gain.In [11] joint relay selection and power allocation is investigated in a CR system to maximize system throughput by developing an optimal approach.Also a suboptimal approach is proposed which reduces complexity while maintaining performance.A relay selection in cooperative cognitive network is proposed in [12], in which a secondary node relays PU information to reduce the power consumption of PU, and simultaneously decrease interference from primary traffic which exposes spectrum opportunities for secondary system.In [13] a two-user AF relaying cognitive network is considered, and jointly allocation of sensing time and power is developed to maximize throughput of secondary network.
It can be widely seen that novel wireless communications tend to use multi-carrier OFDM method because of its flexibility in resource allocation and the benefits on minimizing the impact of multi-path fading.To be in harmony with this trend, our interest in this study lies on an OFDM-based cooperative cognitive radio network.Among three different paradigms of cognitive radios, we choose underlay transmission where both the primary and cognitive users transmit over a same spectrum.Two pairs of Source-Destination (SD) nodes are considered where one pair is the primary user and the other is the secondary user.Also we assume that the secondary user's transmission is assisted with a cognitive DF relay node.In underlay scenario, it is assumed that the secondary network is aware about channel gains from its transmitters to the PU.SU's and the relay's power are limited in this case, because simultaneous transmission with licensed users is available while cognitive transmitters (secondary source and the relay) have not violated an interference temperature which is constrained by primary network.
Our goal is to maximize the throughput of secondary network with interference power constraint defined by primary network, and total power constraint of secondary transmitter and the relay.A joint optimization of power allocation and subcarrier pairing for secondary transmitter and the DF relay node is developed.We formulate the joint power allocation and subcarrier pairing optimization problem, and then solve it using continuous relaxation, to reform the problem to dual problem.An algorithm to achieve feasible subcarrier pairing and power allocation is proposed based on the solution results.Extensive simulation results show that the proposed algorithm almost achieve the optimal weighted sum rate, and outperform existing methods in secondary throughput improvement.
The rest of this paper is organized as follows.In Section 2, we define our cooperative cognitive radio network model and formulate the Weighted Sum Rate (WSR) maximization problem by introducing our objective and constraints.Section 3 solves the optimization problem and proposes an algorithm to achieve feasible solution of the problem.In Section 4, we benchmark the proposed method through extensive simulations.Section 5 concludes the paper.

System Model and Problem Formulation
In this section, we will first define a cooperative cognitive radio system model, and then formulate the WSR maximization problem by introducing the objective function, problem constraints, and the optimization variables based on the system model.

Cooperative Cognitive Radio System Model
We consider an OFDM-based multicarrier cooperative cognitive radio system, in which both primary and secondary users share the same bandwidth for their transmissions.The secondary system model contains a base station (BS S ), a single user (SU) and a DF relay node (R), which transmits simultaneously with a primary system also containing a base station (BS P ) and one user (PU).The secondary transmission consists of two hops.In first hop (time slot) BS S transmits over subcarrier k, while both the relay and SU receive the information.In second hop the half-duplex relay retransmits the information which has been received in the first time slot.The relay uses a Decode-and-Forward (DF) method, and transmits over a different subcarrier m, in second hop.Therefore the secondary system transmission lasts two time slots using a Subcarrier Pair denoted as SP (k, m).So if we divide the total frequency band to M subcarriers we will have M subcarrier pairs for secondary system transmission.The secondary system's downlink destination (SU) receives the information from BS S and the relay node using a Maximal Ratio Combiner (MRC) receiver to exploit spatial diversity.
As shown in Figure 1 downlink transmissions are considered, and the corresponding downlink channel gains are defined.Let k SD g denote the secondary BS (Source) to SU (Destination) downlink channel gain on subcarrier k, k SR g and m RD g denote channel gains of BS S to Relay on subcarrier k, and relay to SU on subcarrier m respectively.Also let f h denote the primary system downlink channel gain on subcarrier f.The channel gains from BS S and the relay to PU are defined as and k S q m R q on subcarriers k and m respectively.As shown in Figure 1, and f h are the desired channel gains (solid lines) while and k S q m R q causes unwanted interference to PU (dashed line).Also let f P S I  denote interfering channel gain from primary BS (BS P ) to SU.It's assumed that all channel gains are independent from each other, and remain constant in a two-slot period.Moreover, all the desired channel gains of the secondary system and the interfering channel gains from secondary to primary system which are q respectively, are assumed to be perfectly known in BS S .The BS S performs subcarrier pairing and power allocation, and informs the relay node and the SU of the corresponding parameters through a control signaling prior to data transmission.In next section we will formulate our optimization problem based on the system model shown in Figure 1.

Joint Subcarrier Pairing and Power Allocation Problem Formulation
The aim of this paper is to maximize the downlink throughput of secondary system.To formulate our optimization problem we have to define our objective function, constraints, and optimization variables.Similar to [14], in our cooperative cognitive model the objective function is the achievable weighted rate of secondary system, which is denoted as where is the channel gain from secondary BS to SU on subcarrier k and is the variance of Additive White Gaussian Noise (AWGN) on the corresponding subcarrier.
Since a relay node is used in secondary transmissions, the subcarrier pair SP (k, m), may works in two different modes depending on relay's activity: direct-link mode and relay mode.In the relay mode the DF relay will forward the message in the second time slot on subcarrier m, while in the direct-link mode the message will be sent only via subcarrier k of SD link in the first time slot, and the relay is not used since it will not improve the throughput.
Now that the parameters are defined, we can formulate the achievable weighted rate of secondary downlink for SP (k, m) as [14,15] where is a weighting factor to represent Quality-of-Service (QoS) requirements on subcarrier k.Also the rate is scaled by 0 k w  1 2 since the transmission lasts two time slots.
Introducing total power constraint in secondary system as , , , for SP (k, m), it is shown [16] that using relay is advantageous and improves the achievable rate if , and .
Similar to [14], by letting following expressions for and and defining the equivalent secondary channel gain the expression in ( 2) can be re-written in a unified format as [17]   , log 1 2 Also we define which indicates whether SP (k, m) is selected or not.Now that we have defined our objective function and the required variables, the WSR optimization problem can be formulated over variables p and t as follows P1 max log 1 2 , 0, , , m a x , , where p and t are M M  matrices with entries is the maximum tolerable interference vector for primary user over each subcarrier.Constraints ( 8)-( 12) are corresponded to secondary system limitations, while constraints ( 13) and ( 14) are the primary user interference temperature thresholds. 1 , ,

Weighted Sum Rate Maximization
In this section, we will first solve the WSR maximization problem in the dual domain, and then propose an algorithm to achieve feasible solution of the problem based on the optimization results.

Solution to Optimization Problem [P1]
The defined problem [P1] is a Mixed Integer Programming (MIP) problem because of the constraint (11).Due to the fact that the constraint is a set of discrete points, MIP problem is hard to solve.The approach to this kind of problem is to approximate the problem by its continuous relaxation as in [14].In this technique, instead of constraint (11) which forces to be an integer (either 0 or 1), it will be relaxed to which is a convex constraint set and is more easier to deal with.
Also in [P1], constraints ( 13) and ( 14) are power limitations on the secondary source and relay respectively, while the optimization variable is , k m as the total power constraint.To satisfy constraints ( 13) and ( 14) and make it in correspondence with , k m , let denote the maximum total allowed power over subcarrier pair SP (k, m).Using expressions in (4) and combining constraints ( 13) and ( 14 Using continuous relaxation and maximum total power constraint in [P1], the relaxed problem becomes s.t. ( 8)-( 10), and , 0, , The objective function in (17) is the same as [P1] for The new relaxed problem is a convex optimization problem since its objective function is concave, and the constraints are convex sets [18], and therefore has a unique global solution.But the optimal value of , k m may not be an integer.Using the same approach in [14], we solve the problem in the dual domain and obtain the optimum value of , k m which satisfies (11).This solution will be an upper bound for problem (17).

t
The Lagrangian of problem (12) by using constraints (8) and ( 9) will be obtained as , , , where 0 , , , , applying constraint (19), the optimal value of will be where i is the iteration index, and  is the step size in the ith iteration.Note that the step size of ( 26) and ( 27 The original WSR problem [P1] has an integer valued constraint which causes non-convexity of the constraints set.Hence, the curve obtained for WSR optimum value by varying the total power constraint may have discrete changes in its slope.These sudden jumps in the slope of the curve are due to the changes of the optimal SP as the total power varies.This will cause the solution to the relaxed dual problem to be an upper bound for the optimum value of [P1].The gap is known as the duality gap.Recently, it is declared that the duality gap will be zero if the optimum answer of the optimization problem is a concave function of the constraints [20], and therefore the solution of the dual problem is the optimum solution to the original problem.Also it is shown in [21] and used in [14]-same approach as our problem in an OFDM downlink resource allocation problem-that the optimum value of the problem tends to be concave over constraints if the number of subcarriers (M) is large enough.Hence for practical OFDM networks (with ) [14], we obtain the optimum solution to the original problem [P1] by solving the dual problem [P2].

M 
where the following notation is used Equation ( 23) is known as "constrained water-filling" [19] since it has maximum power constraint .Also called as "multi-level water-filling" because of weighting factor w k .
Lagrange dual variables  and m  to minimize (21) with the values of , k m and , k m t obtained in the first part.In section 3.1 a closed form solution for the first part is derived in Equations ( 23) and (25).While for the second part a sub-gradient method is discussed to find the optimum values of p  and m  using iterative Equations ( 26) and ( 27).In Algorithm 1, we will first initialize the required variables.In each iteration, we will first calculate the optimum value of , k m and , k m using Equations ( 23) and (25), and then find new values for p t  and m  with the corresponding equations.The loop will continue until the dual Lagrange parameters converge to their optimum values which minimize the constraints cost.Since constraint ( 9) is guaranteed by equation (25), while constraint (10) may be violated, an amendment is used to avoid selecting the same relay subcarrier for more than one source subcarriers, which results in more than one "1"s in a column of t and violate constraint (10) [14].This amendment is applied at the last 10% of the algorithm iterations.At the last step a constrained water-filling is applied to find the optimum power allocation using to maximize WSR for secondary system while the interference threshold of the primary network is not violated.The optimum subcarrier pairing scheme and power allocation which maximize WSR of secondary system while preserving QoS of primary network in an OFDM based cooperative cognitive network is obtained by using this algorithm.2. Do while

Numerical Results and Simulations
In this section the proposed method for our Cooperative Cognitive Radio Network, as shown in Figure 1, is benchmarked through extensive simulations.As we mentioned all the channel gains q , f h , and f P S I  are independent of each other, and remain constant in a two-slot period.Also it is assumed that the channel gains are independent and identically distributed (i.i.d.) over subcarriers and follow a Rayleigh distribution with parameter b differs from one link to another corresponding to the mean square of subcarrier channel gains in different scenarios.The variance of Additive White Gaussian Noise (AWGN) is assumed to be one.
In Algorithm 1, the initial value of dual Lagrange multipliers are set as 1 1     , and .The convergence criterion is , and the step size for the subgradient method is set to be 0.05 , where i is the iteration index.In Un-Weighted Sum Rate (UWSR) scenario, the weighting factor k is set to be all "1" for different subcarriers, and in Weighted Sum Rate (WSR) scenario, it obeys the equation The results are obtained based on 1000 randomly generated channel realizations.
Two different system configurations corresponding to the position of the relay node are investigated.In the first scenario it is assumed that the relay node is placed between the secondary source and destination.The mean square of channel gains in this scenario, are set as Figure 2 shows the secondary system's WSR through different number of subcarriers.The total power con-straint for secondary system is assumed to be budget 5 P  [Watts], and the maximum tolerable interference for primary user is equal to max 1 M  I over all subcarriers.In this figure the WSR of the Fixed Subcarrier Pairing (FSP) method-where for every subcarrier pair -is compared to the Joint Subcarrier Pairing and Power Allocation (JSPPA) approach.In JSPPA1 we only apply the maximum power constraint for Constrained Water-Filling in the last part of Algorithm 1, while in JSPPA2 the is considered in both subcarrier pairing and power allocation to find optimum values and respectively.As it is shown in Figure 2, a 17% improvement in WSR is acquired with JSPPA2 when we have subcarriers.Figure 3 shows the UWSR for secondary system, exploiting FSP, JSPPA1, and JSPPA2 methods with same assumptions.In the second scenario, the relay node is placed close to the secondary source.The mean square of channel gains in this scenario, are set as All other assumptions are the same as the previous scenario.Figures 4 and 5 show the secondary system's WSR and UWSR respectively through different number of subcarriers.In this scenario the improvement in WSR for 64 M  subcarriers is 9%.When the relay node is between the secondary source and destination, in Figure 6 we have shown the WSR of secondary system by varying the total available power budget P for 64 M  subcarriers.And Figure 7 shows   the WSR of secondary system by varying the tolerable interference for subcarriers and fixed total power budget [Watts].Both Figures 6 and 7 show that the JSPPA method outperforms the FSP.

Conclusion
In this study we consider an OFDM based cooperative cognitive network with a pair of Source-Destination nodes as the primary user (PU), and a pair of Source-Destination nodes-which is assisted with a relay-as the secondary user SU.In a shared spectrum underlay paradigm, downlink transmission of the secondary system using a half-duplex relay with DF method under total power constraint is considered.Under defined constraints set, a joint subcarrier pairing and power allocation is proposed for SU to maximize its WSR.The problem is transformed to a convex optimization problem using continuous relaxation method, and solved in the dual domain.In large enough number of subcarriers (M) the solution to dual problem will tend to the optimal WSR solution.An algorithm to achieve feasible solutions is used based on the dual optimization results.The algorithm has two parts, the first part is to find the optimum values for optimization variables which are power allocation and subcarrier pair scheme, and second part is to minimize the Lagrangian function over Lagrange dual variables.Through extensive simulations, the spectrum utilization of the proposed approach is benchmarked, and compared with the existing ones.Interestingly it is shown that the proposed method improves the weighted sum rate of SU, and the simulations endorse our mathematical proofs.Simulations show that in best conditions the proposed method JSPPA2 improves the WSR of SU by 17% in respect to the fixed subcarrier pairing (FSP) method.
Cooperation in multiuser primary and secondary systems with more relay nodes can be considered in future researches by using the same optimization approach.Also using the proposed method for cooperative cognitive networks which obey overlay paradigm, where cognitive relay nodes can help even primary users to enhance their data transmission, can be a subject of interest in future studies.
) can be different.In each iteration the values of , k m and can be updated using the new values obtained for p , k m t  and m  in (23) and (25).As the iteration advances, (26) and (27) will converge and the optimum dual variables will be obtained.(23) rate for selected subcarrier pair SP(k,m) including cost of selecting subcarrier m in second time slot ( m  ) and price of power consumption (the last term in , k m X ).Hence, the optimum pair for subcarrier k in first time slot is subcarrier m in second time slot which maximizes , ) We apply our solution to [P2] using Algorithm 1.The problem [P2] has two parts, first part is maximizing the Lagrangian function (21) over variables , value for dual Lagrange variables, and second part is finding the optimum values for the p Assuming  and m  are given we derived , k m and , k m in (23) and (25) respectively, and the interior part of our dual optimization problem is solved.The next step is *

Algorithm 1 .
Joint Subcarrier Pairing and Power Allocation for SU WSR Maximization in Cooperative Cognitive Network 1. Initialize parameters: 1

Figure 2 .
Figure 2. WSR of secondary system, when the relay is between the secondary source and destination.

Figure 3 .
Figure 3. UWSR of secondary system, when the relay is between the secondary source and destination.

Figure 4 .
Figure 4. WSR of secondary system, when the relay is close to the secondary source.

Figure 5 .
Figure 5. UWSR of secondary system, when the relay is close to the secondary source.

Figure 6 .Figure 7 .
Figure 6.WSR of secondary system through different values of , for subcarriers.budget P 64  M