Antenna Allocation in MIMO Radar with Widely Separated Antennas for Multi-Target Detection
Abstract
: In this paper, we explore a new resource called multi-target diversity to optimize the performance of multiple input multiple output (MIMO) radar with widely separated antennas for detecting multiple targets. In particular, we allocate antennas of the MIMO radar to probe different targets simultaneously in a flexible manner based on the performance metric of relative entropy. Two antenna allocation schemes are proposed. In the first scheme, each antenna is allocated to illuminate a proper target over the entire illumination time, so that the detection performance of each target is guaranteed. The problem is formulated as a minimum makespan scheduling problem in the combinatorial optimization framework. Antenna allocation is implemented through a branch-and-bound algorithm and an enhanced factor 2 algorithm. In the second scheme, called antenna-time allocation, each antenna is allocated to illuminate different targets with different illumination time. Both antenna allocation and time allocation are optimized based on illumination probabilities. Over a large range of transmitted power, target fluctuations and target numbers, both of the proposed antenna allocation schemes outperform the scheme without antenna allocation. Moreover, the antenna-time allocation scheme achieves a more robust detection performance than branch-and-bound algorithm and the enhanced factor 2 algorithm when the target number changes.1. Introduction
Multiple input multiple output (MIMO) radar has received considerable attention in recent years. Unlike a standard phased-array radar, where different antennas transmit the scaled and phase-shifted version of the same waveform, MIMO radar can transmit different waveforms at different antennas [1]. In general, MIMO radars can be classified into two categories—MIMO radar with colocated antennas [1] and MIMO radar with widely separated antennas [2]. MIMO radar with colocated antennas improves parameter estimation and beamforming performance by having more effective spatial degrees of freedom [3,4], since its transmit antennas and receive antennas are close enough to observe coherent signals reflected from the target. MIMO radar with widely separated antennas, also referred to as statistical MIMO radar, improves detection and estimation resolution by exploiting the diversity of the propagation path [5,6]. In this paper, we are concerned with a MIMO radar with widely separated antennas, which can be viewed as a type of multistatic radar. Each antenna of the MIMO radar can steer its beam independently towards any direction as an independent transmitter. However, this system differs from the multistatic radar by emphasizing the joint processing of signals for transmission and reception [2].
There are many resource optimization problems we could consider to improve the performance of MIMO radar. Waveform design is one of the most interesting resource optimization problems [7–17]. Optimal illumination cooperation from MIMO radar waveform design can further enhance the capabilities of the radar system. With limited total transmit power, power allocation is regarded as another important radar resource optimization problem [6,18]. In addition, antenna allocation is also analyzed as a resource optimization problem for parameter estimation [19,20]. In order to select an appropriate antenna subset for single-target localization, antenna allocation was formulated as a knapsack problem in [19]. Other formulations of selecting sensors could also be found in [20–22]. However, these researches focus on the optimization of the resource for a single target. In fact, there is another type of resource that we could use to further optimize the performance of MIMO radar: multi-target diversity. We consider a multi-static radar, where its transmit and receive antennas are spaced far away. Furthermore, we assume each transmit antenna is able to steer its beam independently towards different targets and that the receive antennas are able to receive all of the signals reflected from different targets. Therefore, the entire set of antennas form a statistical MIMO radar. We expect that, by steering each antenna to illuminate the “most proper” target, we could improve the overall detection performance. Current literature on multiple targets is not related to the idea of multi-target diversity in the statistical MIMO radar. In [23] and [24], energy allocation for detecting multiple targets in mono-static radar was considered. Additionally, multi-target localization in MIMO radar was investigated in [25,26]. To the best of our knowledge, antenna allocation of a statistical MIMO radar for multi-target detection has not been addressed so far.
Metrics have been introduced to measure the performance of MIMO radar, and these can be used as criteria to optimize the resources. Information theoretic criterion, which was first introduced to radar receiver design in [27], has received considerable attention. Mutual information criterion was proposed for radar waveform design in [8–14]. By maximizing the mutual information between target impulse response and target echoes, radar systems could improve estimation and detection performance. The minimum mean square error (MMSE) criterion was also introduced in [9], which led to the same estimation result as the mutual information criterion. Focusing the transmit power in the target direction was another way to improve estimation accuracy [3,15,16]. Minimizing the probability of miss detection for a given probability of false alarm and maximizing the signal-to-interference were proposed in [11] and [17], respectively. In [14] and [28], another information theoretic criterion, i.e., relative entropy, was used to study the detection performance. It originated from Stein's lemma [29]. Relative entropy is capable of approximately characterizing the probability of a miss in target detection, while mutual information is incapable. Therefore, we will employ relative entropy to analyze radar detection performance in this paper.
Based on relative entropy, the multi-target detection problem can be formulated as a minimum makespan scheduling problem in the combinatorial optimization framework [30]. Minimum makespan scheduling problems are classified as non-deterministic polynomial-time complete ( -complete) in the strong sense [31]. Exhaustive search can be used to solve these problems. This algorithm is simple and always leads to a solution. However, the computational complexity will exponentially grow as the problem size increases. A solution with a given error to the minimum makespan scheduling problem was shown in [30]. Other polynomial-time approximation schemes were also proposed with little penalty compared to an exhaustive search. In [31], the author proposed a polynomial-time approximation algorithm in identical machines, and in [32], the author considered related machines rather than identical machines. However, the “machines” in our problem may be totally different. Therefore, new algorithms will be used in our problem.
The main contributions of this paper are as follows.
Antenna allocation is introduced to exploit the “multi-target diversity” to enhance the detection performance. To implement antenna allocation, a statistical MIMO radar is employed to illuminate multiple targets simultaneously. Two comparative schemes—time allocation and uniform allocation—are also used in the experiments. These two schemes both sequentially illuminate targets one by one.
The antenna allocation problem is formulated as a minimum makespan scheduling problem in a combinatorial optimization framework. In the antenna allocation problem, the detection performance is characterized by relative entropy. The contributions of antennas can be effectively modeled as “processing times” in makespan, and target cells can be effectively modeled as “machines”, which are not identical. A branch-and-bound algorithm is used to achieve the optimal antenna allocation result through a simple transformation of the original problem. Moreover, a new antenna allocation algorithm, called the enhanced factor 2 algorithm, is proposed. This heuristic algorithm employs a greedy strategy to allocate antennas for multiple targets. Both of the branch-and-bound algorithm and the enhanced factor 2 algorithm just allocate each antenna to illuminate a certain fixed target over the entire illumination time. Therefore, the branch-and-bound algorithm and enhanced factor 2 algorithm can be regarded as an antenna-only allocation scheme.
More importantly, we propose an antenna-time allocation scheme as an alternative to the -hard antenna allocation problem. Using the antenna-time allocation scheme, each transmit antenna is selected to illuminate different targets with different illumination time. Therefore, each target has an opportunity to be illuminated by all transmit antennas, so that a new proper antenna allocation result can be explored. In fact, the antenna-time allocation scheme can be regarded as an extension of the antenna-only allocation scheme, since it optimizes both antenna allocation and illumination time allocation.
The rest of this paper is organized as follows. The signal model of MIMO radar and the multi-target detection problem are introduced in Section 2. Antenna allocation schemes are proposed in Section 3, including the antenna-only allocation scheme and the antenna-time allocation scheme. Numerical experiments and analysis are shown in Section 4. Finally, Section 5 concludes the paper.
2. Problem Formulation
2.1. Signal Model
Assume that the entire search area of the space is divided into K target cells, where each cell contains either one or zero targets. Whether there is a target in a certain cell is independent of that in other cells. We consider a MIMO radar consisting of M transmit antennas and N receive antennas. Each transmit antenna is able to steer its beam independently to illuminate any of the target cells. Denote the entire set of transmit antennas by a set = {1,…, M}. Split into K non-overlapping subsets of antennas, and let k ⊆ denote the k-th subset, which illuminates target k. Let Mk denote the cardinality of the set k. We assume that each antenna illuminates one target cell, so that . In realistic scenarios, target number K can be very large in the entire surveillance region. Therefore, we have to detect several times to cover the entire region. For each time of detection, we assume that the number of targets is less than the number of transmit antennas, i.e., M > K. Denote the waveform of the m-th transmit antenna by a vector sm. Thus, the transmitted waveform for target cell k is denoted by a matrix Sk = [sk1, sk2, ⋯, skMk], where the subscript k1 ∈ k denotes the first transmit antenna in k steered towards target k. The transmitted waveform is assumed to be narrowband. Therefore, the waveform that arrives at the n-th receiver can be modeled as [2,7]:
We consider a statistical MIMO radar with homogeneous and sufficiently separated receivers. Thus, the spatial correlations between the columns of Hk can be ignored [7,8]. Likewise, we assume that the noise at different receivers are independent. Specifically, we make the following three assumptions:
Assumption 1 (white noise): The columns of W are identically and independently distributed (i.i.d.) with distribution wn ∼ (0, Rw) [7,8]. Furthermore, since the temporally white noise can be achieved by prefiltering [4,6], Rw can be considered as a diagonal matrix .
The spatial independence between the columns of W is satisfied since the thermal noise at different receivers should be independent [7,8]. In the colored noise case, after a whitening filter, a predistortion will be incorporated for transmitted waveform matrix design, and it can be treated as a matrix transformation. In [8], waveform design in colored noise is considered to optimize the performance of detecting a single target.
Assumption 2 (orthogonal waveforms): The signals transmitted by different antennas are orthogonal [6]:
Waveforms transmitted by different antennas at different frequency bands can approximately satisfy Assumption 2. In [6], it is also shown that a specific set of frequency spread signals can be orthogonal. In fact, to separate the signals from different targets at the receivers, we do not need to assume orthogonality among all signals. We only need to assume that the waveforms from different k are orthogonal. By this, we can apply waveform optimization to each set of k for each target k. For the moment, we make Assumption 2 for simplicity. Extension to waveform optimization for each set of k will be addressed in future work.
Assumption 3 (target scattering matrix): The columns of Hk are i.i.d. with distribution hk,n ∼ (0, RHk). Furthermore, considering sufficiently separated transmitters and receivers, the target scattering coefficients are different and independent. Thus, RHk could be considered as a diagonal matrix [7].
In a statistical MIMO radar, antennas are sufficiently separated. Thus, hm,k,n's can be assumed as independent random variables. Mathematically, the target scattering matrix Hk is a full rank random matrix [2]. The diversity of the propagation paths in a statistical MIMO radar can be characterized by this target scattering matrix in Assumption 3. A specific Hk can be simulated based on the statistic characteristics of Hk. For the moment, we ignore the correlations between the target scattering coefficients for different propagation paths. Extension to a general RHk case may be addressed in future work.
2.2. Multi-Target Detection
Our objective is to employ a MIMO radar to detect multiple targets. To this end, we need to establish a detection model for all targets. Next, we are going to introduce a performance metric on the detection performance of each target.
By Assumption 2, the waveforms transmitted by different antennas are orthogonal to each other. Therefore, each receive antenna can apply sk1 to the received signal given by Equation (3) to separate the signal from transmit antenna k1:
Let collect , , ⋯, , and thus, we have:
Let , and the binary hypotheses, 0 target-absent and 1 target-present, are given by:
It can be verified that Zk is a sufficient statistics for the optimal Neyman-Pearson detection based on the definition of the sufficient statistics. Here, we omit the proof for brevity.
Next, we investigate the detection performance of all targets. In order to measure the detection performance, we introduce a performance metric of relative entropy, i.e., Kullback-Leibler divergence. It originates from Stein's lemma [29] and approximately characterizes the probability of miss detection. Its advantages have received detailed study in [28]. The relative entropy between p0(Zk) and p1(Zk) is defined as:
For a given scenario, some transmitters contribute more to the detection performance than others, since they have lower path losses, better angular views of the target and more advantageous wave lengths and polarizations. If each target cell can be illuminated by the proper transmit antennas, which can maximize the received energy, antenna allocation will achieve much more benefit. Therefore, our objective is to select an appropriate transmit antenna subset for each target cell, i.e., to find an optimal allocation scheme to partition the entire transmit antenna set. In order to guarantee the detection performance of each target cell, our objective function is to maximize the minimum relative entropy of all targets. Thus, the target with the minimum relative entropy is detectable, and other targets are also detectable. Here, we assume that the value of low relative entropy is not too much smaller than the value of high relative entropy. If not, we could waste resources on a very faint target, which remains undetectable. Due to the waste regarding the faint target, the optimization could leave other targets as undetectable. Therefore, in realistic scenarios, we can omit the too low relative entropy from the optimization. Here, we also need to be aware that there are other objective functions that we can use. In the case that the detection performance of all target cells needs to be guaranteed, maximizing the minimum relative entropy is one of the effective choices. Other choices may be studied in future work. If the transmitted energy per transmit antenna is fixed, the optimization problem can be formulated as:
3. Antenna Allocation Schemes
The analysis in this section is to derive a mechanism for the appropriate allocation of transmit antennas.
3.1. Reduced Expression of Relative Entropy
Before we proceed to antenna allocation schemes, we first derive the closed form expression for the relative entropy given in Equation (8) and study its properties. Based on Assumptions 1 and 3, the probability density functions p0(Zk) and p1(Zk) in Equation (8) are written as Equations (10) and (11), respectively.
Substituting Equations (10) and (11) into Equation (8) leads to:
Assuming that the transmit waveform Sk consists of L identical signal pulses, Sk can be expressed as a Kronecker product. We define:
In fact, Equation (14) is simple, since we make those independence assumptions on noise and target scattering coefficients. However, the process of simplification is not the major concern in this paper, and how to simplify a more general model without the assumptions may be addressed in a future study. Next, we can solve the optimization problem (9) via antenna allocation based on Equation (14). Here, we need to be aware that the target scattering coefficients are not known, or we do not need to allocate antennas to improve the detection performance. However, in this problem, the antenna allocation is a dynamic allocation of a set of available antennas. These antennas are allocated to use at each time during a measure period in order to optimize the detection performance. Measure time is partitioned into a sequence of epochs, and one antenna allocation result is employed in one epoch. Our antenna allocation problem refers to the closed-loop solutions to most sensor management problems, i.e., the next allocation is determined while the MIMO radar system is in operation and in view of the results obtained from prior radar measurements in prior epochs [33]. From this point, the antenna allocation optimizes its decision as to how to allocate antennas for the next measurement. We just focus on antenna allocation in one epoch. Particularly, the detection process can be that all transmitters probe each target with some certain number of signal pulses L̂ first. Thus, based on the received signals, we can obtain a coarse in the first epoch. Next, we can allocate antennas more flexibly in each epoch to further make a more accurate decision for the detection based on the previous epoch. Thus, we assume that Equation (14) is obtained in the previous epoch from now on [25,33].
3.2. Antenna-Only Allocation Scheme
In this scheme, we only allocate antennas to improve the detection performance. Considering that Dk(p0‖p1) is a monotonic nondecreasing function of , the optimal solution of (9) will be achieved when , for m ∈ and k ∈ . Substituting into Equation (14) leads to:
Therefore, the problem (9) can be written as:
This type of problem is essentially a minimum makespan scheduling problem in combinatorial optimization algorithms [30]. In [30], the minimum makespan scheduling problem is shown as: given processing times for n jobs, p1,p2, …, pn and an integer m, find an assignment of the jobs to m identical machines, so that the completion time, also referred to as the makespan, is minimized. In the antenna allocation problem, the performance is characterized by the relative entropy. The objective is to maximize the minimum relative entropy. The contributions of antennas are effectively modeled as “processing times” in the makespan scheduling problem, and target cells are effectively modeled as “machines”, which are not identical.
The minimum makespan scheduling problems are classified as -complete in the strong sense [31]. We use a branch-and-bound algorithm and an enhanced factor 2 algorithm to solve this problem, respectively. For a particular scenario, both of these two algorithms only implement antenna allocation over the entire illumination time. Therefore, we put these two algorithms into the same category, called the antenna-only allocation scheme.
3.2.1. Branch-and-Bound Algorithm
This is a classical algorithm for finding optimal solutions of various optimization problems, which was first proposed by Land and Doig in [34]. If computational burden is not considered, a straightforward approach to solve the antenna allocation problem is to enumerate all possible combinations of antennas. Thus, relative entropy function values of all combinations are computed and compared. The optimal solution is the antenna allocation, which leads to the largest function value. However, for problems involving large numbers of antennas, the number of antenna combinations will grow exponentially, which will make the computational burden overwhelming.
The branch and bound algorithm also consists of an enumeration of all candidate antenna allocation solutions. The search process for all antenna combinations is a tree structure, whose nodes are the antenna subsets of all combinations. The step that splits an antenna set into two or more smaller antenna sets is called branching. The step that computes upper and lower bounds for the relative entropy value within a given antenna set is called bounding. A tree node will be discarded, if its upper bound is less than some other node's lower bound. This step is called pruning [34-36]. Using this algorithm, a large number of fruitless candidate solutions can be discarded, which will ease the computational burden to some extent.
To apply the branch-and-bound algorithm in our problem, the problem (18) can be written as the integer programming problem:
Therefore, a branch-and-bound algorithm can be used for this integer programming problem. We need to be aware that branch-and-bound algorithm has exponential worst case complexity, but luckily, it can work with much less effort in some cases. We obtain the optimal antenna allocation result through the branch-and-bound algorithm as a performance bound and compare it with the performance of the proposed enhanced factor 2 algorithm in the experiments.
3.2.2. Enhanced Factor 2 Algorithm
The enhanced factor 2 algorithm is proposed based on the factor 2 algorithm in [30]. The proposed algorithm employs a greedy strategy to arrange the next antenna to illuminate the target cell that has been assigned the least relative entropy so far. The order of arranging antennas is based on the contribution of antennas, i.e., the antenna that can contribute more to the detection performance must have priority. The heuristic scheme is summarized as Algorithm 1.
The relative entropy of each target cell is initialized ℚk = 0, k ∈ . For each Qk,m corresponding to target k and transmit antenna m, calculate as the significance of antenna m to target k. Out of k ∈ and m ∈ , select a pair (k′,m′) with the maximal rk,m, and ℚk′ = ℚk′ + Qk′,m′. Meanwhile, the antenna is discarded from the remaining antenna set, = \m′. Next, we consider the targets, except k′, and thus, we set rk′,m = − 1(∀m ∈ ). Repeat this step till all targets have been selected once. After that, out of k ∈ , select the one with the minimal ℚk. For this cell k″, select a transmit antenna m″ from the remaining antenna set with the maximal Qk″,m, and ℚk″ = ℚk″ + Qk″,m″. Next, transmit antenna m″ is also discarded. Repeat this step till all transmit antennas have been selected, and max min Dk (p0 ‖ p1) = min ℚk.
Though an exhaustive search algorithm can always achieve a solution if it exists, it is implemented with exponential complexity. The branch-and-bound algorithm has exponential worst case complexity. However, the proposed heuristic enhanced factor 2 algorithm offers significantly reduced complexity, and the performance can be seen in the experiments. Here, we need to note that after max min Dk(p0 ‖ p1) is achieved, there can be some antenna m̃ illuminating the target k̃ that is not assigned the least relative entropy. Thus, we could reduce the power of these antennas, i.e., to save resources. This will not influence the antenna allocation result and max min Dk (p0‖p1), but it is useful in realistic applications.
Algorithm 1 Enhanced factor 2 algorithm. | |
1 | Initialize ℚk = 0, k ∈ . |
2 | Calculate , k ∈ , m ∈ . |
3 | for i = 1 : K |
select k′ ∈ and m′∈ s.t. max rk,m; | |
ℚk′ = ℚk′ + Qk′,m′; | |
rk′,m = −1 (∀m ∈ ), = \ m′; | |
end. | |
4 | while ≠ null |
select k″, s.t. min ℚk; | |
select m″ ∈ , s.t. max Qk″,m; | |
ℚk″ = ℚk″ + Qk″,m″, = \ m″; | |
end. | |
5 | max min Dk (p0‖p1) = min ℚk. |
3.3. Antenna-Time Allocation Scheme
Other valid approximation algorithms with high precision for the minimum makespan scheduling problem can be exploited, but we do not focus on these in this paper. We consider that if each target has an opportunity to be illuminated by all transmit antennas, a new appropriate antenna allocation result can be explored. Therefore, we propose an antenna-time allocation scheme to solve the problem from another perspective. Using the antenna-time allocation scheme, antennas are allocated to illuminate different targets with different illumination times.
In order to achieve a proper antenna-time allocation, we propose a method to partition the illumination time based on the illumination probabilities. In particular, the optimal illumination probabilities are obtained first; next, the illumination probabilities can be converted to the illumination time.
Assume that the total illumination time is divided into identical units of time. Denote the illumination from antenna m to target k at the l-th unit time by:
The optimal solution of pk,m(l) can be obtained by solving the following linear optimization problem:
Through linear programming, we can easily obtain the optimal numerical solution of pk,m(l). Observing the problem (22), we can find that pk,m(l) remains constant for l = 1, 2,…, , since Qk,m is a constant for given k and m. In fact, we assume the target features remain the same over the entire illumination time. Therefore, pk,m(l) can be written as:
Lemma 1
The entire illumination process from antenna m to target k is a stochastic process with ergodicity:
Therefore, when the number of unit time is sufficiently large, the time average of the sequence Ik,m(l) is the same as the ensemble average. Denote the illumination time from antenna m to target k by tk,m. According to Lemma 1, tk,m can be obtained:
Next, in order to show an intuitive explanation of the antenna-time allocation scheme, we compare the antenna-time allocation scheme with the antenna-only allocation scheme proposed before. Uniform allocation and time allocation are also used for comparison, both of which employ all antennas to sequentially detect targets one by one. Time allocation optimizes the division of the entire time, while uniform allocation divides the entire time into equal parts for all targets without any optimization.
Figure 1 shows an intuitive explanation of four different illumination schemes in a two-target detection case. Four transmit antennas are utilized in total, and the illumination time is the length of ten signal pulses. Figure 1a shows the illumination via uniform allocation. The uniform allocation is a general method that employs all four antennas to sequentially illuminate targets one by one. It can just be regarded as the conventional single-target detection scheme. The time for Target 1 equals that for Target 2. Figure 1b shows the illumination via time allocation. It is similar to the uniform allocation. However, the illumination time is divided into two proper parts instead of two equal parts. Figure 1c shows the illumination via antenna-only allocation. We divide all antennas into two proper parts to illuminate two targets simultaneously. Antennas 1 and 2 illuminate Target 1; meanwhile, Antenna 3 and Antenna 4 illuminate Target 2. Figure 1d shows the illumination via antenna-time allocation. Each antenna steers its beam independently towards different targets with different times. The illumination time of Antennas 1–4 for Target 1 occupies 80%, 60%, 20% and 10% of the total time, respectively. Therefore, over most illumination time, we still employ Antennas 1 and 2 to illuminate Target 1; meanwhile, we employ Antennas 3 and 4 to illuminate Target 2. However, it is still possible for both Target 1 and Target 2 to be illuminated by all antennas. In fact, via the antenna-time allocation, we optimize both antenna allocation and illumination time allocation in a more flexible manner.
4. Numerical Results and Analysis
In this section, we set scenarios to investigate the detection performance via four different illumination schemes, including a uniform allocation scheme, a time allocation scheme, an antenna-only allocation scheme and an antenna-time allocation scheme. Moreover, the antenna-only allocation scheme contains the branch-and-bound algorithm and the enhanced factor 2 algorithm. Therefore, we use five methods in total to make a comparative analysis. Before simulations, we need to note that in order to allocate transmit antennas for each target cell, we need M > K. However, the number of targets can be very large in the entire surveillance region. Therefore, we have to detect several times to cover the entire region. For each detection, the number of targets to be illuminated simultaneously should not exceed the number of transmitters, i.e., M > K.
4.1. Experiment Results
We set M = 4, N = 4 and target number K = 2. Consider a typical scenario in which RH,1 = diag(5, 0.5, 0.2, 1), RH,2 = diag(0.01, 2, 1, 0.1), which means some antennas are more suitable to illuminate a certain target than the others. Particularly, Target 1 is more sensitive to the illumination from Antennas 1 and 4, while Target 2 is more sensitive to the illumination from Antennas 2 and 3. On the whole, the target scattering intensity of Target 1 is higher than that of Target 2. The different sensitivities can be caused by different path losses, angular views, wave lengths and polarizations. Figure 2 shows the minimum relative entropy of all targets versus 0–25 dB transmitted power via different schemes. It can be observed that the relative entropy of the schemes that we propose outperforms that of uniform allocation and time allocation. It is obvious that uniform allocation is the worst, since it does not use any optimization. When the transmitted power grows, the relative entropy of the antenna-time allocation scheme gradually outperforms that of the antenna-only allocation scheme, i.e., branch-and-bound algorithm and enhanced factor 2 algorithm. In addition, we can find that branch-and-bound algorithm and enhanced factor 2 algorithm have the same relative entropy value for each transmitted power.
In fact, the branch-and-bound algorithm and enhanced factor 2 algorithm lead to the same antenna allocation result for each point on the curves in Figure 2. Figure 3 shows the antenna allocation results via the branch-and-bound algorithm. When the transmitted power ranges from 0 dB to 12 dB, the antenna allocation remains constant, which is shown in Figure 3a. It can be seen that only Antenna 1 is employed to illuminate Target 1. When the transmitted power ranges from 13 dB to 25 dB, the antenna allocation is shown in Figure 3b. It can be seen that both of Antenna 1 and Antenna 4 are employed to illuminate Target 1. The antenna allocation result via the enhanced factor 2 algorithm is exactly the same as that via the branch-and-bound algorithm. Therefore, we can just use Figure 3 to show the antenna allocation results via branch-and-bound algorithm and the enhanced factor 2 algorithm.
Figure 4 shows the number of nodes searched via the branch-and-bound algorithm. It can be observed that the number of nodes is relative large at the transmitted power of 13 dB. Additionally, it has already been shown in Figure 3 that the antenna allocation also changes at this transmitted power. For the computational complexity, this case is worse than the others. Moreover, we can see that in most cases, the branch-and-bound algorithm needs to search relative small number of nodes.
Figure 5 shows the antenna allocation result in the scenario of Figure 2 via the antenna-time allocation scheme. Figure 5a–d shows the antenna allocation results at the transmitted power of 5 dB, 10 dB, 15 dB and 20 dB, respectively. It can be seen that as the transmitted power grows, the illumination time of Antenna 1 for Target 1 increases unceasingly till it occupies the total time. Antenna 4 does not illuminate Target 1 in the low transmitted power region. As the transmitted power grows, Antenna 4 illuminates Target 1 with ever-increasing time.
Figure 6 shows the minimum detection probabilities of targets versus transmitted power via different schemes. The parameters used for simulation in Figure 6 are the same as those in Figure 2. The probability of false alarm is kept constant as Pf = 0.0001. The transmitted power per antenna ranges from 0 dB to 25 dB. In order to observe the detection performance, 104 Monte Carlo trials are executed to achieve the detection probability of each point on the curves. We can find that the branch-and-bound algorithm, the enhanced factor 2 algorithm and the antenna-time allocation scheme almost have the same detection probability for each transmitted power point. All of them can improve detection performance significantly compared with time allocation and uniform allocation.
In order to investigate the influence of different targets on the proposed schemes, we use all illumination schemes to detect different targets in different scenarios. Figure 7 shows the minimum detection probabilities of targets versus -12-12 dB radar cross-section (RCS) fluctuations of targets. The statistics of RCS fluctuations can be given by the log-normal distribution model. The variance is from -12 dB to 12 dB. In Figure 7, 106 Monte Carlo trials are executed to achieve the detection probability for each RCS fluctuation variance. It can be observed that two antenna-only allocation algorithms, i.e., the branch-and-bound algorithm and the enhanced factor 2 algorithm, outperform the time allocation scheme.
Over a large range of RCS fluctuations, the antenna-time allocation scheme is the best, while the uniform allocation is the worst.
Next, we investigate the influence of the target number on the proposed schemes. We set more antennas: M = 9, N = 9. The target number varies from two to nine. Figure 8 shows the minimum detection probabilities of targets versus target number via different schemes. It can be observed that the proposed antenna-only allocation scheme and antenna-time allocation scheme still have advantages over time allocation and uniform allocation with different target numbers. Moreover, the antenna-time allocation scheme achieves the most robust detection performance when the number of targets is larger.
4.2. Analysis
In the experiments above, we have found that the idea of antenna allocation has advantages over the idea without antenna allocation for multi-target detection. In this part, we will show some theoretical analysis.
In the relative entropy Equation (14), we notice that is the ratio of the received signal energy and the noise level at receiver, i.e., SNR. Let denote SNR for one signal pulse; we obtain:
Observing Equation (28), we notice that Dk (p0‖p1) increases linearly with the number of receive antennas N. Therefore, when the other parameters are fixed, the number of receivers should be as large as possible to increase Dk (p0 ‖ p1).
Next, we fix N to study the influence of Mk and Lk. Applying the Taylor expansion of log(·) in the small ρ region, Equation (28) can be approximately written as:
For multi-target detection, K denotes the target number. In the scheme without antenna allocation, i.e., time allocation or uniform allocation, we employ all antennas to illuminate targets sequentially Thus, Mk = M and Lk = L/K, where M and L are the total number of transmit antennas and the total number of signal pulses, respectively Therefore, we have:
In antenna allocation scheme, we partition all antennas to illuminate K targets simultaneously Thus, the average value of Mk is M/K, and Lk = L. We still obtain Equation (31). Therefore, we can see that for each target, the product of the antenna number and the pulse number remains constant, no matter which illumination scheme we use. Thus, Equation (31) can be written as:
Substituting Equation (32) into Equation (30) leads to:
5. Conclusions
In this paper, a multi-target detection problem for MIMO radar with widely separated antennas is investigated based on relative entropy. In order to explore multi-target diversity, we propose the antenna-only allocation scheme and the antenna-time allocation scheme to implement multi-target detection simultaneously. In the antenna-only allocation scheme, we use a branch-and-bound algorithm and an enhanced factor 2 algorithm, respectively. In contrast to sequential illumination without antenna allocation, simultaneous illumination via antenna allocation can increase illumination time at the expense of decreasing the number of illumination antennas for each target. It is shown that proper antenna allocation outperforms time allocation and uniform allocation over a large range of transmitted power, RCS fluctuations and target numbers. Moreover, the antenna-time allocation scheme can achieve a more robust performance than branch-and-bound algorithm and enhanced factor 2 algorithm for more targets.
Acknowledgments
This work was supported by National Natural Science Foundation of China under Grant No. 61371079 and a Postdoctoral Science Foundation funded project.
Author Contributions
Each author contributed extensively to the preparation of this manuscript. Hao Gao conceived the work, designed the algorithms and wrote the manuscript. Jian Wang analyzed the work and designed the framework. Chunxiao Jiang analyzed the experiment results and made suggestions. Xudong Zhang conceived the work, made suggestions and commented on the manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
Appendix
In this Appendix, we show the proof of Lemma 1.
The mean and the variance of can be obtained by Equations (Al) and (A2), respectively:
Based on Equations (A1) and (A2), equals pk,m with probability one.
References
- Li, J.; Stoica, P. MIMO radar with colocated antennas: Review of some recent work. IEEE Signal Process. Mag. 2007, 24, 106–114. [Google Scholar]
- 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]
- Stoica, P.; Li, J.; Xie, Y. On probing signal design for MIMO radar. IEEE Trans. Signal Process 2007, 55, 4151–4161. [Google Scholar]
- Bekkerman, I.; Tabrikian, J. Target detection and localization using MIMO radars and sonars. IEEE Trans. Signal Process. 2006, 54, 3873–3883. [Google Scholar]
- 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]
- 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]
- Song, X.; Willett, P.; Zhou, S.; Luh, P.B. The MIMO radar and jammer games. IEEE Trans. Signal Process. 2012, 60, 687–699. [Google Scholar]
- Tang, B.; Tang, J.; Peng, Y. MIMO radar waveform design in colored noise based on information theory. IEEE Trans. Signal Process. 2010, 58, 4684–4697. [Google Scholar]
- Yang, Y.; Blum, R.S. MIMO radar waveform design based on mutual information and minimum mean-square error estimation. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 330–343. [Google Scholar]
- Yang, Y.; Blum, R.S. Minimax robust MIMO radar waveform design. IEEE J. Sel. Topics Signal Process. 2007, 1, 147–155. [Google Scholar]
- De Maio, A.; Lops, M. Design principles of MIMO radar detectors. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 886–898. [Google Scholar]
- Leshem, A.; Naparstek, O.; Nehorai, A. Information theoretic adaptive radar waveform design for multiple extended targets. IEEE J. Sel. Topics Signal Process. 2007, 1, 42–55. [Google Scholar]
- Bell, M.R. Information theory and radar waveform design. IEEE Trans. Inf. Theory 1993, 39, 1578–1597. [Google Scholar]
- Kay, S. Waveform design for multistatic radar detection. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 1153–1166. [Google Scholar]
- Fuhrmann, D.R.; San Antonio, G. Transmit beamforming for MIMO radar systems using signal cross-correlation. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 171–186. [Google Scholar]
- Li, J.; Xu, L.; Stoica, P.; Forsythe, K.W.; Bliss, D.W. Range compression and waveform optimization for MIMO radar: A Cramér-Rao bound based study. IEEE Trans. Signal Process. 2008, 56, 218–232. [Google Scholar]
- Friedlander, B. Waveform design for MIMO radars. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1227–1238. [Google Scholar]
- Godrich, H.; Petropulu, A.; Poor, H.V. Power allocation strategies for target localization in distributed multiple-radar architectures. IEEE Trans. Signal Process. 2011, 59, 3226–3240. [Google Scholar]
- Godrich, H.; Petropulu, A.; 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]
- Wang, H.; Yao, K.; Pottie, G.; Estrin, D. Entropy-Based Sensor Selection Heuristic for Target Localization. Proceedings of the 3rd International Symposium on Information Processing in Sensor Networks; Berkeley, CA, USA, April 2004; pp. 36–45. [Google Scholar]
- Joshi, S.; Boyd, S. Sensor selection via convex optimization. IEEE Trans. Signal Process. 2009, 57, 451–462. [Google Scholar]
- Gore, D.A.; Paulraj, A.J. MIMO antenna subset selection with space-time coding. IEEE Trans. Signal Process. 2002, 50, 2580–2588. [Google Scholar]
- Fuhrmann, D.R. Active-Testing Surveillance Systems, or, Playing Twenty Questions with a Radar. Proceedings of the 11th Annual Adaptive Sensor and Array Processing (ASAP) Workshop, Lexington, MA, USA; March 2003; pp. 11–13. [Google Scholar]
- Boggio, L.A.; Fuhrmann, D.R. Active-Testing Surveillance for Multiple Target Detection with Composite Hypotheses. Proceedings of 2003 IEEE Workshop on Statistical Signal Processing, St. Louis, MS, USA; September 2003; pp. 641–644. [Google Scholar]
- Garcia, N.; Haimovich, A.M.; Coulon, M.; Lops, M. Resource Allocation in MIMO Radar With Multiple Targets for Non-Coherent Localization. IEEE Trans. Signal Process. 2014, 62, 2656–2666. [Google Scholar]
- Gorji, A.A.; Tharmarasa, R.; Blair, W.D.; Kirubarajan, T. Multiple unresolved target localization and tracking using colocated MIMO radars. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 2498–2517. [Google Scholar]
- Woodward, P.M. Inormation theory and the design of radar receivers. Proc. IRE 1951, 39, 1521–1524. [Google Scholar]
- Tang, J.; Li, N.; Wu, Y.; Peng, Y. On detection performance of MIMO radar: A relative entropy-based study. IEEE Signal Process. Lett. 2009, 16, 184–187. [Google Scholar]
- Cover, T.M.; Thomas, J. Elements of Information Theory; Wiley: New York, NY, USA, 2006. [Google Scholar]
- Vazirani, V.V. Approximation Algorithms; Springer-Verlag: Berlin, Germany, 2001. [Google Scholar]
- Woeginger, G.J. A polynomial-time approximation scheme for maximizing the minimum machine completion time. Oper. Res. Lett. 1997, 20, 149–154. [Google Scholar]
- Azar, Y.; Epstein, L. Approximation schemes for covering and scheduling in related machines. Approx. Algorithms Comb. Optim. Lect. Notes Comput. Sci. 1998, 1444, 39–47. [Google Scholar]
- Hero, A.O.; Cochran, D. Sensor management: Past, present, and future. IEEE Sens. J. 2011, 11, 3064–3075. [Google Scholar]
- Land, A.H.; Doig, A.G. An automatic method of solving discrete programming problems. Econometrica 1960, 28, 497–520. [Google Scholar]
- Zhao, Y.; Chen, J.; Goldsmith, A.; Poor, H.V. Identification of outages in power systems with uncertain states and optimal sensor locations. IEEE J. Sel. Topics Signal Process accepted, to appear. 2014. [Google Scholar]
- Zhao, Y.; Goldsmith, A.; Poor, H.V. On PMU Location Selection for Line Outage Detection in Wide-area Transmission Networks. Proceedings of the IEEE Power and Energy Society General Meeting, San Diego, CA, USA; July 2012; pp. 1–8. [Google Scholar]
© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license ( http://creativecommons.org/licenses/by/3.0/).
Share and Cite
Gao, H.; Wang, J.; Jiang, C.; Zhang, X. Antenna Allocation in MIMO Radar with Widely Separated Antennas for Multi-Target Detection. Sensors 2014, 14, 20165-20187. https://doi.org/10.3390/s141120165
Gao H, Wang J, Jiang C, Zhang X. Antenna Allocation in MIMO Radar with Widely Separated Antennas for Multi-Target Detection. Sensors. 2014; 14(11):20165-20187. https://doi.org/10.3390/s141120165
Chicago/Turabian StyleGao, Hao, Jian Wang, Chunxiao Jiang, and Xudong Zhang. 2014. "Antenna Allocation in MIMO Radar with Widely Separated Antennas for Multi-Target Detection" Sensors 14, no. 11: 20165-20187. https://doi.org/10.3390/s141120165
APA StyleGao, H., Wang, J., Jiang, C., & Zhang, X. (2014). Antenna Allocation in MIMO Radar with Widely Separated Antennas for Multi-Target Detection. Sensors, 14(11), 20165-20187. https://doi.org/10.3390/s141120165