Robust Energy Efficient Beamforming in MISOME-SWIPT Systems With Proportional Secrecy Rate

The current research on SWIPT systems focus on three types of receivers: 1) time-switching receiver [15, 30]; 2) power-splitting receiver [23, 24, 30]; and 3) antenna-separating receiver [25, 26, 29, 30, 31]. Based on the receiver setups, significant research efforts have been made to investigate the system power minimization [24], system throughput maximization [30, 31] and harvested energy maximization (HEM) [25, 26, 30, 31]. For examples, the authors in [25] and [26] studied the HEM for the antenna-separating receiver with linear energy harvester and non-linear energy harvester, respectively.

Though the broadcast characteristic of wireless channels enables one-to-many energy delivery, it also increases the probability of legitimate messages to be eavesdropped in the SWIPT systems [32, 33, 34, 35, 36]. Therefore, the combination of physical layer security with SWIPT systems becomes an emerging research topic. For examples, the authors in [32] investigated the secrecy rate maximization (SRM) and HEM in the PtP systems with multiple energy-harvesting nodes (EHNs) co-located with eavesdroppers (EVEs). Combining Charnes-Cooper transformation [37] with one-dimensional search method, the authors in [32] proposed optimal solutions to the SRM and HEM problems. Then, the authors in [33] studied the SRM in single-stream and multiple-stream MIMO-PtP systems with perfect channel state information (CSI). Using a semidefinite relaxation (SDR) technique and an inexact block-coordinate-descent method, the optimal and near-optimal solutions are obtained for the single-stream and multiple-stream cases. The authors in [38] studied the secure communication in the renewable resource powered distributed antenna systems with SWIPT capability, which allows each access point to exchange energy with the central processor. By jointly designing the beamforming vector, artificial noise (AN) covariance matrix and energy exchange variables, the system power consumption is minimized in [38].

Due to the tremendous growth of energy consumption, the wireless operators require new approaches to reduce their energy bills or maximize the energy utilization efficiency. Hence, the EE optimization is one of the most important research issues for the future generation wireless communications [39, 40]. When multiple EVEs exist in the wireless communication systems, a natural extension to the traditional EE optimization becomes secrecy energy efficiency (SEE) optimization. The SEE optimization is more important than traditional EE optimization since the traditional EE optimization does not jointly consider the constructive and detrimental effects of artificial noise in energy efficiency. Although there a number of works considered EE optimization in SWIPT systems [23, 41], the SEE optimization in SWIPT systems has become an active research area recently [42]. Moreover, the SEE optimization induces a better energy utilization efficiency compared with the system throughput maximization or system power minimization. With perfect CSI, the authors in [42] investigated the SEE optimization problem in a MIMO-PtP system with multiple EVEs. Using the Dinkelbach method, they proposed an iterative algorithm to obtain the optimal SEE. However, their proposed algorithm cannot be applied to the SWIPT systems with multiple legitimate users (LUEs). When multiple LUEs exist in the MISO-SWIPT systems, the fairness issue among the LUEs becomes an important topic, which has yet been investigated.

Since the perfect CSI is challenging to obtain in MISO and MIMO systems¹¹ 1 This situation is especially true when the users are passive or roaming to a new place and do not interchange information with the local BST., the authors in [34] studied the system power minimization in the multiuser multiple-input-single-output-multiple-eavesdropper (MISOME)-SWIPT systems with imperfect CSI. Specifically, they leveraged the energy signals and AN to secure the legitimate communication and satisfy the energy requirement. Then, the authors extended the work in [34] and jointly considered the transmission power, RF-EH efficiency and interference power leakage as a multiobjective optimization problem with imperfect CSI [35]. In [36], the authors studied the impact of non-linear energy harvester on system power consumption in a multiuser MISOME-SWIPT system. To the authors’ best knowledge, the robust SEE optimization problem in the multiuser MISOME-SWIPT systems with proportional-secrecy-rate (PSR) constraints has not been reported.

In this work, we consider a multiuser MISOME-SWIPT system where a BST transmits information signals and jamming signals to the LUEs and EVEs, respectively. Meanwhile, the EHNs scavenge energy from the information signals and jamming signals. Different from [42], we study the SEE optimization via joint design of beamforming vector and AN covariance matrix under the imperfect CSI of EVEs and EHNs. Moreover, we include the RF-EH constraints and PSR constraints to guarantee the harvested energy at the EHNs and ensure the fairness among LUEs as in [43, 44], respectively. The formulated SEE optimization problem is more complicated than the traditional EE optimization problem [23, 41] due to the non-convex PSR constraints and non-convex SEE function. In order to develop tractable algorithms, we first simplify the SEE optimization problem via exploiting the problem structure. Based on the introduced parameter, we leverage the SDR technique coupled with a one-dimensional search method to address the simplified SEE optimization problem. The major contributions of this work can be summarized as follows.

• •

  We formulate the SEE maximization problem for a multiuser MISOME-SWIPT system via the joint design of beamforming vector and AN covariance matrix. In order to guarantee the fairness among LUEs, we explicitly confine the secrecy rates of LUEs to a predefined ratio such that the secrecy rates of LUEs are proportional. Since the EVEs and EHNs are either roaming users or passive devices, they may not exchange CSI as frequently as that of LUEs. Hence, we consider the bounded-error CSI of EVEs and EHNs.

• •

  The formulated SEE maximization problem in the MISOME-SWIPT system is non-convex in SEE function and PSR constraints. Therefore, it is difficult and ineffective to solve the non-convex SEE maximization problem via standard methods. In order to develop effective algorithms, we exploit the structure of non-convex PSR constraints such that a tractable form of SEE maximization problem is obtained where the strong duality holds. We theoretically prove that: 1) the non-convex PSR constraints can be relaxed into a finite set of convex constraints; and 2) the relaxations are tight. Based on these facts, we propose a two-stage algorithm such that an optimal solution is obtained via solving a series of convex optimization problems and a one-dimensional search method.

• •

  In order to trade the performance for computational complexity, we propose two suboptimal algorithms based on several heuristic beamforming techniques, e.g., maximal ratio transmission (MRT) and zero-forcing beamforming (ZFBF). The computational complexity of the proposed optimal and suboptimal algorithms are quantitatively compared.

Numerical results are used to verify the performance of the proposed optimal and suboptimal algorithms.

The remainder of this paper is organized as follows. The system model and the SEE maximization problem are presented with imperfect CSI of EVEs and EHNs in Section II. The optimal and suboptimal algorithms are proposed in Section III and Section IV, respectively. Numerical results are used to verify the effectiveness of proposed solutions and provide several practical insights in Section V. Finally, Section VI concludes the work.

Notations: Vectors and matrices are shown in bold lowercase letters and bold uppercase letters, respectively. $\mathbb{C}^{N\times M}$ denotes the $N\times M$ dimension complex-value matrices. $\left\|\cdot\right\|_{\mbox{\tiny F}}$ and $\left|\cdot\right|$ are the Frobenius norm and absolute value, respectively. $\sim$ stands for “distributed as”. $\bm{I}_{N}$ and $\bm{0}_{N\times M}$ denote, respectively, an $N$ dimensional identity matrix and a zero matrix with $N$ rows and $M$ columns. The expectation of a random variable is denoted as $\mathds{E}\left[\cdot\right]$. $\mbox{vec}\left[\bm{W}\right]$ obtains a vector by stacking the columns of $\bm{W}$ under the other. $\mbox{diag}\{\bm{W}_{1},\bm{W}_{2},\ldots,\bm{W}_{N}\}$ denotes a diagonal matrix with the diagonal elements given by the matrices $\bm{W}_{1},\bm{W}_{2},\ldots,\bm{W}_{N}$. $\left\{\bm{w}_{n}\right\}_{n\in\mathcal{N}}$ represents the set made of $\bm{w}_{n}$, $n\in\mathcal{N}$. For a square matrix $\bm{W}$, $\bm{W}^{\mbox{\tiny H}}$ and $\mbox{Tr}\left(\bm{W}\right)$ denote its conjugate transpose and trace, respectively. $\bm{W}\succeq\bm{0}$ and $\bm{W}\succ\bm{0}$ respectively denote that $\bm{W}$ is a positive semidefinite and $\bm{W}$ is a positive definite matrix.

We consider downlink transmission of a MISOME-SWIPT system, which consists of one single BST, a set of LUEs, a set of EVEs and a set of EHNs. Let $\mathcal{N}=\left\{1,2,\ldots,N\right\}$, $\mathcal{M}=\left\{1,2,\ldots,M\right\}$ and $\mathcal{I}=\left\{1,2,\ldots,I\right\}$ denote the set of LUEs, set of EVEs and set of EHNs, respectively. The BST is equipped with $N_{t}$ antennas for simultaneous transmission of information and artificial noise over the same single carrier-frequency band. Therefore, the LUEs and EVEs decode information via beamforming. Each LUE is equipped with one antenna to receive the legitimate information. Each EVE, which is also equipped with one single antenna, can be roaming user from other wireless communication systems and searching for the services from the local BST²² 2 This is due to the fact that the malicious roaming users have the potential to intercept signals intended for LUEs under the coverage of the same BST.. Therefore, we leverage the AN to interrupt EVEs and guarantee RF-EH requirement of EHNs. Here, the EHNs are equipped with one antenna and can be the passive sensors to monitor the environment status of the MISOME-SWIPT system.

At the BST, the information-bearing signal for the $n$-th LUE is given as $s_{n}$ with $\mathds{E}\left[|s_{n}|^{2}\right]=1$. Therefore, the transmission signal for LUEs is denoted by

where $\bm{w}_{n}\in\mathbb{C}^{N_{t}\times 1}$ is the beamforming vector for the $n$-th LUE. In order to secure the BST-LUE links, the BST needs to transmit an AN vector $\bm{q}$ to reduce information leakage to EVEs. Hence, the transmission signal at the BST is given as

where the AN vector $\bm{q}\in\mathbb{C}^{N_{t}\times 1}$ is modeled as a complex Gaussian random vector with mean zero and covariance $\bm{Q}\in\mathbb{C}^{N_{t}\times N_{t}}$. In particular, it is assumed that the information-bearing signals $\left\{s_{n}\right\}_{n\in\cal N}$ and the AN vector $\bm{q}$ are statistically independent.

The frame-based frequency non-selective fading channels with unit duration for each frame is considered; therefore, the words “energy” and “power” can be used interchangeably. In each frame, the BST broadcasts information-bearing signals and jamming signals over the same single carrier-frequency band. Considering the Rayleigh fading channels, the channel-coefficient vector for the $n$-th BST-LUE link $\bm{h}_{n}\in\mathbb{C}^{N_{t}\times 1}$ is circularly symmetric complex Gaussian (CSCG) distributed as $\bm{h}_{n}\sim{\cal CN}\left(\bm{0},\Omega_{u,n}^{-1}\bm{I}_{N_{t}}\right)$. Here, $\Omega_{u,n}^{-1}$ captures the joint effect of link distance and carrier frequency for the $n$-th BST-LUE link [46, 47]. As a result, the received signals of the $n$-th LUE is expressed as

where $z_{u,n}\sim{\cal CN}\left(0,\sigma^{2}_{u,n}\right)$ is the additive white Gaussian noise (AWGN) with mean zero and variance $\sigma^{2}_{u,n}$ of the $n$-th LUE.

Based on the received signals in (3), the data rate of the $n$-th LUE is denoted as

where

with $\bm{H}_{n}\triangleq\bm{h}_{n}\bm{h}_{n}^{\mbox{\tiny H}}$ and $\bm{W}_{n}\triangleq\bm{w}_{n}\bm{w}_{n}^{\mbox{\tiny H}}$, $n\in\cal N$. In addition, the rank of each beamforming matrix $\bm{W}_{n}$ is upper bounded as $\mbox{Rank}\left(\bm{W}_{n}\right)\leq 1$. The term ${\cal BW}$ denotes the bandwidth of the system.

Besides, the received signal at the $m$-th EVE is denoted as

where $z_{e,m}\sim{\cal CN}\left(0,\sigma^{2}_{e,m}\right)$ is the AWGN with mean zero and variance $\sigma^{2}_{e,m}$ at the $m$-th EVE, $m\in\cal M$; $\bm{g}_{e,m}\in\mathbb{C}^{N_{t}\times 1}$ is the channel-coefficient vector for the $m$-th BST-EVE link.

Based on the signal model in (6), the information leakage to the $m$-th EVE of the $n$-th LUE by treating the interference as noise is denoted as [46]

where

where $\bm{G}_{e,m}\triangleq\bm{g}_{e,m}\bm{g}_{e,m}^{\mbox{\tiny H}}$.

Denote the channel-coefficient vector of the $i$-th BST-EHN link as $\bm{g}_{h,i}\in\mathbb{C}^{N_{t}\times 1}$, the received baseband signal of the $i$-th EHN is denoted as

where $z_{h,i}\sim{\cal CN}\left(0,\sigma_{h,i}^{2}\right)$ is the AWGN at the $i$-th EHN, and $\bm{g}_{h,i}$ is associated with the pathloss and carrier frequency of the $i$-th BST-LUE link.

Ignoring the energy that can be harvested from the AWGN, the amount of harvested energy of the $i$-th EHN is denoted as

where $\bm{G}_{h,i}\triangleq\bm{g}_{h,i}\bm{g}_{h,i}^{\mbox{\tiny H}}$, and $\xi_{h,i}\in\left(0,1\right]$ is the energy conversion efficiency.

Since we consider strongly secure communication, the information-leakage rate of the $n$-th BST-LUE link is confined to be smaller than or equal to the rate of auxiliary information $R_{e,n}^{\mbox{\tiny REQ}}$. Hence, the secrecy rate of the $n$-th BST-LUE link is written as [45, 48]

with

where $\left(\cdot\right)^{+}=\max\left(\cdot,0\right)$.

Furthermore, the power consumption of BST in the MISOME-SWIPT system is comprised of the power of downlink beamforming, power of AN and power of circuit. Therefore, the expression of power consumption is given as

where the constant circuit power consumption $P^{\mbox{\tiny CIR}}$ is given as $P^{\mbox{\tiny CIR}}=P^{\mbox{\tiny SP}}\left(0.87+0.1N_{t}+0.03N_{t}^{2}\right)$ since the circuit power consumption is positively correlated to the number of antennas [49]. Here, the term $P^{\mbox{\tiny SP}}$ denotes the baseband processing power consumption.

In this work, the SEE is defined as the ratio of the sum secrecy rates of BST-LUE links over the power consumption of BST. Based on (11) and (13), the SEE of MISOME-SWIPT system is given as

whose unit is Nats/joule.

We assume that the MISOME-SWIPT system operates in time-division duplex mode. At the beginning of each frame, each LUE sends a pilot signal to the BST. After receiving the pilot signals, the BST estimates the channels associated with each LUE that sends the pilot signal. Since each LUE facilitates each uplink transmission with pilot signals, the BST can exploit the uplink reciprocity to periodically update the downlink CSI for BST-LUE links. Therefore, we assume that the BST can perfectly estimate the downlink CSI for the BST-LUE links. However, the EVEs and EHNs can be the roaming users and/or low-power nodes and do not exchange the pilot signals with the BST as frequently as the LUEs. The obtained CSI for the BST-EVE links can be imperfect subject to channel estimation errors and/or quantization errors. Hence, we are motivated to use the bounded-error CSI model to formulate the CSI for the BST-EVE and BST-EHNs links as [50, 51]

where $\bar{\bm{g}}_{e,m}\in\mathbb{C}^{N_{t}\times 1}$ and $\Delta\bm{g}_{e,m}\in\mathbb{C}^{N_{t}\times 1}$ are, respectively, the estimated channel-coefficient vector and channel uncertainty of the $m$-th BST-EVE link. Similarly, $\bm{g}_{h,i}$ and $\Delta\bm{g}_{h,i}$ are the estimated channel-coefficient vector and channel uncertainty vector of the $i$-th EHN. The positive constants $\Theta_{e,m}$ and $\Theta_{h,i}$ denote the radii of uncertainty region of $\bm{g}_{e,m}$ and $\bm{g}_{h,i}$. Here, the vectors $\bar{\bm{g}}_{e,m}\sim{\cal CN}\left(\bm{0},\Omega_{e,m}^{-1}\bm{I}_{N_{t}}\right)$ and $\bar{\bm{g}}_{h,i}\sim{\cal CN}\left(\bm{0},\Omega_{h,i}^{-1}\bm{I}_{N_{t}}\right)$ are CSCG distributed having the terms $\Omega_{e,m}$ and $\Omega_{h,i}$ that take into account of the joint effect of pathloss and carrier frequency of the $m$-th BST-EVE link. In addition, the channel uncertainty vectors $\Delta\bm{g}_{e,m}$ and $\Delta{\bm{g}}_{h,i}$ capture the joint effect of estimation errors and time-varying characteristics of wireless channels.

Our objective is to maximize SEE via the joint design of beamforming vector and AN covariance matrix in the MISOME-SWIPT system. At the beginning of each frame, the BST performs beamforming and jamming in a centralized way. The SEE maximization problem is formulated as

where the constraints in (17b) are used to guarantee the proportional fairness on the secrecy rates of LUEs with $\sum\nolimits_{n=1}^{N}\varphi_{n}=1$ and $\varphi_{n}\geq 0$ [43, 44]; and the information-leakage constraints and the harvested-power constraints are respectively defined in (17c) and (17d). Here, the term $P_{h,i}^{\mbox{\tiny REQ}}$ in (17d) denotes the energy requirement of the $i$-th EHN; and the constant $P^{\max}$ in the constraints in (17e) is the power budget of BST due to the circuit limitation.