CSI compared to the other newly proposed resource allocation schemes. Index TermsâOFDMA, relay selection, subcarrier alloca- tion,power allocation, imperfect CSI and cooperative commu- nications. ...... New York: Academic, 2007.
Resource Allocation for Cooperative Relay-assisted OFDMA Networks with Imperfect CSI ∗ Department
Zheng Chang∗† and Tapani Ristaniemi∗ of Mathematical Information Technology, University of Jyväskylä, P.O.Box 35, FIN-40014 Jyväskylä, Finland † Magister Solutions Ltd, Hannikaisenkatu 41, FIN-40100 Jyväskylä, Finland Email: {zheng.chang, tapani.ristaniemi}@jyu.fi
Abstract—This work addresses the radio resource allocation (RRA) problem for cooperative relay assisted OFDMA wireless network. The relays adopt the decode-and-forward protocol and can cooperatively assist the transmission from source to destination. The RRA scheme addresses practical implementation issues of resource allocation in OFDMA networks: the inaccuracy of channel-state information (CSI) available to the source. Instead, the source only knows estimated channel status and distributions of related estimation errors. The objective is to maximize the system throughput of the source-to-destination link under various constraints. Since the optimization problem is known as NPhard, we divide the original problem to three subproblems including relay selection, subcarrier and power allocations. We derive theoretical expressions for the solutions and illustrate them through simulations. Results validate clearly that our proposed algorithm can enhance the performance of system with imperfect CSI compared to the other newly proposed resource allocation schemes. Index Terms—OFDMA, relay selection, subcarrier allocation,power allocation, imperfect CSI and cooperative communications.
I. I NTRODUCTION Orthogonal Frequency Division Multiple Access (OFDMA) is an effective technique that exploits the benefits of Orthogonal Frequency Division Multiplexing (OFDM) for combating against channel noise and multipath effects and finally enables high data rate transmissions over fading channels. Meanwhile, cooperative communication has emerged as one of the main trends to reach even better system performance in terms of throughput, energy efficiency or cell coverage. Therefore, the incorporation of OFDMA and cooperative relays is foreseen to result in a promising structure that offers a possibility to reach many desirable objectives for the future wireless networks. However, a combination of a conventional one-to-many (single hop) OFDMA system and a relay network calls for a careful design of the radio resource allocation (RRA) principles. This means a carefully design and coordination of the power and subcarrier allocation, selection of relay(s) across different hops and optimizing the resource between the hops. The RRA algorithm plays an important role in the developments of both conventional and relay-aided OFDMA systems. The related works have been widely done in several different areas [1]-[6] when assuming perfect channel state information (CSI) is known to the source. A cross-layer optimization algorithm for resource allocation in conventional OFDMA network has been presented in [1] without considering relaying. An iter-
ative algorithm is proposed to solve the subcarrier assignment together with relay selection in [2]. Then, the power allocation problem can be solved by another iterative method based on waterfilling algorithm. Authors in [3] introduces closed-form solution for radio resource allocation for multihop cooperative relay network. However, the per-tone power constraint is used which is not practical. Scheme used in [4] considers fairness constraints when selecting relays. In [5], a threshold method is used to solve two subproblems, subcarrier allocation and power allocation. Although the performance is improved comparing to some other algorithms, the total power constraint is considered, which is not a realistic case since each node has its own power limitation. The work in [6] also proposed a subcarrier and relay pairing algorithm to solve the existing RRA problem, which requires high complexity. [7]-[9] present the work about RRA with imperfect CSI. [7] consider the RRA algorithm for conventional OFDMA networks. [8] investigates the issue of joint RRA and relay selection with imperfect CSI. However, authors use mean rate to characterize the CSI uncertainty which results in different interpretations. Recent work about RRA for OFDMA relay networks with imperfect CSI is introduced in [9], where only one relay is selected for assisting the transmission. As we can see, the existed algorithms have their major drawbacks which need to be improved. In this paper, we investigate a new resource allocation scheme for OFDMA cooperative network with imperfect CSI, which can effectively solve the problems of joint relay selection, subcarrier and power allocation and thus, enhance the system throughput when estimated error existed. In this work, relays are deployed for extending the cell coverage, so we do not consider direct link from source to destination. Conditional expected throughput is considered as our performance evaluation metric. We propose a relay selection and subcarrier allocation schemes, where one set of relays that can obtain the best link data rate is selected. Power is allocated to the source and relays under pernode constraints, which is more realistic than the scheme , e.g., in [5] where only whole system power summit is considered. To the best of our knowledge, such joint optimization for assumed two hop OFDMA network with imperfect CSI has not been studied in the literatures and is important for achieving the better overall system performance. The rest of this paper is organized as follows. Section II describes our relay-assisted OFDMA cooperative wireless
networks and formulates the problem. We consider downlink only in this work, but it can be extended further to the uplink case. In Section III, the proposed resource allocation scheme is presented. We demonstrate the benefits of our proposed algorithm in section IV and finally conclude the paper in Section V. II. S YSTEM M ODEL AND P ROBLEM F ORMULATION A. System model and Assumptions We consider our system as a two-hop downlink OFDMA relay network. The whole system consists of source(i.e. access point, AP), destination node(i.e. mobile terminal, MT) and several relays. The first hop is so called broadcast phase, where AP broadcasts information to a cluster of decode-and-forward (DF) relays. In the second hop, relays cooperate to transmit the information data to the MT, so that, e.g., spatial diversity gain can be achieved (relays are assumed to be far enough to each other). The estimated CSI is assumed to be known at MT through estimator (e.g., minimum mean square error estimator (MMSE)) and then feed back to the transmitter perfectly. We also assume that channel estimation error pertains to the amplitude of the correct channel gain, while the phase of the channel gain can be perfectly obtained. As a result, estimated channel gain with an estimation error is available to both the transmitter and the receiver. The AP acts as a central controller to carry out all resource allocation related operations based on the feedbacks from the MT. In this work, we assume there are total Z relays in the networks, and the selected relay cluster K contains K potential half-duplex relays. The presented relay-assisted cooperative OFDMA network is as shown in Fig. 1 when K = 3.
mean Gaussian random variable with variance σh˜2 . Thus, the ˆ σ 2 ). n is the imperfect CSI h is assumed to follow CN (h, ˜ h additive noise which can be also modeled as complex Gaussian random variable with variance σn2 . Therefore, the square of imperfect CSI h follows a noncentral chi-square probability density function (PDF) given by √ ˆ G+G ( ˆ ) 1 − σ2˜ GG ˆ h J f (G|G) = 2 e (3) 0 2 σh˜ σh˜4 ˆ 2 . J0 is the 0th order ˆ = |h| where we denote G = |h|2 ,G modified Bessel Function of the first kind. In our proposed system model, we suppose hi is the channel transfer function from transmitter to receiver and we assume ˆ i means the the channel is static in a time slot. For example, h s,k channel estimate from AP s to relay node (RN) k over OFDM ˆ j means the channel estimate from RN k to subcarrier i and h k,d destination d over OFDM subcarrier j. We have channel gain ˆ i |2 and second hop G ˆ j |2 . ˆ i = |h ˆ j = |h of the first hop G s,k s,k k,d k,d L is the path loss factor and the noise variance for two hops are σk2 and σd2 . The variance of related estimation error for 2 2 2 two hops are σh,k and σh,d and we assume σh˜2 = σh,k = ˜ ˜ ˜ 2 σh,d . We denote the transmit power assigned to subcarrier i ˜ i for transmitting data as P . In this work, we do not consider the direct link from AP to MT due to distance or obstacles. This assumption is practical for the case that RNs are deployed for cell extension. One RN k occupies subcarrier i in the first hop and j in the second hop. Therefore, at first hop, the data rate of the broadcast phase is determined by the minimum rate of each link between AP and selected RNs. Since transmitter only knows the CSI conditioned on the feedback of receiver, we could obtain the expected achievable throughput of the first hop as follows: M { [ ]} ∑ I i i i i |ˆ i Rs,K = min Eγs,k log(1 + ωs,k ρk Ps,k γs,k ) , (4) γs,k k∈K
i=1 Ls,k Gis,k 2 σk
Figure 1.
Wireless cooperative relay networks
B. Problem Formulation Let x be the transmit data from transmitter to receiver and P is the transmit power gain. So regardless of the path loss, the received data after estimator at receiver is √ y = h P x + n, (1)
We also define ω is the indicator whether certain subcarrier is assigned to RN k, for example,
and we have ˆ + h, ˜ h=h
(2)
ˆ is the estimated channel function and h ˜ is the where h independent estimation error which can be modeled as zero
ˆi Ls,k G
s,k i i where γs,k = and γˆs,k = . The notation 2 σk i |ˆ i Eγs,k means expectation with respect to γ conditioned s,k γs,k i on γˆs,k . M is the subcarrier set of the system that contains M subcarriers. I is the subcarrier set which contains the subcarriers that are allocated to the selected RNs at first hop. We further refer the link throughput and its lower bound interchangeably for simplicity. ρk indicates that whether RN k is chosen for subcarrier allocation, so we obtain { 1 if k is chosen for relaying, ρk = 0 otherwise.
{ i ωs,k
=
1 if i is assigned to k at first hop, 0 otherwise.
For the second hop, it is assumed that the RNs are perfectly
synchronized and transmitted at the same time. Therefore, the second hop can be viewed as a virtual MISO link. The expected throughput can be expressed as J RK,d = Eγ j
[
j γk,d k,d |ˆ
(
log 1 +
M ∑ K ∑
] j j j ) ωk,d ρk Pk,d γk,d ,
( ) I J L(P, ω, ρ, λ, µ) = Rs,K + RK,d −λs (
(5)
j=1 k=1 Lk,d Gj
Therefore, the total achieved end-to-end throughput of source s to destination d through RN set K is [10] } 1{ I J Rsd = min Rs,K , RK,d . (6) 2 Then, we can formulate our problem as max Rsd ,
(7)
subject to M ∑ K ∑
k=1 K ∑
i i ωs,k Ps,k − Ps,max )
i=1 k=1 M K ∑ ∑
λk,d (
k=1
(11) j j Pk,d − Pk,max ) ωk,d
j=1
I J −µ(Rs,K − RK,d ), j i where P = {Ps,k , Pk,d } is the set of power allocation, ω = j i {ωs,k , ωk,d } denotes the subcarrier allocation, and ρ = {ρk } is the relay assignment. The λ = {λs , λk,d } and µ = {µ} are Lagrange multipliers. Then it can be derived that λs , λk,d ≥ 0 and µ ∈ (−1, 1). The Lagrange dual function can be written as:
g(λ, µ) = max L(P, ω, ρ, λ, µ).
i ωs,k
=
(8) ∈ {0, 1}
j j ωk,d = 1, ωk,d ∈ {0, 1}
k=1
where Ps,max is the maximum transmit power of AP and Pk,max is the maximum power of RN. Therefore, our goal is to find the optimal solutions of relay, subcarrier and power allocations which satisfy the problem (7). It can be deduced that (6) can achieve maximum only when J I Rs,K = RK,d . Thus, (7) can be modified to ( ) I J arg max Rs,K + RK,d , (9) subject to conditions in (8) and I J Rs,K = RK,d .
(13)
A. Evaluating Dual Variable
≤ Pk,max
i 1, ωs,k
(12)
Since we assume the number of subcarrier is sufficient large, so that the duality gap between primal problem and dual function can be negligible [12]. consequently, we can solve the problem (7) by minimizing the dual function min g(λ, µ).
i i ≤ Ps,max Ps,k ωs,k
i=1 k=1 M ∑ j j ωk,d Pk,d j=1 K ∑
−
ˆj Lk,d G
j j k,d k,d where γk,d = and γˆk,d = . J is the σd2 σd2 subcarrier set which contains the subcarriers that are allocated j , we also to the selected RNs at second hop. For indicator ωk,d have { 1 if j is assigned to k at second hop, j ωk,d = 0 otherwise.
K M ∑ ∑
Since a dual function is always convex [11], then for example, two methods can be used to minimize g(λ, µ) with guaranteed convergence, which are subgradient method and ellipsoid method [12]. We follow the subgradient method in [12] to derive the subgradient g(λ, µ) with the optimal power allocation p∗ that will be presented in the following subsection. Algorithm 1 Evaluating Dual Variable 1: Initialize λ0 and µ0 2: 3: 4: 5:
while (!Convergance) do Obtain g(λa , µa ) at the ath iteration; Update a subgradient for λa+1 and µa+1 , by λa+1 = λa + υ a △λ and µa+1 = µa + υ a △µ; end while
where △λ = {△λs , △λ1,d , ...△λK,d }, △λs , △λk,d and △µ can be expressed as
(10)
△λs = Ps,max −
i (Ps,k )∗
i=1 k=1 M ∑
III. R ESOURCE A LLOCATION S CHEME In this section, we introduce the adaptive algorithms to solve existing problem (9). Although the resource allocation problem is combinatorial in nature with nonconvex structure, as the number of subcarriers becomes sufficient large, the dual gap tends to be zero [12]. Therefore, it can be solved in dual domain. The Lagrangian of problem (9) is [11]
M ∑ K ∑
△λk,d = Pk,max − △µ = a
I (Rs,K )∗
−
j (Pk,d )∗
(14)
j=1 J (RK,d )∗ .
Here υ is the stepsize and a is the number of iterations. The subgradient algorithm in Algorithm 1 is guaranteed to con-
verge to the optimal λ and µ. The computational complexity of Algorithm 1 is polynomial in the number of dual variable K + 1 [12]. Since (12) can be viewed as nonlinear integer programming problem, whose optimal solution requires high computational cost. Therefore, we are aiming to solve the optimization problem by solving three subproblems, which are relay selection, subcarriers and power allocation. We firstly introduce power allocation scheme. B. Power Allocation Scheme In order to obtain the optimal solution of power allocation, we are aiming to solve the problem solving problem (11) j i over variables Ps,k and Pk,d . However, from (5) and (4), we see that problem (11) involves the conditional expectation of achievable throughput with respect to estimated CSI. Applying Karush-Kuhn-Tucker (KKT) conditions [11] and equations in [14], we could obtain the optimal power allocation schemes by solving following equation numerically.
Figure 2. One example of power allocation as function of estimated channel SNR for various value of σ˜2 . h
M { } ∑ i i i i |ˆ i [log(1 + L(P, ω, ρ, λ) = min Eγs,k ωs,k ρk Ps,k γs,k )] γs,k k∈K
2 i
σ β i i ( σk2 βs,k )αis,k L k Ps,k αs,k i e s,k s,k i i Ps,k Ls,k Ps,k
( Γ −
i ) σk2 βs,k i , αs,k i Ls,k Ps,k
λs = . 1−µ (15) i where Γ(a, b) is the incomplete Gamma function. αs,k = i i (ηs,k + 1)2 /(2ηs,k + 1) is the Gamma shape parameter with i i i i 2 ˆ i + σ 2 ) is Gamma ˆ = αs,k /(G ηs,k = Gs,k /σh˜ and βs,k ˜ s,k h PDF rate parameter. Similarly, for the cooperation phase, the optimal RN power allocation is obtained by solving: j αk,d ( j Pk,d
( j j )αk,d c1 βk,d e Γ c1 βk,d j
−
j j ) αk,d , c1 βk,d
λk,d = , (16) 1+µ
j j j j where αk,d = (ηk,d + 1)2 /(2ηk,d + 1) with ηk,d = j j j j 2 2 ˆ /σ and β ˆ G = α /( G ). We have c + σ 1 = ˜ ˜ k,d k,d k,d k,d h h σd2 +
∑K
m=1,m̸=k
Pm,d Lm,d Gm,d
j Lk,d Pk,d
. Pm,d and Gm,d is the power
allocation and channel gain from relay m to MT d. By using approximation method, e.g. , in [15], we are able to obtain the power allocation with imperfect CSI. One example can be found in Fig. 2 where different value of σh˜2 is considered. We can see that when estimated error is relative small, the power allocation achieved by imperfect CSI is very close to the one when perfect CSI is assumed at AP.
i=1
+Eγ j
j γk,d k,d |ˆ
[log(1 +
j j i ωs,k ρk Pk,d γk,d )]
i=1 k=1
(
M { } ∑ i i i i |ˆ i [log(1 + −µ min Eγs,k ωs,k ρk Ps,k γs,k )] γs,k k∈K
−Eγ j
j γk,d k,d |ˆ
i=1
[log(1 +
M ∑ K ∑
) j j i ωs,k ρk Pk,d γk,d )]
i=1 k=1
+λs ( −
M ∑ K ∑
i i ωs,k Ps,k − Ps,max )
i=1 k=1 M ∑
K ∑
k=1
λk,d (
j j ωk,d Pk,d − Pk,max ).
j=1
(17) By applying KKT condition, the RN is selected according to the following rule, { [ P γ ]} s,k s,k i |ˆ i K∗ = arg max min (1 − µ∗ )Eγs,k γs,k k∈K 1 + Ps,k γs,k k [ ] P γ k,d k,d , +(1 + µ∗ )Eγ j |ˆγ j ∑K k,d k,d 1 + k Pk,d γk,d (18)
C. Optimal Relay Selection (ORS) We consider ORS in this work, unlike some traditional single relay selection algorithms in [6] and [13], as the multiple RNs selection. The proposed algorithm is to select K RNs to form a cluster that can maximize the achieved throughput in (6) based on the imperfect CSI. Assuming the subcarrier and power allocation is done, we can rewrite (11) as
M ∑ K ∑
ˆi Ls,k G
s,k i = . The channel SNR Since we know that γˆs,k 2 σk i γs,k conditioned on γˆs,k is also a non-central Chi-squared distributed random variable with PDF:
i i f (γs,k |ˆ γs,k )=
−
1 e i νs,k
i +γ i γ ˆs,k s,k νi s,k
( J0 2
√
γ ˆi γi s,k s,k (ν i )2 s,k
) (19)
j i |ˆ γk,d )= f (γs,k
1 j νk,d
j j γ ˆ +γ − k,dj k,d ν k,d
( J0
{
v ) u j u γˆ γ j 2t k,d k,d
e
(ν
j )2 k,d
ωj = (20)
j i where νs,k = σk2 /σh˜2 and νk,d = σd2 /σh˜2 . Following same procedure as power allocation, we obtain i |ˆ i Eγs,k γs,k
[ P i γi s,k s,k
]
i γi 1 + Ps,k s,k
i =ψs,k
i i ( θs,k )ψs,k i Ps,k
e
i θs,k Pi s,k
i ) ( σk2 θs,k i Γ − ψs,k , , i Ps,k
Eγ j
j γk,d k,d |ˆ
[
(21)
j i ] j Pk,d γk,d j j ( j )ψk,d c2 θk,d =ψk,d c2 θk,d e ∑K j j 1 + k Pk,d γk,d ( j j ) Γ − ψk,d , c2 θk,d , (22)
i i i i i i where ψs,k = (ζs,k + 1)2 /(2ζs,k + 1) with ζs,k = γˆs,k /νs,k j j j i i i i and θs,k = ζs,k /(ˆ γs,k + νs,k ). ψk,d = (ζk,d + 1)2 /(2ζk,d + 1) j j j j j j j with ζk,d = γˆk,d /νk,d and θk,d = ζk,d /(ˆ γk,d + νk,d ). We have ∑K j c2 = (1 + m=1,m̸=k Pm,d γm,d )/Pk,d . Pm,d and γm,d are the power allocation and channel SNR from relay m to MT. Optimal value of P can be given in (15) and (16). Thus, (18) can be viewed as multi-objective optimization problem, which aims at obtaining the trade-off of the throughputs of first hop and second hop. (18) is also the termination criteria for the whole RRA scheme. Therefore, the relay selection strategy is { 1 if k ∈ K∗ , ρk = 0 otherwise.
1 0
if j ∈ J ∗ , otherwise.
IV. P ERFORMANCE E VALUATION Simulations are presented to evaluate the performance of proposed algorithms in this section. It is assumed that five RNs are located between AP and MT, and MT is 1.8km away from AP. The Stanford University SUI-3 channel model is employed [16], in which the central frequency is 1.9GHz. Channel is assumed to be 3-tap channel and signal fading follows Rician distribution. We choose number of subcarriers N to be 32, so the duality gap can be ignored [6]. Flat quasi-static fading channels are considered, hence the channel coefficients are assumed to be constant during a complete frame, and can vary from a frame to another independently. The noise variance of the two hops are set to be 1 for simplicity. The path loss factor varies according to the different distances from RNs to AP and MT. If distance between RN and AP or RN and MT is shorter than a break point dBP = 100m, the exponent is fixed to 2, otherwise it is 3.5. The maximum transmit power of AP and RN are 40 dBm and 20 dBm respectively. An accurate MMSE estimator with estimated variance σh˜2 = 0.02 is assumed at the receiver. We demonstrate our results comparing with the performance of some other schemes: 1) Equal power allocation combined with proposed subcarrier allocation scheme and relay selection(EPA); 2) Waterfilling power allocation and proposed subcarrier allocation scheme with relay selection (Waterfilling); 3) Modified proportional allocation scheme in [5] with fairness consideration(Fairness SA);
D. Optimal Subcarrier Allocation (OSA) The goal of subcarrier allocation strategy is to assign subcarriers to a given RN that can obtain best throughput performance. Following the same procedure as the relay selection, we could obtain subcarrier allocation criteria as follows: I ∗ = arg max
∗
{
{
{ ]} [ P i γi s,k s,k i |ˆ i min (1 − µ∗ )Eγs,k γs,k i i k∈K 1 + Ps,k γs,k (23) ∗
[
J = arg max (1 + µ )E
j j γk,d |ˆ γk,d
Ls,k Gi
j j ]} Pk,d γk,d , ∑K j j 1 + k Pk,d γk,d (24) Lk,d Gj
j s,k k,d i = where channel SNR γs,k and γk,d = . 2 σk σd2 Therefore, the OSA indicator for the first hop and second hop can be expressed as { 1 if i ∈ I ∗ , ωi = 0 otherwise.
Figure 3. Impact of maximum transmit power Ps,max on system bandwidth Efficiency
Fig. 3 demonstrates the impact of maximum transmit power of AP on the system bandwidth efficiency. We denote Ds,d as the distance between AP and MT, and Ds,k as the distance between AP and RNs. In Fig. 3, we have Ds,d = 1800 m and Ds,k from 1500 m to 1600 m. The considered channel SNR at the RN k is varied from γs,k = −20dB to γs,k = −30dB and at MT d it is varied from γk,d = −15dB to γk,d = −25dB. It can be seen that the proposed scheme achieves the best
performance. The performance gain over other methods in comparison is up to 20%. It can also be noticed that if Waterfilling is used as the power allocation scheme (instead of our proposed scheme), the throughput performance is comparable with Fairness SA. Another performance gain can be seen in power consumption. We can see that with a fixed data rate requirement, our proposed scheme provides a clear power saving gain. For instance, at the level of 1.0 bit/s/Hz bandwidth efficiency our proposed scheme can reach a power saving around 3 dBm compared to the other schemes. Fig. 4 shows the impact of distance between AP and RN on the system throughput. The distance between AP and RN is normalized to distance between AP and MT and vary from 0.1 to 0.9. In Fig. 4, we set maximum AP power Ps,max = 40 dBm and maximum RN power is Pk,max = 20 dBm. From Fig. 4, we can see that the proposed algorithm obtain the highest system capacity when distance is less than 0.9 in Fig. 4. When the average normalized distance between AP and RN is around 0.9, we can find that performance difference between EPA and proposed scheme are smaller. This may due to the fact the some RNs are very close to the MT so that the achieved SNR is rather high. It can be concluded that the proposed algorithm can provide better performance gain over other existed algorithms.
Figure 4. Impact of distance between AP and RN on the system bandwidth efficiency, maximum AP power is 40 dBm, maximum RN power is 20 dBm
Fig. 5 illustrates the convergence speed of the proposed algorithm and the Fairness SA scheme. The considered Ps,max is fixed to 35 dBm andPr,max is 17.5 dBm. Our proposed algorithm reaches the steady state after several iterations, which demonstrates fast convergence speed. V. C ONCLUSION In this paper we investigated the problem of resource allocation for cooperative multi-relay assisted OFDMA networks with assumption of Imperfect CSI. The joint optimization problem for radio resource allocation was solved by addressing three subproblems including optimal selection of collaborative relays, subcarriers and power with the objective of maximizing the expected system throughput. Theoretical expressions were derived for the optimal selections. It was shown that by designing RRA scheme with imperfect CSI for different hops,
Figure 5. efficiency
Impact of the number of iteration on the system bandwidth
it is possible to reach a noticeable gain in the cell-edge throughput. In addition, the results support our theoretical analysis that the proposed scheme obtain power allocation as close to the one when perfect knowledge of CSI is considered. These were also illustrated with simulation examples. R EFERENCES [1] G. Song and Y. Li, "Cross-layer Optimization for OFDM Wireless Networks-Part 1: Theoretical Framework," IEEE Tran. on Wir. Comm., Vol.4, No.2, pp.614-624, Mar. 2005. [2] T.C-Y. Ng and W. Yu, "Joint Optimization of Relay Strategies and Resource Allocation in Cooperative Cellular Networks," IEEE Jour. on Sel. Area in Comm., Vol.25, No.2, pp.328-339 ,Feb. 2007. [3] S-K. Kim, X. Wang, and M. Madihian, "Optimal Resource Allocation in Multi-hop OFDMA Wireless Networks with Cooperative Relay," IEEE Tran. on Wir. Comm., Vol.7, No.5, pp.1833-1838, May 2008. [4] H. Jeong and J. H. Lee, "Adaptive Relay Selection for Regenerative OFDMA Relay Networks with Fairness Constraints," Proceedings of IEEE VTC 2008-fall, Calgary, Canada, Sep. 2008. [5] M. Tang, et al, "Dynamic Resource Allocation with Threshold in OFDMA-based Relay Networks," Proceedings of IEEE VTC 2010-Spring, Taipei, Taiwan, May 2010. [6] W. Dang, et al, "Subcarrier-Pair Based Resource Allocation for Cooperative Multi-Relay OFDM systems," IEEE Tran. on Wir. Comm., Vol.9, No.5, pp. 1640-1649, May 2010. [7] M. K. Awad, et al, "A Dual-Decomposition-Based Resource Allocation OFDMA Networks with Impefect CSI," IEEE Tran. on Veh. Tech., Vol.59, No.5, pp. 2394-2403, May 2010. [8] M. Ahmad and M. Assaad, "Joint Resource Optimization and Relay Selection for cooperative Cellular Networks with Imperfect Channel Knowledge," Proceedings of 11th SPAWC 2010, Morocco, June 2010. [9] J-H. Wang, et al, "Resource Allocation and Relay Selection for cooperative OFDMA Networks with Imperfect Channel Information," Proceedings of 12th IEEE SPAWC 2011, San Francisco, CA, June 2011. [10] Z. Chang and T. Ristaniemi, "Radio Reource Allocation for Cooperative Relay-assisted OFDMA Wireless Networks," Proceedings of 3rd IWCLD, Rennes, France, Nov. 2011. [11] S. Boyd, and L. Vandenberghe, Convex Optimization , Cambridge UK: Cambridge University Press, 2004. [12] W. Yu and R. Lui, "Dual Methods for Nonconvex Spectrum Optimization in Multicarrier System," IEEE Tran. on Comm., Vol.54, No.7, pp.13101322, Jul. 2006. [13] A. S. Ibrahim, et al, "Cooperative Communications with RelaySelection: When to Cooperate and Whom to Cooperate With," IEEE Transaction on Communications, Vol.7, No.7, pp.2814-2827, July 2008. [14] I. Gradshtein, I. Ryzhik, D. Zwillinger and A. Jeffrey, Table of Intergrals, Serie, and Products. New York: Academic, 2007. [15] A. Lauschke, Convergence Acceleration with Canonical Contractions of Continued Fractions: Incomplete Gamma Function [Online]. available: http://216.80.120.13:8080/webMathematica/LC/gamma.jsp [16] R. Jain, Channel Models A Tutorial1 v1.0 , Feb. 2007.