Sensors 2013, 13, 10306-10332; doi:10.3390/s130810306

OPEN ACCESS

sensors

ISSN 1424-8220

www.mdpi.com/j ournal/sensors

Article

A Subcarrier-Pair Based Resource Allocation Scheme Using Proportional Fairness for Cooperative OFDM-Based Cognitive Radio Networks

Yongtao Ma, Liuji Zhou * and Kaihua Liu

School of Electronic Information Engineering, Tianjin University, 92 Weijin Road, Nankai District,

Tianjin 300072, China; E-Mails: mayongtao@tju.edu.cn (Y.M.); liukaihua@tju.edu.cn (K.L.)

* Author to whom correspondence should be addressed; E-Mail: honglin929@tju.edu.cn;

Tel.: +86-158-2245-5869.

Received: 28 June 2013; in revisedform: 5 August 2013 /Accepted: 6 August 2013 /

Published: 9 August 2013

Abstract: The paper presents a joint subcarrier-pair based resource allocation algorithm in order to improve the efficiency and fairness of cooperative multiuser orthogonal frequency division multiplexing (MU-OFDM) cognitive radio (CR) systems. A communication model where one source node communicates with one destination node assisted by one half-duplex decode-and-forward (DF) relay is considered in the paper. An interference-limited environment is considered, with the constraint of transmitted sum-power over all channels and aggregate average interference towards multiple primary users (PUs). The proposed resource allocation algorithm is capable of maximizing both the system transmission efficiency and fairness among secondary users (SUs). Besides, the proposed algorithm can also keep the interference introduced to the PU bands below a threshold. A proportional fairness constraint is used to assure that each SU can achieve a required data rate, with quality of service guarantees. Moreover, we extend the analysis to the scenario where each cooperative SU has no channel state information (CSI) about non-adjacent links. We analyzed the throughput and fairness tradeoff in CR system. A detailed analysis of the performance of the proposed algorithm is presented with the simulation results.

Keywords: cognitive radio; cooperative communication; resource allocation; proportional fairness; spectrum sharing

1. Introduction

Cognitive radio technology (CR) has been proposed as a relatively new concept for improving the overall utilization of spectrum bands. This promising technology can allow the unlicensed secondary users (SUs, also referred to as CR users or CRUs) to access those frequency bands which are not currently being used by licensed primary users (PUs) in a given geographical area [1,2]. Cooperative communication technology [3] allows network nodes with single antennas to use other network nodes' antennas to transmit data, which can generate a virtual multiple-input multiple-output (MIMO) system. Cooperative spectrum sensing is a viable sensing technique to enhance spectral utilization efficiency of secondary users while ensuring the quality of service (QoS) of primary users [4]. In a CR network, SUs are allowed to transmit over the frequency bands of PUs as long as the resulting aggregate interference is kept below a certain threshold. This threshold is known as interference temperature constraint or interference power constraint [1]. As SUs can design power and subcarrier allocation strategies subject to such interference power constraints, the interference introduced to PUs is effectively controlled. A great deal of resource allocation algorithms and interference control strategies has been proposed for spectrum-sharing CR networks. For example, the optimal power allocation strategies to maximize the transmitted data rate of the secondary user with an effective protection of the primary user were studied in [5,6] for spectrum-sharing CR networks.

Orthogonal frequency division multiplexing (OFDM) is an attractive modulation scheme for users in a CR system due to its flexibility in allocating resources among SUs. Since both SUs and PUs may exist in side-by-side bands, yet have different access technologies, mutual interference is the limiting factor for the performance of both networks. Thus, using of the classical subcarrier allocation and power loading algorithms, such as uniform power but variable rate and water-filling algorithms maximizing the transmission capacity of an OFDM-based conventional wireless network may result in higher mutual interference in the PUs' band. There is only one group of users in such a wireless network, i.e., PUs, for a CR system.

According to the latest literature on resource allocation in cooperative communication [7-17], the relay users in the system do not transmit their own data and merely help other non-relay users transmit data. In some wireless applications such as cellular networks, however, each user has its own data to transmit so that it should allocate its total constrained power and subcarriers properly in transmitting its own data and relaying other users' data [18,19]. Tourki [19] focused on efficiency issues by studying how to maximize the total transmitted data rate in non-orthogonal amplify-and-forward (AF) cooperative scheme, which ignores the fairness among the cooperative users. According to [20,21], equal power allocation (EPA) among subcarriers was proposed to separate the user selection from the power of subcarrier. With EPA, the EPA-PRG (proportional rate greedy) [22] algorithm is proposed to maximize the system throughput while keeping the fairness. However, cooperative transmission technology isn't applied in this algorithm. In [23], a linear water-filling scheme (LWF-PI) was proposed. This algorithm maximized the overall transmitted data rate of the CR system while keeping the interference introduced to the PU bands below a threshold. However, the fairness among users was ignored. Chandrashekar et al. [24] proposed an algorithm which is capable of maximizing the total transmitted data rate and achieving a high proportional fairness index. However, this algorithm cannot be applied to the CR network where we must adjust the interference introduced to the PU bands below

a threshold. Tan [25] proposed a joint subcarrier and power algorithm based on Blotto games. This algorithm can achieve a good trade-off performance between fairness and efficiency in OFDMA-based cognitive radio network (CRN), but it cannot obtain the effectiveness of multiuser diversity for the SUs without ability to generate a virtual MIMO system.

A novel scheme was presented in [26] for the allocation of subcarriers, rates, and power in orthogonal frequency-division multiple-access (OFDMA) networks. The resource-allocation problem was solved by decomposing it into a hierarchy of sub-problems. A joint subcarrier and power allocation algorithm was presented in [27] for cooperative MU-OFDM CR systems. In [28], a survey of resource allocation and scheduling schemes in OFDMA wireless networks was presented. Nader et al. in [29] considered the practical case in which only partial CSI for the wireless channel between the secondary base station and SUs is available at the secondary base station. They formulated the resource allocation problem in the secondary network as an optimization problem in which the objective was to maximize the weighted sum rate of the secondary users. A novel sub-channel and transmission power allocation scheme was proposed in [30] for multi-cell OFDMA networks with CR functionality. Tianxiang et al. in [31] discussed optimization over the relay assignment, subcarrier allocation, per node power control, and heterogeneous quality-of-service (QoS) provisioning. Sabit et al. in [32] investigated the performance of an OFDM-based CR spectrum sharing communication system that assumed random allocation and absence of the PU channel occupation information. Hong Xu et al. in [33] formulated a unifying optimization framework based on Nash bargaining solutions to fairly and efficiently allocate resources between primary and secondary networks, in both decentralized and centralized settings. As the optimal resource allocation scheme was highly complex, G. B. et al. [34] proposed a low complexity suboptimal subcarrier and power allocation scheme. They also proposed a suboptimal subcarrier allocation scheme that can guarantee a certain level of fairness among CR users. Naeem et al. introduced in [35] a hybrid heuristic algorithm for the relay assignment and power allocation problem which is a non-convex mixed-integer non-linear optimization problem, and this problem is generally non-deterministic polynomial-time (NP)-hard.

In this paper, a joint subcarrier-pair based resource allocation algorithm in order to improve both efficiency and fairness index is presented first. The definition of fairness is borrowed from the networking literature. In contrast with [36], where large channel fluctuations are intentionally created with "dumb" antennas for long-term proportional fairness resource allocation, this paper proposes a subcarrier-pair based resource allocation algorithm to maintain proportional rates among SUs for each channel realization, which ensures the rates of different SUs to be proportional in any time scale of interest. By formulating the resource allocation and pairing problem in this way, it will be shown that a high transmitted data rate for all SUs (even those with poor channel gains) can be achieved with low computational complexity. Moreover, we extend the analysis to the case in which each SU can only have access to CSI of its adjacent links. This is a more realistic scenario when network nodes are mobile and the timely CSI cannot be exchanged between cooperative users. Consequently, each user can only have access to statistical CSI of non-adjacent links. It is shown that the system performance deteriorates due to limited CSI but still outperforms that of equal power allocation scheme. The key contributions of this work are:

1. It is considered that SUs need to transmit their own data directly to the destination, and in the next phase they also help their partner forward the data received in previous phase to the destination. Simulation results show that in the same situations the system transmitted data rate by proposed algorithm is the highest than that by LWF-PI algorithm [23], EPA algorithm [21] and the Optimal Scheme [37].

2. The proposed subcarrier-pair based resource allocation algorithm ensures the rates of different SUs to be proportional in any time scale of interest, simulation results shown that a high transmitted data rate for all SUs (even those with poor channel gains) can be achieved.

3. It is considered that SU has no CSI about non-adjacent link. In this case, we take full advantage of the statistical information of the non-adjacent links.

Notation: In this paper, a variable with "underline" denotes the temporary optimal values within each iteration process, the "double underline" «denotes the optimal value, and the "bar" i denotes the statistical average value. E ( « ) denotes the expectation operator, and ■ denotes the optimal value when only partial CSI can be obtained by SUs.

2. System Model and Problem Formulation

We consider a hybrid network consisting of a primary network (PRN) and a cognitive radio network (CRN) as shown in Figure 1. The CRN consists of a CR access point (AP) and 2K SUs. The PRN and CRN co-exist within the same geographical area. The access mechanism/modulation format in SUs' band is OFDM. Our focus is mainly on the uplink radio resource allocation in the CRN. The SUs are trying to find the opportunity to access to the AP.

According to [37,38], we also consider that the frequency bands of bandwidth B1, B2, ..., BL which have been occupied by L PUs are sensed by the CR system and known to SU transmitters. Every two SUs form a cooperative partner and they are relay node for each other. As shown in Figure 2, the kth (1 < k < K) cooperative partner consists of two SU transmitters, k1 and k2. As is assumed in [23,37,38], we consider the same side-by-side CR radio access model. The unoccupied bandwidth sensed by SUs

Figure 1. A cooperative MU-OFDM CR uplink system.

for opportunistic spectrum access is located on each side of L PU bands as shown in Figure 3. The available bandwidth for CR transmission is divided into N subcarriers based on OFDM system. It is considered that the access mechanism/modulation format in PUs' band is not known to the CR system and the bandwidth for each CR subcarriers is A/Hz. Some symbols are shown in Table 1.

Figure 2. Model for cooperative transmission.

Figure 3. spectrum access model of cognitive radio system.

In general, there are three instantaneous fading gains in the uplink transmission scenario shown in Figure 1:

(1) The gains between the SU's transmitter and SU's receiver or AP for the nth subcarrier denoted as h^, h^'l, respectively.

(2) The gains between the SU's transmitter and Ith PU's receiver, denoted asht

sp,n ki, pl '

(3) The gains between the Ith PU's transmitter and the SU's receiver or AP, denoted as hp;:kn, hp;:on, respectively.

The channel gains are modeled as independent zero-mean complex Gaussian random variables, where ki denotes the ith SU in kth cooperation partner and pl denotes the lth PU band. According to [39], it is considered that these instantaneous fading gains are perfectly known at the SU's transmitter. Specifically, we assume that the SU's receiver can estimate channel gains h^^'^.y and Zip ^ and report to the CR transmitter. In Section 3.2, we will study the case where the instantaneous fading gains of the non-adjacent links are not perfectly known at the SU transmitter but the statistical information of the non-adjacent links are known at the SU transmitter. Moreover, it is assumed that primary receiver can estimate the channel which is reported to the SU transmitter through a

common control channel.

Table 1. Table of Symbols.

Symbol

Definition

Number of cooperative partners in the CR network

Number of subcarriers

Bandwidth of a subcarrier

Length of a slot

i,SP1(n)

Transmission power of k1th SU on subcarrier n in the first frame

i,SP2(n)

Transmission power of k1th SU on subcarrier n in the second frame

i,SP1(n)

Transmission power of k2th SU on subcarrier n in the first frame

i,SP2(n)

hkskn (i * j)

Transmission power of k2th SU on subcarrier n in the second frame

the channel gain of the communication link from the kth SU to the kj th SU user on the nth subcarrier

the channel gain of the communication link from the k'th SU to AP on the nth subcarrier

hsp,n hki, pl

the channel gain of the interference link

from /th PU to k'th SU user receiver on the nth subcarrier

UPs,n "pl,ki

the channel gain of the interference link

from /th PU to k'th SU user receiver on the nth subcarrier

UPs,n hpl,0

the channel gain of the interference link from /th PU to AP receiver on the nth subcarrier

{z®, z®, z?), i = 1,2,3,4}

the additive noises at the corresponding node

i^k'i' ' ¿>0' '' ~the interference introduced by the PUs into corresponding node

2.1. Cooperative Transmission among SUs

The scenario of a three-node DF diversity model is considered, where one source communicates with one destination assisted by one half-duplex relay, as shown in Figure 2. One transmission period is divided into two consecutive frames. Communication takes place in two phases (listening phase T1 and relaying phase T2, the definition is according to the working state of relay user) for each frame. The power allocation scheme for kth cooperative partner on subcarrier n is shown in Table 2. The source node broadcasts its signal to relay and AP in T1, whereas the relay and AP listen. The relay decodes the signal and forwards it to AP in T2. It is denoted that the subcarrier n in T1 is pairing with subcarrier SP1(n) in T2 for first frame, and pairing with subcarrier SP2(n) for second frame. In the first frame, k2th SU receives data in this time slot while k1th SU transmits a symbol xk1(t) with power level p^^on nth subcarrier in T1. The symbol is received by node 0 (AP) and overheard by

k2th SU as:

^,0 (t ) = (t ) + z01) (t ) + ^0(1) (t )

I--(1)

yk1k 2 (t ) = h^PkuPKn) *« (t ) + zS (t ) + (t )

Table 2. Power allocation scheme for kth cooperative partner on subcarrier n.

T1 in First Frame T2 in First Frame T1 in Second Frame T2 in Second Frame

SUki pn,SP1(n) Pk1,1 0 0 pn,SP2(n)

SUK 0 pn,SP1(n) Pk 2,1 pn,SP2(n) Pk 2,2 0

During this interval, the k2th SU decodes its overheard signal as x*(t) and transmits it to the AP on SP1(n) subcarrier in T2 with the power level P^P'^ . Then the AP receives the signal as:

A 2,0 (t) = hT-yPk?iP1(n) % (t) + ^02) (t) + ^ (t) (2)

In the second frame, the roles of kith SU and k2th SU are reversed. Similarly, k2th SU transmits a symbol xk2(t)with power level P"22P2(")on the nth subcarrier in T1. The symbol is received by node 0

(AP) and overheard by kith SU as:

i — Uss," I p«,SP2(n),

vk2 ( ) ^ z0 ( ) "^0 (l)

^k 2,0 (t)=h;rypS2(-) Xk 2 (t)+z03)(t)+^03)(t)

Jk 2,ki (t)=KT^Pg!2(n) Xk 2 (t)+zki (t)+^k(3) (t)

In T2 of second frame, the AP node receives the noisy signal which is relayed by kith SU with the power levelPk"f2(n), i.e.:

W0 = h^y[P^xl2(t) + z04)(t) + ^0(4)(t) (4)

2.2. Mutual Interference between PUBands and CR Users

In the MU-OFDM CR system, due to the coexistence of PUs and SUs in side by side bands, it is necessary to consider the mutual interference between PUs and SUs. There are two types of interference in the system. One is introduced by the PUs into the SUs band, and the other is introduced by the SUs into the PUs' band. In what follows, we provide brief description and mathematical models for interference between SUs and PUs.

2.2.1. The Interference Introduced into PUs by SUs

CR interference is introduced into the PU spectrum by CR out-of-band (OOB) emissions. OOB emissions arise as a result of transmit pulse shaping such that a portion of the CR radiated power in a vacant subcarrier is leaked into neighboring bands occupied by the PUs. According to [23], the interference factor which is the integration of the power density spectrum of the nth subcarrier across the lth PU band, and can be written as:

Si Pl (di) =| |2 Ts r +Bf (^^f2)2 df (5)

Jdni —z % fTs

where Ts denotes the symbol duration, dnl denotes the distance in frequency between the nth subcarriers of SU band and 1th PU band, and Bl represents occupied bandwidth by 1th PU. It can be

shown from Equation (5) that the interference to PU band is related to the distance between SU band and PU band.

2.2.2. The Interference Introduced into SUs by PUs

The interference introduced into kith SU and AP node transmitting in nth subcarrier by Tth PU can be denoted as Jpli0, respectively. According to [37], the interference value Jpli0 can be

written as:

Jnpl„(dnl,Ppu) =\h% |2 C+f E{In(w^U(eiw))}dw

J"P:Mi, Ppu ) =1 h^ |2 \ \2 E{In (wAu (ejw))}dw (6)

J" (d P ) =| hsp,n I2 {dnl + 2

',0 (dnl,1 PU ) =| hpl,0 1 I Af ElJN dnl 2

E{In (w.Au (e'w)} = -L: f $ru (e-X^^M^ df

2nMisin(w - g) / 2

where w represents the frequency normalized to the sampling frequency, E{IN( • )} is the power density spectrum of the PU signal after M-fast Fourier transform (FFT) processing, 0PU (eJw) is the power density spectrum of the PU signal. The PU signal has been taken to be an elliptically filtered white noise process with amplitude PPU.

According to [40], using a relay is advantageous when:

I Uss," I2 I /,ss,Spl(n) I2 I Kss," I2

r- 1 • / 1 hk1,k2 1 1 hk2,0 1 N . 1 hk1,0 1

framel: min(-^-,---) >-^

k1,k2 I I "k2,0 I n > I k1,0_

L , L ) > L

v +T,Jh2 v2+I4P0(n) v2+]TJ",0 l=1 l=1 l=1

I 7„ss,« |2 I 7 ss,SP2(«) |2 I 7„ss,« |2 (7)

1 hk2,k1 1 1 K,S ) 1 w 1 hk2,0 1

k2,k1 k1,0 frame 2: min(-f-,-°-) > L

'pl ,0

v2 Jp",k1 v2 JSP20( n) v2 J

l=1 l=1 l=1

in selective DF mode. It is considered that the link of source node^ relay node and link of relay node^ destination node are better than that of source node^ destination node. Otherwise, the relay keeps idle on subcarrier n in the relaying phase for x*k2(t) orx*(t). In this paper, we just consider the case in which the relays keep working on each subcarrier, i.e., the Equation (7) is always true.

3. Optimization Problem Formulation

In this section, we analyze the joint optimization of subcarrier-pair based resource allocation algorithm for OFDM-DF based on full CSI and partial CSI, respectively. We are interested in how each SU allocates its power properly across its own data and its relayed data so as to maximize the system transmitted data rate while maintaining reasonable fairness between SUs. The optimization problem is formulated firstly and then solved in the dual domain. It is assumed that the PUs have a constant-rate, constant-power transmission, while the SUs are capable to adjust transmit power over different fading states based on the CSI of the CR network. We study a type of constraint imposed over the secondary transmission to protect the PUs by limiting the interference introduced to the PUs below a threshold.

3.1. Resource Allocation and Subcarrier Pairing Scheme Based on the OFDM-DF

The CR AP combines the received signals from the source node in T1 and the relay node in T2 through the maximal ratio combining. The transmit power is adjusted in each SU's transmitter. According to [27] and [41], when the link of source node->relay node transmission is successful for entire DF process, the transmission rate of kith SU and k2th SU at n subcarrier in relaying mode, which is connected via the Shannon capacity formula, can be shown as I"1SP1(n) and /"2SP2("), respectively:

.f pn,SP1(n) pn,SP1(n) pn,SP1(n)

ikniSP1(n) log2{i+min(|hk- |2 —^-,|hZ|2 —^-+|h;2;S0P1(n)I2 —k2L-)}

4 +1 JnPl ,k 2 +1 JnPl ,o +1J:^P0(n)

i=i i=i i=i

.„ pn,SP2(n) pn,SP2(n) pn,SP2(n)

/;f2(n) = + mindh- |2 -,|h-12 —^-+1 h;;,0SP2(n) |2 —^-)}

4 +1 Km +Z J"Pl ,0 +Z JSPT

l=1 l=1 l=1

where o2 denotes the Additive White Gaussian Noise(AWGN) variance. Here, it is assumed that all the channel gains are constant during two frames and the link between cooperative partners are symmetric, i.e., hsks{l2 = hsks2"j for all k. The factor 1/4 in Equation (8) results from the fact that the transmission takes four slots in the cooperative scheme.

Let P^f?1 ^ + P^Y1 ^ = P^SP 1 ^ . i = 1. 2 where i denotes the ith frame. This formula means

that the average of the transmit power of the source node Pk"1'SPi(n) and that of the relay node P^J1 is

constrained to be pj,SPi("), which is the allocated power on subcarrier n at the source node for direct

transmission. According to [42], the solution to this problem is the transmitted data rate and it is maximized when:

p",SP1(") P",SP1(") pn,SP1(n)

I lss," |2 Pk1,1_ _ I in," i2 Pk1,1__I J ss,SP1(n) |2 Pk 2,1_

| hk1,k2 | L =| hk1,0 | L + | hk2,0 | L

2 i ^ T" ^2 . ^ T" ^2 . ^ T"SP1(")

^ +Z J"Plkk2 ^ +Z JP,0 +Z J

pl,0 ^ 1 pl,0 l=1 l =1 l =1

,SP2( ) ,SP2( ) ,SP2( )

I iss,n |2 Pk2,2_ _ I iss," |2 Pk2,2__i 7 ss,SP2(n) |2 Pk1,2__(9)

| hk2,k1 | L =| hk2,0 | L + | hk1,0 | L

^ +]LJ"l,k1 ^ + 1^,0 ^ +]LJSP20(")

l=1 l=1 l=1

for all 1 < n < < k < ^

I Uss" I2 I Uss" I2 I Uss,n I2 I Uss,n I2

n | hk1k2 | n | "k1,0 | n | "k2,0 | n | kk2,k1 |

Yk1 ,k 2 L -,Гк1 ,0 L -,Гк 2 ,0 L -,Vk 2 ,k1 L

L k1,0 L k2,0 L k2,k1 L

+jk, ^+£.ъ* v2+z j;, ^,k1 (10)

l=1 l=1 l=1 l=1

Together with P^^ + P£fin) = P"fi(n\ i = 1,2we can obtain that:

VSP1(") v" - v"

p",SP1(") __vk2,0_ p",SP1(") p",SP1(") _ v k1,k2 v k1,0 p",SP1(")

Pk1,1 = v" _v" , VSP1(") Pk1 , Pk2,1 = v" _v" , VSP1(") Pk1

v k1,k 2 v k1,0 + v k 2,0 v k1,k 2 v k1,0 +vk 2,0 V" - V" VSP2(n)

p",SP2(") __/ k2,k1 / k2,0 p",SP2(") p",SP2(") _ _/k1,0_ p",SP2(")

Pk1,2 = " " , SP2(") Pk2 , Pk2,2 = " _ " SP2(") Pk2

Vk 2,k1 vk 2,0 + vk1,0 vk 2,k1 vk 2,0 + vk1,0

Denote %iSP1(n),%2SP2(n) as the equivalent channel gain given by:

vn SP1(n) n SP2( n)

n,SP1(n) __' k1,k2/ k2,0_ „n,SP2(n) _ ' k2,k1' ¿1,0_

™ ~ vn _vn , vSP1(n) 2 '¿2 - n n SP2(n) (12)

/k1,k2 / k1,0 + l k2,0 /k2,k1 /k2,0 + l k1,0

By now, we can unify the transmitted data rate as:

I^)(P^)) = f l0g2{1 Pk?SP1(n)>

/knf2(n)(PkfP2(n)) = ^ l0g2{1 + ^knf2(n)PkfP2(n)} (l3)

The proportional fairness is used as the optimized objective to develop subcarrier-pair based resource allocation algorithm in order to maintain a balance between system efficiency and fairness. According to [43,44] and the theorem in [45], we can readily deduce that there exists one unique proportionally fair allocation which can be attained by maximizing the objective function

Y = ££®kphl(in>SP1(n) + i«,sP2(n)) over the feasible set, where cok is the weighting factor to make the K

k=1 n=1

cooperative partners achieve the desirable transmitted data rate. Besides, we should keep the instantaneous interference introduced to the PUs below a certain threshold. The constraints include the aspects of satisfying the maximum power and interference constraints as well as the minimum rate requirements. Therefore, the resource allocation problem can be formulated mathematically as given in Equation (14). Constraint C1 corresponds to the subcarrier allocation constraint that each subcarrier n only can be allocated to one cooperative partner. C2 and C3 define that the sum of all the transmission powers of a particular SU on different subcarriers can't be greater than the maximum allowed limit for that particular SU. C4 ensures the cumulative interference from all SUs and through all subcarriers on a particular PU should not be greater than the interference limit set. C5 ensure that each SU can obtain the minimum rate requirements. This constraint precludes the possibility of multiple SUs simultaneously transmitting at the same subcarrier:

max Y = £ £ ^ Pn (Ikn1SP1( n) + Iknf2(n))

k=1 n=1 K

subject to C1: £pkn < 1,pkn e{0,1}, Vn

N p N p

C2: £PnPk"r(n) < p,C3: £< p

n=1 2 n=1 2 (14)

C4: yy p [pn,SP1( n) sn + pn,SP2(n) sn + pn,SP2( n) s SP2(n) + pn,SP1( n) s SP1(n)] < ^l)

^^^^ ^^^^ kn L k 1,1 k1, pt k 2,2 k 2, pt k1,2 k1, pt k 2,1 k 2, pt J

k=1 n=1

C5: £pJI"kr(n) + Iknf2(n)} *Rk

for all l = {1,2,..., L},k = {1,2,..., K}

where K denotes the number of cooperative partner, k denotes the kth cooperative partner, L denotes the number of PUs, l denotes the 1th PU, pkn is the allocation indicator that equals 1 when the nth subcarrier is allocated to the kth cooperative partner and 0 otherwise, Pt is the transmit power budget

for each cooperative partner, the factor of 1/2 for the terms of Pt/2 which results from the fact that it is a normalization for the transmissions within the duration of a frame, Ith(T) denotes the maximum allowable interference level at the lth PU receiver, Rk is the minimum transmitted data rate for kth cooperative partners.

The optimal solution to Equation (14) can be found by performing an exhaustive search with computational complexity O(KNZ) [46], where KN is the number of possible subcarrier allocations and Z is the complexity of a power allocation algorithm for each subcarrier allocation. To reduce the exponential computational complexity, a suboptimum resource allocation algorithm with less computational complexity is developed in the following. The dual decomposition approach is used to solve the problem. The dual problem of Equation (14) can be formulated as:

k n n k k p n

mm max ^ = £ £ ^ (4*1™ + ^2™) + I ^ "I Pn) + I "I prf

X(0),X(1), P,p k=1 n=1 n=1 k=1 k=1 2 n=1

X(2) ^(3) ^(4)

k u n l k n

|-pn,SP1(n)ç> n , ¡r,«,SP2(«)c n

+!42){p-Ip^™})-EEPkn[Pk1f1(n) Skni, Pl + Pk?2P2(n) Sk-2,

k=1 2 n=1 l=1 k=1 n=1

k1,S2P2(n) SkSPpn + Pkn2S1P1(n) *©?]}+1 w p (4fP1(n) + Ikn2SP2(n)) - Rk ]

k=1 n=1

subjectto Fnfi (n) > 0, X = {if, if, if, if, if} > 0 n e {1,2,...,N},{i, j} e {1,2},k e{1,2...,K}

where i and j denote the ith and jth SU of the kth cooperative partner, respectively. The values of {iflflfl^lf} are the introduced Lagrange multipliers. In the future using, we denote a vector shown in Equation (16):

1(0) = [^X-■40) f X» = [ , ■■ • ■•, xf f, 1(2) = [ > ■•• ■•> ■f

The Equation (15) can be decomposed into two layers of sub-problems. In the lower layer, we can get K sub-problems:

Ok(X(0), X(1), X(2), X(3), X(4)) = max £ZfP1(n),SP2(n)p

subjectto p e{0,1},Vng{1,2,...,N},Vke{1,2,...,K}

n-fr, \ p_/pn,SP1(n) pn,SP2(n) pn,SP1( n) pn,SP2( n) ï > n p = {pn }, P = {Pk1,1 , Pk1,2 , Pk2,1 , Pk2,2 } - 0

where (17)

rn,SP1(n),SP2(n) / pn,SP1(n) pn,SP2(n)\_ /j-n,SP1(n) , j-n,SP2(n)\_^l(3)r pn,SP1(n) on , pn,SP2(n) on

Lk (Pk1 , Pk 2 ) = œk (lk1 + Ik 2 ) ^^Z [Pk1,1 Sk1, pl + Pk 2,2 Sk 2, pl

pn,SP2(n) oSP2(n) pn,SP1(n) r»SP1(n) n n (0) ;i(1)pn,SP1(n) 0(2) pn,SP2(n) 3 (4V /-n,SP1(n) rn,SP2(n) \ +Pk1,2 Sk1, pl + Pk 2,1 Sk 2, pl ] Ak Pk1 Ak Pk 2 +Ak (lk1 + Ik 2 )

We suppose that Yk is the maximum value of the objective function in the lower layer. The master problem in the upper layer could be expressed as:

K N i K i K L K

HP) = !гк +£^ + + i^P +£

Д,(3),,(4) к-1 n-1 2 к-1 2 к-1 1-1 к-1 (18)

subject to {^f^UfU^}>о

Since a dual function is always optimized by first optimizing some variables and then optimizing the remaining ones. We define a subcarrier pairing parameter ¡n m e {0,1} that takes 1 if the nth subcarrier in T1 is pairing to mth subcarrier in T2 and 0 otherwise. We first optimize the primal variables with the assumption that dual variables {â(°),âJ(1),â£2),â1(3),âJ(4)} are given. The resource allocation and subcarrier pairing process can be divided into four stages:

(a) Allocating the optimal power factor Pf1^, Рк{^(п\ РгГ^, Р;™} for SUs. Pi'iP1(") and P"2s2P2(") imply the power used for self-data transmission, respectively. Pti2P2(") and P^fi4") imply the power used for partner-data transmission, respectively.

4. Allocating the optimal set of subcarriers 0.к for kth cooperation partner, i.e., obtaining the optimal subcarrier allocation factor p^.

5. Optimal pairing process for the subcarriers which are allocated to 0.к, i.e., allocating the optimal subcarrier pairing factor fînm.

6. After the temporary optimal primal variables have been obtained in each iteration process, we would find the temporary optimal dual variables ur^^u^,^}, which can minimize the objective function H(P)as shown in

Equation (18). 3.1.1. Power Allocation Algorithm

Let R^x^wp-^w,p;2SP2(n)) = /^w + /^O), if we make pn = 1, the power allocation can be determined in a water-filling fashion. Taking derivatives of Ll,SP1(n),SP2(n) with respect to

pn,SP1(n) pn,SP2(n) Рк 1 , Рк 2 .

•3i-n,SP1(n),SP2(n)/ pn,SP1(n) pn,SP2(n)\ ^pn,SP1(n),SP2(n)

___(Рк1 , Рк2 ) _(„ , 1(4)4 ___ i« V ^ГЛn,SP1(n^n , Fi,SP1(n^SP1(n)-| _ n

-PDn,SP1(n)- ~(0}к + Лк ) к 7)n,SP1(n)--Лк ~^Л1 Ик1,1 + Рк2,1 Sк2,pi J - 0

°Рк1 °Рк1 i-1

' pn,SP1(n) pn,SP2(nK -)pn,SP1(n),SP2(n) L

к_( Рк1 , Рк 2 ) _( , ;(4К °Лк__j(2) 3(3)ГП n,SP2(n^SP2( n) ^ Rn,SP2(n) о n

>n,SP2(n) -(Шк + Лк ) ^n,S

к2 к 2

pj-n,SP1(n),SP2(n) / pn,SP1(n) pn,SP2(n)\ ppn,SP1(n),SP2(n) L

°'Ьк_( Рк1 , Рк 2 ) , l(4K °Кк_ 1(2) V j(3)m n,SP2(n) ç<SP2( n). nn,SP2(n) n П _ n

-ôp„,SP2(n)- -(®к ) ôp„,SP2(n)--Д LDH,2 Sn,pi + Вк2,2 <\2,piJ - 0

vSP1(n) vSP2(n)

^n,SP1(n) __Ук2,0_ nn,SP2(n) __Ук1,0_

,1 _ Vn — v^-l- vSP1(n) , к2,2 _ n _ n SP2(n)

Ук1,к 2 Гк1,0 + /к 2,0 Ук 2,к1 /к 2,0 + Ук1,0

,n,SP2(n) __У к 2,к1 Ук 2 ,0 Tpn, SP1(n) __Ук1,к2 Ук1,0

SP2(n) , к2,1 _ „ „ ,,,SP1(n)

r\n,SP2(n) __/ к2,к1 / к2,0_ 7^n,SP1(n) _

,2 _ 3 _ n SP2(n) , к2,1 _ n _ „ sp

У к 2,к1 Ук 2,0 + Ук1,0 Ук1,к 2 Ук1,0 + Ук 2,0

Taking derivatives of R,SP1(l),SP2(l) = + pf^^with respect to pfP1(l), p;2SP2(i) , and according

to Equation (19), we can get:

SRn,SP1(n),SP2( n) gpn,SP1(n)

n,SP1(n) '!k 1

n,SP1(n) k ,1

4ln2 1 + 'fP1(n) pfP1(n)

^n,SP1(n) n(1) . V1 5(3W A n,SP1(n)on . 77n,SP1(n)pSP1(n)-| ©k ,1 = Ak + ^Al [Ak1,1 Sk1,pl + Fk 2,1 Sk 2, pl J

ÖRn,SP1(n),SP2( n)

n,SP2( n) ''k 2

n,SP2(n) k,2

n,SP2(n)

n,SP2(n) n,SP2(n)

4ln2 1 + 'SP2( n) P

©n,2P2( n) = +Z^i(3)[^;iS2P2( n) SSn) + B2SP2( n) Sfpi J

71 SPI(yi) tl SP2(n)

Together with the constraint P1' , P2' > 0, the temporary optimal solution can be obtained:

n,SP1(n) _

pn2SP2( n) = max

f ',SP1(n)H+44)) -©n,iSP1(n

4ln2_,

fi\n,SP1( n) n,SP1( n) ©k ,1 'lk1

A/ n,SP2( n)( i(4K /^n,SP2(

2 (®k + Ak ) ~©k ,2

4 ln 2 k2 k k k,2

n,SP2(n) n,SP2(n)

©k ,2 'Ik 2

The temporary SU's transmit power can be obtained for the given dual variables:

n,SP1(n) _ >n,SP1(n) pn,SP1(n)

pn,SP1(n Pk1,1

_ ^n,SP1(n) pn

= Ak1,1 Pk1

pn,SP1(n) _ pn,SP1(n) pn,SP1(n)

P2,1 _ Fk2,1 Pk1

n,SP2(n) k1,2

n,SP2(n) n,SP2(n)

n,SP2(n) k2,2

n,SP2(n) n,SP2(n) P

3.1.2. Subcarrier Allocation Algorithm

The subcarrier allocation constraint is that each subcarrier is allocated to no more than one SU cooperative partner, which prevents mutual interference among SUs. According to Section 3.1.1, we can get a temporary optimum power {^f1^,PtT^,PiJ2^,PS*00} . We substitute this

temporary optimum power vector into the objective functionLfP1(n),SP2(n)(pfP1(n),p«2SP2(n)) and the objective functionR^w^Mp^w,p«2SP2(n)) to obtain the temporary max valueLfP1(n),SP2(n) and R«,SP1(n),SP2(n), respectively. Taking account of resource fairness, we can formulate the optimization problem of subcarrier allocation as:

k * t p«Rri(nmf

max h (p) = £ £ PnLT^ + -

p ^ k£ (£ PnRr1^n))2

k=1 n=1 (23)

subject to Pkn ej0,1}, £p« <1 Vk,«,

where p = {p},1 < k < K,1 < n < N

where the value Z is the weighting factor to balance the total transmitted data rate and fairness index among SUs. The bigger the value Z is, the greater the fairness can be obtained, and otherwise the greater the transmitted data rate is. It is a mixed binary integer programming problem that is difficult to

solve. To reduce the exponential computational complexity, a suboptimum subcarrier allocation algorithm with less computational complexity is developed. The pseudo-code of subcarrier allocation algorithm can be described as follows:

(1) Initialization

(i) Make pM = 0, = 0, Vk, n

(2) For n = 1 to N

(i) For k = 1 to K

a) Make pkn = 1

b) Assign Vnk = h(p) according to Equation (23)

c) Make pkn = 0.

(ii) Assign k * = argmax Vnk

(iii) Assign pk,n =1 nk. = nk. Pn = 0, Vk * k *

According to Section 3.1.2, we can obtain a temporary subcarrier allocation vector P = {p,1 - k < K,1 < n < N} for the given dual variables.

3.1.3. Subcarrier Pairing Algorithm

The pairing constraint is that each subcarrier m in listening phase only pairs with at most one subcarrier n in the relaying phase. We assume that the pairing for deferent frames is not the same. The pairing process of the subcarrier allocated to kth cooperation partner can be expressed as:

Frame 1 : max £ £ Pnm I? ; Frame2:max £ £ p^ itf

neQk meQk neQk m'eQk

subject to C6 : £ p,m < 1,C7: £ P,m < 1,Vn e Qk

meQ^ m gQ^

C8 : £ P,m - 1,C9: £ p,m < 1,Vm e Qk ,Vm ' e Qk

neQ neQ^

C10:p, m e {0,1}, C11: p, „e {0,1}

Constraints C6, C8 and C10 correspond to the pairing constraint that each subcarrier n in listening phase only pairs with one subcarrier m in the relaying phase in the first frame. Constraint C7, C9 and C11 correspond to the pairing constraint that each subcarrier n in listening phase only pairs with one subcarrier m' in the relaying phase for the second frame. We can obtain the temporary optimal m, m' for any n as:

m = argmax {I£[m }, m = argmax {I^™ } /95^

mm ' ^ '

That is:

Frame 1 : pnmm = 1, SP1(n) = m, Pn,m = 0 Vm * m

Frame2 : Pnm, = 1, SP2(n) = m, Pnm, = 0 Vm ' * m (26)

The subcarrier pairing scheme can be shown as following: For k = 1 to K

(a) A = Q ,, A' = Q ,, p, = 0, p, = 0, Vn, m, m' eQt

(b) While A*0, A'*0

(i) n gA, n ' eA '

(ii) Find m gA, m' gA ' satisfying Inklm > Inkl, pf ' > I;-7 ', VI gA, Vl ' e A '

(iii) Assign

SP1(n) = m, Pn m = 1, A = A - {n} - {m} SP2(n ') = m', pn' m- = 1, A ' = A' - {n '} - {m '}

Through Section 3.1.3, we can obtain a temporary optimal subcarrier pairing vector /?n,m, (3n,m>, 1 < {n, m, m'} < N for the given dual variables.

3.1.4. Optimizing the Dual Variables

The optimal values of dual variables can be achieved iteratively by the sub-gradient method as follows:

^n°)J+1 = {¿n?)J-<[1 -I Pn ]}+ k=1

AT1 = {a™c p-IPk,Pkn1,SP1(n)]}

42>,i+1 = A™-<[P-IPk„Pkn2SP2(n)]}+ (27)

A^1 = W-c^ -IIPn(PktSP1(n)s;,p1 + Pkn2S2P2(n)skpl + P;^sspp« + Pkn2SP1(n)SkSp;^n))]}+

k=1 n=1

at1={Ak4X--c«[^ Pn (i:rn) + i;2sp2( n)) - Rk ]}

where{c((iJ),c\'i,c(2i,c()),c(lyis a small positive step size for the ith iteration, i is the iteration number

and( ^)+ = max( ^ ,0), with the appropriate step sizes, the iterations are convergent. The remaining issue is how to determine the step sizeC^, Cfk ,a(2')k,a(3'],a((!)k}. Clearly, performing a line search at each iteration process perform well. For a given current iteration H = {A(0)'i ,A2)'' ,A(3)'' ,A4),i} and a search direction d(i), we compute step size ak'(i) by:

(cS'Ci'cS-CSl-CS} = (0 afg^i^ (o0^

c0,n'a1,k 'a2,k ,a3,l ,a4,k

A* +cfM'), Ak2),' +c2M) ,AP' +c;(')L,Ak4),' -cgd«}

K N K (28)

K rk+I A;o)+I ^

k=1 n=1 k=1

where QU^A^AfUf), A^) = 2^ +

+I42)P+IAfWl -2AfR, Vk,n k=1 1=1 k=1

Within each iteration process, the power allocation vectors can be updated respectively by Equation (22), the subcarrier allocation vectors can be updated respectively by subcarrier allocation algorithm shown in Section 3.1.2, the subcarrier pairing vectors can be updated respectively by subcarrier pairing algorithm which is shown in Section 3.1.3, with the updated value {Af ,Ak2)'1 ,A(3),i}. Therefore, the dual variable 1 ' will converge to the dual optimum A as

i ^ œ and the temporary primal optimum variable will also converge to the primal optimum value after several iterations, i.e., temporary vector {P™, P;2S1P1("), PTj'^, PCf^} will converge to optimal

vector {P^1(n),pk1;2P2(n),P2r(n),P2S2P2(n)} , temporary optimal vector p = {p, 1 < k < K,1 < n < N}will

converge to the optimal vector p = {pfa ,1 < k < K,1 < n < N}, temporary vector /3nm, Pnm, will converge

to the optimal vector 3nm, 3n m, .

3.2. Resource Allocation with Partial CSI

If full CSI can be achieved at the k1th and k2th (1 < k < K) SU transmitter, the optimal subcarrier allocation vector p , subcarrier pairing vector /?n,m, /?n,m, , the power allocation vector

{Pk«1Spi(n),pnf2(n),Pkn2SP1(n),Pkn2,S2P2(n)} can be determined simply by the proposed algorithm which is

shown in Section 3.1. However, the practical case in which only partial CSI of the wireless channel between the secondary base station and SUs is available have to be considered. The CSI of non-adjacent link may be undesirable and even unavailable when the SUs are mobile. In this section, we investigate the optimal resource allocation in SU cooperation network with partial CSI at each transmitter. Specifically, the kith SU transmitter has full CSI of its adjacent links h^f^(i ^ j), ti^Q

but only statistical CSI of non-adjacent linkhp^Q . We assume that the SUs can knowhpsn,hpp]'l, i.e., hpi m,hpi 'o are still available at SU transmitter, and the link between cooperative users is symmetric for simplicity, i.e., hpSk; = h*]! .

Under these assumptions and according to [47], the objective function:

*=EI ^k p; {ikn1SP1(n) (PkîSP1(n) )+i;2SP2(;) ( Pkn2,SP2(n))} (29)

k=1 n=1

can be rewritten as:

* = II^kP;{lnI"k^iP1(n)(P;SP1(n)) + ln I^^p;^21^)} (30)

where I;1SP1(;) and I;2SP2(;) are defined as:

Ik1 (P1 ) = Ehss,n [Ik1 (Pk1 )]

t;,SP2(;V pn,SP2(n^^_ p r t;,SP2(; V pn,SP2(n)\-i Ik 2 (P2 ) = Ehss,n [*k2 (Pk2 )]

In order to seek the optimal power allocation solution, we derive the explicit expressions for

I;,SP1(;) and Il,SP2(n) described as:

Leth^n = )vxti0,i = 1,2, wheredfa.Odenotes the distance between the kith SU and the AP, v is

the path-loss exponent, xki,0, (i = 1,2) is the normalized complex Gaussian random variable distributed as CN (0, 1), then at high signal-to-noise ratio (SNR):

- Af Pn,SP1(n)vn W + W Pn,SP1(n) vn W + W

jn,SP1(n) _ Af (rk1 fk1,k2' 1 + ' 1 n _ Y k1,k2W1 + W1 x _

^ , n n V n n '

}n,SP1(n) n

4ln2 Yk1,k 2 Yk1,0 Yk1,k 2 Yk1,0 Yk1,k 2 ^1,0 Yk1,k 2 Yk1,0

Af P.f^Y^W + W PT^Y^kW + W2

jn,SP2(n) _ L-y i k2 i k2,k1 2 k2 k2,k1 2 2 \

k 2 /i i n n ^ n n ^

;>n,SP2(n) n

4lnz Ykn2,k1 "Ykn2,(

Yk 2,k1 Yk 2,

Yk 2,k1 Yk 2,0 Yk 2,k1 Yk2:

where Wi, W2 in Equation (32) are defined as:

W = [(a2 +tjlPr)d2ki0T\W2 = [(a2 +]^Ddk21V0]-1

Proof: It is sufficient to shown, and Inklvl{r>n can be derived in exactly the same way, then

7-n,sp1(n) _f fafu0n +„",SP1(")pn,sp1(n)n

Ik1 = phss_,n [1k1 ] = Ehss_,n { A log2(1 + ^k1 Pk1 )}

= Ess," !og2(1 +

hk 2,0 4

hk 2,0 4

Yn ysp1(«) Y k1,k 2Yk 2,0

hk 2,0

v" -yn +ysp1(n) k1

Yk1,k 2 "¿1,0 + 2,0

n Ihgr^l2 Yk1,k 2

3«,sp1(n)

Af !og2(1 +

jsp1(n) 2 Jpl ,0

/=1 p^spkn)

, , ss,sp1(n) ,2 k1

vn vn + 1 h 2,0 I

7k1,k2 - Yk1,0 + l

jsp1(n)

2 J pi ,0

Yki1,k2[(a2 +1 Jnpl,0)d2v,0]-11 xk2,0 I2 pkr1(n)) log2 (1 +---;-~-

(Yk1,k2 - Ykt1,0 + "

(a2+1 j sTV^

The third equality holds since we assume each transmission block is long enough to undergo different channel realizations as argued in [48]. Next, let:

W = [(a2 +2 J^TX«,!"1,y =I Xk2,0 I2 then according to [49], we can obtain:

rn,sp1(n) _ E rrn,sp1(n)1 _fE I ,

1k 1 _ Ehss,n [1k 1 1 _ , Ehss,n \ log2

hk 2,0 4 hk 2,0

rhk2 - Yk1,0 + W1y

1, rkn1,kWyP^1" rkn1,k2 - rk1,0 +w1 y

hf f+- e"yPkn1,SP1(n)rkn1,k2W1(rkn1,k2 " rkn1,0)dy

4ln2 J0 (rkn1,k2 - rkn1,0 + W1 y)(rkn1,k2 - rkn1,0 + W1 y + Pkn1,SP1(n)rkn1,k2W1 y)

hf r+TO

P{-------} dy (36)

4ln2'0 r,n1,k2 -rk1,0 r,n1,k2 - rk:1,0

y W y Pn,SP1(n)vn w+w

W1 Pk1 rk1,k 2W1 + W1

Yk\,k2-ïk\,0 r" n-v".

hf f1^^dV + e^VnV2W1+W1 fPkn1,SP1(n)rn1,k2W1+W1 y

J-TO 1; J-TO 1;

41n2 J-TO y J-TO y

Yk\,k 2-rk1,0 n n _rk1,k 2-rk1,0__n n

-hf W r~,rn1,0 -rn1,k2. P/fP1(nV"1k2W1+W1rv rm,0 -rn1,k2 „ i-e 1 Ei( , _ , )+e k1 rk1,k2 1 1 Ei( spi ' ' )}

41nr ^ ^ ^"VA,^ +^1'

where A/'(I) denotes the exponential integral function defined as Ei(t) = exx dx and Ii,(t) can be

J—CO

expanded asymptotically as:

e n' k e* 11

E'(t) - ^ (1 + £ ^f1-) - ^ [1 +1 + o(1)1 (37)

t x t t t

Then using the assumption that SNR is high, we readily obtain Equation (32). This completes the proof.

We can use the resource allocation and subcarrier pairing algorithm which is proposed in Section 3.1 to solve the corresponding optimization Equation (14) once again. When partial CSI can be achieved by SUs, the optimal subcarrier allocation vector p, subcarrier pairing vector (3nm, <

and opti can get

and optimal power allocation vector P^SP 1 ^, plSP2n can be obtained. According to Equation (11), we

pn,SPl(ri) _ E ,__ p

Ukl,k2 Yk 1,0 Yic2,0 )

- f+œ_Î^Z_'pnJpiCn) -y ,

Jo rhv-rZi.° + y n

'pn,SPl(nj f L__Vkl,k2 ~ Yk 1,0 I -y ,

kl Jo I rki,k2-rk o + Me dy

SPl(n) A ~Yk2,0_ pn,SPl(n) (

Yki, o~ Yki, k 2 (38)

via ec

= ^IR+ Yko) J .

~Hkl + w± 61 { W± /

_ 'pn,SPl(rij

-l,n _ 17n fcl

Kfcl,fc2 Kfcl.O

Then we also can get:

pn,5Pl(n) _ / 1__\ pn,5Pl(n) pn,SP2(ri)

k.2,1 — I 1 ...n _ ...n I "fcl ' kl,2 \ Vkl,k2 Yk 1,0/

(1__\ p7i,SP2(n) pn,SP2(n) _ pn,SP2(nj

_ k2 k2,2 ~ 17n _ 17n fc2

K k2,0/ Vk2,kl Vk2,0

Intuitively, the system performance would be degraded due to limited CSI though this scheme which does not require each cooperative node to have full CSI of nonadjacent link, as verified in the numerical simulation.

3.3. Comparison with Classical Resource Allocation Algorithms

Several existing schemes are compared with the proposed subcarrier-pair based resource allocation algorithm in terms of the system transmitted rate and fairness, respectively. These existing schemes include the following:

LWF-PI-without-SP: The subcarrier is allocated according to the channel gain. The messages transmitted on subcarrier m at the source node will be retransmitted on the same subcarrier m at the relay node. The power is allocated according to LWF-PI algorithm [23] on each subcarrier.

EPA-without-SP: The subcarrier is allocated according to the channel gain. The messages transmitted on subcarrier m at the source node will be retransmitted on the same subcarrier m at the relay node. The power is allocated equally on each subcarrier.

LWF-PI-with-SP: The subcarrier is allocated according to the channel gain. The messages transmitted on subcarrier m at the source node will be retransmitted on subcarrier n, which is selected by proposed subcarrier pairing algorithm, at the relay node. The power is allocated according to LWF-PI algorithm [23] on each subcarrier.

EPA-with-SP: The subcarrier is allocated according to the channel gain. The messages transmitted on subcarrier m at the source node will be retransmitted on subcarrier n, which is selected by proposed subcarrier pairing algorithm at the relay node. The power is allocated equally on each subcarrier.

Optimal-Scheme-with-SP: The subcarrier is allocated according to the channel gain. The messages transmitted on subcarrier m at the source node will be retransmitted on subcarrier n, which is selected by proposed subcarrier pairing algorithm at the relay node. The power is allocated according to Optimal Scheme [37] on each subcarrier.

Partial CSI: the difference between the previous five algorithms and Partial CSI scheme is that the SU has only imperfect CSI of non-adjacent link under partial CSI scheme. In this case, the objective function is formulated by statistical methods.

4. Simulation Results

We have studied asymmetric or linear network with all SUs of kth cooperative partner located in the same line. Specifically, kith SU and the destination AP are fixed at (0, 0), (1, 0) respectively, and k2th SU is located at (d, 0), 0 < d < 1, without loss of generality. The results for path loss exponent v = 2 are presented and all channels are modeled as Rayleigh flat fading with AWGN. Some simulation parameters are shown in Table 3.

Table 3. Simulation Parameters.

Parameter Value

Number of cooperative partners K 4

Number of PUs L 2

Number of subcarriers N 20

Length of a slot Ts 4u s

value of amplitude PPU 10 x 10-3 W

Ith(l) (l = 1,2) 2.7 x 10-3 W

average channel power gain 10 dB

△f 0.315 MHz

B2 1 MHz

B2 2 MHz

Simulation results are presented in this section to verify the performance of the proposed subcarrier-pair based resource allocation algorithm. In our simulations, the CVX, a package for specifying and solving convex programs, is used to solve formulated optimization resource allocation problems.

The channel gainsh^,h&o,hppi,Kit,K^n used in this paper are assumed to be Rayleigh fading,

since the channel fading gains for different realizations of channel gain can be different, an average transmission capacity of 10,000 independent simulation runs is considered. And individual fairness index is defined as [46]:

{E E PkR",SP1( n),SP2(n)}2

fair = —lT^N--(40)

K E (E PknRk",SP1(n),SP2( n) )2

k=1 n=1 -

The fairness ranges between 0 and 1. The higher the value fairness is, the more fair the throughput distribution among SUs is.

4.1. The System Transmitted Data Rate Obtained under Resource Allocation Algorithms

As shown in Figure 4, the achievable maximum CR system transmitted data rate is plotted versus the power budget Pt. The upper curve denotes the transmitted data rate by proposed subcarrier-pair based resource allocation algorithm with full CSI. It can be noted that the proposed subcarrier-pair based resource allocation algorithm achieves the highest transmitted data rate under a given total power constraint. The transmitted data rate achieved using the proposed algorithm is the highest among that using LWF-PI-with-SP, Optimal-Scheme-with-SP algorithm, LWF-PI-without-SP algorithm and EPA algorithm. The main reason is that the proposed algorithms can make full use of the entire available interference threshold, while the LWF-PI algorithms can only guarantee that the total interference is under the interference threshold as shown in Figure 5. We observe that only the proposed algorithm outperforms the compared existing resource allocation algorithms. Therefore, we can conclude that the proposed subcarrier-pair based resource allocation algorithm makes valuable contribution to system transmitted data rate.

Figure 4. Maximum transmitted data rate versus Power Budget.

4.2. The Fairness Index Obtained under Resource Allocation Algorithms

Both the proposed resource allocation algorithm and EPA-with-SP algorithm exhibit the best fairness performance and LWF-PI-without-SP algorithm shows the least fairness as shown in Figure 6. Moreover, the fairness loss of EPA-with-SP allocation algorithm compared to the proposed subcarrier-pair based resource allocation algorithm is acceptable. In the LWF-PI-without-SP algorithm, LWF-PI-with-SP algorithm and Optimal-Scheme-with-SP algorithm, most power and subcarrier will be assigned to SUs with good channel conditions to improve system efficiency. Unlike the proposed subcarrier-pair based resource allocation algorithm and EPA schemes, these two schemes achieve significant performance improvement to ensure fairness among SUs. The fairness obtained by EPA-with-SP effectively approaches to 0.98. Due to multiuser diversity, the fairness attained by the proposed algorithm is above 0.95, while the achievable system transmitted data rate is higher than that of LWF-PI and EPA-without-SP algorithms.

Figure 6. Power Budget versus Fairness index.

4.3. The Transmitted Data Rate of Each SU for the Resource Allocation Schemes

As shown in Figure 7, the achievable transmitted data rate of each SU is plotted. The transmitted data rate of each SU, under the proposed subcarrier-pair based resource allocation algorithm, is the most balanced and stable among EPA, LWF-PI and Optimal Scheme [37]. Especially under the LWF-PI algorithm and Optimal Scheme [37], some SUs with bad channel conditions have lower transmitted data rates, which would result in unsuccessful communication with destination node. We can conclude that the proposed subcarrier-pair based resource allocation algorithm makes valuable contribution to balance transmitted data rates among SUs.

Figure 7. Maximum transmitted data rate versus each SU (The power budget Pt is 20 x 10-3 W).

4.4. The System Transmitted Data Rate Obtained under Full CSI and Partial CSI

As shown in Figure 8, the performance of proposed subcarrier-pair based resource allocation algorithm is presented for the case that each SU has no CSI of non-adjacent link. It is shown that the

system transmitted data rate decreases due to limited CSI compared to that of full CSI scenario. However, this scheme does not require the full CSI of non-adjacent link at each cooperative SU. Moreover, we show that the transmitted data rate of the proposed subcarrier-pair based resource allocation algorithm with partial CSI significantly outperforms that of EPA. The transmitted data rate is close to that of subcarrier-pair based resource allocation algorithm with full CSI.

Figure 8. Power Budget versus sum transmitted data rate.

4.5. The Transmitted Data Rate Obtained by SUs under Different Distance

Last but not least, we have also studied how the distance between cooperative SUs impacts on the system performance. Under the proposed subcarrier-pair based resource allocation algorithm with proportional fairness, Figure 9 shows that the system has comparatively better performance when d e [0.4,0.6]. This result provides a guideline for grouping and partner selection in user cooperative networks.

Figure 9. Maximum transmitted data rate versus distance between two cooperative SUs.

5. Conclusions

In this paper, we have developed a novel subcarrier-pair based resource allocation algorithm that maximizes the transmission data rate while the interference introduced to the PUs remains within a given limit. Using the proportional fairness as the optimized objective function, we can improve proportional fairness of resource allocation and achieve substantial transmitted data rate gains. The sum of power constraint for the source and relay nodes are considered. Moreover, we extended the analysis to the case that the CSI of nonadjacent link is not available at cooperative SUs, and found that even in this case the proposed schemes perform better than the classical schemes.

Compared to the existing resource allocation algorithms which are introduced in the paper, our algorithm considers that each secondary relay user has its own data to be transmitted. Simulation results have shown that, either in improving the system throughput or in improving fairness of resource allocation, the proposed subcarrier-pair based resource allocation algorithm offers the best performance conditions among several existing compared resource allocation algorithms under various power budgets, while keeping the interference introduced to PU bands below a specified threshold. Besides, the transmitted data rate of proposed subcarrier-pair based resource allocation algorithm obviously outperforms that of EPA scheme when only partial CSI can be obtained by SUs. The contribution by taking full advantage of the statistical non-adjacent links channel information is demonstrated clearly in the simulation results. In addition, we notice that the system efficiency loss of partial CSI scheme compared to that of full CSI scheme is acceptable.

Acknowledgments

This research was supported by the key project of Tianjin Key Technology R&D Program (09ZCKFGX29200) of China.

Conflict of Interest

The authors declare no conflict of interest. References

1. Simon, H. Cognitive radio brain-empowered wireless communications. IEEE J. Sel. Areas Commun. 2005, 23, 201-220.

2. Liang, Y.C.; Chen, K.C.; Li, G.Y. Cognitive radio networking and communications: An overview. IEEE Trans. Veh. Technol. 2011, 60, 3386-3407.

3. Laneman, J.N.; Tse, D.N.C.; Wornell, G.W. Cooperative diversity in wireless networks: Efficient protocols and outage behavior. IEEE Trans. Inf. Theory 2004, 50, 3062-3080.

4. Choi, Y.J.; Pak, W.; Xin, Y.; Rangarajan, S. Throughput analysis of cooperative spectrum sensing in Rayleigh-faded cognitive radio systems. IET Commun. 2012, 6, 1104-1110.

5. Kang, X.; Liang, Y.-C.; Nallanathan, A. Optimal power allocation for fading channels in cognitive radio networks: Ergodic capacity and outage capacity. IEEE Trans. Wirel. Commun. 2009, 8, 940-950.

6. Kang, X.; Zhang, R.; Liang, Y.-C. Optimal power allocation strategies for fading cognitive radio channels with primary users outage constraint. IEEE J. Sel. Areas Commun. 2011, 29, 374-383.

7. Liu, J.; Shroff, N.B.; Sherali, H.D. Optimal power allocation in multi-relay MIMO cooperative networks: Theory and algorithms. IEEE J. Sel. Areas Commun. 2012, 30, 331-340.

8. Kwangyul, K.; Yoan, S. An Improved Power Allocation Scheme Using Particle Swarm Optimization in Cooperative Wireless Communication Systems. In Proceedings of the 17th Asia-Pacific Conference on Communications (APCC), Sabah, Malaysia, 2-5 October 2011; pp. 654-658.

9. Oberoi, J.S.; Pareek, U.; Naeem, M. EDA-Based Joint Power, Subcarrier Allocation and Relay Assignment Scheme for Multiuser Relaying in OFDMA-Based Cognitive Radio Systems. In Proceedings of the 2011 5th International Conference on Signal Processing and Communication Systems (ICSPCS), Honolulu, HI, USA, 12-14 December 2011; pp. 1-7.

10. Ding, Z.G.; Chin, W.H.; Leung, K.K. Distributed beamforming and power allocation for cooperative networks. IEEE Trans. Wirel. Commun. 2008, 7, 1817-1822.

11. Luo, J.; Blum, R.S.; Cimini, L.J. Decode-and-forward cooperative diversity with power allocation in wireless networks. IEEE Trans. Wirel. Commun. 2007, 6, 793-799.

12. Pan, Y.W.; Nix, A.; Beach, M. Distributed resource allocation for OFDMA-based relay networks. IEEE Trans. Veh. Technol. 2011, 60, 919-931.

13. Shen, Y.Y.; Feng, G.; Yang, B. Fair resource allocation and admission control in wireless multiuser amplify-and-forward relay networks. IEEE Trans. Veh. Technol. 2012, 61, 1383-1397.

14. Salem, M.; Adinoyi, A.; Rahman, M.; Yanikomeroglu, H.; Falconer, D.; Kim, Y.D. Fairness-aware radio resource management in downlink OFDMA cellular relay networks. IEEE Trans. Wirel. Commun. 2010, 9, 1628-1639.

15. Vu, H.N.; Kong, H.-Y. Joint subcarrier matching and power allocation in OFDM two-way relay systems. J. Commun. Netw. 2012, 14, 257-266.

16. Boostanimehr, H.; Bhargava, V.K. Selective subcarrier pairing and power allocation for DF OFDM relay systems with perfect and partial CSI. IEEE Trans. Wirel. Commun. 2011, 10, 4057-4067.

17. Dang, W.B.; Tao, M.X.; Mu, H. Subcarrier-pair based resource allocation for cooperative multi-relay OFDM systems. IEEE Trans. Wirel. Commun. 2010, 9, 1640-1649.

18. Sendonaris, A.; Erkip, E.; Aazhang, B. User cooperation diversity part I: System description. IEEE Trans. Commun. 2003, 51, 1927-1938.

19. Tourki, K.; Gesbert, D.; Deneire, L. Cooperative Diversity Using Per-User Power Control in the Multiuser MAC Channel. In Proceedings of the 2007 IEEE International Symposium on Information Theory, Nice, France, 24-29 June 2007; pp. 1911-1915.

20. Maciel, T.F.; Klein, A. On the performance, complexity, and fairness of suboptimal resource allocation for multiuser MIMO-OFDMA systems. IEEE Trans. Veh. Technol. 2010, 59, 406-419.

21. Papoutsis, V.D.; Fraimis, I.G.; Kotsopoulos, S.A. User selection and resource allocation algorithm with fairness in MISO-OFDMA. IEEE Commun. Lett. 2010, 14, 411-413.

22. Lu, W.S.; Ji, F.; Yu, H. A general resource allocation algorithm with fairness for SDMA/MISO/OFDMA systems. IEEE Commun. Lett. 2011, 15, 1072-1074.

23. Wu, J.; Yang, L.X.; Liu, X. Subcarrier and Power allocation in OFDM Based Cognitive Radio Systems. In Proceedings of the 2011 Fourth International Conference on Intelligent Computation Technology and Automation (ICICTA), Shenzhen, China, 28-29 March 2011; pp. 728-731.

24. Mohanram, C.; Bhashyam, S. A sub-optimal joint subcarrier and power allocation algorithm for multiuser OFDM. IEEE Commun. Lett. 2005, 9, 685-687.

25. Tan, C.K.; Chuah, T.C.; Tan, S.W. Fair subcarrier and power allocation for multiuser orthogonal frequency-division multiple access cognitive radio networks using a Colonel Blotto game. IET Commun. 2011, 5, 1607-1618.

26. Mohamad, K.A.; Veluppillai, M.; Mehri, M. A dual-decomposition-based resource allocation for OFDMA networks with imperfect CSI. IEEE Trans. Veh. Technol. 2010, 59, 2394-2403.

27. Chen, C.-H.; Wang, C.-L.; Chen, C.-T. A resource allocation scheme for cooperative multiuser OFDM-based cognitive radio systems. IEEE Trans. Commun. 2011, 59, 3204-3215.

28. Yaacoub, E.; Dawy, Z. A survey on uplink resource allocation in OFDMA wireless networks. IEEE Commun. Surv. Tutor. 2012, 14, 322-337.

29. Nader, M.; Keivan, N.; Mohammad, G.K. Downlink radio resource allocation in OFDMA spectrum sharing environment with partial channel state information. IEEE Trans. Wirel. Commun. 2011, 10, 3482-3495.

30. Choi, K.W.; Hossain, E.; Kim, D.I. Downlink subchannel and power allocation in multi-cell OFDMA cognitive radio networks. IEEE Trans. Wirel. Commun. 2011, 10, 2259-2271.

31. Tianxiang, L.; Feifei, G.; Xian-Da, Z. Joint resource scheduling for relay-assisted broadband cognitive radio networks. IEEE Trans. Wirel. Commun. 2012, 11, 3090-4000.

32. Sabit, E.; Mohamed, M.A.; Khalid, A.Q. Random subcarrier allocation in OFDM-based cognitive radio networks. IEEE Trans. Signal Process. 2012, 60, 4758-4774.

33. Xu, H.; Li, B.C. Resource allocation with flexible channel cooperation in cognitive radio networks. IEEE Trans. Mob. Comput. 2013, 12, 957-970.

34. Bansal, G.; Hossain, M.D.J.; Bhargava, V.K. Subcarrier and power allocation for OFDMA-based cognitive radio systems with joint overlay and underlay spectrum access mechanism. IEEE Trans. Veh. Technol. 2013, 62, 1111-1122.

35. Naeem, M.; Pareek, U.; Lee, D.C.; Anpalagan, A. Estimation of distribution algorithm for resource allocation in green cooperative cognitive radio sensor networks. Sensors 2013, 13, 4884-4905.

36. Viswanath, P.; Tse, D.N.C.; Laroia, R. Opportunistic beamforming using dumb antennas. IEEE Trans. Inf. Theory 2002, 48, 1277-1294.

37. Bansal, G.; Hossain, M.D.J.; Bhargava, V.K. Optimal and suboptimal power allocation schemes for OFDM-based cognitive radio systems. IEEE Trans. Wirel. Commun. 2008, 7, 4710-4718.

38. Sami, M.A.; Gordon, L.S. Interference-aware radio resource allocation in OFDMA-based cognitive radio networks. IEEE Trans. Veh. Technol. 2011, 60, 1699-1713.

39. Kang, X.; Garg, H.K.; Liang, Y.-C. Optimal power allocation for OFDM-based cognitive radio with new primary transmission protection criteria. IEEE Trans. Wirel. Commun. 2010, 9, 2066-2075.

40. Liu, Y.; Chen, W. Limited-feedback-based adaptive power allocation and subcarrier pairing for OFDM DF relay. IEEE Trans. Veh. Technol. 2012, 61, 2559-2571.

41. Sharma, S.; Shi, Y.; Liu, J. Network coding in cooperative communications: Friend or foe? IEEE Trans. Mob. Comput. 2012, 11, 1073-1085.

42. Zhang, Q.; Zhang, J.M.; Shao, C.J. Power Allocation for Regenerative Relay Channel with Rayleigh Fading. In Proceedings of IEEE 59th Vehicular Technology Conference (VTC 2004), Milan, Italy, 17-19 May 2004; pp. 1167-1171.

43. Le, L.B.; Hossain, E. Resource allocation for spectrum underlay in cognitive radio networks. IEEE Trans. Wirel. Commun. 2008, 7, 5306-5315.

44. Kelly, F.P.; Maulloo, A.K.; Tan, D.K.H. Rate control in communication networks: Shadow prices, proportional fairness and stability. J. Oper. Res. Soc. 1998, 49, 237-252.

45. Kim, H.; Keunyoung, K.; Youngnam, H. A Proportional Fair Scheduling for Multicarrier Transmission Systems. In Proceedings of IEEE 60th Vehicular Technology Conference (VTC2004-Fall), Los Angeles, CA, USA, 26-29 September 2004; pp. 409-413.

46. Shen, Z.; Andrews, J.G.; Evans, B.L. Adaptive resource allocation in multiuser OFDM systems with proportional rate constraints. IEEE Trans. Wirel. Commun. 2005, 4, 2726-2737.

47. Ma, Y.-T.; Lin, X.Q. Power Allocation with Proportional Fairness in User Cooperation Networks. In Proceedings of the 19th Annual Wireless and Optical Communications Conference (WOCC), Shanghai, China, 14-15 May 2010; pp. 1-5.

48. Host-Madsen, A.; Zhang, J.S. Capacity bounds and power allocation for wireless relay channels. IEEE Trans. Inf. Theory 2005, 51, 2020-2040.

49. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Academic Press: San Diego, CA, USA, 2007.

© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Copyright of Sensors (14248220) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.