Next Article in Journal
Magnetic-based Soft Tactile Sensors with Deformable Continuous Force Transfer Medium for Resolving Contact Locations in Robotic Grasping and Manipulation
Previous Article in Journal
Kinetic Energy Harvesting for Wearable Medical Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precoding Design for Energy Efficiency Maximization in MIMO Half-Duplex Wireless Sensor Networks with SWIPT

1
School of Information and Electrical Engineering, Hebei University of Engineering, Handan 056038, China
2
School of Mechanical and Electrical Engineering, Handan University, Handan 056038, China
3
Key Laboratory of Trustworthy Distributed Computing and Service, Ministry of Education, School of Information and Communication Engineering, Beijing University of Posts and Telecommunications (BUPT), Beijing 100876, China
4
Key Laboratory of System Control and Information Processing, Ministry of Education, Shanghai Jiaotong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(22), 4923; https://doi.org/10.3390/s19224923
Submission received: 20 October 2019 / Revised: 7 November 2019 / Accepted: 8 November 2019 / Published: 12 November 2019
(This article belongs to the Section Electronic Sensors)

Abstract

:
This paper explores the energy efficiency (EE) maximization problem in single-hop multiple-input multiple-output (MIMO) half-duplex wireless sensor networks (WSNs) with simultaneous wireless information and power transfer (SWIPT). Such an energy efficiency maximization problem is considered in two different scenarios, in which the number of energy-harvesting (EH) sensor nodes are different. In the scenario where the single energy-harvesting sensor node is applied, the modeled network consists of two multiple-antenna transceivers, of which the energy-constrained energy-harvesting sensor node harvests energy from the signals transmitted from the source by a power splitting (PS) scheme. In the scenario of multiple EH sensor nodes, K energy-constrained sensor nodes are applied and the same quantity of antennas are equiped on each of them. The optimization problem is formulated to maximize the energy efficiency by jointly designing the transceivers’ precoding matrices and the PS factor of the energy-harvesting sensor node. The considered constraints are the required harvested energy, the transmission power limit and the requirement on the data rate. The joint design of the precoding matrices and the PS factor can be formulated as an optimization problem, which can be transformed into two sub-problems. An alternating algorithm based on Dinkelbach is proposed to solve the two sub-problems. The convergence of the proposed alternating algorithm, the solution optimality and the computational complexity are analyzed in the paper. Simulation results demonstrate the convergence and effectiveness of our proposed algorithm for realizing the maximum energy efficiency.

1. Introduction

The emergence of tens of billions of intelligent Internet-of-Things (IoT) devices have attracted serious concerns about their-energy waste and the radiation pollution. Green communications technology opens a new road for energy exploitations and gradually promotes its popularity nowadays [1]. Simultaneous wireless information and power transfer (SWIPT) was first proposed in [2], and can be regarded as the candidate technology that meets the idea of green communications because it aims at realizing a high data rate and at enabling energy reuse. It enables the energy-constrained nodes to harvest recyclable energy and to decode information on the same radio frequency (RF) signals. The tendency inspires a feasible scheme to extend the lifetime of the energy constrained nodes in wireless sensor networks. In SWIPT, the EH receiver and information-decoding (ID) receiver can be functioned in two modes, i.e., co-located and separated [3]. For the co-located receiver, Zhang, R. et al. in [4] proposed two efficient information transmission mechanisms, known as power-splitting (PS) and time-switching (TS). In the PS mechanism, the co-located receiver divides the received signal power into two power streams in an appropriate proportion, one for information decoding and the other for energy harvesting. In the TS mechanism, the co-located receiver is designed to switch its modes from energy harvesting to information decoding in different time slots.
Following these pioneering works, SWIPT began its applications in the single-input single-output (SISO) system [5] and is now developed to be used in scenarios, e.g., the multiple-input single-output (MISO) systems [6,7,8], multiple-input multiple-output (MIMO) systems [9,10,11], broadcast systems [12], orthogonal frequency-division multiplexing (OFDM) relay systems [13], etc. In [8], the two-way multi-antenna relay networks are considered, where the beamforming of relay antennas is designed to maximize the sum rate under the energy constraint. Jiang, R. et al. in [12] studied the broadcast networks system, where the transmitting power of a hybrid access point (H-AP) equipped with multiple antennas is minimized under constraints such as the minimum signal-to-interference-plus-noise ratio (SINR) and the power harvested at each user with a single antenna. Shen, Y. et al. in [13] investigated the resource optimization problem in TS relay networks, which take into account a realistic energy consumption model. It is worth noting that MIMO technology has been widely used in modern mobile communication networks due to its considerable diversity gain and effective resistance to path loss. The combination of MIMO and SWIPT is the cutting edge research direction at present. Following this research trend, Zhang, R. et al. in [4] studied the rate-energy trade-off region under the constraint of nonlinear energy harvesting. By giving an actual EH model, the optimal trade-off between the harvested energy and data rate was explored in [4]. By optimizing the precoding at the source and the relay, the achievable transmission rate with TS or with PS schemes, respectively, in a two-hop decode-and-forward (AF) relay network are investigated in [9]. In [11], the two-user MIMO interference channel is proposed and secure communication is realized by designing optimal precoding on transmitters.
However, few of the above studies were concerned with the energy efficiency, but energy efficiency has developed to be the performance metric of energy-constrained networks due to the requirement on energy management [14,15,16,17]. Energy efficiency is defined as the ratio of the network sum rate to the total energy consumption. In [10], all of the desired receivers decode information and harvest energy in a point-to-point communication system with a PS scheme. The desired receivers can achieve the optimal energy efficiency by jointly optimizing the transmission precoding matrix, the artificial noise covariance matrix and the power splitting ratios. Not only that, the probability of eavesdropping is also reduced. In [14], Ng, D.W.K. et al. study the resource allocation algorithm design for OFDM downlink systems to maximize the energy efficiency. An alternating iteration algorithm is developed in [15] to infer the optimal precoding matrix at the relay and transmitters in the two-way relay networks. Zhou, X. et al. in [16] studied a two-way amplify-and-forward relay networks. An iterative algorithm based on the constrained concave convex procedure (CCCP) is proposed to jointly design the sources’ and relay’s precoding matrices and the PS ratios. The iterative algorithm maximizes the energy efficiency of the network. Similarly, the energy efficiency is maximized in [17], in which a distributed antenna system is studied. Through adjusting the PS ratios, this enables the IoT devices to coordinate the energy harvesting and information decoding. The users of the wireless powered communication networks in [18] can harvest energy from dedicated power station, which ensures that the information can be transmitted to the information reception station. Then the authors Wu, Q. et al. in [18] maximize the energy efficiency of the system via the joint optimization of slot allocation and power control. A utility function is proposed in a MISO system [19] with the purpose of maximizing the ratio of the achieved utility to the total power consumption. A MISO downlink data transmission link is considered in [20], in which the zero-forcing (ZF) algorithm was applied by the base station. The energy efficiency of the system is maximized through joint optimization of beamforming and power splitting.
The realizations of the above works all benefited from the development of EH technologies. They paved the way for filling up the deficiencies of the energy constrained networks. In recent years, EH technologies continue to be developed, i.e., recyclable energy can be harvested from either the surrounding environment or the dedicated RF signals and the lifetime of wireless networks can be improved [21]. Gui, L. et al. in [22] propose a new secure system model with both fixed and mobile jammers to guarantee secrecy in the worst-case scenario, where all jammers are thought to have EH capability. In [23], sensor nodes harvest energy from mobile phones’ RF signals and, in order to acquire available power use from the ambient mobile phones, a contract theory based incentive scheme is investigated. With the linear EH model, the harvested energy in unit time is usually proportional to the power of the received RF signals. By analyzing the experimental results, Boshkovska, E. et al. in [24] validate that for the EH device there exists a nonlinear relationship between the power of output direct current (DC) and the power of input RF signals. However, as in the conclusion given in [17], the linear EH model is still applicable for the low input power devices. The main reasons are described as follows. Firstly, the nonlinear EH model can be approximated to be a piecewise linear model when the EH devices are assigned with relatively low or high input power, separately. With such assumption, the linear EH model has acceptable accuracy. Secondly, the intelligent IoT devices are usually exposed in an electromagnetic surrounding filled with low power RF signals.
In this paper, the maximization of energy efficiency is studied in MIMO half-duplex SWIPT wireless sensor networks with a linear EH model. In fact, some efforts have been made to achieve highly efficient utilization of network resources in WSNs. Guo, W. et al. in [25] propose a soft real-time fault-tolerant task allocation algorithm (FTAOA) for WSNs to support the fault tolerance mechanism to improve resource utilization. To reduce the total traffic cost, Yao, Y. et al. in [26] develop the data collection protocol EDAL. In [27], a node scheduling strategy is proposed to reduce the number of active nodes and the amount of messages in the network. The strategy proposed in [27] reduces the energy consumption and improves the network lifetime. Moreover, in order to ensure the communication reliability in wireless networks, the proposal in [28] studies the distributed scheduling problem, by which the scheduling cycle length is to be minimized.
This paper incorporates SWIPT with WSNs, in which the energy efficiency maximization problem is studied. The main contributions of this paper are summarized as follows:
  • The EE maximization problem in the SWIPT-based MIMO wireless sensor network is considered by employing the PS mechanism. Our goal is to maximize the energy efficiency by jointly optimizing the power splitting factor and the precoding matrices of transceivers. Moreover, the upper limit on the power supply and the required data rate for quality of service (QoS) guarantee is bounded as the constraints.
  • For the proposed optimization problem, the objective function is non-convex. The eigenvalue decomposition and singular value decomposition are applied to reduce the dimensionality. Meanwhile, the original coupling problem can be transformed to two sub-problems. For these two sub-problems, an alternating algorithm based on Dinkelbach is used to guarantee the convergence.
  • Excluding the convergence analysis, the computational complexity is given in the paper and the feasibility of our proposal is also verified by simulation results. The performance comparisons between our proposal and the other selected baseline algorithms validate the effectiveness of our proposal.
The rest of the paper is organized as follows: Section 2 introduces the system model. In Section 3, the optimization problem of power allocation for EE maximization is formulated and then an alternating algorithm based on Dinkelbach is introduced to solve the problem. The convergence and computational complexity of the algorithm are also analyzed in Section 3. Section 4 shows the extended energy efficiency maximization problem in the scenario of the multiple EH sensor nodes used. Section 5 shows the simulation results and performance comparisons. Section 6 concludes the paper.

2. Single EH Sensor Nodes System Model

2.1. Notations

This paper specifies that the bold italics uppercase and the bold italics lowercase letters denote matrices and vectors, respectively. The operators, i.e., A H , A 1 , A + , T r ( A ) , A and A represent conjugate transpose, inverse, pseudo-inverse, trace, determinant and Frobenius norm of the matrix A , respectively. A 0 means that A is positive semidefinite. C N × M is defined as the space of N × M matrices with complex entries. I N denotes the N × N identity matrix and 0 denotes zero matrix.

2.2. System Model

The MIMO half-duplex wireless sensor network is considered in Figure 1, which contains two transceivers, i.e., the source and the EH sensor node, equipped with N 1 and N 2 antennas, respectively. The source node has a persistent energy supplement, whereas the the EH sensor node is an energy-constrained node. To maintain the EH sensor node’s survival, the PS scheme is adopted at the sensor node to harvest energy from the source’s RF signals, which are divided into two power flows by the power splitting scheme. α ( 0 α 1 ) is defined as the power splitting factor, by which the α portion of the received signals is used for the ID receiver and the ( 1 α ) portion of the received signals is used for the RF energy conversion circuit. Following the same channel models given in [12,17], the frequency flat slow block fading channel model that obeys Rayleigh distribution is adopted in this paper. The channel state information (CSI) of all the links are also known beforehand according to the historical statistics.

2.2.1. Source—EH Sensor Node Communication

The two-way data transmission between the source and the sensor node occurs on two consecutive equal time slots, just as shown in Figure 1. During the time slot 1, the source transmits its signals to the sensor node and the received signals at the sensor node are given as
y 2 = HF s 1 + n 2 ,
where s 1 C N 1 × 1 ( s 1 C N ( 0 , 1 ) ) is the symbol vector generated at the source, F C N 1 × N 1 denotes the precoding matrix of the source, H C N 2 × N 1 denotes the channel coefficient matrix from the source to the sensor node. The channel can be assumed to be quasi-static [29], so that the channel coefficient matrix keeps unvaried in the two time slots. n 2 C N 2 × 1 is the zero-mean additive white Gaussian noises (AWGN) at the sensor node with the covariance matrix σ 2 2 I N 2 .
By the power splitting factor α , the received signals at the sensor node can be split into two flows. In the paper, each antenna of the sensor node is assigned with the same power splitting factor for problem simplification. Therefore, the signals for ID at the sensor node can be given as
y 2 I D = α HF s 1 + n 2 + n 2 I D ,
where n 2 I D C N ( 0 , σ I D 2 I N 2 ) is the AWGN generated by the baseband signal processing of the ID receiver, which converts the RF signals into baseband signals. Thus the achievable rate in bits/s/Hz at the sensor node can be given by
R 2 = log 2 I N 2 + α HF F H H H ( α σ 2 2 + σ I D 2 ) I N 2 .
The data rate R 2 in Equation (3) can be achieved due to a perfect CSI being known, which can be obtained by training at either the uplink or the downlink using pilot signals or by using the 1-bit feedback information [17].
On the one hand, the signals for EH are shown as in Equation (4)
y 2 E H = 1 α HF s 1 + n 2 .
Based on Equation (4), the harvested energy at the sensor node can be deduced as in Equation (5) [17]
E E H = η ( 1 α ) T r HF F H H H + σ 2 2 I N 2 ,
where η 0 , 1 is the EH efficiency that reflects the proportion of the harvested energy from the received RF signals. The considered MIMO sensor network always operates at low input power at the sensor node, the linear EH model can be still applicable to the MIMO sensor network.
During time slot 2, the sensor node is powered by the harvested energy and transmits signals to the source. The received signal at the source can be expressed as
y 1 = H H D s 2 + n 1 ,
where s 2 C N 2 × 1 denotes the unit-norm information symbol vector generated by the sensor node, D C N 2 × N 2 denotes the precoding matrix at the source and n 1 C N ( 0 , σ 1 2 I N 1 ) represents the AWGN with n 1 C N 1 × 1 . Similar to the achievable rate shown as Equation (3), the achievable rate in bits/s/Hz at the source can be expressed as Equation (7).
R 1 = log 2 I N 1 + 1 σ 1 2 H H D D H H .
The rate R 1 can be obtained just following the data rate given in Equation (3). Assume that the the time slot for transmission at the source or the sensor node occupies a 1 / 2 time frame; thus, the total amount of the transmitted data equals T 2 ( R 2 + R 1 ) in bits/Hz.

2.2.2. Energy Efficiency

The energy efficiency (bits/Hz/J) in the paper is defined as the ratio of total number of bits conveyed in networks to the networks’ energy consumption. According to Equations (1) and (6), the transmission power of the source and the sensor node denoted as P 1 and P 2 can be obtained as follows
P 1 = F 2 = T r ( Q 1 ) ,
P 2 = D 2 = T r ( Q 2 ) ,
where Q 1 = F F H , Q 2 = D D H . Therefore, the networks’ energy consumption within the considered time slots can be given by Equation (10) [16].
E ( α , Q 1 , Q 2 ) = T 2 ( P 1 + P 2 ) ( P 1 + P 2 ) ς ς + P c ,
where ς ( 0 , 1 ] is the drain efficiency for the power amplifiers of transceivers. Assume that the source and the sensor node have the same drain efficiency, and are set as 100%, i.e., ς = 1 . The factor T 2 is due to the power consumption, i.e., P 1 and P 2 separately occurred in two time slots. The period of each time slot is 1 / 2 unit time frame. P c denotes the sum of circuit power consumption at the source and the sensor node and P c = P c 1 + P c 2 , where P c 1 and P c 2 represent the circuit power consumption of the source and the sensor node, respectively, The circuit power is used to enable transceivers to transmit and receive the signals.
The harvested energy during the time slot should maintain the sensor node’s survival and hence is required to meet the minimum requirement for the power needs, therefore the constraint in Equation (11) can be established as
E E H P 2 + P c 2 .
According to the definitions given above, the EE in the paper can be formulated as follows
E E = T 2 ( R 2 + R 1 ) E ( α , Q 1 , Q 2 ) .

3. Problem Formulation and Solution

3.1. Optimization Problem Formulation

This paper aims at maximizing the network energy efficiency. It considers the joint optimization of the power splitting factor α and the precoding matrices Q 1 and Q 2 . The constraints such as the upper limit on the transmission power, the required minimum data rate for the QoS guarantee and the semi-definite attribute of the precoding matrices are taken into account in the problem formulation. To enable the passive sensor node to receive sustaining energy for its operation, the harvested energy during the two time slots should cover the sensor node’s lowest power needs. Therefore, the problem for maximizing energy efficiency in the considered network can be formulated as P1,
P 1 : max α , Q 1 , Q 2 T 2 ( R 2 + R 1 ) E ( α , Q 1 , Q 2 )
s . t .   E E H P 2 + P c 2
T r Q 1 P ¯ max 1
T r Q 2 P ¯ max 2
R 2 r ¯ min 2
R 1 r ¯ min 1
0 α 1 , Q 1 0 , Q 2 0 ,
where P ¯ max 1 and P ¯ max 2 stand for the upper limits on the transmission power of the source and the sensor node, respectively, r ¯ min 1 and r ¯ min 2 represent data rate threshold for the source and the sensor node, respectively, Q 1 0 and Q 2 0 imply that the precoding matrices Q 1 and Q 2 are positive semi-definite. The constraint of Equation (13b) guarantees that the energy harvested at the sensor node can cover its lowest power needs.

3.2. Proposed Optimal Solution

The formulated problem P 1 is non-convex due to the parameters α and precoding matrix Q 1 being highly coupled in both the objective function and the constraints of Equations (13b) and (13e). Therefore, it is generally quite difficult to solve the non-convex optimization problem. As far as we know, there exist no practical methods that can guarantee to converge to the global optimal solution. In this section, a diagonalization method with low complexity is used to solve the problem P 1 .
The precoding matrix Q 1 can be decomposed as Q 1 = U Q 1 Λ Q 1 U Q 1 H through eigenvalue decomposition (EVD). Similarly, the precoding matrix Q 2 can also be decomposed as Q 2 = U Q 2 Λ Q 2 U Q 2 H by EVD. If singular value decomposition (SVD) is applied to the channel coefficient matrix H , then H can be decomposed as H = U H Λ H 1 2 V H H . Based on these decomposition expressions, the problem P 1 can be transformed to problem P 2 shown as at the top of this page,
P 2 : max α , Λ Q 1 , Λ Q 2 log 2 I N 2 + α α σ 2 2 + σ I D 2 U H Λ H 1 2 V H H U Q 1 Λ Q 1 U Q 1 H V H Λ H 1 2 H U H H T r ( U Q 1 Λ Q 1 U Q 1 H ) + T r ( U Q 2 Λ Q 2 U Q 2 H ) + P c
+ log 2 I N 1 + 1 σ 1 2 V H Λ H 1 2 H U H H U Q 2 Λ Q 2 U Q 2 H U H Λ H 1 2 V H H T r ( U Q 1 Λ Q 1 U Q 1 H ) + T r ( U Q 2 Λ Q 2 U Q 2 H ) + P c
s . t .   η ( 1 α ) T r U H Λ H 1 2 V H H U Q 1 Λ Q 1 U Q 1 H V H Λ H 1 2 H U H H + σ 2 2 I N 2 T r ( U Q 2 Λ Q 2 U Q 2 H ) + P c 2
T r U Q 1 Λ Q 1 U Q 1 H P ¯ max 1
T r U Q 2 Λ Q 2 U Q 2 H P ¯ max 2
log 2 I N 2 + α α σ 2 2 + σ I D 2 U H Λ H 1 2 V H H U Q 1 Λ Q 1 U Q 1 H V H Λ H 1 2 H U H H r ¯ min 2
log 2 I N 1 + 1 σ 1 2 V H Λ H 1 2 H U H H U Q 2 Λ Q 2 U Q 2 H U H Λ H 1 2 V H H r ¯ min 1
0 α 1 ,
where Λ Q i = d i a g λ 1 , Q i , λ 2 , Q i , λ 3 , Q i , λ N i , Q i with λ 1 , Q i λ 2 , Q i λ 3 , Q i λ N i , Q i 0 ( i = 1 , 2 ) and Λ H 1 2 = d i a g λ 1 , H , λ 2 , H , λ 3 , H , λ min { N 1 , N 2 } , H with λ 1 , H λ 2 , H λ 3 , H λ min { N 1 , N 2 } , H 0 . U H C N 2 × min { N 1 , N 2 } and V H C N 1 × min { N 1 , N 2 } are all the unitary matrices, which are constituted by the orthogonal unit norm column vectors.
As to the maximization of network energy efficiency, the best precoding matrices Q 1 * and Q 2 * make the denominator of the objective function Equation (14a) a minimum and if they simultaneously maximize the numerator, then the matrices Q 1 * and Q 2 * can also achieve the maximum of Equation (14a). Inspired by the form of this fraction, if precoding matrices Q 1 and Q 2 are individually quantified with a definite value, the denominator of EE can correspondingly obtain a definite value. The numerator is the sum rates of two transceivers and no coupling relations among Q 1 and Q 2 appear in both expressions of R 1 and R 2 . When Q 1 and Q 2 are definite, the maximum of the objective function of P 2 can be obtained only if the function value of the numerator is maximized. The sum rates, i.e., R 1 plus R 2 given in the numerator, can be maximized by separating the sum rates into two individual items, i.e., R 1 and R 2 , and then by maximizing each of them. As the fundamentals for solving P 2 , two basic operation laws should be introduced at first.
Operation law 1: For any two matrices A C m × n and B C n × m , and det I + AB = det I + BA is always satisfied.
Operation law 2: For any two matrices A C m × n and B C n × m , and T r ( AB ) = T r ( BA ) is always true.
According to the operation laws 1 and 2, the objective function of P 2 can be restructured as in Equation (15),
E E = log 2 I N 2 + α α σ 2 2 + σ I D 2 V H Λ H V H H U Q 1 Λ Q 1 U Q 1 H + log 2 I N 1 + 1 σ 1 2 U H Λ H U H H U Q 2 Λ Q 2 U Q 2 H T r ( Λ Q 1 ) + T r ( Λ Q 2 ) + P c .
Denote Π 1 = Θ 1 Λ Q 1 Θ 1 H with Θ 1 = V H H U Q 1 and Π 2 = Θ 2 Λ Q 2 Θ 2 H with Θ 2 = U H H U Q 2 . According to the operation law 1, the energy efficiency shown as in Equation (15) can be further transformed to Equation (16), which can be regarded as the objective function of a newly generated optimization problem P 3 .
E E = log 2 I N 2 + α α σ 2 2 + σ I D 2 Λ H Π 1 + log 2 I N 1 + 1 σ 1 2 Λ H Π 2 T r ( Π 1 ) + T r ( Π 2 ) + P c .
P 3 : m a x Π 1 , Π 2 , α   log 2 I N 2 + α α σ 2 2 + σ I D 2 Λ H Π 1 + log 2 I N 1 + 1 σ 1 2 Λ H Π 2 T r ( Π 1 ) + T r ( Π 2 ) + P c
s . t .   η ( 1 α ) T r Λ H Π 1 + σ 2 2 I N 2 T r ( Π 2 ) + P c 2
T r Π 1 P ¯ max 1
T r Π 2 P ¯ max 2
log 2 I N 2 + α α σ 2 2 + σ I D 2 Λ H Π 1 r ¯ min 2
log 2 I N 1 + 1 σ 1 2 Λ H Π 2 r ¯ min 1
0 α 1 .
According to the result of the Hadamard inequality [30], given that Π 1 = Θ 1 Λ Q 1 Θ 1 H and Π 2 = Θ 2 Λ Q 2 Θ 2 H , the optimization variables Π 1 and Π 2 must be diagonal matrices. For better conveniences, their optimums in P 3 can be given as Π 1 * and Π 2 * . Furthermore, Π 1 and Π 2 are similar to the diagonal matrices Λ Q 1 and Λ Q 2 , respectively. Therefore, Π 1 * and Q 1 * should obtain the same eigenvalues and such a conclusion applies equally well to Π 2 * and Q 2 * .
Theorem 1.
For any two commutative matrices A and B , if A is positive definite Hermitian matrices and B is n × n positive semidefinite Hermitian matrices that have eigenvalues denoted as ω i i = 1 n and τ i i = 1 n , then the function value of A + B can be bounded by sorting their eigenvalues in multiplication calculations. For better clarity, Theorem 1 can be illustrated by a set of inequalities as in Equation (18),
i = 1 n ( 1 + τ ¯ i ω ¯ i ) I n + A 1 B i = 1 n ( 1 + τ ¯ i i ) ,
where the overline and underline represent the ascending and descending order of the eigenvalues, respectively.
Proof. 
Please refer to [15] for the proof of Property 2. □
Corollary 1.
For any two n × n matrices A and B , if A is positive semidefinite diagonal Hermitian matrices and B is positive semidefinite Hermitian matrices, the following inequality is always true
i = 1 n ( 1 + ω ̲ i + τ ¯ i ) I n + A + B i = 1 n ( 1 + ω ¯ i + τ ¯ i ) ,
where, ω i + i = 1 n represents the eigenvalues of A + .
Based on the result of Theorem 1 and the corollary, the EE in P 3 can be maximized such that Π 1 = Λ Q 1 , Π 2 = Λ Q 2 . In order to achieve above equation relations, U Q 1 = V H and U Q 2 = U H must be required. Meanwhile, regarding Λ H as A + ; therefore, the problem P 3 can be simplified to the problem P 4 as follows
P 4 : max α , λ Q 1 , λ Q 2 i = 1 N 1 log 2 1 + α λ i , H λ i , Q 1 α σ 2 2 + σ I D 2 + i = 1 N 2 log 2 1 + λ i , H λ i , Q 2 σ 1 2 i = 1 N 1 λ i , Q 1 + i = 1 N 2 λ i , Q 2 + P c
s . t .   η ( 1 α ) ( i = 1 N 1 λ i , H λ i , Q 1 + σ 2 2 ) i = 1 N 2 λ i , Q 2 + P c 2
i = 1 N 1 λ i , Q 1 P ¯ max 1
i = 1 N 2 λ i , Q 2 P ¯ max 2
i = 1 N 1 log 2 1 + α λ i , H λ i , Q 1 α σ 2 2 + σ I D 2 r ¯ min 2
i = 1 N 2 log 2 1 + λ i , H λ i , Q 2 σ 1 2 r ¯ min 1
0 α 1 , λ Q 1 0 , λ Q 2 0 ,
where λ Q 1 and λ Q 2 are the set of eigenvalues for Q 1 and Q 2 , respectively, i.e., λ Q 1 = λ i , Q 1 i = 1 , 2 , N 1 , λ Q 2 = λ i , Q 2 i = 1 , 2 , N 2 .
Due to the power splitting factor α and the eigenvalues λ i , Q 1 are coupled with each other in both the objective function and the constraint conditions of Equations (20b) and (20e), problem P 4 is still a non-convex optimization problem. However, problem P 4 can be further divided into two correlated sub-problems and their optimal solutions can still converge to the global optimum in the sequel.
Giving a function f ( x ) = g 1 ( x ) g 2 ( x ) , if g 1 ( x ) is a concave function and g 2 ( x ) is an affine function, then f ( x ) is defined as a pseudo-concave function. As to the objective function of problem P 4 , as shown in Equation (20a), its numerator can be regarded as a concave function with respect to λ Q 1 and λ Q 2 , while its denominator is an affine function of λ Q 1 and λ Q 2 , the objective function is hence a pseudo-concave function with respect to variables λ Q 1 and λ Q 2 . To verify the concavity, the Hessian matrix in terms of the variable α is given as in Equation (21), where the parameters in Equation (21) are all non-negative. The second derivative of α is negative and thus this validates that the objective function Equation (20a) is also a concave function with respect to α .
α 2 = σ I D 2 ln 2 i = 1 N 1 2 α σ 2 4 + 2 σ 2 2 σ I D 2 + 2 α σ 2 2 λ i , H λ i , Q 1 + σ I D 2 λ i , H λ i , Q 1 ( α 2 σ 2 4 + σ I D 4 + 2 α σ 2 2 σ I D 2 + α 2 σ 2 2 λ i , H λ i , Q 1 + α σ I D 2 λ i , H λ i , Q 1 ) 2 < 0 .
The concavity of the objective function about α indicates that there must be a closed-form solution for α . For each fixed α , there is always a set of feasible solutions for λ Q 1 , λ Q 2 . The optimal eigenvalues of the two precoding matrices will be obtained by calculating the closed-form solution through circular iteration. The analysis and resolving solutions are summarized as follows:
(1) Optimal λ Q 1 * , λ Q 2 * for Given α
For a given α , problem P 4 with respect to λ Q 1 , λ Q 2 can be reduced to P 5 :
P 5 : max λ Q 1 , λ Q 2 Ξ ( λ Q 1 , λ Q 2 ) Ξ ( λ Q 1 , λ Q 2 ) Ψ ( λ Q 1 , λ Q 2 ) Ψ ( λ Q 1 , λ Q 2 )
s . t .   η ( 1 α ) ( i = 1 N 1 λ i , H λ i , Q 1 + σ 2 2 ) i = 1 N 2 λ i , Q 2 + P c 2
i = 1 N 1 log 2 1 + α λ i , H λ i , Q 1 α σ 2 2 + σ I D 2 r ¯ min 2
λ Q 1 0 , λ Q 2 0
( 20 c ) ( 20 d ) ( 20 f ) ,
where Ξ ( λ Q 1 , λ Q 2 ) = i = 1 N 1 log 2 1 + α λ i , H λ i , Q 1 α σ 2 2 + σ I D 2 + i = 1 N 2 log 2 1 + λ i , H λ i , Q 2 σ 1 2 and Ψ ( λ Q 1 , λ Q 2 ) = i = 1 N 1 λ i , Q 1 + i = 1 N 2 λ i , Q 2 + P c . The constraint of Equation (22b) is affine with respect to λ Q 1 , λ Q 2 , while the constraint of Equation (22e) can be seen as a concave set of λ Q 1 , λ Q 2 . Meanwhile, the objective function of P 5 is a non-convex fractional form. Therefore, the Dinkelbach method [31] and CVX can be applied to solve this fractional programming problem P 5 .
The Dinkelbach method is considered to be an effective iterative algorithm to solve the fractional programming problem with multi-parameters. The Dinkelbach algorithm can convert the objective function from its fractional form to the subtractive form and can guarantee the convergence to the global optimal solution. Herein, the optimal energy efficiency can be redefined as q * :
q * = Ξ ( λ Q 1 * , λ Q 2 * ) Ξ ( λ Q 1 * , λ Q 2 * ) Ψ ( λ Q 1 * , λ Q 2 * ) Ψ ( λ Q 1 * , λ Q 2 * ) = max { λ Q 1 , λ Q 2 } F Ξ ( λ Q 1 , λ Q 2 ) Ξ ( λ Q 1 , λ Q 2 ) Ψ ( λ Q 1 , λ Q 2 ) Ψ ( λ Q 1 , λ Q 2 ) .
In Equation (23), F is a compact feasible set, with F = λ Q 1 , λ Q 2 ( 22 b ) ( 22 c ) ( 22 d ) ( 22 e ) and λ Q 1 * , λ Q 2 * represents the optimal value of λ Q 1 , λ Q 2 to P 5 . Following the equivalent transformation shown in Equation (23), Theorem 2 can be given as follows.
Theorem 2.
The maximum energy efficiency q * can be achieved if and only if
max { λ Q 1 , λ Q 2 } F Ξ ( λ Q 1 , λ Q 2 ) q * Ψ ( λ Q 1 , λ Q 2 ) = Ξ ( λ Q 1 * , λ Q 2 * ) q * Ψ ( λ Q 1 * , λ Q 2 * ) = 0
with Ξ ( λ Q 1 , λ Q 2 ) 0 , Ψ ( λ Q 1 , λ Q 2 ) > 0 .
Proof. 
Assume that λ Q 1 * , λ Q 2 * F is the optimal precoding strategy for problem P 5 and q * is the optimal energy efficiency. For any feasible solutions, they satisfy the following inequality, Equation (25),
Ξ ( λ Q 1 , λ Q 2 ) Ψ ( λ Q 1 , λ Q 2 ) q * = Ξ ( λ Q 1 * , λ Q 2 * ) Ψ ( λ Q 1 * , λ Q 2 * ) .
Equation (25) can be transformed to the following two inequalities, Equations (26) and (27),
Ξ ( λ Q 1 , λ Q 2 ) q * Ψ ( λ Q 1 , λ Q 2 ) 0 ,
Ξ ( λ Q 1 * , λ Q 2 * ) q * Ψ ( λ Q 1 * , λ Q 2 * ) = 0 .
Equations (26) and (27) imply that maximizing the energy efficiency q * is equivalent to finding the optimization variables λ Q 1 * , λ Q 2 * that satisfy max { λ Q 1 , λ Q 2 } F Ξ ( λ Q 1 , λ Q 2 ) q * Ψ ( λ Q 1 , λ Q 2 ) = 0 . □
Theorem 2 validates the theoretical operability for the Dinkelbach method. It enables the fractional programming problem to be transformed to the subtractive-form via introducing q. The optimal solution of the fractional programming problem can be indirectly obtained by solving the subtractive-form optimization problem. The problem P 6 , which is converted from problem P 5 , is hence obtained and shown as follows
P 6 : max { λ Q 1 , λ Q 2 } F Ξ ( λ Q 1 , λ Q 2 ) q Ψ ( λ Q 1 , λ Q 2 ) .
Based on the Dinkelbach method, q can be iterated from an initial point and the convergence condition proposed in Theorem 2 should be verified in each iteration. In other words, substituting q of each iteration into P 6 verifies whether the objective function of P 6 can be converged to the maximum tolerance. If the algorithm converges, the optimal solution is returned; otherwise, q is updated in the following iteration until the maximum tolerance is satisfied.
(2) Optimal α * for Given λ Q 1 , λ Q 2
Given a feasible λ Q 1 , λ Q 2 , problem P 4 with respect to α can be reduced to P 7
P 7 : max α i = 1 N 1 log 2 1 + α λ i , H λ i , Q 1 α σ 2 2 + σ I D 2 + i = 1 N 2 log 2 1 + λ i , H λ i , Q 2 σ 1 2 i = 1 N 1 λ i , Q 1 + i = 1 N 2 λ i , Q 2 + P c
s . t .   η ( 1 α ) ( i = 1 N 1 λ i , H λ i , Q 1 + σ 2 2 ) i = 1 N 2 λ i , Q 2 + P c 2
i = 1 N 1 log 2 1 + α λ i , H λ i , Q 1 α σ 2 2 + σ I D 2 r ¯ min 2
0 α 1 .
Due to the objective function of P 7 being concave with respect to α , the optimization problem P 7 has a closed-form solution for it. When taking the EE function’s first derivative about α , the inequality EE α > 0 is obtained. This means that EE is a monotonically increasing function with respect to α . Then, the constraint of Equation (29b) can be processed by equivalent transformation and Equation (30) can be obtained as follows
α 1 i = 1 N 2 λ i , Q 2 + P c 2 η ( i = 1 N 1 λ i , H λ i , Q 1 + σ 2 2 ) .
The optimal variable alpha in P 7 must achieve the upper bound 1 i = 1 N 2 λ i , Q 2 + P c 2 η ( i = 1 N 1 λ i , H λ i , Q 1 + σ 2 2 ) constrained by Equation (30). Otherwise, the objective function of P 7 can obtain a larger value via increasing α instead of violating any constraints. Thus, if the optimization variables λ Q 1 , λ Q 2 are maintained as constants, the optimal α can be obtained as
α * = 1 i = 1 N 2 λ i , Q 2 + P c 2 η ( i = 1 N 1 λ i , H λ i , Q 1 + σ 2 2 ) .
Remark 1.
In Equation (31), all parameters in the second term on the right side of the equal sign are positive. This means that the optimal value of α is in 0 , 1 . The optimization variable α cannot be negative or the problem P 7 will be infeasible.
(3) Joint optimization for λ Q 1 , λ Q 2 and α
With the Dinkelbach method, λ Q 1 , λ Q 2 can converge to their optimums by fixing α . Consequently, an alternating algorithm based on the Dinkelbach method can also be proposed to guarantee that both λ Q 1 * , λ Q 2 * and α * converge to their optimums synchronously. The procedures of the proposed alternating algorithm based on the Dinkelbach method are explained using the pseudo-code. Here, m and n denote the iteration numbers of α and λ Q 1 , λ Q 2 , respectively.
In the pseudo-code of the proposed algorithm (Algorithm 1), with an initialized α ( m = 0 ) , the Dinkelbach method is adopted in order to find λ Q 1 ( n ) , λ Q 2 ( n ) . Then, with the known λ Q 1 ( n ) , λ Q 2 ( n ) , α ( m + 1 ) is updated according to Equation (31). The iterations will terminate until the convergence tolerance is met and then the optimal variables can be obtained.
Algorithm 1: Alternating algorithm based on Dinkelbach
1. Initialization
Set maximum tolerance ε 1 , ε 2 and iterations n = 0 , m = 0 ;
Set initial points α ( 0 ) , q = 0 and q p r e = 0 ;
2. Repeat
3.  Repeat
4.   Substitute α ( m ) into Equation (28) and use CVX to calculate λ Q 1 ( n ) , λ Q 2 ( n ) ;
5.   If Ξ ( α ( m ) , λ Q 1 ( n ) , λ Q 2 ( n ) ) q Ψ ( λ Q 1 ( n ) , λ Q 2 ( n ) ) < ε 1
6.    Convergence = true;
7.    Return ( α * , λ Q 1 * , λ Q 2 * ) = ( α ( m ) , λ Q 1 ( n ) , λ Q 2 ( n ) ) and q * = Ξ ( α * , λ Q 1 * , λ Q 2 * ) Ξ ( α * , λ Q 1 * , λ Q 2 * ) Ψ ( λ Q 1 * , λ Q 2 * ) Ψ ( λ Q 1 * , λ Q 2 * ) ;
8.   Else q = Ξ ( α ( m ) , λ Q 1 ( n ) , λ Q 2 ( n ) ) Ξ ( α ( m ) , λ Q 1 ( n ) , λ Q 2 ( n ) ) Ψ ( λ Q 1 ( n ) , λ Q 2 ( n ) ) Ψ ( λ Q 1 ( n ) , λ Q 2 ( n ) ) and n = n + 1 ;
9.    Convergence = false;
10.  End
11. Until Convergence = true
12.  If q * q p r e ε 2
13.    Convergence = true;
14.  Else Calculated α ( m + 1 ) according to Equation (31), m = m + 1 and α ( m ) = α ( m + 1 ) ;
15.    q p r e = q * ;
16.   Convergence = false;
17.  End
18. Until Convergence = true
Convergence and Optimality Analysis:The objective function of P 4 is joint pseudo-concave with λ Q 1 , λ Q 2 for each fixed α . By observing the compact feasible set F = λ Q 1 , λ Q 2 ( 22 b ) ( 22 c ) ( 22 d ) ( 22 e ) , it can conclude that there is no coupling between the optimized variables. Therefore, the set F is a strictly convex feasible set and can ensure the convergence to a Karush–Kuhn–Tucker (KKT) point [32]. Meanwhile, according to the Dinkelbach method, we transform the objective function of P 5 to the subtractive form in P 6 and combine P 6 with the convex feasible set F to obtain a convex optimization problem. So we can get the optimal value of λ Q 1 * , λ Q 2 * by the internal loop of the proposed algorithm.
In the outer loop, the proposed algorithm generates an increasing objective function at each iteration, which means the optimization problem is upper-bounded and hence the convergence of P 7 is guaranteed. According to the derivative that EE α > 0 , it can be concluded that, when fixing the other variables except α , the objective function of P 7 is a monotonically increasing function with respect to α . In the proposed algorithm, for a given λ Q 1 * , λ Q 2 * , α is optimized to increase EE, which means that the value of α will increase with iterations. Therefore, if α is initialized to be extremely small, the proposed algorithm will converge.
Complexity Analysis: Now we will analyze the computational complexity of the proposed alternating algorithm based on Dinkelbach. The computational complexity is inferred based on the optimization problem P 4 . It must be known that the Dinkelbach algorithm has the characteristic of super-linear convergence with K 1 iterations. In addition, the computational complexity can be represented by a polynomial composed of the number of variables and constraints contained in the algorithm’s internal loop due to the optimization problem P 4 being convex, i.e., O ( K 1 ( N 1 + N 2 ) ) . The computational complexity for searching α is O ( 1 ) . Thus, the total complexity of the proposed algorithm is O ( K 2 ( K 1 ( N 1 + N 2 ) + 1 ) ) , where K 2 is the number of iterations within the outer loop.

4. The Multiple EH Sensor Nodes Scenario

The energy efficiency maximization problem in the extended scenario with one source and K EH sensor nodes being used is considered. Similar to the settings in the single node scenario, the source is equipped with N 1 antennas and each sensor node is equipped with N 2 antennas.
The received signals at the k-th EH sensor node can be given as in Equation (32),
y 2 , k = k = 1 K H k F k ϑ 1 , k + n 2 , k = H k F k ϑ 1 , k D e s i r e d s i g n a l + l k K H k F l ϑ 1 , l I n t e r f e r e n c e s i g n a l + n 2 , k A W G N ,
where ϑ 1 , k C N 1 × 1 is the symbol vector for the k-th EH sensor node. In particular, it can be assumed that EH sensor nodes can remove all interference signals before decoding the information of a desired signal [33]. Thus, the achievable sum-rate in bits/s/Hz of K EH sensor nodes can be given as in Equation (33),
R 2 , s e n s o r s = k = 1 K log 2 I N 2 + α H k Q 1 , k H k H ( α σ 2 , k 2 + σ I D , k 2 ) I N 2 ,
where Q 1 , k = F k F k H . Meanwhile, the amount of totally harvested energy of K EH sensor nodes is shown as follows
E E H , tol = η ( 1 α ) T r k = 1 K ( H k Q 1 , k H k H + σ 2 , k 2 I N 2 ) .
When the sensor nodes transmit information to the source by utilizing the harvested energy, the received signals at the source can be deduced as in Equation (35),
y 1 , s o u r c e = k = 1 K H k H D k ϑ 2 , k + n 1 .
According to Equation (35), the achievable data rate of the source is given by Equation (36),
R 1 , s o u r c e = log 2 I N 1 + k = 1 K H k H Q 2 , k H k σ 1 2 I N 1 ,
where Q 2 , k = D k D k H .
Thus, the energy efficiency of the MIMO half-duplex wireless sensor networks with multiple EH sensor nodes can be denoted as
E E m u l = R 2 , s e n s o r s + R 1 , s o u r c e P tol ,
where P tol = P ˜ 1 + P ˜ 2 + P ˜ c , P ˜ 1 = k = 1 K T r ( Q 1 , k ) , P ˜ 2 = k = 1 K T r ( Q 2 , k ) and P ˜ c = P c 1 + K P c 2 . The circuit power consumption of an EH sensor node is assumed to be the same for each other and hence K P c 2 denotes the total circuit power consumption of all sensor nodes.
For the formulation of the energy efficiency maximization problem, this paper adopts the same constraints as initially considered in the single node scenario. Therefore, the optimization problem can be formulated as P8,
P 8 : max α , Q 1 , k , Q 2 , k R 2 , s e n s o r s + R 1 , s o u r c e P tol
s . t . E E H , tol P ˜ 2 + K P c 2
k = 1 K T r ( Q 1 , k ) P ˜ max 1
k = 1 K T r ( Q 2 , k ) P ˜ max 2
R 2 , s e n s o r s r ˜ min 2
R 1 , s o u r c e r ˜ min 1
0 α 1 , Q 1 , k 0 , Q 2 , k 0 ,
where Equation (38b) is the harvested energy constraint, P ˜ max 1 is the upper limit of the source’s transmitting power and P ˜ max 2 is the total available power of all sensor nodes. r ˜ min 2 and r ˜ min 1 are the required rate threshold of K nodes and the source, respectively.
By applying the proposed optimal solution in Section 3.2 to problem P8, the optimization problem can be simplified as P9.
P 9 : max α , λ Q 1 , k , λ Q 2 , k k = 1 K i = 1 N 1 log 2 1 + α λ i , H k λ i , Q 1 , k α σ 2 , k 2 + σ I D , k 2 + i = 1 N 2 log 2 1 + k = 1 K λ i , H k λ i , Q 2 , k σ 1 2 k = 1 K i = 1 N 1 λ i , Q 1 , k + k = 1 K i = 1 N 2 λ i , Q 2 , k + P ˜ c
s . t . η ( 1 α ) k = 1 K i = 1 N 1 λ i , H k λ i , Q 1 , k + σ 2 , k 2 k = 1 K i = 1 N 2 λ i , Q 2 , k + K P c 2
k = 1 K i = 1 N 1 λ i , Q 1 , k P ˜ max 1
k = 1 K i = 1 N 2 λ i , Q 2 , k P ˜ max 2
k = 1 K i = 1 N 1 log 2 1 + α λ i , H k λ i , Q 1 , k α σ 2 , k 2 + σ I D , k 2 r ˜ min 2
i = 1 N 2 log 2 1 + k = 1 K λ i , H k λ i , Q 2 , k σ 1 2 r ˜ min 1
0 α 1 , λ Q 1 , k 0 , λ Q 2 , k 0 ,
where λ Q 1 , k and λ Q 2 , k are the set of eigenvalues of Q 1 , k and Q 2 , k , respectively, i.e., λ Q 1 , k = λ i , Q 1 , k i = 1 , 2 , N 1 , k and λ Q 2 , k = λ i , Q 2 , k i = 1 , 2 , N 2 , k .
The problem P9 can be solved by the alternating algorithm based on Dinkelbach, and the simulation results are described in Section 5.2.

5. Simulation Results

5.1. Single EH Sensor Node Scenario

In this subsection, numerical results are exhibited to show the performances of the proposed alternating algorithm based on the Dinkelbach method. In the simulations, an EE maximization precoding design in a SWIPT-based MIMO half-duplex wireless sensor network is presented by setting N 1 = N 2 , i.e., the number of the equipped antennas at the source is equal to that of the sensor node. The block fading channels modeled as complex Gaussian random variables are adopted in this paper. The elements of the channel coefficient matrix H are generated complying with the distribution C N ( 0 , 1 ) . The maximal transmitting power and the noise power at the source and the sensor node are assumed to be the same, i.e., P ¯ max 1 = P ¯ max 2 = P max and σ 1 2 = σ 2 2 = σ I D 2 = σ 2 = 0.2 [15]. The total circuit power consumption can be set as P c = 3 W . The source and the sensor node have the same requirement on the data rate threshold, i.e., r ¯ min 1 = r ¯ min 2 = r min . The energy conversion efficiency is simplified as η = 1 [10]. In other simulating scenarios, as different η is used, the alternating algorithm maintains the same procedures just as that of η = 1 . All the numerical simulation results are averaged over 500 independent randomly generated channel realizations unless otherwise specified.
H = 0.2772 + 1.1585 i 0.2124 0.2818 i 0.4562 0.0982 i 0.1645 + 0.4383 i 0.4567 0.4329 i 0.3644 0.4480 i 0.0737 + 1.1462 i 1.3267 + 0.4853 i 1.4030 0.0115 i 0.7624 0.7654 i 0.5862 + 0.4177 i 0.0699 0.2621 i 0.2058 0.1429 i 0.3878 + 0.6187 i 0.1906 + 0.1061 i 0.4486 1.6408 i
Figure 2 presents the energy efficiency achieved by the proposed algorithm versus the number of iterations. The results are obtained using five different initial points by setting the maximal transmitting power P max = 8 W and the data rate threshold r min = 1 bps/Hz. The channel coefficient matrix H is generated complying with the distribution C N ( 0 , 1 ) . The source and the sensor node are equipped with the same amount of antennas by setting N 1 = N 2 = 4 . Simulation results in Figure 2 validate that the proposed algorithm enables the energy efficiency to converge to its optimum value 2.7504 bits/Hz/J with no more than three iteration times when iterating from five different points.
In Figure 3, the performance comparisons among the alternating algorithm and the other two baseline schemes, i.e., the non-EH restriction scheme and the non-precoding design for the source scheme, are conducted by setting the other parameters as N 1 = N 2 = 4 and r min = 1 bps/Hz. The non-EH restriction scheme means that both the source and the sensor node are served with the precoding design and are supplied with sufficient energy. It is unnecessary for the sensor node to harvest the RF energy from the source. The received signals at the sensor node are decoded as the data information. In addition, the non-precoding design for the source scheme means that each diagonal element of the source’s eigenvalue matrix is equal to P max P max N 1 N 1 . The simulation results show that the energy efficiency of the non-EH restriction scheme is the highest among all the baseline algorithms. This is because all of the received signals at the sensor node are used to transmit information, then the sensor node can achieve a higher transmission rate at the same energy consumption in the non-EH restriction scheme. However, in the scheme of a non-precoding design for the source, the transmission power is distributed evenly across each antenna of the source. The energy efficiency of this scheme is always the lowest among all baseline algorithms. The energy efficiency without precoding even decreases if P max continues to improve. This is because the transmitting power of the source cannot be properly configured according to the channel states between antenna pairs, so the information rate of the sensor node is lower at the same energy consumption. By contrast, the energy efficiency in the proposed alternating algorithm based on Dinkelbach can rapidly converge to an optimal value. Unlike the non-precoding design scheme, the energy efficiency will not descend with the increase of P max . To sum up, the precoding design scheme we designed can not only achieve a higher energy efficiency than the non-precoding scheme, but also rapidly converges as P max increases.
Figure 4 presents the performance comparison between our proposed alternating algorithm and the non-EH restriction scheme by varying the number of antennas at both the source and the sensor node. In the simulation, the required rate is set as r min = 1 bps/Hz and the influences on energy efficiency as varying the quantity of antennas are compared by gradually increasing P max in the two schemes. It can be intuitively observed from Figure 4 that the more antennas are equipped, the higher the energy efficiency obtained in each scheme. This is because the increase in the number of antennas can enhance the spatial freedom degree and the antennas’ diversity gain. Meanwhile, it can be concluded from Figure 4 that the energy efficiency in the proposed “alternating algorithm based on Dinkelbach” scheme approximates to that of the non-EH restriction scheme at two scenarios with a different number of antennas.
In Figure 5, the energy efficiency performance is studied by changing the required minimum data rate. The antenna numbers in subgraphs (a) and (b) are configured as N 1 = N 2 = 4 and N 1 = N 2 = 8 , respectively. The transmission power constraints of subgraphs (a) and (b) are all set as P max = 8 W . It can be observed from Figure 5 that the highest energy efficiencies shown in the two subgraphs are obtained at two distinct required minimum data rates, i.e., 3 bps/Hz and 6 bps/Hz, respectively. As the required data rate is improved, the energy efficiency improves and then it drops with the rate increment. This is because, by lifting the data rate threshold, the actual transmission power of antennas needs to be increased synchronously. According to the objective function of P 4 shown as the fraction in Equation (20a), the energy efficiency is simultaneously affected by the joint contribution of its numerator and denominator. The variation of transmission power influences the fraction value. When the energy efficiency achieves its maximum, continuing to increase the data rate threshold will cause the energy efficiency to decrease inversely.
In Figure 6, the probability of feasibility under different P max settings is investigated by increasing the required minimum data rate. If the achievable data rate fails the rate threshold, the maximization problem P 4 will be infeasible, i.e., no optimal solutions to P 4 exist. In order to verify the probability feasibility under different data rate thresholds, the number of antennas at transceivers is set as N 1 = N 2 = 4 . When the amount of feasible solutions to problem P 4 is less than half of the total simulation results, the problem is called ergodically infeasible [34]. The result of Figure 6 shows that the probability of feasibility decreases with the increase of the minimum required data rate. The higher the transmission power P max that can be obtained, the larger probability of feasibility achieved.

5.2. Multiple EH Sensor Nodes Scenario

In this subsection, the energy efficiency performance of the proposed algorithm is evaluated in the networks with multiple EH sensor nodes. In the scenario, three EH sensor nodes, i.e., K = 3 are simulated in the scenario and the quantity of the equipped antennas at the source equals that of each sensor node, such that N 1 = N 2 = 4 . The maximal transmitting power and the data rate threshold for both the source and each EH sensor node are given as P ˜ max 1 = P ˜ max 2 = P max and r ˜ min 1 = r ˜ min 2 = 3 bps / Hz , respectively. The other parameters are kept the same as those used in the single EH sensor node scenario.
In Figure 7, the energy efficiency performance is obtained by varying the maximal transmitting power P max . The two baseline algorithms for comparison are designed by assuming each sensor node with indefinite power supply (i.e., non-EH restriction) and by omitting the source precoding design (i.e., non-source precoding), respectively. Distinct from the single EH sensor node scenario, the source eigenvalue matrix in the baseline algorithm of non-source precoding can be obtained as P max K N 1 I N 1 . This is because the source has to allocate power for each sensor node. By observing Figure 7, it can be seen that the energy efficiency of the proposed alternating algorithm based on Dinkelbach outperforms that of the non-precoding design algorithm. As the maximal transmitting power P max ranges from 1 W to 1.8 W, the energy efficiency of the non-source precoding algorithm falls to zero and the value range of P max can be deemed as the outage range, which is due to fact that the evenly distributed source power P max on each antenna cannot support the required sum rate, i.e., the sum rate cannot meet the rate threshold. In contrast, the maximal transmitting power P max at the source can be optimally assigned to each antenna and hence be capable of ensuring the optimal performance. The energy efficiency of the non-EH restriction algorithm is always the highest among all the candidates, since a much higher achievable data rate at sensor nodes can be obtained when all sensor nodes are free from the constraint of power use.
Figure 8 presents the performance comparison between the alternating algorithm and the fixed power-splitting scheme. In the fixed power-splitting scheme, the power splitting factor α is fixed. In the numerical experiment, the power splitting factor was taken as 0.2, 0.4, 0.6 and 0.8, separately. Figure 8 shows that with the increase of maximal transmitting power, the energy efficiency of all schemes exhibits the monotonous non-decreasing tendency. Meanwhile, it is observed from Figure 8 that the alternating algorithm proposed in the paper can always achieve higher energy efficiency. This is because the fixed power-splitting scheme is unable to guarantee the optimized trade-off of power use between information decoding and energy harvesting, thus resulting in lowering the achievable data rate and consuming more energy. When the maximal transmitting power ranges from 1 W to 2 W, the energy efficiency of the fixed power-splitting scheme falls to 0 as α = 0.8 . This is because a greater proportion of source power is spent on information decoding and correspondingly the sensor node cannot harvest enough RF energy to satisfy the data rate constraint. The optimization problem P 9 is essentially infeasible. In the numerical experiment, the proposed alternating algorithm can always achieve better performance than its counterpart, i.e., a fixed power-splitting scheme.

6. Conclusions

In this paper, energy efficiency optimization is investigated in a single-hop MIMO half-duplex wireless sensor network considering a SWIPT scheme. The perfect channel state information is introduced in the energy efficiency maximization problem, where the precoding matrix of the source and the EH sensor node as well as the power splitting factor are jointly designed. The formulated non-convex problem can be equivalently converted by the eigenvalue and singular value decomposition strategy. The alternating algorithm based on Dinkelbach is proposed to solve the converted non-convex optimization problem. The proposed algorithm is validated to achieve the optimality. Simulation results are analyzed to verify the performance of different scenarios for the baseline algorithms. The simulation results evaluate the factors influencing energy efficiency such as the antenna number and data rate threshold. It is shown that the precoding design scheme proposed in the paper can significantly improve energy efficiency. Compared with the non-source precoding scheme, the precoding scheme proposed in the paper shows rapid convergence behavior. These results may provide some insights for better understanding energy efficiency in the multi-antenna SWIPT system. In our future work, the energy efficiency maximization problem will be investigated in the half-duplex SWIPT system with more transceivers.

Author Contributions

Conceptualization and methodology, L.X. and J.-L.W.; investigation and formal analysis, J.L.; validation, L.X.; writing-original draft preparation, L.X. and J.-L.W.; software, J.-L.W.; writing-review and editing, Y.-L.W.; supervision, X.-P.G.; project administration and funding acquisition, L.X.

Funding

This work is partially supported by National Natural Science Foundation of China (NSFC) under Grants 61304131, 61673117, the Natural Science Foundation of Hebei Province under Grant F2016402054, the Scientific Research Plan Project of Hebei Education Department under Grants BJ2014019, ZD2015087, QN2015046, the Research Program of Talent Cultivation Project in Hebei Province under Grant A2016002023.

Acknowledgments

The authors would like to thank Xinping Guan, IEEE Fellow, at Shanghai Jiao Tong University, Yanyan Shen at Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, and Zhezhuang Xu at Fuzhou University for their constructive comments. The authors would also like to acknowledge Na Li, Jin Wang, Lin Zuo, Datian Liu, Yao Ma, Chenyang Xu, and Wanqiang Qin at Hebei University of Engineering, and Yanlong Wang at Beijing University of Posts and Telecommunications for their input into the completion of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ismail, M.; Zhuang, W. Network cooperation for energy saving in green radio communications. IEEE Wirel. Commun. 2011, 18, 76–81. [Google Scholar] [CrossRef]
  2. Varshney, L.R. Transporting information and energy simultaneously. In Proceedings of the 2008 IEEE International Symposium on Information Theory, Toronto, ON, Canada, 6–11 July 2008; pp. 1612–1616. [Google Scholar]
  3. Xiong, K.; Wang, B.; Liu, K.R. Rate-energy region of SWIPT for MIMO broadcasting under nonlinear energy harvesting model. IEEE Trans. Wireless Commun. 2017, 16, 5147–5161. [Google Scholar] [CrossRef]
  4. Zhang, R.; Ho, C.K. MIMO broadcasting for simultaneous wireless information and power transfer. IEEE Trans. Wireless Commun. 2013, 12, 1989–2001. [Google Scholar] [CrossRef]
  5. Li, L.; Wang, H.; Wang, Z.; Paulraj, A. Simultaneous wireless information and power transfer in multi-user interference SISO system. In Proceedings of the IEEE International Conference on Ubiquitous Wireless Broadband (ICUWB), Montreal, QC, Canada, 4–7 October 2015; pp. 1–5. [Google Scholar]
  6. Shi, Q.; Liu, L.; Xu, W.; Zhang, R. Joint transmit beamforming and receive power splitting for MISO SWIPT systems. IEEE Trans. Wireless Commun. 2014, 13, 3269–3280. [Google Scholar] [CrossRef]
  7. Wang, F.; Peng, T.; Huang, Y.; Wang, X. Robust transceiver optimization for power-splitting based downlink MISO SWIPT systems. IEEE Signal Process. Lett. 2015, 22, 1492–1496. [Google Scholar] [CrossRef]
  8. Li, Q.; Zhang, Q.; Qin, J. Beamforming in non-regenerative two-way multi-antenna relay networks for simultaneous wireless information and power transfer. IEEE Trans. Wireless Commun. 2014, 13, 5509–5520. [Google Scholar] [CrossRef]
  9. Benkhelifa, F.; Salem, A.S.; Alouini, M.S. Rate maximization in MIMO decode-and-forward communications with an EH relay and possibly imperfect CSI. IEEE Trans. Commun. 2016, 64, 4534–4549. [Google Scholar] [CrossRef]
  10. Niu, H.; Guo, D.; Huang, Y.; Zhang, B. Robust energy efficiency optimization for secure MIMO SWIPT systems with non-linear EH model. IEEE Commun. Lett. 2017, 21, 2610–2613. [Google Scholar] [CrossRef]
  11. Gong, S.; Xing, C.; Fei, Z.; Kuang, J. Secure communications for SWIPT over MIMO interference channel. In Proceedings of the 2016 IEEE 27th Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC), Valencia, Spain, 4–8 September 2016; pp. 1–5. [Google Scholar]
  12. Jiang, R.; Xiong, K.; Fan, P.; Zhong, S.; Zhong, Z. Optimal beamforming and power splitting design for SWIPT under non-linear energy harvesting model. In Proceedings of the GLOBECOM 2017-2017 IEEE Global Communications Conference, Singapore, 4–8 December 2017; pp. 1–6. [Google Scholar]
  13. Shen, Y.; Huang, X.; Kwak, K.S.; Yang, B.; Wang, S. Subcarrier-pairing-based resource optimization for OFDM wireless powered relay transmissions with time switching scheme. IEEE Trans. Signal Process. 2017, 65, 1130–1145. [Google Scholar] [CrossRef]
  14. Ng, D.W.K.; Lo, E.S.; Schober, R. Energy-efficient resource allocation in multiuser OFDM systems with wireless information and power transfer. In Proceedings of the 2013 IEEE Wireless communications and networking conference (WCNC), Shanghai, China, 7–10 April 2013; pp. 3823–3828. [Google Scholar]
  15. Rostampoor, J.; Razavizadeh, S.M.; Lee, I. Energy efficient precoding design for SWIPT in MIMO two-way relay networks. IEEE Trans. Veh. Technol. 2017, 66, 7888–7896. [Google Scholar] [CrossRef]
  16. Zhou, X.; Li, Q. Energy efficiency for SWIPT in MIMO two-way amplify-and-forward relay networks. IEEE Trans. Veh. Technol. 2018, 67, 4910–4924. [Google Scholar] [CrossRef]
  17. Huang, Y.; Liu, M.; Liu, Y. Energy-efficient SWIPT in IoT distributed antenna systems. IEEE Internet Things J. 2018, 5, 2646–2656. [Google Scholar] [CrossRef]
  18. Wu, Q.; Tao, M.; Ng, D.W.K.; Chen, W.; Schober, R. Energy-efficient resource allocation for wireless powered communication networks. IEEE Trans. Wirel. Commun. 2015, 15, 2312–2327. [Google Scholar] [CrossRef]
  19. Vu, Q.-D.; Tran, L.-N.; Farrell, R.; Hong, E.-K. An efficiency maximization design for SWIPT. IEEE Signal Process. Lett. 2015, 22, 2189–2193. [Google Scholar] [CrossRef]
  20. Shi, Q.; Peng, C.; Xu, W.; Hong, M.; Cai, Y. Energy efficiency optimization for MISO SWIPT systems with zero-forcing beamforming. IEEE Trans. Signal Process. 2015, 64, 842–854. [Google Scholar] [CrossRef]
  21. Xiong, C.; Lu, L.; Li, G.Y. Energy efficiency tradeoff in downlink and uplink TDD OFDMA with simultaneous wireless information and power transfer. In Proceedings of the 2014 IEEE International Conference on Communications (ICC), Sydney, NSW, Australia, 10–14 June 2014; pp. 5383–5388. [Google Scholar]
  22. Gui, L.; Bao, F.; Zhou, X.; Yu, C.; Li, J. Energy-efficient wireless powered secure transmission with cooperative jamming for public transportation. IEEE Trans. Green Commun. Netw. 2019, 1. [Google Scholar] [CrossRef]
  23. Gui, L.; Shi, Y.; Cai, W.; Shu, F.; Zhou, X.; Liu, T. Design of incentive scheme using contract theory in energy-harvesting enabled sensor networks. Phys. Commun. 2018, 28, 166–175. [Google Scholar] [CrossRef]
  24. Boshkovska, E.; Ng, D.W.K.; Zlatanov, N.; Schober, R. Practical non-linear energy harvesting model and resource allocation for SWIPT systems. IEEE Commun. Lett. 2015, 19, 2082–2085. [Google Scholar] [CrossRef]
  25. Guo, W.; Li, J.; Chen, G.; Niu, Y.; Chen, C. A pso-optimized real-time fault-tolerant task allocation algorithm in wireless sensor networks. IEEE Trans. Parall. Distr. Syst. 2014, 26, 3236–3249. [Google Scholar] [CrossRef]
  26. Yao, Y.; Cao, Q.; Vasilakos, A.V. EDAL: An energy-efficient, delay-aware and lifetime-balancing data collection protocol for heterogeneous wireless sensor networks. IEEE/ACM Trans. Netw. 2015, 23, 810–823. [Google Scholar] [CrossRef]
  27. Cheng, H.; Su, Z.; Xiong, N.; Xiao, Y. Energy-efficient node scheduling algorithms for wireless sensor networks using markov random field model. Inf. Sci. 2016, 329, 461–477. [Google Scholar] [CrossRef]
  28. Cheng, H.; Xiong, N.; Yang, L.T.; Jeong, Y.-S. Distributed scheduling algorithms for channel access in TDMA wireless mesh networks. J. Supercomput. 2013, 63, 407–430. [Google Scholar] [CrossRef]
  29. Nguyen, K.-G.; Vu, Q.-D.; Tran, L.-N.; Juntti, M. Energy efficiency fairness for multi-pair wireless-powered relaying systems. IEEE J. Sel. Area. Comm. 2018, 37, 357–373. [Google Scholar] [CrossRef] [Green Version]
  30. Cover, T.M.; Thomas, J.A. Elements of Information Theory; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012. [Google Scholar]
  31. Dinkelbach, W. On nonlinear fractional programming. Manag. Sci. 1967, 13, 492–498. [Google Scholar] [CrossRef]
  32. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  33. Sun, Y.; Ng, D.W.K.; Zhu, J.; Schober, R. Multi-objective optimization for robust power efficient and secure full-duplex wireless communication systems. IEEE Trans. Wirel. Commun. 2016, 15, 5511–5526. [Google Scholar] [CrossRef]
  34. Chen, H.; Gershman, A.B.; Shahbazpanahi, S. Filter-and-forward distributed beamforming in relay networks with frequency selective fading. IEEE Trans. Signal Process. 2009, 58, 1251–1262. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The system model for simultaneous wireless information and power transfer (SWIPT) in a multiple-input multiple-output (MIMO) half-duplex wireless sensor network.
Figure 1. The system model for simultaneous wireless information and power transfer (SWIPT) in a multiple-input multiple-output (MIMO) half-duplex wireless sensor network.
Sensors 19 04923 g001
Figure 2. Energy efficiency of the proposed algorithm with five different initial points versus the number of iterations.
Figure 2. Energy efficiency of the proposed algorithm with five different initial points versus the number of iterations.
Sensors 19 04923 g002
Figure 3. Energy efficiency of the proposed algorithm versus the maximal transmitting power.
Figure 3. Energy efficiency of the proposed algorithm versus the maximal transmitting power.
Sensors 19 04923 g003
Figure 4. Energy efficiency versus maximum transmitting power when varying the number of antennas.
Figure 4. Energy efficiency versus maximum transmitting power when varying the number of antennas.
Sensors 19 04923 g004
Figure 5. The performance comparison of energy efficiency by varying the required minimum data rate and the quantity of equipped antennas.
Figure 5. The performance comparison of energy efficiency by varying the required minimum data rate and the quantity of equipped antennas.
Sensors 19 04923 g005
Figure 6. Probability of feasibility versus the minimum data rate requirement.
Figure 6. Probability of feasibility versus the minimum data rate requirement.
Sensors 19 04923 g006
Figure 7. Energy efficiency versus maximal transmitting power in multiple EH sensor nodes scenario.
Figure 7. Energy efficiency versus maximal transmitting power in multiple EH sensor nodes scenario.
Sensors 19 04923 g007
Figure 8. The comparison of energy efficiency between the proposed algorithm and the fixed power-splitting scheme.
Figure 8. The comparison of energy efficiency between the proposed algorithm and the fixed power-splitting scheme.
Sensors 19 04923 g008

Share and Cite

MDPI and ACS Style

Xue, L.; Wang, J.-L.; Li, J.; Wang, Y.-L.; Guan, X.-P. Precoding Design for Energy Efficiency Maximization in MIMO Half-Duplex Wireless Sensor Networks with SWIPT. Sensors 2019, 19, 4923. https://doi.org/10.3390/s19224923

AMA Style

Xue L, Wang J-L, Li J, Wang Y-L, Guan X-P. Precoding Design for Energy Efficiency Maximization in MIMO Half-Duplex Wireless Sensor Networks with SWIPT. Sensors. 2019; 19(22):4923. https://doi.org/10.3390/s19224923

Chicago/Turabian Style

Xue, Liang, Jin-Long Wang, Jie Li, Yan-Long Wang, and Xin-Ping Guan. 2019. "Precoding Design for Energy Efficiency Maximization in MIMO Half-Duplex Wireless Sensor Networks with SWIPT" Sensors 19, no. 22: 4923. https://doi.org/10.3390/s19224923

APA Style

Xue, L., Wang, J. -L., Li, J., Wang, Y. -L., & Guan, X. -P. (2019). Precoding Design for Energy Efficiency Maximization in MIMO Half-Duplex Wireless Sensor Networks with SWIPT. Sensors, 19(22), 4923. https://doi.org/10.3390/s19224923

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop