Collaborative Relay Beamforming Strategies for Multiple Destinations with Guaranteed QoS in Wireless Machine-to-Machine Networks

Abstract

The Machine-to-Machine (M2M) communications allow information exchange between machine devices, which can be carried out without any human interaction. However, there are a large number of small and low-power machine devices in the wireless M2M networks. To guarantee the quality of service (QoS) requirements of the destination devices, we study the amplify-and-forward (AF) relay beamforming, where multiple relay M2M devices can transmit signals from the source M2M device to multiple destination M2M devices. In this paper, we propose two iterative strategies to jointly optimize the source antenna selection and the collaborative relay beamforming weights with the aid of perfect channel state information (CSI). The aim of the proposed strategies is to maximize the worst-case received signal-to-interference-and-noise ratio (SINR) under two different types of relay power constraints, which are the total relay power constraint and individual relay power constraints, respectively. Using the semidefinite relaxation (SDR) technique, the optimization problem of collaborative relay beamforming can be formulated as a semidefinite programming (SDP) problem, which can be optimally solved. Simulation results validate our theoretical analysis and demonstrate that after several iterations, the performance of the proposed iterative strategies can obtain near-optimal performance.

1. Introduction

Machine devices, such as sensors, communicate with each other in machine-to-machine (M2M) networks, which do not need intervention from human beings [1]. A wireless M2M device may act as a source point and communicate data packets on behalf of the other M2M destination devices [2]. However, the wireless M2M networks are composed of a large number of small and power-limited machine devices, which can only transmit low-power signals. Due to the propagation loss between the source devices and the destination devices, it is very challenging to guarantee the quality of service (QoS) requirements of the destination devices. Recently, many increasing researches showed interest in the relay devices to deal with this problem [3]. The sources transmit signals through one or more relay devices towards the destination devices in wireless M2M networks. Relay devices can be designed to work collaboratively to guarantee the QoS requirements by forming a virtual beamforming system [4, 5].

In general, there are three types of relay schemes, including the amplify-and-forward (AF) [6], decode-and-forward (DF) [7], and compress-and-forward (CF) [8] schemes. In the AF scheme, relay nodes receive the signals transmitted from the source and forward the amplified signals to the destinations. In the DF scheme, relay nodes decode and reencod the received signals to the destinations. In the CF scheme, relay nodes send quantized and compressed signals to the destinations. Among them, the AF scheme is the most attractive for the collaborative relay beamforming due to its low implementation complexity. With the aid of channel state information (CSI), the relay nodes can work collaboratively in a way similar to a multiple antenna systems to construct virtual beams towards the destinations [9].

In this paper, we study the collaborative relay beamforming for multiple destinations in the multiple relay M2M networks, where multiple relay M2M devices form virtual beams to transmit signals from a source M2M device with multiple antennas towards multiple destination M2M devices. Firstly, we present two methods to select the source antenna to transmit the destination's signals. In Method $1$ , we use the greedy algorithm to select the best available channel for each destination until all the destinations are allocated properly with the source antennas; in Method $2$ , we balance the channel condition for each destination and its interference to the other destinations based on Method $1$ . Secondly, we optimize the relay beamforming weights with the aid of perfect CSI to maximize the worst-case received SINR under two different types of relay power constraints, which are the total relay power constraint and individual relay power constraints, respectively. Using the semidefinite relaxation (SDR) technique [10], we show that the optimization problem can be formulated as a semidefinite programming (SDP) problem, which can be optimally solved using interior-point methods [11]. Finally, according to the source antenna selection methods and relay beamforming schemes, we propose two iterative strategies to jointly optimize the source antenna selection and the relay beamforming weights. We show that after several iterations, the performance of the proposed iterative strategies can obtain near optimal performance.

According to our proposed iterative strategies, the source M2M device must know all the CSI of the wireless M2M networks to jointly optimize the source antenna selection and the relay beamforming weights. However in practice, the feedback channels are always limited and could be overloaded if the number of relays and destinations is large. To deal with this practical problem, Choi [12, 13] studied the distributed AF relay beamforming schemes for the single destination where the relays can work cooperatively to form a virtual beam with local CSI (The local CSI means that each relay node knows its instantaneous incoming and outgoing channel coefficients (from relay nodes’ point of view). However, the CSI of the other relay nodes is assumed to be not known.) when a consensus algorithm was employed in [12], and different MMSE criteria were proposed in [13]. D. H. N. Nguyen and H. H. Nguyen [14] proposed two relay power allocation schemes in a distributed manner to minimize the sum relay power with guaranteed SNR at the destinations and maximize the SNR margin subject to individual relay power constraints, respectively, under the orthogonal transmissions of each source-destination pair. In the end of this paper, we discuss the distributed multirelay beamforming strategy to maximize the worst-case received SINR for multiple destinations, where each relay may independently calculate its beamforming weight with the local CSI.

The rest of this paper is organized as follows. In Section 2, we present a brief discussion of related works on relay beamforming which motivated our work. Section 3 describes the multiple relay M2M networks with multiple destinations. Section 4 presents two methods to select the source antenna for each destination. In Section 5, we optimize the multi-relay weights that achieve the maximum worst-case SINR under two different types of relay power constraints. Then, two iterative strategies to jointly select the source antennas and optimize the multirelay weights are proposed, respectively. Section 6 presents and analyzes the simulation results. The distributed multi-relay beamforming strategy is discussed in Section 7. Conclusions are drawn in Section 8.

Notation.

Vectors are written in boldface lowercase letters, for example, x, while matrices are denoted by boldface uppercase letters for example, X. The superscripts *, T, and H stand for the conjugate, transposition, and hermitian of a complex vector or matrix, respectively. $Tr (S)$ represents the trace of a matrix S, and $diag (s_{1}, s_{2}, \dots, s_{k})$ denotes a diagonal square matrix with $s_{1}, s_{2}, \dots, s_{k}$ as the diagonal elements. $𝔼 [z]$ calculates the expectation of a random entity z, while $| z |$ takes the modulus of a complex number z. The real and complex number fields are denoted by ℝ and ℂ, respectively. In addition, $A ≽ 0$ means that A is a hermitian semidefinite matrix, and O means the upper bound of the computational complexity.

2. Related Works

Andreev et al. [2] studied a typical smart metering M2M application and proposed a client relay scheme to improve the performance of M2M devices with poor communication link. Khajehnouri and Sayed [3] designed relay strategies in sensor networks with minimum mean square error (MMSE) performance subject to local and global relay power constraints, respectively. Elkheir et al. [15] presented an M2M communication scheme, where the source communicates with the corresponding destination node with the help of one or more relay nodes. The objective in [15] was to minimize the total relay power under specific performance constraints and the individual power of each relay. Zheng et al. [16, 17] proposed the relay beamforming schemes that achieve the maximum received signal-to-noise ratio (SNR) for a single destination under both total and individual power constraints with the aid of perfect CSI [16] and imperfect CSI [17]. In addition, distributed beamforming for AF relay nodes was investigated when a consensus algorithm was employed for cooperative beamforming in [12], and the minimum mean square error (MMSE) criteria were used in [13].

However, existing approaches focused on relay strategies with the single or multiple relay beamforming for only one destination in wireless M2M networks or sensor networks. The relay M2M devices can be used to support multiple destinations for the M2M communications. Nguyen et al. [18, 19] studied the optimal relay beamforming schemes in a multiple AF relay cooperating network with multiple source-destination pairs. The aim in [18] was to minimize the total relay power under the signal-to-interference-and-noise ratio (SINR) requirement at each destination. In [19], Nguyen et al. minimized the sum relay power with the guaranteed QoS in terms of SNR at the destinations in orthogonal channels. Liu and Petropulu [20] studied the relay beamforming to minimize the sum transmit power at the relays with multiple pairs of the source and destination, assisted by multiple relays. Zheng and Blostein [21] considered the collaborative relay beamforming system to minimize the sum power at the relays and the transmit power at the source, where a single source with multiple antennas transmits signals to multiple users through multiple single-antenna relays. Phan et al. [22] studied the joint power allocation and AF relay selection problem to maximize the destination rates and minimize the total transmit power at the relays, respectively, for multiple destinations in multiple orthogonal source-destination pair networks.

The main difference between our work and the works in [18–22] is that we present two methods to select the source antenna and optimize the relay beamforming weights to maximize the worst-case received SINR under two different types of relay power constraints, which are the total relay power constraint and individual relay power constraints, respectively, while the proposals in [18–21] were to minimize the total relay transmit power under the QoS of destinations, and Phan et al. [22] maximized the destination rates and minimized the total transmit power at the relays, respectively.

3. System Model

We consider a multiple relay M2M network consists of a source with L antennas, M relay nodes, and K destinations each with a single antenna as shown in Figure 1, where $K \leq L$ . (If $K > L$ , the source can select destinations to transmit the signals to exploit selection diversity among destinations [23].) Since we consider the small and power-limited M2M devices, there is no direct link between the source and destinations. The AF relay nodes form virtual beams to transmit the source signals towards multiple destinations on the same frequency band. Each destination may receive interference signals of other destinations from relay nodes. We assume that the coefficients of the channels from the source to the mth relay and the mth relay to the destination k denote $g_{l, m} \in ℂ$ and $h_{k, m} \in ℂ$ , respectively. We consider a flat fading channel model, where all the channel coefficients, including $g_{m} = [g_{1, m}, g_{2, m}, \dots, g_{L, m}]$ , $m = 1,2, \dots, M$ and $h_{k} = [h_{k, 1}, h_{k, 2}, \dots, h_{k, M}]$ , $k = 1,2, \dots, K$ are independent in our multiple relay M2M networks. For convenience, the propagation losses are also contained in $g_{m}$ and $h_{k}$ to reflect the relative locations of the relays and the destinations.

Figure 1

The multiple relay M2M networks with multiple destinations.

It is assumed that the AF relay protocol has two steps. During the first step, the source selects the best available antenna to broadcast the kth destination's signal $s_{k}$ to the relay nodes. Two methods to select the best antenna for the destination k will be presented in Section 4. We assume that the destinations’ signals are independent with $𝔼 [s_{k}] = 0$ and $𝔼 [| s_{k} |^{2}] = P_{k}$ . The received signal at the relay node m is given by

\begin{matrix} r_{m} = g_{m} \sum_{k = 1}^{K} ‍ d_{k} s_{k} + w_{m}, \end{matrix}

(1)

where $w_{m}$ is independent white Gaussian noise at relay node m with zero mean and variance of $N_{m}$ . The column vector $d_{k}$ is the vector that represents the source antenna selection for the kth destination.

Then, during the second step, the relay nodes forward the received signals to the destinations using complex weights $v_{m}$ , $m = 1,2, \dots, M$ for all destinations for the collaborative relay beamforming. Thus, the transmitted signal vector from the mth relay node is

\begin{matrix} x_{m} = v_{m} r_{m} . \end{matrix}

(2)

The destination k receives the additive signals from all the relay nodes, which is given by

\begin{array}{l} y_{k} = \sum_{m = 1}^{M} ‍ h_{k, m} v_{m} r_{m} + n_{k} \\ = \sum_{m = 1}^{M} ‍ h_{k, m} v_{m} g_{m} d_{k} s_{k} + \sum_{m = 1}^{M} ‍ h_{k, m} v_{m} g_{m} \sum_{j = 1, j \neq k}^{K} ‍ d_{j} s_{j} \\ + \sum_{m = 1}^{M} ‍ h_{k, m} v_{m} w_{m} + n_{k}, \end{array}

(3)

where $n_{k}$ is independent white Gaussian noise at the kth destination with zero mean and variance of $N_{k}$ . Thus, the received SINR of the destination k is given by

\begin{array}{l} SIN R_{k} = \frac{P_{k} {| \sum_{m = 1}^{M} h_{k, m} v_{m} g_{m} d_{k} |}^{2}}{\sum_{j = 1, j \neq k}^{K} P_{j} {| \sum_{m = 1}^{M} h_{k, m} v_{m} g_{m} d_{j} |}^{2} + P_{n}} \\ = \frac{P_{k} {| b_{k}^{T} d_{k} |}^{2}}{\sum_{j = 1, j \neq k}^{K} P_{j} {| b_{k}^{T} d_{j} |}^{2} + P_{n}}, \end{array}

(4)

where

\begin{matrix} P_{n} = \sum_{m = 1}^{M} ‍ {| h_{k, m} |}^{2} {| v_{m} |}^{2} N_{m} + N_{k}, \\ b_{k} = {[\sum_{m = 1}^{M} ‍ h_{k, m} v_{m} g_{1, m}, \dots, \sum_{m = 1}^{M} ‍ h_{k, m} v_{m} g_{L, m}]}^{T} . \end{matrix}

(5)

4. Source Antenna Selection for Each Destination

In this section, we present two methods to optimize the source-antenna selecting vector $d_{k}$ for the kth destination with fixed relay weights: (1) a greedy algorithm, that is, select the best available channel for each destination until all the destinations are allocated properly with the source antennas; (2) balance the channel condition for each destination k and its interference to the other destinations.

In Method $1$ , we first choose the element with the maximum magnitude denoted as $G_{l_{k}, k}$ from the $L \times K$ composite channel matrix $G = [b_{1}, b_{2}, \dots, b_{K}]$ . Then the antenna $l_{k}$ is selected to transmit signals to the destination k. The channel coefficients between the rest antennas and destinations form a new $(L - 1) \times (K - 1)$ channel matrix. We can repeat the above processes K times until each destination is allocated with one antenna.

In Method $2$ , when we select the $l_{k}$ th antenna to transmit signals towards the kth destination using Method $1$ , we consider its interference to the other destinations, that is, $\sum_{m = 1}^{M} ‍ h_{j, m} v_{m} g_{l_{k}, m}$ , $j \neq k$ , which are the $l_{k}$ th row elements of the composite channel matrix G. If all the magnitudes of the interference are less than or equal to the interference threshold ξ, we keep the channel and allocate the antenna for the corresponding destination on this channel. When one of the magnitudes of the interference is larger than the threshold ξ, we drop the selected channel and choose the second largest magnitude element from the composite channel matrix. We repeat the above processes until all the magnitudes of the interference are less than or equal to the threshold ξ. If the feasible solution can not be obtained under the interference threshold ξ, then we may double the threshold value and repeat the above processes to select another antenna for the destination.

In our methods, the computational complexity of selecting one antenna scales linearly with the dimension of channel matrix G, that is, $K L$ . Since we need to select K antennas for K destinations from the channel matrixes whose dimensions decrease step by step, the computational complexity of our methods is at most $O (K! (L! / (L - K)!))$ .

5. Iterative Strategies for Max-Min SINR with Relay Power Constraints

In this section, we assume that the source-antenna selection vector $d_{k}$ for each destination k has been already chosen by the methods in Section 4. Thus, we assume that the $l_{k}$ th antenna transmits $s_{k}$ to the kth destination. We can rewrite the expression of the received SINR at the destination k in (4) as

\begin{array}{l} SIN R_{k} = \frac{P_{k} {| \sum_{m = 1}^{M} h_{k, m} g_{l_{k}, m} v_{m} |}^{2}}{\sum_{j = 1, j \neq k}^{K} P_{j} {| \sum_{m = 1}^{M} h_{k, m} g_{l_{j}, m} v_{m} |}^{2} + P_{n}} \\ = \frac{P_{k} {| f_{l_{k}}^{T} v |}^{2}}{\sum_{j = 1, j \neq k}^{K} P_{j} {| f_{l_{j}}^{T} v |}^{2} + v^{H} H_{k} v + N_{k}}, \end{array}

(6)

where

\begin{matrix} v = {[v_{1}, v_{2}, \dots, v_{M}]}^{T}, \\ f_{l_{k}} = {[h_{k, 1} g_{l_{k}, 1}, h_{k, 2} g_{l_{k}, 2}, \dots, h_{k, M} g_{l_{k}, M}]}^{T}, \\ H_{k} = diag ({| h_{k, 1} |}^{2} N_{1}, {| h_{k, 2} |}^{2} N_{2}, \dots, {| h_{k, M} |}^{2} N_{M}) . \end{matrix}

(7)

In the following, we optimize the multi-relay beamforming weights to maximize the worst-case received SINR for multiple destinations and propose the iterative strategies under two different types of relay power constrained in the multiple relay M2M network. We first investigate the case where the total relay power is limited, and then study the scenario where the individual relay transmit power is constrained.

5.1. Total Relay Power Constraint

In this subsection, we optimize the multi-relay beamforming weights for multiple destinations under a total relay power constraint $P_{total}$ . Then we present the iterative strategy to jointly optimize the source-antenna selection and the relay beamforming weights under the total relay power constraint in the multiple relay M2M networks.

For the optimum, we can formulate the problem as

\begin{array}{l} \max_{v} \min_{k} SIN R_{k} \\ s.t. P_{t} \leq P_{total}, \end{array}

(8)

where $P_{t}$ is the total transmission power of relay nodes, it can be presented as

\begin{array}{l} P_{t} = \sum_{m = 1}^{M} ‍ {| v_{m} |}^{2} E [{| r_{m} |}^{2}] \\ = \sum_{m = 1}^{M} ‍ {| v_{m} |}^{2} (\sum_{k = 1}^{K} ‍ {| g_{l_{k}, m} |}^{2} P_{k} + N_{m}) \\ = v^{H} E v, \end{array}

(9)

\begin{array}{l} E = diag (\sum_{k = 1}^{K} ‍ {| g_{l_{k}, 1} |}^{2} P_{k} + N_{1}, \dots, \sum_{k = 1}^{K} ‍ {| g_{l_{k}, M} |}^{2} P_{k} + N_{M}) . \end{array}

(10)

In order to guarantee the QoS for multiple destinations, we consider the maximization of the worst SINR among all destinations in (8). However, the object function of the optimization problem (8) is quasi-convex. Directly solving Problem (8) is usually difficult due to its nonconvexity. We could adopt the bisection search algorithm to efficiently find the optimal solution of Problem (8). Specifically, we could solve a sequence of the following feasibility problems for each given $Γ$ [11].

\begin{array}{l} \min_{v} P_{t} \\ s.t. {\begin{cases} P_{t} \leq P_{total}, \\ SIN R_{k} \geq Γ, & \begin{array}{l} k = 1,2, \dots, K . \end{array} \end{cases} \end{array}

(11)

Denote the optimal solution of (8) by $γ_{opt}$ . If $Γ = γ_{opt}$ in (11), then the optimal solution v of (11) is also optimal for (8). Thus we can repeatedly search for the optimal $Γ$ using the bisection search method [11]. In the following, we focus on solving the problem in (11) for a given $Γ$ .

With the received SINR in (6) and the total relay transmit power in (9), we can express Problem (11) as

\begin{array}{l} \min_{v} v^{H} E v \\ s.t. {\begin{cases} v^{H} E v \leq P_{total} \\ \frac{P_{k}}{Γ} {| f_{l_{k}}^{T} v |}^{2} - \sum_{j = 1, j \neq k}^{K} ‍ P_{j} {| f_{l_{j}}^{T} v |}^{2} \\ - v^{H} H_{k} v \geq N_{k}, & k = 1,2, \dots, K . \end{cases} \end{array}

(12)

Using the definition $V ≜ v v^{H}$ , we can reformulate Problem (12) as

\begin{array}{l} \min_{V} Tr (E V) \\ s.t. {\begin{cases} V ⪰ 0, \\ rank (V) = 1, \\ Tr (E V) \leq P_{total}, \\ \frac{1}{Γ} Tr (F_{k} V) - \sum_{j = 1, j \neq k}^{K} ‍ Tr (F_{j} V) \\ - Tr (H_{k} V) \geq N_{k}, & k = 1,2, \dots, K, \end{cases} \end{array}

(13)

where $F_{k} = P_{k} (f_{l_{k}}^{*} f_{l_{k}}^{T})$ and $V ≽ 0$ means that V is a hermitian semidefinite matrix. Note that the rank constraint $rank (V) = 1$ comes from $V ≜ v v^{H}$ , which is not convex. Thus, the optimization problem (13) is not convex. We can use the SDR technique to relax Problem (13). Here, we drop the rank-one constraint and obtain the following problem.

\begin{array}{l} \min_{V} Tr (E V) \\ s.t. {\begin{cases} V ⪰ 0 \\ Tr (E V) \leq P_{total}, \\ \frac{1}{Γ} Tr (F_{k} V) - \sum_{j = 1, j \neq k}^{K} ‍ Tr (F_{j} V) \\ - Tr (H_{k} V) \geq N_{k}, & k = 1,2, \dots, K . \end{cases} \end{array}

(14)

Note that Problem (14) becomes an SDP problem, which can be optimally solved using interior-point methods, whose computational complexity is polynomial [11]. As shown in [24], for a given solution accuracy ϵ, the interior-point methods need at most $O (β (α^{3} + β α^{2} + β^{2}) \sqrt{α} \log (1 / ϵ))$ arithmetic operations, where α is the dimension of the matrix variable, and β is the number of constraints. In our SDP problem in (14), the dimension of the matrix variable is $α = M$ , and the number of constraints is $β = K$ . Therefore, the computational complexity of Problem (14) is $O (K^{3} M^{3.5} \log (1 / ϵ))$ .

The solution of the SDP problem in (14) provides a lower bound of the objective function in the original problem in (13) due to excluding the rank-one constraint. However, in our extensive simulations, we have never observed that the optimal solution of Problem (14) has a rank higher than one. The similar observation was also reported in [18] and [25] to design the optimal beamforming schemes. Thus, the original problem in (13) can be optimally solved. For the cases where Problem (14) has a rank higher than one, several randomization techniques [26] can be used to provide a good approximation solution $\tilde{v}$ from the optimal solution V of Problem (14).

Based on the above results, we propose an iterative strategy (see, Algorithm 1) to jointly select the source antennas and optimize the multi-relay beamforming weights with the total relay power constraint for multiple destinations.

Algorithm 1: Iterative Strategy under the total relay power constraint.

1: Initialize $γ_{L B}$ , $γ_{U B}$ and $P_{total}$ . For some given ε > 0

2: while $γ_{U B}$ – $γ_{L B}$ > εdo

3: Initialize the multirelay beamforming weights $v_{m}$ =

$c_{m} e^{j θ_{m}}$ , $m = 1,2, \dots, M$ , where $c_{m} = P_{total} / M$ . $θ_{m}$

is a uniform random variable chosen from the interval

[0, 2π].

4: Apply Method 1 or Method 2 in Section 4 to select

the source antenna for each destination with the current

fixed relay weights.

5: Optimally solve Problem (14) with $Γ = (γ_{L B} + γ_{U B}) / 2$ .

6: If Problem (14) is feasible, then obtain the optimal

solution $V_{opt}$ .

7: if rank( $V_{opt}$ ) = 1 then

8: obtain the multirelay beamforming weights vector v,

which is the optimal solution of the original problem

(11).

9: else

10: obtain the multirelay beamforming weights vector v

using the randomization techniques.

11: end if

12: If the total relay power is sufficiently close to a fixed

point, or the iteration number exceeds a predetermined

number, then stop the iteration and go to Step 13.

Otherwise go back to Step 4.

13: if Problem (14) is feasible then

14: $γ_{L B} = Γ$

15: else

16: $γ_{U B} = Γ$

17: end if

18: end while

19: return the maximization of the worst-case SINR $Γ$ and

the multirelay beamforming weights vector v.

According to [11], the bisection search method requires $⌈ \log_{2} ((γ_{U B} - γ_{L B}) / ε) ⌉$ iterations. Thus, the complexity of Algorithm 1 is $O (K! (L! / (L - K)!) K^{3} M^{3.5} \log (1 / ϵ) ⌈ \log_{2} ((γ_{U B} - γ_{L B}) / ε) ⌉)$ for each iteration. Since the best available channel is selected for each destination in Step $4$ , and the worst-case SINR is maximized from Step $5$ to Step $18$ , after several iterations, Algorithm 1 will converge to a fixed point. We can observe this convergence in the simulations.

5.2. Individual Relay Power Constraints

In this subsection, we study the scenario where each relay power is constrained by $P_{m}^{i}$ , $m = 1,2, \dots, M$ . Note that the individual relay power constraints are important when the battery lifetime of the relay nodes is limited. Then we present the iterative strategy under individual relay power constraints in the multiple relay M2M networks.

We can formulate the problem as

\begin{array}{l} \max_{v} \min_{k} SIN R_{k} \\ s.t. P_{m} \leq P_{m}^{i}, m = 1,2, \dots, M, \end{array}

(15)

where $P_{m}$ is the transmission power of relay node m, it can be given by

\begin{matrix} P_{m} = {| v_{m} |}^{2} E [{| r_{m} |}^{2}] = {| v_{m} |}^{2} (\sum_{k = 1}^{K} ‍ {| g_{l_{k}, m} |}^{2} P_{k} + N_{m}) . \end{matrix}

(16)

Since the optimization problem in (15) is quasi-convex, we can also obtain the optimal solution of (15) by repeatedly solving the following problem using the bisection search method to search for the $Γ$ [11].

\begin{array}{l} \min_{v} P_{t} \\ s.t. {\begin{cases} P_{m} \leq P_{m}^{i}, & m = 1,2, \dots, M, \\ SIN R_{k} \geq Γ, & k = 1,2, \dots, K . \end{cases} \end{array}

(17)

With the received SINR in (6) and the mth relay transmit power in (16), we can express Problem (17) as

\begin{array}{l} \min_{v} v^{H} E v \\ s.t. {\begin{cases} {| v_{m} |}^{2} \\ \times (\sum_{k = 1}^{K} ‍ {| g_{l_{k}, m} |}^{2} P_{k} + N_{m}) \leq P_{m}^{i}, \\ \frac{P_{k}}{Γ} {| f_{l_{k}}^{T} v |}^{2} - \sum_{j = 1, j \neq k}^{K} ‍ P_{j} {| f_{l_{j}}^{T} v |}^{2} \\ - v^{H} H_{k} v \geq N_{k}, & m = 1,2, \dots, M, \\ \begin{array}{l} k = 1,2, \dots, K . \end{array} \end{cases} \end{array}

(18)

Using the definition $V ≜ v v^{H}$ , $F_{k} = P_{k} (f_{l_{k}}^{*} f_{l_{k}}^{T})$ and $A = (\sum_{k = 1}^{K} ‍ | g_{l_{k}, m} |^{2} P_{k} + N_{m})$ , we can reformulate Problem (18) as

\begin{array}{l} \min_{V} Tr (E V) \\ s.t. {\begin{cases} V ⪰ 0, \\ rank (V) = 1, \\ {[V]}_{m, m} \leq \frac{P_{m}^{i}}{A}, & m = 1,2, \dots, M, \\ \frac{1}{Γ} Tr (F_{k} V) - \sum_{j = 1, j \neq k}^{K} ‍ Tr (F_{j} V) \\ - Tr (H_{k} V) \geq N_{k}, & \begin{array}{l} k = 1,2, \dots, K, \end{array} \end{cases} \end{array}

(19)

where $[V]_{m, m}$ denotes the $(m, m)$ th element of V and $V ≽ 0$ means that V is a hermitian semidefinite matrix. Note that the constraint $rank (V) = 1$ is not convex. As mentioned in Section 5.1, we also use the SDR technique to relax Problem (19) and obtain the following SDP problem.

\begin{array}{l} \min_{V} Tr (E V) \\ s.t. {\begin{cases} V ⪰ 0, \\ {[V]}_{m, m} \leq \frac{P_{m}^{i}}{A}, & m = 1,2, \dots, M, \\ \frac{1}{Γ} Tr (F_{k} V) - \sum_{j = 1, j \neq k}^{K} ‍ Tr (F_{j} V) \\ - Tr (H_{k} V) \geq N_{k}, & \begin{array}{l} k = 1,2, \dots, K . \end{array} \end{cases} \end{array}

(20)

Note that Problem (20) can be optimally solved with a polynomial time complexity [11]. As shown in Section 5.1, the computational complexity of Problem (20) is the same as that of Problem (14), which is $O (K^{3} M^{3.5} \log (1 / ϵ))$ . In our extensive simulations, we have never observed that the optimal solution of SDP relaxation problem in (20) has a rank higher than one. Thus, the original problem in (19) can be optimally solved. For the cases where Problem (20) has a rank higher than one, we can also use randomization techniques [26] to provide a good approximation solution $\tilde{v}$ .

In the following, we propose another iterative strategy to jointly optimize the source-antenna selection and the multi-relay beamforming weights with the individual relay power constraints for multiple destinations.

According to [11], the complexity of Algorithm 2 is $O (K! (L! / (L - K)!) K^{3} M^{3.5} \log (1 / ϵ) ⌈ \log_{2} ((γ_{U B} - γ_{L B}) / ε) ⌉)$ for each iteration. We can observe that Algorithm 2 will converge to a fixed point after several iterations in the simulations.

Algorithm 2: Iterative Strategy under the individual relay power constraints.

1: Initialize $γ_{L B}$ , $γ_{U B}$ and $P_{m}^{i}$ . For some given ε > 0

2: while $γ_{U B} - γ_{L B} > ε$ do

3: Initialize the multirelay beamforming weights $v_{m} =$

$c_{m} e^{j θ_{m}}$ , $m = 1,2, \dots, M$ , where $c_{m} = P_{m}^{i}$ . $θ_{m}$ is

a uniform random variable chosen from the interval

[0, 2π].

4: Apply Method 1 or Method 2 in Section 4 to select

the source antenna for each destination with the current

fixed relay weights.

5: Optimally solve Problem (20) with $Γ = (γ_{L B} + γ_{U B}) / 2$ .

6: If Problem (20) is feasible, then obtain the optimal

solution $V_{opt}$ .

7: As shown in Algorithm 1, obtain the multirelay beam-

forming weights vector v from $V_{opt}$ .

8: If the total relay power is sufficiently close to a fixed

point, or the iteration number exceeds a predetermined

number, then stop the iteration and go to Step 9.

Otherwise go back to Step 4.

9: if Problem (20) is feasible then

10: $γ_{L B} = Γ$

11: else

12: $γ_{U B} = Γ$

13: end if

14: end while

15: return the maximization of the worst-case SINR $Γ$ and

the multirelay beamforming weights vector v.

6. Simulation Results

In this section, we present simulation results to validate the performance of our iteration strategies for multiple destinations in the multiple relay M2M networks. It is assumed that the flat channel coefficients $g_{m} = [g_{1, m}, \dots, g_{L, m}]$ , $m = 1,2, \dots, M$ and $h_{k} = [h_{k, 1}, \dots, h_{k, M}]$ , $k = 1,2, \dots, K$ are independent Rayleigh random variables with zero mean and unit variance. For convenience, it is assumed that the noise at all nodes has the same power density, that is, $N_{m} = N_{k} = N_{0}$ . In the scenario of individual relay power constraints, we assume that each relay has the same power constraint, that is, $P_{m}^{i} = P_{total} / M$ . The source transmission SNR is set to $P_{k} / N_{0} = 10 (dB)$ , $k = 1,2, \dots, K$ , which is the same for all destinations. We use CVX [27] to numerically solve the SDP problems.

In Figure 2, we compare the performance of the proposed iteration strategies for different values of the total relay SNR of $P_{total} / N_{0}$ . The multiple relay M2M network has one source with $8$ antennas, $20$ relays, and $5$ destinations. We run simulations for our proposed strategies using Method $1$ and Method $2$ mentioned in Section 4 to select the antenna for each destination in Figures 2(a) and 2(b), respectively. The results in Figures 2(a) and 2(b) show that the performance increases as the number of iterations increases and after $3$ iterations, the proposed iteration strategies can obtain near-optimal performance. Thus, $3$ iterations could be enough for the proposed iteration strategies in various multiple destination scenarios. The results indicate that the performance of the proposed iteration strategy with the total power constraint (solid lines) achieves 1~2 (dB) higher SINR than that of the iteration strategy with the individual power constraints (dashed lines). In addition, comparing between Figures 2(a) and 2(b), the performance of our proposed strategies using Method $2$ is slightly superior to that of our proposed strategies using Method $1$ , as Method $2$ considers the interference to the other destinations from a selected channel based on Method $1$ .

Figure 2

Expected worst-case received SINR of the proposed iteration strategies with the total relay power constraint (solid lines) and the individual relay power constraints (dashed lines) for different numbers of iterations.

Figure 3 demonstrates the worst-case received SINR results for the proposed iteration strategies with different values of the total relay SNR of $P_{total} / N_{0}$ . In Figure 3, the curves with the square symbols represent that the multiple relay M2M network has $8$ antennas, $20$ relays, and $5$ destinations; the curves with the circle symbols represent that the multiple relay M2M network has $6$ antennas, $10$ relays, and $5$ destinations; the curves with the star symbols represent that the multiple relay M2M network has $6$ antennas, $5$ relays, and $5$ destinations. We run simulations for the proposed iteration strategies using Method $1$ to select the source antenna for each destination with $3$ iterations in Figure 3. As shown in the figure, in order to improve the performance of the multiple relay M2M network, we can increase the number of antennas or relays with the same relay power constraint. It is also confirmed that the performance of the proposed iteration strategy with the total power constraint (solid lines) is superior than that of the iteration strategy with the individual power constraints (dashed lines). The reason is that in the case of total power constraints, more power can be allocated to the relay nodes with good channel gains by optimization.

Figure 3

Expected worst-case received SINR of the proposed iteration strategies with the total relay power constraint (solid lines) and the individual relay power constraints (dashed lines).

7. Discussion

In the above results, the source M2M device must calculate the collaborative relay beamforming weights with the aid of perfect CSI and then broadcast them to the relays, respectively. Thus, the perfect CSI from the source to relays and from relays to the destinations needs to provide feedback to the source M2M device. The source can also know the noise level of relays and destinations from feedback channels. However in practice, wireless M2M networks, the feedback channels are always limited and could be overloaded if the number of relays and destinations is large. To deal with this issue, we can develop the distributed relay beamforming strategies to reduce the CSI exchange and allow each relay to learn its beamforming weight based on the local CSI. It is assumed that each relay node can learn the local CSI from the source by training and from destinations by feedback, respectively, and measure its noise level. In [16], a distributed strategy with local CSI is proposed for maximizing the received SNR at the single destination based on the Karush-Kuhn-Tucker (KKT) analysis. For the optimization problem of the received SINR at destinations, most previous works have focused on reducing the CSI exchange overhead. Instead of the instant CSI feedback, statistical channel information and partial CSI feedback are used in [28]. In [29], the CSI exchange overhead is drastically reduced by partitioning relays into two clusters which have different needs for CSI for calculating the relay beamforming weights. However, distributed relay beamforming strategies with the local CSI for maximizing the worst-case received SINR at multiple destinations have not been intensively discussed. In order to derive distributed beamforming strategies for the received SINR problem in the nonorthogonal systems, we can consider another criterion, that is, MMSE criterion to construct distributed relay beamforming strategies in [13] and [30]. It has been proved to be adopted for the single destination in the multiple relay systems to develop distributed beamforming strategies with the local CSI [13] and can be expected to be suitable for multiple destinations [30]. It is demonstrated that the MMSE minimization problems can be converted into the maximization of the worst-case SINR problems. Thus, in the future work, we can develop the maximization of the worst-case SINR distributed relay beamforming strategies by using the MMSE distributed beamforming strategies.

8. Conclusion

In this paper, we proposed two iterative strategies to jointly optimize the source-antenna selection and the AF relay beamforming weights with the guaranteed QoS in terms of SINR for multiple destinations in wireless M2M networks. We presented two methods to select the source antenna for each destination: (1) use the greedy algorithm to select the best available channel for each destination; (2) balance the channel condition for each destination and its interference to the other destinations based on Method $1$ . The maximization of the worst-case received SINR for collaborative relay beamforming was studied under two different types of power constraints, which are the total relay power constraint and the individual relay power constraints, respectively. To obtain the solution of the optimal relay beamforming weights, an SDR technique was employed to formulate the optimization problem as an SDP problem, which can be optimally solved. Simulation results demonstrated that after $3$ iterations, the performance of the proposed iterative strategies can obtain near-optimal performance. In the end of this paper, we discussed the distributed relay beamforming strategy to maximize the worst-case received SINR for multiple destinations in the practical M2M networks.

Footnotes

Acknowledgment

This work has been supported by the National Basic Research Program of China (973 Program, no. 2009CB320403).

References

Zhang

Shan

Yang

Mobile cellular networks and wireless sensor networks: toward convergence

IEEE Communications Magazine 2012 50 3 164 169

Andreev

Galinina

Koucheryavy

Energy-efficient client relay scheme for machine-to-machine communication

Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '11)

December 2011

Houston, Tex, USA

1 5

Khajehnouri

Sayed

A. H.

Distributed MMSE relay strategies for wireless sensor networks

IEEE Transactions on Signal Processing 2007 55 7 I 3336 3348

2-s2.0-34347381098

10.1109/TSP.2007.894410

Sadek

A. K.

Liu

K. J. R.

Multinode cooperative communications in wireless networks

IEEE Transactions on Signal Processing 2007 55 1 341 355

2-s2.0-33947376466

10.1109/TSP.2006.885773

Alatossava

Taparugssanagorn

Holappa

V. M.

Ylitalo

Measurement based capacity of distributed MIMO antenna system in Urban microcellular environment at 5.25 GHz

Proceedings of the IEEE 67th Vehicular Technology Conference-Spring (VTC '08)

May 2008

Marina Bay, Singapore

430 434

2-s2.0-47749154498

10.1109/VETECS.2008.101

Liang

Y. W.

Schober

Amplify-and-forward multi-antenna beamforming with joint source-relay power constraint

Proceedings of the IEEE 72nd Vehicular Technology Conference Fall (VTC '10)

September 2010

Ottawa, Canada

1 5

2-s2.0-78649398655

10.1109/VETECF.2010.5594311

Laneman

J. N.

Tse

D. N. C.

Wornell

G. W.

Cooperative diversity in wireless networks: efficient protocols and outage behavior

IEEE Transactions on Information Theory 2004 50 12 3062 3080

2-s2.0-5044252003

10.1109/TIT.2004.838089

Kramer

Gastpar

Gupta

Cooperative strategies and capacity theorems for relay networks

IEEE Transactions on Information Theory 2005 51 9 3037 3063

2-s2.0-26444619804

10.1109/TIT.2005.853304

Luo

Yue

Optimized distributed MIMO for cooperative relay networks

IEEE Communications Letters 2010 14 1 9 11

2-s2.0-75449115026

10.1109/LCOMM.2010.01.091457

10.

Luo

Z. Q.

W. K.

Zhang

Semidefinite relaxation of quadratic optimization problems

IEEE Signal Processing Magazine 2010 27 3 20 34

2-s2.0-77951427135

10.1109/MSP.2010.936019

11.

Boyd

Vandenberghe

Convex Optimization 2004

Cambridge, UK

Cambridge University Press

12.

Choi

Distributed beamforming using a consensus algorithm for cooperative relay networks

IEEE Communications Letters 2011 15 4 368 370

2-s2.0-79954630945

10.1109/LCOMM.2011.012511.102156

13.

Choi

MMSE-based distributed beamforming in cooperative relay networks

IEEE Transactions on Communications 2011 59 5 1346 1356

2-s2.0-79956258275

10.1109/TCOMM.2011.030411.100158

14.

Nguyen

D. H. N.

Nguyen

H. H.

Power allocation in wireless multiuser multi-relay networks with distributed beamforming

IET Communications 2011 5 14 2040 2051

15.

Elkheir

G. A.

Lioumpas

A. S.

Alexiou

Energy efficient AF relaying under error performance constraints with application to M2M networks

Proceedings of the IEEE 22nd International Symposium on Personal Indoor and Mobile Radio Communications (PIMRC '11)

September 2011

Toronto, Canada

56 60

16.

Zheng

Wong

K. K.

Paulraj

Ottersten

Collaborative-relay beamforming with perfect CSI: optimum and distributed implementation

IEEE Signal Processing Letters 2009 16 4 257 260

17.

Zheng

Wong

K. K.

Paulraj

Ottersten

Robust collaborative-relay beamforming

IEEE Transactions on Signal Processing 2009 57 8 3130 3143

2-s2.0-68249162098

10.1109/TSP.2009.2021456

18.

Nguyen

D. H. N.

Nguyen

H. H.

Tuan

H. D.

Distributed beamforming in relay-assisted multiuser communications

Proceedings of the IEEE International Conference on Communications (ICC '09)

June 2009

Dresden, Germany

1 5

2-s2.0-70449461367

10.1109/ICC.2009.5198922

19.

Nguyen

D. H. N.

Nguyen

H. H.

Pham

T. T.

Distributed beamforming in multiuser multi-relay networks with guaranteed QoS

Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '09)

December 2009

Honolulu, Hawaii, USA

1 6

2-s2.0-77951573371

10.1109/GLOCOM.2009.5425247

20.

Liu

Petropulu

A. P.

Cooperative beamforming in multi-source multi-destination relay systems with SINR constraints

Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '10)

March 2010

Dallas, Tex, USA

2870 2873

2-s2.0-78049380185

10.1109/ICASSP.2010.5496181

21.

Zheng

Blostein

Downlink distributed beamforming through relay networks

Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '09)

December 2009

Honolulu, Hawaii, USA

1 6

2-s2.0-77951520921

10.1109/GLOCOM.2009.5426225

22.

Phan

K. T.

Nguyen

D. H. N.

Le-Ngoc

Joint power allocation and relay selection in cooperative networks

Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '09)

December 2009

Honolulu, Hawaii, USA

1 5

2-s2.0-77951527947

10.1109/GLOCOM.2009.5425752

23.

Yoo

Jindal

Goldsmith

Multi-antenna downlink channels with limited feedback and user selection

IEEE Journal on Selected Areas in Communications 2007 25 7 1478 1491

2-s2.0-34547499778

10.1109/JSAC.2007.070920

24.

Mobasher

Taherzadeh

Sotirov

Khandani

A. K.

A Near-maximum-likelihood decoding algorithm for MIMO systems based on semi-definite programming

IEEE Transactions on Information Theory 2007 53 11 3869 3886

2-s2.0-36348978204

10.1109/TIT.2007.907472

25.

Bengtsson

Ottersten

Optimal and suboptimal transmit beamforming

Handbook of Antennas in Wireless Communications 2001

CRC Press

26.

Sidiropoulos

N. D.

Davidson

T. N.

Luo

Z. Q.

Transmit beamforming for physical-layer multicasting

IEEE Transactions on Signal Processing 2006 54 6 I 2239 2251

2-s2.0-33744542350

10.1109/TSP.2006.872578

27.

Grant

Boyd

CVX users guide for CVX version 1. 21 (build 808)

2011, http://cvxr.com/

28.

Hammarwall

Bengtsson

Ottersten

Utilizing the spatial information provided by channel norm feedback in SDMA systems

IEEE Transactions on Signal Processing 2008 56 7 3278 3293

2-s2.0-46649117659

10.1109/TSP.2008.920484

29.

Esli

Wittneben

Distributed multiuser cooperative network with heterogenous relay clusters

Proceedings of the 41st Asilomar Conference on Signals, Systems and Computers (ACSSC '07)

November 2007

Pacific Grove, Calif, USA

829 834

2-s2.0-50249173983

10.1109/ACSSC.2007.4487333

30.

Verdu

Multiuser Detection 1998

Cambridge, UK

Cambridge University Press