Next Article in Journal
Spatiotemporal Variations in Energy Consumption and Their Influencing Factors in China Based on the Integration of the DMSP-OLS and NPP-VIIRS Nighttime Light Datasets
Next Article in Special Issue
Performance Modeling Framework for IoT-over-Satellite Using Shared Radio Spectrum
Previous Article in Journal
Automatic Windthrow Detection Using Very-High-Resolution Satellite Imagery and Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Heterogeneous Satellites for Passive Detection of Moving Aerial Target

1
State Key Laboratory of Integrated Service Networks, Xidian University, Xi’an 710071, China
2
Bradley Department of Electrical and Computer Engineering, Virginia Tech, Blacksburg, VA 24060, USA
3
School of Engineering, University of Warwick, Coventry CV4 7AL, UK
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(7), 1150; https://doi.org/10.3390/rs12071150
Submission received: 20 February 2020 / Revised: 1 April 2020 / Accepted: 2 April 2020 / Published: 3 April 2020
(This article belongs to the Special Issue Satellite Communication)

Abstract

:
Passive detection of a moving aerial target is critical for intelligent surveillance. Its implementation can use signals transmitted from satellites. Nowadays, various types of satellites co-exist which can be used for passive detection. As a result, a satellite signal receiver may receive signals from multiple heterogeneous satellites, causing difficult in echo signal detection. In this paper, a passive moving aerial target detection method leveraging signals from multiple heterogeneous satellites is proposed. In the proposed method, a plurality of direct wave signals is separated in a reference channel first. Then, an adaptive filter with normalized least-mean-square (NLMS) is adopted to suppress direct-path interference (DPI) and multi-path interference (MPI) in a surveillance channel. Next, the maximum values of the cross ambiguity function (CAF) and the fourth order cyclic cumulants cross ambiguity function (FOCCCAF) correspond into each separated direct wave signal and echo signal will be utilized as the detection statistic of each distributed sensor. Finally, final detection probabilities are calculated by decision fusion based on results from distributed sensors. To evaluate the performance of the proposed method, extensive simulation studies are conducted. The corresponding simulation results show that the proposed fusion detection method can significantly improve the reliability of moving aerial target detection using multiple heterogeneous satellites. Moveover, we also show that the proposed detection method is able to significantly improve the detection performance by using multiple collaborative heterogeneous satellites.

Graphical Abstract

1. Introduction

Satellite communications have attracted tremendous attention in both industries and academia because of their significant advantages, such as all-weather operations, wide overage, strong viability, and high reliability [1,2]. To achieve globalization and performance enhancement, satellites have evolved from a single satellite model to a complicated system with multi-satellite concurrent compatibility. Accordingly, in addition to navigation in the traditional application of satellite systems, wireless communications and positioning can also be supported by current satellite systems. Recently, representative satellite systems, such as the global positioning system (GPS), the digital satellite broadcasting system (DVB-S) and the international maritime satellite (INMARSAT) have been widely used in military, industries, and human life [3].
Passive target detection using satellite illumination is an important application in satellite systems, which utilizes signals to detect moving targets through signal processing technologies [4,5]. Most existing works on passive target detection makes a strong assumption of a single satellite, in which only signals for a single satellite need to be processed [6,7,8,9,10]. However, in practice, a receiver in passive target detection is likely to receive illumination signals from different satellites. Thus, strong interference caused by undesired satellite illumination may greatly compromise the performance of traditional passive target detection methods designed for a single satellite illumination [9,11].
In order to adapt to the new multiple heterogeneous satellites environments, in this paper, we focus on developing the effective passive target detection method with multiple satellite illumination. Many technical challenges have to be addressed in this regard. For example, the first obstacle is interference suppression and expected signal purification. Specifically, a near-Earth orbit target may be simultaneously covered by multiple satellite beams and reference signals transmitted in a reference channel may be contaminated by unexpected signals, which makes the estimation of reference signals very difficult [12]. In addition, interference caused by other satellite signals may have negative impact on direct-path interference (DPI)/multi-path interference (MPI) suppression, resulting in poor subsequent echo detection performance in a surveillance channel [13]. Thus, multiple satellite signals received in the reference channel and the DPI/MPI in the surveillance channel should be purified and suppressed, respectively, for effective passive target detection.
Another challenge in passive target detection is that satellite signals received by a satellite signal receiver is generally weak. The huge propagation loss caused by the long-distance propagation from satellites to satellite signal receivers leads to a low received power [14,15]. Furthermore, satellite signals would be further weakened after being absorbed and reflected by a moving target, so-called target echo. Therefore, a long coherent accumulation time need to be taken to detect extremely weak satellite echo. To improve target detection probabilities, the viable option of developing a collaborative detection system using multiple heterogeneous satellite signals will be studied in this paper.
In our earlier work, a passive method leveraging signals from multiple homogeneous satellites was proposed for moving aerial target detection [16]. In this paper, a moving aerial target detection method utilizing multiple collaborative heterogeneous satellites is proposed. the main contributions of the paper are described as follows:
  • For signal purification, parallel band-pass filters are employed to separate the direct wave signals in the reference channel. Moreover, an adaptive filter based on normalized least mean square (NLMS) is applied to suppress direct-path interference (DPI) and multi-path interference (MPI) in the surveillance channel.
  • The maximum value of the cross ambiguity function (CAF) and the fourth order cyclic cumulants cross ambiguity function (FOCCCAF) of each separated direct wave signal and the echo signal are used to improve the signal-to-noise ratio (SNR) of moving aerial target detection.
  • Detection results are obtained by fusing the decisions of multiple sensors to enhance the reliability of moving aerial target detection.
The reminder of this paper is organized as follows. In Section 2, the system model used in this paper is presented. The separation of direct wave in the reference channel and suppression of DPI and MPI in the surveillance channel are described in Section 3. Detection statistics construction and the detector design are proposed in Section 4. In Section 5, the performance of the moving aerial target detection is analyzed. Section 6 shows the numerical results to verify the classification performance. Finally, Section 7 concludes the whole paper.

2. System Model

In this paper, a passive detection system, consisting of one moving aerial target, a satellite signal receiver, and multiple collaborative heterogeneous satellites is considered in the paper, as shown in Figure 1. The satellite signal receiver is comprised of a reference channel antenna and a surveillance channel antenna. The antenna dedicated to the reference channel receives direct wave signals from satellites, including GPS, DVB-S, and INMARSAT satellites. On the other hand, the antenna of the surveillance channel is used to receive target echo signals with the direct wave and multi-path interference.
In Figure 1, the signal x ( t ) in the reference channel is given by [17]
x ( t ) = i = 1 M r i s i ( t ) + n t ,
where M denotes the number of different satellites, n ( t ) is a statistically independent stationary Gaussian white noise with zero mean, s i ( t ) represents different direct wave signal from the i-th satellite, and r i stands for the amplitude of the different direct wave signal in the reference channel.
The signal z ( t ) in the surveillance channel is expressed as
z ( t ) = i = 1 M r i s i ( t D i ) e j 2 π f d i t + i = 1 M Ω i s i ( t ) + i = 1 M j = 1 H ω i j s i ( t τ i j ) + n ( t ) ,
where r i is the amplitude of the i-th echo signal, D i represents the delay of the echo wave relative to the corresponding different direct wave, f d i denotes the doppler shift of the echo of different direct wave signals, Ω i is the amplitude of the direct wave signal in the surveillance channel, ω i j is the amplitude of the direct signal of the i-th satellite passing the j-th multipath channel. τ i j stands for the direct signal of the i-th satellite after the j-th multipath channel. H represents the diameter of the multipath channel in the surveillance channel, n ( t ) is a statistically independent stationary Gaussian white noise with zero mean.

3. Interference Suppression

3.1. Direct Wave Signal Separation in the Reference Channel

As the receiver in the reference channel receives multiple direct wave signals simultaneously, it is necessary to separate direct wave signals in the reference channel to facilitate suppression of DPI and MPI in the surveillance channel and subsequently process. Practically, GPS, DVB-S and INMARSAT satellites use different frequency bands to transmit signals. Therefore, different band-pass filters need to be used to separate various direct wave signals. The band-pass filter is determined by [18]    
H ( e j ω ) = r = 0 Q b r e j ω r 1 + k = 1 S a k e j ω k ,
where a k and b r are the coefficient for the band-pass filter. ω p denotes pass-band upper boundary frequency of the band-pass filter. ω p u represents pass-band lower boundary frequency of the band-pass filter. ω s stands for stop-band upper boundary frequency. ω s u , α p , and α s denote the stop-band lower boundary frequency, the pass-band maximum attenuation and the stop-band minimum attenuation, respectively. The parameters of different filters are set according to different spectral ranges, which determine the values of Q, S, a k , b r of different band-pass filter system functions. Then, a plurality of band-pass filters are connected in parallel to separate direct wave signals at the same time to obtain different direct wave signals. After the separation of direct wave signals, the time-domain signals overlapped in the reference channel have become M independent direct wave signals, where the i-th direct wave signal x i ( t ) can be expressed as
x i t = r i s i ( t ) + n ( t ) .
The time-domain overlapped signals mixed with GPS, DVB-S and INMARSAT satellite signals in the reference channel are used for verification. The spectrogram of mixed direct wave signals is shown in Figure 2, while those of these three signals after separated by the band-pass filter are shown in Figure 3. From Figure 2 and Figure 3, it can be seen that after passing through the band-pass filter, the three satellite signals in time overlapped are effectively separated.

3.2. DPI and MPI Suppression Based on NLMS Adaptive Filter

The echo signals in the surveillance channel are extremely weak, while the power of DPI and MPI are much larger than that of echo signals. The gain obtained by increasing the coherent accumulation time can not satisfy the requirement of echo detection. Therefore, it is necessary to suppress the DPI and MPI in the surveillance channel.
The cascade suppression of DPI and MPI in the surveillance channel is achieved by using an adaptive filter [19,20,21]. The adaptive filter employs NLMS algorithm following [22]
e i ( n ) = z i ( n ) x i T ( n ) w i n ,
w i ( n + 1 ) = w i ( n ) + μ λ + x i ( n ) 2 e i ( n ) x i ( n ) ,
where x i ( n ) denotes the separation direct wave signal in the reference channel, z i ( n ) represents the received signal of the surveillance channel, w i n is the weight vector of the filter, and e i ( n ) is the error signal. In order to avoid a very small x i ( n ) 2 that can cause divergence, λ is set to a small positive constant and μ is the fixed normalized step factor with 0 < μ < 2 .
The surveillance channel may contain a number of different direct wave signals and multi-path interference. Therefore, it is necessary to suppress the DPI and MPI one by one. The process of DPI and MPI suppression is as follows: the direct wave signal of each satellite is taken as an input signal x i ( n ) , firstly. Secondly, the output signal e i ( n ) of the upper stage after the interference suppression by the adaptive filter is used as the training signal z i ( n ) of the adaptive filter of the next stage. This process is repeated M times to complete cascade suppression of the DPI and MPI. Let y ( t ) denote the signal in which the DPI and MPI are suppressed in the surveillance channel as
y ( t ) = z ( t ) i = 1 M j = 0 H w i j . x i ( t j Δ ) ,
where w i j represents the optimal weight of the j-th trajectory of the i-th satellite, and j Δ stands for the delay of the j-th path. Thus, y ( t ) contains only the echo signal of each satellite signal and Gaussian noise in the ideal case.

4. Detection Statistics Construction

4.1. Detection Statistic Construction Based on CAF

After the cascade NLMS adaptive filtering, there are only echo signals and Gaussian white noise in the surveillance channel. In order to detect the echo signals, the direct wave signal x i ( t ) after separation are correlated with the signal y t of the surveillance channel after DPI and MPI suppression to obtain CAF, which is expressed as [23]
χ i ( τ , f d ) = 0 T F x i ( t τ ) y ( t ) e j 2 π f d t d t ,
where τ denotes time delay, f d is the Doppler shift, T F stands for coherent accumulation time. In this paper, several different doppler-time-delay spectra χ i ( τ , f d ) are obtained by Equation (8). Then we can construct the following detection statistic Λ as follows
Λ = χ i ( τ , f d ) .

4.2. Detection Statistic Construction Based on FOCCCAF

After the cascade NLMS adaptive filtering, the separated direct wave signal x i ( t ) and the signal y t after the DPI and MPI suppression in the surveillance channel are processed based on FOCCCAF, which is expressed as
χ y , x i α i f , α i ( u , f ) = + C x i x i x i y α i f ( τ ) C 4 x i α i ( τ u ) e j π f τ d τ ,
where C 4 x i α i ( τ ) represents the fourth-order self-circulating cumulant of the direct wave x i ( t ) , and the cyclic frequency is α i . C x i x i x i y α i f ( τ ) denotes the fourth-order mutual-circulating cumulant of x i ( t ) and y ( t ) , the cyclic frequency is α i f . C 4 x i α i ( τ ) is expressed as
C 4 x i α i ( τ ) = r i 4 M 4 s i α i ( τ ) 3 A r i 2 M 2 s i α i ( τ ) ,
where A is the mean squared value of s i ( t ) . M 4 s i α i ( τ ) and M 2 s i α i ( τ ) are the fourth-order cyclic moment and the second-order cyclic moment of s i ( t ) at α i , which can be expressed as
M 4 s i α i ( τ ) = lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t ) s i ( t ) s i ( t + τ ) ) e j 8 π α i t ,
M 2 s i α i ( τ ) = lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t + τ ) ) e j 4 π α i t ,
where T is the length of time. C x i x i x i y α i f ( τ ) is the fourth-order mutual-circulating cumulant of x i ( t ) and y ( t ) , which can be expressed as [24,25]
C x i x i x i y α i f ( τ ) = r i 3 r i e j π ( α i f + f d i ) D i M 4 s i α i f + f d i ( τ D i ) + 3 r i 2 r i e j π f d i τ e j π ( α i f + f d i ) D i M s i s i n s i α i f + f d i ( τ D i ) 3 B r i r i e j π f d i τ e j π ( α i f + f d i ) D i R s i α i f + f d i ( τ D i ) + r i 3 M s i s i n s i α i f ( τ ) ,
where B = E [ x i ( t ) y ( t ) ] , α i f denotes the fourth-order mutual-cyclic frequency of x i ( t ) and y ( t ) . R s i α i f + f d i ( τ D i ) is the cyclic autocorrelation of the i t h signal s i ( t ) , which is given by
R s i α i f + f d i ( τ D i ) = lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t + τ ) ) e j 2 π α i t .
M 4 s i α i f + f d i ( τ D i ) is the fourth-order mutual-circulating moment of s i ( t ) and n ( t ) in the reference channel, which can be expressed as
M s i s i n s i α i f + f d i ( τ D i ) = lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t ) n ( t ) s i ( t + τ D i ) ) e j 2 π ( α i f + f d i ) t .
M s i s i n s i α i f ( τ ) is the fourth-order mutual-circulating moment of s i ( t ) and n ( t ) , which can be expressed as
M s i s i n s i α i f + f d i ( τ D i ) = lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t ) n ( t ) s i ( t + τ D i ) ) e j 2 π ( α i f + f d i ) t .
Substituting Equation (11) and Equation (14) into Equation (10), we can obtain
χ y , x i α i f , α i ( u , f ) = + C x i x i x i y α i f ( C 4 x i α i ( τ u ) ) e j π f τ d τ = r i 7 r i e j π ( α i f + f d i ) D + M 4 s i α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ + 9 A B r i 3 r i e j π ( α i f + f d i ) D + R s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) ) e j π ( f f d i ) τ d τ 3 B r i 5 r i e j π ( α i f + f d i ) D + ( M 4 s i α i ( τ u ) ) R s i α i f + f d i ( τ D i ) e j π ( f f d i ) τ d τ 3 A r i 5 r i e j π ( α i f + f d i ) D + M 4 s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ 9 A r i 4 r i e j π ( α i f + f d i ) D + M s i s i n s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ + 3 r i 6 r i e j π ( α i f + f d i ) D + M s i s i n s i α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ 3 A r i 5 + M s i s i s i n α i f ( τ ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) e j π f τ d τ + r i 7 + M s i s i s i n α i f ( τ ) ( M 4 s i α i ( τ u ) ) e j π f τ d τ .
In Equation (18), the first four items do not contain any noise, which is used to detect the echo signal, and constructs the following detection statistic Λ as
Λ = r i 7 r i e j π ( α i f + f d i ) D + M 4 s i α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ + 9 A B r i 3 r i e j π ( α i f + f d i ) D + R s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) ) e j π ( f f d i ) τ d τ 3 B r i 5 r i e j π ( α i f + f d i ) D + ( M 4 s i α i ( τ u ) ) R s i α i f + f d i ( τ D i ) e j π ( f f d i ) τ d τ 3 A r i 5 r i e j π ( α i f + f d i ) D + M 4 s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ .
In Equation (19), let ρ M R denote the correlation coefficient of the fourth-order self-circulating moment M 4 s i α i ( τ ) of the signal s i ( t ) and the cyclic autocorrelation function R s i α i ( τ ) of s i ( t ) as
ρ M R = cov ( M 4 s i α i ( τ ) , R s i α i ( τ ) ) var ( M 4 s i α i ( τ ) ) var ( R s i α i ( τ ) ) ,
where var ( M 4 s i α i ( τ ) ) represents the variance of M 4 s i α i ( τ ) , var ( R s i α i ( τ ) ) denotes the variance of R s i α i ( τ ) , and cov ( M 4 s i α i ( τ ) , R s i α i ( τ ) ) is the covariance of M 4 s i α i ( τ ) and R s i α i ( τ ) , which can be expressed as
cov ( M 4 s i α i ( τ ) , R s i α i ( τ ) ) = E [ M 4 s i α i ( τ ) R s i α i ( τ ) ] E [ M 4 s i α i ( τ ) ] E [ R s i α i ( τ ) ] .
From Equation (12) and Equation (15), we can obtain
E [ M 4 s i α i ( τ ) R s i α i ( τ ) ] = E [ lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t ) s i ( t ) s i ( t + τ ) ) e j 8 π α i t 1 T t = 0 T 1 ( s i ( t ) s i ( t + τ ) ) e j 2 π α i t ] = E [ lim T 1 T t = 0 T 1 ( s i ( t ) s i ( t ) s i ( t ) s i ( t ) ) e j 8 π α i t lim T 1 T t = 0 T 1 ( s i ( t + τ ) s i ( t + τ ) ) e j 2 π α i t ] = E [ M 4 s i α i ( τ ) ] E [ R s i α i ( τ ) ] .
From the above analysis, we can obtain
cov ( M 4 s i α i ( τ ) , R s i α i ( τ ) ) = E [ M 4 s i α i ( τ ) R s i α i ( τ ) ] E [ M 4 s i α i ( τ ) ] E [ R s i α i ( τ ) ] = 0 .
Using Equation (23) in Equation (20), we can also obtain
ρ M R = 0 ,
and Equation (19) can be simplified as
Λ = r i 7 r i e j π ( α i f + f d i ) D + M 4 s i α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ + 9 A B r i 3 r i e j π ( α i f + f d i ) D + R s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) ) e j π ( f f d i ) τ d τ .

4.3. Detector Design

In order to improve the reliability of the final detection results, distributed sensors are used for collaborative detection. Each sensor is designed as a binary detector, whose binary hypothesis can be described as: assumption H 1 shows that the surveillance channel signal y ( t ) contains the echo signal r i s i ( t D i ) e j 2 π f d i t corresponding to the i-th direct wave signal r i s i ( t ) ; assumption H 0 shows that the surveillance channel y ( t ) does not contain the echo signal r i s i ( t D i ) e j 2 π f d i t corresponding to the i-th direct wave r i s i ( t ) ), which can be expressed as
H 1 : y ( t ) = η = 1 η = i M r η s η ( t D η ) e j 2 π f d η t + n ( t ) , H 0 : y ( t ) = η = 1 η i M r η s η ( t D η ) e j 2 π f d η t + n ( t ) .
Based on the binary hypothesis in Equation (26), the binary detector is
Λ > H 1 < H 0 T d i ,
where T d i is the detection threshold. Each sensor uses detector in Equation (27) to make the decision. Each sensor decision u i is
u i = 0 e l s e , 1 Λ T d i .
When output u i of the sensor is one, target is detected, and H 1 is true. When output u i of the sensor is zero, no target is detected, and H 0 is true. The performance of the i t h sensor can be denoted as P ( u i / H j ) , j = 0 , 1 . The false alarm probability of the i t h sensor can be denoted as P i ( 1 / H 0 ) = P f i , the detection probability of the sensor can be denoted as P i ( 1 / H 1 ) = P d i . According to Equation (28), the fusion center of the distributed sensor fuses the decision result of each sensor according to the criterion of "h in k". It means that if there are more than or equal k sensors that detect the target in h sensors, the overall decision is H 1 . The decision of each sensor is independent of each other, the fusion center gives the decision output u as
u = 1 i = 1 h u i k , 0 e l s e .
After the fusion decision of the fusion center, the final detection probability P d is expressed as
P d = P ( u = 1 / H 1 ) = P ( u = 1 / u 1 , u 2 , u h , H 1 ) P ( u 1 , u 2 , u h / H 1 ) = P ( u = 1 / u 1 , u 2 , u h , H 1 ) i = 1 h P ( u i / H 1 ) = P ( u = 1 / u 1 , u 2 , u h , H 1 ) i S 1 P ( u i = 1 / H 1 ) i S 0 P ( u i = 0 / H 1 ) = P ( u = 1 / u 1 , u 2 , u h , H 1 ) i S 1 P d i i S 0 ( 1 P d i ) = i = k h C i / h { P d 1 , P d 2 , P d h } = i = k h u i = i i P d i u i ( 1 P d i ) 1 u i ,
where S 1 = i u i = 1 , i = 1 , , h , S 0 = i u i = 0 , i = 1 , , h . The final false alarm probability P f is expressed as
P f = P ( u = 1 / H 0 ) = P ( u = 1 / u 1 , u 2 , u h , H 0 ) i S 1 P f i i S 0 ( 1 P f i ) = i = k h C i / h P f 1 , P f 2 , , P f h = i = k h u i = i i P f i u i ( 1 P f i ) 1 u i .
From the above, the procedure of the moving aerial target detection method based on multiple collaborative heterogeneous satellites is summarized in Algorithm 1.
Algorithm 1 The procedure of the proposed passive detection method.
1:
Different direct wave signals in the reference channel are parallelly separated by Equation (3);  
2:
DPI and MPI are in the surveillance channel cascadely suppressed by Equation (5) and Equation (6);  
3:
FOCCCAF and CAF are introduced by Equation (10) and Equation (8) to construct the detection statistic Λ ;  
4:
Fuse the binary decision results of each sensor among the distributed sensors to get the final target detection results by Equation (29).

5. Moving Aerial Target Detection Performance Analysis

5.1. Performance Analysis of Detection Based on CAF

Suppose that the real and imaginary parts of Gaussian noise n ( t ) and Gaussian noise n ( t ) are distributed as N ( 0 , σ n 2 ) and N ( 0 , σ n 2 ) , respectively. The separated reference signal and the echo signal from which the direct wave signal and multipath interference are suppressed in the surveillance channel are calculated CAF and discrete, which can be expressed as
χ i ( l T s , f d ) = m = 0 N 1 x i ( m T s l T s ) y ( m T s ) e j 2 π f d m T s ,
where T s represents the sample interval, N is the number of sampling points and T s stands for the sampling time, N = T F T F T s T s . Equation (32) can also be expressed as
χ i ( l , f d ) = m = 0 N 1 x i ( m l ) y ( m ) e j 2 π f d m = m = 0 N 1 Re { y ( m ) x i ( m l ) e j 2 π f d m } + j m = 0 N 1 Im { y ( m ) x i ( m l ) e j 2 π f d m } .
Lemma 1.
The distribution of the detection based on CAF under the H 1 hypothesis is
( Λ / H 1 ) C N ( r i r i χ s i s i ( l ζ , f d f d i ) , 2 N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) ) ,
and the distribution of the detection based on CAF under the H 0 hypothesis is
( Λ / H 0 ) C N ( 0 , 2 N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 ) ) ,
where C N ( . ) represents the complex Gaussian process.
Proof. 
See Appendix A.1. □
Form the above analysis, the detection probability of each sensor is
p d i = T d i + ρ ( Λ / H 1 ) d Λ = Q 1 N r i r i 2 N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) , T d i N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) ,
where ρ ( Λ / H 1 ) represents the probability density function of ( Λ / H 1 ) , Q 1 ( · , · ) denotes the Marcum Q function [26,27,28]. The false alarm probability for each sensor is
P f i = T d i + ρ ( Λ / H 0 ) d Λ = exp T d i 2 N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 ,
where ρ ( Λ / H 0 ) represents the probability density function of ( Λ / H 0 ) . Substituting Equation (36) and Equation (37) into Equation (30) and Equation (31), the overall detection probability and false alarm probability after decision fusion are
P d = i = k h u i = i i Q 1 N r i r i 2 N ( r i 2 σ n 2 + η = 1 M ( r η ) 2 σ n 2 ) , T d i N ( r i 2 σ n 2 + η = 1 M ( r η ) 2 σ n 2 ) u i · 1 Q 1 N r i r i 2 N ( r i 2 σ n 2 + η = 1 M ( r η ) 2 σ n 2 ) , T d i N ( r i 2 σ n 2 + η = 1 M ( r η ) 2 σ n 2 ) 1 u i ,
and
p f = i = k h u i = i i exp T d i 2 N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 u i 1 exp T d i 2 N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 1 u i .
The echo signal-to-noise ratio S N R η of the surveillance channel in the η -th detection channel is defined as
S N R η = 10 lg ( r η ) 2 σ n 2 .
From Equation (40), we can obtain
10 S N R η 10 = ( r η ) 2 σ n 2 .
In this case when η is equal i, by substituting Equation (41) into Equation (36), we can get the relationship between the detection probability p d i of a single sensor and the echo signal-to-noise ratio S N R i as
p d i ( S N R i ) = Q 1 N r i r i N ( r i 2 ( r i ) 2 10 S N R i 10 + η = 1 M ( r η ) 2 σ n 2 ) , T d i N ( r i 2 ( r i ) 2 10 S N R s i 10 + η = 1 M ( r η ) 2 σ n 2 ) ,
and the relationship between the final detection probability p d after the fusion and the echo signal-to-noise ratio S N R is
p d = i = k M u i = i i = 1 Q 1 N r i r i 2 N ( r i 2 r i 2 10 S N R i 10 + η = 1 M r η 2 σ n 2 ) , T d i N ( r i 2 r i 2 10 S N R i 10 + η = 1 M r η 2 σ n 2 ) u i · 1 Q 1 N r i r i 2 N ( r i 2 r i 2 10 S N R i 10 + η = 1 M r η 2 σ n 2 ) , T d i N ( r i 2 r i 2 10 S N R i 10 + η = 1 M r η 2 σ n 2 ) 1 u i .
From Equation (43), it can be seen that the detection probability p d is related to the number of sampling points N, the power ratio of the direct echo r i 2 r i 2 r i 2 r i 2 , the detection threshold T d i and the noise power σ n 2 in the reference channel with different echo signal-to-noise S N R i .

5.2. Performance Analysis of Detection Based on FOCCCAF

It is assumed that the noises n ( t ) and n ( t ) obey the Gaussian distribution of N ( 0 , σ n 2 ) and N ( 0 , σ n 2 ) , respectively. We can get the FOCCCAF from the separated reference signal and the echo signal after the DPI and MPI suppression in the supervise channel by using Equation (10).
Lemma 2.
The distribution of the detection based on FOCCCAF under the H 1 hypothesis is
Λ / H 1 C N ( r i 7 r i e j π ( α i f f d i ) D i χ 4 s i α i f f d i , α i ( u , f ) + 9 A B r i 3 r i e j π ( α i f f d i ) D i χ 2 s i α i f f d i , α i ( u , f ) , σ s i 10 N 3 ( 81 A 2 r i 8 ( r i ) 2 σ n 2 + 9 r i 12 ( r i ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) ) .
and the distribution of the detection based on FOCCCAF under the H 0 hypothesis is
Λ / H 0 C N ( 0 , σ s i 10 N 3 ( η = 1 η i M 81 A 2 r i 8 ( r η ) 2 σ n 2 + η = 1 η i M 9 r i 12 ( r η ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) ) .
Proof. 
See Appendix A.2. □
From the above analysis, the detection probability for each sensor is
p d i = T d i + ρ ( Λ / H 1 ) d Λ = Q 1 N 2 r i 7 r i + 9 A B r i 3 r i 2 Y , T d i Y ,
where ρ ( Λ / H 1 ) represents the probability density function of ( Λ / H 1 ) and Y = σ s i 10 N 3 ( 81 A 2 r i 8 ( r i ) 2 σ n 2 + 9 r i 12 ( r i ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) . The false alarm probability for each sensor is
P f i = T d i + ρ ( Λ / H 0 ) d Λ = exp T d i σ s i 10 N 3 ( η = 1 η i M 81 A 2 r i 8 ( r η ) 2 σ n 2 + η = 1 η i M 9 r i 12 ( r η ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) ,
where ρ ( Λ / H 0 ) represents the probability density function of ( Λ / H 0 ) . Substituting Equation (46) and Equation (47) into Equation (30) and Equation (31), the detection probability and false alarm probability after distributed sensor fusion are
P d = i = k h u i = i i Q 1 N 2 r i 7 r i + 9 A B r i 3 r i 2 Y , T d i Y u i · 1 Q 1 N 2 r i 7 r i + 9 A B r i 3 r i 2 Y , T d i Y 1 u i ,
and
p f = i = k h u i = i i exp ( T d i Δ ) u i ( 1 exp ( T d i Δ ) ) 1 u i ,
where Δ = σ s i 10 N 3 ( η = 1 η i M 81 A 2 r i 8 ( r η ) 2 σ n 2 + η = 1 η i M 9 r i 12 ( r η ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) . In the supervise channel, the i t h ( i = η ) echo signal-to-noise ratio S N R i is defined as
S N R i = 10 lg ( r i ) 2 σ n 2 .
From Equation (50), we can obtain
10 S N R i 10 = ( r i ) 2 σ n 2 .
Substituting Equation (51) into Equation (46), the relationship between the detection probability of a single sensor P d i and the echo signal-to-noise ratio S N R i is expressed as follows
p d i = T d i + ρ ( Λ / H 1 ) d Λ = Q 1 N 2 r i 7 r i + 9 A B r i 3 r i 2 Ξ , T d i Ξ ,
where Ξ = σ s i 10 N 3 ( 81 A 2 r i 8 ( r i ) 2 σ n 2 + 9 r i 12 ( r i ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 9 A 2 r i 10 ( r i ) 2 10 S N R i 10 + r i 14 ( r i ) 2 10 S N R i 10 σ s i 4 ) . The relationship between the final detection probability after fusion P d and the echo signal-to-noise ratio SNR is expressed as follows:
P d = i = k h u i = i i Q 1 N 2 r i 7 r i + 9 A B r i 3 r i 2 Ξ , T d i Ξ u i · 1 Q 1 N 2 r i 7 r i + 9 A B r i 3 r i 2 Ξ , T d i Ξ 1 u i .
From Equation (53), it can be seen that the detection probability P d is related to the number of sampling points N and the power ratio of the direct echo r i 2 r i 2 r i 2 r i 2 , the power of s i ( t ) and the noise power in the reference channel with different echo signal-to-noise.

6. Numerical Results and Discussion

To verify the effectiveness of the proposed passive detection method and investigate and the influence of different parameters on the detection performance, a series of simulation experiments are conducted using MATLAB (9.5.0.944444 (R2018b), MathWorks Company, Natick, MA, USA). Simulation parameters are set according to [29,30]. GPS satellite, DVB-S satellite and INMARSAT satellite signals are considered in the simulation with Gaussian white noise. The ratio of echo signal to direct wave power is defined as
D S R = P D P S ,
where N represents the number of sampling points, P S and P D are the power of the direct wave and the echo, respectively. The detection accuracy rate is used to evaluate the detection algorithm as
δ H = N R N × 100 % ,
where δ H and N R are the number of corrected detection and the total number of detection in the simulation, respectively.
In order to evaluate the influence of sampling points on the detection performance of weak echo signals, three satellite signals GPS, DVB-S and INMARSAT are used for simulation experiments. The carrier frequencies of the three signals are: 1.57 GHz, 12.38 GHz, 4.2 GHz; the chip-rates are: 1.023 MHz, 22.425 MHz, 2.2 MHz; the delay is 1 us, 2 us, 3 us; the doppler shift are: 100 Hz, 200 Hz, 300 Hz; the intensity of the three signal direct waves are: −130.1 dBw, −111.83 dBw, −120.61 dBw, and the direct wave and its corresponding echo power difference are 40 dBw, P f a = 10 4 , coherent accumulation points are: 10 7 , 10 8 , 10 9 .2000 Monte Carlo simulation experiments are carried out and the simulation results are shown in Figure 4.
From Figure 4, it can be seen that the detection performance of the weak echo gradually increases with the increase of the accumulated points under the same SNR. This is because, as the number of accumulated points increases, the peaks at the corresponding delay and Doppler frequency shift become more prominent, which is equivalent to increasing the echo SNR gain, thereby improving detection performance. When the cumulative number of points is 10 9 and the SNR is −75 dB, the detection probability of the echo signal reaches 99%. For the same cumulative number of points, as the SNR becomes larger, the detection performance improves. The simulation results agree with the theoretical results, and the feasibility of detecting weak echo signals based on FOCCCAF is verified.
In order to compare the proposed method based on FOCCCAF with the method based on CAF, three satellite signals GPS, DVB-S and INMARSAT satellite signals are used for simulation experiments. P f a = 10 4 , coherent accumulation point is 10 9 , 2000 Monte Carlo simulation experiments are carried out and the simulation results are shown in Figure 5.
It can be seen from Figure 5 that the detection performance based on FOCCCAF is better than that on CAF, and the detection performance is improved by about 6 dB. This is because both FOCCCAF and CAF perform echo detection by detecting the delay of the corresponding echo signal and the peak at the Doppler shift, and the fourth-order cyclic cumulant can suppress noise, which is equivalent to increasing the detected peak value, so that the FOCCCAF-based peak under the same echo SNR is more prominent, which is beneficial to the detection of echo signals. The theoretical and actual detection performance are consistent, verifying the superiority of FOCCCAF-based detection performance to CAF.
In order to verify the impact of false alarm probability on the performance of the proposed method, the detection performance of the method under different false alarm probabilities is simulated. GPS, DVB-S and INMARSAT satellite signals are used for simulation experiments. coherent accumulation point is 10 9 , 2000 Monte Carlo simulation experiments are carried out and the simulation results are shown in Figure 6.
From Figure 6, it can be seen that, under the same SNR, the greater the false alarm probability, the greater the detection probability. The simulation gives the relationship between the echo SNR and the detection probability of the false alarm probability at 10 4 , 10 5 and 10 6 . When the echo SNR is around −75 dB and the accumulated point reaches 10 9 , the detection probability of the echo signal reaches 1. When the echo SNR is −95 dB, the detection probability is 0. This is because the echo SNR is too low, even if the accumulated points reach 10 9 , the detected peak cannot be highlighted. To detect weak echo signals at this SNR, we need more points to accumulate.
In order to verify the impact of SDR on detection performance, GPS, DVB-S and INMARSAT satellite signals are used for simulation experiments. P f a = 10 4 , coherent accumulation point is 10 9 , 2000 Monte Carlo simulation experiments are carried out and the simulation results are shown in Figure 7.
In order to verify the impact of SDR on detection performance, GPS, DVB-S and INMARSAT satellite signals are used for simulation experiments. P f a = 10 4 , coherent accumulation point is 10 9 , 2000 Monte Carlo simulation experiments are carried out and the simulation results are shown in Figure 7.
It can be seen from Figure 7 that the detection performance of the echo signal is becoming better and better with the increase of the echo direct wave power ratio. This is because under the condition that the direct wave power is constant, as the echo direct wave power ratio increases, the echo signal is enhanced, and the detected peak can be highlighted. This is equivalent to improving the SNR, so the detection of echo signals is promoted. The correctness of the theory is verified by actual simulation.
In order to compare the detection performance with different methods, GPS, DVB-S and INMARSAT satellite signals are used for simulation experiments. P f a = 10 4 , coherent accumulation point is 10 9 , 2000 Monte Carlo simulation experiments are carried out and the simulation results are shown in Figure 8. From Figure 8, it can be seen that under the same conditions, the detection probability based on FOCCCAF which is proposed in this paper is better than that of CAF and CCAF in [16].

7. Conclusions

In this paper, a novel method of the passibe moving aerial target detection has been proposed under collaborative heterogeneous multiple satellites, in which parallel band-pass filters and an NLMS-based adaptive filter based on NLMS are used to separate direct wave signals in the reference channel and suppress DPI/MPI in the surveillance channel, respectively. Then, the detection statistic of weak echo signals from multiple heterogeneous satellites was obtained by using the anti-interference property of FOOCCCAF. Through distributed multi-sensors fusion, multiple detection results have been superimposed to obtain final detection results. Finally, depending on the theoretical analysis of the detection statistics, the relationship between the detection performance and a series of key parameters is studied. Numerical results have demonstarted that the proposed method is able to precisely detect a moving target in multiple heterogeneous satellite environments. This indicated that in a complicated environment with multiple satellites, the proposed method is more applicable than the methods relying on a single satellite as the illumination for passive detection.

Author Contributions

M.L. conceptualized and performed the algorithm and wrote the paper; K.L. and X.G. analyzed the experiment data; H.S. provided technical assistance to the research; Y.C. helped improve the language; F.G. is the research supervisor. The manuscript was discussed by all co-authors. 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 (61501348 and 61801363), the Shaanxi Provincial Key Research and Development Program (2019GY-043), the Joint Fund of Ministry of Education of the People’s Republic of China (6141A02022338), the China Postdoctoral Science Foundation (2017M611912), the 111 Project (B08038), and the China Scholarship Council (201806965031).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Proof of Lemma 1

In Equation (33), the real part Re χ i ( l , f d ) of the CAF under assumption H 1 is expressed as
Re χ i ( l , f d ) = m = 0 N 1 Re { y ( m ) x i ( m l ) e j 2 π f d m } = m = 0 N 1 Re { η = 1 η = i M r η s η ( m ζ ) r i s i ( m l ) e j 2 π f d η m e j 2 π f d m } + m = 0 N 1 Re { n ( m l ) η = 1 η = i M r η s η ( m ζ ) e j 2 π f d η m e j 2 π f d m } + m = 0 N 1 Re { n ( m ) r i s i ( m l ) e j 2 π f d m } + m = 0 N 1 Re { n ( m ) n ( m l ) e j 2 π f d m } .
The first term of Equation (A1) belongs to the real part of the CAF of the reference signal and the corresponding echo signal. It belongs to the determined detection statistics Re ( Λ ) and can be expressed as follows
m = 0 N 1 Re r i r η Re { η = 1 η = i M s η ( m ζ ) s η ( m l ) e j 2 π f d i m e j 2 π f d m } = Re r i r i Re { χ s i s i ( l ζ , f d f d i ) } ,
where χ s i s i ( l ζ , f d f d i ) represents the self-ambiguity function of the radiation source signal, when ζ is equal zero and f d i is equal zero, we can get
χ s i s i ( l , f d ) = m = 0 N 1 s i ( m ) s i ( m l ) e j 2 π f d m .
Equation (A3) is a self-ambiguity function of s i ( m ) and s i ( m ) is the signal of the normalized amplitude, so the amplitude of s i ( m ) is equal to one and E s i ( m ) is equal to zero. Therefore, when ζ is equal zero and f d i is equal zero, the self-ambiguity function obtains the maximum value, which is expressed as
χ s i s i ( 0 , 0 ) = N .
The second and third terms of Equation (A1) represent the real part of the CAF between the noise and the satellite illumination signal, which is equivalent that Gaussian noise is linearly transformed, so the second term still subject to Gaussian distribution of the following statistics characteristic
m = 0 N 1 Re { n ( m l ) η = 1 η = i M r η s η ( m ζ ) e j 2 π f d η m e j 2 π f d m } N ( 0 , η = 1 M N r η 2 σ n 2 ) ,
and
m = 0 N 1 Re { n ( m ) r i s i ( m l ) e j 2 π f d m } N ( 0 , N r i 2 σ n 2 ) .
The fourth term of Equation (A1) represents the cross correlation between the reference channel and the noise of the surveillance channel. Since they two are statistically independent, so this term is zero under the circumstance that N is large sufficiently. To sum up, under the assumption H 1 , Re χ i ( l , f d ) subject to the Gaussian distribution of the following statistical characteristic
Re χ i ( l , f d ) N ( Re r i r i Re { χ s i s i ( l ζ , f d f d i ) } , N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) ) .
In Equation (33), the imaginary part Im χ i ( l , f d ) of the CAF under the assumption H 1 is assumed as
Im χ i ( l , f d ) = m = 0 N 1 Im { y ( m ) x i ( m l ) e j 2 π f d m } .
From Equation (A8), we can obtain
Im χ i ( l , f d ) = m = 0 N 1 Im { y ( m ) x i ( m l ) e j 2 π f d m } = m = 0 N 1 Im { η = 1 η = i M r η s η ( m ζ ) r i s i ( m l ) e j 2 π f d η m e j 2 π f d m } + m = 0 N 1 Im { n ( m l ) η = 1 η = i M r η s η ( m ζ ) e j 2 π f d η m e j 2 π f d m } + m = 0 N 1 Im { n ( m ) r i s i ( m l ) e j 2 π f d m } + m = 0 N 1 Im { n ( m ) n ( m l ) e j 2 π f d m } .
The first term of Equation (A9) belongs to the imaginary part of the CAF of the reference signal and its corresponding echo signal, which is a determined detection statistics Im ( Λ ) and it can be expressed as follows
m = 0 N 1 Im r i r η Im { η = 1 η = i M s η ( m ζ ) s η ( m l ) e j 2 π f d i m e j 2 π f d m } = Im r i r i Im { χ s i s i ( l ζ , f d f d i ) } .
The second and third terms of Equation (A9) represent the imaginary part of the CAF between the noise and the source signal, which is equivalent that Gaussian noise is linearly transformed, so the second term still subject to the Gaussian distribution of following statistics characteristic
m = 0 N 1 Im { n ( m l ) η = 1 η = i M r η s η ( m ζ ) e j 2 π f d η m e j 2 π f d m } N ( 0 , η = 1 M N r η 2 σ n 2 ) ,
and
m = 0 N 1 Im { n ( m ) r i s i ( m l ) e j 2 π f d m } N ( 0 , N r i 2 σ n 2 ) .
The fourth term of Equation (A9) represents the cross correlation between the reference channel and the noise of the surveillance channel, since they two are statistically independent, so this term is zero when N is large enough. To sum up, under the assumption H 1 , Im χ i ( l , f d ) subjects to the Gaussian distribution of the following statistical characteristics
Im χ i ( l , f d ) N ( Im r i r i Im { χ s i s i ( l ζ , f d f d i ) } , N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) ) .
From Equation (A7) and Equation (A13), under the assumption H 1 ,the test statistic Λ obeys the complex Gaussian distribution with mean μ 1 = r i r i χ s i s i ( l ζ , f d f d i ) and variance σ 1 2 = 2 N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) , which can be expressed as
( Λ / H 1 ) C N ( r i r i χ s i s i ( l ζ , f d f d i ) , 2 N ( r i 2 σ n 2 + η = 1 M r η 2 σ n 2 ) ) ,
where CN (.) represents complex Gaussian distribution.
Under the assumption H 0 , the real part of the CAF is given by
Re χ i ( l , f d ) = m = 0 N 1 Re { y ( m ) x i ( m l ) e j 2 π f d m } = m = 0 N 1 Re { η = 1 η i M r η s η ( m ζ ) r i s i ( m l ) e j 2 π f d η m e j 2 π f d m } + m = 0 N 1 Re { n ( m l ) η = 1 η i M r η s η ( m ζ ) e j 2 π f d η m e j 2 π f d m } + m = 0 N 1 Re { n ( m ) r i s i ( m l ) e j 2 π f d m } + m = 0 N 1 Re { n ( m ) n ( m l ) e j 2 π f d m }
Under the hypothesis H 0 , the different satellite signals are uncorrelated, so the first term of Equation (A15) can be ignored, the second and third terms of Equation (A15) represent of the imaginary part of the CAF between Gaussian noise and radiation signal. The second term and the third term still subject to the Gaussian distribution of the following statistical characteristic
m = 0 N 1 Re { n ( m l ) η = 1 η i M r η s η ( m ζ ) e j 2 π f d η m e j 2 π f d m } N ( 0 , N η = 1 η i M r η 2 σ n 2 ) ,
and
m = 0 N 1 Re { n ( m ) r i s i ( m l ) e j 2 π f d m } N ( 0 , N r i 2 σ n 2 ) .
The fourth term of Equation (A15) represents the cross correlation between the reference channel and the noise of the surveillance channel. Since they two are statistically independent, this term is zero when N is large enough. To sum up, under the assumption H 0 , Re χ i ( l , f d ) subjects to the Gaussian distribution of the following statistical characteristic
Re χ i ( l , f d ) N ( 0 , N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 ) ) .
Similarly, the imaginary part Im χ i ( l , f d ) under H 0 obeys the Gaussian distribution of the following statistical properties
Im χ i ( l , f d ) N ( 0 , N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 ) ) .
Under the assumption H 0 , the detection statistic obeys the following complex Gaussian distribution
( Λ / H 0 ) C N ( 0 , 2 N ( r i 2 σ n 2 + η = 1 η i M r η 2 σ n 2 ) ) .

Appendix A.2. Proof of Lemma 2

In the Equation (18), the expectation of M s i s i n s i α i f + f d i ( τ ) equals zero, and the variance of M s i s i n s i α i f + f d i ( τ ) is
var M s i s i n s i α i f + f d i ( τ ) = E M s i s i n s i α i f + f d i ( τ ) 2 E M s i s i n s i α i f + f d i ( τ ) 2 = E M s i s i n s i α i f + f d i ( τ ) 2 = 1 N 2 t 1 = 0 N 1 t 2 = 0 N 1 E τ = 0 N 1 s i 3 ( t 1 ) s i 3 ( t 2 ) n ( t 1 + τ ) n ( t 2 + τ ) e j 2 π ( α i f + f d i ) ( t 2 t 1 ) = 1 N σ n 2 σ s i 6 ,
where σ n 2 represents the variance of n ( t ) , σ s i 2 denotes the variance of s i ( t ) . From Equation (A21), the M s i s i n s i α i f + f d i ( τ D i ) in the fifth term and the M s i s i s i n α i f ( τ ) in the seventh term obey the Gaussian distribution of the following statistical characteristics
M s i s i n s i α i f + f d i ( τ D i ) N ( 0 , 1 N σ n 2 σ s i 6 ) ,
and
M s i s i s i n α i f ( τ ) N ( 0 , 1 N σ n 2 σ s i 6 ) ,
where σ n 2 is the variance of n ( t ) . From Equation (18), the R s i α i ( τ u ) in the fifth and seventh terms and the M 4 s i α i ( τ u ) in the sixth and eighth terms obey the Gaussian distribution of the following statistical characteristics
R s i α i ( τ u ) N ( 0 , 1 N σ s i 4 ) ,
and
M 4 s i α i ( τ u ) N ( 0 , 1 N σ s i 8 ) .
Let
Ω 1 = 9 A r i 4 r i e j π ( α i f + f d i ) D + M s i s i n s i α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ ,
using the above analysis results, we can know that E [ Ω 1 ] = 0 and
var [ Ω 1 ] = E { [ Ω 1 ] 2 } E { [ Ω 1 ] } 2 = E { [ Ω 1 ] 2 } = 81 A 2 r i 8 ( r i ) 2 E 1 N τ 1 = 0 N 1 n ( τ 1 v ; α ) M 3 s i ( τ 1 v ; α ) R s i ( τ 1 v ; α ) e j 2 π τ 1 · 1 N τ 2 = 0 N 1 n ( τ 2 v ; α ) M 3 s i ( τ 2 v ; α ) R s i ( τ 2 v ; α ) e j 2 π τ 2 = 81 A 2 r i 8 ( r i ) 2 N 2 E τ 1 = 0 N 1 τ 2 = 0 N 1 n ( τ 1 v ; α ) M 3 s i ( τ 1 v ; α ) R s i ( τ 1 v ; α ) e j 2 π τ 1 · n ( τ 2 v ; α ) M 3 s i ( τ 2 v ; α ) R s i ( τ 2 v ; α ) e j 2 π τ 2 = 81 A 2 r i 8 ( r i ) 2 N 2 τ 1 = 0 N 1 τ 2 = 0 N 1 E V n ( τ 1 v ; α ) n ( τ 2 v ; α ) R s i ( τ 1 v ; α ) R s i ( τ 2 v ; α ) · M 3 s i ( τ 1 v ; α ) M 3 s i ( τ 2 v ; α ) e j 2 π ( τ 1 τ 2 ) = 81 A 2 r i 8 ( r i ) 2 N 3 σ n 2 σ s i 10 ,
and Ω 1 obeys the Gaussian distribution of the following statistical characteristics
Ω 1 N 0 , 81 A 2 r i 8 ( r i ) 2 N 3 σ n 2 σ s i 10 .
Let
Ω 2 = 3 r i 6 r i e j π ( α i f + f d i ) D + M s i s i n s i α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ .
Similarly, Ω 2 obeys the Gaussian distribution of the following statistical characteristics
Ω 2 N 0 , 9 r i 12 ( r i ) 2 N 3 σ n 2 σ s i 14 .
Let
Ω 3 = 3 A r i 5 + M s i s i s i n α i f ( τ ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) e j π f τ d τ .
Similarly, Ω 3 obeys the Gaussian distribution of the following statistical characteristics
Ω 3 N 0 , 9 A 2 r i 10 N 3 σ n 2 σ s i 10 .
Let
Ω 4 = r i 7 + M s i s i s i n α i f ( τ ) ( M 4 s i α i ( τ u ) ) e j π f τ d τ .
Similarly, Ω 4 obeys the Gaussian distribution of the following statistical characteristics
Ω 4 N 0 , r i 14 N 3 σ n 2 σ s i 14 .
In summary, under the assumption H 1 , the detection statistic Λ is obeyed the complex Gaussian distribution, which is expressed as
Λ / H 1 C N ( r i 7 r i e j π ( α i f f d i ) D i χ 4 s i α i f f d i , α i ( u , f ) + 9 A B r i 3 r i e j π ( α i f f d i ) D i χ 2 s i α i f f d i , α i ( u , f ) , σ s i 10 N 3 ( 81 A 2 r i 8 ( r i ) 2 σ n 2 + 9 r i 12 ( r i ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) ) .
Under the assumption H 0 ,
χ y , x i α i f , α i ( u , f ) = + C x i x i x i y α i f ( C 4 x i α i ( τ u ) ) e j π f τ d τ = r i 7 r η e j π ( α i f + f d i ) D i + M s i s i s i s η α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ + 9 A B r i 3 r η e j π ( α i f + f d i ) D i + R s i s η α i f + f d i ( τ D i ) ( R s i α i ( τ u ) ) e j π ( f f d i ) τ d τ 3 B r i 5 r η e j π ( α i f + f d i ) D i + ( M s i s i s i s η α i ( τ u ) ) R s i α i f + f d i ( τ D i ) e j π ( f f d i ) τ d τ 3 A r i 5 r η e j π ( α i f + f d i ) D i + M s i s i s i s η α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ η = 1 η i M 9 A r i 4 r η e j π ( α i f + f d i ) D i + M s i s i n s η α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ + η = 1 η i M 3 r i 6 r η e j π ( α i f + f d i ) D + M s i s i n s η α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ 3 A r i 5 + M s i s i s i n α i f ( τ ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) e j π f τ d τ + r i 7 + M s i s i s i n α i f ( τ ) ( M 4 s i α i ( τ u ) ) e j π f τ d τ .
Since there is no echo corresponding to the i t h direct wave in the supervise channel, that is η i . At this time, there is no correlation between M s i s i s i s η α i f + f d i ( τ D i ) and M 4 s i α i ( τ u ) , R s i s η α i f + f d i ( τ D i ) and R s i α i ( τ u ) , R s i α i f + f d i ( τ D i ) and M s i s i s i s η α i ( τ u ) , M s i s i s i s η α i f + f d i ( τ D i ) and R s i α i ( τ u ) . Let
Ω 1 = η = 1 η i M 9 A r i 4 r η e j π ( α i f + f d i ) D i + M s i s i n s η α i f + f d i ( τ D i ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) d τ ,
using the above analysis results, we can know that E [ Ω 1 ] = 0 , and
var [ Ω 1 ] = E { [ Ω 1 ] 2 } E { [ Ω 1 ] } 2 = E { [ Ω 1 ] 2 } = η = 1 η i M 81 A 2 r i 8 ( r η ) 2 E 1 N τ 1 = 0 N 1 n ( τ 1 v ; α ) M 3 s i ( τ 1 v ; α ) R s i ( τ 1 v ; α ) e j 2 π τ 1 · 1 N τ 2 = 0 N 1 n ( τ 2 v ; α ) M 3 s i ( τ 2 v ; α ) R s i ( τ 2 v ; α ) e j 2 π τ 2 = η = 1 η i M 81 A 2 r i 8 ( r η ) 2 N 2 E τ 1 = 0 N 1 τ 2 = 0 N 1 n ( τ 1 v ; α ) M 3 s i ( τ 1 v ; α ) R s i ( τ 1 v ; α ) e j 2 π τ 1 · n ( τ 2 v ; α ) M 3 s i ( τ 2 v ; α ) R s i ( τ 2 v ; α ) e j 2 π τ 2 = η = 1 η i M 81 A 2 r i 8 ( r η ) 2 N 2 τ 1 = 0 N 1 τ 2 = 0 N 1 E V n ( τ 1 v ; α ) n ( τ 2 v ; α ) R s i ( τ 1 v ; α ) R s i ( τ 2 v ; α ) = η = 1 η i M 81 A 2 r i 8 ( r η ) 2 N 2 N σ n 2 1 N σ s i 6 1 N σ s i 4 = η = 1 η i M 81 A 2 r i 8 ( r η ) 2 N 3 σ n 2 σ s i 10 ,
and Ω 1 obeys the Gaussian distribution of the following statistical characteristics
Ω 1 N 0 , η = 1 η i M 81 A 2 r i 8 ( r η ) 2 N 3 σ n 2 σ s i 10 .
Let
Ω 2 = η = 1 η i M 3 r i 6 r η e j π ( α i f + f d i ) D + M s i s i n s η α i f + f d i ( τ D i ) ( M 4 s i α i ( τ u ) ) e j π ( f f d i ) τ d τ .
Similarly, Ω 2 obeys the Gaussian distribution of the following statistical characteristics
Ω 2 N 0 , η = 1 η i M 9 r i 12 ( r η ) 2 N 3 σ n 2 σ s i 14 .
Let
Ω 3 = 3 A r i 5 + M s i s i s i n α i f ( τ ) ( R s i α i ( τ u ) e j π ( f f d i ) τ ) e j π f τ d τ .
Similarly, Ω 3 obeys the Gaussian distribution of the following statistical characteristics
Ω 3 N 0 , 9 A 2 r i 10 N 3 σ n 2 σ s i 10 .
Let
Ω 4 = r i 7 + M s i s i s i n α i f ( τ ) ( M 4 s i α i ( τ u ) ) e j π f τ d τ .
Similarly, Ω 4 obeys the Gaussian distribution of the following statistical characteristics
Ω 4 N 0 , r i 14 N 3 σ n 2 σ s i 14 .
In summary, under the assumption H 0 , the detection statistic Λ is obeyed the complex Gaussian distribution, which is expressed as
Λ / H 0 C N ( 0 , σ s i 10 N 3 ( η = 1 η i M 81 A 2 r i 8 ( r η ) 2 σ n 2 + η = 1 η i M 9 r i 12 ( r η ) 2 σ n 2 σ s i 4 + 9 A 2 r i 10 σ n 2 + r i 14 σ n 2 σ s i 4 ) ) .

References

  1. Hashim, I.; Al-Hourani, A.; Rowe, W.; Scott, J. Adaptive X-Band Satellite Antenna for Internet-of-Things (IoT) over Satellite Applications. In Proceedings of the 13th International Conference on Signal Processing and Communication Systems, Surfers Paradise, Australia, 16–18 December 2019; pp. 1–7. [Google Scholar]
  2. Li, F.; Lam, K.; Zhao, N.; Liu, X.; Zhao, K.; Wang, L. Spectrum Trading for Satellite Communication Systems with Dynamic Bargaining. IEEE Trans. Commun. 2018, 66, 4680–4693. [Google Scholar] [CrossRef]
  3. Ilyushin, Y.; Padokhin, A.; Smolov, V. Global Navigational Satellite System Phase Altimetry of the Sea Level: Systematic Bias Effect Caused by Sea Surface Waves. In Proceedings of the 2019 PhotonIcs & Electromagnetics Research Symposium—Spring (PIERS-Spring), Rome, Italy, 17–20 June 2019; pp. 1618–1627. [Google Scholar]
  4. Qiao, J.; Chen, W.; Ji, S.; Weng, D. Accurate and Rapid Broadcast Ephemerides for Beidou-Maneuvered Satellites. Remote. Sens. 2019, 11, 787. [Google Scholar] [CrossRef] [Green Version]
  5. Powell, S.; Akos, D. GNSS Reflectrometry Using the L5 and E5a Signals for Remote Sensing Applications. In Proceedings of the 2013 US National Committee of URSI National Radio Science Meeting, Boulder, CO, USA, 9–12 January 2013; p. 1. [Google Scholar]
  6. Clarizia, M.; Chotiros, N.; Vaccaro, M. A GPS-Reflectometry Simulator for Target Detection Over Oceans. In Proceedings of the IGARSS 2018–2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 23–27 July 2018; pp. 450–451. [Google Scholar]
  7. Garvanov, I.; Kabakchiev, C.; Behar, V.; Garvanova, M. Target Detection Using a GPS Forward-Scattering Radar. In Proceedings of the 2015 International Conference on Engineering and Telecommunication, Mosocw, Russia, 18–19 November 2015; pp. 29–33. [Google Scholar]
  8. Clarizia, M.; Braca, P.; Ruf, C.S.; Willett, P. Target Detection Using GPS Signals of Opportunity. In Proceedings of the 2015 18th International Conference on Information Fusion, Washington, DC, USA, 6–9 July 2015; pp. 1429–1436. [Google Scholar]
  9. Gronowski, K.; Samczynski, P.; Stasiak, K.; Kulpa, K. First, Results of Air Target Detection Using Single Channel Passive Radar Utilizing GPS Illumination. In Proceedings of the 2019 IEEE Radar Conference, Boston, MA, USA, 22–26 April 2019; pp. 1–6. [Google Scholar]
  10. Mojarrabi, B.; Homer, J.; Kubik, K. Power Budget Study for Passive Target Detection and Imaging Using Secondary Applications of GPS Signals in Bistatic Radar Systems. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002; pp. 449–451. [Google Scholar]
  11. Zeng, H.; Chen, J.; Wang, P.; Yang, W.; Liu, W. 2D Coherent Integration Processing and Detecting of Aircrafts Using GNSS-Based Passive Radar. Remote. Sens. 2018, 10, 1164. [Google Scholar] [CrossRef] [Green Version]
  12. AlJewari, Y.; Ahmad, R.; AlRawi, A. Impact of Multipath Interference and Change of Velocity on the Reliability and Precision of GPS. In Proceedings of the 2014 2nd International Conference on Electronic Design, Penang, Malaysia, 19–21 August 2014; pp. 427–430. [Google Scholar]
  13. Wu, X.; Gong, P.; Zhou, J.; Liu, Z. The Applied Research on Anti-multipath Interference GPS Signal Based on Narrow-related. In Proceedings of the 2014 IEEE 5th International Conference on Software Engineering and Service Science, Beijing, China, 27–29 June 2014; pp. 771–774. [Google Scholar]
  14. Suberviola, I.; Mayordomo, I.; Mendizabal, J. Experimental Results of Air Target Detection with a GPS Forward-Scattering Radar. IEEE Geosci. Remote Sens. Lett. 2011, 9, 47–51. [Google Scholar] [CrossRef]
  15. Lashley, M.; Bevly, D.; Hung, J. Performance Analysis of Vector Tracking Algorithms for Weak GPS Signals in High Dynamics. IEEE J. Sel. Top. Signal Process. 2009, 3, 661–673. [Google Scholar] [CrossRef]
  16. Liu, M.; Gao, Z.; Chen, Y.; Song, H.; Li, Y.; Gong, F. Passive Detection of Moving Aerial Target Based on Multiple Collaborative GPS Satellites. Remote. Sens. 2020, 12, 263. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, M.; Yi, F.; Liu, P.; Li, B. Cramer-Rao Lower Bounds of TDOA and FDOA Estimation Based on Satellite Signals. In Proceedings of the 2018 14th IEEE International Conference on Signal Processing, Beijing, China, 12–16 August 2018; pp. 1–4. [Google Scholar]
  18. Masjedi, S.; Moddares-Hashemi, M. Theoretical Approach for Target Detectionand Interference Cancellation in Passive Radar. IET Radar Sonar Navig. 2013, 7, 205–216. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Xi, S. Application of New LMS Adaptive Filtering Algorithm with Variable Step Size in Adaptive Echo Cancellation. In Proceedings of the 2017 IEEE 17th International Conference on Communication Technology, Chengdu, China, 27–30 October 2017; pp. 1715–1719. [Google Scholar]
  20. Niranjan, D.; Ashwini, B. Noise Cancellation in Musical Signals Using Adaptive Filtering Algorithms. In Proceedings of the 2017 International Conference on Innovative Mechanisms for Industry Applications, Bangalore, India, 21–23 February 2017; pp. 82–86. [Google Scholar]
  21. Mugdha, A.; Rawnaque, F.; Ahmed, M. A Study of Recursive Least Squares (RLS) Adaptive Filter Algorithm in Noise Removal from ECG Signals. In Proceedings of the 2015 International Conference on Informatics, Electronics Vision, Fukuoka, Japan, 15–18 June 2015; pp. 1–6. [Google Scholar]
  22. Cardinali, R.; Colone, F.; Ferretti, C.; Lombardo, P. Comparison of Clutter and Multipath Cancellation Techniques for passive Radar. In Proceedings of the 2007 IEEE Radar Conference, Boston, MA, USA, 17–20 April 2007; pp. 469–474. [Google Scholar]
  23. Liu, M.; Zhang, J.; Li, B. Feasibility Analysis of OFDM/OQAM Signals as Illuminator of Opportunity for Passive Detection. In Proceedings of the 2018 14th IEEE International Conference on Signal Processing, Beijing, China, 12–16 August 2018; pp. 1–4. [Google Scholar]
  24. Hong, J.; Wang, S.X.; Lu, H.J. An Effective Direction Estimation Algorithm in Multipath Environment Based on Fourth-order Cyclic Cumulants. In Proceedings of the IEEE 5th Workshop on Signal Processing Advances in Wireless Communications, Lisbon, Portugal, 11–14 July 2004; pp. 263–267. [Google Scholar]
  25. Yan, X.; Wang, S.; Wang, K.; Jiang, H. Localization of Near field Cyclostationary Source Based on Fourth-order Cyclic Cumulant. In Proceedings of the 2008 9th International Conference on Signal Processing, Beijing, China, 26–29 October 2008; pp. 1629–1632. [Google Scholar]
  26. Jun, H. A proper integral representation of Marcum Q-Function. In Proceedings of the 2014 XXXIth URSI General Assembly and Scientific Symposium, Beijing, China, 16–23 August 2014; pp. 1–3. [Google Scholar]
  27. Bocus, M.; Dettmann, C.; Coon, J. An Approximation of the First Order Marcum Q-Function with Application to Network Connectivity Analysis. IEEE Commun. Lett. 2013, 17, 499–502. [Google Scholar] [CrossRef] [Green Version]
  28. Ermolova, N.; Tirkkonen, O. Laplace Transform of Product of Generalized Marcum Q, Bessel I, and Power Functions With Applications. IEEE Trans. Signal Process. 2014, 62, 2938–2944. [Google Scholar]
  29. Hu, B.; Liu, M.; Yi, F.; Song, H.; Jiang, F.; Gong, F.; Zhao, N. DOA Robust Estimation of Echo Signals Based on Deep Learning Networks with Multiple Type Illuminators of Opportunity. IEEE Access 2020, 8, 14809–14819. [Google Scholar] [CrossRef]
  30. Liu, M.; Zhang, J.; Tang, J.; Jiang, F.; Liu, P.; Gong, F.; Zhao, N. 2D DOA Robust Estimation of Echo Signals Based on Multiple Satellites Passive Radar in the Presence of Alpha Stable Distribution Noise. IEEE Access 2019, 7, 2169–3536. [Google Scholar]
Figure 1. Passive detection system model based on multiple collaborative heterogeneous satellites.
Figure 1. Passive detection system model based on multiple collaborative heterogeneous satellites.
Remotesensing 12 01150 g001
Figure 2. Spectrogram of three time-domain overlapped satellite signals.
Figure 2. Spectrogram of three time-domain overlapped satellite signals.
Remotesensing 12 01150 g002
Figure 3. Spectrograms of three time-domain overlapped satellite signals after separation. (a) Spectrogram of GPS satellite signal; (b) Spectrogram of INMARSAT satellite signal; (c) Spectrogram of DVB-S satellite signal.
Figure 3. Spectrograms of three time-domain overlapped satellite signals after separation. (a) Spectrogram of GPS satellite signal; (b) Spectrogram of INMARSAT satellite signal; (c) Spectrogram of DVB-S satellite signal.
Remotesensing 12 01150 g003
Figure 4. Relationship between coherent accumulation points and detection performance.
Figure 4. Relationship between coherent accumulation points and detection performance.
Remotesensing 12 01150 g004
Figure 5. Detection performance of the proposed method.
Figure 5. Detection performance of the proposed method.
Remotesensing 12 01150 g005
Figure 6. Detection performance with different false alarm probabilities.
Figure 6. Detection performance with different false alarm probabilities.
Remotesensing 12 01150 g006
Figure 7. Relationship between SDR and detection performance.
Figure 7. Relationship between SDR and detection performance.
Remotesensing 12 01150 g007
Figure 8. Detection performance comparison with different methods.
Figure 8. Detection performance comparison with different methods.
Remotesensing 12 01150 g008

Share and Cite

MDPI and ACS Style

Liu, M.; Li, K.; Song, H.; Chen, Y.; Gao, X.; Gong, F. Using Heterogeneous Satellites for Passive Detection of Moving Aerial Target. Remote Sens. 2020, 12, 1150. https://doi.org/10.3390/rs12071150

AMA Style

Liu M, Li K, Song H, Chen Y, Gao X, Gong F. Using Heterogeneous Satellites for Passive Detection of Moving Aerial Target. Remote Sensing. 2020; 12(7):1150. https://doi.org/10.3390/rs12071150

Chicago/Turabian Style

Liu, Mingqian, Kunming Li, Hao Song, Yunfei Chen, Xiuhui Gao, and Fengkui Gong. 2020. "Using Heterogeneous Satellites for Passive Detection of Moving Aerial Target" Remote Sensing 12, no. 7: 1150. https://doi.org/10.3390/rs12071150

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