Next Article in Journal
Defect Detection and Segmentation Framework for Remote Field Eddy Current Sensor Data
Next Article in Special Issue
A Strategic Bargaining Game for a Spectrum Sharing Scheme in Cognitive Radio-Based Heterogeneous Wireless Sensor Networks
Previous Article in Journal
Detecting Steps Walking at very Low Speeds Combining Outlier Detection, Transition Matrices and Autoencoders from Acceleration Patterns
Previous Article in Special Issue
Medium Access Control Protocols for Cognitive Radio Ad Hoc Networks: A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Weighted Sum Harvested Energy Maximization for SWIPT Cognitive Radio Networks Based on Particle Swarm Optimization

1
School of Electrical and Computer Engineering, University of Ulsan, Ulsan 680-749, Korea
2
Faculty of Physics, University of Education, Hue University, 34 Le Loi Str., Hue City 530000, Vietnam
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(10), 2275; https://doi.org/10.3390/s17102275
Submission received: 31 July 2017 / Revised: 27 September 2017 / Accepted: 2 October 2017 / Published: 6 October 2017
(This article belongs to the Special Issue Cognitive Radio Sensing and Sensor Networks)

Abstract

:
In this paper, we consider multiuser simultaneous wireless information and power transfer (SWIPT) for cognitive radio systems where a secondary transmitter (ST) with an antenna array provides information and energy to multiple single-antenna secondary receivers (SRs) equipped with a power splitting (PS) receiving scheme when multiple primary users (PUs) exist. The main objective of the paper is to maximize weighted sum harvested energy for SRs while satisfying their minimum required signal-to-interference-plus-noise ratio (SINR), the limited transmission power at the ST, and the interference threshold of each PU. For the perfect channel state information (CSI), the optimal beamforming vectors and PS ratios are achieved by the proposed PSO-SDR in which semidefinite relaxation (SDR) and particle swarm optimization (PSO) methods are jointly combined. We prove that SDR always has a rank-1 solution, and is indeed tight. For the imperfect CSI with bounded channel vector errors, the upper bound of weighted sum harvested energy (WSHE) is also obtained through the S-Procedure. Finally, simulation results demonstrate that the proposed PSO-SDR has fast convergence and better performance as compared to the other baseline schemes.

1. Introduction

Radio frequency energy harvesting is emerging as an active research area in the fields of both academics and industry, due to the fact that it can solve the bottle-neck of battery-powered wireless devices. It is especially important and useful in wireless sensor networks since sensor nodes have a limited amount of energy. In this new paradigm, the RF can bring both data and power from transmitters to receivers [1,2,3,4,5,6,7,8], referred to as “simultaneous wireless information and power transfer (SWIPT)”. Thus, new pre-coding techniques for transmitters and optimizing design techniques for receivers need to be investigated. The receiver operates in two modes where it switches between information decoding (ID) and energy harvesting (EH), i.e., time-switching (TS) mode, or shares the incoming signal into EH and ID parts, i.e., power-splitting (PS) [1]. Joint power-splitting SWIPT and beamforming are applied to some conventional networks such as multiuser multi-input single-output (MISO) systems [3,4] where transmission energy and energy efficiency are optimized, respectively. These research issues are also considered in systems with multi-antenna transmitters and receivers [5], cooperative networks [6], and interference channels [7], as well as robust secure transmission [8]. In addition, the authors in [9,10,11] studied the transmit power minimization problem in multiuser SWIPT MISO systems under the imperfect channel side information (CSI) with bounded channel vector errors and stochastic channel vector errors. As an important criterion, the sum harvested energy maximization in multiuser power-splitting SWIPT MISO system is investigated under the perfect channel side information of channels from transmitters to receivers [12]. However, the solution of [12] based on successive second-order cone programming cannot be applied to cognitive radio networks where interference threshold constraints of primary users should be considered with the case of imperfect CSI.
Similar to the works introduced by Goldsmith et al. [13], the applications of RF-powered techniques to harvest energy and transfer data in cognitive radio networks (CRNs) were summarized by Mohjazi et al. [14]. Other researchers studied SWIPT for different scenarios in order to provide energy to receivers while ensuring quality of service. In [15], the secondary network exploited both spectrum and energy in primary networks while assisting primary data transmission. In [16], Ng et al. studied one secondary data link in the presence of multiple energy harvesting receivers and primary users. Yang et al. in [17] and Lee et al. in [18] also investigated SWIPT in cognitive relay and cognitive wirelessly powered networks, respectively. However, CRNs with power-splitting SWIPT and beamforming designs in multiuser scenarios have not been well studied so far. The studies on power-splitting SWIPT in CRNs have many potential applications, such as wirelessly powered cognitive sensor networks and cognitive cellular networks, where users need to both receive information and obtain energy. In the our prior work [19], we considered the SWIPT cognitive radio network scenarios, and designed the system by considering the following important two criteria: one is “max–min harvested energy of cognitive users” and the other is “the worst-user trade-off between harvested energy of cognitive users and interference power of primary users” under the perfect CSI. Another important criterion in SWIPT cognitive radio networks is “weighted sum harvested energy maximization” under both perfect CSI and imperfect CSI cases, which have not been fully investigated yet.
To the best of our knowledge, this is the first work that investigates a multiuser power-splitting SWIPT for CRN in which one secondary transmitter (ST) equipped with a multi-antenna will transmit information and energy to multiple single-antenna secondary receivers (SRs) that have a PS structure, in the existence of multiple single-antenna primary users (PUs). The primary goal of this paper is to maximize the weighted sum harvested energy (WSHE) of all SRs by jointly optimizing the transmit beamforming vectors and the PS ratios while satisfying the minimum requirement of each SR’s signal-to-interference-plus-noise ratio (SINR), the ST’s limited transmission power, and each PU’s specified interference. It is noteworthy that the harvested energy received by an SR can be controlled by adjusting the weighted factors at the ST. The user-centric energy criteria such as max–min fairness harvested energy or the trade-off between EH and interference are studied in our prior work [19] where the more energy is transmitted over the poor channel to combat the channel attenuation for the worst user. Therefore, the worst user is always guaranteed in both information and energy. Unlike [19], in this paper the system-centric WSHE criterion is investigated where the total harvested energy of system is maximized regardless of the harvested energy of worst user. Thus, more energy is transmitted over the better channel, while the information rate is still guaranteed for every cognitive user. In addition, in [19], we only studied the problems for the ideal case with the perfect CSI. Unlike [19], we in this paper also study the WSHE maximization problem in the practical scenario under the imperfect CSI with errors bound. It is noted that this realistic case frequently occurs in wireless cognitive sensor networks with low complexity sensor nodes.
More specifically, in this paper we solve two completely novel research issues that in prior works have not been studied. The first one is WSHE maximization in multi-user SWIPT cognitive radio networks under the perfect CSI. The other issue is consideration of the WSHE problem under the imperfect CSI. The main contributions of this paper can be summarized as follows.
  • For the perfect CSI, we formulate the non-convex optimization problem for the WSHE of all secondary receivers including constraints to limited transmission power at the ST, minimum required SINR at each SR, and interference threshold at each PU. The objective function is non-convex due to the coupled design variables of both the transmit beamforming vectors at the ST and the PS ratios at the SRs. Therefore, we solve it in two steps based on two variable groups of beamforming vectors and PS ratios. In the first step, we fix the PS ratios and obtain the optimal beamforming vectors by applying a semidefinite relaxation (SDR) technique. Interestingly, we can show that SDR is tight for our problem. In the second step, we propose an algorithm based on particle swarm optimization (PSO) to find the approximate optimal PS ratios. The baseline schemes are considered for performance comparison, in which the zero-forcing beamforming vectors at the ST and the equal power splitting ratios at the SRs are applied.
  • For the imperfect CSI, the S-Procedure and SDR technique are applied to recast the robust WSHE maximization problem as a semidefinite programming (SDP) problem. Then, the similar PSO-based method is used to solve the formulated problem. However, we only obtain the upper bound value of WSHE since the optimal solutions are unlikely to satisfy the rank-1 constraints.
  • Finally, simulation results show that the proposed PSO-SDR has fast convergence and better performance compared to the other baseline schemes. In addition, the the proposed PSO-SDR converges to the optimal value, achieved by the brute-force search (BFS) method, while obtaining significantly lower computational complexity.
The remainder of the paper is organized as follows. The system description is presented in Section 2. The weighted sum harvested energy maximization problem in the perfect CSI case is solved in Section 3 by the joint SDR technique and PSO method. Section 4 presents the solution of the robust WSHE maximization problem in the imperfect CSI case. The simulation results are provided in Section 5, followed by conclusions in Section 6.
Notation 1.
Vectors and matrices are indicated by boldface lower case and capital letters, respectively. X * and X H represent the conjugate and Hermitian transpose of matrix X , respectively. We use Tr X , and rank X to indicate trace and rank of matrix X , respectively. X 0 and X 0 indicate that matrix X is positive semidefinite or positive definite, respectively. The Euclidean norm of a complex vector and the absolute value of a complex scalar are represented by · and · , respectively. The space of m × n matrices with complex entries is denoted by C m × n . H N denotes the space of N × N Hermitian matrices. I denotes the identity matrix with appropriate size. C N μ , σ 2 represents the distribution of a circularly symmetric complex Gaussian (CSCG) random variable with mean μ and variance σ 2 , and ’∼’ means “distributed as". Moreover, we summarize the important abbreviations and explain their meanings in Table 1.

2. System Description

Figure 1a shows a considered MISO cognitive radio networks with SWIPT capability that includes one ST, M SRs, and L PUs. M information messages are simultaneously sent by the ST equipped with N antennas to M SRs equipped with a single antenna at the same frequency as the PU’s communication when L single-antenna PUs exist. We denote h i C N × 1 , i 1 , . . . , M and g l C N × 1 , l 1 , . . . , L as the baseband equivalent channels from the ST to S R i , and P U l , respectively. The ST is supposed to obtain these channel vectors perfectly with their elements as independent and identically distributed (i.i.d) circularly symmetric complex Gaussian (CSCG) variables.
The secondary transmitter can practically achieve channel side information (CSI) vectors from the ST to the SRs by utilizing the conventional channel estimation methods. First, the pilots are transmitted by the ST, and then the ST obtains feedback on the results of channel estimation from the SRs. In addition, when there is cooperation between licensed and cognitive systems, the channel estimation can be sent from the PUs directly to the ST via the pilot signals from the cognitive network or indirectly through a frequency manager [23]. In an other way, the ST can achieve link information by periodically detecting the signal sent by the PUs since the licensed network utilizes time division duplex (TDD) for communication. However, it is difficult for the ST to obtain perfect instantaneous CSI for both cross and direct links. Thus, we also consider the harvested energy of the SRs for imperfect CSI case in the paper.
The ST sends the communication signal, which can be represented as follows
x = i = 1 M w i s i ,
where s i C is the symbol which carries information addressed to S R i with E s i 2 = 1 , and the corresponding the precoding beamforming vector w i C N × 1 . The average transmission power is calculated as P = i = 1 M w i 2 . The over-the-air transmit power is limited by the maximum transmit power P max [24,25]. Although the power amplifier efficiency and the constant circuit power consumption of the secondary transmitter accounting for antenna circuits, transmit filter, mixer, frequency synthesizer, and digital-to-analog converter, etc., [24,25] are important parameters in considering energy system consumption, these factors will not affect our main objective maximizing the harvested energy at secondary users under the limited transmit power. Therefore, we neglect these parameters when evaluating the sum harvested energy problem. The baseband signal received at S R i is represented as follows
r i = h i H w i s i + h i H j = 1 , j i M w j s j + n i , i ,
where n i C N 0 , σ n 2 is is the antenna noise added at S R i .
In the system, the SRs equipped with the PS receive structure shown in Figure 1b can simultaneously decode information and harvest energy from incoming signals. Under this PS scheme, S R i splits the incoming signal into two streams, where one stream with power ratio θ i 0 , 1 is utilized for ID, and the other with power ratio ( 1 θ i ) is utilized for EH. Figure 1b shows that the incoming signal for ID is r i , ID = θ i r i + v i where v i C N 0 , σ v 2 is the circuit noise appended by the ID of S R i . Thus, the SINR at S R i is calculated as
SINR i = θ i h i H w i 2 θ i j = 1 , j i M h i H w j 2 + σ n 2 + σ v 2 , i .
Moreover, the incoming signal for EH of S R i is r i , EH = 1 θ i r i . The initial energy can exist at the beginning of each time block and we assume that the size of battery is large enough to store both initial energy and harvested energy. We here focus on jointly optimizing precoding beamforming vectors at the transmitter and power splitting ratios at the secondary receivers so that the SRs can harvest maximum energy.
Then, the energy harvested by the EH of S R i is computed as follows
E H i = η i 1 θ i j = 1 M h i H w j 2 + σ n 2 , i ,
where the efficiency of S R i ’s energy harvester is denoted by η i 0 , 1 . Furthermore, the power of interference at P U l caused by the ST is computed as
I T l = i = 1 M g l H w i 2 , l .

3. Problem Formulation and Solution

In this paper, the main objective is to maximize the WSHE of all SRs by jointly finding the received PS ratios θ i and the beamforming vectors w i subject to the required SINR at each SR, the maximum transmit power at the ST, and the interference threshold for each PU. Therefore, the WSHE problem can be expressed as follows
max w i , θ i i = 1 M λ i η i 1 θ i j = 1 M h i H w j 2 + σ n 2
s . t . i = 1 M w i 2 P max
θ i h i H w i 2 θ i j = 1 , j i M h i H w j 2 + σ n 2 + σ v 2 a i , i , i
i = 1 M g l H w i 2 I t , l
0 < θ i < 1 , i ,
where we set i = 1 M λ i = 1 , and the weight factor λ i > 0 emphasizes the different priority for harvesting energy at S R i . Constraint Equation (6b) corresponds to the power constraint of the ST. Constraint Equation (6c) is corresponds to the SINR constraint such that SIN R i at S R i is larger than the minimum required SINR, a i . Moreover, constraint Equation (6d) shows that the interference power at which the ST interferes with P U l should be lower than I t in the underlay cognitive mode. Since the problem is non-convex due to the objective function and constraint Equation (6c), it is quite challenging to solve the WSHE problem directly. Therefore, in the next sub-section we propose an algorithm to solve the WSHE problem based on semidefinite relaxation and particle swarm optimization.

3.1. SDR Approach with Fixed PS Ratios

The idea for solving problem (6) is to separate the variables into two groups including the precoding vectors w i and the power splitting ratios θ i . To do this, we transform the objective function in Equation (6a) into:
max 0 < θ i < 1 , i max w i i = 1 M λ i E H i .
After that, let us consider problem (7) with fixed PS ratios θ i and only beamforming vector variables as follows:
max w i i = 1 M λ i E H i
s . t . Equations ( 6 a ) , ( 6 c ) and ( 6 d ) .
Using the semidefinite relaxation method [20], we convert Equation (7) into standard semidefinite programming (SDP) [26] which can be effectively resolved via general-purpose numerical solver like CVX [21]. For this, let us denote W i = w i w i H , H i = h i h i H , i and G l = g l g l H , l . Based on the formulas u = Tr ( u ) , u 2 = u H u , and Tr ( UV ) = Tr VU , the following results are obtained as follows:
w i 2 = w i H w i = Tr w i H w i = Tr W i
h i H w j 2 = Tr w j H h i h i H w j = Tr H i W j .
Similarly, g l H w i 2 = Tr G l W i . We also have the property as W i = w i w i H W i 0 and rank W i = 1 . Therefore, the problem (7) can be recast as follows:
min W i i = 1 M λ i η i 1 θ i j = 1 M Tr H i W j + σ n 2
s . t . i = 1 M Tr W i P max 0
Tr H i W i a i + j i M Tr H i W j + σ n 2 + σ v 2 θ i 0 , i
i = 1 M Tr G l W i I t 0 , l
W i 0 , i
rank W i = 1 , i .
The above optimization problem is non-convex due to constraint (8f). According to the SDR technique, we remove constraint (8f) to make problem (8) a standard semidefinite programming (SDP) problem. Here, we define problem (8) without constraint (8f) as the problem (8)-SDR. Then, we solve the problem (8)-SDR via a numerical solver for the convex optimization problem. Note that we refer to Remark 2 for convex problem with respect to complex-valued variables. However, the general SDP problem does not always give the rank-1 solutions and there is also no general method to prove the rank-1 solutions. In [27], the uplink-downlink duality is exploited to obtain the rank-1 solution for the SDP problem with M separable matrix variables and M linear constraints. The idea is to convert the downlink problem into the virtual uplink problem and then construct the rank-1 solution. However, this method can not be applied to our problem due to the different objective function of weighted sum harvested energy and the appearance of the constraints of total transmit power limit and interference threshold at primary users. In [28], in particular, the rank-1 solution of the SDP problem is always obtained whenever M L + 2 with the number of matrix variables M and the number of constraints L. However, this result can not be applied to our proposed SDP problem since the number of constraints and the number of variables in problem (8)-SDR do not satisfy the above rank-1 conditions whenever the number of primary users is larger than 1. Interestingly, the optimal solution to problem (8)-SDR can be proven as rank-1. Therefore, it is also the optimal solution of the original problem (8).
L W i , α , β i , γ l , E i = i = 1 M λ i η i 1 θ i j = 1 M Tr H i W j + σ n 2 + α i = 1 M Tr W i P max + i = 1 M β i 1 1 a i a i Tr H i W i + j i M Tr H i W j + σ n 2 + σ v 2 θ i + l = 1 L γ l i = 1 M Tr G l W i I t i = 1 M Tr E i W i .
Lemma 1.
The optimal solution, W i , of problem (8)-SDR is rank-1, i.e., rank W i = 1 , i = 1 , . . . , M .
Proof. 
The Karush–Kuhn–Tucker (KKT) optimality conditions are applied to the proof of Lemma 1. The Lagrangian function of problem (8)-SDR is written as seen in (9) where α 0 , β i 0 , i , γ l 0 , l , and E i 0 , i are the dual variables associated with constraints (8b), (8c), (8d), and (8e), respectively. Since W i 0 , H i = h i h i H , and G l = g l g l H , we derive that Tr W i , Tr H j W i and Tr G l W i are real-valued for all i , j , l . Therefore, the objective function (8a), the constraint functions in (8b), (8c), (8d) and then Lagrangian (9) are real-valued functions with complex-valued Hermitian matrix variables. At the optimal points, the gradient of Lagrangian (9) must vanish at all independent real variables, i.e., all real and imaginary variables of complex-valued matrix variables W i . This is equivalent to the gradient of real-valued Lagrangian vanishing at W i according to Theorem 2 in ([29], Section 4). ☐
Following Lemma 1 in ([30], Section 4B), we obtain the gradient of Tr AX as X Tr AX = Tr AX X * = A where A and X are Hermitian matrices. Since W i , H j , and G l are Hermitian matrices for all i , j , l , we derive the gradient of Lagrangian in (10). As a result, we obtain the KKT optimality conditions used for the proof as follows:
j = 1 M λ j η j 1 θ j H j + α I β i a i H i + j i M β j H j + l = 1 L γ l G l E i = 0 , i
E i W i = 0 , i
α , β i , γ l 0 ; E i , W i 0 , i , l ,
where (10) and (11), respectively, are calculated from W i L = 0 and Tr E i W i = 0 with E i , W i 0 . In addition, I C N × N is an identity matrix. From (10), we derive
E i = A β i + β i a i H i , i ,
where
A = α I + j = 1 M β j H j + l = 1 L γ l G l j = 1 M λ j η j 1 θ j H j .
First, α > 0 is proved by contradiction. Assuming that α = 0 , let us introduce
D i = [ h 1 , . . . , h i 1 , h i + 1 , . . . , h M , g 1 , . . . , g L ] , i .
Since D i C N × ( M + L 1 ) , we set a basic of Null D i H as L i C N × ( N M L + 1 ) . In this paper, the number of antennas N is assumed to be larger than M + L 1 . Thus, we derive H j L i = 0 , j i , and G l L i = 0 , i , l . After that, with x C ( N M L + 1 ) × 1 and y = L i x , we obtain
y H E i y = λ i η i 1 θ i + β i a i y H H i y
y H E i y = λ i η i 1 θ i + β i a i h i H L i x 2 0 .
Owing to E i 0 and (16), we obtain the following:
y H E i y 0 , y h i H L i x = 0 , x h i H L i = 0 .
From the fact that the transmission channels h i , g l , i , l , are independent and random, we can assume that h i Range D i . Thus, we derive h i H L i 0 , which contradicts (17). Subsequently, α > 0 is satisfied. Moreover, since G l 0 , we have
α I + l = 1 L γ l G l 0 .
Second, we also prove that A 0 by contradiction. We assume that A 0 , i.e., there exists at least a vector z 0 such that z H Az 0 . From (13), we derive
z H E i z = z H Az β i + β i a i z H H i z , i .
With E i , H i 0 , we have z H E i z 0 , and z H H i z 0 . Thus, from (19) we derive z H A z = 0 , and z H H i z = 0 , i . From (14), we derive z H α I + l = 1 L γ l G l z = 0 , which contradicts (18). It thus follows that A 0 .
Last, we prove that rank W i = 1 , i . We know that rank ( X Y ) rank ( X ) rank ( Y ) is a basic rank inequality property. Thus, according to (13), we derive the following:
rank E i rank A rank β i + β i a i H i = N 1
since rank A = N and rank H i = 1 . From (11), we have Range W i Null E i , so rank W i N rank E i . Then, from (20) we derive rank W i 1 . Since W i = 0 does not satisfy constraint (8c), we can finally conclude that rank W i = 1 . Thus, Lemma 1 is completely proven.
Remark 1.
Inspired by the idea in [3], we can propose another approach to prove Lemma 1 by considering the minimization of the partial Lagrangian and the dual problem. First, we have the partial Lagrangian of the primal problem (8)-SDR expressed as follows:
L W i , α , β i , γ l = i = 1 M λ i η i 1 θ i j = 1 M Tr H i W j + σ n 2 + α i = 1 M Tr W i P max + i = 1 M β i 1 1 a i a i Tr H i W i + j i M Tr H i W j + σ n 2 + σ v 2 θ i + l = 1 L γ l i = 1 M Tr G l W i I t
where α 0 , β i 0 , i , and γ l 0 , l are the dual variables associated with constraints (8b), (8c), and (8d), respectively. The dual function of dual problem is expressed as
g α , β i , γ l = min W i 0 , i L W i , α , β i , γ l .
The dual problem is expressed as follows:
max α , β i , γ l g α , β i , γ l
s . t . α 0 , β i 0 , γ l 0 , i , l .
We rewrite the partial Lagrangian function and equivalently obtain the minimization problem as:
min W i 0 , i i = 1 M Tr j = 1 M λ j η j 1 θ j H j + α I β i a i H i + j i M β j H j + l = 1 L γ l G l W i .
That is equivalent to the optimization problem for separate variables W i :
min W i 0 Tr F i W i , i
where we denote F i = j = 1 M λ j η j 1 θ j H j + α I β i a i H i + j i M β j H j + l = 1 L γ l G l , i . For simple notation, we denote α , β i , γ l as the optimal dual solution of (23). Then, we derive F i 0 by contradiction. Suppose that F i is not a positive semidefinite matrix, i.e., there exists at least a vector x 0 such that x H F i x < 0 . Therefore, if we choose W i = t x x H where t > 0 , then Tr F i W i = Tr F i t x x H = t x H F i x < 0 is unbounded below when t + . Therefore, we can not obtain the bounded optimal dual value. Hence, F i 0 . From (25), we derive Tr F i W i = 0 and then F i W i = 0 since F i 0 and W i 0 , i . The rest of the rank-1 proof is similar to the KKT method from Equation (13) where F i has the same role as E i in (13).
Remark 2.
The real-valued objective and constraint functions with complex-valued variables in the problem (8)-SDR exist with partial derivatives according to Definition 2.2 ([31], Chapter 2) or Theorem 1 ([32], Section 2). These functions are affine and convex in complex variables W i according to Remark 3.15 ([33], Chapter 3). Thus, the problem (8)-SDR is a convex optimization problem with respect to complex-valued matrix variables. Moreover, KKT conditions can be applied to the problem (8)-SDR with real-valued functions and complex-valued variables from Remark 9.13 ([33], Chapter 9). Another approach to solving problem (8)-SDR is to convert problems with complex-valued variables to equivalent problem with real-valued variables according to Remark 1.20 ([33], Chapter 1). However, the proof of rank-1 constraints which we have not obtained yet, is still much more complicated than that of the problem (8)-SDR with complex-valued variables, although the same optimal values are achieved and the rank-1 constraints are satisfied in numerical experiments.

3.2. PSO-SDR Approach to Maximizing WSHE

We can apply the brute-force method to search the optimal PS ratios over all possible collections of θ i . Since the brute-force search method has very high computational complexity, i.e., in the M-dimension search, it is very slow to obtain the optimal results. To avoid this difficulty, a PSO-based algorithm [22,34,35] with low complexity, high convergence rate, and high rigor is exploited to search the best power-splitting ratios in the paper.
A main description of the PSO-based algorithm is presented as follows. The limited number of iterations and the number of elements in a swarm are denoted by T max and N S , respectively. We assign each element’s position to a set of M PS ratios θ 1 , . . . , θ M . By observing the SINR constraint (8c) in problem (8), we do not obtain this constraint if θ i is extremely small. As a result, we set a lowest value for θ i , represented by θ min . Therefore, θ min , 1 is the search interval of each PS ratio θ i . For the n-th particle, x n , v n , and p b , n , are represented for its position, velocity and local optimal position, respectively. The global best position denoted as g b is obtained by collecting information from all the particles. With the set of PS ratios, x n , we solve problem (8) with the SDR technique to achieve the maximum WSHE value denoted as f x n . For the update step, the new velocity is affected by the previous velocity with the inertia weight, i w , and the local and global best positions with the cognitive and social factors, denoted as c 1 and c 2 , respectively. Finally, problem (6) is solved by the proposed PSO-based algorithm described in detail in Table 2.
We can calculate the computational complexity of the PSO-based algorithm according to T max , N S , and the complexity in solving the SDP problem (8)-SDR. In the PSO method, the problem (8)-SDR is solved within T max N S times. In the problem (8)-SDR, there exist M + L + 1 linear constraints for M matrix variables, and each matrix variable has a size N × N . Therefore, the complexity of solving problem (8)-SDR is O M N M 3 N 6 + M + L + 1 M N 2 log 1 1 ζ ζ [36]. As a result, the computational complexity of the proposed algorithm is O T max N S M N M 3 N 6 + M + L + 1 M N 2 log 1 1 ζ ζ , for a given solution accuracy ζ .
The proposed solutions are not only suitable for the conventional cognitive radio (CR) network but also CR sensor networks where the secondary receivers are the wireless sensor nodes with a limited energy and low complexity due to the following reason. After the channel side information (CSI) from the secondary transmitter to the secondary receivers and primary users is obtained, the proposed SDR-PSO algorithm in Table 2 is performed at the secondary transmitter, which is usually a base station equipped with powerful computational processors. After that, the optimal power-splitting ratios are sent to the secondary receivers for setting the PS ratios by the secondary transmitter. As a result, the secondary receivers can split the incoming signal by optimal PS ratios without running the algorithms. Therefore, the proposed algorithms can be applied to low resources-constrained CR sensor networks.

3.3. Zero-Forcing Beamforming (ZFBF) Problem

In the ZFBF problem, we choose the weighted vector for the information message to the S R i which are orthogonal to the channels from the ST to the other S R j , j i , and the PUs, i.e., w i N u l l B i H where B i = [ h 1 , . . . , h i 1 , h i + 1 , . . . , h M , g 1 , . . . , g L ] C N × ( M + L 1 ) . The null-space of B i H has an orthonormal basic Q i C N × N ( M + L 1 ) where Q i H Q i = I . As a result, w i can be expressed w i = Q i u i where u i C ( N ( M + L 1 ) ) × 1 . Let us introduce U i = u i u i H . Some calculations are then performed as follows:
Tr W i = Tr Q i u i u i H Q i H = Tr Q i H Q i u i u i H = Tr U i , i
Tr H i W i = Tr H i Q i u i u i H Q i H = Tr Q i H H i Q i U i , i
Tr H j W i = 0 , Tr G l W i = 0 , j i , l .
As a result, the optimization problem with ZFBF beamforming is formulated as follows:
max U i , θ i i = 1 M λ i η i 1 θ i Tr Q i H H i Q i U i + σ n 2
s . t . i = 1 M Tr U i P max 0
Tr Q i H H i Q i U i a i + σ n 2 + σ v 2 θ i 0 , i
U i 0 , 0 < θ i < 1 , i
rank U i = 1 .
The ZFBF problem is solved by using the similar solution of the general problem (6). The optimal value of this ZFBF problem is compared with the optimal WSHE value of (6) in the numerical section.

3.4. Equal Power Splitting (EPS) Problem

In the special case of EPS design, the power-splitting factors θ i are simply assigned to be 0 . 5 , i . It means that the secondary receivers share a half received signal energy for decoding information and a half for harvesting power. The secondary transmitter can control the direction and power of beamforming vectors to achieve the WSHE maximization without join to the secondary receivers. Thus, the EPS problem is formulated as follows:
min W i 0 . 5 i = 1 M λ i η i j = 1 M Tr H i W j + σ n 2
s . t . i = 1 M Tr U i P max 0
Tr H i W i a i + j i M Tr H i W j + σ n 2 + σ v 2 0 . 5 0 , i
i = 1 M Tr G l W i I t 0 , l
W i 0 , i
rank W i = 1 , i .
For the EPS problem with fixed PS factors, we only have the variables of beamforming vectors in the optimization problem. Therefore, we can obtain the optimal solution by SDR technique and CVX solver. The optimal WSHE value in the EPS problem is considered as a special case in comparision with the general WSHE problem (6).

4. Robust WSHE Maximization with Imperfect Channels

The perfect CSI assumption is not practical due to many reasons such as the mobility of users, the estimation errors, limited CSI feedback quantization, etc. Therefore, we propose a robust design for the imperfect CSI of the channels. First, the models of imperfect CSI from the ST to the SRs and the PUs, are respectively considered as follows:
h i = h ^ i + Δ h i , Δ h i C N × 1 , Δ h i ε i
g l = g ^ l + Δ g l , Δ g l C N × 1 , Δ g l δ l ,
where h i and g l are the actual channels; h ^ i and g ^ l are the estimated channels at the ST; and Δ h i and Δ g l are the CSI errors. The error bounds, ε i and δ l , represent the radius of the uncertainty region of the estimated CSI channels. When these radiuses go to zero, the estimated channels become the perfect channels. For the imperfect CSI case, the optimization problem (6) becomes non-convex and has infinite number of constraints. The idea to solve this problem is to apply the vector inequalities and the S-Procedure [26,33] in order to obtain a tractable SDP problem. The robust WSHE maximization problem can be expressed as problem (29).
max w i , θ i min Δ h i ε i , i i = 1 M λ i η i 1 θ i j = 1 M h ^ i H + Δ h i H w j 2 + σ n 2
s . t . i = 1 M Tr W i P max
θ i h ^ i H + Δ h i H w i 2 θ i j = 1 , j i M h ^ i H + Δ h i H w j 2 + σ n 2 + σ v 2 a i , Δ h i ε i , i
i = 1 M g ^ l H + Δ g l H w i 2 I t , Δ g l δ l , l
0 < θ i < 1 , i .
First, we consider the objective function (29a) as follows
min Δ h i ε i , i i = 1 M λ i η i 1 θ i j = 1 M h ^ i H + Δ h i H w j 2 + σ n 2
From the inequalities x + y x y and u H v u v , we obtain the results as follows
h ^ i H + Δ h i H w j h ^ i H w j Δ h i H w j h ^ i H w j ε i w j .
Subsequently, we have
h ^ i H + Δ h i H w j 2 h ^ i H w j ε i w j 2
h ^ i H w j ε i w j 2 = h ^ i H w j 2 + ε i 2 w j 2 2 ε i h ^ i H w j w j
h ^ i H w j ε i w j 2 w j H h ^ i h ^ i H w j + ε i 2 w j H w j 2 ε i h ^ i w j 2
h ^ i H + Δ h i H w j 2 w j H H ˜ i w j ,
where H ˜ i = h ^ i h ^ i H + ε i 2 2 ε i h ^ i I . Thus, the objective function is recast as follows
max w i , θ i i = 1 M λ i η i 1 θ i j = 1 M w j H H ˜ i w j + σ n 2 .
Next, we apply S-Procedure to convert the SINR and PU constraints to linear matrix inequality (LMI) constraints.
Lemma 2.
(S-Procedure) We have the function
f k x = x H A k x + 2 Re b k H x + c k ,
where k = 1 , 2 , A k H N , b k C N × 1 , and c k C . Then, the implication f 1 x 0 f 2 x 0 holds if and only if there exists a ξ 0 such that
ξ A 1 b 1 b 1 H c 1 A 2 b 2 b 2 H c 2 0 ,
provided that there exists a point x ¯ such that f 1 x ¯ < 0 .
We first analyze the SINR constraints (29c) at the SRs as follows:
f 1 Δ h i = Δ h i H I Δ h i + 2 Re 0 H Δ h i + ε i 2 0
f 2 Δ h i = h ^ i H + Δ h i H w i w i H h ^ i + Δ h i a i + h ^ i H + Δ h i H j i M w j w j H h ^ i + Δ h i + σ n 2 + σ v 2 θ i 0 .
We can rewrite f 2 Δ h i as follows
f 2 Δ h i = Δ h i H A i Δ h i + 2 Re A i h ^ i H Δ h i + h ^ i H A i h ^ i + σ n 2 + σ v 2 θ i 0 ,
where A i = j i M w j w j H w i w i H a i . Note that A i , i , are Hermitian matrices.
According to Lemma 2, there exists α i 0 satisfying
α i I 0 0 ε i 2 A i A i h ^ i h ^ i H A i h ^ i H A i h ^ i + σ n 2 + σ v 2 θ i 0
α i I A i A i h ^ i h ^ i H A i α i ε i 2 h ^ i H A i h ^ i + σ n 2 + σ v 2 θ i 0 , i .
Finally, we consider the interference constraints (29d) at the PUs as follows:
f 1 Δ g l = Δ g l H I Δ g l + 2 Re 0 H Δ g l + δ l 2 0
f 2 Δ g l = Δ g l H B Δ g l + 2 Re B g ^ l H Δ g l + g ^ l H B g ^ l I t 0 ,
where B = i = 1 M w i w i H . Note that B is Hermitian matrix.
According to Lemma 2, there exists β l 0 satisfying
β l I 0 0 δ l 2 B B g ^ l g ^ l H B g ^ l H B g ^ l I t 0 , l
β l I B B g ^ l g ^ l H B β l δ l 2 g ^ l H B g ^ l + I t 0 , l .
From (30), (32), and (33), we obtain the robust WSHE problem with imperfect CSI channels as follows
max w i , θ i , α i , β l i = 1 M λ i η i 1 θ i j = 1 M w j H H ˜ i w j + σ n 2
s . t . ( 29 b ) , ( 32 ) , ( 33 ) , ( 29 e )
α i 0 , β l 0 , i , l .
We denote W i = w i w i H , and transform the robust problem to the equivalent problem as follows
max W i , θ i , α i , β l i = 1 M λ i η i 1 θ i j = 1 M Tr H ˜ i W j + σ n 2
s . t . i = 1 M Tr W i P max 0
( 32 ) , A i = j i M W j W i a i , i
( 33 ) , B = i = 1 M W i , l
0 < θ i < 1 , α i 0 , β l 0 , i , l
rank W i = 1 , i .
This optimization problem is solved by an approach similar to the WHSE problem (6) with perfect channels. For fixed power splitting ratios θ i , by relaxing the rank-1 constraints, the robust problem becomes the SDP form and can be solved by numerical solvers. By using the PSO method, we can obtain the optimal value of the robust problem by relaxing the rank-1 constraints. However, since the solutions do not always guarantee the properties of rank-1, we only obtain the upper bound of the robust WSHE value. The rank-1 suboptimal solutions have not found yet in the paper, and are left for future work where the randomization method [20] can be applied.

5. Simulation Results

In this section, the Monte Carlo simulations are performed to evaluate the numerical results of the proposed algorithms. The simulation parameters are given by: the number of SRs and PUs, M = 2 , L = 2 , respectively; the maximum transmit power of ST, P max = 30 dBW ; the variance of noises, σ n 2 = 75 dBm , σ v 2 = 55 dBm ; and the energy harvesting efficiency of all SRs, η i = 1 , i . The weighted factors for two SRs are given as λ 1 , λ 2 = 0 . 2 , 0 . 8 . We assume that the power attenuation of the channels from the ST to the SRs and the PUs are identical and equal, at 50 dB . Each entry of channel vectors is randomly generated from i.i.d. Rayleigh flat-fading according to the above power attenuation.
We define the normalized maximum channel estimation error of secondary receivers and primary users as ε i nor = ε i h i , δ l nor = δ l g l , i , l . We assume that the normalized maximum channel estimation errors ε i nor = δ l nor = 0 . 001 , i , l . For the PSO method in Table 2, according to the experiments in [34,35], the simulation parameters are set as follows: i w = 0 . 7 , c 1 = 1 . 494 , and c 2 = 1 . 494 , with which the good convergence can be obtained in the proposed algorithm. Furthermore, the swarm size, the maximum number of iterations, and the minimum value of θ i , are set as N S = 10 , T max = 30 , and θ min = 0 . 001 , respectively.
Figure 2 and Figure 3 show the convergence characteristics of the WSHE of the proposed scheme according to the number of iterations in both cases of perfect and imperfect CSI, respectively, when the number of antennas, N, is 10. We can see that the fitness function (the WSHE) can quickly converge in less than 30 iterations. As a result, we use the maximum iteration number T max of 30 for futher simulations. Furthermore, it is observed that the increase of the minimum required SINR corresponds to the decreases of the WSHE of which reasons will be explained in Figure 4 in more detail.
Figure 4 shows average weighted sum harvested power according to the minimum required SINR when the interference threshold I t = 45 dBm , and the number of antennas N = 10 . Two suboptimal zero-forcing beamforming (ZFBF) and equal power splitting (EPS) ratios are also considered as the baseline schemes in comparison with the optimal design in perfect CSI case. The design parameters of the ZFBF is obtained by the solution of problem (27) in Section 3.3. In case of the EPS design, the SRs equally divides the received signal for the ID part and the EH part, i.e., θ i = 0 . 5 , i . Therefore, the optimization problem (28) has only the beamforming vectors variables in the EPS design. Also, the results of problem (6) without the PU constraints are also provided. In addition, the brute-force search (BFS) method with M-dimension search for θ i , i = 1 , . . . , M , is also considered to check the optimal value obtained by the proposed PSO-based method.
As shown in Figure 4, the maximum weighted sum harvested energy decreases as the required minimum SINR increases in all cases except the ZFBF scheme. The reason is that the secondary receivers should increase the power splitting ratios and adjust the beamforming vectors in order to get more energy in information decoding parts and further to satisfy the increasing SINR constraints. Hence, the secondary receivers harvest less energy in energy harvesting parts. The ZFBF scheme is not sensitive to the SINR value because the SINR constraints in (27c) are easily satisfied with small noise variances. Also, Figure 4 shows that the optimal solution (6) with perfect CSI achieves better WSHE than ZFBF, EPS, and robust imperfect CSI (35) since both power splitting ratios and beamforming vectors are optimized. The case without primary users gives better WSHE since no PU interference constraints are considered. Furthermore, Figure 4 shows that the proposed PSO-based method achieves the same optimal value as BFS method in both perfect and robust imperfect CSI cases.
Figure 5 shows the average weighted sum harvested energy according to the PU interference threshold when the minimum required SINR a i = 5 dB , and the number of antennas N = 10 . We observe that the WSHE increases with the increasing interference threshold in all schemes except the no PUs and ZFBF schemes. The reason is that when the PU interference threshold increases, the beamforming vectors of the secondary transmitter are more flexible to transmit information and energy simultaneously. The proposed scheme with perfect CSI gives better performance than ZFBF, EPS and the proposed robust schemes. Obviously, the scheme without the constraints of PU interference provides the highest performances. Moreover, when the interference threshold gets smaller, the gap between the proposed perfect CSI and no PUs schemes gets larger because the beamforming vectors have to limit the interference leakage to primary users. On the other hand, when the interference threshold is high, the WSHE of the proposed perfect CSI solution is similar to that of the no PUs scheme because the PU interference constraints are easily obtained. Also, Figure 5 shows the proposed PSO-based method achieves the same performance as the BFS-based method in both perfect and robust imperfect CSI schemes.
Finally, Figure 6 shows the weighted sum harvested energy according to the number of antennas N when the minimum required SINR a i = 5 dB , and the interference threshold I t = 45 dBm . It is observed that the weighted sum harvested energy increases as the number of antennas increases, in all schemes. When the secondary transmitter has more antennas, it can exploit the extra degrees of freedom in optimizing the direction of the beamforming vectors. Thus, the transmitted signal can be more exactly steered towards the secondary receivers. Consequently, we can obtain more energy harvesting at the secondary receivers and less interference to primary users. Similar to Figure 4 and Figure 5, the performances of the proposed scheme with perfect CSI are better than those of ZFBF, EPS and the proposed robust solution with imperfect CSI. Moreover, the proposed method gives the same harvested energy as that of BFS method in both perfect and imperfect CSI cases while requiring much lower computational complexity.

6. Conclusions

We have investigated a multiuser SWIPT cognitive radio network to maximize the weighted sum harvested energy of the secondary receivers under the constraints on total transmit power at the secondary transmitter, the required SINR at the secondary receivers with power-splitting structure, and the interference level at the primary users. The efficient algorithms based on SDR and PSO methods are proposed to find the optimal beamforming at the secondary transmitter and the optimal power-splitting ratios at the secondary receivers for both cases of the perfect and imperfect CSIs. Through simulation, it was shown that the proposed PSO-based algorithm achieves a fast convergence within 30 iterations. It was also shown that the obtained maximum WSHEs in both cases of perfect and imperfect CSI are increased when the allowable interference level at primary users as well as the number of antenna at secondary transmitter get larger.
In the paper, however, the upper bound of WSHE is only archived for the robust imperfect CSI case, and subsequently, the optimal solution satisfying the rank-1 constraints needs to be investigated more in future. In addition, the SDR-based solution has large computational complexity since the number of variables is squared when vector variables are changed to matrix variables. Therefore, the SDR-based approach may not be appropriate to the large-scale system where the antenna number N is very large and the complexity of SDR method is O N 6 . 5 . Thus, fast solutions should be investigated in future works such as with respect to the alternative direction method of multipliers [37].

Acknowledgments

This work was supported by the National Research Foundation of Korea funded by the MEST under Grant NRF 2015R1A2A1A15053452 and Grant NRF 2015R1D1A1A09057077.

Author Contributions

The first and second authors conceived the idea. The first author performed the simulation experiments and wrote the paper. The second author supervised the research, critically reviewed and revised the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bi, S.; Ho, C.K.; Zhang, R. Wireless powered communication: Opportunities and challenges. IEEE Commun. Mag. 2015, 53, 117–125. [Google Scholar] [CrossRef]
  2. Ding, Z.; Ng, D.W.K.; Peng, M.; Suraweera, H.A.; Schober, R.; Poor, H.V. Application of smart antenna technologies in simultaneous wireless information and power transfer. IEEE Commun. Mag. 2015, 53, 86–93. [Google Scholar] [CrossRef]
  3. Shi, Q.; Liu, L.; Xu, W.; Zhang, R. Joint transmit beamforming and receive power splitting for MISO SWIPT systems. IEEE Trans. Wirel. Commun. 2014, 13, 3269–3280. [Google Scholar] [CrossRef]
  4. Vu, Q.D.; Tran, L.N.; Farrel, R.; Hong, E.K. An Efficiency Maximization Design for SWIPT. IEEE Signal Process. Lett. 2015, 22, 2189–2193. [Google Scholar] [CrossRef]
  5. Zhang, R.; Ho, C.K. MIMO Broadcasting for Simultaneous Wireless Information and Power Transfer. IEEE Trans. Wirel. Commun. 2013, 12, 1989–2001. [Google Scholar] [CrossRef]
  6. Ding, Z.; Krikidis, I.; Sharif, B.; Poor, H.V. Wireless information and power transfer in cooperative networks with spatially random relays. IEEE Trans. Wirel. Commun. 2014, 13, 4440–4453. [Google Scholar] [CrossRef]
  7. Timotheou, S.; Krikidis, I.; Zheng, G.; Ottersten, B. Beamforming for MISO Interference Channels with QoS and RF Energy Transfer. IEEE Trans. Wirel. Commun. 2014, 13, 2646–2658. [Google Scholar]
  8. Feng, R.; Li, Q.; Zhang, Q.; Qin, J. Robust Secure Transmission in MISO Simultaneous Wireless Information and Power Transfer System. IEEE Trans. Veh. Technol. 2015, 64, 400–405. [Google Scholar] [CrossRef]
  9. Liao, J.; Khandaker, M.R.A.; Wong, K.K. Robust power-splitting SWIPT beamforming for broadcast channels. IEEE Commun. Lett. 2016, 20, 181–184. [Google Scholar] [CrossRef]
  10. 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]
  11. Chu, Z.; Zhu, Z.; Xiang, W.; Hussein, J. Robust beamforming and power splitting design in MISO SWIPT downlink system. IET Commun. 2016, 10, 691–698. [Google Scholar] [CrossRef]
  12. Nasir, A.A.; Tuan, H.D.; Ngo, D.T.; Durrani, S.; Kim, D.I. Path-Following Algorithms for Beamforming and Signal Splitting in RF Energy Harvesting Networks. IEEE Commun. Lett. 2016, 20, 1687–1690. [Google Scholar] [CrossRef]
  13. Goldsmith, A.; Jafar, S.A.; Maric, I.; Srinivasa, S. Breaking spectrum gridlock with cognitive radios: An information theoretic perspective. Proc. IEEE 2009, 97, 894–914. [Google Scholar] [CrossRef]
  14. Mohjazi, L.; Dianati, M.; Karagiannidis, G.K.; Muhaidat, S. RF-powered cognitive radio networks: Technical challenges and limitations. IEEE Commun. Mag. 2015, 53, 94–100. [Google Scholar] [CrossRef]
  15. Zheng, G.; Ho, Z.; Jorswieck, E.A.; Ottersten, B. Information and Energy Cooperation in Cognitive Radio Networks. IEEE Trans. Wirel. Commun. 2014, 62, 2290–2303. [Google Scholar] [CrossRef]
  16. Ng, D.W.K.; Lo, E.S.; Schober, R. Multi-Objective Resource Allocation for Secure Communication in Cognitive Radio Networks with Wireless Information and Power Transfer. IEEE Trans. Veh. Technol. 2016, 65, 3166–3184. [Google Scholar] [CrossRef]
  17. Yang, Z.; Ding, Z.; Fan, P.; Karagiannidis, G.K. Outage performance of cognitive relay networks with wireless information and power transfer. IEEE Trans. Veh. Technol. 2016, 65, 3828–3833. [Google Scholar] [CrossRef]
  18. Lee, S.; Zhang, R. Cognitive wireless powered network: Spectrum sharing models and throughput maximization. IEEE Trans. Cogn. Commun. Netw. 2015, 1, 335–346. [Google Scholar] [CrossRef]
  19. Tuan, P.V.; Koo, I. Optimal Multiuser MISO Beamforming for Power-Splitting SWIPT Cognitive Radio Networks. IEEE Access 2017, 5, 14141–14153. [Google Scholar] [CrossRef]
  20. Luo, Z.Q.; Ma, W.K.; So, A.M.C.; Ye, Y.; Zhang, S. Semidefinite relaxation of quadratic optimization problems. IEEE Signal Process. Mag. 2010, 27, 20–34. [Google Scholar] [CrossRef]
  21. Grant, M.; Boyd, S. CVX: Matlab Software for Disciplined Convex Programming Version 2.1. Available online: http://cvxr.com/cvx (accessed on 22 February 2017).
  22. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the IEEE International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  23. Ghasemi, A.; Sousa, E.S. Fundamental limits of spectrum sharing in fading environment. IEEE Trans. Wirel. Commun. 2007, 6, 649–658. [Google Scholar] [CrossRef]
  24. 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. 2016, 15, 2312–2327. [Google Scholar] [CrossRef]
  25. Wu, Q.; Chen, W.; Ng, D.W.K.; Chen, W.; Schober, R. User-Centric Energy Efficiency Maximization for Wireless Powered Communication Networks. IEEE Trans. Wirel. Commun. 2016, 15, 6898–6912. [Google Scholar] [CrossRef]
  26. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  27. Bengtsson, M.; Ottersten, B. Optimal and suboptimal transmit beamforming. In Handbook of Antennas in Wireless Communications; CRC Press: Boca Raton, FL, USA, 2001; Chapter 18. [Google Scholar]
  28. Huang, Y.; Palomar, D.P. Rank-constrained separable semidefinite programming with applications to optimal beamforming. IEEE Trans. Signal Process. 2010, 58, 664–678. [Google Scholar] [CrossRef]
  29. Hjorungnes, A.; Gesbert, D. Complex-Valued Matrix Differentiation: Techniques and Key Results. IEEE Trans. Signal Process. 2007, 55, 2740–2746. [Google Scholar] [CrossRef]
  30. Zheng, G.; Wong, K.K.; Ottersten, B. Robust cognitive beamforming with bounded channel uncertainties. IEEE Trans. Signal Process. 2009, 57, 4871–4881. [Google Scholar] [CrossRef]
  31. Hjorungnes, A. Complex-Valued Matrix Derivatives with Applications; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  32. Brandwood, D.H. A complex gradient operator and its application in adaptive array theory. Proc. Inst. Elect. Eng. 1983, 130, 11–16. [Google Scholar]
  33. Chi, C.Y.; Li, W.C.; Lin, C.H. Convex Optimization for Signal Processing and Communications: From Fundamentals to Applications; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  34. Robinson, J.; Samii, Y.R. Particle swarm optimization in electromagnetics. IEEE Trans. Antennas Propag. 2004, 52, 397–407. [Google Scholar] [CrossRef]
  35. Tuan, P.V.; Koo, I. Throughput maximisation by optimising detection thresholds in full-duplex cognitive radio networks. IET Commun. 2016, 10, 1355–1364. [Google Scholar] [CrossRef]
  36. Karipidis, E.; Sidiropoulos, N.D.; Luo, Z.Q. Quality of service and max-min-fair transmit beamforming to multiple co-channel multicast groups. IEEE Trans. Signal Process. 2008, 56, 1268–1279. [Google Scholar] [CrossRef]
  37. Huang, K.; Sidiropoulos, N.D. Consensus-ADMM for general quadratically constrained quadratic programming. IEEE Trans. Signal Process. 2016, 64, 5297–5310. [Google Scholar] [CrossRef]
Figure 1. (a) The multi-input single-output (MISO) cognitive radio network being considered in the paper, and (b) the power-splitting receive model at SR i . PU: primary user; SR: secondary receiver; PS: power splitting; ID: information decoding; EH: energy harvesting.
Figure 1. (a) The multi-input single-output (MISO) cognitive radio network being considered in the paper, and (b) the power-splitting receive model at SR i . PU: primary user; SR: secondary receiver; PS: power splitting; ID: information decoding; EH: energy harvesting.
Sensors 17 02275 g001
Figure 2. The iteration behavior of the proposed PSO-based algorithm when the required signal-to-interference-plus-noise ratios (SINRs) are given as 2, 4, 6, 8, and 10 dB, respectively and a random perfect CSI is considered.
Figure 2. The iteration behavior of the proposed PSO-based algorithm when the required signal-to-interference-plus-noise ratios (SINRs) are given as 2, 4, 6, 8, and 10 dB, respectively and a random perfect CSI is considered.
Sensors 17 02275 g002
Figure 3. The iteration behavior of the proposed PSO-based algorithm when the required signal-to-interference-plus-noise ratios (SINRs) are given as 2, 4, 6, 8, and 10 dB, respectively, and a random imperfect CSI is considered.
Figure 3. The iteration behavior of the proposed PSO-based algorithm when the required signal-to-interference-plus-noise ratios (SINRs) are given as 2, 4, 6, 8, and 10 dB, respectively, and a random imperfect CSI is considered.
Sensors 17 02275 g003
Figure 4. Average weighted sum harvested energy according to minimum required SINR. BFS: brute-force search.
Figure 4. Average weighted sum harvested energy according to minimum required SINR. BFS: brute-force search.
Sensors 17 02275 g004
Figure 5. Average weighted sum harvested energy according to the interference threshold.
Figure 5. Average weighted sum harvested energy according to the interference threshold.
Sensors 17 02275 g005
Figure 6. Average weighted sum harvested energy according to number of antennas.
Figure 6. Average weighted sum harvested energy according to number of antennas.
Sensors 17 02275 g006
Table 1. Abbreviations and meanings.
Table 1. Abbreviations and meanings.
Abbre.Full FormMeanings
CRNsCognitive radio networks [13]The CRNs use licensed bands of primary networks for communication under overlay or underlay access modes. In this paper, the underlay access mode is applied where the interference from secondary transmitter to primary users is lower than the prescribed threshold.
SWIPTSimultaneous wireless information and power transfer[1]In this paper, the transmitter simultaneously sends both information and power to the receivers which are equipped with an information decoder and an energy harvester.
WSHEWeighted sum harvested energyWSHE is the sum of all harvested energy of secondary receivers where each value has one weight factor.
CSIChannel state informationCSI is the complex value of the baseband channel between the transmitter and the receiver.
SDRSemidefinite relaxation [20]The fundamental idea of SDR is based on the properties as follows: X = x x H X 0 and rank X = 1 . Thus, we can change the variable from x to X by adding the constraints of X 0 and rank X = 1 . Then, SDR technique removesthe constraint of rank X = 1 to obtain the semidefinite programming problem solved by the numerical solver CVX [21]. This optimal solution needs to to be checked for satisfying the rank-1 constraint, or not.
PSOParticle swarm optimization [22]PSO is a heuristic evolutional search algorithm which is based on simulating a swarm of particles (birds, fish …) which share information of positions and fitness values with each other. For each iteration, PSO basically performs three steps, called velocity update, position update, and evaluation of the fitness function. Both deterministic and probabilistic rules are used to search overall design space. Each particle moves to a new position with likely improvement of fitness value.
ZFBFZero-forcing beamformingIn the ZFBF case, the transmitter designs the beamforming vector, bringing the intended message for one user and being orthogonal to the CSI of other users. Thus, ZFBF design does not cause interference to other users.
EPSEqual power splittingThe receivers will share half the power of theincoming signal for information decoder, and half for the energy harvester.
Table 2. The PSO-based algorithm for solving the WSHE problem (6).
Table 2. The PSO-based algorithm for solving the WSHE problem (6).
1:input parameters: T max , N S , θ min , v max , i w , c 1 , c 2 ,
and variables x n , n .
2:initialization:
2:Assign the iteration index of PSO loop: m = 1 .
3:The x n ’s values of elements, n , are
assigned randomly in θ min , 1 , then obtain f x n by solving (8).
4:Assign the global maximum value: f g b = max 1 n N S f x n .
5:Set the best position of particle: p b , n = x n , n .
Set the velocity of particle: v n = 0 , n .
6:repeat(PSO loop)
7:for n = 1 : N S do
8:  Calculate particle’s new velocity:
   v n i w v n + c 1 π 1 , n p b , n x n + c 2 π 2 , n g b x n ,
  ⊙ denotes the Hadamard product, and the vectors π 1 , n , π 2 , n
  have independent uniformly distributed elements in 0 , 1 .
9:  Restrict vector v n ’s each element in v max , v max .
10:  Calculate particle’s new position: x n x n + v n .
11:  Restrict each element of vector x n in θ min , 1 .
12:  Accessment: Compute f x n and the optimal beamforming
  variables w i n , i from the solution of problem (8)-SDR
  according to the PS ratios set, x n .
13:  Update the new best position of particle:
  if f x n > f p b , n then
   Assign: p b , n x n .
  end if
14:  Update particle’s new global best position:
  if f x n > f g b then
   Assign: g b x n , w i * w i n , i .
  end if
15:end for
16: Update iteration index: m m + 1 .
17:until m > T max (end PSO loop)
18:final results: the global best value f g b is the optimal value of
WSHE problem (6) according to the optimal PS ratios θ i * = g b ,
and the optimal beamforming vectors w i * , i .

Share and Cite

MDPI and ACS Style

Tuan, P.V.; Koo, I. Robust Weighted Sum Harvested Energy Maximization for SWIPT Cognitive Radio Networks Based on Particle Swarm Optimization. Sensors 2017, 17, 2275. https://doi.org/10.3390/s17102275

AMA Style

Tuan PV, Koo I. Robust Weighted Sum Harvested Energy Maximization for SWIPT Cognitive Radio Networks Based on Particle Swarm Optimization. Sensors. 2017; 17(10):2275. https://doi.org/10.3390/s17102275

Chicago/Turabian Style

Tuan, Pham Viet, and Insoo Koo. 2017. "Robust Weighted Sum Harvested Energy Maximization for SWIPT Cognitive Radio Networks Based on Particle Swarm Optimization" Sensors 17, no. 10: 2275. https://doi.org/10.3390/s17102275

APA Style

Tuan, P. V., & Koo, I. (2017). Robust Weighted Sum Harvested Energy Maximization for SWIPT Cognitive Radio Networks Based on Particle Swarm Optimization. Sensors, 17(10), 2275. https://doi.org/10.3390/s17102275

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