Next Article in Journal
An Open Data Approach for Estimating Vegetation Gross Primary Production at Fine Spatial Resolution
Next Article in Special Issue
Reflective Tomography Lidar Image Reconstruction for Long Distance Non-Cooperative Target
Previous Article in Journal
Worldwide Statistical Correlation of Eight Years of Swarm Satellite Data with M5.5+ Earthquakes: New Hints about the Preseismic Phenomena from Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Antenna Placement and Power Allocation for Target Detection in a Distributed MIMO Radar Network

Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(11), 2650; https://doi.org/10.3390/rs14112650
Submission received: 14 April 2022 / Revised: 21 May 2022 / Accepted: 30 May 2022 / Published: 1 June 2022

Abstract

:
Radar network configuration and power allocation are of great importance in military applications, where the entire surveillance area needs to be searched under resource budget constraints. To pursue the joint antenna placement and power allocation (JAPPA) optimization, this paper develops a JAPPA strategy to improve target detection performance in a widely distributed multiple-input and multiple-output (MIMO) radar network. First, the three variables of the problem are incorporated into the Neyman–Pearson (NP) detector by using the antenna placement optimization and the Lagrange power allocation method. Further, an improved iterative greedy dropping heuristic method based on a two-stage local search is proposed to solve the NP-hard issues of high-dimensional non-linear integer programming. Then, the sum of the weighted logarithmic likelihood ratio test (LRT) function is constructed as optimization criteria for the JAPPA approach. Numerical simulations and the theoretical analysis confirm the superiority of the proposed algorithm in terms of achieving effective overall detection performance.

1. Introduction

1.1. Background and Related Studies

The detection theory development and system optimization utilization of distributed MIMO radars is one of the most important directions in metrology for radar. Multi-Input Multi-Output (MIMO) radar networks have recently received significant attention due to their superior detection and tracking performance compared to traditional monostatic radar systems [1,2,3]. They also possess superior anti-stealth, anti-destroy, and anti-interference performance by cooperating multiple subarrays with transmitting independent and orthogonal waveforms [4,5]. As an important category of the MIMO radar network [6], the distributed MIMO radar has enhanced detection and high localization precision because its multiple transmitters and multiple receivers are sufficiently far from each other [7,8,9,10].
Among the research topics of MIMO radars, radar resource allocation is an important pot in the military field. It plays a great promotion role in the rapid deployment of self-propelled radar equipment and air defense planning in key areas, which are critical military needs with significant applications. In particular, it is utilized in air defense operations to enhance combat efficiency within the restrictions of complex battle characteristics, such as electromagnetics, topography, and hostile situations [11,12].
Antenna position deployment and transmitting power allocation are two factors concerned by the radar resource allocation, and they determine the performance of MIMO radars. However, it is a vital challenge to maximize resource consumption efficiency and radar network performance. Various research has been devoted to them. According to the focus of the method, the related studies can be classified into three categories: the system configurations [13,14,15], the transmitting power allocation [16,17,18], and the integration of the former two [19,20,21,22].
The first category is to optimize the system configuration to enhance the performance of the radar system. Considering that the antenna placement problem is NP hard [23], traditional methods for the solution are the exhaustive search and heuristic search algorithms. However, they are either computationally intensive or only performed at a local–optimal level. As a result, the greedy algorithm [13] and intelligent algorithm [14,15] are often adopted to solve the problem. Specifically, ref. [24] proposes simulation-based discrete stochastic optimization algorithms to adaptively select a better receiving antenna subset using target tracking resolution. The effects of optimizing the remaining antennas’ placement on the radar performance are studied when fixing receiving antennas and transmitting antennas separately. When optimization is extended to consider both transmitting and receiving antennas, ref. [25] investigates the target velocity estimation performance of the distributed MIMO radar using the Cramér–Rao bound and finds that increasing the number of antennas and better placement can lead to better speed estimation performance. A combinatorial optimization model for the joint placement of transmitters and receivers is established in [13], and the proposed algorithm utilizes convex relaxation to obtain the approximation of the original optimal solution. Further, [26] proposes an antenna placement scheme for compressed-sensing-based collocated MIMO radars to improve the target detection performance by minimizing the sensing matrix’s coherence. An optimization algorithm is proposed for iterative placement between transmitters and receivers to obtain better detection performance. Therefore, geometry configuration is crucial to the optimization on distributed MIMO radar systems that achieves a superior performance.
The second category is to adaptively optimize the transmitting power allocation. In principle, the optimal power allocation strategy is to obtain a solution based on the convex optimization theory. For example, an alternating global search algorithm is designed in ref. [16] to minimize the non-convex Cramér–Rao low bound (CRLB) of target position estimation through power allocation, and ultimately, enhance the power utilization of the distributed radar systems. Similarly, ref. [17] presents a deterministic convex power allocation scheme and transforms it into a non-linear quality model with the Karush–Kuhn–Tucker conditions. Ref. [18] proposes a novel low probability of intercept power resource allocation (LPI-PRA) scheme that delivers superior low probability of intercept (LPI) performance in terms of minimizing the total power consumption by exploiting the semi-definite programming (SDP) and Karush–Kuhn–Tucker (KKT) conditions. Further, ref. [27] studies the robust power allocation problem for multi-target tracking in the quality of service framework and proposes an iterative parallel search algorithm to solve such a non-convex optimization problem. Hence, the flexibility and excellent performance of power allocation are shown in the superior performance of the improved radar network.
The last category is joint optimization of the former two items. The research on the joint antenna placement and power allocation in MIMO radars also yielded fruitful results. For instance, ref. [28] proposes a suboptimal, but computationally efficient, method for scheduling and power allocation based on greedy programming in multi-target tracking. It improves tracking accuracy by adopting the suitable subset of antennas to be employed in each tracking interval, as well as the power transmitted by these antennas. In addition, ref. [19] establishes the SNR criteria for the antenna placement and power allocation optimization. It discretizes the radar placement area into small grids, and then a sequential exhaustive enumeration geometric method and the water-filling power distribution algorithm are used to optimize variables, respectively. Different from [19], ref. [20] minimizes the total transmit power by the Lagrangian method and domain decompositions. Both [19] and [20] enhance the detection performance of the radar system; however, the antenna placement and power allocation in optimization processes are often optimized separately according to different objective functions, resulting in deviations from the global optimum.
For this reason, some collaborative methods may be suitable for solving such joint optimization problems. For instance, the combination of distributed MIMO radar nodes selection and power allocation based on target localization accuracy is investigated in ref. [21]. It utilizes the performance matrix as the single criterion for a two-step optimization method and converts it into a second-order cone problem to solve. Further, ref. [22] studies the joint optimization problem of tracking multiple targets and it uses predicted conditional Cramér–Rao lower bound (PC-CRLB) as an optimization criterion for a two-step semidefinite-programming-based solution to solve this problem. Apart from the previous works, ref. [29] formulates the cost function by minimizing the sum of weighted PC-CRLB, integrating the subarray selection with power allocation optimization by a two-stage local search-based algorithm. More recently, ref. [30] develops a collaborative radar node selection and transmitter resource allocation algorithm for target tracking applications in multiple radar architectures. An efficient two-stage-based solution methodology that incorporates the interior point method and cyclic minimization framework makes it converge to the optimal solution.

1.2. Contributions

In light of the aforementioned literature, joint optimization has outstanding effectiveness in improving the radar system performance. However, few studies that concentrate on the problem of joint radar antenna placement and transmitter resource allocation (JAPPA) with the objective of target detection have been investigated. This gap promotes this research.
In the current paper, the JAPPA under a collaboratively unified objective function is investigated by considering radar detection performance optimization. Our primary contributions are summarized as follows.
  • We formulate the joint power allocation and antenna placement problem as an optimization model subject to the resource budget and area priority level. In the proposed JAPPA model, binary composite hypothesis testing is established to design the Neyman–Pearson (NP) detector for the whole surveillance area with the targets Radar Cross Section (RCS) obeying Rayleigh distribution. Compared to previous oversimplified works [19,20], the weighted NP-based logarithm likelihood ratio test (LRT) function is specified as the utility function that combines all optimization factors into a single objective function. In addition, we evaluate the average utility function values after changing the single stationary target’s position in the whole region to describe the global target detection performance of the radar network.
  • We propose an efficient JAPPA strategy that incorporates antenna placement with power allocation to optimize the target detection performance in the radar network. Different from the 0–1 programming in [21,31], we choose suitable antenna deployment positions in the regional grid points to establish the non-linear mixed integer programming problem. To devise computationally feasible methods for practical applications, a two-stage local-search-based algorithm is proposed to split the coupled joint optimization. Herein, it isolates integer programming from continuous variable optimization and provides equivalent performance while requiring less computing effort.
  • We develop a joint optimization closed-loop scheme for the joint antenna placement and power allocation optimization [30]. In our scheme, the optimal position results in the current cycle are used for guiding the power allocation. The detection performance will be further improved through this link, which renders a closed-loop scheme to repeat the iteration of transmitting and receiving antenna placement using the optimal power allocation scheme.

1.3. Organization

The remainder of the paper is organized as follows: Section 2 establishes the system model and NP-based detection model. Section 3 develops the optimization model and proposes an efficient solver. Based on the solver, a two-stage closed-loop method for resource allocation is built. Simulation results and contrast test data are presented in Section 4 to verify the effectiveness and superiority of the proposed method. Finally, Section 5 concludes the paper’s results and achievements.
Notations: For clarity, some notations used in this paper are listed as follows, and others will be explained when they first appear: The boldface is used for vectors s (lower case) and matrices A (upper case). Superscript ( ) T denotes the transpose operators. N ( a , b ) represents the Gaussian normal distribution where the mean is a , and the variance is b . P T , k m indicates the allocated power to the m-th transmitting antenna in the k-th iteration loop. Moreover, L T , k and L R , k denote the transmitting antenna and receiving antenna locations in the k-th iteration loop, respectively.

2. System Model and Detection Formulation

2.1. Signal Model

Note that the detection performance of a distributed MIMO radar is an essential indicator for evaluating radar network surveillance performance. This paper discusses the joint antenna placement and optimal power allocation method in a two-dimensional Cartesian coordinate system. We consider a distributed MIMO radar network with M transmitters and N receivers. The m-th transmitter is positioned at the coordinate ( x T m , y T m ) for m = 1 , , M , the n-th receiver is positioned at the coordinate ( x R n , y R n ) for n = 1 , , N , respectively. Besides, the position of the target is ( x t , y t ) . Figure 1 illustrates the considered distributed MIMO radar network configuration sketch with M = N = 3.
The MIMO radar transmitting array elements simultaneously transmit a set of normalized orthogonal signals to the far-field target, and the precoding waveform vector satisfies
s m ( t ) s m H ( t τ ) d t = { 1 for m = m , τ = 0 0 for otherwise
where s m ( t ) denotes the corresponding baseband signal, f c is the signal carrier frequency, and τ m is the duration time of the m-th signal.
As such, the baseband signal transmitted from the m-th transmitting antenna to the n-th receiving antenna can be expressed as:
r n m = P T m L n m C n m s m ( t τ n m ) e j 2 π f d , n m t + n n ( t )
P T m L n m = [ P T 1 L 11 P T 2 L 12 P T M L 1 M P T N L N 1 P T 2 L N 2 P T M L N M ] N × M
C n m = [ c 11 c 12 c 1 N c M 1 c M 2 c M N ] M × N
where L n m is the propagation that includes antenna gains, P T m is the radar transmit power, C n m is scattering amplitude that includes any unknown oscillator phase terms, s m is the signal transmitted by the m-th transmitting antenna, τ n m is the time-delay, f d , n m is the Doppler shift caused by target movement. Moreover, n n ( t ) presents the additive complex white Gaussian noise vector at the receiver.
With the assumptions above, the baseband signal received by the n-th receiving antenna of the distributed MIMO radar can be written as:
r n ( t ) = m = 1 M r m n ( t ) = m = 1 M E n m c n m s m ( t τ n m ) e j 2 π f d , n m t + n n ( t )
where n n ( t ) is the zero-mean white complex noise at n-th receiver, and
E n m E n = [ P T 1 L n 1 , P T 2 L n 2 , , P T M L n M ] 1 × M
c n m C n = [ c n 1 , c n 2 , , c n M ] 1 × M

2.2. Neyman–Pearson Detection Model

The decision of target existence in a distributed MIMO radar network is made without any prior probability. Therefore, the NP criterion is suitable for radar detection [31]. Herein, we develop the NP detector for a distributed MIMO radar network in this section, seeking the optimal antenna position and power allocation scheme.
The NP detector establishes the LRT and compares it with the threshold η specified by the P f a . It yields two likelihood hypotheses, as shown in (8): the null hypothesis (signal-absence) H 0 and the alternative hypothesis (signal-presence) H 1 . It is noticeable that distributed MIMO radar network detection hypotheses are given for each receiver.
{ H 0 : r n = n n ( t )   H 1 : r n = m = 1 M E n m c n m s m ( t τ n m ) e j 2 π f d , n m t + n n ( t )   , n = 1 , , N
In this paper, the target comprises multiple randomly distributed independent scatterers of roughly equal areas. Hence, the target RCS is assumed under the Rayleigh scattering model, and the Swerling I model is defined as [32]:
f σ t = x σ 0 2 e x 2 2 σ 0 2
where σ 0 2 is the average RCS of the object. According to the free space propagation loss equation analysis, the amplitude of the signal received by the n-th receiver can be cast as (10) [32].
α n m = σ t R m t R t n P T m G t G r I p λ 2 ( 4 π ) 3 L c L r
where P T m is the transmitted power, σ t is RCS of the target, G t and G r are the transmitter and receiver gains, respectively, I p is the receiving processing gain, λ is radar signal wavelength, L c and L r are the scattering loss and receiving loss, respectively. R m t and R t n are the distance between the m-th transmitter to target and target to the n-th receiver. Due to the widely separated antennas, different α n m s are regarded as independent and the echo amplitude is proportional to the RCS, resulting in a Rayleigh distribution for amplitude α n m . Substituting (10) in (9), it can be deduced that
f α n m = x σ n m 2 e x 2 2 σ n m 2
where
σ n m 2 = P T m G r G t I p σ 0 2 λ 2 ( 4 π ) 3 R m t 2 R t n 2 L c L r
Therefore, we have
α n m = E n m c n m e j 2 π f d , n m t
The received signal r n can be redefined as
r n = m = 1 M α n m s m ( t τ n m ) e j 2 π f d , n m t + n n ( t )
Then, a likelihood ratio detection session is performed after the echo signal passes through the filter bank and sampler. Thus, a binary hypothesis can be expressed as follows:
y n = { m = 1 M a n m s m ( t τ n m ) + n n , H 1 n n , H 0
where y n = [ y 1 , y 1 , , y n ] T represents the received signal test data vector, s m is the sample version signal to be detected, n n is the noise in the test. In addition, we assume the noise model is a Gaussian distributed model that satisfies n n ~ N ( 0 , σ 2 ) , then the joint probability density function (PDF) of the test data y n is defined as
{ f Y n ( y n 0 ) N ( 0 , σ 2 I N ) f Y n ( y n 1 , α n m ) N ( r 0 n , σ 2 I N )
In the above equation, r 0 n m = 1 M a n m s m ( t τ n m ) , and σ 2 I N denotes the noise variance and identity matrix I N = d i a g ( 1 , , 1 ) 1 × N . Then, the LRT for the NP detector is usually given as follows [33]
L ( y n ) = f Y n ( y n | 1 ) f Y n ( y n | 0 ) > 1 < 0 η
where η is a threshold determined by P f a . Based on this derivation, hypothesis 0 and 1 probability density function models are, respectively, defined as:
f Y n ( y n | 0 ) = n = 1 N 1 2 π σ 2 e x p { 1 2 ( y n 2 σ 2 ) }
f Y n ( y n | 1 , α m n ) = n = 1 N 1 2 π σ 2 e x p { 1 2 ( y n r 0 n 2 σ 2 ) }
According to the specific PDF in (16) and substituting from (11), (19) can be derived as:
f Y n ( y n | 1 ) = 0 ( M × N 2 ) 0 f Y n ( y n | 1 , α n m ) m = 1 M N = 1 N f α n m ( α n m )   d α n m = 0 ( M × N 2 ) 0 { 1 ( t τ n m ) N / 2 e x p n = 1 N { 1 2 ( y n 2 2 y n T r 0 n + r 0 n 2 σ 2 ) } } m = 1 M N = 1 N f α n m ( α n m )   d α n m
Based on the previous Formulas (1) and (10)–(14), we can make a substitution due to the orthogonality of the distributed radar signals that
r 0 n 2 = r 0 n T r 0 n = m = 1 M m = 1 M α n m s T ( t τ n m ) α m n s ( t τ n m ) = m = 1 M α m n 2
and
y n 2 = y n T y n = 1 M m = 1 M y n 2
2 y n T r 0 n = 2 α m n m = 1 M y n T s m ( t τ n m )
Substituting (21)–(23) into (20), it is easy to find that
f Y n ( y n | 1 ) = 0 ( M × N 2 ) 0 { 1 ( 2 π σ 2 ) N / 2 e x p { 1 2 ( 1 M n = 1 N m = 1 M y n 2 2 α m n y n T s m ( t τ n m ) + α m n 2 σ 2 ) } } × m = 1 M n = 1 N f α m n ( α m n )   d α m n = 1 ( 2 π σ 2 ) N / 2 m = 1 M n = 1 N ( 0 e x p ( n = 1 N m = 1 M 1 M y n 2 2 α n m y n T s m ( t τ n m ) + α n m 2 2 σ 2 ) f α n m ( α n m )   d α n m ) = 1 ( 2 π σ 2 ) N / 2 m = 1 M n = 1 N K n m
where
K n m = 0 e x p ( 1 M y n 2 2 x y n T s m ( t τ n m ) + α n m 2 2 σ 2 ) f α n m ( α n m )   d α n m = 0 e x p ( 1 M y n 2 2 x y n T s m ( t τ n m ) + x 2 2 σ 2 ) x σ n m 2 e x 2 2 σ n m 2   d x = 1 σ n m 2 e x p ( y n 2 2 M σ 2 ) 0 x e x p ( x y n T s m ( t τ n m ) σ 2 x 2 ( 1 2 σ n m 2 + 1 2 σ 2 ) )   d x
Substituting 1 2 σ n m 2 + 1 2 σ 2 = b 2 and y n T s m ( t τ n m ) σ 2 = c into (25), then, K n m can be rewritten as
K n m = 1 σ n m 2 e x p ( y n 2 2 M σ 2 ) 0 x e x p ( x c x 2 b 2 )   d x = 1 2 b 2 σ n m 2 e x p ( y n 2 2 M σ 2 ) ( 1 + c b e x p ( c 2 4 b 2 ) c 2 b e x p ( p 2 ) d p   )
where p = b x c 2 b . Therefore, the corresponding LRT function can be calculated as:
L ( y n ) = f Y n ( y n | 1 ) f Y n ( y n | 0 ) = m = 1 M n = 1 N 1 1 + σ n m 2 σ 2 ( 1 + π c b e x p ( c 2 4 b 2 ) c 2 b e x p ( p 2 ) d p )
To alleviate the detector calculation complexity, the radar network high-SNR condition, σ n m 2 σ 2 , is typically utilized to simplify the formula. Thus, b 2 can be redefined as:
b 2 1 2 σ n m 2 + 1 2 σ 2 = 1 2 σ 2
Using   c 2 b e x p ( p 2 ) d p 1 1 π ( c / b ) e x p ( c 2 / 4 b 2 ) , σ n m / σ 1 , (27) can be rewritten as
L ( y n ) = m = 1 M n = 1 N 1 1 + σ m n 2 σ 2 ( 1 + π c b e x p ( c 2 4 b 2 )   c 2 b e x p ( p 2 ) d p ) = m = 1 M n = 1 N ( π c b e x p ( c 2 4 b 2 ) 1 + σ n m 2 σ 2 )
Herein, we have tests statistics of the NP detector,
T MIMO = m = 1 M n = 1 N ( π c b e x p ( c 2 4 b 2 ) 1 + σ n m 2 σ 2 ) > 1 < 0 η
where
c b = y n T s m ( t τ n m ) σ 2 1 2 σ 2 = 2 y n T s m ( t τ n m ) σ = 2 E n m c n m σ
Since the test statistic T MIMO . is always positive, we can substitute (30) into (31) and take the logarithm of the two ends to obtain the logarithmic LRT function (log-LRT) shown in (32). Finally, Ln ( T MIMO ) can be compared with the new threshold η which equals Ln ( η ) .
Ln ( T MIMO ) = m = 1 M n = 1 N [ Ln ( 2 π E n m c n m σ e x p ( ( E n m c n m ) 2 2 σ 2 ) ) Ln ( 1 + σ n m 2 σ 2 ) ] > 1 < 0 η
Remark: The simplification of expression (28) can be performed only under high-SRN conditions. As for the case of low-SNR conditions, it will be considered in our future works, which will further enrich the radar system application scenarios in the future.

3. Antenna Placement and Power Allocation Strategy

Aiming at solving the JAPPA problem, we proposed an efficient JAPPA strategy for target detection by using a distributed MIMO radar system. First, the combat area needs to be discretized into grid areas to optimize antenna placement using the numerical method. Then, as shown in the blue dashed box, the antenna location deployment and transmitter power allocation can be optimized by the strategy based on the NP criterion. Finally, the jointly optimal solution can be obtained by a closed-loop iterative approach, shown in red. The whole method structure is illustrated in Figure 2.

3.1. Optimization Formulation

Mathematically, the JAPPA strategy incorporates the antenna placement with power allocation to achieve optimal target detection performance based on the regional characteristics and predetermined power budget. The NP-based detection criterion deduced in Section 3 is appropriate to construct the objective function.
First, we define the following vectors: the transmitter antenna position vector is L T k = [ L T , k 1 , , L T , k m , , L T , k M ] T , L T , k m = ( x T , k m , y T , k m ) . The receiver antenna position vector is L R k = [ L R , k 1 , , L R , k n , , L R , k N ] T , L R , k n = ( x R , k n , y R , k n ) . The transmitter power allocation vector is P T k = [ P T 1 , k , , L T m , k , , L T M , k ] T . Then, the log-LRT can be rewritten as
Ln { T MIMO ( L T k , L R k , P T k ) } = m = 1 M n = 1 N [ Ln ( 2 π E n m c n m σ e x p ( ( E n m c n m ) 2 2 σ 2 ) ) L n ( 1 + σ n m 2 σ 2 ) ] > 1 < 0 η
Considering a realistic confrontation, we set up an area with a priority level using the historical statistic of non-cooperative targets (enemy targets) violations. The position with a higher threat level has higher priority and higher detection probability requirements, and vice versa. In this way, the priority level is utilized as a weighting coefficient multiplying with the NP-based log-LRT at each location. Therefore, it is preferred to use the sum of weighted NP-based log-LRT as the objective function.
F ( L T k , L R k , P T k ) = E { ζ q [ T MIMO ( L T k , L R k , P T k ) ] }
where E { } is expectation operation, T MIMO is the NP-based log-LRT statistic. ζ q is the priority level coefficient, which can be determined by the comprehensive consideration of the military objective location, operational missions, and other parameters. Although the numerical computation of ζ q is beyond the scope of this paper, just for simplicity, we pre-define the priority levels for each regional grid point directly in the simulation.
Many constraints should be taken into account. Firstly, the optimal strategy for the radar network should be subject to acceptable target detection performance.
Ln { T MIMO ( L T k , L R k , P T k ) } > η
Then, the sum of the transmitters’ consumption power is bound by a power budget P T O T A L .
1 M T P T k P T O T A L
where 1 M T = [ 1 , 1 , , 1 ] 1 × M . Next, each transmitter sets the minimum and maximum power values as the constraints [34].
{ P T k m i n P T m , k P T m , k P T k m a x
Therefore, the JAPPA optimization model can be formulated as
P :   max   F ( L T k , L R k , P T k ) s . t . C 1 :   1 M T P T k P T O T A L C 2 :   0 P T k m i n P T m , k , m : { 1 , , M } C 3 :   0 P T m , k P T k m a x , m : { 1 , , M } C 4 :   x T m , y T m N T , m : { 1 , , M } C 5 :   x R n , y R n   N T , n N : { 1 , , N }
where and N represent the sets of antenna labels, and N T represents the antenna candidate points set. Besides, C 1 is the total power consumption constraint, C 2 And C 3 represents the upper and lower power bounds for each transmitter antenna, respectively. C 4 and C 5 ensure the antennas are placed on the coordinate grid point.

3.2. Solution Technique

Due to the presence of the radar transceiver antenna term L T k and L R k in the objective function, (38) is essentially a non-convex non-linear mixed integer programming problem, and problem-solving is non-deterministic polynomial-time hard. A direct method is to partition the discrete coordinate position variables and continuous power variables to search for a solution [22,35].
Here, we propose a novel algorithm based on a two-stage local search, a mixed-integer global optimization algorithm based on greedy take-away and Lagrange. The discrete variable coordinates and the continuous variables are divided using the cutting plane to obtain the integer variable linear approximation of the non-linear programming model. Once the value of the integer variable is determined, the corresponding continuous variable is solved by the Lagrangian function [36]. Finally, integer variables and continuous variables together constitute a feasible solution to the original model. The steps are shown as follows.

3.2.1. Local Search for Antenna Placement

In this case, we define the power of each transmitter as uniformly allocated, i.e., P T , a v e r a g e = P T O T A L / M . The transmitting and receiving antenna deployment at the points from the grid area are known to be NP hard. Simultaneous completion of transmitting and receiving antennas is an arduous computational task; although the exhaustive search method is feasible, the computational effort grows exponentially [27]. Therefore, it is necessary to separate the two antenna arrangements, and the optimization equations for transmitter and receiver antenna placement are, respectively, given below.
T r a n s m i t t e r : P :   max F ( L T , k + 1 , L R , k , P T , a v e r a g e ) s . t . C 1 : x T , k + 1 m , y T , k + 1 m N T , m : { 1 , , M } C 2 : x R , k n , y R , k n N T , n N : { 1 , , N }
R e c e i v e r : P : max F ( L T , k + 1 , L R , k + 1 , P T ,   a v e r a g e   ) s . t . C 1 : x T , k + 1 m , y T , k + 1 m N T , m : { 1 , , M } C 2 : x R , k + 1 n , y R , k + 1 n N T , n N : { 1 , , N }
Initially, a dense uniformed grid with different priority levels over the antenna placement area is assumed, each of the grid points of which serves as a potential antenna placement location. Based on this structure, we propose a site selection algorithm based on a two-stage greedy dropping heuristic to optimize the antenna positions of the transmitter and receiver successively. The purpose is to choose M positions to place the transmitting antennas and then select N positions to place the receiving antennas from the candidate positions to optimize the objective function. Meanwhile, a taboo list has been set up to accelerate the convergence speed. The summary of this algorithm is given in Algorithm 1 (Local search for antenna placement).
Algorithm 1: Two-stage greedy dropping heuristic local search for antenna placement
Input: M ,   N ,   N T ,   ζ q
Output: optimal L T , k = [ L T , k 1 , , L T , k m , , L T , k M ] T and L R , k = [ L R , k 1 , , L R , k n , , L R , k N ] T
1 Initialize k = 1 , q = 1 , L T , 1 1 , , L T , 1 M and L R , 1 1 , , L R , 1 N ;
2 double T e m p 0 , 1 = T e m p 1 , 0 = 0 , ε = 1 e 6 ;
3 while   F ( L T , k , L R , q , P T , a v e r a g e ) T e m p k 1 , q ε  do
4  Replaces T e m p k , q = F ( L T , k , L R , q , P T , a v e r a g e )
5  Recombines L T , k : substitute L T , k m with a new coordinate L ,   L L T , k L T , 0
6  if    F ( L T , k , L R , q , P T , a v e r a g e ) T e m p q , k  then
7    L T , k + 1 = L T , k
8   Removes the L from the set of candidate points
9  end if
10   k = k + 1
11 end while
12 while F ( L T , k , L R , q , P T , a v e r a g e ) T e m p k , q 1 ε  do
13  Replaces T e m p k , q = F ( L T , k , L R , q , P T , a v e r a g e )
14  Recombines L R , q : substitute L R , k n with a new coordinate L ,   L L R , k L R , 0
15  if    F ( L T , k , L R , q , P T , a v e r a g e ) T e m p k , q  then
16    L R , q = L R , q
17   Removes the L from the set of candidate points
18  end if
19   q = q + 1
20 end while
21 return L T , k + 1 , L R , k + 1

3.2.2. Local-Search-Based Lagrange-KKT for Power Allocation

The premise of power allocation is antenna placement; therefore, the optimization model can be given by
P :   max   F ( P T , k + 1 ) C 1 :   1 M T P T , k + 1 = P T O T A L C 2 :   0 P T m i n P T , k m , m : { 1 , , M } C 3 :   0 P T , k m P T m a x , m : { 1 , , M }
where the only variable is the transmitter power vector P T , k + 1 and the problem has been transformed into a univariate optimization model. Thus, the corresponding continuous variable power values can be determined by solving the sub-problems (41). Herein, a Lagrangian multiplier method and the Karush–Kuhn–Tucker (KKT) condition are used to solve the optimal power allocation among transmitters [37], and the Lagrangian function is written in the following
= F ( P T , k + 1 ) λ 1 ( P T m i n P T m , k + 1 ) λ 2 ( P T m , k + 1 P T m a x ) μ ( m = 1 M P T m , k + 1 P T O T A L )
Appendix A proves that (41) is a non-linear convex optimization problem and concludes that the approximate task utility function is unimodal [38]. From this, the optimal power allocation strategy can be obtained, and we have
P T , k + 1 ( P T , k + 1 , λ ¯ , μ ¯ ) = 0
where P T , k + 1 is the local optimum of the problem (41), λ ¯ = [ λ 1 , λ 1 ] and μ ¯ are multiplier vectors; the KKT conditions are formulated as
{ P T , k + 1 = n = 1 N ( 1 2 P T m , k + 1 + L n m c n m 2 2 σ 2 L n m c n m 2 σ 2 + L n m c n m 2 P T m , k + 1 ) + λ 1 λ 2 μ = 0 λ 1 = P T m i n P T m , k + 1 = 0 λ 2 = P T m , k + 1 P T m a x = 0 μ = P T O T A L m = 1 M P T m , k + 1 = 0 λ ¯ 0
Considering σ 2 P T , k + 1 m L m n C m n 2 and all transmitting antennas have to operate, i.e., the inequality constraint cannot fetch an equal sign. Therefore, we simplify the constraint conditions that λ 1 = λ 2 = 0 . The optimal transmitter power allocation strategy can be calculated as follows
{ P T , k + 1 m = N 2 ( μ + n = 1 N L n m c n m 2 2 σ 2 ) m = 1 M P T m , k + 1 = P T O T A L
According to the properties of the convex function, the local optimum P T satisfying the KKT condition is the global optimum of (42), which is the optimal transmitter power allocation scheme. At this time, we have obtained the optimal transmitter power allocation scheme corresponding to the integer variable ( L T , k + 1 , L R , k + 1 ) that maximizes the target detection capability of the MIMO radar network based on the existing antenna deployment.

3.3. Closed-Loop Resource Allocation Scheme

Further on, we exploit the previous results in the last optimization cycle as the feedback information to implement joint antenna placement and power allocation optimization in distributed MIMO radar networks by successive optimization. Overall, we iteratively repeat transmitting and receiving antenna placement and power allocation. In such an operation, first, the updated power allocation scheme is used to obtain the optimal antenna placement scheme. Then, the placement scheme is employed to repeat the power allocation optimization. Finally, the JAPPA strategy incorporates the prior information to achieve the optimal antenna placement and power allocation, which results in superior detection performance in the next round of the period.
In this paper, we combine the idea of cycle minimization (CM) to separate coordinate deployment variables and power variables for cyclic optimization [29,30]. Before the convergence properties analysis, we have
Lemma 1. 
For a given antenna position vector ( L T k , L R k ) , the power allocation problem can be rewritten as a convex optimization problem.
Proof. 
Same as Appendix A. □
Lemma 2. 
Based on the given power allocation scheme P T k , the antenna deployment solution of the local search is the optimal scheme [19].
The closed-loop scheme consists of the two local search algorithms described above and optimizes the ( L T k , L R k ) and P T k in each loop, respectively. In each iteration, the local search algorithm for antenna placement is based on the optimal power solution from the last cycle. Instead, its current placement scheme serves as an input parameter for power allocation, obtaining the initial power value for the next cycle by convex optimization solving.
In summary, the optimization is accomplished by selecting the locations of transmitting and receiving antennas, and then completing the optimal transmit power allocation. On this basis, the joint optimization problem can be solved by using a two-staged algorithm and cyclic optimization, and the cycle terminates to the optimum once satisfying convergence conditions [29,30,39]. This scheme ultimately improves the overall network target detection performance, so that a closed loop is generated. The solution steps summary is given in Figure 3.

3.4. Computational Complexity Analysis

The proposed JAPPA strategy mainly contains two phases that are as follows:
(1)
The local search for antenna placement;
(2)
The Lagrange-KKT for Power Allocation.
Strictly speaking, the computational complexity of the proposed optimization model solving algorithm is mainly determined by the number of iterations of the antenna position selection problem and in the proposed local search method, the exploration for a better antenna placement can be via swapping in and out more points. Therefore, the proposed antenna placement method is with a worst-case computational complexity of O ( ( M + N ) Q 3 ) , wherein Q is the number of selectable coordinate points in the area. Additionally, when the antenna position vectors in (34) are determined, they become convex. At that time, with the Lagrangian-KKT method having a complexity of O ( M N ) in each optimization round [37], the JAPPA method computational complexity can be evaluated as O ( ( M + N ) Q 3 + M N ) . For comparison, the exhaustive search algorithm [23] is with the worst-case complexity of O ( ( C Q M + N / ( M ! N ! ) ) M N ) O ( Q M + N M N ) , which possesses a much heavier computational burden and the optimization method TWANPA, proposed in [19], has a polynomial computational burden of O ( M Q 2 ( N Q + 1 ) + N 2 ) . Therefore, we list the calculation complexity comparison given in Table 1.

4. Simulation Results and Discussion

In this section, the simulation results and numerical analysis are presented to demonstrate the effectiveness of the proposed JAPPA strategy. Herein, in our simulation, a 3 × 3 MIMO radar network ( M = 3 transmitters and N = 3 receivers) is assumed to be placed in a square-shaped plane with sides equal to 20 km. We ignore the electromagnetic coupling between widely separated antennas, and we assume the Swerling I model for the target RCS fluctuation. In addition, to simulate the combat scenario of realistic air defense confrontation, we define a priority level ranging from approximately 0 to 1 for the whole square from the perspective of target detection. As Figure 4 shows, we randomly set an irregular priority level at each grid point, which will be utilized as the weighting coefficient to measure the target detection capability. In addition, the relevant radar parameters used in the simulation are given in Table 2.
The echo SNR is an essential evaluation indicator of a radar system to evaluate its target detection performance. We calculate the average SNR of the radar system by calculating the SNR of each grid point and the echo SNR is given as
S N R = E { y n | 1 } E { y n | 0 }
Moreover, the corresponding detection probability ( P D ) of the radar network at the k-th iteration is defined as [40].
P D k 0.5 × e r f c ( ln P f a S N R k + 0.5 ) P D k 0.5 × e r f c ( ln P f a S N R k + 0.5 )
where complementary error function e r f c is formulated as
e r f c ( z ) = 1 1 π 0 z e v 2 d v
Subsequently, to tackle the initialization of the optimization circulation, we firstly fix the transmitters’ position parameters to optimize the remaining two parameters. Three transmitter antennas are placed at the lowest system radiation power density point one after another [19]. Therefore, the initial transmitting antenna positions are shown in Table 3 and our simulations are run in MATLAB 2019a on a computer with a 2.60 GHz CPU and 16.0 GB RAM. 16.0 GB.
At this point, we utilize the initial transmitting antennas position scheme for the first iterative to start the iteration optimization loop. The positions are chosen one by one according to the weighted NP-based log-LRT and converge to the optimal strategy after six iterations of optimization. The result of receiver and transmitter coordinate positions is shown in Table 4 and Table 5, respectively. Further, the optimal power allocation scheme and antenna geometry are shown in Table 6 and Figure 5, respectively.
The resulting transmitting power allocated to the transmitters over 15 iterations is present in Figure 6. The power allocation values remain stable after the five iterations, which verified the algorithm’s convergence inferred in the previous paper. The transmitting antenna #2 has the highest priority level of the plane, so it receives more power resources than the other two antennas. The relatively large distance between the transmitting antenna #1 and the rest of the radars and its inferior priority position results in a lower allocated power level.

4.1. Factors That Affect the Results of Iterative Optimization

4.1.1. Case 1: Optimization Parameters’ Degree of Freedom

In principle, expanding the freedom of optimizing parameters can effectively improve the operational performance of radar systems. Here, we discuss the effect of optimization parameters’ degree of freedom (DOF). There are three DOFs regarding antenna placement and power allocation: the positions of the transmitting antenna, the positions of the receiving antennas, and the transmitting power. In this case, separate optimization and combinatorial optimization of antenna placements and power allocation are carried out to analyze the different detection performance improvements, respectively. The following groups of different optimization objects will be investigated separately.
  • Single-parameter optimization:
    (1)
    Transmitting antenna position optimization (TAO): The approach simply adjusts the location of the transmitting antenna to enhance the radar network’s target detection capacity.
    (2)
    Receiving antenna position optimization (RAO): The radar network transmitter sites remain constant with evenly distributed transmitting power, and only the placements of the receiving antennas are modified to improve the radar network’s detecting capability.
    (3)
    Transmit power allocation optimization (TPO): The transmitting power is efficiently allocated to increase target detection capabilities based on the random deployment of the radar network antenna space architecture.
  • Two-element optimization:
    (1)
    Joint antenna placement optimization (JAP): In this experiment, the transmitting power is uniformly distributed and the antenna deployment positions are optimized using a two-stage method, in which transmitting and receiving antenna positions are iteratively optimized. The iteration completes when the increment reaches the given error limit.
    (2)
    Joint transmitter placement and power allocation optimization (JTPPA): We jointly design the transmit power allocation and transmitters’ placement to optimize radar performance criteria.
Figure 7 compares the radar target detection performance by the receiver-operating characteristic (ROC) curves of different optimization methods. The larger the area under the curve, the higher the SNR of the radar network, which implies better target detection performance. It can be observed that the optimal transmitting power allocation with random antenna placement disregards the importance of location and the priority level and, thus, TPO provides the worst detection performance. By contrast, the optimization considering the antenna placement with the same power has a slightly better effect. The JAP and JTPPA ROC curves indicate that addressing the power allocation and antenna deployment simultaneously can significantly improve the performance of the radar network. The suboptimal two-element optimization method has an overall higher ability to improve radar performance than the single-parameter optimization method. Furthermore, it can be seen that the proposed JAPPA method optimizes three parameters simultaneously to improve detection ability, thus, providing the most superior detection performance.

4.1.2. Case 2: Iterative Loop Containment Range

This case investigates the effect of iterative loop range on resource allocation, and further analyzes the joint optimization problem of antenna placement and power allocation in distributed MIMO radar networks, focusing on the collaborative optimization of optimizing power allocation and antenna location problem.
Two experiments are carried out: Iterative optimization of antenna placement (IOAP), an algorithm that first iteratively optimizes the antenna position and finally optimizes the power allocation, and iterative loop of antenna position and power distribution (JAPPA). The optimization difference between the two iterative loops is shown in Figure 8a. The first separation of the two dotted lines is due to the JAPPA method performing the optimal power allocation, while IOAP executes the next round of antenna position optimization deployment. Finally, both the curves terminate the optimization process after satisfying the convergence conditions. Figure 8b gives an SNR comparison of the two algorithms, demonstrating that JAPPA can more fully use the limited resources to enhance the system.
The IOAP antenna placement scheme is given in Table 7 and Table 8. Figure 9a illustrates the resulting antenna geometry distribution in the square-shaped plane and Figure 9b presents the solution for the power allocation problem and compares the power allocation schemes of the IOPA and JAAPA methods. The transmitting antenna with the highest power allocation is located at the same location compared to the JAPPA scheme. The essence is that the point has the highest priority, and its results also justify the optimization results.

4.1.3. Case 3: Numbers of Antennas

Obviously, improvements such as raising the transmitting power, improving the antenna gain, and increasing the amount of system equipment can enhance detection performance [9,25]. This simulation is extended to the different numbers of the transmitting and receiving antennas, and the effect of different antenna quantities to complete the target detection task under the constraint of power resource budget has been investigated.
The curves in Figure 10 represent different numbers of transceiver stations, and the detection probability varies with the probability of false alarms. The results indicate that with the increase in the antenna amounts, the radar network possesses a superior detection performance due to the SNR improvement brought by the diversity gain. However, the target detection performance of radar networks under power constraints is not always positively correlated with the number of antennas. Since the radar network detection performance of M = 3 , N = 2 is superior to that of M = 3 , N = 3 , the optimal antennas quantity scheme for monitoring in a particular area is not superior in number. It deserves further research since it has significant guiding relevance for practical applications.

4.2. The Efficiency of the Iterative Two-Stage Closed-Loop Solver

In this section, the time efficiency of the proposed method is illustrated using the exhaustive search-based algorithm as a benchmark. In Figure 11a, it can be observed that the SNR performance of JAPPA is close to that of the exhaustive method, which proves the effectiveness of JAPPA. The threshold ε = 10 6 is set for each iteration increment to fairly compare the computational complexity. Then, the running time of the algorithm is calculated by Monte Carlo experiments, and a visual comparison is shown in Figure 11b. Specifically, it separately takes the exhaustive search method and the TWANPA method about 10 5 and 3 × 10 2   s to complete optimization. By contrast, the proposed two-staged JAPPA algorithm only takes about 2 × 10 2   s in problem solving. That implies the proposed JAPPA method can provide comparable performance to the exhaustive search method while less time-consuming than the TWANPA and conforms to the complex conclusion in Table 1.

4.3. Effectiveness of the Proposed Method

In this part, the effectiveness of the proposed JAPPA technique is analyzed with three comparative tests by evaluating the degree of optimization of the algorithm on the radar target recognition performance.
  • Random antenna placement with uniform power allocation (RAP-UPA): This method randomly selects the positions of the antennas with the uniform transmitting power resource allocation.
  • Random antenna placement withnon-uniform power allocation (RAP-UUPA): The realistic scenario of non-uniform power allocation and random antenna placement will be considered in this simulation.
  • Optimal antenna placement with uniform power allocation (OAP-UPA) [13]: In this scenario, the transmitters and receivers are optimally placed sequentially with the power consumption uniformly allocated to the transmitting antennas.
  • Three-stage water-filling-type antenna placement and power allocation strategy (TWANPA) [19]: This method improves the radar detection performance through successive optimization under the total power constraint and designs a suboptimal method to effectively locate transmitting and receiving antennas under the water-filling power allocation. It completes the transmitting antenna placement using the power density criterion and positioning the receiving antenna with the SNR criterion. Subsequently, a closed-loop optimization mechanism is established by repeating the antennas’ positions using the allocated powers.
The simulation results in Figure 12 present the ROC comparison between the proposed JAPPA strategy and the other four benchmarks. Therefore, it is apparent from a brief look at Figure 12 that the superior performance of the proposed JAPPA strategy is very significant to improve the SNR of the radar system. Among these four benchmarks, OAP-UPA outperforms RAP-UPA and RAP-UUPA due to the optimal antenna positions at both the transmitter and receiver end. TWANPA shows advantages over the OAP-UPA and achieves superior detection performance, since it uses an optimal power allocation scheme based on the water-filling-type method. The TWANPA method selects the optimal antenna placement set based on two different objective functions and iteratively converges to a suboptimal solution, while the proposed JAPPA method completes the optimization based on a unified log-LRT function. The result verifies the effectiveness of the established optimization model.

5. Conclusions

In this paper, the joint resource allocation problem of a distributed MIMO radar network based on environmental information is investigated and an optimization scheme JAPPA is developed to maximize the detection performance of the radar system for the whole surveillance area. The intractable mixed-integer programming problem is modified into a two-stage optimization problem, which is solved by the greedy take-away method and the Lagrangian-KKT method, respectively. Finally, the jointly optimal solution is guaranteed using a closed-loop scheme and the simulation results demonstrate the superiority of the proposed JAPPA method and show the following aspects.
  • The log-LRT function unifies the antenna position variables and transmitter power allocation variables into a single collaborative objective function, which is suitable for guiding the joint resource allocation optimization.
  • The joint antenna and power allocation optimization can more effectively improve the target detection performance of the radar system.
  • The proposed strategy is effective and efficient in solving the JAPPA problem.
Future research could consider more environmental factors, such as terrain, clutter distribution, etc., to get closer to the actual battlefield environment [41,42], and extend to dynamic target detection and joint detection and tracking [43,44]. Beyond that, game theory [45] and machine learning may be taken into account.

Author Contributions

Conceptualization, C.Q.; methodology, C.Q.; validation, C.Q.; formal analysis, C.Q.; investigation, C.Q.; resources, C.Q. and J.X.; data curation, C.Q.; writing—original draft preparation, C.Q.; writing—review and editing, J.X. and H.Z.; visualization, C.Q. and H.Z.; supervision, J.X. and H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 62001506.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Associate Editor and all the reviewers for their valuable and insightful comments and suggestions during the revision of this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Prove that Formula (42) is a convex optimization model.
The constraint in the Formula (42) is linear and convex. The objective function is rewritten as follows:
F ( P T ) = m = 1 M n = 1 N [ Ln ( 2 π E n m c n m σ e x p ( E n m c n m 2 2 σ 2 ) ) L n ( 1 + E n m c n m 2 σ 2 ) ]
To simplify the calculation process, we have
G = Ln ( 2 π c n m L n m P T m σ e x p ( c n m 2 L n m P T m 2 σ 2 ) ) L n ( 1 + c n m 2 L n m P T m σ 2 )
It is necessary to derive the second derivative of G
P m P m G = 2 P T m 2 ( c n m 2 L n m σ 2 ) 2 ( 1 + C n m 2 L n m σ 2 P T m ) 2 2 P T m 2 ( 1 + c n m 2 L n m σ 2 P T m ) 2
where P m P m denotes the second partial derivative with regard to P m . From the literature [36], if the second derivative of a one-variable function is always non-negative in the feasible region, the function is convex. At this point, only the sign of the numerator needs to be considered and it is easy to prove that in the case of the actual radar power consumption order of magnitude, P m P m G is always positive and F ( P T ) can be regarded as the affine transformation of G with respect to P m . Therefore, F ( P T ) is convex. That is, problem (45) is a convex optimization problem [23].
The certificate is completed.

References

  1. Fishler, E.; Haimovich, A.; Blum, R.; Chizhik, D.; Cimini, L.; Valenzuela, R. MIMO radar: An idea whose time has come. In Proceedings of the 2004 IEEE Radar Conference (IEEE Cat. No.04CH37509), Philadelphia, PA, USA, 29 April 2004; pp. 71–78. [Google Scholar]
  2. Yang, J.; Aubry, A.; Maio, A.D.; Yu, X.; Cui, G. Design of Constant Modulus Discrete Phase Radar Waveforms Subject to Multi-Spectral Constraints. IEEE Signal Process. Lett. 2020, 27, 875–879. [Google Scholar] [CrossRef]
  3. Yao, Y.; Liu, H.T.; Miao, P.; Wu, L.N. MIMO Radar Design for Extended Target Detection in a Spectrally Crowded Environment. IEEE Trans. Intell. Transp. Syst. 2021, 1–10. [Google Scholar] [CrossRef]
  4. Stoica, P. MIMO Radar Signal Processing; Wiley-IEEE Press: Hoboken, NJ, USA, 2009. [Google Scholar]
  5. Li, X.; Chen, H.; Zhuang, Z.; Yang, J.; Qiu, Z. Effects of transmitting correlated waveforms for co-located multi-input multi-output radar with target detection and localisation. Signal Process. IET 2013, 7, 897–910. [Google Scholar]
  6. Haimovich, A.M.; Blum, R.S.; Cimini, L.J. MIMO Radar with Widely Separated Antennas. IEEE Signal Process. Mag. 2008, 25, 116–129. [Google Scholar] [CrossRef]
  7. Gogineni, S.; Nehorai, A. Target Estimation Using Sparse Modeling for Distributed MIMO Radar. IEEE Trans. Signal Process. 2011, 59, 5315–5325. [Google Scholar] [CrossRef]
  8. Fishler, E.; Haimovich, A.; Blum, R.S.; Cimini, L.J.; Chizhik, D.; Valenzuela, R.A. Spatial Diversity in Radars—Models and Detection Performance. IEEE Trans. Signal Process. 2006, 54, 823–838. [Google Scholar] [CrossRef]
  9. He, Q.; Blum, R.S.; Haimovich, A.M. Noncoherent MIMO Radar for Location and Velocity Estimation: More Antennas Means Better Performance. IEEE Trans. Signal Process. 2010, 58, 3661–3680. [Google Scholar] [CrossRef]
  10. Godrich, H.; Haimovich, A.M.; Blum, R.S. Cramer Rao bound on target localization estimation in MIMO radar systems. In Proceedings of the 2008 42nd Annual Conference on Information Sciences and Systems, Princeton, NJ, USA, 19–21 March 2008; pp. 134–139. [Google Scholar]
  11. Gacho, J.F.A.L. Researching UAV Threat—New Challenges. In Proceedings of the 2021 International Conference on Military Technologies (ICMT), Brno, Czech Republic, 8–11 June 2021; pp. 1–6. [Google Scholar]
  12. Yao, Y.; Wu, L.; Liu, H.T. Robust Transceiver Design in the Presence of Eclipsing Loss for Spectrally Dense Environments. IEEE Syst. J. 2021, 15, 4334–4345. [Google Scholar] [CrossRef]
  13. Yi, J.; Wan, X.; Leung, H.; Lü, M. Joint Placement of Transmitters and Receivers for Distributed MIMO Radars. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 122–134. [Google Scholar] [CrossRef]
  14. Chitgarha, M.M.; Radmard, M.; Majd, M.N.; Nayebi, M.M. Improving MIMO radar’s performance through receivers’ positioning. IET Signal Process. 2017, 11, 622–630. [Google Scholar] [CrossRef]
  15. Wang, Y.; Yang, S.; Zhou, T.; Li, N. Geometric Optimization of Distributed MIMO Radar Systems with Spatial Distance Constraints. IEEE Access 2020, 8, 199227–199241. [Google Scholar] [CrossRef]
  16. Feng, H.Z.; Liu, H.W.; Yan, J.K.; Dai, F.Z.; Fang, M. A fast efficient power allocation algorithm for target localization in cognitive distributed multiple radar systems. Signal Process. 2016, 127, 100–116. [Google Scholar] [CrossRef]
  17. Yan, J.; Pu, W.; Liu, H.; Jiu, B.; Bao, Z. Robust Chance Constrained Power Allocation Scheme for Multiple Target Localization in Colocated MIMO Radar System. IEEE Trans. Signal Process. 2018, 66, 3946–3957. [Google Scholar] [CrossRef]
  18. Shi, C.G.; Wang, Y.J.; Wang, F.; Salous, S.; Zhou, J.J. Power resource allocation scheme for distributed MIMO dual-function radar-communication system based on low probability of intercept. Digit. Signal Prog. 2020, 106, 12. [Google Scholar] [CrossRef]
  19. Radmard, M.; Chitgarha, M.M.; Majd, M.N.; Nayebi, M.M. Antenna placement and power allocation optimization in MIMO detection. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 1468–1478. [Google Scholar] [CrossRef]
  20. Jebali, S.; Keshavarz, H.; Allahdadi, M. Joint Power allocation and target detection in distributed MIMO radars. IET Radar Sonar Navig. 2021, 15, 1433–1447. [Google Scholar] [CrossRef]
  21. Ma, B.; Chen, H.; Sun, B.; Xiao, H. A Joint Scheme of Antenna Selection and Power Allocation for Localization in MIMO Radar Sensor Networks. IEEE Commun. Lett. 2014, 18, 2225–2228. [Google Scholar] [CrossRef]
  22. Xie, M.C.; Yi, W.; Kirubarajan, T.; Kong, L.J. Joint Node Selection and Power Allocation Strategy for Multitarget Tracking in Decentralized Radar Networks. IEEE Trans. Signal Process. 2018, 66, 729–743. [Google Scholar] [CrossRef]
  23. Godrich, H.; Petropulu, A.P.; Poor, H.V. Sensor Selection in Distributed Multiple-Radar Architectures for Localization: A Knapsack Problem Formulation. IEEE Trans. Signal Process. 2012, 60, 247–260. [Google Scholar] [CrossRef]
  24. Berenguer, I.; Xiaodong, W.; Krishnamurthy, V. Adaptive MIMO antenna selection via discrete stochastic optimization. IEEE Trans. Signal Process. 2005, 53, 4315–4329. [Google Scholar] [CrossRef] [Green Version]
  25. He, Q.; Blum, R.S.; Haimovich, A.M.; He, Z. Antenna placement for velocity estimation using MIMO radar. In Proceedings of the 2008 42nd Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 26–29 October 2008; pp. 619–623. [Google Scholar]
  26. Ajorloo, A.; Amini, A.; Tohidi, E.; Bastani, M.H.; Leus, G. Antenna Placement in a Compressive Sensing-Based Colocated MIMO Radar. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 4606–4614. [Google Scholar] [CrossRef]
  27. Yuan, Y.; Yi, W.; Hoseinnezhad, R.; Varshney, P.K. Robust Power Allocation for Resource-Aware Multi-Target Tracking with Colocated MIMO Radars. IEEE Trans. Signal Process. 2021, 69, 443–458. [Google Scholar] [CrossRef]
  28. Chavali, P.; Nehorai, A. Scheduling and Power Allocation in a Cognitive Radar Network for Multiple-Target Tracking. IEEE Trans. Signal Process. 2012, 60, 715–729. [Google Scholar] [CrossRef]
  29. Zhang, H.W.; Liu, W.J.; Xie, J.W.; Zhang, Z.J.; Lu, W.L. Joint Subarray Selection and Power Allocation for Cognitive Target Tracking in Large-Scale MIMO Radar Networks. IEEE Syst. J. 2020, 14, 2569–2580. [Google Scholar] [CrossRef]
  30. Shia, C.; Ding, L.; Wang, F.; Zhou, J.; Yan, J. Collaborative radar node selection and transmitter resource allocation algorithm for target tracking applications in multiple radar architectures. Digit. Signal Prog. 2022, 121, 103318. [Google Scholar] [CrossRef]
  31. Zhang, H.; Xie, J.; Shi, J.; Zhang, Z.; Fu, X. Sensor Scheduling and Resource Allocation in Distributed MIMO Radar for Joint Target Tracking and Detection. IEEE Access 2019, 7, 62387–62400. [Google Scholar] [CrossRef]
  32. Skolnik, M.I. Radar Handbook; McGraw-Hill Education: New York, NY, USA, 1990. [Google Scholar]
  33. News, S.; News, S. Detection, Estimation and Modulation Theory; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1968. [Google Scholar]
  34. Yuan, Y.; Yi, W.; Kirubarajan, T.; Kong, L.J. Scaled accuracy based power allocation for multi-target tracking with colocated MIMO radars. Signal Process. 2019, 158, 227–240. [Google Scholar] [CrossRef]
  35. Yan, J.K.; Liu, H.W.; Jiu, B.; Chen, B.; Liu, Z.; Bao, Z. Simultaneous Multibeam Resource Allocation Scheme for Multiple Target Tracking. IEEE Trans. Signal Process. 2015, 63, 3110–3122. [Google Scholar] [CrossRef]
  36. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  37. Boyd, S.; Vandenberghe, L.; Faybusovich, L. Convex Optimization. IEEE Trans. Autom. Control. 2006, 51, 1859. [Google Scholar]
  38. Bertsekas, D.P.; Nedic, A.; Ozdaglar, A.E. Convex Analysis and Optimization; Pitman Advanced Pub. Program; MIT: Cambridge, MA, USA, 2003. [Google Scholar]
  39. Stoica, P.; Selen, Y. Cyclic minimizers, majorization techniques, and the expectation-maximization algorithm: A refresher. IEEE Signal Process. Mag. 2004, 21, 112–114. [Google Scholar] [CrossRef]
  40. Mahafza, B.R.; Elsherbeni, A. MATLAB Simulations for Radar Systems Design; Chapman and Hall/CRC: London, UK, 2003. [Google Scholar]
  41. Liang, J.; Zhang, T.; Yang, Y.; Cui, G.; Yang, J. An efficient antenna placement method for MIMO radar under the situation of multiple interference regions. In Proceedings of the 2017 20th International Conference on Information Fusion (Fusion), Xi’an, China, 10–13 July 2017. [Google Scholar]
  42. Zhang, H.; Liu, W.; Zhang, Z.; Lu, W.; Xie, J. Joint Target Assignment and Power Allocation in Multiple Distributed MIMO Radar Networks. IEEE Syst. J. 2020, 15, 694–704. [Google Scholar] [CrossRef]
  43. Ristic, B.; Rosenberg, L.; Du, Y.K.; Wang, X.; Williams, J. A Bernoulli track-before-detect filter for maritime radar. IET Radar Sonar Navig. 2019, 14, 356–363. [Google Scholar] [CrossRef]
  44. Zhang, H.; Liu, W.; Shi, J.; Fei, T.; Zong, B. Joint detection threshold optimization and illumination time allocation strategy for cognitive tracking in a networked radar system. IEEE Trans. Signal Process. 2022; to be published. [Google Scholar]
  45. He, B.; Su, H.; Huang, J. Joint Beamforming and Power Allocation Between A Multistatic MIMO Radar Network and Multiple Targets Using Game Theoretic Analysis. Digit. Signal Prog. 2021, 115, 103085. [Google Scholar] [CrossRef]
Figure 1. System configuration for distributed MIMO radar network.
Figure 1. System configuration for distributed MIMO radar network.
Remotesensing 14 02650 g001
Figure 2. Illustration of the proposed JAPPA method.
Figure 2. Illustration of the proposed JAPPA method.
Remotesensing 14 02650 g002
Figure 3. Block diagram of the closed-loop scheme.
Figure 3. Block diagram of the closed-loop scheme.
Remotesensing 14 02650 g003
Figure 4. Priority level.
Figure 4. Priority level.
Remotesensing 14 02650 g004
Figure 5. Optimal antennas’ geometry distribution.
Figure 5. Optimal antennas’ geometry distribution.
Remotesensing 14 02650 g005
Figure 6. Power allocation results over 15 iterations.
Figure 6. Power allocation results over 15 iterations.
Remotesensing 14 02650 g006
Figure 7. Average ROC curves for power allocation method.
Figure 7. Average ROC curves for power allocation method.
Remotesensing 14 02650 g007
Figure 8. Simulation results comparison (a) optimum iterative process; (b) average ROC curves for IOAP and JAPPA.
Figure 8. Simulation results comparison (a) optimum iterative process; (b) average ROC curves for IOAP and JAPPA.
Remotesensing 14 02650 g008
Figure 9. Optimal scheme based on IOAP. (a) IOAP Antennas’ geometry distribution; (b) IOAP power assigned to the transmitters.
Figure 9. Optimal scheme based on IOAP. (a) IOAP Antennas’ geometry distribution; (b) IOAP power assigned to the transmitters.
Remotesensing 14 02650 g009
Figure 10. Average ROC curves for different antenna amounts.
Figure 10. Average ROC curves for different antenna amounts.
Remotesensing 14 02650 g010
Figure 11. Performance comparison of TWANPA, the exhaustive search, and the JAPPA: (a) average ROC curves; (b) runtime.
Figure 11. Performance comparison of TWANPA, the exhaustive search, and the JAPPA: (a) average ROC curves; (b) runtime.
Remotesensing 14 02650 g011
Figure 12. Average ROC curves for five power allocation methods.
Figure 12. Average ROC curves for five power allocation methods.
Remotesensing 14 02650 g012
Table 1. Computational complexity comparison.
Table 1. Computational complexity comparison.
AlgorithmExhaustive SearchTWANPA MethodJAPPA Strategy
Computational complexity O ( Q M + N M N ) O ( M Q 2 ( N Q + 1 ) + N 2 ) O ( ( M + N ) Q 3 + M N )
Table 2. System parameters.
Table 2. System parameters.
NamesSymbolsSettings
Transmitting Antenna Gain G m 30 dB
Receiving Antenna Gain G n 30 dB
Processing Gain At Receiver I p 20 dB
Carrier Frequency f c 10 GHz
Target Average RCS σ 0 2 2 m2
Radar Network Scattering Loss L c 0 dB
Radar Network Receiving Loss L r 0 dB
Noise Factor F n 4 dB
Bandwidth B 5 MHz
Transmit Power P T O T A L 30 KW
Minimum Transmit Power P T min 0 KW
Maximum Transmit Power P T max 30 KW
Table 3. Initial transmitting antennas positions.
Table 3. Initial transmitting antennas positions.
Transmitterx (km)y (km)
#11410
#2214
#3418
Table 4. Transmitting antennas positions.
Table 4. Transmitting antennas positions.
Transmitterx (km)y (km)
#142
#2106
#31410
Table 5. Receiving antennas positions.
Table 5. Receiving antennas positions.
Transmitterx (km)y (km)
#1168
#21212
#384
Table 6. Receiving antennas positions.
Table 6. Receiving antennas positions.
Transmitter#1#2#3
Power Allocation (KW)0.1335.74424.123
Table 7. Transmitting antennas positions.
Table 7. Transmitting antennas positions.
Transmitterx (km)y (km)
#11814
#2106
#31410
Table 8. Receiving antennas positions.
Table 8. Receiving antennas positions.
Transmitterx (km)y (km)
#11612
#2128
#384
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qi, C.; Xie, J.; Zhang, H. Joint Antenna Placement and Power Allocation for Target Detection in a Distributed MIMO Radar Network. Remote Sens. 2022, 14, 2650. https://doi.org/10.3390/rs14112650

AMA Style

Qi C, Xie J, Zhang H. Joint Antenna Placement and Power Allocation for Target Detection in a Distributed MIMO Radar Network. Remote Sensing. 2022; 14(11):2650. https://doi.org/10.3390/rs14112650

Chicago/Turabian Style

Qi, Cheng, Junwei Xie, and Haowei Zhang. 2022. "Joint Antenna Placement and Power Allocation for Target Detection in a Distributed MIMO Radar Network" Remote Sensing 14, no. 11: 2650. https://doi.org/10.3390/rs14112650

APA Style

Qi, C., Xie, J., & Zhang, H. (2022). Joint Antenna Placement and Power Allocation for Target Detection in a Distributed MIMO Radar Network. Remote Sensing, 14(11), 2650. https://doi.org/10.3390/rs14112650

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