Selection based inter-relay power allocation for multiple AF single-antenna relays

This article deals with inter-relay power allocation for multiple-input multiple-output systems with multiple singleantenna relays. It is difficult to derive an optimal solution in a closed form for the case of multiple multi-antenna relays as well as the case of multiple single-antenna relays. In this article, we propose an approximate solution which is effective only for the case of multiple single-antenna relays. A key contribution of this article is a low complexity inter-relay power allocation method which is based on relay selection. This approach also reduces the feedback information of the gain factors of all the relays. An incremental greedy search algorithm is also proposed to further reduce the complexity of the relay selection process with negligible performance degradation. Simulations indicate that the performance is comparable to the optimal exhaustive search algorithm.


Introduction
The relay channel was first studied by Meulen [1] and Cover and Gamal [2], which ignited research interest on the relay channel. Using a relay node, the transmission coverage can be extended with low cost, and a source node can transmit the signal more reliably through increased diversity. These advantages come from the cooperation between the source, the geographically distributed relays, and well developed combining techniques at the destination. Many articles have considered different forwarding schemes at the relay node to use the system resource efficiently [3,4]. They are usually classified by regenerative and non-regenerative methods depending on whether or not the source signal is regenerated at the relay node. The well known decode-andforward (DF) and amplify-and-forward (AF) schemes are examples of the two cases [5]. This general area is called cooperative communication, but some authors call it user cooperation by considering the case where each user can act as a source and a relay [6].
Power allocation methods for cooperative relay channel have been considered in [7,8]. However they all focused on the scenario of three nodes, each of which is equipped with a single antenna. The solutions can be obtained with a relatively straightforward manner. Recently, the scenario was extended to two different cases. One is the scenario of multiple relay nodes with a single antenna. In this case, the optimal power allocation has been derived in [9,10]. The other is the scenario of a single relay with multiple antennas while there are still three nodes in the entire system, and the source and the destination also have multiple antenna. The optimal power allocation for the latter case was proposed in [11,12], which minimizes mean square error (MSE). Pairwise error probability (PEP) using distributed space-time codes was investigated in [13]. Overall channel capacity was considered in [14,15]. They all compute the singular value decomposition of the first hop and the second hop channels, and distribute the total power among the eigenmodes to optimize the performance criteria, MSE or channel capacity.
Recently the most general extension, i.e., a multipleantenna source and a multiple-antenna destination aided by multiple relays with multiple antennas, is considered in [16,17]. They computed the precoding matrix which is used for forwarding in the second hop. In [16], even though the proposed approach has a closed form for a diagonal matrix denoted by F, the final relay precoding matrix should be chosen in order to be a block diagonal matrix with minimum norm. That is, the solution does not have a closed form. Several relaying schemes for each relay are proposed in [17], but they are not capacity-optimal solutions. In [18,19], the authors considered a similar system model with relaying through distributed spacetime codes, and showed that the power allocation solution results in the best relay selection. In fact, it is difficult to obtain a closed form solution for this general case with the block diagonal constraint. As a simpler case, if multiple relays with a single antenna are distributed geographically in a network, the matrix formed by all relay gains should be a diagonal matrix when full cooperation is not possible. However, few related results exist even for this simple case. Recently, [20] considered this simple system model, and derive an asymptotical capacity scaling law. However, they did not deal with the relay power allocation problem as in this article. Thus we focus on this simplified system where there are multiple relays with a single antenna while the source and the destination have multiple antennas. This may also be a practical scenario in the sense that the computational complexity of the relay nodes needs to be minimized. In order to focus on the effect of relay power gain, and compare with the performance of conventional non-regenerative relaying schemes, we only consider the AF scheme and real valued power gain a in this article. For the aforementioned scenario, we propose a relay selection based approach as an approximation to solve for the power allocation problem. It will be shown by simulations that the performance of the proposed method is close to that of the exhaustive (optimal) search approach. We then apply an incremental greedy search algorithm to this proposed method in order to further reduce the computational complexity of the selection based approach. This article deals only with the AF scheme, but a similar approach can be applied to the DF scheme with a modified object function instead of channel capacity.
The remainder of the article is organized as follows. In Section 2, we describe the system model and formulate the problem. We consider the existing inter-relay power allocation in Section 3, and propose the relay selection based approach for capacity maximization in Section 4. In Section 5, we discuss the feedback overhead and the computational complexity. Simulation results are given in Section 6 to assess the performance of the proposed method. Finally, we provide the conclusion in Section 7.
Notation: Boldface lowercase and uppercase letters represent vectors and matrices, respectively. (·)* denotes the complex conjugate transpose of a matrix. |·| means the determinant or the Euclidean (l 2 -) norm depending on the operand (a matrix or a vector).

System model and problem formulation
The source node and the destination node have M transmit antennas and N receive antennas, respectively. In the network, there are R relays with a single antenna. To satisfy the rank requirement of a multiplexing system, we assume M ≤ N and M ≤ R. However, we do not constrain the relationship between R and N. First, the source broadcasts M symbols to all relays in the first hop, and each relay multiplies the received noisy signal with a power gain. After that, all relays transmit these signals to the destination simultaneously in the second hop assuming perfect synchronization among relays. As depicted in Figure 1, the equivalent discrete time signal model can be summarized as follows, where x is the (M × 1) transmitted signal vector from the source node with covariance matrix R x , and y r (R × 1) and y d (N × 1) denote the received signal at all the relays and the destination. H 1 and H 2 are (R × M) and (N × R) channel coefficient matrices from the source to all the relays, and from all the relays to the destination, which have independent and complex Gaussian distributed elements with zero mean and unit variance, i.e., Rayleigh fading. We assume that the source-relay and the relay-destination channels have an identical Rayleigh distribution, which means that all the relays are equidistant from the source and the destination. n 1 and n 2 are additive complex Gaussian noise with covariance matrices Rn 1 and Rn 2 , and G is a diagonal matrix made by relay power gains g = (g 1, . . . , g R ). H = H 2 GH 1 is the equivalent channel matrix, n = H 2 Gn 1 + n 2 is the compounded noise vector with covariance R n = H 2 GR n 1 G * H * 2 + R n 2 . The general capacity formula or the mutual information can be written as in [21] where the factor 1/2 is multiplied because we use two hops to transmit the signal. Our goal is to maximize this value with the individual power constraint allocated on each relay, is the maximum value of the kth relay power. In addition, there is the average power constraint on the relay output such as tr{R x r } ≤ E r , where R x r = GH 1 R x H * 1 G * + GR n 1 G * and E r is the total power of all the relays.
Problem formulation: In the objective function, the diagonal matrix G depends on the first hop and the second hop channels, so it is difficult to obtain a closed form solution.

Inter-relay power allocation
Since the problem may have no closed form solution, we try to find a solution for optimal power allocation by numerical search. Power allocation is performed at the destination node, and the power gain values are fed back to all the relays. We assume that the destination node knows all the channel state information. We also Since

Uniform
A simple method is to use a uniform gain across all the relays. If the values of g k are a constant (b), we can set G = bI where b is chosen to satisfy the constraints. In this case, the capacity can reach the maximum when the equality of the constraint holds. From the average power constraint, where h T 1k means the kth row vector of H 1 . From the individual power constraint, this value should be less than or equal to (g k ) max for each k. Therefore, where we use the notation [x] u l = max (l, min (x, u)) . If the individual constraint is sufficiently large, all g k 's will be the same as b.

Conventional power allocation
In [4], the authors proposed a so called conventional power gain. It is given by the following equation It should be noted that the conventional power allocation depends only on its own first hop channel while the uniform gain allocation depends on the first hop channels of all the relays.

Proposed relay selection based approach
The utility function to be maximized is the determinant of K = I M + R x H * R −1 n H. substituting the equivalent channel and the compounded noise covariance, it is given by where g i is the ith diagonal element of G, h 2i is the ith column of H 2 , and h T 1i is the ith row of H 1 . It was observed that the optimal solution for G may never use certain relays in a poor channel condition, which motivates us to develop a relay selection based approach. In the selection based approach, the solution can be found by where r is a subset of relays, and g r is a vector with zeros at the position of r C . For example, if R = 5 and r = [1 2 4], g r = [g 1 , g 2 , 0, g 4 , 0]. To reduce the search complexity, the optimization over g r can be skipped. A sub-optimal solution can be obtained by finding a solution for the lower bound of the original problem as where g Unif. r denotes the gain vector for uniform gain allocation. The solution is obtained as if only the relays indicated by r participates in the relaying step. In (10), we have the inequality because the optimization over g r is skipped. The approach in (10) turns out to be a relay selection technique with uniform gain allocation among only selected relays to fully satisfy the original constraint. We chose to use the uniform gain allocation instead of the conventional power allocation because the former has better performance than the latter in the terms of capacity, which will be discussed in the simulation section. As will be shown in the simulation results, even with the sub-optimal approach which maximizes the lower bound, the proposed method can achieve performance similar to the optimal solution obtained by exhaustive search. The proposed approach obviously has lower complexity than the exhaustive search.

Feedback requirement
Let us consider the required amount of feedback bits from the destination to the relays. In the optimal solution by the exhaustive search, the power gain values need to be fed back to the corresponding relays. To use uniform gain allocation, we need to feed back the common power gain b to the relays, which can be broadcasted by the destination. The conventional power allocation depends only on its own first hop channel from the source to each relay which is assumed available at the relay. Thus relay power can be calculated in each relay. Assuming that the total power E r and the number of relays R are known at the relays, there is no need for feedback in this scheme.
The relay selection based approach needs to feed back the selection indicator bits as well as the power gain values of the selected relays. As for the selection indicator bits, one bit (0 or 1) for each relay can be used to indicate whether a given relay is selected or not. Note that the feedback overhead of the proposed method will be significantly less than that of the optimal exhaustive search. When uniform gain allocation is used with the relay selection, only a single value (the common power gain) can be broadcasted by the destination. If conventional power allocation is used with the relay selection, the number of selected relays need to be broadcasted by the destination. As mentioned before, we employ uniform gain allocation in this article because it has slightly better performance than the conventional power allocation when no relay selection is used, which will be shown in the simulation results.

Computational complexity
The full search solution of proposed method gives us information on which relays are selected. However, the full search needs the computation of the utility function for every possible relay subset, and comparison of these values. When the number of relays R is large, the number of possible relay subsets is which increases exponentially as R increases, O(2 R ). Thus the complexity of the selection process itself may become infeasible. In order to reduce this selection complexity, we use an incremental greedy search algorithm. This algorithm first selects the relay with maximum utility function assuming that there is only one relay in the network. At the second stage, the algorithm searches for another relay (among the remaining R-1 relays) which maximizes the utility function when combined with previously selected relay. These operations are repeated until the utility function decreases or all the relays are selected. This greedy algorithm is summarized as in Algorithm 1. We use the notations where END From simulations, it was observed that the greedy algorithm usually stops at a value much smaller than R. Even in the worst case where the iteration continues until m = R, this only requires R(R+1) 2 computations and comparisons which have quadratic growth rate, O(R 2 ). Therefore, if R ≥ 3, the greedy algorithm becomes advantageous over the full search in terms of computational complexity. This gain will be huge when R is large. For example, when R = 8 which is used in our simulations, the full search considers 255 combinations while the greedy algorithm considers 36 combinations even in the worst case.
Let us compare the overall complexity roughly. We denote the complexity of the optimal numerical search by F opt , and that of our selection based method by F sel . If the grid size (interval) used by optimal numerical search is Δ, the total number of grid points is On the other hand, the selection process even for full search has the relay combination number which is pro- Thus it is obvious that we have F opt F sel even if we use the well-developed optimization tools for optimal solution. This complexity gap will be larger when the grid is sufficiently dense and/or R is large. Furthermore, if we combine the selection process with a greedy algorithm, the complexity reduction can be much greater.

Simulation results
In this section, we compare the proposed power allocation methods with the optimal (ex-haustive) method and two existing power allocation methods. In all simulations, we use E s = E r , and the noise variance is σ 2 R = σ 2 D = 1. 'Exhaustive' means the power gains optimally found by a numerical method. 'Uniform' and 'Conventional' use the power allocation methods described in the section of inter-relay power allocation. The above three methods use all the relays since relay selection is not used. 'FullSearch' and 'Greedy' are the proposed relay selection based approaches. Figures 2 and 3 shows the performance of various power allocation methods in terms of average capacity. Figure 2 is in linear scale for the y-axis, and a log scale is used for the y-axis in Figure 3 to see the difference clearly in the low SNR regime, It is observed that the selection based methods perform similarly to the optimal (exhaustive) method. When the solution for power allocation problem can not be easily derived in a closed form solution, high SNR assumption or an asymptotical approach is usually used in the literature. These tend to deviate from the optimal solution especially in the low SNR regime. However, it is interesting to note that our proposed selection techniques which are 'FullSearch' and 'Greedy' have comparable performance in all SNR regimes.
Let us consider the performance difference between the uniform gain allocation g Unif. and the conventional power allocation G Conv. . To simplify the analysis, we focus on g 1 . The gain factors are then If it is assumed that h T 11 2 is the largest among all rows, that means the first hop channel to the first relay is very good. In this case, it is intuitively better to allocate more power for this relay. However, the conventional case allocates less power than the uniform case since R · |h 11 | 2 ≥ i |h 1i | 2 . Therefore, the uniform method performs slightly better than the conventional method. This is why we use the uniform gain allocation instead of the conventional power allocation in the proposed relay selection based approach. In order to show the results for more general antenna configurations, we show the performance comparison with respect to the number of relays in Figure 4, where we use E r = 20 dB to show the performance difference more clearly. We fix M = 3 and N = 6. The proposed methods have similar performance to the optimal (exhaustive) case for all R values. The performance of the two proposed methods begin to deviate from the exhaustive method slightly as the number of relays gets  larger. Note that the deviation is less than 0.1 bps/Hz at R = 8. But the reduction of feedback overhead and computational complexity will be more significant as R gets larger, so the proposed methods may still be attractive in practice. Figure 4 also shows that if we remove the relays in bad condition and allocate more power to the other remaining relays, the performance can be improved when uniform gain allocation is used. It is also observed that the incremental greedy search algorithm has negligible performance loss compared to the full search case.

Conclusions
In this article, we proposed a relay selection based approach for inter-relay power allocation problem. Since it is difficult to obtain a closed form solution for the problem, we propose a low-complexity technique based on relay selection using uniform gain allocation, which has near-optimal performance in terms of capacity. We also combine with the proposed method with an incremental greedy search algorithm in order to further reduce the search complexity for the relay selection. Simulations show that the proposed methods for relay power allocation have performance close to the optimal (exhaustive) power allocation in terms of average capacity. The proposed methods appear to be promising in terms of computational complexity and capacity performance.
Endnote a If we consider the complex valued gain, the performance can be improved by adjusting the phase with extra processing. When a complex number is treated as two real values (magnitude and phase), the feedback overhead (as will be discussed) may overshadow the benefit of power allocation. Thus we focus on the real valued power allocation problem in this article.