Next Article in Journal
Predictor-Based Motion Tracking Control for Cloud Robotic Systems with Delayed Measurements
Previous Article in Journal
A Novel Arc Fault Detection Method Integrated Random Forest, Improved Multi-scale Permutation Entropy and Wavelet Packet Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stackelberg Game-Theoretic Low Probability of Intercept Performance Optimization for Multistatic Radar System

1
Key Laboratory of Radar Imaging and Microwave Photonics, Nanjing University of Aeronautics and Astronautics, Ministry of Education, Nanjing 210016, China
2
Science and Technology on Electro-Optic Control Laboratory, Luoyang 471009, China
3
Key Laboratory of Dynamic Cognitive System of Electromagnetic Spectrum Space Nanjing University of Aeronautics and Astronautics, Ministry of Industry and Information Technology, Nanjing 210016, China
4
School of Engineering and Computing Sciences, Durham University, Durham DH1 3DE, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Electronics 2019, 8(4), 397; https://doi.org/10.3390/electronics8040397
Submission received: 12 March 2019 / Revised: 26 March 2019 / Accepted: 29 March 2019 / Published: 2 April 2019
(This article belongs to the Section Circuit and Signal Processing)

Abstract

:
In this paper, the problem of Stackelberg game-theoretic low probability of intercept (LPI) performance optimization in multistatic radar system is investigated. The goal of the proposed LPI optimization strategy is to minimize the transmitted power of each radar while satisfying a predetermined signal-to-interference-plus-noise ratio (SINR) requirement for target detection. Firstly, a single-leader multi-follower Stackelberg game is adopted to formulate the LPI optimization problem of multistatic radar system. In the considered game model, the hostile intercept receiver plays a role of leader, who decides the prices of power resource first through the maximization of its own utility function. The multiple radars are followers to compete with each other in a non-cooperative game according to the imposed prices from the intercept receiver subsequently. Then, the Nash equilibrium (NE) for the considered game model is derived, and the existence and uniqueness of the NE are analytically proved. Furthermore, a pricing-based distributed iterative power control algorithm is proposed. Finally, some simulation examples are provided to demonstrate that the proposed scheme has remarkable potential to enhance the LPI performance of the multistatic radar system.

1. Introduction

With the recent technology development of advanced passive intercept systems, it is highly important to find a solution to the problem of low probability of intercept (LPI) performance enhancement for different radar systems caused by this trend [1]. In theory, low radiated power [2,3,4], short dwell time [5,6], large revisit interval [7], adaptive beamforming [8,9,10], and waveform optimization [11] will result in improved LPI performance. In recent years, researchers have developed a great number of techniques to satisfy the LPI performance requirement. For example, in [2], the revisit interval, transmit power, and waveform parameters are jointly optimized to improve the LPI performance for radar networks. In [4], an LPI-bsaed joint transmitter selection and resource management scheme is proposed for single target tracking in radar network, where the LPI performance criterion of radar network is minimized by optimizing the revisit interval, dwell time, transmitter selection, and transmit power while maintaining a specified target tracking accuracy. In [6], Wang X.L. et al. propose a joint revisit and dwell time management approach for target tracking in phased array radar system, in which the time resource consumption is minimized for a given target tracking performance. Lawrence D.E. in [8] introduces a novel transmit array beamforming algorithm, which offers a better LPI performance for surveillance radars exploiting phased array antennas. Reference [10] presents a cognitive LPI-based transmit beamforming algorithm by utilizing frequency diverse array (FDA) and multiple-input multiple-output (MIMO) hybrid array antenna, which minimizes the beam power at the target location while maximizing the power at the radar receiver without degrading the target detection performance. The work in [11] addresses the problem of power minimization-based robust orthogonal frequency division multiplexing (OFDM) radar waveform design for radar and communication systems in spectral coexistence. It is shown that the LPI performance of the radar system can be efficiently strengthened by employing the communication waveforms scattered off the target at radar receiver.
Game theory provides a natural and efficient tool in modeling the interactions among independent players [12,13,14]. A lot of work has been developed for radar systems and made significant progress [15,16,17,18,19,20]. In [16], the authors model the interaction between a smart target and a smart MIMO radar as a two-person zero-sum game. The non-cooperative game-based code design approach in radar networks is proposed in [17], whose purpose is to maximize the signal-to-interference-plus-noise ratio (SINR) of each radar. The work in [18] studies the game-theoretic power allocation for a distributed MIMO radar network. Shi C.G. et al. in [19] propose a novel non-cooperative game-theoretic power allocation strategy for multistatic radar in a spectrum sharing environment. As an extension, a cooperative game-theoretic framework is proposed for power control in multistatic radar underlaying a communication system [20]. However, the above game-theoretic resource management protocols ignore the presence of hostile intercept receiver. The Stackelberg game-theoretic model has been exploited in several research [21,22,23,24] to analyze the hierarchical competition with different optimization purposes. In [25], the authors formulate a MIMO radar and target Stackelberg game model in the presence of clutter, and various optimization criteria at Stackelberg equilibrium (SE) of target dominant and radar dominant are obtained, respectively. However, in view of the above studies, there are still no published references that investigate this hierarchical interactions between the hostile intercept receiver and multistatic radar system. Therefore, the problem of Stackelberg game-theoretic LPI performance optimization for multistatic radar system should be addressed.
Specifically, the major contributions of this study are as follows:
(1)
The problem of Stackelberg game-theoretic LPI performance optimization strategy for multistatic radar system is investigated. Mathematically, the LPI optimization strategy can be formulated as a problem of minimizing the radiated power of each radar for a specified target detection performance. In earlier literature, although both non-cooperative game [19] and cooperative game [20] have been utilized to control the transmit power of multistatic radar system, as aforementioned, the hostile intercept receiver is not considered in such a scenario. Therefore, we take the hierarchical interactions between intercept receiver and multiple radars into consideration and formulate the LPI performance optimization between them as a single-leader multiple-follower Stackelberg game. In the underlying game model, the hostile intercept receiver plays a role of leader, who decides the prices of unit power resource first through the maximization of its own utility. The multiple radars are followers to compete with each other in a non-cooperative game according to the imposed prices from the interceptor subsequently.
(2)
We incorporate the total received power at intercept receiver, the unit power prices, the specified SINR requirement, and the transmit power of each radar to define the novel utility functions for the single leader and multiple followers. Then, we analyze the followers’ non-cooperative game model with the released prices from the leader, and the Nash equilibrium (NE) solution for the considered game model is derived. Additionally, the existence and uniqueness of the NE solution are strictly proved.
(3)
A pricing-based distributed iterative power control method is presented to solve the resulting optimization problem, which guarantees the convergence to the Stackelberg equilibrium (SE) points.
(4)
Some numerical examples are provided to confirm the convergence of the approach to the unique SE solution and verify the effectiveness of our proposed strategy in terms of LPI performance enhancement.
The rest of this paper is structured as follows: Section 2 provides the system model and assumptions. In Section 3, the Stackelberg game model for the LPI performance optimization is formulated. Section 4 provides the numerical results and analyses to demonstrate the proposed strategy. Finally, Section 5 concludes this paper.

2. System Model and Assumptions

In this work, we consider a multistatic radar system composed of Q R radars, whose purpose is to minimize the radiated power of each radar while guaranteeing a predetermined SINR requirement for target detection. As illustrated in Figure 1, the ith radar receives the echoes from the target due to its emitted waveforms as well as the waveforms from the other radars, both scattered off the target and through a direct path. The signals transmitted from different radars might not be orthogonal due to various reasons, including the absence of radar transmission synchronization [18], which could induce significant mutual interference. It is supposed that successive interference cancellation (SIC) technique is employed at each radar receiver to remove both direct and target scattered communication signals from the observed signal [26].
In the considered multistatic system, each radar performs target detection autonomously and sends its local decision to the fusion center, which takes a global decision once the data coming from all the radars is collected. It is also assumed that each radar can determine the presence of a target by employing a binary hypothesis testing on the received signal based on the generalized likelihood ratio test (GLRT) [18,20,26]. Thus, the M time-domain samples of the received signals for the ith radar, with H 0 corresponding to the target absence hypothesis and H 1 corresponding to the target presence hypothesis, can be expressed by:
H 0 : r i = j = 1 , j i Q R ξ i , j P j s j + n i , H 1 : r i = ϱ i P i s i + j = 1 , j i Q R ξ i , j P j s j + n i ,
where s i = ϕ i b i denotes the transmitted waveform from radar i, b i = [ 1 , e j 2 π f D , i , , e j 2 π ( M 1 ) f D , i ] T C M × 1 denotes the Doppler steering vector of radar i with respect to the target, f D , i is the Doppler shift associated with the radar i, M is the number of received pulses during the dwell time, and ϕ i is the predesigned waveform emitted from radar i. ϱ i represents the channel gain at the direction of the target, P i is the transmit power of radar i, ξ i , j stands for the cross gain between radar i and j, and n i denotes a zero-mean white Gaussian noise with variance σ n 2 . It is assumed that ϱ i CN ( 0 , h i , i T ) , ξ i , j CN ( 0 , c i , j ( h i , j T + h i , j D ) ) and n i CN ( 0 , σ n 2 ) , where h i , i t represents the variance of the channel gain for the radar i-target-radar i path, c i , j h i , j T represents the variance of the channel gain for the radar i-target-radar j path, c i , j h i , j D represents the variance of the channel gain for the direct radar i-radar j path, and c i , j denotes the cross correlation coefficient between the ith radar and jth radar.
The propagation gains of the corresponding paths are defined as follows:
h i , i T = G T G R σ i , i RCS λ 2 ( 4 π ) 3 R i 4 , h i , j T = G T G R σ i , j RCS λ 2 ( 4 π ) 3 R i 2 R j 2 , h i , j D = G T G R λ 2 ( 4 π ) 2 d i , j 2 , g i D = G T G I λ 2 ( 4 π ) 2 d i 2 ,
where h i , i T represents the propagation gain for the radar i-target-radar i path, h i , j T represents the propagation gain for the radar i-target-radar j path, h i , j D represents the direct radar i-radar j path, g i D represents the direct radar i-intercept receiver path. G T is the radar main-lobe transmitting antenna gain, G R is the radar main-lobe receiving antenna gain, G T is the radar side-lobe transmitting antenna gain, G R is the radar side-lobe receiving antenna gain, and G I is the interceptor receiving antenna gain. σ i , i RCS is the RCS of the target with respect to the ith radar, σ i , j RCS is the RCS of the target from radar i to radar j, λ denotes the wavelength, R i denotes the distance from radar i to the target, R j denotes the distance from radar j to the target, d i , j denotes the distance between radar i and radar j, and d i denotes the distance between radar i and intercept receiver. It is supposed that all the path propagation gains are fixed during observation period.

3. Problem Formulation

3.1. Stackelberg Game Formulation

Stackelberg game is a strategic game, which is composed of a single leader and multiple followers competing with each other [21,22]. In this study, the hostile intercept receiver plays a role of leader, who decides the prices on per unit of received power from different radars first through the maximization of its own profit. While the radars are multiple followers, which move subsequently and compete selfishly in a non-cooperative Nash game in view of the assigned prices. In the considered LPI optimization problem, multiple radars in the multistatic system are selfish, which act solely according to their own strategies. From the hostile interceptor’s point of view, those selfish moves lead to inefficient power resource utilization and degradation of the LPI performance of multistatic radar system. In this sequel, our propose is to formulate a Stackelberg game-theoretic power contorl strategy among different radars subject to a specified SINR requirement for target detection.

3.1.1. Leader-Level Game

Under the Stackelberg game model, the principal aim of the intercept receiver is to maximize its own utility from selling the received power to different radars. Mathematically, the utility of the intercept receiver can be written as:
u INT ( φ , P ) = i = 1 Q R φ i P i g i D · u S min i = 1 Q R P i g i d ,
where φ is the price vector with φ = [ φ 1 , , φ Q R ] T R Q R × 1 , P is the transmit power vector with P = [ P 1 , , P Q R ] T R Q R × 1 , φ i is the unit power price for the ith radar, u ( · ) denotes step function, and S min denotes the sensitivity of the interceptor. It should be pointed out that the transmit power of the ith radar P i depends on the imposed price φ i in the Stackelberg game formulation. To this end, the intercept receiver needs to find the best price φ i to maximize its own revenue. Hence, the optimization problem of the intercept receiver can be expressed as:
P 3.1 : max φ u INT ( φ , P ) , subject to : φ 0 .
where 0 = [ 0 , , 0 ] T R Q R × 1 .

3.1.2. Follower-Level Game

At the multistatic radar system’s side, the GLRT is adopted to determine the appropriate detector [26]. The probabilities of detection p D , i ( α i , γ i ) and false alarm p FA , i ( α i ) can be written as:
p D , i ( α i , γ i ) = 1 + α i 1 α i · 1 1 + M γ i 1 M , p FA , i ( α i ) = ( 1 α i ) M 1 ,
where α i is the detection threshold, M is the total number of received pulses during the dwell time. γ i denotes the SINR received at the ith radar expressed by:
γ i = h i , i t P i j = 1 , j i Q R c i , j h i , j D + h i , j T P j + σ n 2 = h i , i T P i I i ,
where σ n 2 represents the background noise at the ith radar, I i denotes the total interference and noise received at the ith radar, that is,
I i = j = 1 , j i Q R c i , j h i , j D + h i , j T P j + σ n 2 .
As mentioned before, the multiple radars take the part of followers to maximize their individual utilities through power control in a non-cooperative Nash game. Thus, the utility of each radar can be defined as:
u MRS , i ( P i , P i , φ i ) = ln ( γ i γ min ) Profit Part + P i max P i φ i P i g i D Cos t Part ,
where P i = [ P 1 , , P i 1 , P i + 1 , , P Q R ] T R ( Q R 1 ) × 1 denotes the vector of power allocation for all radars apart from radar i, γ min is the desired SINR threshold for target detection, and P i max is the maximum transmit power of radar i. From (5), it is noticed that the utility function of each radar is composed of profit part and cost part. If each radar increases its transmit power, the target detection performance improves, and thus the profit goes up. However, as the transmit power increases, the power consumption and the received power at intercept receiver are increased, and so does the cost [21,22]. Hence, there exists a trade-off between the profit and cost in each radar. Mathematically, for radar i, the optimization problem can be expressed by:
P 3.2 : max P u MRS , i ( P i , P i , φ i ) , subject to : γ i γ min , i , 0 P i P i max , i .
Therefore, the sub-problems P 3.1 and P 3.2 together form a Stackelberg game-theoretic model for the considered problem scenario. The purpose of the Stackelberg game is to find the SE points, which is studied in the following sub-section.

3.2. Analysis of the Proposed Stackelberg Game

For our formulated Stackelberg game model, the SE is defined as follows:
Definition 1.
The point ( φ * , P * ) is the SE for the considered Stackelberg game model if for any ( φ , P ) , the following conditions are satisfied [21,22]:
u INT ( φ * , P * ) u INT ( φ , P ) ,
u MRS , i ( φ * , P * ) u MRS , i ( φ , P ) ,
where φ * denotes the optimal solution for P 3.1 , and P * denotes the optimal solution for P 3.2 .
Generally speaking, the SE point for a Stackelberg game model can be achieved by obtaining the NE for the follower-level sub-game [21]. It has been shown that the different radars in multistatic radar system compete in a non-cooperative Nash game. On the other hand, the best response function of the intercept receiver can be obtained by solving the sub-problem P 3.1 . To achieve the SE point, the best response functions for the different radars (followers) should be obtained first, while the intercept receiver (leader) derives its best response function based on those of the radars subsequently.
Lemma 1.
For a given price of power resource φ i , the optimal solution for sub-problem P 3.2 can be expressed by:
P i ( k + 1 ) = γ m i n P i ( k ) γ i ( k ) + 1 1 2 P i m a x P i ( k ) + φ i ( k ) g i D 0 P i m a x ,
where { x } a b = m a x { m i n { x , b } , a } , and k denotes iteration index.
Proof. 
Taking the first derivative of u MRS , i ( P i , P i , φ i ) with respect to P i , we can obtain:
P i u MRS , i ( P i , P i , φ i ) = 1 γ i γ min h i , i T I i 1 2 P i max P i φ i g i D ,
Let u i ( p i , P i ) p i = 0 and rearrange terms, we have:
1 γ i γ min h i , i T I i = 1 2 P i max P i φ i g i D ,
thus,
γ i = γ min + h i , i T I i 1 1 2 P i max P i + φ i g i D .
Since γ i = P i h i , i T I i , then we can obtain:
P i = I i h i , i T γ min + 1 1 2 P i max P i + φ i g i D .
Finally, the following equation can be obtained to achieve the NE through iterations:
P i ( k + 1 ) = γ min P i ( k ) γ i ( k ) + 1 1 2 P i max P i ( k ) + φ i ( k ) g i D 0 P i max ,
which completes the proof. □
Proposition 1.
The non-cooperative Nash game model P 3.2 has at least one NE.
Proof. 
According to [19,26], the conditions for the existence of NE are listed as follows:
(i) 
P i is a non-null, convex and tight subset in a finite Euclidean space;
(ii) 
u MRS , i ( P i , P i , φ i ) is continuous and quasi-concave with P i .
Each power element is limited between 0 and P i max , and P i is a convex tight subset. Thus, the condition (i) is easily met. For condition (ii), we take the second order derivative of u MRS , i ( P i , P i , φ i ) with respect to P i and obtain
2 P i 2 u MRS , i ( P i , P i , φ i ) = h i , i T I i 2 2 1 ( γ i γ min ) 3 1 4 ( P i max P i ) 3 < 0 .
Thus, u MRS , i ( P i , P i , φ i ) is a concave function of P i . Both conditions (i) and (ii) hold. Therefore, there exists at least one NE in P 3.2 , which completes the NE existence proof. □
Proposition 2.
The NE of the non-cooperative Nash game model P 3.2 is unique.
Proof. 
In order to show that the NE of the game model P 3.2 is unique, we need to prove that the ith radar’s best response strategy function
y ( P i ) = γ min P i γ i + 1 1 2 P i max P i + φ i g i D
should be standard, which satisfies the following conditions [26]:
(i) 
Positivity: For i , y ( P i ) > 0 ;
(ii) 
Monotonicity: If P i m > P i n , then y ( P i m ) > y ( P i n ) ;
(iii) 
Scalability: For ζ > 1 , ζ y ( P i ) > y ( ζ P i ) .
For Condition (i), since for i , it is obvious to obtain
y ( P i ) = γ min P i γ i + 1 1 2 P i max P i + φ i g i D > 0 ,
thus Condition (i) is satisfied.
For Condition (ii), if P i m > P i n , we can obtain:
y ( P i m ) y ( P i n ) = γ min γ i P i max P i n + P i max P i m W · P i max P i n + P i max P i m ,
where
W 1 1 2 P i max P i m + φ i g i D 1 2 P i max P i n + φ i g i D × 1 2 P i max P i m P i max P i n .
Since we have i = 1 Q R P i g i D S min , that is:
P i S min j = 1 , j i Q R P j g j D / g i D η .
Thus:
P i max P i n + P i max P i m 2 P i max η .
When γ min γ i · 2 P i max η > W , we can obtain:
y ( P i m ) > y ( P i n ) .
In such a case, Condition (ii) is satisfied.
For Condition (iii),
ζ y ( P i ) y ( ζ P i ) = V × ζ P i max P i P i max ζ P i 2 P i max P i P i max ζ P i + ζ φ i g i D ,
where
V 1 1 2 P i max P i + φ i g i D 1 2 P i max ζ P i + φ i g i D .
Owing to ζ > 1 , it is apparent that:
ζ P i max P i P i max ζ P i > 0 .
Therefore, we can obtain:
ζ y ( P i ) y ( ζ P i ) > 0 ,
and Condition (iii) is satisfied.
As a result, all the above conditions are met, which completes the NE uniqueness proof. □

3.3. Distributed Approach for Calculating Stackelberg Equilibrium

Based on the above theoretical derivations and analyses, we develop a distributed iterative power control method to perform LPI performance optimization for multistatic radar system. As aforementioned, to achieve the SE point, the sub-problem P 3.2 must be solved first for a given price of power resource φ i . Then, we solve the sub-problem P 3.1 for the optimal price φ i * with the calculated transmit power P i * of multiple radars. The detailed steps of the pricing-based distributed iterative power control approach is provided in Algorithm 1, where ( · ) * is the SE solution, and Δ φ is the step size.
Algorithm 1: Pricing-Based Distributed Iterative Power Control Approach
Electronics 08 00397 i001

4. Numerical Examples and Performance Evaluation

In this section, some numerical examples are presented to evaluate the performance of our proposed LPI optimization strategy. For ease of exposition, we consider a target detection scenario similar to Figure 1, which consists of a multistatic radar system, a target, and a intercept receiver. It is assumed that the multistatic radar system consists of Q R = 6 widely deployed radar nodes, which are located at ( 50 , 0 ) km , ( 25 , 25 / 3 ) km , ( 25 , 25 / 3 ) km , ( 50 , 0 ) km , ( 25 , 25 / 3 ) km , and ( 25 , 25 / 3 ) km , respectively. The position of the intercept receiver is ( 30 , 30 ) km . In each time slot, each radar receives M = 512 pulses. The probabilities of target detection and false alarm are set as p D , i = 0.9973 ( i ) and p FA , i = 10 6 ( i ) , respectively. Thus, the corresponding detection threshold α i ( i ) and the SINR threshold γ min can be computed as 0.0267 and 10 dB, respectively. The other system parameters are given in Table 1.
To evaluate the influence of target reflectivity on the power allocation results, we consider two different target RCS models. In the first case, σ RCS , 1 = [ 1 , 1 , 1 , 1 , 1 , 1 ] m 2 . In the second case, σ RCS , 2 = [ 0.5 , 20 , 5 , 0.25 , 16 , 50 ] m 2 . On the other hand, in order to reveal the effect of system configuration on resource allocation results, we consider two different target locations, that is, [ 0 , 0 ] km and [ 30 , 40 ] km .
Figure 2 depicts the convergence behavior for the transmit power of each radar in the first case with different initial power allocation values P ( 1 ) = [ 3600 , 800 , 300 , 1500 , 600 , 4000 ] W and P ( 1 ) = [ 100 , 2300 , 1500 , 280 , 4600 , 900 ] W , respectively. Similarly, Figure 3 shows the convergence behavior for the transmit power of each radar in the second case with P ( 1 ) = [ 2500 , 2500 , 2500 , 2500 , 2500 , 2500 ] W and P ( 1 ) = [ 1000 , 1500 , 1000 , 1500 , 1000 , 1500 ] W , respectively. The proposed algorithm stops when u INT ( k + 1 ) u INT ( k ) < ϵ is within the desired accuracy. It can be observed that our presented Stackelberg game-based LPI performance optimization strategy converges fast to the unique NE points for all initial values of transmit power.
In Figure 4 and Figure 5, we illustrate the convergence process for the transmit power allocation in both cases. Here, we define the ratio of transmit power as κ i = P i / i = 1 Q R P i [21,22]. One can see from Figure 3b that the proposed LPI optimization scheme would like to assign more transmit power to Radar 1, Radar 2 and Radar 6, which are farther from the target. In addition, from Figure 3a, it can be noticed that more power resource is allocated to Radar 4 and Radar 1, whose RCS with respect to the target is much smaller than other radars. Hence, we can conclude that the transmit power allocation results depend on the system geometry between target and multistatic radar and target reflectivity [26]. Furthermore, as shown in Figure 3b, the proposed optimization algorithm tends to distribute less transmit power to the radars with larger path propagation gains.
Figure 6 and Figure 7 show the convergence behavior for the achieved SINR of each radar in both cases. As can be interpreted in these figures, by employing our proposed strategy, the achieved SINR value of each radar tends to converge to the predefined SINR threshold after 6–8 iterations, and thus the desired target detection requirement can be satisfied.
Moreover, the convergence performance for the normalized utility function of intercept receiver in both scenarios is illustrated in Figure 8 and Figure 9. The results indicate that, as the number of iterations increase, the revenues of intercept receiver eventually converge to the SE points in both cases.
Figure 10 and Figure 11 show the convergence behavior for the received power at intercept receiver. It is evident that the total transmit power received at intercept receiver from the multistatic radar system is below the sensitivity of interceptor S min . This is due to the fact that the proposed strategy can coordinate the radiated power from the radar transmissions through updating the unit power prices. In such a case, the LPI performance of the multistatic radar system can be guaranteed by minimizing the power consumption of each radar. To conclude, our proposed Stackelberg game-theoretic LPI optimization scheme is effective to enhance the LPI performance of the multistatic radar system.

5. Conclusions

This paper studies the problem of Stackelberg game-theoretic LPI performance optimization in multistatic radar system, whose purpose is to minimize the radiated power of each radar for a specified target detection performance. A single-leader multi-follower Stackelberg game is established to systematically formulate the intercept receiver and the multiple radars’ behaviors, where the hostile intercept receiver acts as a leader and the multiple radars are followers, respectively. The Stackelberg game model jointly investigates the utility maximization of the interceptor and multiple radars. Based on our theoretical findings on the existence and uniqueness of the NE in the game model, we present a pricing-based distributed iterative power control method, and its convergence to the NE is verified by numerical simulations. Moreover, the results of this work are useful to practically strengthen the LPI performance for target detection in multistatic radar system. In future work, the derivations and simulation results will be extended to the cooperative-theoretic case.

Author Contributions

C.S., W.Q. and F.W. conceived and designed the experiments; C.S., W.Q. and F.W. performed the experiments; S.S. and J.Z. analyzed the data; C.S. wrote the paper; S.S. and J.Z. contributed to data analysis revision; S.S. and J.Z. contributed to English language correction. All authors of article provided substantive comments.

Funding

This research has recieved no external funding.

Acknowledgments

This work is supported in part by the National Natural Science Foundation of China (Grant No. 61801212, No. 61671239), in part by the Natural Science Foundation of Jiangsu Province (Grant No. BK20180423), in part by the National Aerospace Science Foundation of China (Grant No. 20172752019, No. 2017ZC52036), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PADA) and Key Laboratory of Radar Imaging and Microwave Photonics (Nanjing Univ. Aeronaut. Astronaut.), Ministry of Education, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China, and in part by the Key Laboratory of Dynamic Cognitive System of Electromagnetic Spectrum Space (Nanjing Univ. Aeronaut. Astronaut.), Ministry of Industry and Information Technology, Nanjing 210016, China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pace, P.E.; Tan, C.K.; Ong, C.K. Microwave-photonics direction finding system for interception of low probability of intercept radio frequency signals. Opt. Eng. 2018, 57, 1–8. [Google Scholar] [CrossRef]
  2. Zhang, Z.K.; Tian, Y.B. A novel resource scheduling method of netted radars based on Markov decision process during target tracking in clutter. EURASIP J. Adv. Signal Process. 2016, 16, 1–9. [Google Scholar] [CrossRef]
  3. Zhang, Z.K.; Salous, S.; Li, H.L.; Tian, Y.B. Optimal coordination method of opportunistic array radars for multi-target-tracking-based radio frequency stealth in clutter. Radio Sci. 2015, 50, 1187–1196. [Google Scholar] [CrossRef] [Green Version]
  4. Shi, C.G.; Wang, F.; Salous, S.; Zhou, J.J. Joint transmitter selection and resource management strategy based on low probability of intercept optimization for distributed radar networks. Radio Sci. 2018, 53, 1108–1134. [Google Scholar] [CrossRef]
  5. Wang, X.L.; Yi, W.; Xie, M.C.; Kong, L.J. A joint beam and dwell time allocation strategy for multiple target tracking based on phase array radar system. In Proceedings of the 20th International Conference on Information Fusion (Fusion), Xi’an, China, 10–13 July 2017; pp. 1–5. [Google Scholar]
  6. Wang, X.L.; Yi, W.; Xie, M.C.; Zhai, B.W.; Kong, L.J. Time management for target tracking based on the predicted Bayesian Cramer-Rao lower bound in phase array radar system. In Proceedings of the 20th International Conference on Information Fusion (Fusion), Xi’an, China, 10–13 July 2017; pp. 1–5. [Google Scholar]
  7. Cheng, T.; Zou, D.Q.; He, Z.S. Adaptive waveform and sampling interval tracking based on estimation accuracy for Doppler radar. In Proceedings of the 2013 IET International Radar Conference, Xi’an, China, 14–16 April 2013; pp. 1–4. [Google Scholar]
  8. Lawrence, D.E. Low probability of intercept antenna array beamforming. IEEE Trans. Antennas Propag. 2010, 58, 2858–2865. [Google Scholar] [CrossRef]
  9. Zhou, C.W.; Gu, Y.J.; He, S.B.; Shi, Z.G. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2018, 67, 1099–1112. [Google Scholar] [CrossRef]
  10. Xiong, J.; Wang, W.Q.; Cui, C.; Gao, K.D. Cognitive FDA-MIMO radar for LPI transmit beamforming. IET Radar Sonar Navig. 2017, 11, 1574–1580. [Google Scholar] [CrossRef]
  11. Shi, C.G.; Wang, F.; Sellathurai, M.; Zhou, J.J.; Salous, S. Power minimization-based robust OFDM radar waveform design for radar and communication systems in coexistence. IEEE Trans. Signal Process. 2018, 66, 1316–1330. [Google Scholar] [CrossRef]
  12. Fang, H.; Xu, L.; Wang, X.B. Coordinated multiple-relays based physical-layer security improvement: A single-leader multiple-followers Stackelberg game scheme. IEEE Trans. Inf. Forensics Secur. 2018, 13, 197–209. [Google Scholar] [CrossRef]
  13. Du, Z.Y.; Wu, Q.H.; Jiang, B.; Xu, Y.H.; Qin, Z.Q. Interference-aware spectrum access self-organization: A weighted graph game perspective. IEEE Syst. J. 2018, 12, 3250–3259. [Google Scholar] [CrossRef]
  14. Zhang, Y.L.; Xu, Y.H.; Wu, Q.H. Group buying based on social aware in D2D networks: A game theoretic approach. In Proceedings of the IEEE/CIC International Conference on Communications in China (ICCC), Qingdao, China, 22–24 October 2017; pp. 1–6. [Google Scholar]
  15. Sun, B.; Chen, H.W.; Wei, X.Z.; Wang, H.Q.; Li, X. Power allocation for range-only localisation in distributed multiple-input multiple-output radar networks-a cooperative game approach. IET Radar Sonar Navig. 2014, 8, 708–718. [Google Scholar] [CrossRef]
  16. Song, X.F.; Willett, P.; Zhou, S.L.; Luh, P.B. The MIMO radar and jammer games. IEEE Trans. Signal Process. 2012, 60, 687–699. [Google Scholar] [CrossRef]
  17. Piezzo, M.; Aubry, A.; Buzzi, S.; Maio, A.D.; Farina, A. Non-cooperative code design in radar networks: A game-theoretic approach. EURASIP J. Adv. Signal Process. 2013, 2013, 63. [Google Scholar] [CrossRef]
  18. Deligiannis, A.; Panoui, A.; Lambotharan, S.; Chambers, J.A. Game-theoretic power allocation and the Nash equilibrium analysis for a multistatic MIMO radar network. IEEE Trans. Signal Process. 2017, 65, 6397–6408. [Google Scholar] [CrossRef]
  19. Shi, C.G.; Wang, F.; Sellathurai, M.; Zhou, J.J. Non-cooperative game theoretic power allocation strategy for distributed multiple-radar architecture in a spectrum sharing environment. IEEE Access 2018, 6, 17787–17800. [Google Scholar] [CrossRef]
  20. Shi, C.G.; Wang, F.; Salous, S.; Zhou, J.J.; Hu, Z.T. Nash bargaining game-theoretic framework for power control in distributed multiple-radar architecture underlying wireless communication system. Entropy 2018, 20, 267. [Google Scholar] [CrossRef]
  21. Kang, X.; Zhang, R.; Motani, M. Price-based resource allocation for spectrum-sharing femtocell networks: A Stackelberg game approach. IEEE J. Sel. Areas Commun. 2012, 30, 538–549. [Google Scholar] [CrossRef]
  22. Wang, W.B.; Hoang, D.T.; Niyato, D.; Wang, P.; Kim, D.I. Stackelberg game for distributed time scheduling in RF-powered backscatter cognitive radio networks. IEEE Trans. Wirel. Commun. 2018, 17, 5606–5622. [Google Scholar] [CrossRef]
  23. Yin, R.; Zhong, C.J.; Yu, G.D.; Zhang, Z.Y.; Wong, K.K.; Chen, X.M. Joint spectrum and power allocation for D2D communications underlying cellular networks. IEEE Trans. Veh. Technol. 2016, 65, 2182–2195. [Google Scholar] [CrossRef]
  24. Yuan, Y.L.; Liang, C.C.; Kaneko, M.; Chen, X.; Hogrefe, D. Topology control for energy-efficient localization in mobile underwater sensor networks using Stackelberg game. IEEE Trans. Veh. Technol. 2019, 68, 1487–1500. [Google Scholar] [CrossRef]
  25. Lan, X.; Li, W.; Wang, X.L.; Yan, J.B.; Jiang, M.R. MIMO radar and target Stackelberg game in the presence of clutter. IEEE Sens. J. 2015, 15, 6912–6920. [Google Scholar] [CrossRef]
  26. Shi, C.G.; Wang, F.; Sellathurai, M.; Zhou, J.J. Non-cooperative game-theoretic distributed power control technique for radar network based on low probability of intercept. IET Signal Process. 2018, 12, 983–991. [Google Scholar] [CrossRef]
Figure 1. Illustration of the system model.
Figure 1. Illustration of the system model.
Electronics 08 00397 g001
Figure 2. Convergence behavior for the transmit power of each radar in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 2. Convergence behavior for the transmit power of each radar in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g002
Figure 3. Convergence behavior for the transmit power of each radar in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 3. Convergence behavior for the transmit power of each radar in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g003
Figure 4. Convergence behavior for the transmit power allocation in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 4. Convergence behavior for the transmit power allocation in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g004
Figure 5. Convergence behavior for the transmit power allocation in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 5. Convergence behavior for the transmit power allocation in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g005
Figure 6. Convergence behavior for the achieved SINR of each radar in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 6. Convergence behavior for the achieved SINR of each radar in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g006
Figure 7. Convergence behavior for the achieved SINR of each radar in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 7. Convergence behavior for the achieved SINR of each radar in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g007
Figure 8. Convergence behavior for the normalized utility of intercept receiver in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 8. Convergence behavior for the normalized utility of intercept receiver in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g008
Figure 9. Convergence behavior for the normalized utility of intercept receiver in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 9. Convergence behavior for the normalized utility of intercept receiver in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g009
Figure 10. Convergence behavior for the received power at intercept receiver in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 10. Convergence behavior for the received power at intercept receiver in the first case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g010
Figure 11. Convergence behavior for the received power at intercept receiver in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Figure 11. Convergence behavior for the received power at intercept receiver in the second case: (a) Target position is [ 0 , 0 ] km ; (b) Target position is [ 30 , 40 ] km .
Electronics 08 00397 g011
Table 1. Coexisting DMRS and Wireless Communication System Parameters.
Table 1. Coexisting DMRS and Wireless Communication System Parameters.
ParameterValueParameterValue
G T 30 dB G R 30 dB
G T 30 dB G R 30 dB
G I 0 dB σ n 2 10 18 W
S min 102 dBmW c i , j 0.01
λ 0.10 m P i , max ( i ) 5000 W

Share and Cite

MDPI and ACS Style

Shi, C.; Qiu, W.; Wang, F.; Salous, S.; Zhou, J. Stackelberg Game-Theoretic Low Probability of Intercept Performance Optimization for Multistatic Radar System. Electronics 2019, 8, 397. https://doi.org/10.3390/electronics8040397

AMA Style

Shi C, Qiu W, Wang F, Salous S, Zhou J. Stackelberg Game-Theoretic Low Probability of Intercept Performance Optimization for Multistatic Radar System. Electronics. 2019; 8(4):397. https://doi.org/10.3390/electronics8040397

Chicago/Turabian Style

Shi, Chenguang, Wei Qiu, Fei Wang, Sana Salous, and Jianjiang Zhou. 2019. "Stackelberg Game-Theoretic Low Probability of Intercept Performance Optimization for Multistatic Radar System" Electronics 8, no. 4: 397. https://doi.org/10.3390/electronics8040397

APA Style

Shi, C., Qiu, W., Wang, F., Salous, S., & Zhou, J. (2019). Stackelberg Game-Theoretic Low Probability of Intercept Performance Optimization for Multistatic Radar System. Electronics, 8(4), 397. https://doi.org/10.3390/electronics8040397

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