arxiv: v1 [cs.it] 14 Apr 2016

Transcription

1 1 Beamforming for Simultaneous Wireless Information and Power Transfer in Two-Way Relay Channels Wei Wang Rui Wang and Hani Mehrpouyan Department of Communication Engineering Nantong University Jiangsu China. Department of Information and Communications Tongji University Shanghai China. Department of Electrical and Computer Engineering Boise State University Boise USA. s: arxiv: v1 [cs.it] 14 Apr 016 Abstract This paper studies simultaneous wireless information and power transfer (SWIPT) systems in two-way relaying (TWR) channels. Here two source nodes receive information and energy simultaneously via power splitting (PS) from the signals sent by a multi-antenna relay node. Our objective is to imize the weighted sum of the harvested energy at two source nodes subject to quality of service (QoS) constraints and the relay power constraints. Three well-known and practical two-way relay strategies are considered i.e. amplify-andforward (AF) bit level XOR based decode-and-forward (DF- XOR) and symbol level superposition coding based DF (DF- SUP). For each relaying strategy we formulate the joint energy transmit beamforming and PS ratios optimization as a nonconvex quadratically constrained problem. To find a closed-form solution of the formulated problem we decouple the primal problem into two subproblems. In the first problem we intend to optimize beamforming vector for a given PS ratio. In the second subproblem we optimize the PS ratio with a given beamforming vector. It is worth noting that although the corresponding subproblem are nonconvex the optimal solution of each subproblem can still be found by using certain techniques. We provide numerical results that demonstrate the advantage of adapting the different relaying strategies and weighted factors to harvest energy in two-way relaying channel. Index Terms Beamforming energy harvesting simultaneous wireless information and power transfer (SWIPT) two-way relaying (TWR) power splitting (PS). I. INTRODUCTION Energy harvesting (EH) from surrounding environments is an emerging solution to prolong the operational time of energy-constrained nodes in wireless networks [1] []. Compared with conventional energy sources radio frequency (RF) signals can carry both information and energy simultaneously. Simultaneous wireless information and power transfer (SWIPT) has recently drawn significant attention where SWIPT has been investigated for various wireless channels [3] [11]. For example a point-to-point single-antenna additive white Gaussian noise (AWGN) channel was first studied in [3] where the authors used a capacity-energy function to study the fundamental performance tradeoff for simultaneous information and power transfer. Later on SWIPT was extended to a frequency selective channels in [4]. The authors in [5] studied This work is supported by the National Natural Science Foundation of China under grant the Natural Science Foundation of Jiangsu Province under grant BK the Project of Flagship-Major Construction of Jiangsu Higher Education Institutions of China under grant /00 and NASA Early Investigation Grant. SWIPT for fading channels subject to time-varying co-channel interference. In [6] [7] SWIPT schemes for multiple-inputmultiple-output (MIMO) channels were considered. The transmit beamforming design was studied for SWIPT in multipleinput-single-output (MISO) broadcast channels in [8] [9]. Moreover SWIPT has been investigated in other physical layer setups such as the OFDM and more. in [10] [11]. Besides the above studies related to one-hop transmission SWIPT technique has also been extended to wireless relay networks [1] [0]. For the one-way single-antenna relay channel two protocols namely time switching (TS) and power splitting (PS) are proposed for amplify-and-forward (AF) relay networks in [1] [13]. Later on SWIPT was extended to a full-duplex wireless-powered one-way relay channel in [14] [15] where the data and energy queues of the relay are updated simultaneously in every time slot. Because twoway relaying (TWR) is able to simultaneously enlarge wireless coverage and enhance spectral efficiency the SWIPT protocols for TWR channel recently have attracted much attention. In [16] the authors provided a SWIPT protocol in two-way AF relaying channels where two sources exchange information via an energy harvesting relay node. The authors investigated the sum-rate imization problem in two-way AF relaying channels in [17] where two source nodes harvest energy from multiple relay nodes. In [18] The authors studied the relay beamforming design problem for SWIPT in a non-regenerative two-way multi-antenna relay network. The authors investigated a compute-and-forward (CF) relay networks optimal beamforming design problem in [19] where two source nodes harvest energy by SWIPT from relay nodes. Moreover for different relaying strategies the authors in [0] studied the sum-throughput imization problem in a two-way AWGN relay channel where all nodes are powered by EH. A. Motivation So far most studies on SWIPT in relay networks focused on energy-constrained relay nodes [1] [16] [0]. As a matter of fact the sensor nodes or other low-power devices often have very limited battery storage and require an external charging sources to remain active in wireless cooperative or sensor networks. Although replacing or recharging batteries provides a solution to this problem it may incur a high cost and sometimes even be unavailable due to some physical or economic limitations. As show in Fig. 1 when the sensor node

2 Fig. 1. Example applications of SWIPT in TWR systems with battery-limited sensor nodes. that is inside the body or embedded in a building structure is depleted of energy it cannot fulfill its role any longer unless the source of energy is replenished. Therefore EH in such kind of scenarios where the relay node serves as the energy source is particularly important as it can provide a much safer and much more convenient solution. Hence in this paper we consider a TWR SWIPT system with batterylimited source nodes and a relay node that acts also as a source of energy (as described with more detail in Section II). Under this setup the source nodes receive information and energy simultaneously from the signals sent by a relay node. Furthermore to enhance bandwidth efficiency and power transfer we consider a scenario where the relay node is equipped with multiple antennas. This setup applies to lots of practical wireless transmission scenarios. Since TS can be regarded as a special case of PS with only binary split power ratios [5] [6] we focus our study on PS receivers instead of TS receivers. B. Related Works To the best of our knowledge only three papers e.g. [17] [18] and [19] thus far studied multi-antenna TWR SWIPT systems with the battery-limited source nodes. In [17] and [18] the authors considered TWR SWIPT systems and studied the optimal relay beamforming problem which imizes the weighted sum rate subject to the transmit power constraint at relay and the EH constraint at source nodes. However the authors in [17] assumed that the source node is able to decode information and extract power simultaneously which as explained in [6] may not hold in practice. In contrast to [17] the authors in [18] considered the case of separated EH and information decoding (ID) receivers. The energy receiver could harvest energy from the signals transmitted by both the source and the relay and the source node can only receive the information from the signals forwarded by the relay. Furthermore in their works only AF relaying strategy was considered. Authors in [19] considered a similar overlay protocol as ours. However it focused on the relay transmit power minimization problem for a TWR SWIPT network with CF relaying strategy. In this paper we study a TWR based PS-SWIPT system where TWR consists of a multi-antenna relay node and two single-antenna source nodes. Here two source nodes receive information and energy simultaneously via PS from the signals sent by the relay node. In particular we consider three types relaying strategies: AF bit level XOR based decode-and-forward (DF-XOR) and symbol level superposition coding based DF (DF-SUP). Different from [17] [19] our objective is to imize the weighted sum of the harvested energy at two source nodes subject to a given minimum signal-to-interference-and-noise ratio (SINR) constraint at source nodes and a imum transmit power constraint at the relay node. This scenario is of particular interest in supporting sensor networks that empower the next set of applications such as internet of things and also medical monitoring as outlined in Fig. 1. In these applications a large number of sensors will be operating in close vicinity. Hence SINR is a important metric for maintaining a given throughput while imizing energy transfer of the sensors or source nodes by the relay. The latter imizes the operational time of the sensors which can be a important metric in the scenarios presented in Fig. 1. To the authors best knowledge the joint beamforming and PS optimization for this new setup has not been studied in existing works. C. Our Contributions Under the above setup a TWR based PS-SWIPT system is considered in this paper. Different from existing works we assume that source nodes can receive information and energy simultaneously via PS from the signals sent by the relay node. Moreover various two-way relaying strategies may result in different transmit signals at the relay node. The impact of various relaying strategies on the amount of harvest energy has not been considered in existing studies. Besides for another challenging doubly-near-far problem [] in the TWR SWIPT system which refers to a node far away from the relay harvests much lower energy but consumes more to transmit data than a node near the relay could be mitigated effectively by setting different EH priorities for different source nodes. The main contributions of this work are summarized as follows. Firstly this is the first work to investigate joint beamforming and PS optimization for a TWR SWIPT system with battery-limited source nodes. Here two source nodes receive information and energy simultaneously via PS from the signals sent by a multi-antenna relay node. To achieve this goal we propose a two-phase PS-based relaying protocol. Secondly based on the above system setup we present different transmit signals at the relay node by considering three practical twoway relaying strategies due to their implementation simplicity [1] [] i.e. AF DF-XOR and DF-SUP. To explore the performance limit of the system for each relaying strategy we formulate the joint energy transmit beamforming and PS ratios optimization as a nonconvex quadratically constrained problem. Thirdly for each nonconvex relay beamforming optimization problem we find a solution by decoupling the primal problem into two subproblems. The first subproblem

3 3 only optimizes the beamforming vectors. We solve this nonconvex problem by applying the technique of semidefinite programming (SDP) [3]. The second subproblem only includes the PS ratios. We propose an novel algorithm to find the optimal closed-form solutions by separating the latter nonconvex subproblem into eight cases. Then a near optimal solution of the original optimization problem is found based on a two-tie iterative algorithm. Finally we provide numerical results for each relaying scheme to evaluate the performance of the proposed optimal beamforming designs. It is shown that when the priority and the distance of two source nodes are symmetric the DF-XOR relaying strategy performs better than the other two strategies. While the distances of two source nodes are asymmetric for all three considered relaying strategies the furthest node can harvest more energy when its energy weight factor is set to a larger value which can provide an effective solution to the doubly-near-far problem []. D. Organization The remainder of this paper is organized as follows. The TWR SWIPT system model is described in Section II. In Section III the weighted sum-power harvested imization problems are formulated for different relaying strategies. The solutions for the associated optimization problems by using suitable optimization tools are presented in Section IV. In Section V numerical simulation results are provided. Finally the paper is concluded in Section VI. Notations: Boldface lowercase and uppercase letters denote vectors and matrices respectively. For a square matrixaa T A A H Tr(A) Rank(A) and A denote its transpose conjugate conjugate transpose trace rank and Frobenius norm respectively. A 0 indicates that A is a positive semidefinite matrix. vec(a) denotes the vectorization operation by stacking the columns of A into a single vector a. E( ) denotes the statistical expectation. denotes the Kronecker product. denotes the XOR operator. 0 and I denote the zero and identity matrix respectively. The distribution of a circular symmetric complex Gaussian vector with mean vector x and covariance matrix Σ is denoted by CN(xΣ). C x y denotes the x y domain of complex matrices. II. SYSTEM MODEL Consider a half-duplex TWR system where two singleantenna source nodes S 1 and S exchange information with each other through an N-antenna relay node R as shown in Fig.. The channel matrices from S 1 and S to the relay are denoted by h 1 and h respectively and the channel matrices from the relay to S 1 and S are denoted by g 1 and g respectively. To further improve the spectral efficiency the two-time slot TWR model is used to realize bidirectional communication. Throughout this paper the following set of assumptions are made: The source nodes cannot communication with each other directly. Hence all messages are sent through the relay. This Fig.. A two-time slot TWR system where each source node coordinates information decoding and energy harvesting. occurs when the direct link is blocked due to long-distance path loss or obstacles [4] [5]. The relay is connected to the power grid which implies that it has access to reliable power at all times. However the source nodes are powered by the energy limited batteries or capacitors and need to replenish their energy by wireless power transfer. Amongst the different relaying protocols AF DF-XOR and DF-SUP schemes are applied at the relay node due to their implementation simplicity [1] []. Quasi-static block fading channels are assumed here i.e. channels are unchanged in a time-slot of T but change from time slot to time slot. The use of such channels is motivated by prior research in this field [5] [9] [15] [17] and practical consideration. As shown in Fig. 3 we propose a two-phase PS-based protocol for the TWR system. In the first phase of duration T/ two source nodes S 1 and S deliver their information to the relay node R simultaneously. In the second phase with the remaining time duration T/ the received information signal at R is processed by the aforementioned relaying strategies and then forwarded to the source nodes. Note that here by assuming a PS ratio ρ the transmit signal from the relay is used to simultaneously achieve information and power transfer. Based on the above system setup the received signal at the relay after the first phase i.e. the multiple access (MAC) phase is given by y R = h 1 x 1 +h x +n R (1) where x i for i {1} represents the transmit signal from node S i h i C N 1 is the channel vector from node S i to the relay node R and n R denotes the additive complex Gaussian noise vector at the relay following CN(0σ r I N). Each transmit signalx i is assumed to satisfy an average power constraint i.e. E( x i ) = P i. Upon receiving y R the relay node performs certain processing and then forwards its signal to the source nodes in the second phase also referred as broadcast (BC) phase. Let the transmit signal from the relay be denoted by x R = x 1 +x () where x 1 is the combined signal consisting of the messages from two nodes by using physical-layer network coding (PLNC). Note that here besides x 1 we also include a new

4 4 be the power splitting ratio meaning that 1 ρỹ i is used for information processing. As a result after converting the received signal to baseband and performing self-interference cancelation the obtained signal is denoted as Fig. 3. Energy harvesting and information processing relaying protocol based on PS with splitting ratio ρ. signal x which provides us with more degrees of freedom to optimize power transfer from relay to the source nodes. It is worth noting that various two-way relaying strategies may result in different transmit signal x R. The three relaying strategies we considered namely AF DF-XOR and DF-SUP are all favorable for practical implementation and the precoding designs based on these strategies are mathematically tractable [1] []. The primary focus of this work is to imize the weighted sum of the harvested power based on practical two-way relay strategies while meeting a minimum SINR for each source node. III. RELAYING STRATEGIES AND OPTIMIZATION PROBLEM FORMULATIONS Based on the channel setup described in Section II in this section we shall present different transmit signals x R for the TWR SWIPT system by considering three practical two-way relay strategies. Moreover to explore the system performance limit we also formulate three optimization problems for these approaches in this section. A. AF Relay Strategy With the AF relaying strategy the relay transmit signal x R in () can be expressed as x R = x 1 +x = Wh 1 x 1 +Wh x +Wn R +x (3) where W represents the precoding matrix used at the relay. In addition we assume that the relay node has the imum transmit power P r i.e. Tr{E(x R x H R )} P r which is equivalent to P 1 Wh 1 +P Wh +Tr(Q x)+σ r W F P r (4) where Q x = E(xx H ) is the covariance matrix of x and P 1 and P are the transmit powers used at nodes S 1 and S respectively. Then the received signals at the two nodes in the second T/ time interval are given by ỹ i = g T i Wh i x i +gt i Wh i x i +g T i x+gt i Wñ R +n id (5) where i = if i = 1 and i = 1 if i =. Note that here ỹ i x i x and ñ R denote the signals in the RF band and n id is the additive Gaussian noise due to the receiving antenna that follows CN(0σid ) for i {1}. Specifically as shown in Fig. 3 the received signal ỹ i at each end node is split into two portions for EH and information processing. Let ρ (01) y i = 1 ρ(g T i Wh i x i +g T i x+g T i Wñ R +n id )+n ic (6) where n ic CN(0σic ) is the additive Gaussian noise introduced by the signal conversion from RF band to baseband. Accordingly the SINR at the node S i is given by SINR AF i = P i g T i Wh i g T i Q xg i +σ r g T i W +σ id + σ ic 1 ρ. (7) Moreover the other portion of the received signal ρỹ i is used for EH. Since the background noise at the EH receiver is negligible and thus can be ignored [6] the harvested energy E i during EH time T/ is given by E i = ηt ρ( gt i Wh i P i + g T i Wh i P i +g T i Q x g i) (8) where η is the energy conversion efficiency with 0 < η < 1 which depends on the rectification process and the EH circuitry [6]. Note that in (8) the self-interference can be used for EH which is different from information processing. Our design goal is to imize the weighted sum power harvested at two EH nodes which is defined as the harvested energy minus the consumed energy. The corresponding optimization problem can be formulated as P 1P ρwq x 0 s.t. SINR AF i τ i i = 1 P i P i i = 1 Tr{E(x R x H R )} P r. α(e 1 P 1T )+β(e P T ) In (9) α and β correspond to the given energy weights for the two EH receivers S 1 and S respectively where a larger weight value indicates a higher priority of transferring energy to the corresponding EH receiver as compared to other EH receiver. τ i and P i are the SINR requirement and the imum transmit power at node S i respectively. B. DF-XOR Relay Strategy If the relay node adopts the DF relaying strategy it needs to decode the messages sent from both source nodes in the MAC phase and then transmit a function of the two messages in the BC phase. The rate region of the MAC channel is characterized by [6] [7] C MAC ( R 1 R ) = R 1 log (1+ P1 h1 σ ) r R log (1+ P h σ ) r R 1 + R log det(i+ P1 σ r h 1 h H 1 + P σ rh h H ) (9) (10) where R 1 and R are the transmit rates at nodes S 1 and S respectively.

5 5 We assume that the messages sent from the nodes in the MAC phase can be successfully decoded at the relay node. Let b i denote the decoded bit sequence from S i for i {1 }. With the DF-XOR relaying strategy the combined bit sequence is yielded as b 1 = b 1 b. Then the transmit signal in the second time interval of T/ denoted by x R in () can be expressed as x R = s 1 +x (11) where s 1 is the modulated signal of bit sequence b 1. The corresponding relay power constraint for signal x R is denoted as Tr{E(x R x H R )} = Tr(Q s)+tr(q x ) P r (1) where Q s = E(s 1 s H 1 ) is the covariance matrix of s 1. Then the received signals at node S i in the RF band is given by ỹ i = gi T s 1 +gi T x+n idi = 1. (13) Using power splitting 1 ρỹi is used for information processing at the end nodes. After baseband conversion and self-interference cancelation the obtained signal is denoted as y i = 1 ρ(g T i s 1 +g T i x+n id)+n ic (14) Subsequently the SINR at the node S i can be determined as SINR XOR i = g T i Q sg i g T i Q xg i +σ id + σ ic 1 ρ. (15) On the other hand ρỹ i is used for EH at the node S i. The harvested energy E i is given as [6] E i = ηt ρ(gt i Q s g i +g T i Q x g i)i = 1. (16) Similarly our aim is to imize the weighted sum power harvested at the two source nodes subject to a given minimum SINR constraint at each source node and a imum transmit power constraint at the relay node. The corresponding optimization problem can be formulated as ρq s 0Q x 0 s.t. SINR XOR i τ i i = 1 Tr(Q s )+Tr(Q x ) P r. α(e 1 P 1T )+β(e P T ) (17) Note that here different from (9) it is not necessary to optimize P 1 and P as they are determined via the constraints presented in (10). C. DF-SUP Relay Strategy In this subsection we consider a case where the relay uses the DF-SUP relaying strategy. Again we assume that the messages sent from the nodes in the MAC phase can be successfully decoded at the relay node. Then by applying DF- SUP relaying strategy the transmit signal in the BC phase denoted by x R can be expressed as x R = s 1 +s +x (18) where s i is the modulated signal of bit sequence of node S i. The corresponding relay power constraint for signal x R is denoted as Tr{E(x R x H R )} = Tr(Q s1)+tr(q s )+Tr(Q x ) P r (19) where Q si = E(s i s H i ) is the covariance matrix of s i. Then the received signals at node S i in the RF band is given by ỹ i = g T i s 1 +g T i s +g T i x+n idi = 1. (0) Let us assume that 1 ρỹ i portion of the received signal is used for information processing at the end nodes. After converting this signal to baseband and performing selfinterference cancelation the obtained signal is given by y i = 1 ρ(g T i s i +gt i x+n id)+n ic (1) Then the SINR at the node S i can be denoted as SINR SUP i = g T i Q si g i g T i Q xg i +σ id + σ ic 1 ρ. () On the other hand ρỹ i is used for EH at the node S i. The harvested energy E i is given as [6] E i = ηt ρ(gt i Q s1gi +gt i Q sgi +gt i Q xgi ) (3) Also our design goal is to imize the weighted sum power harvested at the two EH nodes which is defined as the harvested energy minus the consumed energy. The corresponding optimization problem can be formulated as ρq s1 0Q s 0Q x 0 α(e 1 P 1T ) +β(e P T ) (4) s.t. SINR SUP i τ i i = 1 Tr(Q s1 )+Tr(Q s )+Tr(Q x ) P r. Here similar to (17) we are not necessary to optimize P 1 and P as they are determined by (10). For different relaying strategies we have formulated three weighted sum-power harvested imization problems in (9) (17) and (4). In the following sections we will propose three algorithms to solve these optimization problems. IV. OPTIMAL DESIGN OF THREE MAXIMIZATION PROBLEMS For the AF relaying strategy the optimization problem in (9) is nonconvex due to not only the coupled beamforming vectors {WQ x } and the remaining parameters {P i ρ} in both the SINR and transmitted power constraints but also all the quadratic terms involving W. In general it is difficult or even intractable to obtain the global optimal solution to a nonconvex problem [9] [18]. However it is well known that a function can be imized by first imizing over some of the variables and then imizing over the remaining ones [8 Sec 4.1.3]. Thus when P i and ρ are first fixed the resulting beamforming optimization problem reduces to that of a conventional nonconvex problem with a

6 6 rank-one constraint. The latter can be efficiently solved by using some rank relaxation techniques [3]. Moreover we note that when the beamforming vectors {WQ x } are fixed the resulting beamforming optimization problem over {P i ρ} is still a nonconvex problem. However as show later the optimal solution can be obtained in closed-form by separating this subproblem into eight cases. In the following we first decouple problem (9) into two subproblems that can be solved separately and then propose a two-tie iterative algorithm to obtain the near optimal solution of the original optimization problem. Finally similarly we decouple problems (17) and (4) into two subproblems. Note that here different from (9) two subproblems from (17) and (4) only involve the beamforming vectors and PS ratios. A. Joint Beamforming and PS Optimization for AF relaying strategy Let us solve the two subproblems stemming from (9). In the first subproblem we try to find the solutions of W and Q x for fixed P 1 P and ρ values. Then we update the values of P 1 P and ρ by fixing the remaining parameters. These subproblems are updated in an alternating manner. 1) Optimize W and Q x for fixed P 1 P and ρ: Note that when fixing P 1 P and ρ the problem of optimizing variables W and Q x is equivalent to WQ x 0 αρ( g T 1 Wh P + g T 1Wh 1 P 1 +g T 1 Q x g 1)+βρ( g T Wh 1 P 1 + g T Wh P +g T Q x g ) s.t. SINR AF i τ i i = 1. P 1 Wh 1 +P Wh +Tr(Q x ) (5) +σr W F P r Although problem (5) has a simpler form than the original problem in (9) it is still a nonconvex problem. To find the optimal solution of problem (5) we conduct some further transformations on (5). To be specific we transform g1 TWh and Wh 1 into their equivalent forms as and similarly g T 1Wh = Tr(g T 1Wh h H W H g 1) (6a) = Tr(g 1 gt 1 Wh h H WH ) (6b) = w H (h ht g 1 gt 1 )w (6c) = Tr((h ht g 1 gt 1 )wwh ) (6d) Wh 1 = Tr(Wh 1 h H 1 W H ) (7a) = Tr(IWh 1 h H 1 W H ) (7b) = w H (h 1 ht 1 I)w (7c) = Tr((h 1 ht 1 I)wwH ). (7d) where w = vec(w). In obtaining (6c) and (7c) we have used the identity Tr(ABCD) = (vec(d T )) T (C T A)vec(B). (8) Similar to (6) and (7) we apply the above transformations to other terms in (5). Let W ww H (5) can be rewritten as W 0Q x 0 Tr(A 1 W)+Tr(B1 Q x ) s.t. Tr(C i 1 W) Tr(τ i g ig T i Q x ) D i 1 i = 1. Tr(E 1 W)+Tr(Qx ) P r Rank( W) = 1. (9) where A 1 (P h h T +P 1 h 1h T 1) (αρg 1g T 1 + βρg gt ) B 1 (αρg 1 gt 1 +βρg gt ) Ci 1 (P i h i ht τ i i σr I) g i gt i Di 1 (σid + σ ic 1 ρ )τ i and E 1 (P 1 h 1 ht 1 +P h ht +σ ri) I. Due to the rank-one constraint finding the optimal solution of (9) is difficult. We therefore resort to relaxing it by deleting the rank-one constraint namely W 0Q x 0 Tr(A 1 W)+Tr(B1 Q x ) s.t. Tr(C i 1 W) Tr(τ i g ig T i Q x ) D i 1 i = 1. Tr(E 1 W)+Tr(Qx ) P r (30) It is noted that if the problem (30) has a rank-one optimal solution of W i.e. the problem (9) has a rank-one optimal solution of W the problem (9) is equivalent to the problem (5). Fortunately we have the following lemma. Lemma 1: The rank-one optimal solution of the problem (30) always exists. Proof: Please refer to Appendix A. By acquiring the optimal rank-one solution of (30) we can further get the optimal solution of (9) and then the optimal solution of (5). ) Optimize P 1 P and ρ for fixed W and Q x : In the second step we need to optimize the power P 1 P and the power ratio ρ with the remaining variables fixed. The corresponding optimization problem can be formulated as P 1P ρ α(e 1 P 1T )+β(e P T ) s.t. SINR AF i τ i i = 1. P 1 Wh 1 +P Wh +Tr(Q x ) +σ r W F P r 0 < P i P i i = 1. 0 < ρ < 1 (31) Similar to (6) and (7) we apply the transformations in (31). The problem of optimizing the variables P 1 P and

7 7 ρ is equivalent to P 1P ρ A ρp +B ρp 1 αp 1 βp +C ρ (3a) s.t. (E P D )(1 ρ) τ 1 σ 1c (3b) (G P 1 F )(1 ρ) τ σ c (3c) P 1 J +P K P r L (3d) 0 < P 1 P 1 (3e) 0 < P P (3f) 0 < ρ < 1. (3g) where A αηt gt 1 Wh + βηt gt Wh B αηt gt 1 Wh 1 + βηt gt Wh 1 C αηt gt 1 Q xg1 + βηt gt Q xg D (g1 TQ xg1 + σ r gt 1 W + σ 1d )τ 1 E g1 TWh F (g TQ xg +σ r g TW +σ d )τ G g T Wh J Wh 1 K Wh and L Tr(Q x )+σr W F. Problem (3) is still nonconvex in its current form since both the SINR and transmit power constraints involve coupled transmit power P i and power ratio ρ. To find the optimal solution of (3) we give the following lemma. Lemma : Let {P1P ρ } denote an optimal solution of problem (3) we have (1) for the optimal solution {P1P ρ } either the SINR or the transmit power constraint must hold with equality; () the optimal solution {P1 P ρ } can be obtained in closed-form by comparing the following eight cases: When the two SINR constraints (3b) and (3c) hold with equality the optimal solution {P 1P ρ } are given by τ σ c P1 1 ρ +F = P = G ρ a1 +a a 3 = 1. a 1 τ 1σ 1c 1 ρ +D E (33) where a 1 (A D G + B E F + C E G ) a A E τ 1 σ 1c + B E τ σ c + A D G and a 3 αe τ σ c +βg τ 1 σ 1c. When the constraints (3b) and (3d) hold with equality the optimal solution {P 1P ρ } are given by P1 = P r L P K P J = ρ = 1 b 1 +b. J E b 3 τ 1σ 1c 1 ρ +D E (34) where b 1 (A J B K )τ 1 σ 1c b (αk βj )τ 1 σ 1c and b 3 (P re L E D K )B +(A D +C E )J J E. When the constraints (3b) and (3e) hold with equality the optimal solution {P 1P ρ } are given by τ 1σ1c P1 = P 1 P = 1 ρ +D E ρ c c 1 = 1. E c 3 (35) where c 1 A τ 1 σ 1c c βτ 1 σ 1c and c 3 A D +B E P 1+C E E. When the constraints (3c) and (3d) hold with equality the optimal solution {P 1P ρ } are given by τ σc P1 = 1 ρ +F P G = P r L P1 J K ρ = 1 d 1 +d. K G d 3 (36) where d 1 (B K A J )τ σ c d (βj αk )τ σ c and d 3 (P rg L G F J )A +(B F +C G )K G K. When the constraints (3c) and (3f) hold with equality the optimal solution {P 1P ρ } are given by τ σc P1 = 1 ρ +F ρ = 1 G e e 1 G e 3. P = P (37) where e 1 B τ σ c e ατ σ c and e 3 B F +A G P +C G G. When the two transmit power constraints (3d) and (3e) hold with equality the optimal solution {P 1 P ρ } are given by P1 = P 1 P = P r L J P 1 K ρ = min{1 τ 1σ1c E P D 1 τ σ (38) c G P1 F }. When the constraints (3d) and (3f) hold with equality the optimal solution {P1P ρ } are given by P1 = P r L K P P = P J ρ = min{1 τ 1σ1c E P D 1 τ σ (39) c G P1 F }. When the constraints (3e) and (3f) hold with equality the optimal solution {P1P ρ } are given by P1 = P 1 P = P ρ = min{1 τ 1σ1c E P D 1 τ σc (40) G P1 F }. Proof: Please refer to Appendix B. We compare all objective function values by substituting (33) (40) into (3a) and select one {P1 P ρ } as the optimal solution if they lead to the greatest value of the objective function f(ρ ).

8 8 The proposed iterative algorithm to obtain the near optimal solution to problem (9) is summarized in the following Algorithm 1. Algorithm 1 Finding a near optimal solution to problem (9) Initialize P 1 P and ρ; Repeat Update the beamforming matrixes W and Q x for fixed P 1 P and ρ using the following steps: First solve problem (30) by CVX and obtain the optimal solution as W and Q x. Then derive the optimal solution W of (5) by eigenvalue decomposition (EVD) of W ; Update P 1 P and ρ with W and Q x fixed using the following steps: First check whether there exists some feasible solutions {P 1 P ρ } of problem (3). If yes compute all the objective function values and select one {P 1 P ρ } corresponding to the largest value of objective function as the optimal solution to (3). Otherwise exit the algorithm; Until The difference between the value of objective function in (9) from one iteration to another is smaller than a pre-fixed predetermined threshold. B. Joint Beamforming and PS Optimization for DF-XOR relaying strategy In this subsection we consider optimization problem (17) where the relay node adopts the DF-XOR two-way relaying strategy. Similar to problem (9) we decouple problem (17) into two subproblems. It is worth noting that here different from (9) two subproblems from (17) only involve the beamforming vectors and PS ratios where P 1 and P are not necessary to be optimized as they are determined via the constraints presented in (10). 1) Optimize Q s and Q x for fixed ρ: Note that when fixing ρ the problem of optimizing variables Q s and Q x can be equivalent to Q s 0Q x 0 α(g T 1 Q s g 1 +g T 1Q x g 1) s.t. SINR XOR i τ i i = 1. Tr(Q s )+Tr(Q x ) P r which is rewritten as Q s 0Q x 0 +β(g T Q sg +gt Q xg ) Tr(A 3 Q s )+Tr(A 3 Q x ) s.t. Tr(B 3 Q s ) Tr(τ 1 B 3 Q x ) D 3 Tr(C 3 Q s ) Tr(τ C 3 Q x ) E 3 Tr(Q s )+Tr(Q x ) P r (41) (4) where A 3 αg 1 gt 1 + βg gt B 3 g 1 gt 1 C 3 g gt D 3 σ1d + σ 1c 1 ρ τ 1 ande 3 σd + σ c 1 ρ τ. It is easy to verify that (4) is a standard SDP problem. Thus its optimal solution {Q s Q x } can be easily obtained using existing software e.g. CVX [9]. ) Optimize ρ for fixed Q s and Q x : In the second step we need to optimize the PS ratio ρ with the remaining variables fixed. The corresponding optimization problem can be formulated as ρ αηt ρ(gt 1 Q sg 1 +gt 1 Q xg 1 ) αt P 1 + βηt ρ(gt Q sg +gt Q xg ) βt P s.t. SINR XOR i τ i i = 1. which is equivalent to ρ (A 4 +B 4 )ρ T (αp 1 +βp ) s.t. C 4 (1 ρ) τ 1 σ 1c D 4 (1 ρ) τ σ c αηt (43) (44) where A 4 (gt 1 Q sg1 + gt 1 Q xg1 ) B 4 βηt (gt Q sg + gt Q xg ) C 4 g1 TQ sg1 (gt 1 Q xg1 + σ1d )τ 1 and D 4 gq T s g (g Q T x g+σ d )τ. According to the definition of P 1 andp in (10) the simplified PS design problem yields the following ρ (A 4 +B 4 )ρ s.t. ρ 1 τ 1σ 1c C 4 (45) ρ 1 τ σc. D 4 It can be observed that the objective function in (45) achieves a higher value when one of the SINR constraints holds with equality. Hence the optimal solutionρ = min{1 τ1σ 1c C 4 1 τ σc D 4 } can be obtained from problem (45). To summarize Algorithm below summarizes the solution to (17). Note that Algorithm differs from Algorithm 1 in two main aspects: First in step 1) the optimal beamforming matrixes {Q sq x} can be easily obtained due to the absence of rank-one constraint; and second step ) only involves PS ratios ρ P 1 and P are not necessary to be optimized. Algorithm Finding the near optimal solution to problem (17) Initialize ρ; Repeat solve problem (4) by CVX and obtain the optimal solution of Q s and Q x in (41) which are denoted by Q s and Q x respectively; Update PS ratios ρ using (45) for fixed Q s and Q x; Until Termination criterion is satisfied. C. Joint Beamforming and PS Optimization for DF-SUP relaying strategy In this subsection we consider that the DF-SUP relaying strategy is adopted at the relay node. To find the optimal solution of the joint optimization problem in (4) We similarly decouple problem (4) into two subproblems and then propose a two-tie iterative algorithm to obtain a near optimal solution of the original optimization problem. 1) Optimize Q s1 Q s and Q x for fixed ρ: In the first step we need to optimize the beamforming matrices Q s1

9 9 Q s and Q x with the PS ratio ρ fixed. The corresponding optimization problem can be formulated as Q s1 0Q s 0Q x 0 α(g T 1 Q s1 g 1 +g T 1Q s g 1 +g T 1 Q xg 1 )+β(gt Q x1g +g T Q s g +g T Q x g ) s.t. SINR SUP i τ i i = 1. Tr(Q s1 )+Tr(Q s )+Tr(Q x ) P r which is rewritten as Q s1 0Q s 0Q x 0 Tr(A 5 (Q s1 +Q s +Q x )) s.t. Tr(B 5 Q s ) Tr(τ 1 B 5 Q x ) D 5 Tr(C 5 Q s1 ) Tr(τ C 5 Q x ) E 5 Tr(Q s1 )+Tr(Q s )+Tr(Q x ) P r. (46) (47) where A 5 αg 1g T 1 + βg g T B 5 g 1g T 1 C 5 g g T D 5 σ1d + σ 1c 1 ρ τ 1 and E 5 σd + σ c 1 ρ τ. Note that (47) is a standard SDP problem. Thus its optimal solution can be easily obtained [9]. ) Optimize ρ for fixed Q s1 Q s and Q x : In the second step we need to optimize the PS ratio ρ with the remaining variables fixed. The corresponding optimization problem can be formulated as ρ αηt ρ(gt 1 Q s1g 1 +gt 1 Q sg 1 +g T 1Q x g 1) αt P 1 + βηt ρ(gt Q s1 g +g T Q sg +gt Q xg ) βt P s.t. SINR SUP i τ i i = 1. which is equivalent to ρ (A 6 +B 6 )ρ T (αp 1 +βp ) s.t. C 6 (1 ρ) τ 1 σ 1c D 6 (1 ρ) τ σ c (48) (49) where A 6 αηt (gt 1 Q s1g1 +gt 1 Q sg1 +gt 1 Q xg1 ) B 6 βηt (gt Q s1 g +gq T s g +gq T x g) C 6 g1 T Q s g1 (g1 TQ xg1 + σ 1d )τ 1 and D 6 g TQ s1g (gt Q xg + σd )τ. Since P 1 and P are determined based on the first phase problem (49) is simplified as ρ (A 6 +B 6 )ρ s.t. ρ 1 τ 1σ 1c C 6 (50) ρ 1 τ σc. D 6 Similar to the problem (45) the optimal PS solution ρ = min{1 τ1σ 1c C 6 1 τσ c D 6 } can be obtained from problem (50). Similar to Algorithm the proposed iterative algorithm to problem (4) is summarized in Algorithm 3. Algorithm 3 Finding the near optimal solution to problem (4) Initialize ρ; Repeat Update the beamforming matrixes Q s1 Q s and Q x with (47) for fixed ρ; Update PS ratios ρ using (50) for fixed Q s1 Q s and Q x; Until Termination criterion is satisfied. V. SIMULATION RESULTS In this section we numerically evaluate the performance of the proposed energy harvesting schemes. The channel vector h i and g i are set to be Rayleigh fading i.e. the elements of each channel matrix or vector are complex Gaussian random variables with zero mean and unit variance. The channel gain is modeled by the distance path loss model [7] given as g ij = c d n ij where c is an attenuation constant set as 1 n is the path loss exponent and fixed at 3 and d ij denotes the distance between nodes i and j. For simplicity we assume that the noise power at all the destinations are the same i.e. σic = σ id = σ r = σ = 0 dbm i and the energy conversion efficiency η = 50%. Moreover the imum transmit powers at the two sources if not specified are set as P 1 = P = P = 5 dbm. In all simulations the weighted sum power of the relay network is computed by using 1000 randomly generated channel realizations. In Figs. 4-6 we present the weighted sum power harvested for different relaying strategies at different distance of two sources when the relay node is equipped with N = 4 transmit antennas. In Figs. 4 we illustrate the harvested power of the two source nodes when the relay node adopts AF relaying strategy. Specifically the distances of the two source nodes are symmetric i.e. d RS1 = d RS = 1 meter. For comparison two different energy weights with α = β = 1/ and α = 4/5 and β = 1/5 are simulated for this scenario. From simulation results in Fig. 4(a) we find that when S 1 and S have the same priority the two nodes can achieve a fair energy efficiency. However in Fig. 4(b) when S 1 and S have different priorities node S 1 can harvest more energy since its energy weight factor is set to a larger value. Moreover note that when the relay transmit power P r is low the harvested energy at the nodes is negative which implies that the harvested power at nodes from the relay is smaller than the consumed power for signal transmission. In Figs. 5 and 6 we illustrate the harvested power at the two source nodes i.e. P S1 and P S and the shared power ratio P at the S node i.e. S P S1 +P S in different rate requirements R 1 = R = γr and R 1 = R = γ3r where R = 1 log (1 + P σ ). It is noted that in asymmetric scenario i.e. d RS1 = 1 meter and d RS = meters in Fig. 5(a) and with d RS1 = 1 meter and d RS = 3 meters in Fig. 6(a) although two source nodes S 1 and S have same priority the node S still harvests much lower energy in different rate requirements. The main reason is that the location of S is far away from the relay node R which could result in very small channel gain as compared to the near

10 10 (a) (a) (b) Fig. 4. Performance comparison for AF Relaying Strategy when changing P r. (a) Harvested power of the two source nodes with same priority. (b) Harvested power of the two source nodes with different priority. node. This coupled effect is referred to as the doubly-near-far problem []. However when with higher priority i.e. β = 3/4 and β = 6/7 in Fig. 5(b) and with β = /3 and β = 9/10 in Fig. 6(b) we find that node S all can share more power for the harvested total power in different rate requirements. This indicates that under the asymmetric scenario the far node will be able to harvest more energy when its energy weight factor is set to a larger value which can provide an effective solution to the doubly-near-far problem. In addition note that here different the results in Figs. 4(a) and 4(b) the harvested energy at the nodes is positive which implies that the harvested power at nodes from the relay is greater than the consumed power for signal transmission. In Fig. 7 we illustrate the weighted sum power harvested for different relay strategies with different number of antennas at relay when changingp r. For fair comparison the priorities the rate requirements and the distances of two source nodes are set to be the same i.e. α = β = 1/ R 1 = R = 0.1R and d RS1 = d RS = 1 meter and simulated for each scenario. From simulation results we find that the DF- XOR relaying strategy achieves the best performance and the DF-SUP relaying strategy outperforms the AF relaying Fig. 5. Performance comparison for DF-XOR Relaying Strategy at γ = 0.1 when changing P r. (a) Harvested power of the two source nodes with same priority. (b) Power ratio shared by S with different priority. strategy. This indicates that DF relaying strategy has a higher EH efficiency due to the assumption that the relay has enough processing ability to correctly decode the received signals. Moreover combining the information using XOR is better than using superposition since the power of the relay node can be used more efficiently in the DF-XOR relaying strategy. In addition it is observed that when the number of transmit antennas increases (N = 4 8) the all three considered two-way relaying strategies achieve better performance. This demonstrates the significant benefit by applying large or even massive antenna arrays for efficiently implementing TWR SWIPT systems in practice. Finally in Fig. 8 we compare the proposed joint beamforming and PS optimization scheme with the other two schemes i.e. only precoding without power and PS ratio allocation and only power and PS ratio allocation without precoding scheme. For fairness the setting of each node is the same with the one in Fig. 7. From simulation results for three considered twoway relaying strategies we find that the joint beamforming and PS optimization scheme all achieves the best performance as it uses the degrees of the freedom of both power PS ratio allocation and precoding. Moreover Fig. 8 also shows that the (b)

11 11 (a) Fig. 7. Performance comparison for different relay strategies with different number of antennas at relay when changing P r. (b) Fig. 6. Performance comparison for DF-SUP Relaying Strategy at γ = 0.1 when changing P r. (a) Harvested power of the two source nodes with same priority. (b) Power ratio shared by S with different priority. Fig. 8. Performance comparison with different schemes at α = β = 1/ R 1 = R = 0.1R and d RS1 = d RS = 1. precoding only scheme can improves the system performance and it performs much better than the power and PS ratio allocation only scheme. asymmetric for all three considered relaying strategies the furthest node can harvests more energy when its energy weight factor is set to a larger value. VI. CONCLUSIONS This paper has studied the joint energy transmit beamforming and power splitting design for a multi-antenna TWR system based simultaneous wireless information and power transfer (SWIPT). The weighted sum power harvested at all source nodes was imized subject to given SINR constraints at different source nodes and transmitted power constraints at relay node. Considering three different relaying strategies the design problems are formulated as three nonconvex quadratically constrained problem which are decoupled into two subproblems that can be solved separately by applying suitable optimization tools. The performance of three relay strategies were compared and some practical implementation issues were discussed. Simulation results showed when the priority and the distance of two source nodes are symmetric the DF-XOR relaying strategy performs better than the other two strategies. While the distances of two source nodes are APPENDIX A PROOF OF LEMMA 1 Note that Problem (30) is a quasi-convex optimization problem its optimal solution can be easily obtained [9]. Let us denote the optimal solution of W and Qx in (30) by W andq x respectively. It is easy to verify that W is an optimal solution of the following optimization problem W 0 Tr(A 1 W) s.t. Tr(C i W) 1 D1 i +Tr(τ i gig i T Q x) i = 1. (51) Tr(E 1 W) Pr Tr(Q x ) According to Lemma 3.1 in [3] there exists an optimal solution W for the problem (51) such that (Rank( W )) 3 (5)

12 1 It can be verified that W 0. Thus from (5) we have Rank( W ) = 1. The proof of Lemma 1 is thus completed. APPENDIX B PROOF OF LEMMA Suppose problem (3) is feasible and let {P 1 P ρ } and f( ) denote the optimal solution and the objective function respectively. Next we show that for problem (3) with the optimal solution {P 1P ρ } either the SINR constraints or the transmit power constraints must hold with equality. Now we prove this result by contradiction. Namely if the above conditions are not satisfied we can find another solution of (3) denoted by {P 1 P ρ } which achieves a higher value for the objective function f(p 1 P ρ ) than f(p 1P ρ ). First suppose that the two SINR constraints do not hold with equality. In this case note that in problem (3) the two SINR constraints are equivalent to and ρ 1 τ 1σ 1c EP D (53) ρ 1 τ σc GP 1 F. (54) It can be observed that when PS solution ρ = min{1 τ 1σ1c EP 1 τσ c D GP 1 F } i.e. one of the SINR constraints holds with equality the objective functionf(p 1 P ρ ) = (A ρ β)p + (B ρ α)p 1 + C ρ achieves a higher value as the same transmit power constraints. Hence this assumption is not true. Next consider the case that the transmit power constraints do not hold with equality. In this case if we want the PS solution ρ to increase to achieve a higher value of objective function the transmit power solution P i increases. Then for the three transmit power constraints (3d) (3e) and (3f) there exists at least a power constraint that holds with equality. Thus the assumption that the transmit power constraints do not hold with equality also cannot be true. To summarize for the optimal solution {P1 P ρ } either the SINR or transmit power constraint must hold with equality. Moreover given the optimal solution {P1 P ρ } we can prove that two constraints of problem (3) are achieved with equality. First suppose that the optimal solution {P1 P ρ } can be obtained if and only if the SINR constraint (3b) holds with equality. In this case we can easily find another solution of P 1 for (3) while two transmit power constraints (3d) or (3e) hold with equality. We denote this solution of P 1 by P 1. We can verify that this solution achieves a higher value for the objective function f( P 1 P ρ ) than f(p1 P ρ ). Hence this assumption cannot be true. Similarly for all the other assumptions where if and only if a constraint holds with equality we can easily prove these assumptions cannot be true too. In conclusion for the optimal solution {P1P ρ } there exist at least two constraints of problem (3) holding with equality. Based on the above observations we can separate the problem (3) into ten cases by setting the two SINR constraints (3b) (3c) and the three transmit power constraints (3d) (3e) (3f) hold with equality. However in the ten cases when the constraints (3b) and (3f) hold with equality if we want the objective function f( ) to increase the transmit power solution P 1 increases which leads to the conclusion that the constraints (3d) or (3e) are active. Hence this combination is included in (3d) and (3f) or (3e) and (3f) so that this case can be removed. Similarly constraints (3c) and (3e) combination will also be removed since this case was contained in (3d) and (3e) or (3e) and (3f). Therefore the optimal solution {P 1 P ρ } of problem (3) is able to be obtained in closedform by comparing following eight cases. When the two SINR constraints (3b) and (3c) hold with equality we obtain the following two equations τ 1σ1c P = 1 ρ +D P1 E = τ σ c 1 ρ +F G (55) By substituting (55) into (3a) the objective function f(ρ ) can be equivalently written as τ 1σ 1c τ σ f(ρ ) =A ρ 1 ρ +D c +B ρ 1 ρ +F E G τ σ c τ 1σ 1 ρ +F 1c α 1 ρ +D β +C ρ G E which is further equivalent to f(ρ ) = a 1(ρ ) +a ρ a 3 E G (1 ρ ) (56) a 4. (57) where a 1 a a 3 are defined as in (33) and a 4 αf /G + βd /E. Hence problem (3) is simplified as ρ f(ρ ) s.t. P 1J +P K P r L 0 < P 1 P 1 0 < P P 0 < ρ < 1. (58) To proceed to solve (58) we have the following lemma. Lemma 3: The optimal solution ρ = 1 a1+a a 3 a 1 can be obtained from problem (58) while 0 < a a 3 < a 1. Proof: Please refer to Appendix C. Then we obtain the optimal solution in (33). When the constraints (3b) and (3d) hold with equality we obtain the following two equations τ 1σ 1c P 1 ρ +D = P1 = P r L PK (59) E J By substituting (59) into (3a) the objective function f(ρ ) can be equivalently written as f(ρ ) = b 3J E (ρ ) +(b 3 J E +b 1 )ρ +b J E (1 ρ ) b 4 (60)

13 13 where b 1 b b 3 are defined as in (34) and b 4. (60) is equivalent to (P re L E D K )α+d J β J E f(ρ ) = b 3 (1 ρ b 1 +b )+ J E (1 ρ ) b 1 +b 3 b 4. J E (61) Then problem (3) is equivalent to the following problem s.t. ρ f(ρ ) τ σ c 1 ρ +F G P 1 P 1 0 < P P (6) 0 < ρ < 1. Similar to Lemma 3 when b 3 > 0 and b 1 + b < 0 the objective function f(ρ ) must have a imum value which can be further derived from b 3 (1 ρ b ) = 1+b J E (1 ρ ). On the other hand to guarantee the optimal solution ρ satisfying 0 < ρ < 1 we have b 1 +b +J E b3 > 0 which results in an optimal solution ρ of problem (6) given as ρ = 1 b 1 +b. (63) J E b 3 we thus obtain the solution given in (34). When the constraints (3b) and (3e) hold with equality we obtain the following two equations τ 1σ1c P 1 ρ +D = P1 = P 1 (64) E By substituting (64) into (3a) the objective function f(ρ ) can be written as f(ρ ) = c 3E (ρ ) +(c 3 E +c 1 )ρ c E (1 ρ ) c 4 (65) where c 1 c c 3 are defined as in (35) and c 4. (65) is equivalent to αe P 1+βD E f(ρ ) = c 3 (1 ρ )+ c 1 c E (1 ρ ) c 1 +c 3 c 4. (66) E Then problem (3) is simplified as ρ f(ρ ) τ σc s.t. P1 1 ρ +F G P1 J +P K P r L 0 < P P (67) 0 < ρ < 1. Due to the fact c 3 > 0 i.e. c 3 < 0. Similar to Lemma 3 if c1 c E < 0 i.e. c > c 1 the objective function f(ρ ) must have a imum value which can be inferred from a fact that c 3 (1 ρ ) = c1 c E (1 ρ ). Note that to guarantee the optimal solution ρ satisfying 0 < ρ < 1 c must satisfy c < c 1 +E c 3. As a result the optimal solutionρ of problem (67) can be derived as ρ c c 1 = 1. (68) E c 3 Then we obtain the optimal solution in (35). When the constraints (3c) and (3d) hold with equality we obtain the following two equations P 1 = τ σ c 1 ρ +F G P = P r L P 1J K (69) By substituting (69) into (3a) the objective function f(ρ ) can be equivalently written as f(ρ ) = d 3K G (ρ ) +(d 3 K G +d 1 )ρ +d K G (1 ρ ) d 4 (70) where d 1 d d 3 are defined as in (36) and d 4. (70) is equivalent to (P rg L G F J )β+f K α G K f(ρ ) = d 3 (1 ρ )+ Then problem (3) can be rewritten as s.t. d 1 +d K G (1 ρ ) d 1 K G +d 3 d 4. ρ f(ρ ) τ 1σ 1c 1 ρ +D E 0 < P 1 P 1 0 < ρ < 1. P P (71) (7) Similar to Lemma 3 when d 3 > 0 and d 1 + d < 0 the objective function f(ρ ) must have a imum value which can be further derived from d 3 (1 ρ d ) = 1+d K G (1 ρ ). When d 1 + d + K G d 3 > 0 the optimal solution ρ of problem (7) can be derived as ρ = 1 d 1 +d K G d 3. (73) Then we obtain the optimal solution in (36). When the constraints (3c) and (3f) hold with equality we obtain the following two equations τ σc P1 = 1 ρ +F P G = P (74) By substituting (74) into (3a) the objective function f(ρ ) can be equivalently written as f(ρ ) = e 3G (ρ ) +(e 3 G +e 1 )ρ e G (1 ρ ) e 4 (75) where e 1 e e 3 are defined as in (37) and e 4. (75) is further equivalent to βg P +αf G f(ρ ) = e 3 (1 ρ )+ e 1 e G (1 ρ ) e 1 G +e 3 e 4. (76)