Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: wireless powered communication network has arisen as a new system with stable and self - sustainable power supplies in shaping future - generation wireless communications @xcite . the enabling technology , known as simultaneous wireless information and power transfer ( swipt ) , has particularly drawn an upsurge of interests owing to the far - field electromagnetic power carried by radio - frequency ( rf ) signals that affluently exist in wireless communications . with the transmit power , waveforms , and dimensions of resources , etc . , being all fully controllable , swipt promises to prolong the lifetime of wireless devices while delivering the essential communication functionality , as will be important for low - power applications such as rf identification ( rfid ) and wireless sensor networks ( wsns ) ( see @xcite and the references therein ) .    on the other hand , privacy and authentication have increasingly become major concerns for wireless communications and physical ( phy)-layer security has emerged as a new layer of defence to realize perfect secrecy transmission in addition to the costly upper - layer techniques such as cryptography . in this regard , relay - assisted secure transmission was proposed @xcite and phy - layer security enhancements by means of cooperative communications have since attracted much attention @xcite . in particular , cooperative schemes can be mainly classified into three categories : _ decode - and - forward ( df ) _ , _ amplify - and - forward ( af ) _ , and _ cooperative jamming ( cj ) _ @xcite with cj being most relevant to phy - layer security . specifically , coordinated cj refers to the scheme of generating a common jamming signal across all single - antenna relay helpers against eavesdropping @xcite , while uncoordinated cj considers that each relay helper emits independent artificial noise ( an ) to confound the eavesdroppers @xcite . it is expected that in the scenarios where the direct link is broken between the transmitter ( tx ) and the legitimate receiver ( rx ) , some of the relays have to take on their conventional role of forwarding the information while others will perform cj @xcite . a recent paradigm that generalizes all the above - mentioned cooperation strategies is _ cooperative beamforming ( cb ) mixed with cj _ @xcite , where the available power at each relay is split into two parts : one for forwarding the confidential message and the other for cj . however , mixed cb - cj approaches may be prohibitive in applications with low power devices because idle relays with limited battery supplies would likely prefer saving power for their own traffic to assisting others communication . in light of this , swipt provides the incentive for potential helpers to perform dedicated cb mixed with cj at no expense of its own power , but opportunistically earn harvested energy . motivated by this , our work considers secrecy transmission from a tx to a legitimate rx with the aid of a set of single - antenna wireless energy harvesting ( weh)-enabled af - operated relays in the presence of multiple single - antenna eavesdroppers . as a matter of fact , cooperative schemes that involve weh - enabled relays was recently investigated in @xcite . we consider the use of the _ dynamic power splitting ( dps ) _ receiver architecture , initially proposed for swipt in @xcite , which divides the received power with an adjustable ratio for energy harvesting ( eh ) and information receiving ( ir ) . weh - enabled relays using dps receivers have also been considered in @xcite and @xcite , without ( w / o ) and with secrecy consideration , respectively . note that there is also interest in addressing the threat that weh receivers may attempt to intercept the confidential messages in swipt - enabled networks @xcite . nevertheless , we will focus on exploiting the benefits of weh - enabled relays when they are trustful . in particular , motivated by the strong interest in swipt and the vast degree - of - freedom ( dof ) achievable by cooperative relays , this paper aims to maximize the secrecy rate with the aid of weh - enabled af - operated relays , subject to the eh power constraints of individual relays by jointly optimizing the cb of the relays and the cj covariance matrix . in this paper , we assume that there is no direct link between the source and destination nodes , and perfect global channel state information ( csi ) is available for the case of centralized optimization . it is worth pointing out that although our setting may look similar to @xcite , their optimal cb - cj design is not applicable to ours due to the multiplicative nature in beamforming weights incurred by the power splitting ( ps ) ratios that intrinsically poses more intractability to our optimization problem . further , our work also differs from @xcite where an efficient algorithm was proposed to maximize the secrecy rate for the optimization of the ps ratios and af relay beamforming . the difference is twofold . first , an was not considered in the second transmission phase in @xcite and in addition , their algorithm only converged to a local optimum , as opposed to our work that gives the global optimal solutions for cb . the rest of the paper is organized as follows . section  [ sec : system model ] describes two types of weh - enabled rx architecture for the af relays and defines the secrecy rate region of the relay wiretap channel . section  [ sec : problem formulation ] then formulates the secrecy rate maximization problems that jointly optimize the an ( or cj ) and the af - relay cb for the weh - enabled relays operating with the two types of rx . the problems are respectively solved by centralized schemes in section  [ sec : centralized ] and distributed approaches in section  [ sec : distributed ] . section  [ sec : numerical results ] provides simulation results to evaluate the performance of the proposed schemes . finally , section  [ sec : conclusion ] concludes the paper .    _ notations_we use the uppercase boldface letters for matrices and lowercase boldface letters for vectors . the superscripts @xmath0 , @xmath1 , @xmath2 and @xmath3 represent , respectively , the transpose , conjugate , conjugate transpose operations on vectors or matrices , and the optimum . in addition , @xmath4 stands for the trace of a square matrix . moreover , @xmath5_{i , j}$ ] denotes the @xmath6th entry of a matrix , while @xmath7 and @xmath8 represent the euclidean norm and the entry - wise absolute value square of a vector , respectively . also , @xmath9 denotes a diagonal matrix with its diagonal specified by the given vector and @xmath5_{i=1}^n$ ] represents an @xmath10 vector with each element indexed by @xmath11 . furthermore , @xmath12 and @xmath13 stand for product and hadamard product , respectively . @xmath14 denotes the field of complex ( real ) matrices with dimension @xmath15 and @xmath16}$ ] indicates the expectation operation . finally , @xmath17 is short for @xmath18 . in this paper , we consider secrecy transmission in a swipt - enabled wsn , where a tx ( alice ) wants to establish confidential communication with the legitimate rx ( bob ) with no direct link but with the aid of @xmath19 weh - enabled sensors operating as af relays , denoted by @xmath20 , in the presence of multiple eavesdroppers ( eves ) , denoted by @xmath21 , all equipped with single antenna . we assume that there is no direct link from the tx to any of the eves , due to , for instance , severe path loss or shadowing , as illustrated in fig . [ fig : system model ] . we consider a two - hop relaying protocol based on two equal time slots and the duration of one transmit - slot is normalized to be one unit so that the terms `` energy '' and `` power '' are interchangeable with respect to ( w.r.t . ) one transmit - slot .    at the receiver of each af relay , we introduce two types of weh - enabled receiver architecture , namely , _ static power splitting ( sps ) _ ( fig .  [ fig : subfig : sps rx ] ) and _ dps _ ( fig . [ fig : subfig : dps rx ] ) , both of which allow the relay to harvest energy and receive information from the same received signal . specifically , the receiver first splits a portion of @xmath22 , of the received power for eh and the rest @xmath23 for ir , @xmath24 . the @xmath22 portion of harvested power is further divided into two streams with @xmath25 used for generating the an to confound eves and @xmath26 used for amplifying the received signal , where @xmath27 is the @xmath11th element of the received signal @xmath28 , and @xmath29 denotes the eh efficiency . note that dps with adjustable @xmath22 s is presently the most general receiver operation because practical circuits can not directly decode the information from the stream used for eh @xcite and sps is just a special case of dps with @xmath30 , @xmath24 , fixed for the whole transmission duration . however , sps , advocated for its ease of implementation , is introduced separately in the sequel for its simplified relay beamforming design . in the first transmit - slot , the received signal at each individual relay can be expressed as @xmath31 where the transmit signal @xmath32 is a circularly symmetric complex gaussian ( cscg ) random variable with zero mean and unit variance , denoted by @xmath33 , @xmath34 denotes the complex channel from the tx to the @xmath11th relay , @xmath35 is the transmit power at the tx , and @xmath36 is the additive white gaussian noise ( awgn ) introduced by the receiving antenna of the @xmath11th relay , denoted by @xmath37 . as such , the linearly amplified baseband equivalent signal at the output of the @xmath11th relay is given by @xmath38 where @xmath39 denotes the complex af coefficient , and @xmath40 denotes the noise due to signal conversion from the rf band to baseband , denoted by @xmath41 . since @xmath42 is constrained by the portion of the harvested power for forwarding , i.e. , @xmath26 , @xmath39 is accordingly given by @xmath43 where @xmath44 denotes the phase of the af coefficient for the @xmath11th relay . next , we introduce the cj scheme . denote the cj signal generated from @xmath19 relays by @xmath45^t$ ] and define its covariance matrix as @xmath46 $ ] . then the coordinated cj transmission can be uniquely determined by the truncated eigenvalue decomposition ( evd ) of @xmath47 given by @xmath48 , where @xmath49)$ ] is a diagonal matrix with @xmath50 s denoting all the positive eigenvalues of @xmath47 and @xmath51 is the precoding matrix satisfying @xmath52 . note that @xmath53 denotes the rank of @xmath54 which will be designed later . as a result , the cj signal can be expressed as @xmath55 where @xmath56 s are drawn from the columns of @xmath57 , and @xmath58 s are independent and identically distributed ( i.i.d . ) complex gaussian variables denoted by @xmath59 , which is known as the optimal distribution for an @xcite . on the other hand , @xmath60 , @xmath24 , denotes the power constraint for jamming at the @xmath11th relay , which implies that @xmath61 where @xmath62 is a diagonal matrix with its diagonal @xmath63 ( a unit vector with the @xmath11th entry equal to @xmath64 and the rest equal to @xmath65 ) . note that the cj scheme proposed above is of the most general form . for the special case when @xmath66 , i.e. , @xmath67 , each relay transmits a common jamming signal @xmath68 with their respective weight drawn from @xmath69 @xcite . this case is desirable in practice since it has the lowest complexity for implementation . in summary , the transmitted signal at the @xmath11th relay is given by @xmath70    according to together with , , and , the transmit signal from all relays can be expressed in vector form as @xmath71 where @xmath72 and @xmath73 are , respectively , diagonal matrices with their diagonals composed of @xmath74 @xmath75 and @xmath76 . in addition , @xmath77_{i=1}^n$ ] , @xmath78_{i=1}^n$ ] , and @xmath79_{i=1}^n$ ] . @xmath80    in the second transmit - slot , the received signal at the desired receiver , i.e. , bob , is given by @xmath81 where @xmath82_{i=1}^n$ ] comprises complex channels from the @xmath11th relay to the rx and @xmath83 is the corresponding receiving awgn . by substituting into , @xmath84 can be expressed as @xmath85 the received signal at the @xmath86th eve , @xmath87 , is given by @xmath88 where @xmath89_{i=1}^n$ ] denotes the complex channels from the relays to the @xmath86th eve and @xmath90 is the awgn at the @xmath86th eavesdropper . the mutual information for the rx ( bob ) is given by @xmath91 , and that for the @xmath86th eve is @xmath92 , @xmath93 , where @xmath94 and @xmath95 , which denote their respective signal - to - interference - plus - noise ratios ( sinrs ) , are given at the top of next page . next , we define the secrecy rate region that consists of all the achievable secrecy rate for the relay wiretap channel given transmit power @xmath35 , denoted by @xmath96 , which is given by @xcite @xmath97 in this section , we consider the secrecy rate maximization problem by jointly optimizing the an beams , relay beam and their power allocations for weh - enabled af relays operating with sps , i.e. , @xmath30 , @xmath24 , is fixed .    by replacing @xmath39 with , @xmath98 in at the next page becomes @xmath99 and @xmath100_i\triangleq h_{sr_i}h_{r_id}\sqrt{\tfrac{\eta\bar\alpha_i(1-\bar\alpha_i)|h_{sr_i}|^2p_s}{(1-\bar\alpha_i)(|h_{sr_i}|^2p_s + \sigma_{n_a}^2)+\sigma_{n_c}^2}}.\ ] ] in addition , @xmath101 and @xmath102 in can also be combined as @xmath103_{i , i},\ ] ] where @xmath104_{i , i}=\frac{\eta\bar\alpha_ip_s|h_{sr_i}|^2|h_{r_id}|^2((1-\bar\alpha_i)\sigma_{n_a}^2+\sigma_{n_c}^2 ) } { ( 1-\bar\alpha_i)(|h_{sr_i}|^2p_s+\sigma_{n_a}^2)+\sigma_{n_c}^2}.\end{aligned}\ ] ] as a result , @xmath105 can be rewritten as @xmath106 where @xmath107_{i=1}^n$ ] . similarly by letting @xmath108_i\triangleq h_{sr_i}h_{r_ie , k}\sqrt{\frac{\eta\bar\alpha_i(1-\bar\alpha_i)|h_{sr_i}|^2p_s}{(1-\bar\alpha_i)(|h_{sr_i}|^2p_s + \sigma_{n_a}^2)+\sigma_{n_c}^2}}\ ] ] and @xmath109_{i , i}\triangleq\frac{\eta\bar\alpha_ip_s|h_{sr_i}|^2|h_{r_ie , k}|^2((1-\bar\alpha_i)\sigma_{n_a}^2+\sigma_{n_c}^2 ) } { ( 1-\bar\alpha_i)(|h_{sr_i}|^2p_s+\sigma_{n_a}^2)+\sigma_{n_c}^2},\ ] ] we have @xmath110    by some simple manipulation , is reformulated as a per - relay jamming power constraint given by @xmath111    now , the secrecy rate maximization problem w.r.t . @xmath112 s , @xmath44 s and @xmath47 for sps - based relays can be formulated as @xmath113    @xmath114      here , we consider the secrecy rate maximization problem for weh - enabled af relays with adjustable ps ratios @xmath115 by jointly optimizing the an beams , relay beam , weh ps ratios @xmath115 , and an ps ratios @xmath116 . first , consider the following variable transformation : @xmath117 using this , @xmath118 can then be expressed as @xmath119 , where @xmath120_{i=1}^n$ ] and @xmath121_{i=1}^n$ ] . moreover , @xmath122 and @xmath123 can be simplified as @xmath124 and @xmath125 , respectively , where @xmath126_{i=1}^n$ ] with @xmath127 , @xmath24 , and @xmath128_{i=1}^n$ ] . similarly , we have @xmath129 , @xmath130 , and @xmath131 , where @xmath132_{i=1}^n$ ] , @xmath133 . then we apply the above transformation to @xmath134 ( c.f .  ) and @xmath135 ( c.f .  ) , @xmath93 , to get and ( see next page ) . now , we can recast the constraints w.r.t . @xmath47 , @xmath22 s , and @xmath112 s to those w.r.t . the transformed variables @xmath136 s and @xmath137 s . in accordance with , the optimization variables , @xmath22 s and @xmath112 s , can be alternatively given by @xmath138 where @xmath139 . replacing @xmath22 s and @xmath112 s with , is reformulated as @xmath140 on the other hand , since @xmath141 and @xmath142 , @xmath24 , after some simple manipulation , it follows from that @xmath143    as such , the secrecy rate maximization problem for dps - based relays becomes @xmath144 in this section , we resort to centralized approaches to solve problem @xmath145 and @xmath146 , respectively , assuming that there is a central optimizer that is able to collect all csis including @xmath147 , @xmath148 and @xmath149 , perform the optimization , and broadcast to relays their individual optimized parameters .      to start with , we recast @xmath145 into a two - stage problem by introducing a slack variable @xmath150 . first of all , we solve the epigraph reformulation of @xmath145 with a fixed @xmath151 $ ] as @xmath152    defining @xmath153 as the optimum value of @xmath154 and denoting @xmath155 , the objective function of @xmath145 is given by @xmath156 where @xmath157 in the objective function has been omitted and we claim a zero secrecy rate if admits a negative value . as a result , @xmath145 can be equivalently given by @xmath158 note that this single - variable optimization problem allows for simple one - dimension search over @xmath159 $ ] , assuming that @xmath160 is attainable given any @xmath150 in this region . as the physical meaning of @xmath161 in @xmath154 can be interpreted as the maximum permitted sinr for the best eavesdropper s channel , feasibility for a non - zero secrecy rate implies that @xmath162 where cauchy - schwarz inequality has been applied in @xmath163 and @xmath164 follows from @xmath165 , @xmath166 . the above epigraph reformulation of non - convex problems like @xmath145 has been widely employed in the literature @xcite , and @xmath167 admits the same optimal value as @xmath145 while @xmath154 with the optimal @xmath150 provides the corresponding optimal solution to @xmath145 . we summarize the steps for solving @xmath145 here : given any @xmath159 $ ] , solve @xmath154 to obtain @xmath160 ; solve @xmath167 via a one - dimensional search over @xmath150 . before developing solutions to @xmath154 , we have the lemma below . @xmath160 is a concave function of @xmath150 . [ lemma : concavity of h_1(tau ) ]    see appendix  [ appendix : proof of concavity of h_1(tau ) ] .    using lemma [ lemma : concavity of h_1(tau ) ] , it is easy to verify that @xmath168 is also a concave function of @xmath150 according to the composition rule @xcite , which allows for a more effective search for the optimum @xmath150 , e.g. , bi - section method , than the exhaustive search used in @xcite . moreover , although @xmath160 is not differentiable w.r.t . @xmath150 , the bi - section method can still be implemented , the algorithm involving which is similarly applied in solving @xmath146 and thus will be given later in section  [ sec : proposed solutions for dps ] .    in the sequel , we focus on solving @xmath154 . by introducing @xmath169 and ignoring the rank - one constraint on @xmath170 , @xmath171 can be alternatively solved by @xmath172 note that the objective function has been multiplied by @xmath150 compared with that of @xmath154 for ready computation of @xmath160 . although @xmath173 is made easier to solve than @xmath154 by rank relaxation , it is still a quasi - convex problem considering the linear fractional form of the objective function and constraints , for which charnes - cooper transformation @xcite will be applied for equivalent convex reformulation . specifically , by substituting @xmath174 and @xmath175 into @xmath173 , it follows that @xmath176 problem @xmath177 can now be optimally and efficiently solved using interior - point based methods by some off - the - shelf convex optimization toolboxes , e.g. , cvx @xcite . we have the following results :    1 . the optimal solution to @xmath177 satisfies @xmath178 ; 2 . @xmath179 , where @xmath180 is given by @xmath181 in which @xmath182 is given in appendix  [ appendix : proof of rank - one solution with cj for sps ] ; 3 . @xmath183 . [ prop : rank - one solution with cj for sps ]    see appendix  [ appendix : proof of rank - one solution with cj for sps ] . proposition  [ prop : rank - one solution with cj for sps ] implies that the rank - one relaxation of @xmath173 from @xmath154 is tight for an arbitrary given @xmath150 . the @xmath184 s and @xmath185 s can thus be retrieved from the magnitude and angle of @xmath186 , respectively , by applying evd to @xmath187 . similar to section  [ sec : optimal solutions for sps ] , in this section , we aim at solving the two - stage reformulation of @xmath146 by introducing a slack variable @xmath188 $ ] . first , for a given @xmath150 , we solve @xmath189 next , denoting @xmath190 by @xmath191 , where @xmath192 is the optimum value for problem @xmath193 , we solve the following problem that attains the same optimum value as @xmath146 : @xmath194 where @xmath195 is similarly derived as @xmath196 so that we directly arrive at @xmath197 , denoted by @xmath195 . we claim that @xmath198 can be solved by bi - section for @xmath150 over the interval @xmath199 $ ] assuming that @xmath191 is valid for any given @xmath150 ( otherwise a zero secrecy rate , i.e. , @xmath200 , is returned . ) , since @xmath191 has the following property . @xmath191 is a concave function of @xmath150 . [ lemma : concavity of h_2(tau ) ]    the proof is similar to that for lemma  [ lemma : concavity of h_1(tau ) ] , and thus is omitted . it is also seen that how to attain @xmath191 forms the main thrust for solving @xmath146 . however , the constraints in , and are not convex w.r.t . @xmath136 and/or @xmath137 , @xmath24 , due to their high orders and multiplicative structure . @xmath193 thus turns out to be very hard to solve in general . to cope with these non - convex constraints , we introduce the following lemma .    the restricted hyperbolic constraints which have the form @xmath201 , where @xmath202 , @xmath203 , are equivalent to rotated second - order cone ( soc ) constraints as follows . @xmath204[lemma : hyperbolic to soc ]    for convenience , denoting @xmath205 , @xmath206 , @xmath207 by @xmath208 , @xmath209 , and @xmath210 , respectively , @xmath24 , can be rewritten as @xmath211 according to and , it is easily verified that @xmath212 . hence , is eligible for lemma  [ lemma : hyperbolic to soc ] , which is reformulated into the soc constraint : @xmath213 similarly , can be simplified as @xmath214 , and after some manipulation , it is recast into a constraint of the restricted hyperbolic form as @xmath215 is thus , in line with lemma  [ lemma : hyperbolic to soc ] , equivalent to an soc constraint given by @xmath216 at last , is a linear constraint w.r.t . @xmath208 and @xmath209 given by @xmath217    @xmath218    note that  have so far been equivalently transformed into the soc constraints , the linear constraints , as well as , the latter two of which are jointly convex w.r.t . @xmath208 and @xmath209 , @xmath24 . however , is still not convex w.r.t . @xmath210 , @xmath24 , yet . to circumvent this , in the sequel we propose to solve problem @xmath146 by alternating optimization . the upshot of the algorithm is that first we fix @xmath47 by @xmath219 and thus @xmath210 by @xmath220 , @xmath24 , and solve problem @xmath221 ( @xmath193,@xmath198 ) with fixed @xmath47 as @xmath221 ( @xmath222,@xmath223 ) in the sequel . ] to find the optimal @xmath224 , @xmath225 and @xmath226 via @xmath222 and @xmath223 ; then with @xmath227 , @xmath24 , we devise the optimal solution derived in section  [ sec : optimal solutions for sps ] to obtain the optimal cj covariance , viz @xmath228 , and thus @xmath229 , @xmath24 ; finally , by updating @xmath230 and @xmath231 , @xmath24 , problems @xmath221 and @xmath145 are iteratively solved until they converge . the remaining challenges lie in solving problem @xmath222 now that , and are all made convex w.r.t . their variables @xmath208 , @xmath209 , @xmath24 . similar to that for @xmath154 , we introduce @xmath232 and @xmath233 and exempt problem @xmath222 from @xmath234 and @xmath235 as follows : @xmath236 @xmath191 ( c.f .  ) and are given at the top of next page . recalling the procedure to deal with @xmath173 , we now apply charnes - cooper transformation to convert @xmath237 into a convex problem , denoted by @xmath238 , by replacing @xmath239 and @xmath240 with @xmath241 and @xmath242 , respectively . the solution for @xmath238 is tight and characterized by the following proposition . we have the following results :    1 . the optimal solution to @xmath238 satisfies @xmath243 such that @xmath244 ; 2 . @xmath245 is given by @xmath246 where @xmath247 , @xmath248 , and @xmath249 is given in appendix  [ appendix : proof of rank - one solution with cj for dps ] ; 3 . @xmath250 , of which the diagonal entries compose a vector denoted by @xmath251 , can be reconstructed by @xmath252 with rank one . [ prop : rank - one solution with cj for dps ]    see appendix  [ appendix : proof of rank - one solution with cj for dps ] . the @xmath253 s and @xmath112 s are thus attained according to via evd of @xmath254 and @xmath255 . the proposed algorithm for solving @xmath146 is presented in table  [ table : algorithm i ] .    ' '' ''    @xmath228 ; @xmath256 that denotes the optimum value for @xmath145 given @xmath257 , @xmath24 @xmath258 , @xmath259 @xmath260 @xmath261 , @xmath262 , @xmath24 , and solve @xmath221 : @xmath263 , @xmath264 , @xmath265 , @xmath266 , @xmath267 @xmath268 , @xmath269 solve @xmath222 and @xmath191 @xmath270 @xmath271 , where @xmath272 and @xmath273 denotes an arbitrary small value . @xmath274 @xmath275 @xmath276 , @xmath277 , and obtain @xmath278 according to @xmath279 , @xmath24 , and solve @xmath145 via @xmath154 and @xmath167 @xmath280 , @xmath281 , and obtain @xmath282 and @xmath283 according to @xmath284 , @xmath24 update @xmath285 according to @xmath278 , @xmath282 , @xmath283 , and @xmath228    ' '' '' in this section , we investigate heuristic algorithms to solve problems @xmath145 and @xmath146 in a completely distributed fashion . note that different from the paradigm of _ distributed optimization _ that allows for certain amount of information exchange based on which iterative algorithms are developed to gradually improve the system performance , we herein assume that each individual relay can only make decision based on its local csis , namely , @xmath34 , @xmath286 , @xmath287 , @xmath24 , and there is no extra means of information acquisition for ease of implementation . the purpose for such an algorithm is twofold : on one hand , we aim to answer the question that in the least favourable situation , namely , no coordination over the relays , how to improve the achievable secrecy rate of the system ? on the other hand , it provides a lower - bound for the centralized schemes proposed in section  [ sec : centralized ] , which sheds light upon the trade - off achievable between secrecy performance and complexity . besides , we emphasize the jamming scheme that is different from the cj in the centralized schemes . unlike the cj signal coordinately transmitted by all relays , in the distributed implementation , each relay is only able to generate its an locally , i.e. , @xmath288^t$ ] , in which @xmath289 s are i.i.d . an beams , denoted by @xmath290 . this type of cj is known to be uncoordinated with the covariance matrix given by @xmath291)$ ] . in this section , we assume that each relay consumes all of its remaining power from af for an , i.e. , @xmath292 , @xmath166 ( c.f .  ) . hence , the an design solely depends on @xmath22 s and/or @xmath112 s . first , we propose a heuristic scheme for the @xmath11th af relay to decide on @xmath112 , @xmath24 , which is given by @xmath293 where @xmath294 is a constant controlling the relay s level of jamming . for example , a larger @xmath295 indicates that each relay prefers to splitting a larger portion of power for jamming and vice versa . the intuition behind is that if the @xmath11th relay observes that @xmath296 , which means that a nonnegative secrecy rate is achievable even if there is only itself in the system , it will shut down the an ; otherwise , it will split up to @xmath295 portion of the harvested power for jamming . in an extreme case of @xmath297 , probably when an eve is located within the very proximity of this relay , it allocates the maximum permissible portion of power , i.e. , @xmath295 , for an . next , since an individual relay can not evaluate the secrecy performance of the whole system , @xmath44 s are simply chosen to be the optimum for the multi - af relaying without security considerations , i.e. , @xmath298 , @xmath24 .      following the same designs for @xmath112 s and @xmath44 s in section  [ sec : distributed alg . for sps ] , the remaining task for weh - enabled relays operating with dps is to set proper values for @xmath22 s . we choose @xmath22 s that maximize the `` hypothetical sinr '' . this `` hypothetical '' sinr may not be the actual sinr for the destination , but just a criterion calculated based on the `` hypothetical '' received signal seen by the @xmath11th relay , given by @xmath299 the corresponding sinr is thus expressed as @xmath300 where @xmath301 . consequently , the maximization of w.r.t . @xmath22 , @xmath24 , is formulated as @xmath302    the optimal @xmath22 , @xmath24 , to @xmath303 is @xmath304[prop : optimal alpha to distr ]    it is easy to verify that problem @xmath303 is convex and the minimum solution of its objective function derived from the first - order derivative happens to fall within the feasible region of @xmath22 , which is seen in .    with @xmath112 s , @xmath44 s and @xmath22 s set , each af relay is then able to decide its relay weight and an transmission . in this section we compare our proposed schemes operating with sps or dps with some benchmarks . in the centralized case , the optimal solution for sps in section  [ sec : optimal solutions for sps ] is denoted by _ cj - sps _ , while algorithm  [ table : algorithm i ] in section  [ sec : proposed solutions for dps ] is denoted by _ cj - dps_. the distributed schemes in section  [ sec : distributed alg . for sps ] and section  [ sec : distributed alg . for dps ] are referred to as _ distributed - sps _ and _ distributed - dps _ , respectively . to demonstrate the effectiveness of our an - aided secure multi - af relay beamforming algorithms , we also provide three benchmark schemes : _ nocj - sps _ , _ nocj - dps _ and _ random ps_. for _ nocj - sps _ , we solve problem @xmath145 by replacing @xmath47 with @xmath305 . similarly , for _ nocj - dps _ , we initialize @xmath306 and quit the loop in algorithm  [ table : algorithm i ] after the very first time of solving problem @xmath221 . _ random ps _ , on the other hand , picks up i.i.d . @xmath22 and @xmath112 uniformly generated over @xmath307 $ ] , respectively , and co - phases @xmath308 , @xmath24 . consider that @xmath19 weh - enabled af relays and @xmath309 eavesdroppers are located within a circular area of radius @xmath310 . specifically , we assume that their respective radius and radian are drawn from uniform distributions over the interval @xmath311 $ ] and @xmath312 , respectively . we also assume that the channel models consist of both large - scale path loss and small - scale multi - path fading . the unified path loss model is given by @xmath313 where @xmath314 , @xmath315 denotes the relevant distance , @xmath316 m is a reference distance , and @xmath317 is the path loss exponent set to be @xmath318 . @xmath34 , @xmath286 , and @xmath319 , @xmath166 , @xmath133 , are generated from independent rayleigh fading with zero mean and variance specified by . the simulation parameters are set as follows unless otherwise specified : the radius defining the range is @xmath320 m ; the transmit power at the source is @xmath321dbm ; the noise variances are set as @xmath322dbm , @xmath323dbm , @xmath324 , and @xmath325 , @xmath93 ; the eh efficiency is set to @xmath326 . also , numerical results were averages over @xmath327 independent channel realizations .      here , we evaluate the performance of the proposed centralized designs in section  [ sec : centralized ] . the efficiency of the alternating optimization that iteratively attains numerical solution to @xmath146 is studied in fig .  [ fig : alternating sec vs iter ] , which shows the increment of the achievable secrecy rate after each round of the iteration . the most rapid increase is observed after the first iteration , which illustrates that the optimization of the ps ratios , @xmath22 s , accounts for the main factor for the secrecy rate performance gains over a sps scheme that sets @xmath328 . it is seen that the alternating algorithm converges within the relative tolerance set to be @xmath329 , after an average of @xmath330@xmath331 iterations for several channel realizations , which appears reasonable in terms of complexity . dbm , @xmath332 , and @xmath333.,width=288 ]    fig . [ fig : centralized compare vs n ] shows the achievable secrecy rate for the legitimate rx with the increase in the number of af relays by different schemes . as we can see , the secure multi - af relaying schemes assisted by the transmission of an outperforms that without an for both sps and dps . in addition , with the increase in @xmath19 , the role of cj gradually reduces for both schemes of sps and dps this is because as @xmath19 gets larger , the optimal designs tend to suppress the interception at the most capable eavesdropper more effectively with @xmath19 dof , enforcing the numerator of @xmath334 to a relatively low level , which can also be observed from @xmath335 in fig . [ fig : subfig : eve s rate vs n ] , and therefore the optimal amount of power allocated to an beams inclines to be little ; otherwise the jamming yielded will be detrimental to the reception at the legitimate rx . besides , given the same number of af relays , the secrecy performance gains brought by the proposed schemes with cj are more prominent in the presence of more eavesdroppers , since it is hard to reduce all the eavesdroppers channel capacity without resorting to cj properly . db for @xmath332 and @xmath336 , respectively.,width=288 ]    fig . [ fig : centralized compare sec vs eves ] shows the achievable secrecy rate for the legitimate rx versus the number of eavesdroppers by different schemes . first , similar to the results shown in fig . [ fig : centralized compare vs n ] , the proposed an - aided multi - af relaying designs operating with dps - enabled relays , viz _ cj - dps _ , perform best among all the schemes . secondly , as @xmath309 goes up , the an - aided schemes , _ cj - dps _ and _ cj - sps _ , allow the secrecy rate to drop slowly , in other words , more robust against multiple eavesdroppers , while the secrecy rate of their nocj counterparts almost goes down linearly with @xmath309 . furthermore , with @xmath309 increasing , for example , more than @xmath337 , the increase in the number of relays , from @xmath332 to @xmath338 , can not replace the role of cj as shown in fig . [ fig : centralized compare vs n ] and [ fig : centralized compare sec vs eves ] , since in the presence of many eavesdroppers , more relays may also result in improved eavesdroppers decoding ability w / o the assistance of cj . it is also noteworthy that with @xmath339 , there is little use of cj by the centralized schemes , which was also observed in @xcite .     for @xmath332 and @xmath336 , respectively.,width=288 ] [ fig : centralized compare sec vs ps ] provides simulation results of different schemes by varying the source transmit power . it is seen that with more power available at the source , the advantage of cj is more pronounced , since given other variables fixed , larger @xmath35 indicates larger feasible regions for @xmath145 and @xmath146 . furthermore , as similarly seen in fig . [ fig : centralized compare vs n ] , with a mild number of eavesdroppers ( @xmath333 ) , subject to the same @xmath35 , a large number of cooperative relays enables more dof in designing the optimal @xmath22 s and @xmath44 s , which alleviates the dependence on an beams to combat eves .      here , we study the performance of the distributed schemes , namely , _ distributed - sps _ and _ distributed - dps _ in section  [ sec : distributed ] . as mentioned earlier , these heuristics are provided as benchmarks to demonstrate what can be done under the extreme `` no - coordination '' circumstance , in comparison with _ random ps_. note that any other distributed schemes with certain level of cooperation among relays are supposed to increase the secrecy performance up to the proposed centralized algorithms , namely , _ cj - dps _ and _ cj - sps _ , at the expense of extra computational complexity and system overhead . db and @xmath333.,width=288 ]    fig . [ fig : distributed compare sec vs n ] provides the results for the achievable secrecy rate of various schemes versus the number of relays . _ distributed - sps _ and _ distributed - dps _ , are observed to be outperformed by their centralized counterparts though , they are considerably superior to _ random ps_. it is also seen that the performance gap between the centralized and distributed approaches is enlarged as @xmath19 increases , which is expected , since larger @xmath19 yields more dof for cooperation that is exclusively beneficial for the centralized schemes . furthermore , compared with the centralized schemes , the distributed ones are more vulnerable to the increase in the eavesdroppers number . db and @xmath19=8.,width=288 ]    in fig . [ fig : distributed compare sec vs eves ] , we investigate the relationship between the secrecy rate performance and the number of eavesdroppers by different methods . as can be observed , compared with the centralized schemes , the secrecy rates achieved by _ distributed - sps _ and _ distributed - dps _ both reduce more drastically with the increase in @xmath309 due to the lack of effective cooperation . also , the advantage of dps over sps for the distributed schemes is compromised since @xmath22 s are not jointly designed with other parameters . at last , a similar observation has been made as that for fig . [ fig : distributed compare sec vs n ] , that is , larger @xmath19 yields more visible performance gap between the centralized and distributed approaches .     and @xmath333.,width=288 ]    in fig . [ fig : distributed compare sec vs ps ] , we examine the effect of increasing the transmit power at the source on the secrecy performance of different schemes under the same settings as those in fig . [ fig : centralized compare sec vs ps ] . among all the presented designs , _ cj - dps _ still achieves the best secrecy rate as observed in other examples . also , the fact that larger @xmath19 benefits more from cooperative designs is corroborated again due to the same reason as that for fig . [ fig : distributed compare sec vs n ] and [ fig : distributed compare sec vs eves ] . furthermore , the secrecy rate of _ distributed - sps _ or _ distributed - dps _ is quickly saturated when @xmath340db while that for their centralized counterparts still rises at fast speed . this paper studied secure communications using multiple single - antenna weh - enabled af relays assisted by an via cj for a swipt network with multiple single - antenna eavesdroppers . using ps at the relays , the achievable secrecy rates for the relay wiretap channel were maximized by jointly optimizing the cb and the cj covariance matrix along with the ps ratios for relays operating with , respectively , sps and dps . for dps , reformulating the constraints into restricted hyperbolic forms enabled us to develop convex optimization - based solutions . further , we proposed an information - exchange - free distributed algorithm that outperforms the random decision . as the optimal solution to @xmath177 is proved to be optimal to problem @xmath154 , the optimum value for @xmath177 is @xmath341 . hence , we identify the property of @xmath160 by investigating @xmath177 , with its lagrangian given by @xmath342 where @xmath343 denotes a tuple consisting of all the primal and dual variables . specifically , @xmath344 , @xmath345 and @xmath346 are lagrangian multipliers associated with @xmath347 , @xmath348 and the first constraint of @xmath177 , respectively ; @xmath349 are the dual variables associated with the sinr constraint for the @xmath86th eve , respectively ; @xmath350_{i=1}^n)$ ] with each diagonal entry @xmath351 denoting the dual variable associated with the per - relay power constraint ; @xmath352 is the lagrangian multiplier associated with @xmath353 . in addition , @xmath354_{i=1}^n)$ ] . the karush - kuhn - tucker ( kkt ) conditions for are partially listed as follows :      the associated dual problem is accordingly given by @xmath356 where @xmath344 , @xmath345 , and @xmath352 are given by  , respectively . since it is easily verified that @xmath177 satisfies the slater condition , the strong duality holds @xcite . this implies that the dual optimum value in @xmath357 is @xmath160 , which turns out to be a point - wise minimum of a family of affine functions and thus concave over @xmath159 $ ] @xcite . pre- and post - multiplying @xmath359 with the left hand side ( lhs ) and right hand side ( rhs ) of , respectively , and substituting @xmath360 into , @xmath361 can be rewritten as @xmath362 introducing the notation of @xmath5_{\rm offd}$ ] to represent a square matrix with its diagonal entries removed , it follows from   that @xmath363_{\rm offd}=\boldsymbol 0.\ ] ] by subtracting from , @xmath344 can be rewritten as @xmath364_{\rm d}\notag\\ & + \left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol y_2\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}+\sum_{k=1}^k\theta_k p_s\tilde{\boldsymbol h}_{se , k}^\dag\tilde{\boldsymbol h}_{se , k}^t,\label{eq : y1 intermediate2}\end{aligned}\ ] ] where @xmath5_{\rm d}$ ] denotes a square matrix with only the diagonal remained . observing that @xmath365_{\rm d}^{-1}\left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol h_{re , k}^\dag\boldsymbol h_{re , k}^t\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}\nonumber\\ & = { \rm diag}\left([|h_{r_ie , k}|^2\left.\middle/\right.|h_{r_id}|^2]_{i=1}^n\right)\notag\\ & \equiv \boldsymbol r_{ed , k},\label{eq : r_{ed , k}}\end{aligned}\ ] ] @xmath361 can be finally recast as @xmath366 where @xmath367_{\rm d}-\left(\left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol h_{rd}^\dag\boldsymbol h_{rd}^t\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}-\boldsymbol d_{\hat{sd}}\right ) \left(\lambda\boldsymbol i-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)\boldsymbol r_{ed , k}\right).\end{aligned}\ ] ]    in the following , we show that @xmath368 is a positive definite matrix . note that since @xmath369 is a diagonal matrix , its definiteness is only determined by the signs of its diagonal entries , for which we commence with the discussion in three difference cases below . 1 .   * case i * : @xmath370 such that @xmath371_{i , i}<0 $ ] . since @xmath372_{\rm d}-\boldsymbol d_{\hat{sd}}\right]_{i , i}=$ ] + @xmath373 , it follows from that @xmath374_{i , i}>0 $ ] in this case . * case ii * : @xmath370 such that @xmath371_{i , i}>0 $ ] . we have @xmath375_{\rm d}\right]_{i , i}-[\boldsymbol w_0]_{i , i}\vert h_{r_id}\vert^2(\lambda^\ast-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)[\boldsymbol r_{ed , k}]_{i , i})\ge 0 $ ] in accordance with , which implies that @xmath374_{i , i}=\left[\left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol y_2\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}\right]_{i , i}-[\boldsymbol w_0]_{i , i}\vert h_{r_id}\vert^2(\lambda-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)[\boldsymbol r_{ed , k}]_{i , i})+[\boldsymbol d_{\hat{sd}}]_{i , i}(\lambda-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)[\boldsymbol r_{ed , k}]_{i , i})>0 $ ] ( c.f .  ) . * case iii * : @xmath370 such that @xmath371_{i , i}=0 $ ] . in this case , it follows that @xmath374_{i , i}= \left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol y_2\boldsymbol w_0^{\frac{1}{2}}\right]_{i , i}$ ] + @xmath376 . it is noteworthy that the number of @xmath11 s such that @xmath371_{i , i}=0 $ ] can not exceed one . this can be proved by contradiction as follows . if @xmath377 , @xmath378 , such that @xmath371_{i_1,i_1}=0 $ ] and @xmath371_{i_2,i_2}=0 $ ] , it implies that @xmath379_{i_1,i_1}=\sum_{k=1}^k\theta_k[\boldsymbol r_{ed , k}]_{i_2,i_2}$ ] , which contradicts to the fact that for any two independent continuously distributed random variables , the chance that they are equal is zero .    in summary , @xmath374_{i , i}\ge 0 $ ] , @xmath24 . if @xmath374_{i , i } > 0 $ ] , then it is obvious that @xmath380 . next , we show that it still holds true in the case that @xmath381 , such that @xmath374_{i^\prime , i^\prime}=0 $ ] , @xmath382 , by definition . define the null - space of @xmath369 by @xmath383 and multiply @xmath384 and @xmath385 , @xmath386 , on the lhs and rhs of @xmath368 , respectively . if @xmath387 , it is straightforward to obtain @xmath388 ; otherwise , it follows that @xmath389_{i^\prime}\vert^2>0 $ ] , as @xmath390_{i^\prime}\neq 0 $ ] in probability . this completes the proof . as a result , @xmath344 in is shown to always take on a special structure , that is , a full - rank matrix minus a rank - one matrix . note that this observation plays a key role in proving the rank - one property of @xmath391 , which is also identified in ( * ? ? ? * appendix c ) .      finally , multiplying both sides of by @xmath391 , as per , we obtain the following equation : @xmath393 which further implies that @xmath394 . in addition , since the optimality of @xmath177 suggests that @xmath395 , @xmath396 is thus proved . as @xmath391 can be decomposed as @xmath397 by evd , results in @xmath398 , which further implies that @xmath399 therefore , admits a unique solution @xmath182 up to a scaling factor , which is given by @xmath400 consequently , we have @xmath401 , where @xmath402 . on the other hand , by plugging @xmath401 into the equality constraint of @xmath177 , we have @xmath403 , which yields @xmath404    at last , we show @xmath405 of proposition  [ prop : rank - one solution with cj for sps ] . for the case of @xmath406 , it is obvious that @xmath407 . for the case of @xmath408 , only a sketch of the proof is provided here due to the length constraint . according to , first it is provable that @xmath409 is a full - rank matrix when @xmath177 obtains its optimum value ; next , observing that @xmath410 , it follows that @xmath411 as a result of @xmath412 ; then according to , @xmath413 is thus obtained .          for convenience of the proof , the optimum value for @xmath237 and @xmath415 are denoted by @xmath416 and @xmath417 , respectively . assuming that @xmath418 is the optimal solution to @xmath237 , it is easily verified to be feasible for @xmath415 as well , which implies that @xmath419 . on the other hand , if @xmath415 returns an optimal solution of @xmath420 , by defining @xmath421 and @xmath422 , @xmath24 , we show that @xmath423 is also feasible for @xmath237 as follows . as for , @xmath424 where @xmath163 is due to , and @xmath164 comes from @xmath425 and @xmath426 . similarly , @xmath427 + @xmath428 can be proved to satisfy as well . in addition , @xmath429 , @xmath24 , i.e. , holds true . these feasibility implies that @xmath430 . by combining the above two facts , we have @xmath431 and complete the proof .    then , we apply the charnes - cooper transformation again to @xmath415 , the result of which is denoted by @xmath432 , to study the property of @xmath433 and @xmath250 . it is noteworthy that the charnes - cooper transformed constraint of admits the form given by @xmath434 , @xmath24 , where @xmath435 is the column vector inside @xmath7 of the lhs of and @xmath436 indicates the rhs expression . since it is easily checked that @xmath437 as a result of feasibility , we have @xmath438 , which implies that @xmath439 according to schur complement . it thus follows that holds true , @xmath440 : @xmath434 . as such , we can show that can be recast into a constraint without @xmath441 , i.e. , not related to @xmath442 , @xmath443 , following the same procedure as ( * ? ? ? * appendix iii ) by exploiting ( * ? ? ? * _ lemma 2 _ ) . similarly , the charnes - cooper transformed constraint of can also be rewritten without @xmath442 , @xmath443 . hence , the partial lagrangian for @xmath432 in terms of @xmath442 and @xmath443 can be expressed as ( [ eq : lagrangian of ( p2.1-sdp - eqv ) ] ) ( see top of next page ) ,      where @xmath445 denotes a tuple comprising all the associated primal and dual variables : @xmath344 , @xmath345 , and @xmath349 are lagrangian multipliers associated with @xmath442 , @xmath443 , and , @xmath93 , respectively ; @xmath346 is the dual variable associated with the only equality constraint ; @xmath446_{i=1}^n)$ ] and @xmath447_{i=1}^n)$ ] denote those associated with @xmath448 and @xmath449 , @xmath24 , respectively ; finally , the diagonal entry of @xmath450_{i=1}^n)$ ] denotes the dual variable associated with @xmath451 , @xmath24 . the kkt conditions related to are accordingly given by        1 . * case i * : @xmath456 , @xmath133 . in this case , since @xmath457 ( c.f .  ) due to the strong duality , it is easily verified that @xmath458 and therefore @xmath459 . * case ii * : @xmath460 such that @xmath461 . according to , @xmath462 , which is a positive semidefinite diagonal matrix . according to the similar argument made in * case iii * of appendix  [ appendix : proof of rank - one solution with cj for sps ] , it is thus shown that @xmath249 has maximum one zero diagonal entry and the positive definiteness of @xmath463 can be proved by definition without difficulty .    as @xmath344 ( c.f .  ) again complies with the difference between a positive definite matrix and a rank one matrix , i.e. , @xmath464 , it turns out that @xmath465 according to . then , following the same procedure as that in appendix  [ appendix : proof of rank - one solution with cj for sps ] , @xmath466 of proposition  [ prop : rank - one solution with cj for dps ] can be proved ( details omitted for brevity ) . finally , it is verified that @xmath238 is related to @xmath250 merely with its diagonal entries , viz @xmath467_{i , i}$ ] , @xmath166 , ( c.f .  , ) . furthermore , denoting @xmath468_{i , i}^{1/2}]_{i=1}^n$ ] by @xmath251 , it is easily checked that the diagonal entries remain the same after we replace @xmath250 with @xmath252 . combining the above two facts , we arrive at the conclusion that such modification does not affect its optimality while returning a rank - one @xmath250 for @xmath238 , which completes the proof for @xmath405 . e. tekin , s. member , and a. yener , `` the general gaussian multiple - access and two - way wiretap channels : achievable rates and cooperative jamming , '' _ ieee trans . inf . theory _ 54 , no . 6 , pp . 27352751 , jun . 2008 . j. zhang , m. c. gursoy , `` relay beamforming strategies for physical - layer security , '' in _ proc . ieee annual conference on information sciences and systems ( ciss10 ) _ , princeton , usa , mar . 2010 , pp.16 . y. yang , q. li , w .- k . ma , j. ge , and p. c. ching , `` cooperative secure beamforming for af relay networks with multiple eavesdroppers , '' _ ieee signal process . 1 , pp . 3538 , jan . 2013 . h. xing , k .- k wong , z. chu , and a. nallanathan , `` to harvest and jam : a paradigm of self - sustaining friendly jammers for secure af relaying , '' _ to appear in ieee trans . signal process . _ , available on - line at arxiv:1502.07066v3 . z. ding , k. k. leung , d. l. goeckel , and d. towsley , `` opportunistic relaying for secrecy communications : cooperative jamming vs. relay chatting , '' _ ieee trans . 10 , no . 6 , pp . 17251729 , jun . 2011 . y. yang , q. li , w .- k . ma , j. ge , and m. lin , `` optimal joint cooperative beamforming and artificial noise design for secrecy rate maximization in af relay networks , '' in _ proc . ieee workshop on signal processing advances in wireless communications ( spawc ) _ , darmstadt , de , jun . 2013 , pp.360364 . q. li , y. yang , w .- k . ma , m. lin , j. ge , and j. lin ,  robust cooperative beamforming and artificial noise design for physical - layer secrecy in af multi - antenna multi - relay networks , " _ ieee trans . signal process . 1 , pp . 206220 , jan . k. ishibashi , `` dynamic harvest - and - forward : new cooperative diversity with rf energy harvesting , '' in _ proc . ieee international conference on wireless communications and signal processing ( wcsp ) _ , hefei , cn , oct . 2014 , pp . 15 . a. nasir , x. zhou , s. durrani , and r. kennedy , `` relaying protocols for wireless energy harvesting and information processing , '' _ ieee trans.wireless commun . 36223636 , jul . z. ding , s. perlaza , i. esnaola , and h. poor , `` power allocation strategies in energy harvesting wireless cooperative networks , '' _ ieee trans . wireless commun . _ , vol . 13 , no . 2 , pp . 846860 , feb . 2014 . s. timotheou , i. krikidis , g. zheng , and b. ottersten , `` beamforming for miso interference channels with qos and rf energy transfer , '' _ ieee trans . wireless commun . _ , vol.13 , no.5 , pp . 26462658 , may 2014 . m. zhao , x. wang , and s. feng , `` joint power splitting and secure beamforming design in the multiple non - regenerative wireless - powered relay networks , '' _ ieee commun . 9 , pp . 15401543 , sept . 2015 . d. ng , e. lo , and r. schober , `` robust beamforming for secure communication in systems with wireless information and power transfer , '' _ ieee trans . wireless commun . _ , vol . 13 , no . 8 , pp 45994615 , aug . 2014 . q. li , q. zhang , and j. qin , `` secure relay beamforming for simultaneous wireless information and power transfer in non - regenerative relay networks , '' _ ieee trans . 5 , pp . 24622467 , jun . 2014 . q. li and w .- k . ma , `` spatially selective artificial - noise aided transmit optimization for miso multi - eves secrecy rate maximization , '' _ ieee trans . signal process . 27042717 , may 2013 .
Abstract:

Answer: this paper studies secrecy transmission with the aid of a group of wireless energy harvesting ( weh)-enabled amplify - and - forward ( af ) relays performing cooperative jamming ( cj ) and relaying . the source node in the network does _ simultaneous wireless information and power transfer ( swipt ) _ with each relay employing a power splitting ( ps ) receiver in the first phase ; each relay further divides its harvested power for forwarding the received signal and generating artificial noise ( an ) for jamming the eavesdroppers in the second transmission phase . in the centralized case with global channel state information ( csi ) , we provide closed - form expressions for the optimal and/or suboptimal af - relay beamforming vectors to maximize the achievable secrecy rate subject to individual power constraints of the relays , using the technique of semidefinite relaxation ( sdr ) , which is proved to be tight . a fully _ distributed _ algorithm utilizing only local csi at each relay is also proposed as a performance benchmark . simulation results validate the effectiveness of the proposed multi - af relaying with cj over other suboptimal designs .    artificial noise , cooperative jamming , amplify - and - forward relaying , secrecy communication , semidefinite relaxation , wireless energy harvesting .    [ section ] [ section ] [ section ] [ section ] [ proposition ] [ section ] [ section ] [ section ]