Abstract

Direction-of-arrival (DOA) estimation in multipath environment is an important issue for passive bistatic radar (PBR) using frequency agile phased array VHF radar as illuminator of opportunity. Under such scenario, the main focus of this paper is to cope with the closely spaced uncorrelated and coherent signals in low signal-to-noise ratio and limited snapshots. Making full use of the characteristics of moduli of eigenvalues, the DOAs of the uncorrelated signals are firstly estimated. Afterwards, their contributions are eliminated by means of spatial difference technique. Finally, in order to improve resolution and accuracy DOA estimation of remaining coherent signals while avoiding the cross-terms effect, a new beamforming solution based iterative adaptive approach (IAA) is proposed to deal with a reconstructed covariance matrix. The proposed method combines the advantages of both spatial difference method and the IAA algorithm while avoiding their shortcomings. Simulation results validate its effectiveness; meanwhile, the good performances of the proposed method in terms of resolution probability, detection probability, and estimation accuracy are demonstrated by comparison with the existing methods.

1. Introduction

Passive bistatic radar (PBR) exploiting illuminators of opportunity as uncooperative transmitters has become an emerging technology because it allows target detection and localization with advantages such as low cost, covert detection, and low vulnerability to electronic jamming [1, 2]. However, these researches mainly focus on FM radio broadcasting [3, 4], mobile network base station [5], navigation satellites [6], digital video broadcasting, and digital audio broadcasting [79]. Compared with the aforementioned illuminators of opportunity, a dedicated radar usually has high power and an idea ambiguity function. In recent years, frequency agility and phased array technology are widely used in very high frequency (VHF) radar systems for remote sensing and surveillance purpose because it provides flexible beam scanning, improved detection probability, and high valuable counter to stealth technology [10, 11]. The use of frequency agile phased array VHF radar can extend the range of available illuminators of opportunity, while it also poses many challenges in PBR signal processing, such as low signal-to-noise ratio (SNR), limited snapshots, etc. Thus, some complex signal processing algorithms must be employed, which increase the complexity of the systems. Due to multipath propagation, the direct signal of target and its multipath signal are coherent and closely spaced in the mainlobe. Meanwhile, the transmitting signals emitted from the exploited illuminator are received simultaneously, and they may come from the same direction as the direct signal. These problems significantly degrade the direction-of-arrival (DOA) estimation performance and even result in a failure of target detection and localization.

At present, many high resolution algorithms have been developed to cope with the scenario when uncorrelated and coherent signals coexist. The most popular solutions are the forward/backward spatial smoothing techniques [1215]. The problem, however, is that the decorrelation via spatial smoothing technique is acquired at the cost of the reduction in array aperture, which would further increase the width of mainlobe and make the resolution worse. Besides, these methods cannot resolve the uncorrelated and coherent signals that are in the same direction. Spatial difference methods are proposed to separately estimate the DOAs of the uncorrelated and coherent signals in the multipath environment [1625]. These methods firstly estimate the DOAs of the uncorrelated signals and then eliminate their contribution using spatial difference technique. After that, the DOAs of the remaining coherent signals are estimated. However, the method in [16] is difficult in realization because the covariance matrix of the uncorrelated signals should be known as a priori. The method in [17] suffers from rank deficient problem. The required number of snapshots is too high to put [18, 19] into application. Since [2023] eliminate the contributions of uncorrelated signals by constructing a difference matrix, which encounter the problem of cross-terms effect and may yield wrong DOA estimates when the number of snapshots is small, the method introduced in [24, 25] makes full use of the property of the moduli of eigenvalues to resolve the uncorrelated signals from the coherent signals. However, the main disadvantage is that the DOA estimation performance degrades significantly in the low SNR and limited snapshots. Recently, an iterative adaptive approach (IAA) has attracted much attention [2628]. A number of variant IAA-based beamforming solutions have shown a superresolution performance in dealing with closely spaced signals and work well with few numbers of snapshots [29]. Unfortunately, when uncorrelated and coherent signals are in the same direction, these methods fail to work, either. Combining the above analysis, we can get some conclusions that little literature has been published up to now in this field, partly because of the difficulties in PBR signal processing.

Motivated by previous work, we propose a new DOA estimation method in multipath environment for PBR. By making full use of the property of moduli of eigenvalues, the DOAs of the uncorrelated signals are estimated. Then the contribution of the uncorrelated signals is eliminated by using spatial difference technique; that is, only coherent signals remain in the spatial difference matrix. To eliminate the cross-terms effect while improving the resolution and accuracy of DOA estimation of the remaining coherent signals, a new IAA-based beamforming solution is conducted on a reconstructed covariance matrix. Simulation results validate the better performance of the proposed method by comparison with the existing algorithms.

2. Problems Statement

The basic operating scenario of PBR system is illustrated in Figure 1. As shown in the dashed line box, the PBR system consists of a surveillance antenna, a reference antenna, and a receiver. Meanwhile, a frequency agile phased array VHF radar is used as the uncooperative illuminator of opportunity. Refer to analysis in [10]; the frequency agility technology can increase the ability of antijamming and increase the detection probability. The phased array technology provides flexible beam pattern, angular diversity, and more degrees of freedom, which makes it competent for the task of detecting and tracking different targets.

However, the types of illuminator are enriched, while the difficulties for PBR signal processing also come along. In case of PBR exploiting uncooperative illuminator of opportunity as the transmitter, there is a problem that the instantaneous parameters of transmitting signal are unknown at the receiver; thus, a dedicated reference antenna needs to steer toward the uncooperative illuminator to receive the transmitting signal. Meanwhile, the frequency agility technology destroys the coherency between pulses and the ability of rapidly changing beam scanning makes it impossible for PBR to predict the next beam position. Therefore, these problems make the SNR become low and the number of snapshots is finite.

Furthermore, among the characteristics of the environment, the multipath propagation in VHF band is serious. As shown in the solid line box, the received signal is a mixture of transmitting signal of the exploited illuminator, the direct signal of target, and its multipath signal. The direct signal and multipath signal are coherent and closely spaced in the mainlobe. Meanwhile, the transmitting signal is received simultaneously and may come from the same direction as the direct signal.

Therefore, under such scenarios, the problems of DOA estimation in multipath environment for PBR system are summarized as how to deal with the closely spaced uncorrelated and coherent signals with high resolution and accuracy in the low SNR with limited snapshots.

3. Signal Model

Without loss of generality, in PBR scenario, a uniform linear array (ULA) with isotropic sensors is considered as the surveillance antenna. Let the first sensor be the reference, and take as the DOA of an impinging signal. Then, the steering vector can be expressed aswhere is the wavelength of the signal, is the distance between adjacent sensors, and the superscript denotes the transpose operator.

Assume that far-field narrowband signals are received simultaneously and are mixtures of uncorrelated signals and coherent signals with several groups. The signals are uncorrelated, and the signal that comes from direction corresponding to the signal with power . The remaining coherent signals are divided into groups, which come from statistically independent signals with power , and with multipath signals for each signal. That is, the relationship between uncorrelated and coherent signals is In the coherent group, the multipath signal that comes from the direction corresponds to the multipath propagation of the independent signal and its complex reflection coefficient is . It is assumed that the uncorrelated signals and coherent signals in different groups are regarded as uncorrelated with each other. Consequently, the model of the received signals is given bywhere , , and are the array flow pattern. and contain the reflection coefficient of multipath propagation. ; . , is the source vector of the uncorrelated signals, and is the source vector of the coherent groups. is the additive Gaussian white noise with mean zero and variance , which is uncorrelated with each other and the signals. Besides, the operator represents a block diagonal matrix.

After that, the array covariance matrix can be obtainedwhere is the covariance matrix of uncorrelated signals and is the covariance matrix of coherent signals. stands for the identity matrix. Besides, the superscript denotes the transpose operator; the operator represents a diagonal matrix with the diagonal entry .

4. The Proposed DOA Estimation Method

4.1. DOA Estimation of Uncorrelated Signals

For the DOA estimation of the uncorrelated signals, the eigenvalue decomposition (EVD) of the covariance matrix is conducted and the result is givenwhere , , , and . The is signal subspace matrix whose columns are the eigenvectors corresponds to largest eigenvalues; is noise subspace composed of the eigenvectors correspond to smallest eigenvalues. and are the diagonal matrixes constituted by their eigenvalues.

Moreover, the is also spanned by the [24, 25]; thus it is written aswhere is a full-rank matrix. The can be partitioned into two overlapped submatrixes with the size of , which are expressed aswhere and are the first and last rows of . with . Afterwards, a new matrix is constructed aswhere the superscript denotes the matrix pseudo inverse operator and the matrixes and vectors in (9) have the following forms:andBased on the property of the diagonal elements, which is well documented in [24], and we can obtain the following results:Clearly, by performing EVD of the , the moduli values of the eigenvalues corresponding to the uncorrelated signals are equal to 1, and the moduli values of the eigenvalues corresponding to the coherent signals are all less than 1. Therefore, we can choose a threshold to resolve the uncorrelated signals from the coherent signals. In the meanwhile, the number of uncorrelated signals is estimated by sorting the in descending order and substituting them into the following equation:Thus, if the is the first to satisfy , the number of the uncorrelated signals is , and the corresponding DOA is obtained as follows:

4.2. DOA Estimation of Coherent Signals

The DOAs of the remaining signals are then estimated by using a novel IAA-based beamforming solution. Herein, the equation in (4) is rewritten asin (16), the element in the covariance matrix of uncorrelated signals can be expressed asApparently, is a Toeplitz matrix, according to the property that any Toeplitz matrix will satisfy the fact that where is the exchange matrix, which has unity elements along the counter-diagonal and zeros elsewhere. Thus, the uncorrelated signals can be eliminated by means of spatial difference technique [17], which can be expressed asTheoretically, only coherent signals remain in the spatial differencing matrix in this step. However, in the practical PBR application scenario, the noise cannot be eliminated completely. Furthermore, the EVD of would generate cross-terms effect, which would result in a high probability of failure in DOA estimation of closely spaced coherent signals.

To solve the aforementioned problems, we propose a new IAA-based beamforming solution with respect to an improved covariance matrix. Then, we rewrite (4) aswhere and denote the reconstructed array flow pattern and signal power matrix, which can be expressed asandwhere the diagonal value represents the corresponding signal power at each angle.

Afterwards, according the IAA algorithm [27] (see Appendix), the IAA-based beamformer for each DOA of the remaining coherent signals is computed asTherefore, the estimated signal power of corresponding direction can be obtainedIn order to eliminate the cross-terms effect while improving the resolution and accuracy of DOA estimation, a new covariance matrix is constructed by squaring the spatial difference matrix [21]Then, the spatial spectrum of each direction is obtained by replacing with As a result, by performing the IAA algorithm with respect to the reconstructed covariance matrix, the DOAs of the remaining coherent signals are estimated.

4.3. Summary and Discussion of the Proposed Method

In this subsection, the proposed method is summarized and discussed as follows:

(i) Since the pervious analysis of practical PBR application scenario, the covariance matrix of the received signal is unavailable: thus it is replaced by the sample-average estimated array autocovariance matrix.where is the received signal of the snapshot: is the number of snapshot.

(ii) Then, by performing the EVD of the matrix , the signal subspace is obtained, and the and are computed to construct . Thus, the moduli of eigenvalues of are obtained.

(iii) Based on the property of moduli of eigenvalues, the uncorrelated signals are resolved from coherent signals while the DOAs of the uncorrelated signals are estimated

(iv) To eliminate the contribution of uncorrelated signals and improve the resolution and accuracy DOA estimation of coherent signals, the spatial difference matrix is computed and a new covariance matrix is constructed by squaring the spatial difference matrix.

(v) Finally, the DOAs of remaining coherent signals are estimated by performing IAA algorithm on the reconstructed covariance matrix .

It is noteworthy that our proposed method combines the advantages of both spatial difference method and the IAA algorithm while avoiding their disadvantages. When the uncorrelated and coherent signals coexist in multipath environment, especially the uncorrelated signal coming from the same direction as the coherent signals do, the method can still resolve them, whereas the IAA algorithm fails to work. The performance of the spatial difference method degrades significantly in the low SNR with limited snapshots, but our proposed method can still work well. Besides, by exploiting our method, the DOAs of uncorrelated and coherent signals are estimated separately, which can be parallelized easily to speed up the processing in practical application.

5. Simulation Results

In this section, the performance of the proposed method is investigated by simulation experiments. For PBR system, a ULA of omnidirectional antenna sensors spaced half-wavelength interspacing is used as the surveillance antenna and the half power width of mainlobe can be computed as . In the simulations, the threshold is set as to resolve uncorrelated signals from coherent signals, and the scanning size is uniform in the range from to with increment.

In the first simulation, we consider five uncorrelated signals from and two groups of coherent signals from and . The reflection coefficient of the first group of coherent signals is and the other group is . For simplicity, all signals are of equal power. It is worth pointing out that the two uncorrelated signals are set as and , which come from the same direction with the two groups of coherent signals; in the first group of coherent signals, the two signals are closely spaced in the mainlobe. The simulation is carried out in the SNR of and the number of snapshots of . Therefore, the spatial spectrum of the uncorrelated and coherent signals is shown in Figure 2. As is apparent, the proposed method can easily resolve the uncorrelated signals from the coherent signal in same direction, and it also has satisfactory performance in resolving closely spaced coherent signals.

The second simulation studies the superresolution performance of two closely spaced signals with the proposed method. In the simulation, we define the probability of resolution as where and represent the times of successful resolution result and the total trials, respectively.

The angular gap between two closely spaced coherent signals is defined as , and when the condition of is satisfied, the result is considered successful. Herein, we set the DOA of the first signal as , and the other one vary from to . Under such simulation scenario, the coherent signals are always in the mainlobe. Figure 3 shows the curves of probability of resolution with respect to the SNR of and , the number of snapshots of and . It is shown that our proposed method has a little performance degradation in the low SNR, but the overall performance is satisfactory as the SNR is larger. Especially, when the angular gap is , the probability of resolution of is almost equal to 1 in and . Therefore, the superresolution DOA estimation performance of the proposed method is verified.

In the third simulation, we investigate the high accuracy performance of DOA estimation with our proposed method. The detection probability in the simulation is defined aswhere and stand for the times of successful detection result of signals and the total tests, respectively.

The DOAs detection result is considered successful when the difference between the estimated DOA and the true DOA is less than , i.e., subject to . The simulation considers the same scenario as the first one. The DOA estimation results with respect to the changes of SNR and the number of snapshots are plotted in Figures 4 and 5. Apparently, the curves of the detection probability indicate that our proposed method has good performance of DOA estimation when the SNR varies from and , and the number of snapshots is higher than . This is because that our proposed method makes full use of the true signal power when dealing with uncorrelated and coherent signals, which is different from the existing subspace-based high resolution algorithms. These methods use the orthogonality of the signal subspace and the noise subspace in a certain direction to calculate the spatial spectrum, and the performance of DOA estimation is poor when the SNR is low and the number of snapshots is small. Therefore, similar to the previous results, the proposed method has high accuracy DOA estimation performance in dealing with uncorrelated and coherent signals.

In the last simulation, we further examine the performance of the proposed method by comparison. We select the method in [23] and the method in [24] as the comparative methods. The performance is tested by the root mean square error (RMSE), which is defined asIn each simulation scenario, Monte Carlo runs are performed. Figures 6 and 7 show the RMSE curves of DOA estimation with respect to the SNR and the number of snapshots.

It can be seen from the RMSE curves of DOA estimation with different methods that the DOA estimation performance of our method is better than the method in [23] and the method in [24], especially in the low SNR with small number of snapshots. For method in [23], the reason of performance degradation is that it utilizes the whole array to estimate the uncorrelated signals and then uses the reduced array to estimate the coherent signals in a subsequent stage. When the coherent signals are decorrelated by using forward/backward spatial smoothing technique, the reduction in array aperture would broaden the mainlobe, which significantly degrades the resolution and accuracy of DOA estimation. For method in [24], the performance of DOA estimation of uncorrelated signals is similar to our method due to the same process in dealing with uncorrelated signals. However, when dealing with the remaining coherent signals, it also uses forward/backward spatial smoothing technique, and the performance is not satisfactory, either. For our method, when dealing with uncorrelated and coherent signals, we combine the advantages of spatial difference method and IAA algorithm and propose a novel IAA-based beamforming solution with respect to the improved covariance matrix, which can effectively improve the performance of resolution and accuracy in DOA estimation.

6. Conclusions

In this paper, we propose a novel DOA estimation method in the multipath environment for PBR using frequency agile phased array VHF radar as illuminator of opportunity. The proposed method combines the advantages of spatial difference method and the IAA algorithm while avoiding their shortcomings and then performs a new IAA-based beamforming solution on an improved covariance matrix to deal with closely spaced uncorrelated and coherent signals. The simulation results show that it has satisfactory DOA estimation performance in the scenario of low SNR and small number of snapshots. Furthermore, in the proposed method, the DOAs of uncorrelated and coherent signals are estimated separately, which can be parallelized easily to speed up processing, and can be a viable solution for PBR applications. Since the presented PBR system operates in the VHF band, multipath propagation in complex terrain is serious and unavoidable. In order to further improve the performance of target detection and localization, the dynamic characteristics of electromagnetic and geographical environment and antimultipath technology should be taken into consideration in future works.

Appendix

The IAA algorithm is a data-dependent, nonparametric method which is based on a weighted least squares (WLS) approach. Firstly, we define the covariance matrix of noise based on the above sparse signal modelThen, the WLS cost function is given bywhere . In (A.2), we perform partial derivative operation of with respect to and make the result be zero to obtain that Then we can compute the mean square error (MSE) asSuppose , are a column vector and a row vector, and the inverse matrix of exists. Then, we can get the following matrix inequality:Let ; with , then we can rewrite (A.4) asClearly, we can get the conclusion that when , the result of MSE is minimal.

After that, using the matrix inversion lemma to yieldTherefore, the beamformer for each direction of the corresponding signal isSince the in IAA algorithm depends on the unknown signal power, it must be done in an iterative way. The initialization is implemented by conventional beamforming. It is worth noting that the matrix and are updated from the previous iteration, not the data snapshot.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partly supported by the National Natural Science Foundation of China under Grant 61401489.