The Scientific World Journal

Volume 2014 (2014), Article ID 380106, 10 pages

http://dx.doi.org/10.1155/2014/380106

Research Article

Joint Subchannel Pairing and Power Control for Cognitive Radio Networks with Amplify-and-Forward Relaying

Yanyan Shen,¹ Shuqiang Wang,² and Zhiming Wei³

¹Department of Mechanical and Biomedical Engineering, City University of Hong Kong, Kowloon, Hong Kong

²Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong Province 518000, China

³Shandong Academy of Agricultural Machinery Sciences, Jinan, Shandong Province 250000, China

Received 2 December 2013; Revised 30 April 2014; Accepted 20 May 2014; Published 15 June 2014

Academic Editor: Juncheng Jia

Copyright © 2014 Yanyan Shen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Dynamic spectrum sharing has drawn intensive attention in cognitive radio networks. The secondary users are allowed to use the available spectrum to transmit data if the interference to the primary users is maintained at a low level. Cooperative transmission for secondary users can reduce the transmission power and thus improve the performance further. We study the joint subchannel pairing and power allocation problem in relay-based cognitive radio networks. The objective is to maximize the sum rate of the secondary user that is helped by an amplify-and-forward relay. The individual power constraints at the source and the relay, the subchannel pairing constraints, and the interference power constraints are considered. The problem under consideration is formulated as a mixed integer programming problem. By the dual decomposition method, a joint optimal subchannel pairing and power allocation algorithm is proposed. To reduce the computational complexity, two suboptimal algorithms are developed. Simulations have been conducted to verify the performance of the proposed algorithms in terms of sum rate and average running time under different conditions.

1. Introduction

With the emergence of wireless services and applications, the demand for spectrum is increasing. However, experimental results have demonstrated that traditional fixed spectrum allocation schemes lead to inefficient spectrum utilization [1]. To improve the spectrum utilization, cognitive radio networks (CRNs) that allow the unlicensed secondary users (SUs) to share the licensed spectrum with the licensed primary users (PUs) have been proposed.

Based on whether the PUs are aware of the existence of the SUs, the dynamic spectrum sharing in CRNs can be divided into two categories. For the first category, the PUs are not aware of the existence of the SUs. The SUs are able to detect the available spectrums and then use them for their own data transmission if the interference to the PUs is not harmful. For the second category, the PUs leverage parts of its spectrums to the SUs in exchange for appropriate remuneration.

Recently, by exploiting the cooperative diversity, the cooperation technique has been incorporated into the CRNs [2, 3]. Based on the above categories for the dynamic spectrum sharing, the cooperation in CRNs can be divided into two categories as well [4]. For the first category, the SUs are located in different areas and thus have different channel conditions. In this case, some SUs may be selected as relays to help forward other SUs’ data [5–9]. It is named as cooperation among SUs. For the second category, the PUs lease some resources (time or frequency) to the SUs; thus the SUs will act as relays to forward the PUs’ data. When the data of the PUs are completely transmitted or the data rate requirements of the PUs are satisfied, then the SUs can transmit their own data [10–12]. We call it cooperation among PUs and SUs. In this paper, we will consider the cooperation among SUs.

The resource allocation problem is significant for CRNs with cooperation among SUs. A typical feature of CRNs, which differs from other kinds of networks, is that the quality of service (QoS) of the PUs should be guaranteed by imposing interference power constraints on the resource allocation problem [13–15]. For the resource allocation problem in relay networks, the problem of power allocation and relay selection has been investigated in [16, 17]. In particular, if wireless relay networks have multiple channels, for example, orthogonal frequency-division multiplexing (OFDM) relay networks, the resource allocation problem including the subchannel allocation and subchannel pairing has also been studied. The existing literature has shown that appropriate subchannel pairing has improved the network performance for not only the amplify-and-forward (AF) relay strategy but also the decode-and-forward (DF) relay strategy [18–22]. However, there have been very few studies devoted to the resource allocation problem with subchannel pairing in relay-based CRNs.

Some studies have investigated the resource allocation problem for the case that secondary relay helps SU’s transmissions in CRNs. Resource allocation in the downlink transmission over a cluster-based cellular network was analyzed in [23], where relays were fixed and supposed to be able to help forward data from the base station to the SUs in its cluster. The formulated problem was an average weighted sum throughput maximization problem; then the dual decomposition method was used to derive an asymptotically optimal algorithm. The relay selection and power allocation problem has been investigated in [5, 6]. The problems under consideration in [5, 6] were formulated as optimization problems under the interference power constraint and the total power constraint; however, the network models were different. Li et al. [5] considered one channel in the model, but Bharadia et al. [6] considered multiple channels. Although there were multiple channels, the channel pairing problem was not considered in [6]. Meanwhile, the power constraint in [6] was the total source and relay power constraint. This constraint can give some insight for the power allocation but is not very suitable for modeling the practical networks where the source and the relay usually have their individual power supplies. Hence, individual source and relay power constraints are more appropriate for modeling the wireless relay networks in reality, and we will study the resource allocation under individual source and relay power constraint in this paper. Soury et al. in [7] studied the subchannel pairing, power allocation, and relay strategy selection (AF and DF selection) problem in CRNs. Although this paper considered subchannel pairing, the optimal subchannel pairing solution was not given. Sidhu et al. in [8] studied the joint subchannel pairing and power allocation problem for AF relay-based CRNs. Our paper is related to [8]. However, it should be noted that the network models in [8] and our paper are very different. To make the network model simple, the direct transmission link from the secondary source node to the secondary destination node is missing in [8], while both the direct transmission link and the relay transmission link are considered in our paper. From this point of view, the network model in [8] can be regarded as a special case of our model. Hence, the network model in our paper is more general than that in [8]. In addition, the proposed algorithms in [8] cannot be directly applied to solve the problem here. In Section 5.2, the performance comparison of the proposed algorithms in this paper and the optimal algorithm in [8] will be analyzed.

In this paper, we will investigate the joint subchannel pairing and power allocation problem in CRNs where a SU is helped by a secondary relay with AF relay strategy. Both the direct transmission link and the relay transmission link are considered. Our objective is to maximize the SU’s sum rate on all subchannels subject to individual power constraints at the source and the relay, the interference power constraints to the PU, and the subchannel pairing constraint. The problem is formulated as a mixed integer programming problem. We jointly optimize the subchannel pairing and power allocation by the dual decomposition method. To reduce the computational complexity further, two suboptimal algorithms are proposed as well. Simulation results show that the two suboptimal algorithms spend much less running time achieving sum rates that are close to the optimal solution under various simulation conditions.

The remainder of the paper is organized as follows. Section 2 describes the system model and problem formulation. The optimal solution is proposed in Section 3. Two suboptimal algorithms with lower computational complexity are proposed in Section 4. Performance evaluation of the proposed algorithms is shown in Section 5, and Section 6 concludes the work.

2. System Model and Problem Formulation

Consider a network setting where a PU and a SU coexist, which is shown in Figure 1. The SU can adopt the available channels that are licensed to the PU for its own data transmission. It is assumed that the total available bandwidth that is licensed to the PU is Hz. The total bandwidth is divided into nonoverlapping subchannels each with Hz bandwidth. Let denote the set of subchannels. Since the SU may cause interference to the PU, relay-based transmission is adopted by the SU to reduce power; thus it can ease the interference to the PU. For simplicity, we will use source, relay, and destination to denote secondary source, secondary relay, and secondary destination in the following.

Figure 1: System model.

The channel gains of subchannel from the source to the destination and from the source to the relay are denoted by and , respectively, and those of subchannel from the relay to the destination are denoted by . , , and are the variances of the additive white Gaussian noises (AWGN) in the corresponding subchannels, respectively. It is assumed that perfect channel state information of each subchannel is known by a central controller. For notational brevity, let , , and . AF relay strategy is used by the relay. It is assumed that the time domain is equally divided into a number of time slots. The AF cooperative transmission is a two-time slot process. In the first time slot, the secondary source uses power to transmit data on subchannel ; the secondary relay and destination can receive signal from it. In the second time slot, the relay amplifies the received signal and uses power to forward it to the destination. At the end of the second time slot, the destination combines the received signals from both the relay and the source using the maximum ratio combiner. To improve the transmission rate further, subchannel pairing is considered. Let denote the subchannel pairing in the two time slots. If subchannel in the first time slot and subchannel in the second time slot are paired, ; otherwise . It is assumed that the signals transmitted on different subchannel pairings are different. The data rate (bps/Hz) for the cooperative transmission on subchannel pairing is denoted by Li et al. [19]

which is not jointly concave with respect to and . To facilitate the analysis in the following, it can be approximated by

Such approximation has also been used in [19–22], and the approximation gap disappears as the signal-to-noise ratio becomes large.

To protect the PU’s QoS, the interference to the primary destination should not be larger than the given threshold and in the first and second time slot, respectively, which can be expressed by

where and and and are the channel gain from the source to the primary destination on subchannel and that from the relay to the primary destination on subchannel , respectively. and are the variances of AWGN on the corresponding subchannels. It should be noted that and might have the same value in practice if the interference tolerance level of the PU is a constant or the required QoSs of the PU are stable in a given time period. Here, we use and to describe a general relay-based CRN model, and the following results in the paper can be easily extended to the case where .

We desire to study the subchannel pairing and power allocation problem to maximize the transmission rate of the SU under the source power constraint, the relay power constraint, the subchannel pairing constraint, and the interference power constraint. The problem under consideration can be formulated as follows:

where and represent the individual source and relay power constraints; and are the total power at the source and the relay, respectively. indicates that the consumed power on each subchannel at the source and the relay should be nonnegative. ~ are the subchannel pairing constraints. reveals that the subchannel pairing indicator is a binary variable. Combined with , indicates that, for any subchannel in the second time slot, at most one subchannel in the first time slot can be paired with it, and indicates that, for any subchannel in the first time slot, at most one subchannel in the second time slot can be paired with it. is the interference power constraints to the primary destination.

3. Optimal Solution

Problem (5) is a mixed integer programming problem; it is generally difficult to obtain the optimal solution. Fortunately, when the number of subchannels goes to infinity, the dual gap between the primal problem and the dual problem goes to zero [24]; thus the dual decomposition method can be used to solve this problem.

3.1. Dual Problem

By introducing the nonnegative Lagrange multipliers to the constraints , , and , the partial Lagrange function can be expressed as

Thus the dual problem is

where is the dual function and is given by

3.2. Solution to the Dual Problem

The dual function is not differential due to the discontinuity of the subchannel pairing variable , and hence its gradient does not exist. Instead, the subgradient method can be adopted. The subgradient of is

where , , and are the solution to (8). We will discuss how to solve (8) in the following. Based on the subgradient of , the update of the dual variable is given by

where is the iteration number, is the step size at iteration , and . If satisfies the diminishing step size rule, then by using the subgradient method in (10), is guaranteed to converge to the optimal [25].

Next, we will deduce the optimal subchannel pairing and power allocation, which is the solution to problem (8). Through some algebra operations, can be rewritten as

To facilitate the calculation of the optimal power and subchannel pairing, in (11) can be deleted without impact on deriving the optimal power allocation, and thus, can be expressed by

where

Taking the first-order derivatives of with respect to and and making them equal 0, then we get

With (14) and (15), we obtain that at the optimal point

where

To make sure , . Substituting (16) into (15), then we get the optimal and shown as

where

Substituting and into (13), we get

Since the aim of (8) is to get the maximal of , intuitively, the subchannel pairing solution can be determined by

It should be noted that although (22) can guarantee that the constraint holds, it cannot guarantee that the constraint simultaneously holds. That is to say, the obtained from (22) may correspond to a situation where more than one subchannel in the first time slot might be paired with the same subchannel in the second time slot, for example, the subchannel pairing solutions and . As a result, the subchannel pairing obtained from (22) may be infeasible. To guarantee the subchannel pairing constraints , , and are satisfied simultaneously, one naive method is the exhaustive search. However, its computational complexity is , which is prohibitive from practical implementation, especially when the number of subchannels is large. In order to reduce the computational complexity without loss of optimality, we will employ the Hungarian method to solve the subchannel pairing problem.

Define matrix ; then the subchannel pairing problem is equivalent to an assignment problem. The rule is to pick up one and only one entry from each row and each column of to make the sum of the picked up entries maximal. The Hungarian method with the computational complexity can solve this kind of problem efficiently. The Hungarian method for subchannel pairing is described in Algorithm 1. The whole algorithm is shown in Algorithm 2, which is a joint optimal subchannel pairing and power allocation algorithm.

(1) Let matrix , and the following operations are for .

(2) Subtract the smallest entry in each row from all the entries of its row.

(3) Subtract the smallest entry in each column from all the entries of its column.

(4) Draw lines through appropriate rows and columns so that all the zero entries of matrix are

      covered and the minimum number of such lines is used.

(5) Test for optimality: (i) If the minimum number of covering lines is , an optimal assignment

      of zeros is possible and this algorithm finishes. (ii) If the minimum number of covering lines

      is less than , an optimal assignment of zeros is not yet possible, and then go to Step  6.

(6) Determine the smallest entry not covered by any line. Subtract this entry from each uncovered

      row, and then add it to each covered column. Return to Step  3.

Algorithm 1: Hungarian method [27] for subchannel pairing.

(1) Initiate the value of ,

(2) for   to the maximum number of iterations, do

(3)    Calculate transmission power and by (18) and (19), respectively.

(4)    Substitute and into in (13).

(5)    Determine the optimal value of , by the Hungarian method.

(6)    Update by (10).

(7) end for

Algorithm 2: Optimal solution to (5).

Although Algorithm 2 can get the solution to problem (5), it is observed from the simulation results that the source power constraints with (3) and the relay power constraint with (4) are hard to be satisfied with equality. Because (5) is a maximization problem with an objective function that is increasing with respect to and , thus, at the optimal point, at least one of the source power constraints and (3) should be met strictly. And it also holds for the relay power constraints and (4). Hence, the power obtained from Algorithm 2 should be scaled until at least one of the source power constraints and (3)and at least one of the relay power constraints and (4) are satisfied with equality.

4. Suboptimal Solutions

In this section, two suboptimal algorithms with low computational complexity will be proposed.

4.1. Suboptimal Algorithm  1

Suboptimal algorithm 1 is a joint fixed subchannel pairing and optimal power allocation algorithm. The subchannel pairing is predetermined; for example, subchannel in the first time slot is paired with subchannel in the second time slot, and thus, and , and . When the subchannel pairing is determined, the power allocation problem becomes much simpler than problem (6). And the dual decomposition method used in the previous section can be adopted here. It should especially be pointed that from the subchannel pairings and , and , it is known that and , and , and thus only and , need to be calculated in each iteration. The total algorithm is shown in Algorithm 3.

(1) Let , , and and ,

(2) Initiate the value of ,

(3) for   to the maximum number of iterations, do

(4)    Calculate transmission power and by (18) and (19), respectively.

(5)    Update by (10).

(6) end for

Algorithm 3: Suboptimal algorithm  1.

Compared with Algorithm 2, this algorithm has similar iteration process as Algorithm 2 expect for two aspects. First, this algorithm does not need to calculate the subchannel pairing by the Hungarian method, but Algorithm 2 needs this. Second, this algorithm only needs to calculate the power and , , but Algorithm 2 needs to calculate and , .

As in Algorithm 2, the power obtained by this algorithm cannot make the power constraints , , and in (5) be satisfied with equality. Therefore, the same power scaling as Algorithm 2 will be used at the end of this algorithm.

4.2. Suboptimal Algorithm  2

Suboptimal algorithm  2 is a joint optimal subchannel pairing and equal power allocation algorithm. The optimal subchannel pairing is calculated by the Hungarian method. For equal power allocation, both total power constraints and interference power constrains should be satisfied. This is achieved by first ignoring the interference power constraints and allocating the total power equally on each subchannel and then checking the interference power constraints and conducting some operations until the interference power constraints are satisfied. The iteration process of this algorithm is shown as follows, and the entire suboptimal algorithm  2 is described in Algorithm 4.

(1) Equally distribute the total power on each subchannel, that is, .

(2) Do optimal subchannel pairing by the Hungarian method.

(3) Check whether the constraints (3) and (4) are satisfied:

(4) if constraints (3) and (4) are satisfied then

(5)    obtain the final solution, and the algorithm terminates;

(6) else if constraint (3) is violated but constraint (4) is satisfied then

(7)    ;

(8) else if constraint (3) is satisfied but constraint (4) is violated then

(9)    ;

(10) else if both constraints (3) and (4) are violated then

(11)  ;

(12) end if

(13) Go to step 1.

Algorithm 4: Suboptimal algorithm  2.

Step 1. Without the consideration of the interference power constraints (3) and (4), the total power is equally distributed to each subchannel; that is, , , , . Then the Hungarian method is used for the subchannel pairing.

Step 2. Based on the allocated power and subchannel pairing, the interference power constraints (3) and (4) will be verified. If the constraints (3) and (4) are satisfied, stop the iteration of the algorithm and now the solution is the final solution; otherwise, reduce the total power consumption on all the subchannels. For the reduction of the total power consumption, three cases should be considered. (a) If (3) is violated but (4) is satisfied, then reduce the total source power by a small constant   . (b) If (3) is satisfied but (4) is violated, then reduce the total relay power by a small constant . (c) If (3) and (4) are violated, then reduce both the source and relay power and by a small constant .

Step 3. Reallocate the source and the relay power equally on all the subchannels and find the optimal subchannel pairing by the Hungarian method; then go to Step 2.

The comparisons of the optimal algorithm and two suboptimal algorithms are summarized in Table 1. From Table 1, it is evident that both suboptimal algorithms do not fully exploit the potential of different channel conditions for network performance improvement since they carry out fixed subchannel pairing or equal power distribution on each subchannel that does not take into account different channel conditions, and thus they will result in sum rate reduction to some extent.

Table 1: Comparison of the proposed optimal algorithm and suboptimal algorithms.

Algorithm
Optimal algorithm
Suboptimal Algorithm  1
Suboptimal Algorithm  2

Resource allocation
Joint optimal subchannel pairing and power allocation
Given fixed subchannel pairing
Optimal subchannel pairing

Optimal power allocation
Equal power distribution

5. Performance Evaluation

5.1. Complexity Comparison

For the optimal algorithm (Algorithm 2), in each iteration it needs operations to compute and , respectively. Then, it needs operations to calculate . The complexity of the Hungarian method is . Hence, the total complexity of Algorithm 2 is , where is the iteration number of the subgradient method. It has been shown in [24] that is a polynomial function with respect to the number of the dual variables that is four in Algorithm 2.

For suboptimal algorithm  1, it needs to compute and in each iteration; thus, it requires operations. The subgradient algorithm needs iterations to converge; hence its total complexity is .

The suboptimal algorithm  2 needs iterations to calculate the power distribution and subchannel pairing, where is related to the values of , , , , and . It needs operations to calculate the power distribution on each subchannel. The Hungarian method with complexity is used to calculate the subchannel pairing. Therefore, its overall complexity is .

It should be noted that the computational complexity of the optimal algorithm is higher than the suboptimal algorithm  1. Although it is difficult to compare the value of and , simulations in Section 5.2 reveal that suboptimal algorithm  2 has a lower complexity than suboptimal algorithm  1 in most scenarios.

5.2. Numerical Results

In this subsection, simulation results are presented to illustrate the performance of the proposed algorithms compared with the joint optimal algorithm in [8], which is indicated by optimal algorithm in [8] in the following figures and tables.

Consider a simulation model shown in Figure 2, where the source and the destination of the SU are located at () and (, 0). The coordinate of the relay is (). The primary destination is located at (). The total bandwidth is  Hz. The channel gain in any transmission pair contains a large-scale Rayleigh fading component and a large-scale path loss component with path loss factor 4. The small-scale fading is modeled as a frequency selective channel consisting of six independent Rayleigh multipaths. Each multipath component is modeled by Clarke’s flat fading model. It is assumed that the power delay profile is exponentially decaying with , where is the multipath index. Therefore, the relative power of the six multipath components is dB [26]. The total powers at the source and the relay are assumed to be equal for simplicity. For the suboptimal algorithm  2, we set in the following simulations. To show the performance of the algorithms in a fair and efficient manner, all the interested data in the simulations are averaged over 1000 independent channel realizations.

Figure 2: Simulation model.

First, simulations have been conducted to investigate the impact of the location of the relay on the sum rate. The parameters are set to be m,  W, W, and . Figure 3 illustrates the sum rates of the four algorithms versus the location of the relay. As the increase of , the distance between the source and the relay becomes longer, and the distance between the relay and the destination becomes shorter. Figure 3 shows that all the algorithms reach the highest sum rate when the relay is located at the center between the source and the destination. The sum rate of the proposed optimal algorithm is higher than those of the other three algorithms under nearly all the values of . When the relay is close to the source () or the destination (), the gap between the sum rates of the proposed optimal algorithm and the proposed suboptimal algorithm  1 is very small, that is, less than . The optimal algorithm in [8] always achieves the lowest sum rate compared with the proposed three algorithms in this paper. This is mainly because the optimal algorithm in [8] does not consider the direct transmission link while the proposed algorithms consider both the direct transmission link and the relay transmission link.

Figure 3: Sum rate versus the location of the relay.

Figure 4 illustrates the sum rate comparison of the algorithms when the number of the subchannels changes. We set m, ,  W, and  W. As is evident from Figure 4, the sum rates of all the algorithms are increasing as the number of subchannels increases; that is because when the number of subchannels is large, there are big chances to allocate more power to the subchannel with higher channel gain, and thus, the total rate increases. Under different numbers of subchannels, the optimal algorithm always obtains a higher sum rate than the suboptimal algorithms and the optimal algorithm in [8], but the sum rates of the proposed two suboptimal algorithms are always very close to the optimal solutions; that is, the suboptimal algorithms  1 and 2 achieve more than and of the corresponding optimal solution, respectively. The sum rate of the optimal algorithm in [8] is always lower than that of the proposed optimal algorithm under different numbers of subchannels.

Figure 4: Sum rate versus the number of subchannels.

The computational complexity of the proposed algorithms is highly related to the number of subchannels as we analyzed in Section 5.1. To compare the computational complexity, we use Table 2 to show the average running time comparison of the algorithms when the number of subchannels changes. The simulation parameters of Table 2 have the same values as those in Figure 4. It is apparent from Table 2 that the proposed optimal algorithm and the optimal algorithm in [8] always need much longer average running time to obtain the final solution compared with the two suboptimal algorithms, but the running times of the proposed optimal algorithm and the optimal algorithm in [8] are also different. When the number of subchannels is 18, the running time of the optimal algorithm in [8] is nearly twice as much as that of the proposed optimal algorithm. Combined with the sum rate comparison in Figure 4, it is indicated that when the number of subchannels is larger than 10, the proposed optimal algorithm achieves a higher sum rate but with a much less running time compared with the optimal algorithm in [8]. In most cases, suboptimal algorithm  2 spends a little less time than suboptimal algorithm  1. From the simulation results of Figure 4 and Table 2, we conclude that the two proposed suboptimal algorithms can achieve near optimal solution with a significant average running time reduction.

Table 2: Average running time (in seconds) comparison for the four algorithms when the number of subchannels changes,  m, ,   W, and  W.

Suboptimal algorithm 1
Suboptimal algorithm 2
Optimal algorithm
Optimal algorithm in [8]

Figure 5 shows the sum rate of the algorithms versus the interference power threshold . The parameters are m, ,  W, and . It is apparent from Figure 5 that the sum rate is largely impacted by the value of . When is very small (), the sum rates of the four algorithms increase dramatically as, the increase of , that is, because the interference power constraints are the dominating constraints in this case. As increases, the total power constraints gradually become the dominating constraints, and thus, the sum rates keep changing in a small region. The proposed optimal algorithm achieves a higher sum rate compared with the other three algorithms under different values of . The ratios of the average sum rate achieved by suboptimal algorithms 1 and 2 to that achieved by the optimal algorithm are more than and under different values of , respectively.

Figure 5: Sum rate versus .

The sum rate versus the total power is illustrated in Figure 6. The parameters are set to be m,  W, , and . It is seen from Figure 6 that the two suboptimal algorithms 1 and 2 achieve solutions that are more than and of the optimal solution for different values of , respectively. The sum rate of the optimal algorithm [8] is always nearly the same as the suboptimal algorithm  2 when  W, and it becomes much smaller than that of suboptimal algorithm  2 when becomes larger. The sum rate of all these algorithms increases as increases. However, the sum rate increase rate gradually becomes slower as increases. That is because, although increases, the interference power constraints restrict the actual total power consumption on all subchannels.

Figure 6: Sum rate versus the total power.

6. Conclusion

In this work, we have considered the problem of subchannel pairing and power allocation in CRNs with the AF relay. The network model contains a PU and a secondary source-destination pair assisted by a secondary relay. The objective is to maximize the rate of the SU under the individual power constraints at the source and the relay and the subchannel pairing constraints. In order to protect the PU, the interference power constraints to the PU are also considered. The problem under consideration is formulated as a mixed integer programming problem. Via the dual decomposition method, we propose an algorithm that jointly optimizes the subchannel pairing and power allocation. To reduce the computational complexity, two suboptimal algorithms are proposed as well. The simulation result has shown that when the relay is located at the center between the source and the destination, the three proposed algorithms obtain the highest sum rate. It has also shown that the proposed suboptimal algorithms significantly reduce the average running time and achieve sum rates that are close to the optimal solution under different conditions. Therefore, the proposed suboptimal algorithms are very suitable for practical implementations where sum rate reduction can be tolerant to some extent. In the future, we will analyze the theoretical approximate ratio of the suboptimal algorithms to the optimal one and will investigate the performance of our proposed algorithms in the presence of imperfect channel state information.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the editors and the anonymous reviewers, whose invaluable comments have helped to improve the quality of the paper.

References

D. Chen, S. Yin, Q. Zhang, M. Liu, and S. Li, “Mining spectrum usage data: a large-scale spectrum measurement study,” in Proceedings of the 15th Annual ACM International Conference on Mobile Computing and Networking (MobiCom '09), pp. 13–24, Beijing, China, September 2009.

K. B. Letaief and W. Zhang, “Cooperative communications for cognitive radio networks,” Proceedings of the IEEE, vol. 97, no. 5, pp. 878–893, 2009.

Q. Zhang, J. Jia, and J. Zhang, “Cooperative relay to improve diversity in cognitive radio networks,” IEEE Communications Magazine, vol. 47, no. 2, pp. 111–117, 2009.

O. Simeone, J. Gambini, Y. Bar-Ness, and U. Spagnolini, “Cooperation and cognitive radio,” in Proceedings of the IEEE International Conference on Communications (ICC '07), pp. 6511–6515, Glasgow, Scotland, June 2007.

L. Li, X. Zhou, H. Xu, G. Y. Li, D. Wang, and A. Soong, “Simplified relay selection and power allocation in cooperative cognitive radio systems,” IEEE Transactions on Wireless Communications, vol. 10, no. 1, pp. 33–36, 2011.

D. Bharadia, G. Bansal, P. Kaligineedi, and V. K. Bhargava, “Relay and power allocation schemes for OFDM-based cognitive radio systems,” IEEE Transactions on Wireless Communications, vol. 10, no. 9, pp. 2812–2817, 2011.

H. Soury, F. Bader, M. Shaat, and M. Alouini, “Joint sub-carrier pairing and resource allocation for cognitive networks with adaptive relaying,” EURASIP Journal on Wireless Communications and Networking, vol. 259, pp. 1–15, 2013.

G. Sidhu, F. Gao, W. Wang, and W. Chen, “Resource allocation in relay-aided OFDM cognitive radio networks,” IEEE Transactions on Vehicular Technology, vol. 62, no. 8, pp. 3700–3710, 2013.

S. Wang, M. Ge, and C. Wang, “Efficient resource allocation for cognitive radio networks with cooperative relays,” IEEE Journal on Selected Areas in Communications, vol. 31, no. 11, pp. 2432–2441, 2013.

K. Khalil, M. Karaca, O. Ercetin, and E. Ekici, “Optimal scheduling in cooperate-to-join cognitive radio networks,” in Proceedings of the IEEE INFOCOM, pp. 3002–3010, Shanghai, China, April 2011.

T. Elkourdi and O. Simeone, “Spectrum leasing via cooperation with multiple primary users,” IEEE Transactions on Vehicular Technology, vol. 61, no. 2, pp. 820–825, 2012.

S. M. M. Toroujeni, S. M.-S. Sadough, and S. A. Ghorashi, “Spectrum leasing for OFDM-based cognitive radio networks,” IEEE Transactions on Vehicular Technology, vol. 62, no. 5, pp. 2131–2139, 2013.

A. T. Hoang and Y.-C. Liang, “Downlink channel assignment and power control for cognitive radio networks,” IEEE Transactions on Wireless Communications, vol. 7, no. 8, pp. 3106–3117, 2008.

L. Zhang, Y. Xin, Y.-C. Liang, and H. V. Poor, “Cognitive multiple access channels: optimal power allocation for weighted sum rate maximization,” IEEE Transactions on Communications, vol. 57, no. 9, pp. 2754–2762, 2009.

D. T. Ngo and T. Le-Ngoc, “Distributed resource allocation for cognitive radio networks with spectrum-sharing constraints,” IEEE Transactions on Vehicular Technology, vol. 60, no. 7, pp. 3436–3449, 2011.

S. Kadloor and R. Adve, “Relay selection and power allocation in cooperative cellular networks,” IEEE Transactions on Wireless Communications, vol. 9, no. 5, pp. 1676–1685, 2010.

K. Vardhe, D. Reynolds, and B. Woerner, “Joint power allocation and relay selection for multiuser cooperative communication,” IEEE Transactions on Wireless Communications, vol. 9, no. 4, pp. 1255–1260, 2010.

W. Ying, Q. Xin-chun, W. Tong, and B.-L. Liu, “Power allocation and subcarrier pairing algorithm for regenerative OFDM relay system,” in Proceedings of the IEEE 65th Vehicular Technology Conference (VTC '07), pp. 2727–2731, Dublin, Ireland, April 2007.

Y. Li, W. Wang, J. Kong, and M. Peng, “Subcarrier pairing for amplify-and-forward and decode-and-forward OFDM relay links,” IEEE Communications Letters, vol. 13, no. 4, pp. 209–211, 2009.

C.-N. Hsu, H.-J. Su, and P.-H. Lin, “Joint subcarrier pairing and power allocation for OFDM transmission with decode-and-forward relaying,” IEEE Transactions on Signal Processing, vol. 59, no. 1, pp. 399–414, 2011.

M. Hajiaghayi, M. Dong, and B. Liang, “Jointly optimal channel pairing and power allocation for multichannel multihop relaying,” IEEE Transactions on Signal Processing, vol. 59, no. 10, pp. 4998–5012, 2011.

W. Dang, M. Tao, H. Mu, and J. Huang, “Subcarrier-pair based resource allocation for cooperative multi-relay OFDM systems,” IEEE Transactions on Wireless Communications, vol. 9, no. 5, pp. 1640–1649, 2010.

R. Wang, V. K. N. Lau, and Y. Cui, “Decentralized fair scheduling in two-hop relay-assisted cognitive OFDMA systems,” IEEE Journal on Selected Topics in Signal Processing, vol. 5, no. 1, pp. 171–181, 2011.

W. Yu and R. Lui, “Dual methods for nonconvex spectrum optimization of multicarrier systems,” IEEE Transactions on Communications, vol. 54, no. 7, pp. 1310–1322, 2006.

S. Boyd, L. Xiao, and A. Mutapcic, “Subgradient methods,” Notes for EE392o, Stanford University, 2003.

Z. Shen, J. G. Andrews, and B. L. Evans, “Adaptive resource allocation in multiuser OFDM systems with proportional rate constraints,” IEEE Transactions on Wireless Communications, vol. 4, no. 6, pp. 2726–2736, 2005.

D. Bruff, “The assignment problem and the hungarian method,” Notes for Math 20, Harvard Univerisity, 2005.