Abstract

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 the optimality of the proposed linear processing scheme. Through the proposed 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.

1. 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 [46]. 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 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.

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 [811] 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 half-duplex relaying modes can be switched opportunistically to obtain a tradeoff between spectral efficiency and self-interference. 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.

There are also a lot of works on the resource allocation for the multiple-input and multiple-output (MIMO) cooperative communication networks, such as [3, 1521]. 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 [2227]. 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 pseudo-convex 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.

3. System Model

In this paper, for matrix , we let , and denote the trace, conjugate transpose and determinant of , respectively. denote the identity matrix. represents the block diagonal matrix with 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 half-duplex 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 and denote the channel gain matrices between BS and relay and between relay and the user, respectively.

The received signals at the relay in the first time slot are given by , , where is the transmission symbol vector and is the transmission power coefficient matrix from BS to relay , and represents the AWGN at relay , with ~. Next, we let denote the combined relay channel matrix between BS and all relays. Similarly, we let , and . Thus, the combined received signal of all relays can be given as .

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 denote the linear processing matrix of all relays. The processed signals at relay can be represented as . Then, the overall relay signals can be given as , where .

Also, we let represent the combined channel gain matrix between all relays and the user. Then, the received signals at user can be written as:where is the AWGN, with 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 can be represented as , where 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 is the power amplification matrix of all relays, where and means the power amplification factor for the th subcarrier at relay . In this section, we intend to investigate the linear processing matrix with fixed and . 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 , the received signals at the user can be rewritten aswhere and are the functions of matrix , 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, is a non-full rank permutation matrix in which there only rows and columns with nonzero entry.

4. Overall System Throughput

As aforementioned, the end-to-end achievable rate can be given bywhere represents the covariance matrix of the equivalent noise term. The factor 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 and are the maximum transmission power of BS and relay , respectively.

In this section, our objective is to find the optimal to maximize the system achievable , 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 . In this case, 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 . In accordance with , (6) can be rewritten as:Substituting and into (7), we haveLet and , and then (8) can be rewritten asThe th diagonal entries of and are given by: Let and be the received 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 and . Therefore, and are equivalent to and , respectively.

To maximize (9), in the first, we give the following .

Lemma 1. Let and be two diagonal matrices. We let and be the ordered sequences of and in descending order, respectively. The optimal to maximize (9) is as follows: That is, is the permutation matrix that matches and .

Proof. Please refer to Appendix

For the scenario with one relay and one user, 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:

Theorem 4. Let and . Let , where . It is clear that the cardinality of is since there are combinations of , and . Let be a ordered sequence in descending order. Under Constraints 2 and 3, the optimal to maximize (9) is as follows:

Proof. We assume the optimal incoming and outgoing channels obtained from Theorem 4 are and respectively, which correspond to a set of , where , and so on. Without loss of generality, we assume are ordered sequence with descending order. Then, the corresponding permutation matrix satisfies:For any other permutation matrix which also satisfies Constraints 2 and 3, we assume the corresponding subcarrier pair set is . From Lemma 1, there isAccording to the selection criterion of Theorem 4, there is at least one . So, there isTherefore, the permutation matrix 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 s. 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 s.

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 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 combinations.

For example, if the number of subcarriers, relays and users are 2, 2 and 1 respectively, there are 4 “possible” incoming channels: , and 4 “possible” outgoing channels: . There are 16 combinations of the incoming and outgoing channels. Assume , , , , , , , . Then, , , , , , , , , , , , , , , , . Therefore, and are selected because another channels do not satisfy Constraints 2 and 3. Thus, the optimal transmission paths are and , which mean that the subcarrier pair is assigned to user with the help of relay and the subcarrier pair is assigned to user with the help of relay , respectively. The corresponding permutation 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.

5. 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    and , respectively. Then, the combined received signals of all relays can be given by , where . Similarly, the processing matrix of all relays is extended to , and the linear processing matrix . Note that is also a special permutation matrix, in which there are only columns and rows with nonzero entries.

We let denote the combined channel gain matrix from all relays to user , where is the channel gain from relay to user , and let denote the power amplification matrix from all relays to user , where denote the power amplification matrix from relay to user . denotes the power amplification factor from relay to user over subcarrier .

Next, let and 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: , where 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 bywhere and 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 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 . In general, and are in the order of . Therefore, the overall complexity of the proposed linear process scheme is . 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.

6. 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 :

: indicates whether subcarrier pairing is assigned to relay to help user . Note that is dependent on the permutation matrix . 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 relays. In this section, we take into account to optimize the power allocation.

Since both and are diagonal matrices, we can readily convert (19) into scalar form as follows:In accordance with the definition of , can be represented as where is the transmit power of relay to user over subcarrier in the second time slot. Let , and . Substituting , and into (20), and then (20) can be rewritten as:

Our objective is to jointly optimize the subcarrier pairing, relay selection, channel assignment and power allocation to maximize the overall system throughput. Let , andThen, 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 is a mixed-integer nonlinear programming problem. Furthermore, 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 .

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 .

Next, we solve the problem using Lagrange dual domain method. The Lagrangian function of is:where and , are the Lagrangian multipliers associated with the constraint (27) and (28).

The dual function of can be formulated as:

6.1. Optimizing and with Given

In this subsection, we optimize and with given to compute the dual function.

In the first, we optimize with fixed . The optimal can be readily obtained using KKT conditions as follows:where .

In the previous section, we can obtain the optimal with fixed power allocation at the BS and the power amplification at the relays using the proposed linear processing scheme. Substituting into (24), we can obtaine the corresponding optimal . Then, an iterative algorithm can be designed. The details will be shown in the following subsection.

6.2. Optimizing the Dual Function

In the aforementioned subsection, the dual function of the problem has been computed. Now, we minimize the dual problem using subgradient method:

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 and . Firstly, Lagrangian multiplier and power allocation are initialized, and then can be obtained by the proposed linear processing scheme. In accordance with , we can obtain corresponding power allocation . If satisfies constraints (27) and (28), the solution and are global optimal solution. otherwise, iterative algorithm would continue to refine the solution. The details of the iterative algorithm is given in Algorithm 1.

Initialize the Lagrange multiplier and power allocation
repeat
set
For given and , obtain the optimal using the proposed linear process scheme. Then, according to the ,
obtain using the Lagrange dual domain method
Update through , where is the step size at the th iteration.
until The convergence of

7. 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.

7.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.

7.1.2. 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.

7.1.3. 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 , , 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 , and , 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.

7.2. 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.

8. 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.

Appendix

Proof. In accordance with the property of the determinant , the objective function (9) can be rewritten as:Since is not a function of , we only need to consider the second determinant:By the property of determinantwhere and are the eigenvalues of and , respectively, sorted in ascending order, the equality is reached when and are both diagonal with the diagonal entries being inverse-order matched. Then, we haveThe equality is reached when is the permutation matrix such that the entries of the ordered sequences and are one-to-one matched.

Conflicts of Interest

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