Linear Processing Design of Amplify-and-Forward Relays for Maximizing the System Throughput

. In this paper, firstly, we study the linear processing of amplify-and-forward (AF) relays for the multiple relays multiple users scenario. We regard all relays as one special “relay”, and then the subcarrier pairing, relay selection and channel assignment can be seen as a linear processing of the special “relay”. Under fixed power allocation, the linear processing of AF relays can be regarded as a permutation matrix. Employing the partitioned matrix, we propose an optimal linear processing design for AF relays to find the optimal permutation matrix based on the sorting of the received SNR over the subcarriers from BS to relays and from relays to users,respectively. Then, we prove theoptimality of theproposed linear processing scheme. Through theproposed linear processing scheme, we can obtain the optimal subcarrier paring, relay selection and channel assignment under given power allocation in polynomial time. Finally, we propose an iterative algorithm based on the proposed linear processing scheme and Lagrange dual domain method to jointly optimize the joint optimization problem involving the subcarrier paring, relay selection, channel assignment and power allocation. Simulation results illustrate that the proposed algorithm can achieve a perfect performance.


Introduction
Relay-assisted cooperative communication and orthogonal frequency-division multiplexing (OFDM) are core technologies for the next generation mobile communication system which attract tremendous attentions.Relay-assist cooperative communication can effectively extend communication coverage, and reduce the consumption of transmit power.It can also improve overall system throughput due to its exploiting spatial diversity and combating channel fading [1,2].OFDM divides a wideband channel into some orthogonal narrow band subcarriers.Then, the fast data flow can be transformed to a set of slow data flows.Furthermore, inner-symbol interference can be eliminated significantly.OFDM can also increase the flexibility for coding and modulation.Therefore, relay-assisted cooperative communication and OFDM are adopted by a lot of communication standards, such as 3GPP and IEEE 802.16e.
There are two common categories of relay strategies: amplify-and-forward (AF) and decode-and forward (DF).The DF relay can decode and re-encode the received signals, and then retransmit these signals to destination node.Different from DF relay, the AF relay amplifies the received signals linearly and forwards these signals to destination node without decoding.The AF strategy is the most practical relay approach as the result of its low complexity transceiver design [3].Moreover, AF relays are more transparent to adaptive modulation techniques than DF relays.
Due to the independent fading on each subcarrier, the incoming and outgoing subcarrier should be matched carefully to maximize the overall system throughput.This problem is also called subcarrier paring.Recently, subcarrier paring has gained a lot of attentions, such as [4][5][6].Reference [4] proposes an ordered subcarrier pairing (OSP) method, in which it is proved that subcarrier paring according to the channel state is optimality for equal power allocation.The authors in [5] further prove that the OSP is still optimality for optimal power allocation.Reference [6] considers a linear processing design at the relay for amplified-and-forward (AF) relaying communication system without direct path and with direct path under fixed power gain.Different from [4,5], [6] dose not only consider the cooperative mode without direct path, but also study the the cooperative mode with direct path.In [6], the authors separate the processing structure into 2 Journal of Electrical and Computer Engineering two components which are the power amplification matrix and the unitary linear processing matrix, respectively.They show the optimal unitary processing matrix is of permutation structure.
The joint resource allocation for cooperative communication has been also attracted tremendous attentions.For the scenario with multiple relays multiple users, the joint optimization problem involves subcarrier pairing, relay selection, channel assignment and power allocation.Due to the combinational nature of the subcarrier pairing, relay selection and channel assignment, the joint optimization problem is difficult to solve when the number of subcarriers, relays and users are large.Therefore, most previous works only consider a subset of these issues.
In this paper, firstly, we exploit the linear processing scheme for multiple relays multiple users network to maximize the system throughput.We regard all relays as one special "relay", and the subcarrier pairing, channel assignment and relay selection can be seen as the linear processing of the special "relay" under fixed power allocation.It can be shown that the optimal linear processing matrix of the special "relay" is a promotion of the permutation matrix.To maximize the system throughput under fixed power allocation, the optimal linear processing matrix should be found out.To this end, we propose a linear processing scheme for the AF relays under fixed power allocation which can find out the optimal permutation matrix, and then we prove the optimality of this linear processing design.It is interesting that, through the proposed linear processing scheme, we can obtain the optimal subcarrier paring, relay selection and channel assignment with fixed power allocation in accordance with the ordered SNR over the incoming and outgoing channels in polynomial time.To the best of our knowledge, this paper is the first work to investigate the linear processing of the AF relays for the multiple relays multiple users network.Finally, based on the proposed linear process scheme and Lagrange dual domain method, we propose an iterative algorithm to solve the joint optimization problem involving subcarrier pairing, relay selection, channel assignment as well as power allocation.Simulation results illustrate that the proposed iterative algorithm is effective to find out an asymptotically optimal solution in polynomial time.
The rest of the paper is organized as follows.In Section 2, the previous related works are reviewed.We describe the system model in Section 3. The optimal linear processing under fixed power gain is introduced in Section 4. In Section 5, we extend the linear process scheme to the scenario with multiple relays multiple users.In Section 6, we describe the iterative algorithm based on the proposed linear process scheme and Lagrange dual domain method to solve the joint resource allocation problem.In Section 7, simulation results are provided to evaluate the performance of the proposed algorithm.Finally, we conclude this paper in Section 8.

Related Work
The resource allocation for the scenario with single BS, multiple relays and multiple users is a joint optimization problem which involves subcarrier paring, relay selection, channel assignment and power allocation.This joint optimization problem can be formulated as a mixed-integer programming due to the combinatorial nature of the subcarrier paring, relay selection and channel assignment.It becomes more and more intractable with the increase of the number of subcarriers, relays and users because the problem of subcarrier pairing, relay selection and channel assignment is NP hard [7].To decrease the computational complexity, most previous works consider a subset of the joint optimization problems.
References [8][9][10][11] consider the scenario with one BS, one relay and one user.Reference [8] jointly optimizes the subcarrier pairing and power allocation, and a mixed-integer programming problem is formulated.Then, the authors in [8] transform the mixed integer programming problem into a convex optimization through continuous relaxation.[9] proposes a hybrid scheme in which the full-duplex and halfduplex relaying modes can be switched opportunistically to obtain a tradeoff between spectral efficiency and selfinterference.For DF relay strategy, if the ratio of the time allocation in the first transmission phase over the whole period is not large enough, the DF relay can not decode the received signals accurately.It is necessarily to optimize the the ratio of the time allocation in the first transmission phase over the whole period.Therefore, reference [10] optimizes not only the power allocation but also time allocation for DF relay strategy to minimize the outage probability.The perfect channel state information (CSI) at the source is impractical, thus it is reasonable to study the resource allocation problem with limited feedback [11].
The resource allocation for the scenario with one BS, one relay and multiple users is studied by reference [7].The authors in [7] propose an optimal algorithm for the joint resource allocation including subcarrier pairing, channel assignment and power allocation through transforming the original problem into some simple linear programming problems.In [12], a combination of transmit-receive weights and Tomlinson-Harashima precoding is used to cancel the interference, and then a low-complexity power allocation is proposed to achieve data rate fairness among different users.
The cooperative communication scenario with one BS, multiple relays and multiple users is studied by [13,14], and the direct transmission mode is also considered at the same time.The authors in [13] transform the combinational optimization problem of subcarrier paring, channel assignment, relay selection and transmission mode selection into a minimum cost network flow problem, and apply the linear optimal distribution algorithm to solve the minimum cost network flow problem.Then, the Lagrange dual domain method is used to solve the power allocation.Different from [13], reference [14] only considers the relay selection and subcarrier assignment for multiuser cooperative network to decrease the computational complexity.In the first, [14] simplifies the original problem into a new problem through decreasing the number of variables which is easily to handle.Then, branch-and-cut is introduced to optimize the new problem.
Referencee [15] considers the problem of minimizing the total power consumption in a two-hop single-relay MIMO network with QoS requirements.In [15], a nonconvex power allocation problem is approximated with a convex problem, and then the convex problem is computed in closed-form through a multistep procedure.[16] studies the problem of resource allocation in relay-enhanced bidirectional MIMO-OFDM networks to minimize the total power consumption.The authors in [16] propose a green resource allocation scheme to jointly optimize the subchannel assignment, power allocation and phase duration assignment, in which both the separate-downlink (DL)-and-uplink (UL) and mixed-DL-and-UL relaying assignments and the linear block diagonalization (LBD) technique are adopted.Reference [17] study the optimal power allocation structure for the multiple relays network.The authors in [17] show that the power allocation at BS and each relay follow a matching structure.The cooperative communication with virtual MIMO has also attracted a lot of attentions, such as [18].Reference [18] studies a wireless network where multiple users cooperate with each other.The authors in [18] propose a new auction-based power allocation framework with multiple auctioneers and multiple bidders to maximize the weighted sum-rates of the users.
Recently, energy efficient (EE) resource allocation for MIMO cooperative communication has become attractive, such as [22][23][24][25][26][27].Reference [22] considers the energy efficient joint source-relay power allocation problem for MIMO AF relaying system.In [22], the objective of optimization is the number of the bits per second per hertz per Joule with the guarantee of the minimum spectral efficiency, and the objective function and the constraint are not convex.Firstly, [22] transforms the original problem into a pseudoconvex problem by employing the high signal-to-noise ratio (SNR) approximation.Then, the authors in [22] propose a relaxation method according to the Jensen inequality to solve the pseudo-convex problem.[23] studies the energy-spectral efficiency (EE-SE) trade-off of the uplink of a multi-user cellular virtual MIMO system.Finally, a heuristic resource allocation algorithm is proposed to optimize the EE-SE tradeoff in [23].The energy-efficient resource allocation for OFDMA cellular networks with user cooperation is studied by [24].In [24], a mixed-integer nonlinear programming problem is formulated, and then an optimal algorithm is proposed to solve the problem.

System Model
In this paper, for matrix M, we let Tr{M}, M † and |M| denote the trace, conjugate transpose and determinant of M, respectively.I denote the identity matrix.diag{M 1 , . . ., M 푛 } represents the block diagonal matrix with M 1 , . . ., M 푛 on its main diagonal.We assume that the perfect CSI is obtained by the technology of channel estimation, and then the accurate channel gain is known.
In this paper, we intend to exploit the optimal linear processing of relays in OFDMA AF-based relaying network consisting multiple relays multiple users, similar as Figure 1.For easy to exposition, we first study the linear processing scheme for the scenario with multiple relays and single user in this section.Then, we will show that our proposed linear processing scheme can be readily extended to the multiple relays multiple users scenario in the next section.In this work, we assume there are  subcarriers and  AF relays.The halfduplex relay is considered, that is, each relay can not transmit and receive the signals in the same subcarrier at the same time.The transmission duration is divided into two equal time slot.Let H 푖 푠푟 ∈ C 푁×푁 and H 푖 푟푑 ∈ C 푁×푁 denote the channel gain matrices between BS and relay  and between relay  and the user, respectively.
In the second time slot, the AF relays process the received signals linearly and then retransmit these signals to the user.Since the relays cannot share their received signals with each other, we let W = diag{W 1 , . . ., W 푀 } ∈ C 푀푁×푀푁 denote the linear processing matrix of all relays.The processed signals at relay  can be represented as x [푖]  푟 = W 푖 y [푖] 푟 .Then, the overall relay signals can be given as the combined channel gain matrix between all relays and the user.Then, the received signals at user can be written as: where n 푑 ∈ C 푀푁×1 is the AWGN, with n 푑 ∼ CN(0,  2 푑 I) at the user in the second time slot.
In the second time slot, the relays combine the received signals linearly and then retransmit the amplified version of the processed signal to the user.Therefore, the processing matrix W can be represented as W = D 푟 U, where U ∈ C 푀푁×푀푁 is a special permutation matrix in which there are only  rows and  columns with nonzero entry, which determines how to combine the incoming signals linearly and D 푟 = diag{D 1 푟 , . . ., D 푀 푟 } ∈ C 푀푁×푀푁 is the power amplification matrix of all relays, where D 푖 푟 = diag{ 푖,1 , . . .,  푖,푁 } and d 푖,푛 means the power amplification factor for the th subcarrier at relay .In this section, we intend to investigate the linear processing matrix U with fixed D 푠 and D 푟 .That is, our goal is to investigate what form of the linear processing scheme leads to the optimal relay selection and subcarrier pairing.Since we are interesting in the permutation matrix U, the received signals at the user can be rewritten as where are the functions of matrix U, and they can be seen as the equivalent channel matrix and noise term.
In fact, there are  relays and  subcarriers so that there are  incoming subchannels and  outgoing subchannels.To avoid interference, we can only select  incoming subchannels and  outgoing subchannels (each subchannel can only be assigned to one relay node at each time slot).Therefore, U is a non-full rank permutation matrix in which there only  rows and  columns with nonzero entry.

Overall System Throughput
As aforementioned, the end-to-end achievable rate can be given by where R 푛 ≜ E[ñ(U)ñ † (U)] represents the covariance matrix of the equivalent noise term.The factor 1/2 in (3) accounts for the two time slots required to complete a transmission.Due to the maximum transmission power limit at BS and each relay , we have the following peak power constraints: where P 푠 and P 푖 are the maximum transmission power of BS and relay , respectively.
In this section, our objective is to find the optimal U * to maximize the system achievable (U), which can be formulated as the following optimization problem: For the scenario with multiple relays and single user, relay selection and subcarrier pairing can be essentially represented as a kind of U. In this case, U ∈ C 푀푁×푀푁 is a non-full rank permutation matrix in which there are only  rows and  columns with nonzero entries.Next, we give our method to find the optimal matrix U * .In accordance with |I + AB| = |I + BA|, (6) can be rewritten as: Substituting The th diagonal entries of F and T are given by: 푑 be the received SNR from BS to relay and from relay to user respectively, and then there are It is clear that { 푗 } and { 푗 } have the same sorting order with {SNR 푟푑,푗 } and {SNR 푠푟,푗 }.Therefore,  푗 and  푗 are equivalent to SNR 푟푑,푗 and SNR 푠푟,푗 , respectively.To maximize (9), in the first, we give the following Lemma.

Proof. Please refer to Appendix
For the scenario with one relay and one user, U is a full rank permutation matrix since there is only subcarrier pairing that need to be optimized.According to Lemma 1, we can readily obtained the optimal subcarrier pairing with fixed power allocation, which is that the incoming subcarrier with the best SNR is paired with the outgoing subcarrier with the best SNR, and the incoming subcarrier with the worst SNR is paired with the outgoing subcarrier with the worst SNR, and so on.However, when the number of relays is more than one, the best incoming subcarrier and the best outgoing subcarrier may do not belong to the same relay which makes these two subcarrier can not be matched.Therefore, the case with multiple relays differs significantly from the case with single relay.To exploit the optimal linear processing scheme for the scenario with multiple relays, in the following, we first give two constraints to avoid interference and make the linear processing scheme be in line with the reality, and then we will give a Theorem to exploit the optimal linear processing for multiple relays network.
Constraint 2. Only  orthogonal subchannels can be selected at each hop to avoid co-channel interference.Constraint 3. Two subcarriers which are matched to one pair should belong to the same relay node.Indeed, Constraint 2 means that each subcarrier can be only assigned to one relay at each transmission slot to avoid interference.
Next, we give a Theorem as follows: For any other permutation matrix Ũ which also satisfies Constraints 2 and 3, we assume the corresponding subcarrier pair set is { g(1) , g(2) , . . ., g(푁) }.From Lemma 1, there is According to the selection criterion of Theorem 4, there is at least one g(푖 ) ≤  (푖) .So, there is Therefore, the permutation matrix U * is optimal for (9).
Theorem 4 gives an optimal linear processing scheme of relays corresponding to one optimal relay selection and subcarrier pairing scheme for the scenario with multiple relays.In fact, as aforementioned, when there is only one relay, the linear processing matrix "U" is a full rank permutation matrix which only corresponds to the subcarrier pairing.In accordance with Lemma 1, the optimal pairing strategy is based on the sorted received SNRs.However, for the scenario with multiple relays considered in this paper, the processing matrix "U" is a non-full rank matrix since each subcarrier can only be allocated to one relay node to avoid interference.Theorem 4 shows that the optimal relay selection and subcarrier pairing strategy is based on the sorted products of the incoming and outgoing subchannel's SNRs.
In fact, there are  ×  "possible" incoming channels between BS and all relays and  ×  "possible" outgoing channels between all relays and the user, and we can only select  incoming and  outgoing channels.There are () 2 combinations of incoming and outgoing channels.Each combination corresponds to one product of one incoming and one outgoing subchannel's SNR.Theorem 4 means the optimal scheme is that, firstly, these combinations are ordered according to their corresponding products.Then,  best combinations satisfying Constraints 2 and 3 would be selected from the ordered () 2 combinations.
1  21  1 22 = 2,   1  12 ,  1 21 } are selected because another channels do not satisfy Constraints 2 and 3. Thus, the optimal transmission paths are (2, 1, 1, 1) and (1, 2, 1, 1), which mean that the subcarrier pair (2, 1) is assigned to user 1 with the help of relay 1 and the subcarrier pair (1, 2) is assigned to user 1 with the help of relay 1, respectively.The corresponding permutation U is the form as follows: Note that it is possible that there exists tie, for instance, there are two combinations of subcarriers with the same product and both these two combinations satisfy Constraints 2 and 3.For this case, arbitrary tie-breaking will be adopted which does not change the optimality of the algorithm.

The Scenario with Multiple Relays Multiple Users
In this section, we aim to extend the linear processing scheme proposed in the previous section to the scenario with multiple relays multiple users.We assume there are one BS,  relays,  users and  subcarriers.To make the computation feasible, in the first time slot, we denote the channel matrix between BS and all relays, the transmit power coefficient matrix and the transmission symbol vector as H 푠푟 = diag{H 푠푟 , H 푠푟 , . .Note that U is also a special permutation matrix, in which there are only  columns and  rows with nonzero entries.
Next, let 푟푑 , D 2 푟푑 , . . ., D 퐾 푟푑 } ∈ C 푀푁퐾×푀푁퐾 denote the combined channel gain matrix from all relays to all users and the combined power amplification matrix from all relays to all users, respectively.The received signals of all users can be given by: y 푑 = H 푟푑 D 푟푑 s + n 푑 , where n 푑 = [n [1] † 푑 , . . ., n [퐾] † 푑 ] ∈ C 푀푁퐾×1 represents the AWGN at users.Finally, the received signals at all users can be expressed as: The achievable rate in this scenario with multiple relays multiple users can be given by where H(U) and R 푛 have the similar form as that in previous section.So far, it is easy to prove that the Theorem 4 can be extended to the scenario with multiple relays multiple users.Thus, for the multiple relays multiple users network, when the power allocation is fixed, we can obtain the optimal permutation U which corresponds to the optimal subcarrier pairing, relay selection and channel assignment through the linear processing scheme proposed in this paper.
The computational complexity of the proposed linear processing scheme for the scenario with multiple relays multiple users is given by O(() 2  log() 2  + () 2 ).In general,  and  are in the order of O( √ ).Therefore, the overall complexity of the proposed linear process scheme is O((7/2) 3 √  log  +  7/2 ).Note that the computational complexity of the proposed algorithm is less than most of the algorithms which solve the similar problems, such as network flow method.
Remark 5. Reference [20] studies the resource allocation of MIMO-OFDM relaying system and achieves excellent performance.However, it is worth nothing that [20] considers the power allocation and subcarrier assignment of relay system, and studies resource allocation problem from scalar optimization.In this work, we consider the linear processing scheme of AF relay system, which is a signal processing problem.We formulate the problem as a matrix optimization problem, and solve it using the method of matrix analysis.

Joint Optimization with Power Allocation
The relay selection and subcarrier pairing under fixed power allocation have been optimized in the previous section.In this section, we intend to tackle the joint resource allocation involving relay selection, subcarrier pairing, channel assignment as well as power allocation using the Lagrange dual domain method and the proposed linear processing scheme.
Firstly, we give a set of binary assignment variables corresponding to the permutation matrix U: 푚푛푖푘 : indicates whether subcarrier pairing (, ) is assigned to relay  to help user .Note that  푚푛푖푘 is dependent on the permutation matrix U.In the aforementioned section, we have optimized the binary assignment variables  푚푛푖푘 which indicates the subcarrier pairing, relay selection and channel-user assignment under fixed power allocation through optimizing the linear processing of the In this section, we take into account to optimize the power allocation.
Then, the joint optimization problem can be formulated as: where constraint (25) and (26) mean each subcarrier can be assigned at most one node (user or relay) at each time slot to avoid interference.( 27) and ( 28) mean the peak power constraints of BS and relays which are converted from ( 4) and ( 5), respectively.Problem P1 is a mixed-integer nonlinear programming problem.Furthermore, P1 is nonconvex, which is difficult to tackle in general.However, it has been shown by [13] that the duality gap is close to zero as the number of subcarriers goes to infinity.Therefore, the Lagrange dual domain method can be used to solve problem P1.
It is worth nothing that the logarithm function in ( 23) is not concave.To make it concave, we remove the "1" in the denominator.Then, the original logarithm function becomes: It is clear that (30) is an upper bound of (23).The authors in [24] have shown that this upper bound is tight at moderate high SNR.
The dual function of P1 can be formulated as: (U, P, ) . (32) 6.1.Optimizing U and P with Given .In this subsection, we optimize U and P with given  to compute the dual function.
In the first, we optimize P with fixed U. The optimal P can be readily obtained using KKT conditions as follows: where [] + = max{, 0}.
In the previous section, we can obtain the optimal U * with fixed power allocation at the BS and the power amplification at the relays using the proposed linear processing scheme.Substituting U * into (24), we can obtaine the corresponding optimal P.Then, an iterative algorithm can be designed.The details will be shown in the following subsection.

Optimizing the Dual Function.
In the aforementioned subsection, the dual function of the problem P1 has been Initialize the Lagrange multiplier  (0) and power allocation P (0) repeat set  ←  + 1 For given  (푡) and P (푡) , obtain the optimal U (푡) () using the proposed linear process scheme.Then, according to the U (푡) (), obtain P (t+1) () using the Lagrange dual domain method Update  through  푡+1 = [ 푡 − ΔV 푡 ] + , where V 푡 is the step size at the th iteration.until The convergence of min 푡 ( 푡 ) Algorithm 1 computed.Now, we minimize the dual problem using subgradient method: min  () Any subgradient-based method can be adopted to minimize ().Then, a subgradient at the point  can be given by Now, we propose an iterative algorithm to optimize the variables U and P. Firstly, Lagrangian multiplier  and power allocation P are initialized, and then U() can be obtained by the proposed linear processing scheme.In accordance with U(), we can obtain corresponding power allocation P().If P() satisfies constraints ( 27) and ( 28), the solution U() and P() are global optimal solution.otherwise, iterative algorithm would continue to refine the solution.The details of the iterative algorithm is given in Algorithm 1.

Simulation Results
In this section, we evaluate the performance of the proposed algorithm using simulations.We assume that the BS is located in the center of the cell, and all users are placed randomly, and the relays are uniformly distributed on the centered circle.Cost231-Hata model is adopted as the path loss model, and the small scale fading is modeled as Rayleigh fading.The simulation parameters are summarized in Table 1.

Comparing with Benchmarks.
In this subsection, we compare the performance of the proposed iterative algorithm with that of some common benchmarks as follows: 7.1.1.No Pairing.In this scheme, the subcarrier pairing is not considered, that is, the subcarrier in the first hop is the same with that in the second hop.Thus, we should only optimize the channel assignment, relay selection and power allocation.

No PA.
In this scheme, we assume each subcarrier is allocated uniform power.Therefore, the original problem  becomes a combinatorial problem since the original problem only consists of the subcarrier pairing and relay selection.

Max Channel Gain.
In this scheme, there is no subcarrier pairing and power allocation.That is, the subcarrier in the first hop and in the second hop are the same.The channel is assigned according to the channel gain over the second hop.
Figure 2 shows the system throughput versus the number of users.The number of subcarriers and relays are set to  = 16,  = 4, respectively.The number of users changes from 2 to 12.It can be seen from Figure 2 that the system  throughput increases with the increase of the number of users since the subcarriers have higher possibility to be allocated to the users who have better channel gain.Meanwhile, the proposed algorithm achieve 15 percent system throughput improvement than the No Pairing scheme.Figure 3 illustrates the system throughput versus the number of subcarriers.It can be shown that the system throughput increases with the increase of the number of users.The number of subcarriers changes from 8 to 32.This result can be easily explained as that the system throughput is proportional to the system bandwidth.
Figure 4 depicts the system throughput versus the total transmit power.The number of subcarriers, relays and users are set to  = 16,  = 2 and  = 4, respectively.Figure 3 shows that the proposed algorithm achieves more system throughput than these benchmark schemes.The result illustrates that it is necessary to jointly optimize the subcarrier pairing, relays selection and power allocation, and the proposed algorithm is effective.

Convergence Performance and System Capacity.
In this subsection, we evaluate the convergence performance of the proposed algorithm and the system capacity versus the number of subcarriers.
Figure 5 shows the system throughput of the proposed algorithm versus the iteration indices.From Figure 5, it can be seen that, since the proposed linear processing scheme is polynomial time algorithm, the proposed algorithm have satisfactory convergence performance, and it converges within 50 iterations.
Figure 6 depicts the system throughput of the proposed algorithm versus the number of subcarriers with different total transmit power.The result shows that the system throughput can be improved significantly with the increase of the system bandwidth and SNR.

Conclusion
This paper investigates the joint resource allocation problem in multiple relays multiple users MIMO AF relaying network.The joint resource allocation problem involves subcarrier pairing, relay selection, channel assignment and power allocation which is difficult to tackle due to the combinatorial nature of the subcarrier pairing, relay selection and channel assignment.The subcarrier pairing, relay selection and channel assignment can be regarded as a linear processing of the relays which can be represented by a permutation matrix.We exploit the linear processing scheme for the scenario with multiple relays multiple users, and then propose an optimal linear processing design with fixed power gain to find the optimal permutation matrix which corresponds to  the optimal subcarrier pairing, relay selection and channel assignment.Through this linear processing design, we can obtain the optimal subcarrier pairing, relay selection and channel assignment with fixed power allocation.Then, the Lagrange dual domain method is adopted to solve the power allocation.Finally, an iterative algorithm is proposed to solve the combinatorial optimization problem and power allocation alternately.
The massive MIMO technology is one of the key topics of 5G wireless networks [28].An issue with this technique is to eliminate the pilot population.In the future works, we will pay more attention to the resource allocation and pilot population elimination for the relay-enhanced massive MIMO system.

Figure 2 :
Figure 2: Operator's profits versus different number of users.

Figure 5 :
Figure 5: Operator's profits versus different number of users.

Figure 6 :
Figure 6: Operator's profits versus different number of users.