Next Article in Journal
Hands-Free User Interface for AR/VR Devices Exploiting Wearer’s Facial Gestures Using Unsupervised Deep Learning
Next Article in Special Issue
Trajectory Planning for Data Collection of Energy-Constrained Heterogeneous UAVs
Previous Article in Journal
Development of a Cost-Effective Optical Sensor for Continuous Monitoring of Turbidity and Suspended Particulate Matter in Marine Environment
Previous Article in Special Issue
Completion Time Minimization for Multi-UAV Information Collection via Trajectory Planning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aerial Cooperative Jamming for Cellular-Enabled UAV Secure Communication Network: Joint Trajectory and Power Control Design

1
The Glorious Sun School of Business and Management, Donghua University, Shanghai 200051, China
2
College of Information Science & Technology, Chengdu University of Technology, Chengdu 610059, China
3
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
4
School of Electronic Information Engineering, Shanghai Dianji University, Shanghai 201306, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(20), 4440; https://doi.org/10.3390/s19204440
Submission received: 22 September 2019 / Revised: 10 October 2019 / Accepted: 11 October 2019 / Published: 14 October 2019
(This article belongs to the Special Issue Optimization and Communication in UAV Networks)

Abstract

:
To improve the secrecy performance of cellular-enabled unmanned aerial vehicle (UAV) communication networks, this paper proposes an aerial cooperative jamming scheme and studies its optimal design to achieve the maximum average secrecy rate. Specifically, a base station (BS) transmits confidential messages to a UAV and meanwhile another UAV performs the role of an aerial jammer by cooperatively sending jamming signals to oppose multiple suspicious eavesdroppers on the ground. As the UAVs have the advantage of the controllable mobility, the objective is to maximize the worst-case average secrecy rate by the joint optimization of the two UAVs’ trajectories and the BS’s/UAV jammer’s transmit/jamming power over a given mission period. The objective function of the formulated problem is highly non-linear regarding the optimization variables and the problem has non-convex constraints, which is, in general, difficult to achieve a globally optimal solution. Thus, we divide the original problem into four subproblems and then solve them by applying the successive convex approximation (SCA) and block coordinate descent (BCD) methods. Numerical results demonstrate that the significantly better secrecy performance can be obtained by using the proposed algorithm in comparison with benchmark schemes.

1. Introduction

Due to the many advantages of controllable mobility, such as on-demand fast deployment, wide coverage, low cost, and line-of-sight (LoS) transmission that offers good channel capacity, unmanned aerial vehicles (UAVs) have been extensively utilized in different scenarios, e.g., surveillance and monitoring [1,2,3], search and rescue [4,5], cargo transportation [6], data collection [7] and mobile relays [8].
Recently, UAVs have attracted increasing attention in wireless communications, and are anticipated to playing an important role in the next-generation wireless networks [9,10]. Generally, there are two promising solutions to UAV communication applications: cellular-enabled UAV communication (CEUC) and UAV-assisted terrestrial communication (UATC) networks [11]. In UATC, the UAVs are flexibly deployed as aerial base stations (BSs) or mobile relays to assist in providing reliable communication services for terrestrial networks [12,13,14]. By contrast, the UAVs are integrated into the wireless network scenarios as aerial users served by ground BSs in the CEUC system [15]. Owing much to the almost ubiquitous accessibility of the existing LTE (Long Term Evolution) and the forthcoming (beyond) fifth-generation ((B)5G) cellular networks, reliable communications can be supported between UAVs and their corresponding BSs [16,17]. The CEUC is anticipated to have a number of appealing advantages over the existing ground-to-UAV communications, including the ease of monitoring and management, ubiquitous accessibility, robust navigation and enhanced performance, etc. [11]. Despite its merits, the UAV communication based on the future (B)5G cellular networks is more susceptible to suspicious eavesdropping on the ground, which leads to a severe security challenge that is urged to be solved.
Currently, the UAV trajectory design combined with physical-layer security techniques, as a promising solution, has drawn significant attention to safeguard the UAV communication. Specifically, UAV secure communications are studied in [18,19], where the average secrecy rate is significantly improved via optimizing the trajectory of the UAV jointly with the power control for a finite mission duration. As cooperative jamming is one of the important physical-layer techniques that can enhance the secrecy performance, reference [20] proposed to employ a UAV as a friendly mobile jammer, to ensure the secrecy of the ground wiretap channel. In [21], a novel full-duplex operation was applied to the rotary-wing UAV to further improve the energy efficiency (EE) of UAV secrecy communications, and the EE was maximized by the joint optimization of the source transmit/UAV’s jamming power and UAV trajectory. A four-node mobile relay and eavesdropper system is proposed in [22], where the UAV was employed as a mobile relay to assist in terrestrial communications. To cope with the non-convex secrecy rate maximization problem, an alternating optimization algorithm is designed by optimizing the power control and UAV trajectory alternatively. The authors in [23,24] proposed a dual-UAV UATC network to enhance the communication quality and improve the secrecy performance, where the downlink transmission from the UAVs is established by adaptively adjusting the UAVs’ trajectories and transmit powers. Note that most of the above studies only focus on the security issues in UATC systems. However, how to design efficient anti-eavesdropping methods, to protect legitimate BS-to-UAV transmission in the CEUC networks, has not been investigated, and thus remains a challenging problem to address.
In light of the above, we propose an anti-eavesdropping scheme by employing an aerial UAV jammer in the CEUC network, where one UAV flies to receive confidential messages from a BS while the mobile UAV jammer confuses multiple suspicious eavesdroppers on the ground by sending jamming signals. Specifically, we take into account the joint optimization of both the UAVs’ trajectories and the BS’s/UAV jammer’s power allocation, in order to maximize the average worst-case secrecy rate of the UAV receiver for a given finite period. In the proposed scheme, the UAVs are subject to the practical mobility as well as both the average and peak power constraints. In contrast with the above-mentioned existing works, the UAVs’ trajectory design in our proposed CEUC network is particularly important, as the interference from other UAVs cannot be practically cancelled, which causes different objective function and constraints. Therefore, well-designed trajectories of the UAVs can not only avoid severe interference between UAVs, but also provide effective jamming signals to the eavesdroppers, which is expected to notably enhance the secrecy performance. As the formulated optimization problem is non-convex with the objective function as well as its constraints, it is very hard to obtain a globally optimal solution (Since the difficulty of the original problem is NP-hard, it is generally impossible to obtain the globally optimal solution by using the present optimization techniques.) To tackle this challenging problem, we first transform it into a lower bound expression with more tractability. Then, an efficient algorithm is designed by applying the block coordinate descent (BCD) method [25,26]. To be specific, we partition the total optimization variables into four blocks for the two UAVs’ trajectories, BS’s transmit power, and UAV jammer’s jamming power control, respectively. Then, each block is alternatively optimized in each iteration with other blocks being fixed. Although we fix the other three blocks, the corresponding optimization problem remains intractable because of its non-convex. To obtain a high-quality approximately optimal solution, we thus introduce a series of slack variables and apply successive convex approximation (SCA) technique [27,28]. The proposed algorithm has the applicable complexity and guarantees to converge to a locally optimal solution to this problem. To best of knowledge, this is the first work that exploits the anti-eavesdropping UAV trajectory design to solve physical-layer security issues of the CEUC system. The numerical results illustrate that the designed algorithm achieves significantly better secrecy performance than all benchmarks without trajectory or power control design, especially the scheme without the UAV jammer, as in [18].
The rest of this paper is organized as follows. Section 2 gives the system model and problem formulation. In Section 3, a joint optimization algorithm is proposed and its complexity and convergence performance are also analyzed. The simulations are presented in Section 4 to verify the effectiveness of the proposed algorithm. Finally, Section 5 concludes the paper.

2. System Model and Problem Formulation

2.1. System Model

Consider a CEUC network, as shown in Figure 1, where a ground BS transmits confidential messages to a mobile UAV receiver (denoted by U) within a given UAV flight period T, while I malicious eavesdroppers on the ground, denoted by E i for i I { 1 , , I } , intercept the messages from the valid UAV communication. To safeguard the legitimate transmission, the potential eavesdroppers are kept under surveillance by an aerial UAV jammer (denoted by J). The aim of the UAV J is to cooperatively send jamming signals to the eavesdroppers to resist their wiretapping. Notice that if there is no friendly UAV J and only one eavesdropper is considered, the proposed scenario reduces to the goround-to-UAV transmission in [18].
Based on the three-dimensional Cartesian coordinate system, we denote w B = [ x B , y B ] T and w E i = [ x E i , y E i ] T as the horizontal coordinates of the BS and E i , respectively, which are assumed to be fixed and known beforehand to the UAVs. The assumption that w E i is known in the network is proper when E i is an active ground node but untrusted by the UAV [29]. Therefore, E i can be detected by the synthetic aperture radar or optical camera mounted on the UAV [18]. The initial and final locations of the UAVs are assumed to be pre-specified, which are denoted by q k , 0 = [ x k , 0 , y k , 0 ] T and q k , F = [ x k , F , y k , F ] T for k { U , J } , respectively. To make it more manageable, the period T is partitioned into N equal-length time slots, i.e., T = δ t N , where δ t is the length of one time slot. As such, the UAV trajectory in time slot n N can be represented approximately by q k [ n ] [ x k [ n ] , y k [ n ] ] T for k { U , J } , with a fixed altitude H. Let Ω = V max δ t be the maximum horizontal distance that the UAV can travel in a single time slot, where V max is the maximum speed of the UAV. Practically, the UAVs should satisfy the following mobility constraints,
| | q k [ n + 1 ] q k [ n ] | | 2 Ω 2 , n = 0 , , N 1 ,
q k [ 0 ] = q k , 0 , q k [ N ] = q k , F ,
| | q U [ n ] q J [ n ] | | 2 d min 2 , n = 0 , , N ,
where d min is the minimum tolerable distance between the two UAVs that ensures the avoidance of a collision.
We assume that the ground-to-UAV and UAV-to-UAV transmissions are mainly governed by LoS channels [18,20,23,24]. Thus, the corresponding channel power gains in time slot n follow the free-space path loss model. (For the purpose of exposition, it is reasonable to assume that the ground-to-UAV follows the free-space LoS channel model when the UAV is deployed in the rural area with sufficiently high altitude. In this case, the probability of Non-LoS state is negligible compared to the dominant LoS state [30]. However, the proposed design is readily extendable to more general channel models in urban areas with Non-LoS effects, e.g., [30].), which are, respectively, given as below,
h BU [ n ] = ρ 0 d BU 2 [ n ] = ρ 0 ( H H B ) 2 + | | q U [ n ] w B | | 2 ,
h JE i [ n ] = ρ 0 d JE i 2 [ n ] = ρ 0 H 2 + | | q J [ n ] w E i | | 2 ,
h JU [ n ] = ρ 0 d JU 2 [ n ] = ρ 0 | | q U [ n ] q J [ n ] | | 2 ,
where d BU [ n ] , d JE i [ n ] and d JU [ n ] are the distances from the BS to the UAV U, from the UAV J to the eavesdropper E i , and between the two UAVs in time slot n, respectively, ρ 0 is the channel power gain at the reference distance d 0 = 1 m and H B is the altitude of the BS. The ground-to-ground transmission is assumed to follow the Rayleigh fading channel. As such, the channel power gain is denoted by
h BE i = ρ 0 ζ i d B E i κ = ρ 0 ζ i | | w B w E i | | κ ,
where d B E i is the distance between the BS and the eavesdropper E i , ζ i is an exponentially distributed random variable with unit mean representing small-scale Rayleigh fading and κ 2 is the distance-dependent path loss exponent.
Denote by P [ n ] and Q [ n ] the BS’s transmit power and the UAV J’s jamming power in time slot n, respectively. In practice, they should satisfy the respective average power constraint P ¯ or Q ¯ , and peak power constraint P ^ or Q ^ , i.e.,
1 N n = 1 N P [ n ] P ¯ , 0 P [ n ] P ^ ,
1 N n = 1 N Q [ n ] Q ¯ , 0 Q [ n ] Q ^ ,
where P ¯ P ^ and Q ¯ Q ^ . Then, the achievable rate in bits/second/Hertz (bps/Hz) of the UAV U in time slot n is given by
R U [ n ] = log 2 1 + P [ n ] h BU [ n ] Q [ n ] h JU [ n ] + σ 2 ,
where Q [ n ] h JU [ n ] is the jamming interference from the UAV J, and σ 2 is the additive white Gaussian noise power at the receivers. Similarly, the achievable rate of the eavesdropper E i in time slot n can be expressed as
R E i [ n ] = E ζ i log 2 1 + P [ n ] h BE i Q [ n ] h JE i [ n ] + σ 2 log 2 1 + P [ n ] ρ 0 | | w B w E i | | κ Q [ n ] h JE i [ n ] + σ 2 R ^ E i [ n ] ,
where E ζ i [ · ] is the expectation operator with respect to (w.r.t.) ζ i . Note that R E i [ n ] is replaced by R ^ E i [ n ] based on Jensen’ inequality and the concavity of R E i [ n ] w.r.t. ζ i , and R ^ E i [ n ] is the largest rate that E i can achieve. Therefore, in accordance with the theoretical results in [31], the worst-case secrecy rate for each time slot can be lower bounded by
R wcs [ n ] = max ( R U [ n ] max i I R ^ E i [ n ] , 0 ) .
Note that by adaptively setting P [ n ] = 0 , the optimal solution to (12) is at least to be zero for any time slot n, without violating the power constraint (8). Therefore, the maximum operation can be dropped in the following optimization problems.

2.2. Problem Formulation

In this paper, we aim to maximize the average worst-case achievable secrecy rate from the BS to the UAV U over N time slots, by jointly optimizing the BS’s transmit power P { P [ n ] , n N } , the jamming power Q { Q [ n ] , n N } of the UAV J, and the UAV trajectory q k = { q k [ n ] , n N } for k { U , J } . Thus, this optimization problem can be formulated as
max P , Q , q U , q J 1 N n = 1 N ( R U [ n ] R E [ n ] )
s . t . ( 1 ) ( 3 ) , ( 8 ) ( 9 ) .
where we let R E [ n ] = max i I R ^ E i [ n ] , and thus R E [ n ] corresponds to the maximum achievable rate among multiple eavesdroppers in time slot n. Optimally solving problem (13) is difficult, in general, due to the following two main reasons: (1) the objective function is not concave w.r.t the corresponding optimization variables even with fixed variables of other blocks and (2) the constraint in (3) is non-convex w.r.t. the UAVs’ trajectory variables.

3. Joint Trajectory and Power Control Algorithm

In this section, an efficient algorithm is proposed to obtain the sub-optimal solution to problem (13). Specifically, we cope with problem (13) by solving four subproblems iteratively, i.e., the alternative optimization of the transmit power P , jamming power Q , UAV U’s trajectory q U , and UAV J’s trajectory q J , by fixing the other three optimization variables. Furthermore, the overall algorithm is presented, and its complexity and convergence are analyzed rigorously.

3.1. Transmit Power Optimization

For simplicity, let a n = γ 0 d BU 2 [ n ] ( 1 + Q [ n ] γ 0 / ρ 0 d JU 2 [ n ] ) , and b n = γ 0 | | w B w E | | κ 1 + Q [ n ] γ 0 / d JE 2 [ n ] , where γ 0 = ρ 0 / σ 2 is the reference signal-to-noise ratio (SNR), and w E is denoted as the horizontal location of the eavesdropper that achieves the largest rate and d JE 2 [ n ] is the distance from the UAV J to the eavesdropper w E . Thus, with given Q , q U , and q J , problem (13) can be simplified as
max P n = 1 N log 2 1 + a n P [ n ] log 2 1 + b n P [ n ]
s . t . ( 8 ) .
Based on the result in [18], the close-form solution to this problem is given by: P * [ n ] = min ( [ Λ n ] + , P ^ ) if a n > b n ; otherwise P * [ n ] = 0 , where Λ n = ( ( 1 / 2 b n 1 / 2 a n ) 2 + ( 1 / b n 1 / a n ) / ( λ ln 2 ) ) 1 2 1 / 2 a n 1 / 2 b n . The value of λ 0 is a constant that ensures constraint (8) is met, which can be obtained cost-effectively via the bisection algorithm [32]. By obtaining the optimal transmit power variables P , they can be seen as the given input for the jamming power optimization problem in the next subsections.

3.2. Jamming Power Optimization

Let c n = P [ n ] γ 0 d BU 2 [ n ] , d n = γ 0 d JU 2 [ n ] , e n = P [ n ] γ 0 | | w B w E | | κ and f n = γ 0 d JE 2 [ n ] . With given P , q U , and q J , we can reformulate problem (13) as
max Q n = 1 N log 2 1 + c n 1 + d n Q [ n ] log 2 1 + e n 1 + f n Q [ n ]
s . t . ( 9 ) .
Problem (15) is a non-convex problem because of the non-convex objective function, which is actually difficult to solve for general N. However, the first term in (15) is convex w.r.t. Q [ n ] , and thus it can be approximated to a convex function within each iteration by applying the SCA method. It is known that the first-order Taylor expansion can be used to obtain the global under-estimator for any convex function at any point [32]. Thus, denoted by Q l = { Q l [ n ] , n N } , the given local point in the l-th iteration, we have
log 2 1 + c n 1 + d n Q [ n ] A n l + B n l ( Q [ n ] Q l [ n ] )
where A n l = log 2 1 + c n 1 + d n Q l [ n ] and
B n l = c n d n ln 2 ( 1 + d n Q l [ n ] ) ( 1 + c n + d n Q l [ n ] ) .
With (16), problem (15) is lower bounded by the following problem for any given Q l ,
max Q n = 1 N B n l Q [ n ] log 2 1 + e n 1 + f n Q [ n ]
s . t . ( 9 ) .
Observe that this subproblem is concave w.r.t. Q [ n ] and thus can be solved efficiently by the interior-point method [32].After solving problem (17), the obtained jamming power Q serves as the given variables for the trajectory optimization problem of the UAVs.

3.3. Trajectory Optimization of the UAV U

Even with given P , Q , and q J , it is still hard to achieve the optimal solution to problem (13), due to the non-concavity of the objective function w.r.t. q U and the non-convexity of the constraint (3). To tackle this subproblem, we first introduce the slack variables α = { α [ n ] = ( H H B ) 2 + | | q U [ n ] w B | | 2 , n N } and β = { β [ n ] = | | q U [ n ] q J [ n ] | | 2 , n N } . After some simple transformations, solving problem (13) is equivalent to solve the following problem,
max q U , α , β n = 1 N log 2 1 + P [ n ] γ 0 α [ n ] + Q [ n ] γ 0 β [ n ] log 2 1 + Q [ n ] γ 0 β [ n ]
s . t . α [ n ] ( H H B ) 2 + | | q U [ n ] w B | | 2 ,
β [ n ] | | q U [ n ] q J [ n ] | | 2 , ( 1 ) ( 3 ) .
In fact, if α [ n ] ( β [ n ] ) is increased (decreased), the objective value of problem (13) will be decreased, and thus the constraints for α and β must satisfy the equalities. Problem (18) is still non-convex, because of the non-convex objective function in (18), and the constraints in (3) and (20). To tackle this difficulty, an important lemma is provided as below.
Lemma 1.
Given K 1 > 0 and K 2 > 0 , the function f ( x , y ) = log 2 1 + K 1 x + K 2 y is jointly convex w.r.t. x > 0 and y > 0 .
Proof. 
See Appendix A. □
Based on Lemma 1, it is easy to prove the convexity of the first term in problem (18). By using the first-order Taylor expansions of a convex function f ( x , y ) in a neighborhood of ( x , y ) = ( x 0 , y 0 ) , i.e., f ( x , y ) = f ( x 0 , y 0 ) + f x ( x 0 , y 0 ) ( x x 0 ) + f y ( x 0 , y 0 ) ( y y 0 ) , the first term in (18) at given local points denoted by α l = { α l [ n ] , n N } and β l = { β l [ n ] , n N } in the l-th iteration, can be given as follows,
log 2 1 + P [ n ] γ 0 α [ n ] + Q [ n ] γ 0 β [ n ] log 2 C n l D n l C n l ln 2
where C n l = 1 + P [ n ] γ 0 α l [ n ] + Q [ n ] γ 0 β l [ n ] and D n l = P [ n ] γ 0 ( α l [ n ] ) 2 ( α [ n ] α l [ n ] ) + Q [ n ] γ 0 ( β l [ n ] ) 2 ( β [ n ] β l [ n ] ) . Similarly, by using the first-order Taylor expansion at the given local point denoted by q U l = { q U l [ n ] , n N } in the l-th iteration, the convex function | | q U [ n ] q J [ n ] | | 2 , w.r.t. q U [ n ] in problem (3) and in problem (20) can be replaced by their convex lower bounds, i.e.,
| | q U [ n ] q J [ n ] | | 2 | | q U l [ n ] q J [ n ] | | 2 + 2 ( q U l [ n ] q J [ n ] ) T ( q U [ n ] q U l [ n ] ) .
As a result, by applying SCA technique in each iteration, we approximate the original convex functions to more manageable functions at given local points. Therefore, with (21)–(22), we have the following optimization problem
max q U , α , β n = 1 N D n l C n l ln 2 log 2 1 + Q [ n ] γ 0 β [ n ]
s . t . β [ n ] | | q U l [ n ] q J [ n ] | | 2 + 2 ( q U l [ n ] q J [ n ] ) T ( q U [ n ] q U l [ n ] ) ,
d min 2 | | q U l [ n ] q J [ n ] | | 2 + 2 ( q U l [ n ] q J [ n ] ) T ( q U [ n ] q U l [ n ] ) , ( 1 ) , ( 2 ) , ( 19 ) .
It is observed that problem (23) is now convex with all convex constraints. As such, the interior-point method can be used efficiently to solve this problem. Note that the lower bounds obtained by the Taylor expansions suggest that the optimal objective value by solving problem (23) is a lower bound of that of problem (18). In the next subsection, the solved q U is input to the trajectory optimization problem of the UAV J as the given variable.

3.4. Trajectory Optimization of the UAV J

With given P , Q , and q U , we let δ = { δ [ n ] = H 2 + | | q J [ n ] w E | | 2 , n N } , to tackle the non-concavity of the objective function w.r.t. q J . Therefore, problem (13) can be rewritten as
max q J , β , δ n = 1 N [ log 2 β [ n ] + c n β [ n ] + Q [ n ] γ 0 log 2 β [ n ] + Q [ n ] γ 0 log 2 1 + e n δ [ n ] δ [ n ] + Q [ n ] γ 0 ]
s . t . δ [ n ] H 2 + | | q [ n ] w E | | 2 , ( 1 ) ( 3 ) , ( 20 ) .
Also, the constraint for δ holds with equalities, otherwise the objective value of problem (13) will be decreased by increasing δ [ n ] . Similarly, by using the first-order Taylor expansion at given local points denoted by δ l = { δ l [ n ] , n N } , β l = { β l [ n ] , n N } and q J l = { q J l [ n ] , n N } in the l-th iteration, the second and third terms in problem (26), and | | q U [ n ] q J [ n ] | | 2 in (3) and in (20) can be substituted by their respective concave upper and convex lower bounds, i.e.,
log 2 β [ n ] + Q [ n ] γ 0 log 2 β l [ n ] + Q [ n ] γ 0 + β [ n ] β l [ n ] ln 2 ( β l [ n ] + Q [ n ] γ 0 ) ,
log 2 1 + e n δ [ n ] δ [ n ] + Q [ n ] γ 0 E n l + F n l ( δ [ n ] δ l [ n ] ) ,
where E n l = log 2 1 + e n δ l [ n ] δ l [ n ] + Q [ n ] γ 0 ,
F n l = e n γ 0 Q [ n ] ln 2 ( γ 0 Q [ n ] + ( e n + 1 ) δ l [ n ] ) ( γ 0 Q [ n ] + δ l [ n ] ) ,
and
| | q U [ n ] q J [ n ] | | 2 | | q U [ n ] q J l [ n ] | | 2 2 ( q U [ n ] q J l [ n ] ) T ( q J [ n ] q J l [ n ] ) .
With problems (28)–(30), we approximate problem, (26) as the following optimization problem
max q J , β , δ n = 1 N log 2 β [ n ] + c n β [ n ] + Q [ n ] γ 0 F n l δ [ n ] β [ n ] ln 2 ( β l [ n ] + Q [ n ] γ 0 )
s . t . d min 2 | | q U [ n ] q J l [ n ] | | 2 2 ( q U [ n ] q J l [ n ] ) T ( q J [ n ] q J l [ n ] ) ,
β [ n ] | q U [ n ] q J l [ n ] | | 2 2 ( q U [ n ] q J l [ n ] ) T ( q J [ n ] q J l [ n ] ) , ( 1 ) , ( 2 ) , ( 27 ) .
Problem (31) is now a convex optimization problem that can be cost-effectively solved by the interior-point method. Furthermore, the Taylor expansions in problems (28)–(30) indicate that the objective value of problem (26) is at least the same as that by solving problem (31). Note that all the obtained variables P , Q , q U , and q J are utilized as the given variables for the next iteration.

3.5. Overall Algorithm

In summary, the overall algorithm for obtaining the locally optimal solution to problem (13) is computed by the joint optimization of both the BS’s transmit power P and the UAV J’s jamming power Q as well as the two UAVs’ trajectories q U , and q J variables, via alternatively solving subproblems (14), (17), (23) and (31) in an iterative way, respectively. The detailed procedure for solving problem (13) is summarized in Algorithm 1.
In the following, we analyze the computation complexity of Algorithm 1. In each iteration, the BS’s transmit power, UAV J’s jamming power, and the trajectories of UAVs U and J are optimized in sequence, based on the interior-point method by using existing solvers, such as CVX [33]. Therefore, the complexity for solving the four subproblems can be expressed by O ( log N ) , O ( N 3.5 log ( 1 / ϵ ) ) , O ( ( 3 N ) 3.5 log ( 1 / ϵ ) ) , and O ( ( 3 N ) 3.5 log ( 1 / ϵ ) ) , respectively, for the given solution precision of ϵ > 0 [34]. In addition, as the complexity for updating all variables in BCD iterations is in the order of log ( 1 / ϵ ) , the total computation complexity of the proposed algorithm is O ( N 3.5 log 2 ( 1 / ϵ ) ) . Due to the polynomial time complexity, Algorithm 1 is applicable to the aerial cooperative jamming for cellular-enabled UAV networks.
Algorithm 1 Proposed algorithm for solving problem (13)
1: Initial P , Q , q U , q J , α , β and δ . Let l = 0 .
2: repeat
3:  Solve problem (14) with given Q l , q U l , and q J l ,
  and denote by P l + 1 the optimal solution.
4:  Solve problem (17) with given P l , q U l , and q J l ,
  and denote by Q l + 1 the optimal solution.
5:  Solve problem (23) with given P l , Q l , q J l , α l and β l ,
  and denote by q U l + 1 the optimal solution.
6:  Solve problem (31) with given P l , Q l , q U l , β l and δ l ,
  and denote by q J l + 1 the optimal solution.
7:  Update l = l + 1 .
8: until Converge to a pre-specified precision ϵ > 0 .
Next, the convergence of Algorithm 1 is discussed as follows. Let ψ ( P l , Q l , q U l , q J l ) denote the value of the objective function in problem (13) in the l-th iteration. Then, we have
ψ ( P l , Q l , q U l , q J l ) ψ P ( P l + 1 , Q l , q U l , q J l ) ,
where ψ P ( P l + 1 , Q l , q U l , q J l ) is defined as the obtained objective value of problem (14) and P l + 1 is the optimal solution to problem (14). For the optimization of the jamming power Q , the following equations hold,
ψ ( P l + 1 , Q l , q U l , q J l ) = ( j 1 ) ψ Q lb ( P l + 1 , Q l , q U l , q J l ) ( j 2 ) ψ Q lb ( P l + 1 , Q l + 1 , q U l , q J l ) ( j 3 ) ψ ( P l + 1 , Q l + 1 , q U l , q J l ) ,
where ψ Q lb is denoted as the objective value of problem (17), ( j 1 ) holds since the first-order Taylor expansion in (16) is tight at the local point Q l in problem (17), ( j 2 ) satisfies due to the optimal solution Q l + 1 to problem (17), and ( j 3 ) is because the computed objective value of problem (17) is lower bounded by that of problem (15). For the two UAVs’ trajectories optimization, the similar derivation procedure as in (35) can be used, which are given as below,
ψ ( P l + 1 , Q l + 1 , q U l , q J l ) = ψ q U lb ( P l + 1 , Q l + 1 , q U l , q J l ) ψ q U lb ( P l + 1 , Q l + 1 , q U l + 1 , q J l ) ψ ( P l + 1 , Q l + 1 , q U l + 1 , q J l ) ,
ψ ( P l + 1 , Q l + 1 , q U l + 1 , q J l ) = ψ q J lb ( P l + 1 , Q l + 1 , q U l + 1 , q J l ) ψ q J lb ( P l + 1 , Q l + 1 , q U l + 1 , q J l + 1 ) ψ ( P l + 1 , Q l + 1 , q U l + 1 , q J l + 1 ) ,
With (34)–(37), we finally obtain that
ψ ( P l , Q l , q U l , q J l ) ψ ( P l + 1 , Q l + 1 , q U l + 1 , q J l + 1 ) .
As a result, Algorithm 1 ensures that the obtained objective value of problem (13) is non-decreasing over the iterations, and thus it guarantees its convergence to the locally optimal solution to problem (13).

4. Numerical Results

In this section, we verify our joint trajectories and powers optimization (denoted as 2T&P) algorithm through simulations. Three benchmark schemes are taken into account as a comparison:
  • UAVs’ trajectories optimization without power control (denoted as 2T/NP);
  • heuristic UAVs’ trajectories with power control (2HT/P);
  • joint optimization of the UAV U’s trajectory and BS’s power control without aerial cooperative jamming from the UAV J (denoted as 1T&P), which is identical with the algorithm proposed in [18].
Specifically, the 2T/NP scheme sets the powers of the BS and the UAV U as P [ n ] = P ¯ and Q [ n ] = Q ¯ , n , respectively, and the trajectories of the two UAVs are obtained by solving problems (23) and (31) iteratively until convergence. In the 2HT/P scheme, the UAV U flies directly to the top of the BS at its maximum speed, then stays hovering as long as possible, and finally travels directly to its destination at its maximum speed by the end of T. Different from UAV U, UAV J keeps hovering right above the eavesdropper with the largest achievable rate. Given heuristic trajectories in the 2HT/P, the powers P [ n ] and Q [ n ] can be obtained by solving problems (14) and (17), respectively. The initial UAV trajectory for the 2T&P and 2T/NP schemes are constructed by the heuristic UAV trajectories as in 2HT/P. The simulation parameters are specified in Table 1.
We first verify the convergence behaviour of the proposed Algorithm 1 versus the iteration numbers for different T in Figure 2. It is illustrated that the average secrecy rate increases quickly and converges within five iterations, and its performance increases significantly with T. This confirms that a locally optimal solution to problem (13) can be converged by using the proposed algorithm.
Figure 3 illustrates the optimized trajectories of the two UAVs by different schemes when T is sufficiently large, e.g., T = 300 s. It is observed that the hovering locations of all algorithms for the UAV U are directly above the BS. This occurs because the locations of the eavesdroppers are not related to the UAV U’s trajectory due to the ground-to-air transmission, and thus the UAV U can obtain its maximum achievable rate hovering at the location on top of the BS. In addition, the the trajectories of the UAV U in 2T&P and 2T/NP show the curved paths in order to escape from the unintended interference caused by the UAV J. However, the trajectories of the UAV J present significant different. In particular, for our 2T&P scheme in Figure 3a, the UAV first flies along an arc-like path and reaches a certain point close to the eavesdropper E 1 to avoid a collision with the UAV U; then, it keeps static at this hovering location for a permission period, and finally reaches its destination by the end of T, also in an arc-like path to prevent it causing much interference for the UAV U. Notice that the hovering location of the UAV J is closer to E 1 compared to E 2 , as the channel quality of BS-to- E 1 link is much better than that of BS-to- E 2 link. The BS-to- E 2 link can also be degraded if the UAV J can guarantee that the secrecy of the worst-case, i.e., BS-to- E 1 link transmission, by taking advantage of the dominant air-to-ground links. Moreover, at their hovering locations, the UAVs can achieve the better secrecy rate by effectively balancing between enhancing the communication of the ground-to-air link and degrading the quality of the BS-to- E i channel. In contrast with the 2T&P scheme, we can observe that on its way to the final location, the UAV J flies in a big arc path to keep away from the UAV U in the 2T/NP scheme as shown in Figure 3c. This is because the BS’ s transmit power and UAV J’s jamming power in 2T/NP are fixed, and thus the UAV J has to fly as far as possible to avoid severe interference with the UAV U over the whole duration, T.
Note that there is a tradeoff between improving the average achievable secrecy rate of the UAV U and avoiding the interference induced by the UAV J. For 2HT/P in Figure 3b, with the pre-specified UAV trajectories, the BS and the UAV jammer can adjust their power allocations to enhance the secrecy performance. Specifically, the BS gradually increases its transmit power before the UAV U flies to its hovering location, while the UAV J properly decreases its jamming power when it reaches above E 1 to suppress the interference to the UAV U. In contrast, a secure communication-aware UAV trajectory design provides additional flexibility to avoid interference between UAVs in our 2T&P scheme. Thus, the UAV J adaptively adjusts its jamming power and trajectory according to the BS’s transmit power and the location of the UAV U to further achieve the better secrecy rate.
Figure 4 illustrates the average secrecy rate versus T. It is expected that the average secrecy rates obtained by all schemes raise with T, and the proposed 2T&P scheme significantly outperforms other benchmark schemes owing to its joint optimization. Moreover, the proposed 2T&P scheme provides the significant gain as compared to the scheme in [18], i.e., 1T&P. This indicates that the advantage brought by the aerial cooperative jamming is more effective and important on notably improving the average secrecy rate. However, the 2T/NP presents the worst performance, which demonstrates that the power control also plays a key role in avoiding jamming from other UAVs, which is necessary in our cellular-enabled UAV communication networks with aerial cooperative jamming; otherwise, the secrecy rate can be significantly degraded as shown in Figure 4. Note that we expect that the proposed 2T&P algorithm can still achieve the best secrecy performance via the joint design, even if the number of the eavesdropper increases. This is because the joint design guarantees that the eavesdropper with the best channel condition can be effectively jammed by the UAV J; other eavesdroppers cannot wiretap confidential messages from the BS. The obtained results validate the advantages of introducing aerial cooperative jamming, and the joint optimization of UAV trajectories and power allocations.

5. Conclusions

Integrating UAVs into the forthcoming 5G cellular networks faces new security challenges. Thus, a new type of cooperative aerial jamming scheme for the cellular-enabled UAV secure communication networks has been investigated in this paper. In particular, the UAV receiver and the UAV jammer cooperate closely with each other to maximize the worst-case average secrecy rate by jointly optimizing their trajectories and the BS/UAV transmit/jamming power. An efficient iterative solution has been proposed to approximately tackle the secrecy rate maximization problem over a given flight period, by means of the BCD and SCA methods. The proposed algorithm is guaranteed to converge to a locally optimal solution with suitable computational complexity. We have demonstrated, by numerical results, that the friendly UAV jammer provides flexible mobility for interference with the ground eavesdroppers, as well as effective power control of preventing it from jamming the UAV receiver, and thereby improves the system secrecy performance. Furthermore, the proposed scheme significantly outperforms the benchmark schemes with simple heuristic trajectories and pre-configured powers. The current scenario can also be extended to the general case with multiple legitimate UAVs, where optimal communication scheduling between the BS and each UAV should be considered. In this case, the design for UAV trajectories needs to avoid collision between UAVs more effectively, and reconciles a tradeoff between maximizing the minimum secrecy rate among multiple UAVs and suppressing the interference from the UAV jammer, which is an interesting problem to be resolved in the future.

Author Contributions

Conceptualization, methodology, and software, H.S. and B.D.; validation and investigation, Z.W. and X.L.; and writing—original draft preparation, B.D.; and writing—review and editing, H.S. and C.G.

Funding

This work was partially supported by the National Natural Science Foundation of China (Nos. 61701064, 71874027), partially by the China Scholarship Council (No. 201808510008), partially by the Sichuan Science and Technology Program (No. 2018GZ0454), partially by the Chengdu Science and Technology Program (No. 2017-RK00-00363-ZF), and partially by the Chongqing Natural Science Foundation (cstc2019jcyj-msxmX0264).

Acknowledgments

The authors the appreciate suggestions of the editors and reviewers, which are of great value to the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Appendix

We prove Lemma 1 via the definition of convex functions. First, since f ( x , y ) = log 2 1 + K 1 x + K 2 y where x > 0 , y > 0 , K 1 > 0 and K 2 > 0 , the first-order derivatives of f ( x , y ) w.r.t. x and y are given by
f x ( x , y ) = K 1 x 2 G ln 2 , f y ( x , y ) = K 2 y 2 G ln 2 ,
where we let G = 1 + K 1 x + K 2 y for brevity. Then, the Hessian of f ( x , y ) is
2 f ( x , y ) = 2 K 1 x y + 2 K 1 K 2 x + K 1 2 y x 4 y G 2 ln 2 K 1 K 2 x 2 y 2 G 2 ln 2 K 1 K 2 x 2 y 2 G 2 ln 2 2 K 2 x y + 2 K 1 K 2 y + K 2 2 x x y 4 G 2 ln 2
For any t = [ t 1 , t 2 ] T , we have
t T 2 f ( x , y ) t = F 1 ( x , y ) + F 2 ( x , y ) + F 3 ( x , y ) x 4 y 4 G 2 ln 2 0
where F 1 ( x , y ) = 2 K 1 t 1 2 x y 3 ( y + K 2 ) , F 2 ( x , y ) = 2 K 2 t 2 2 x 3 y ( x + K 1 ) and F 3 ( x , y ) = ( t 1 K 1 y 2 t 2 K 2 x 2 ) 2 for x > 0 , y > 0 , K 1 > 0 and K 2 > 0 , which finally leads to the convexity of the function f ( x , y ) .

References

  1. Moon, J.; Lee, H.; Song, C.; Lee, S.; Lee, I. Proactive eavesdropping with full-duplex relay and cooperative jamming. IEEE Trans. Wirel. Commun. 2018, 17, 6707–6719. [Google Scholar] [CrossRef]
  2. Lu, H.; Zhang, H.; Dai, H.; Wu, W.; Wang, B. Proactive eavesdropping in UAV-aided suspicious communication systems. IEEE Trans. Veh. Technol. 2019, 68, 1993–1997. [Google Scholar] [CrossRef]
  3. Li, Y.; Zhang, Y.; Cai, L. Optimal location of supplementary node in UAV surveillance system. J. Netw. Comput. Appl. 2019, 140, 23–39. [Google Scholar] [CrossRef]
  4. Zhao, N.; Lu, W.; Sheng, M.; Chen, Y.; Tang, J.; Yu, F.R.; Wong, K. UAV-assisted emergency networks in disasters. IEEE Wirel. Commun. 2019, 26, 45–51. [Google Scholar] [CrossRef]
  5. Tuna, G.; Nefzi, B.; Conte, G. Unmanned aerial vehicle-aided communications system for disaster recovery. J. Netw. Comput. Appl. 2014, 41, 27–36. [Google Scholar] [CrossRef]
  6. Peng, K.; Du, J.; Lu, F.; Sun, Q.; Dong, Y.; Zhou, P.; Hu, M. A hybrid genetic algorithm on routing and scheduling for vehicle-assisted multi-drone parcel delivery. IEEE Access 2019, 7, 49191–49200. [Google Scholar] [CrossRef]
  7. You, C.; Zhang, R. 3D trajectory optimization in Rician fading for UAV-enabled data harvesting. IEEE Wirel. Commun. 2019, 18, 3192–3207. [Google Scholar] [CrossRef]
  8. Zeng, Y.; Zhang, R.; Lim, T.J. Throughput maximization for UAV-enabled mobile relaying systems. IEEE Wirel. Commun. 2016, 64, 4983–4996. [Google Scholar] [CrossRef]
  9. Zeng, Y.; Wu, Q.; Zhang, R. Accessing from the Sky: A Tutorial on UAV Communications for 5G and Beyond. 2019. Available online: https://arxiv.org/abs/1903.05289 (accessed on 13 March 2019).
  10. Zeng, Y.; Zhang, R.; Lim, T.J. Wireless communications with unmanned aerial vehicles: Opportunities and challenges. IEEE Commun. Mag. 2016, 54, 36–42. [Google Scholar] [CrossRef]
  11. Zeng, Y.; Lyu, J.; Zhang, R. Cellular-connected UAV: Potential, challenges, and promising technologies. IEEE Wirel. Commun. 2019, 26, 120–127. [Google Scholar] [CrossRef]
  12. Wu, Q.; Zhang, R. Common throughput maximization in UAV-enabled OFDMA systems with delay consideration. IEEE Trans. Commun. 2018, 66, 6614–6627. [Google Scholar] [CrossRef]
  13. Wu, Q.; Zeng, Y.; Zhang, R. Joint trajectory and communication design for multi-UAV enabled wireless networks. IEEE Trans. Wirel. Commun. 2018, 17, 2109–2121. [Google Scholar] [CrossRef]
  14. Lyu, J.; Zeng, Y.; Zhang, R.; Lim, T.J. Placement optimization of UAV-mounted mobile base stations. IEEE Commun. Lett. 2017, 21, 604–607. [Google Scholar] [CrossRef]
  15. Zhang, S.; Zeng, Y.; Zhang, R. Cellular-enabled UAV communication: A connectivity-constrained trajectory optimization perspective. IEEE Trans. Commun. 2019, 67, 2580–2604. [Google Scholar] [CrossRef]
  16. Lin, X.; Yajnanarayana, V.; Muruganathan, S.D.; Gao, S.; Asplund, H.; Maattanen, H.; Bergstrom, M.; Euler, S.; Wang, Y.E. The sky is not the Limit: LTE for unmanned aerial vehicles. IEEE Commun. Mag. 2018, 56, 204–210. [Google Scholar] [CrossRef]
  17. Der Bergh, B.V.; Chiumento, A.; Pollin, S. LTE in the sky: Trading off propagation benefits with interference costs for aerial nodes. IEEE Commun. Mag. 2016, 54, 44–50. [Google Scholar] [CrossRef]
  18. Zhang, G.; Wu, Q.; Cui, M.; Zhang, R. Securing UAV communications via joint trajectory and power control. IEEE Trans. Wirel. Commun. 2019, 18, 1376–1389. [Google Scholar] [CrossRef]
  19. Cui, M.; Zhang, G.; Wu, Q.; Ng, D.W.K. Robust trajectory and transmit power design for secure UAV communications. IEEE Trans. Veh. Technol. 2018, 67, 9042–9046. [Google Scholar] [CrossRef]
  20. Li, A.; Wu, Q.; Zhang, R. UAV-enabled cooperative jamming for improving secrecy of ground wiretap channel. IEEE Wirel. Commun. Lett. 2019, 8, 181–184. [Google Scholar] [CrossRef]
  21. Duo, B.; Wu, Q.; Zhang, R. Energy Efficiency Maximization for Full-Duplex UAV Secrecy Communication. 2019. Available online: https://arxiv.org/abs/1906.07346 (accessed on 18 June 2019).
  22. Wang, Q.; Chen, Z.; Li, H.; Li, S. Joint power and trajectory design for physical-layer secrecy in the UAV-aided mobile relaying system. IEEE Access 2018, 6, 62849–62855. [Google Scholar] [CrossRef]
  23. Zhong, C.; Yao, J.; Xu, J. Secure UAV communication with cooperative jamming and trajectory control. IEEE Commun. Lett. 2019, 32, 286–289. [Google Scholar] [CrossRef]
  24. Cai, Y.; Cui, F.; Shi, Q.; Zhao, M.; Li, G.Y. Dual-UAV enabled secure communications: Joint trajectory design and user scheduling. IEEE J. Sel. Areas Commun. 2018, 36, 1972–1985. [Google Scholar] [CrossRef]
  25. Beck, A.; Tetruashvili, L. On the convergence of block coordinate descent type methods. SIAM J. Optim. 2013, 23, 2037–2060. [Google Scholar] [CrossRef]
  26. Hong, M.; Razaviyayn, M.; Luo, Z.; Pang, J. A unified algorithmic framework for block-structured optimization involving big data: With applications in machine learning and signal processing. IEEE Sign. Process. Mag. 2016, 33, 57–77. [Google Scholar] [CrossRef]
  27. Razaviyayn, M.; Hong, M.; Luo, Z. A unified convergence analysis of block successive minimization methods for nonsmooth optimization. SIAM J. Optim. 2013, 23, 1126–1153. [Google Scholar] [CrossRef]
  28. Zeng, Y.; Zhang, R. Energy-efficient UAV communication with trajectory optimization. IEEE Trans. Wirel. Commun. 2017, 16, 3747–3760. [Google Scholar] [CrossRef]
  29. Mukherjee, A.; Fakoorian, S.A.A.; Huang, J.; Swindlehurst, A.L. Principles of physical layer security in multiuser wireless networks: A survey. IEEE Commun. Surv. Tuts 2014, 16, 1550–1573. [Google Scholar] [CrossRef]
  30. You, C.; Zhang, R. Hybrid Offline-Online Design for UAV-Enabled Data Harvesting in Probabilistic LoS Channel. 2019. Available online: https://arxiv.org/abs/1907.06181v1 (accessed on 14 July 2019).
  31. Gopala, P.K.; Lai, L.; El Gamal, H. On the secrecy capacity of fading channels. IEEE Trans. Inf. Theory 2008, 54, 4687–4698. [Google Scholar] [CrossRef]
  32. Boyd, S.; Boyd, S.; Vandenberghe, L.; Press, C.U. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  33. Grant, M.; Boyd, S. CVX: Matlab Software for Disciplined Convex Programming, Version 2.1. 2018. Available online: http://cvxr.com/cvx (accessed on 14 October 2019).
  34. Ben-Tal, A.; Nemirovskiaei, A.S. Lectures on Modern Convex Optimization: Analysis, Algorithms, and Engineering Applications; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2001. [Google Scholar]
Figure 1. Cellular-enabled UAV secure communication network with aerial cooperative jamming.
Figure 1. Cellular-enabled UAV secure communication network with aerial cooperative jamming.
Sensors 19 04440 g001
Figure 2. The convergence performance of the proposed Algorithm 1.
Figure 2. The convergence performance of the proposed Algorithm 1.
Sensors 19 04440 g002
Figure 3. UAVs’ trajectories by different schemes for T = 300 s. All trajectories are sampled every 5 s. The horizontal locations of the BS, eavesdroppers, UAVs’ initial and final locations are marked with ▲, ⚫, ☐ and ◇, respectively.
Figure 3. UAVs’ trajectories by different schemes for T = 300 s. All trajectories are sampled every 5 s. The horizontal locations of the BS, eavesdroppers, UAVs’ initial and final locations are marked with ▲, ⚫, ☐ and ◇, respectively.
Sensors 19 04440 g003
Figure 4. Average secrecy rate versus T with different trajectory and power control designs.
Figure 4. Average secrecy rate versus T with different trajectory and power control designs.
Sensors 19 04440 g004
Table 1. Simulation parameters.
Table 1. Simulation parameters.
NotationPhysical MeaningSimulation Value
q U , 0 Initial horizontal location of the UAV U [ 800 , 200 ] T m
q U , F Final horizontal location of the UAV U [ 800 , 200 ] T m
q J , 0 Initial horizontal location of the the UAV J [ 800 , 200 ] T m
q J , F Final horizontal location of the UAV J [ 800 , 200 ] T m
w B Horizontal location of the BS [ 0 , 0 ] T m
w E 1 Location of the first eavesdropper [ 150 , 0 ] T m,
w E 2 Location of the second eavesdropper [ 200 , 0 ] T m
HAltitude of UAVs100 m
H B Altitude of BS10 m
V max Maximum speed of UAVs40 m/s
d min Minimum safe distance between UAVs20 m
δ t Time slot length1 s
ρ 0 Channel power gain at the reference distance−60 dB
σ 2 Noise power levels−110 dBm
P ¯ and P ^ Average and peak power of UAV U20 dBm and 26 dBm
Q ¯ and Q ^ Average and peak power of UAV J10 dBm and 16 dBm
ϵ Accuracy threshold 10 4

Share and Cite

MDPI and ACS Style

Sun, H.; Duo, B.; Wang, Z.; Lin, X.; Gao, C. Aerial Cooperative Jamming for Cellular-Enabled UAV Secure Communication Network: Joint Trajectory and Power Control Design. Sensors 2019, 19, 4440. https://doi.org/10.3390/s19204440

AMA Style

Sun H, Duo B, Wang Z, Lin X, Gao C. Aerial Cooperative Jamming for Cellular-Enabled UAV Secure Communication Network: Joint Trajectory and Power Control Design. Sensors. 2019; 19(20):4440. https://doi.org/10.3390/s19204440

Chicago/Turabian Style

Sun, Hanming, Bin Duo, Zhengqiang Wang, Xiaochen Lin, and Changchun Gao. 2019. "Aerial Cooperative Jamming for Cellular-Enabled UAV Secure Communication Network: Joint Trajectory and Power Control Design" Sensors 19, no. 20: 4440. https://doi.org/10.3390/s19204440

APA Style

Sun, H., Duo, B., Wang, Z., Lin, X., & Gao, C. (2019). Aerial Cooperative Jamming for Cellular-Enabled UAV Secure Communication Network: Joint Trajectory and Power Control Design. Sensors, 19(20), 4440. https://doi.org/10.3390/s19204440

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