Hostname: page-component-586b7cd67f-2plfb Total loading time: 0 Render date: 2024-12-03T19:02:00.492Z Has data issue: false hasContentIssue false

Reference signal extraction based on multiprocessing method in dual-polarization antenna passive radar

Published online by Cambridge University Press:  08 November 2022

Yucheng Yi*
Affiliation:
School of Information Engineering, East China Jiao Tong University, Nanchang, Jiangxi, China
Lu Zhu
Affiliation:
School of Information Engineering, East China Jiao Tong University, Nanchang, Jiangxi, China
Long Liang
Affiliation:
Southwest China Research Institute of Electronic Equipment, Chengdu, Sichuan, China
Yan Fu
Affiliation:
Wuhan Maritime Communication Research Institute, Wuhan, Hubei, China
*
Author for correspondence: Yucheng Yi, E-mail: ycyi@ecjtu.edu.cn
Rights & Permissions [Opens in a new window]

Abstract

Extraction of reference signal is an indispensable step in the signal processing of polarization diversity passive radar (PDPR) based on a digital television signal. A conventional reference signal extraction method requires an additional reference antenna, which has a certain demand for space. Single dual-polarization antenna passive radar (SDPPR) systems do not require a reference antenna, and the radar station layout is flexible, which is suitable for a large-scale radar network. It is a main research direction of PDPR in future. However, its reference signal extraction needs to rely on the signal reconstruction method. When the signal to interference and noise ratio of the direct-path signal is relatively low, the signal reconstruction method will fail. In this paper, we propose a reference signal extraction method based on sub-carrier processing method, blind adaptive oblique projection technology, and extensive cancelation algorithm to solve the above problem. Experimental results show that the method proposed in this paper is a reasonable alternative after the failure of reference signal reconstruction, and it is an effective supplement to the reference signal extraction technology.

Type
Radar
Copyright
© The Author(s), 2022. Published by Cambridge University Press in association with the European Microwave Association

Introduction

Passive radar or passive coherent location systems have been a hotspot globally [Reference Colone, O'Hagan, Lombardo and Baker1Reference Tan, Sun, Lu, Lesturgle and Chan7]. Over the decades, illuminators of opportunity used in passive radar have been expanded rapidly, using analog signals such as frequency-modulated commercial radios [Reference Wan, Yi, Zhao and Ke2] as well as digital video broadcasting [Reference Poullin4], digital audio broadcasting [Reference Coleman and Yardley5], wireless LAN signals [Reference Colone, Woodbridge, Guo, Mason and Baker6], global systems for mobile communication signals [Reference Tan, Sun, Lu, Lesturgle and Chan7], and china mobile multimedia broadcasting [Reference Wan, Yi, Zhao and Ke2] and even LTE signals [Reference Bartoletti, Conti and Win8]. With the rapid development of 5G technology [Reference Zhao, Ni, Yang, Zhang, Gong and Yu9,Reference Zhao, Liu, Yang, Ai and Ni10], it will inject vitality into the development of passive radar. One difference between active radar and passive radar is whether the reference signal is extracted. Thus, passive radar systems usually have two types of channels: one is used to receive echo signals of targets called surveillance channel. The other one is used to receive a direct-path signal of illuminators of opportunity called reference channel. The pure reference signal can be used to remove disturbance of the direct-path signal and strong multi-path clutter by time-domain filtering [Reference Colone, O'Hagan, Lombardo and Baker1] and improve the target signal to noise ratio (SNR) by match filtering [Reference Wan, Yi, Zhao and Ke2]. Thus, reference signal extraction is a necessary task in passive radar signal processing flow. The reference signal extraction methods in conventional passive radar based on frequency-modulated signals usually adopt multi-path suppression techniques, such as spatial beam-forming, adaptive filtering in the time domain, and so on [Reference Colone, Cardinali, Lombardo, Crognale, Cosmi, Lauri and Bucciarelli11]. But the signal purity is limited by antenna directivity and the correlation between clutters and direct-path signal. According to passive radar based on digital television signals, the authors in [Reference Palmer, Harms, Searle and Davis12Reference Wan, Wang, Hong and Tang16] propose the reference signal reconstruction (RSR) method. The basic signal reconstruction process is shown in Fig. 1.

Fig. 1. Signal reconstruction process in passive radar.

The application of conventional passive radars may be challenged by the attributes of the receiving antenna system in terms of size, weight, and cost. In order to reduce the strength of the direct-path signal in the surveillance channel, the angle of the reference antenna and the surveillance antenna need to be separated. This will increase the space requirement of the antenna system, and makes it impossible to setup a passive radar system on certain space-constrained platforms. Meanwhile, cost control is extremely important to realize the large-scale network of passive radar. A single-antenna digital television-based passive radar (SDPR) system can effectively solve the above problems [Reference Searle, Davis and Palmer17]. In order to further improve the anti-jamming capability of the radar system, we propose a single-antenna scheme of the polarization diversity passive radar [Reference Yi, Wan, Yi and Cao18, Reference Yi, Wan, Yi and Cao19] system on the basis of a slight increase in the cost. Different from the conventional passive radar mode that extract the reference signal from reference channel, the SDPR or the single dual-polarization antenna passive radar (SDPPR) extracts reference signal from the surveillance channel. Fortunately, the RSR method is also applicable to the surveillance channel of passive radar, and the signal can be completely recovered under certain signal to interference and noise ratio (SINR) of the direct-path signal conditions (the signal refers to direct-path signal, the interference include multi-path clutter and target signal). However, due to the complex propagation environment and antenna side-lobe attenuation, sometimes SINR of direct-path signal is not ideal in the surveillance channel. The decrease in the SINR of direct-path signal will reduce the purity of reconstructed signal, and even make the signal unable to be reconstructed. This will bring challenges to the RSR method. Assuming that the RSR in surveillance channel cannot be completed, SDPPR system will fail.

In this paper, a new method based on the orthogonal frequency division multiplex (OFDM) modulation of the digital television signal [Reference Malanowski, Kulpa and Misiurewicz3], blind adaptive oblique projection technology [Reference Behrens and Scharf20Reference Cao, Zhang, Liang, Wen, Jin and Zhang24], and extensive cancelation algorithm (ECA) [Reference Wan, Yi, Zhao and Ke2] is explored to address this difficulty. It is interesting to note that the inter-carrier orthogonality in a multi-path environment means that the zero-Doppler signals are coherent with each other in each sub-carrier. The coherence feature indicates that the direct-path signal and multi-path clutter only occupy one degree of freedom in sub-carrier domain. This feature indicates that the polarization state of composite clutter signal is coherently combined in each sub-carrier. As an effective signal separation method, the largest difference between blind adaptive oblique projection technology and conventional oblique projection in polarization signal processing is that it does not need to obtain the target signal polarization state. It is possible to construct an oblique projection operator only by the composite clutter signal polarization state and achieve signal separation. In order to obtain a more pure reference signal, we use the ECA to filtering the extra multi-path clutter signals and leaving the direct-path signal only. The experimental result shows that the proposed method can extract the reference signal effectively and obtain well target detection effect.

The paper is organized as follows. Section “Signal model” addresses the signal model of SDPPR in the surveillance channel. Section “Theoretical basis” describes the principle of sub-carrier processing, oblique projection, and ECA. Section “Multi-processing method of reference signal extraction” proposes the multi-processing method of the reference signal extraction. In Section “Error analysis,” the simulation analyses prove that the method of reference signal extraction-based multi-processing is effective. Section “Experimental research” shows the comparative experimental results. Section “Conclusion” draws conclusions based on the experimental results.

Signal model

SDPPR is a single-antenna passive radar system with polarization diversity capability. The received signal of SDPPR contains direct-path signal, strong clutter, and target echoes.

In this paper, scalars are represented by italic letters, and vectors are represented by bold letters. Suppose the received signal noise is white Gaussian noise. The signal model of SDPPR can expressed in the Jones vector form:

(1)$$\eqalign{{\boldsymbol S}_J = & \left[{\matrix{ {E_H} \cr {E_V} \cr } } \right] = \left[{\matrix{ {E_{dh}} \cr {E_{dv}} \cr } } \right]c_1d( t) + \sum\limits_{n = 2}^{M_c} {\left[{\matrix{ {E_{nh}} \cr {E_{nv}} \cr } } \right]} c_nd( t-\tau _n) \cr & \!\! + \sum\limits_{m = 1}^{M_t} {\left[{\matrix{ {E_{mh}} \cr {E_{mv}} \cr } } \right]} \alpha _md( t-\tau _m) e^{\,j2\pi f_m^D t_m} + \left[{\matrix{ {w_h} \cr {w_v} \cr } } \right]} $$

where d(t) is the signal from transmitter. Edh, Enh, and Emh represent direct-path signal, clutters, and targets in the horizontal (H) polarization component, respectively. Edv, Env, and Emv represent direct-path signal, clutters, and targets in the vertical (V) polarization component, respectively. wh and wv represent the Gaussian white noise in the H-polarized channel and V-polarized channel, respectively. c 1, cn, and αm represent the amplitude of direct-path signal, clutters, and targets, respectively. τn and τm represent the delay with respect to the direct-path signal. $f_m^D$ represents the Doppler frequency of the m-th target. Mc and Mt represent the number of clutter and target, respectively.

Theoretical basis

This section will elaborate on the basic principles of sub-carrier processing, oblique projection technology, and ECA.

Sub-carrier processing

Taking a digital television signal based on cyclic prefix (CP-OFDM) as an example, the received signal in the useful duration at k-th sub-carrier can be expressed as [Reference Yi, Wan, Yi and Cao18]:

(2)$$\eqalign{S_k & = c_1Q_k^T + \sum\limits_{n = 2}^{M_c} {c_ne^{{-}j2\pi f_k\tau _n^c }} Q_k^T \cr & \quad + \sum\limits_{m = 1}^{M_t} {\alpha _me^{{-}j2\pi ( {\,f_k-f_m^D } ) \tau _m^t }} U_{k, m}^T + W_k} $$

where Qk = [C 0, k, C 1, k,…, CL −1, k]T, L is the OFDM block number used as snapshots, C 1,k represent the normalized modulation code of the l-th OFDM block in k-th subcarrier, and

(3)$$U_{k, l}\approx [ {C_{0, k}, \;C_{1, k}e^{\,j2\pi f_l^D T_s}, \;\ldots , \;C_{L-1, k}e^{\,j2\pi f_l^D T_s( {L-1} ) }} ] ^T$$

where TS is the duration of one OFDM block. It can be seen from (2) that the samples of the direct-path signal and multi-path clutters at one sub-carrier are coherent with each other, whereas target echoes are uncorrelated to the clutter signals due to phase rotation introduced by the Doppler frequency. Therefore, all the zero-Doppler signals as shown in the first term of the right-hand side of (4) are integrated into a single term:

(4)$$\eqalign{S_k = & \left[{c_1 + \sum\limits_{n = 2}^{M_c} {c_ne^{{-}j2\pi f_k\tau_n^c }} } \right]Q_k^T \cr & \quad + \sum\limits_{m = 1}^{M_t} {\alpha _me^{{-}j2\pi ( {\,f_k-f_m^D } ) \tau _m^t }} U_{k, m}^T + W_k} $$

Oblique projection

Oblique projection operator is an extension of the orthogonal projection operators, which has attracted much attention in the field of signal processing [Reference Behrens and Scharf20]. Oblique projection technology has been used to interference and clutter suppression [Reference Hou, Mao and Li25, Reference Yi, Ji, Kirubarajan, Xie and Hu26], direction of arrival estimation [Reference Boyer and Bouleux21, Reference Tannaka and Imai27], target detection [Reference Sun, Xie, Hao and Yi28], and polarization diversity technology [Reference Mao, Liu, Hou, Hong, Guo and Deng22, Reference Cao, Zhang, Liang, Wen, Jin and Zhang24, Reference Wang, Zhang and Cao29]. Assume that S is a complex matrix of size n × m having full column rank and, likewise, that J is a complex matrix of size n × k having full column rank. Assume further that subspace 〈S〉 and 〈J〉 are disjoint, which requires m + k ≤ n. The H represents conjugate transpose symbols. The well-known formula to build a projection with subspace 〈S〉 is

(5)$${\bf P}_S = {\bf S}( {\bf S}^H{\bf S}) ^{{-}1}{\bf S}^H$$

where (•)H represents conjugate transposition. (•)−1 represents inverse operation. The orthogonal projection with subspace 〈S is given by

(6)$${\bf P}_S^\bot = {\bf I}-{\bf P}_S$$

Now examine the matrix (S J) that is composed from the columns of S and of J. Assume that m + k < n, and the column rank of (S J) is less than n. The orthogonal projection onto the linear subspace 〈S J〉 is

(7)$${\bf P}_{SJ} = \left[{\matrix{ {\bf S} & {\bf J} \cr } } \right]\left[{\matrix{ {{\bf S}^H{\bf S}} & {{\bf S}^H{\bf J}} \cr {{\bf J}^H{\bf S}} & {{\bf J}^H{\bf J}} \cr } } \right]^{{-}1}\left[{\matrix{ {{\bf S}^H} \cr {{\bf J}^H} \cr } } \right]$$

We may decompose this orthogonal projection as follows:

(8)$${\bf P}_{SJ} = {\bf E}_J + {\bf E}_S$$
(9)$${\bf E}_J = \left[{\matrix{ {\bf S} & 0 \cr } } \right]\left[{\matrix{ {{\bf S}^H{\bf S}} & {{\bf S}^H{\bf J}} \cr {{\bf J}^H{\bf S}} & {{\bf J}^H{\bf J}} \cr } } \right]^{{-}1}\left[{\matrix{ {{\bf S}^H} \cr {{\bf J}^H} \cr } } \right]$$
(10)$${\bf E}_S = \left[{\matrix{ 0 & {\bf J} \cr } } \right]\left[{\matrix{ {{\bf S}^H{\bf S}} & {{\bf S}^H{\bf J}} \cr {{\bf J}^H{\bf S}} & {{\bf J}^H{\bf J}} \cr } } \right]^{{-}1}\left[{\matrix{ {{\bf S}^H} \cr {{\bf J}^H} \cr } } \right]$$

Formulas (9) and (10) may be simplified as follows:

(11)$${\bf E}_J = {\bf S}( {\bf S}^H{\bf P}_J^\bot {\bf S}) ^{{-}1}{\bf S}^H{\bf P}_J^\bot $$
(12)$${\bf E}_S = {\bf J}( {\bf J}^H{\bf P}_S^\bot {\bf J}) ^{{-}1}{\bf J}^H{\bf P}_S^\bot $$

The subspaces 〈S〉 and 〈J〉 are disjoint, so

(13)$${\bf E}_J\left[{\matrix{ {\bf S} & {\bf J} \cr } } \right] = \left[{\matrix{ {\bf S} & 0 \cr } } \right]$$
(14)$${\bf E}_S\left[{\matrix{ {\bf S} & {\bf J} \cr } } \right] = \left[{\matrix{ 0 & {\bf J} \cr } } \right]$$

As shown in (13) and (14), the oblique projection operator EJ and ES are orthogonal to the subspace 〈J〉 and 〈S〉, respectively. The subspace 〈S〉 and 〈J〉 are not required to be orthogonal.

Extensive cancelation algorithm

ECA is a clutter suppression algorithm based on interference subspace projection [Reference Wan, Yi, Zhao and Ke2]. The algorithm is to project the echo signals of the surveillance channel into the direct-path signal and its time delay spread space to eliminate clutter. Construct a subspace matrix X consisting of the reference signal as follows:

(15)$${\bf X} = \left[{\matrix{ {S_{ref}( 1 ) } & {S_{ref}( 2 ) } & \cdots & {S_{ref}( N ) } \cr } } \right]^T$$

where Sref(n), n = 1, 2, …, N represent reference signal and (•)T represents transposition.

The subspace projection method is used to suppress the direct-path signal and multi-path clutter in the surveillance channel. The orthogonal projection matrix can be represented by the subspace matrix X formed by the reference signal

(16)$${\bf P} = {\bf I}-{\bf X}( {{\bf X}^H{\bf X}} ) ^{{-}1}{\bf X}^H$$

Projecting the echo signal of the surveillance channel to the orthogonal subspace of the reference signal, the residual signal S ECA in the surveillance channel can be expressed as follows:

(17)$$S_{ECA} = [ {{\bf I}-{\bf X}{( {{\bf X}^H{\bf X}} ) }^{{-}1}{\bf X}^H} ] S_{surv} = {\bf P}S_{surv}$$

After the orthogonal projection processing, the direct-path signal and clutter in the surveillance channel are canceled, and the target signal is retained.

Multi-processing method of reference signal extraction

The SDPPR is a special form of bistatic/multistatic radar system. The illuminators cannot be controlled and are given priority to cover urban areas, thus the surveillance signal contain direct-path signal, strong clutter from environment echoes, and target echoes. Affected by the low SINR of direct-path signal in the surveillance channel, the RSR method is likely to fail, and the radar system will not be able to detect target at this time. In this section, we will introduce a multi-processing method of the reference signal extraction.

The ideal reference channel signal model is shown in (18):

(18)$${\bf S}_J^{ref} = \left[{\matrix{ {E_{dh}^{ref} } \cr {E_{dv}^{ref} } \cr } } \right]l_1d( t) + \left[{\matrix{ {w_h^{ref} } \cr {w_v^{ref} } \cr } } \right]$$

Comparing (1) and (18), we observed that the target and multi-path clutter must be suppressed in order to obtain the ideal reference signal in the surveillance channel. A new reference signal extraction method is proposed in two steps.

Target suppression

It is assumed that the surveillance signal only contains direct-path signal, target signal, and Gaussian white noise, and the polarization state vectors of direct-path signal and target signal can be obtained by estimation, and are expressed as follows:

(19)$${\bf D} = \left[{\matrix{ {\cos \gamma_d} \cr {\sin \gamma_de^{\,j\varphi_d}} \cr } } \right]$$
(20)$${\bf T} = \left[{\matrix{ {\cos \gamma_t} \cr {\sin \gamma_te^{\,j\varphi_t}} \cr } } \right]$$

Using the method of oblique projection in Section “Theoretical basis,” the oblique projection operator can be express as

(21)$${\bf E}_{op} = {\bf D}( {{\bf D}^H{\bf P}_T^\bot {\bf D}} ) ^{{-}1}{\bf D}^H{\bf P}_T^\bot $$
(22)$${\bf E}_{op}\left[{\matrix{ {\bf D} & {\bf T} \cr } } \right] = \left[{\matrix{ {\bf D} & 0 \cr } } \right]$$

In (22), the target signal is suppressed. But under actual conditions, the energy of target signal is very low, and its polarization state vector is difficult to estimate. In order to solve this problem, blind adaptive oblique projection technology is used for reference [Reference Behrens and Scharf20Reference Cao, Zhang, Liang, Wen, Jin and Zhang24].

We define the surveillance signal covariance matrix as below:

(23)$${\bf R}_{{\rm sur}v} = {\bf R}_{DT} + \sigma ^2{\bf {\rm I}}{\rm} = \left[{\matrix{ {\bf D} & {\bf T} \cr } } \right]\left[{\matrix{ {R_D} & 0 \cr 0 & {R_T} \cr } } \right]\left[{\matrix{ {\bf D} & {\bf T} \cr } } \right]^H + \sigma ^2{\bf {\rm I}}$$

where ${\bf R}_{{\rm sur}v} = {\bf S}_{surv}{\bf S}_{surv}^H$, ${\bf S}_{surv} = {\bf D}s_D( t ) + {\bf T}s_T( t ) + {\bf W}_{surv}$, s D(t) = cd(t), d(t) is the signal from transmitter, c represents the amplitude of direct-path signal; $s_T( t ) = \alpha d( t-\tau ) e^{j2\pi f_{}^D t}$, α represents the amplitude of target, τ represents the delay with respect to the direct-path signal, $f_{}^D$ represents the Doppler frequency of the target; ${\bf W}_{surv}$ $\left({{\bf W}_{surv} = \left[{\matrix{ {w_H} & {w_V} \cr } } \right]} \right)$ represents the Gaussian white noise of H-polarized channel and V-polarized channel; ${\bf R}_{DT} = [ {{\bf D}s_D( t ) + {\bf T}s_T( t ) } ] [ {{\bf D}s_D( t ) + {\bf T}s_T( t ) } ] ^H$, RD (R D = [s D(t)][s D(t)]H) is the eigenvalue of RDT related to the direct-path signal. RT (R T = [s T(t)][s T(t)]H) is the eigenvalue of RDT related to the target. σ 2 represents noise variance. I represents the identity matrix.

Then, we use the method proposed in [Reference Coon, Sandell, Beach and Mcgeehan30, Reference Boumard31] to estimate the noise variance σ 2, and get RDT as shown below:

(24)$${\bf R}_{DT} = {\bf R}_{surv}-\sigma ^2{\bf {\rm I}} = \left[{\matrix{ {\bf D} & {\bf T} \cr } } \right]\left[{\matrix{ {R_D} & 0 \cr 0 & {R_T} \cr } } \right]\left[{\matrix{ {\bf D} & {\bf T} \cr } } \right]^H$$

According to the principle of blind adaptive oblique projection, the modified oblique projection operator as follows:

(25)$${\bf E}_{b-op}^{} = {\bf D}( {{\bf D}^H{\bf R}_{DT}^ + {\bf D}} ) ^{{-}1}{\bf D}^H{\bf R}_{DT}^ + $$

The symbol (•)+ represents the pseudo-inverse. We define

(26)$${\bf A} = \left[{\matrix{ {\bf D} & {\bf T} \cr } } \right]$$
(27)$${\bf B} = \left[{\matrix{ {R_D} & 0 \cr 0 & {R_T} \cr } } \right]$$

Then, we can obtain the pseudo-inverse of matrix A:

(28)$${\bf A}^ + { = } \left[{\matrix{ {{( {{\bf D}^H{\bf P}_T^\bot {\bf D}} ) }^{{-}1}{\bf D}^H{\bf P}_T^\bot } \cr {{( {{\bf T}^H{\bf P}_D^\bot {\bf T}} ) }^{{-}1}{\bf T}^H{\bf P}_D^\bot } \cr } } \right]$$

where ${\bf P}_T^\bot$ represents the orthogonal projection matrix of T. ${\bf P}_D^\bot$ represents the orthogonal projection matrix of D. The covariance matrix of surveillance channel signal without noise can be simplified as

(29)$${\bf R}_{DT}^ + { = } ( {{\bf AB}{\bf A}^H} ) ^ + $$

In order to remove the target signal and extract the reference signal, we need

(30)$${\bf E}_{b-op}^{} \left[{\matrix{ {\bf D} & {\bf T} \cr } } \right] = \left[{\matrix{ {\bf D} & 0 \cr } } \right]$$

We can plug (28) and (29) into (25), and get

(31)$$\eqalign{& {\bf E}_{b-op} = {\bf D}({\bf D}^H{\bf R}_{DT}^ + {\bf D})^{-1}{\bf D}^H({\bf AB}{\bf A}^H)^ + \cr & \quad \quad \quad = \;{\bf D}({\bf D}^H{\bf R}_{DT}^ + {\bf D})^{-1}({\bf A}^ + {\bf D})^H{\bf B}^{-1}{\bf A}^ + \cr & \quad = {\bf D}({\bf D}^H{\bf R}_{DT}^ + {\bf D})^{-1}\left[ {\matrix{ {{({\bf D}^H{\bf P}_T^\bot {\bf D})}^{-1}{\bf D}^H{\bf P}_T^\bot {\bf D}} \hfill \cr {{({\bf T}^H{\bf P}_D^\bot {\bf T})}^{-1}{\bf T}^H{\bf P}_D^\bot {\bf D}} \hfill \cr } } \right]^H \cr & \quad {\bf B}^{-1}\left[ {\matrix{ {{({\bf D}^H{\bf P}_T^\bot {\bf D})}^{-1}{\bf D}^H{\bf P}_T^\bot } \hfill \cr {{({\bf T}^H{\bf P}_D^\bot {\bf T})}^{-1}{\bf T}^H{\bf P}_D^\bot } \hfill \cr } } \right] \cr & \quad = {\bf D}({\bf D}^H{\bf R}_{DT}^ + {\bf D})^{-1}\left( {\left[ {\matrix{ {\bf I} \hfill \cr 0 \hfill \cr } } \right]} \right)^H{\bf B}^{-1}\left[ {\matrix{ {{({\bf D}^H{\bf P}_T^\bot {\bf D})}^{-1}{\bf D}^H{\bf P}_T^\bot } \hfill \cr {{({\bf T}^H{\bf P}_D^\bot {\bf T})}^{-1}{\bf T}^H{\bf P}_D^\bot } \hfill \cr } } \right]} $$

Then, combined with (30) and (31), we can obtain

(32)$$\eqalign{{\bf E}_{b-op}^{} {\bf T} & = {\bf D}({\bf D}^{H}{\bf R}_{DT}^{+} {\bf D} ) ^{-1} \left({\left[{\matrix{{\bf I} \cr 0 \cr}}\right]}\right)^{H}{\bf B}^{-1} \left[{\matrix{{{({{\bf D}^H{\bf P}_T^\bot {\bf D}})}^{-1}{\bf D}^H{\bf P}_T^\bot {\bf T}} \cr {{( {{\bf T}^H{\bf P}_D^\bot {\bf T}} ) }^{{-}1}{\bf T}^H{\bf P}_D^\bot {\bf T}} \cr }}\right] \cr & \quad \quad \quad \quad \quad \quad \quad \quad = {\bf D}({{\bf D}^H{\bf R}_{DT}^ + {\bf D}} ) ^{{-}1} \left({ \left[{ \matrix{ {\bf I} \cr 0 \cr }}\right]}\right)^H {\bf B}^{{-}1} \left[{ \matrix{ 0 \cr {\bf I} \cr } } \right] \cr & = {\bf D}({{\bf D}^H{\bf R}_{DT}^ + {\bf D}} )^{{-}1} \left({ \left[{ \matrix{ {\bf I} \cr 0 \cr } } \right]} \right)^H \left[{ \matrix{{{( {R_D} ) }^{{-}1}} & 0 \cr 0 & {{( {R_T} ) }^{{-}1}} \cr } } \right] \left[{ \matrix{ 0 \cr {\bf I} \cr } } \right] = 0}$$

As shown in (32), the target signal is completely removed. However, in the actual working process of the SDPPR, it will inevitably be affected by the complex and variable electromagnetic environment. The electromagnetic wave received by the radar will exhibit a partially polarized wave characteristic. As shown in (1), clutter will cause random modulation of amplitude on the direct-path signal [Reference Yi, Wan, Yi and Cao18].

The composite clutter signal that contains both direct-path signal and random scattered clutter signals is regarded as a random amplitude modulation [Reference Yi, Wan, Yi and Cao18, Reference Behrens and Scharf20]. In the presence of clutter, the reception signal model is shown as below:

(33)$$\eqalign{{\bf S} = & \left[{\matrix{ {\cos \gamma_d} \cr {\sin \gamma_de^{\,j\varphi_d}} \cr } } \right]c_1d( t) + \sum\limits_{n = 2}^{M_C} {\left[{\matrix{ {\cos \gamma_c^n } \cr {\sin \gamma_c^n e^{\,j\varphi_c^n }} \cr } } \right]c_nd( t-\tau _n) } \cr & + \sum\limits_{m = 1}^{M_t} {\left[{\matrix{ {\cos \gamma_t^m } \cr {\sin \gamma_t^m e^{\,j\varphi_t^m }} \cr } } \right]\alpha _md( t-\tau _t) e^{\,j2\pi f^Dt}} + \left[{\matrix{ {w_h} \cr {w_v} \cr } } \right]} $$

where γ d and ϕ d represent polarization angle and polarization phase difference of the direct-path signal, respectively. $\gamma _c^n$ and $\varphi _c^n$ represent polarization angle and polarization phase difference of clutter, respectively. $\gamma _t^m$ and $\varphi _t^m$ represent polarization angle and polarization phase difference of target, respectively.

According to (4), we rewrite (33) and get

(34)$$\eqalign{{\bf S} = & \left[{\matrix{ {( {1 + \beta_h} ) \cos \gamma_d} \cr {( {1 + \beta_v} ) \sin \gamma_de^{\,j\varphi_d}} \cr } } \right]\varepsilon d( t) \cr & + \sum\limits_{m = 1}^{M_t} {\left[{\matrix{ {\cos \gamma_t^m } \cr {\sin \gamma_t^m e^{\,j\varphi_t^m }} \cr } } \right]\alpha _md( t-\tau _t) e^{\,j2\pi f^Dt}} + \left[{\matrix{ {w_h} \cr {w_v} \cr } } \right]} $$

where ε represents the amplitude of the composite clutter signal, βh represents amplitude modulation factor of H-polarized channel. βv represents channel amplitude modulation factor of V-polarized. βh and βv follow normal distribution.

In order to avoid the influence of amplitude modulation in (34), we will use sub-carrier processing technology. Convert (34) into sub-carrier expression

(35)$$\eqalign{{\bf S}_k = & \left[{\matrix{ {c_{1, h} + \sum\limits_{n = 2}^{M_c} {c_{n, h}e^{{-}j2\pi f_k\tau_n^{} }} } \cr {c_{1, v} + \sum\limits_{n = 2}^{M_c} {c_{n, v}e^{{-}j2\pi f_k\tau_n^{} }} } \cr } } \right]Q_k^T \cr & + \left[{\matrix{ {\sum\limits_{m = 1}^{M_t} {\alpha_{m, h}e^{{-}j2\pi ( {\,f_k-f_m^D } ) \tau_m^{} }} } \cr {\sum\limits_{m = 1}^{M_t} {\alpha_{m, v}e^{{-}j2\pi ( {\,f_k-f_m^D } ) \tau_m^{} }} } \cr } } \right]U_{k, m}^T + \left[{\matrix{ {W_{k, h}} \cr {W_{k, v}} \cr } } \right]} $$

where c 1,h and c 1,v represent the complex amplitude of direct-path signal in the H-polarized channel and V-polarized channel, respectively. cn ,h and cn ,v represent the complex amplitude of the n-th stationary multi-path clutters in the H-polarized channel and V-polarized channel, respectively. τn represents the delay of the n-th stationary multi-path clutters. αm ,h and αm ,v represent the complex amplitude of target in the H-polarized channel and V-polarized channel, respectively. τm and $f_m^D$ are the delay and the Doppler frequency of m-th target, respectively. Mc and Mt represent the number of multi-path clutter and target, respectively. Wk ,h and Wk ,v represents white Gaussian noise in the H-polarized channel and V-polarized channel, respectively.

$\left[{c_{1, h} + \sum\nolimits_{n = 2}^{M_c} {c_{n, h}e^{{-}j2\pi f_k\tau_n^{} }} } \right]$ and $\left[{c_{1, v} + \sum\nolimits_{n = 2}^{M_c} {c_{n, v}e^{{-}j2\pi f_k\tau_n^{} }} } \right]$ could be viewed as the composite weight of clutters in the H-polarized and V-polarized channels, and in each sub-carrier the composite weight is constant.

Considering single target, (35) may be simplified as

(36)$${\bf S}_k = {\bf D}_k^{cs} \varepsilon _kQ_k^T + {\bf T}_k^t \alpha _ke^{{-}j2\pi ( {\,f_k-f^D} ) \tau ^t}U_k^T + \left[{\matrix{ {W_{k, h}} \cr {W_{k, v}} \cr } } \right]$$

${\bf D}_k^{cs} = \left[{\matrix{ {\cos \gamma_k^c } \cr {\sin \gamma_k^c e^{j\varphi_k^c }} \cr } } \right]$ represents polarization vector of composite clutter signal in the k-th sub-carrier. $\gamma _k^c$ and $\varphi _k^c$ represent polarization angle and polarization phase difference of composite clutter signal. ${\bf T}_k^t = \left[{\matrix{ {\cos \gamma_k^t } \cr {\sin \gamma_k^t e^{j\varphi_k^t }} \cr } } \right]$ represents polarization vector of target in the k-th sub-carrier. $\gamma _k^t$ and $\varphi _k^t$represent polarization angle and polarization phase difference of target. $\varepsilon _k$ and α k represent the amplitude of the composite clutter signal and target in the k-th sub-carrier.

The blind adaptive oblique projection operator in each sub-carrier can be expressed as follows:

(37)$${\bf E}_k^{b-op} = {\bf D}_k^{cs} ( {{( {\bf D}_k^{cs} ) }^H{( {\bf R}_k^{DT} ) }^ + {\bf D}_k^{cs} } ) ^{{-}1}( {\bf D}_k^{cs} ) ^H( {\bf R}_k^{DT} ) ^ + $$
(38)$${\bf R}_k^{DT} = ( {{\bf A}_k{\bf B}_k{\bf A}_k^H } ) ^ + $$

where ${\bf A}_k = \left[{\matrix{ {{\bf D}_k^{DT} } & {{\bf T}_k^t } \cr } } \right]$ and ${\bf B}_k = \left[{\matrix{ {R_{k, D}} & 0 \cr 0 & {R_{k, T}} \cr } } \right]$.

Combining (36) and (37), we can effectively suppress the target signal:

(39)$$\eqalign{{\bf S}_k^{b-op} & = {\bf {\rm E}}_k^{b-op} {\bf S}_k = {\bf D}_k^{cs} ( {{( {\bf D}_k^{cs} ) }^H{( {\bf R}_k^{DT} ) }^ + {\bf D}_k^{cs} } ) ^{{-}1}( {\bf D}_k^{cs} ) ^H( {\bf R}_k^{DT} ) ^ + {\bf S}_k \cr& \quad= {\bf D}_k^{cs} \varepsilon _kQ_k^T + \left[{\matrix{ {\tilde{W}_k^h } \cr {\tilde{W}_k^v } \cr } } \right]} $$

where $\tilde{W}_k^h = E_k^{b-op} W_k^h$ and $\tilde{W}_k^v = E_k^{b-op} W_k^v$. As shown in (39), the composite clutter signal is extracted and the target signal is suppressed.

Reference signal purification

Performing IDFT on the data obtained in the previous step to convert the frequency domain signal into the time domain. At this point, the target signal has been suppressed, leaving only the composite signal containing the direct-path wave and zero-Doppler multi-path clutter with different time delays. In order to obtain a purer reference signal, we need to further suppress multi-path clutter. Therefore, we will study in detail the use of ECA to suppress multi-path clutter in this section.

Assume that the transmitted signal is V-polarized polarization, the energy of direct-path signal received by the V-polarized channel is much stronger than that of the H-polarized channel in general. Thus, we consider extracting the direct-path signal in the V-polarized channel and construct composite signal delay subspace matrix Cv as follows:

(40)$${\bf C}_v = \left[{\matrix{ 0 & {S_{sc}^v ( 1) } & {S_{sc}^v ( 2) } & \cdots & {S_{sc}^v ( N-1) } \cr 0 & 0 & {S_{sc}^v ( 1) } & \cdots & {S_{sc}^v ( N-2) } \cr \vdots & \vdots & \vdots & \ddots & \vdots \cr 0 & 0 & 0 & \cdots & {S_{sc}^v ( N-K + 1) } \cr } } \right]^T$$

where $S_{sc}^v ( n) \,\,\matrix{ {n = 1, \;2, \;\ldots , \;N} \cr }$ represents composite signal in the V-polarized channel. N represents the total data length. K represents the number of cancelation distance cells. In matrix Cv, the first row to the (K − 1)-th row represent delay of one range bin to (K − 1)-th range bin in sample signal, respectively.

According to the orthogonal characteristics of the subspace, the orthogonal projection operator can be obtained as

(41)$${\bf P}_{sc}^v = {\bf I}-{\bf C}_v( {{\bf C}_v^H {\bf C}_v} ) ^{{-}1}{\bf C}_v^H $$

The residual signal $S_{ECA, d}^v$ in the V-polarized channel of surveillance signal can be expressed as follows:

(42)$$S_{ECA, d}^v = [ {{\bf I}-{\bf C}_v{( {{\bf C}_v^H {\bf C}_v} ) }^{{-}1}{\bf C}_v^H } ] S_{sc}^v = {\bf P}_{sc}^v S_{sc}^v $$

where $S_{ECA, d}^v$ is the purified reference signal.

The signal processing flow of SDPPR by using the new method of reference extraction is shown in Fig. 2.

Fig. 2. Signal processing flow of combined processing method.

As shown in Fig. 2, the echo signal has been received in the surveillance channel by the single dual-polarization antenna. After the echo signal is processed by sub-carrier, blind adaptive oblique projection technology and ECA, the pure reference signal is obtained. Then, we can suppress the direct-path signal and multi-path clutters by employing signal processing method of [Reference Yi, Wan, Li and Leung32].

Error analysis

It can be seen from (24) and (32), the estimation precision of noise variance σ 2 has an effect on target signal suppression. Assume that there is error in estimation of noise variance in each sub-carrier. Then, matrix $\tilde{{\bf B}}_k$ can be expressed as

(43)$$\eqalign{{\tilde{\bf B}}_k & = \left[{\matrix{ {R_{k, D}} & 0 \cr 0 & {R_{k, T}} \cr } } \right] + \left[{\matrix{ {{\bf D}_k^{cs} } & {{\bf T}_k^t } \cr } } \right]^{{-}1} \left[{\matrix{ {\Delta \sigma^2} & 0 \cr 0 & {\Delta \sigma^2} \cr } } \right] \cr & \times \left({{\left[{\matrix{ {{\bf D}_k^{cs} } & {{\bf T}_k^t } \cr } } \right]}^H} \right)^{{-}1}}$$

where Δσ 2 represents estimation error of noise variance.

Assume that

(44)$$\left[{\matrix{ {{\bf D}_k^{cs} } & {{\bf T}_k^t } \cr } } \right]^{{-}1} = \left[{\matrix{ {\kappa_x} & {\mu_x} \cr {\kappa_y} & {\mu_y} \cr } } \right]$$
(45)$$\eqalign{{\tilde{{\bf B}}}_k = & \left[{\matrix{ {R_{k, D}} & 0 \cr 0 & {R_{k, T}} \cr } } \right] + \left[{\matrix{ {( {\kappa_x\kappa_x^\ast{ + } \mu_x\mu_x^\ast } ) \Delta \sigma^2} & {( {\kappa_x\kappa_y^\ast{ + } \mu_x\mu_y^\ast } ) \Delta \sigma^2} \cr {( {\kappa_x^\ast \kappa_y + \mu_x^\ast \mu_y} ) \Delta \sigma^2} & {( {\kappa_y\kappa_y^\ast{ + } \mu_y\mu_y^\ast } ) \Delta \sigma^2} \cr } } \right]\cr = & \left[{\matrix{ {R_{k, D} + ( {\kappa_x\kappa_x^\ast{ + } \mu_x\mu_x^\ast } ) \Delta \sigma^2} & {( {\kappa_x\kappa_y^\ast{ + } \mu_x\mu_y^\ast } ) \Delta \sigma^2} \cr {( {\kappa_x^\ast \kappa_y + \mu_x^\ast \mu_y} ) \Delta \sigma^2} & {R_{k, T} + ( {\kappa_y\kappa_y^\ast{ + } \mu_y\mu_y^\ast } ) \Delta \sigma^2} \cr } } \right]} $$

It is assumed that the polarization state estimation of composite clutter signal is unbiased estimation. Then, formula (32) can be expressed as

(46)$${\eqalign{{\bf E}_k^{b-op} {\bf T}_k^t =\;\; & {\bf D}_k^{cs} ( {{( {\bf D}_k^{cs} ) }^H{( {{\bf A}_k{\tilde{{\bf B}}}_k{\bf A}_k^H } ) }^ + {\bf D}_k^{cs} } ) ^{{-}1}\left({\left[{\matrix{ {\bf I} \cr 0 \cr } } \right]} \right)^H\quad \cr & \tilde{{\bf B}}_k^{{-}1} \left[{\matrix{ {{( {{( {\bf D}_k^{cs} ) }^H{\bf P}_T^\bot {\bf D}_k^{cs} } ) }^{{-}1}{( {\bf D}_k^{cs} ) }^H{\bf P}_T^\bot } \cr {{( {{( {\bf T}_k^t ) }^H{\bf P}_D^\bot {\bf T}_k^t } ) }^{{-}1}{( {\bf T}_k^t ) }^H{\bf P}_D^\bot } \cr } } \right]{\bf T}_k^t \cr \;=\;\; & {\bf D}_k^{cs} ( {{( {\bf D}_k^{cs} ) }^H{( {{\bf A}_k{\tilde{{\bf B}}}_k{\bf A}_k^H } ) }^ + {\bf D}_k^{cs} } ) ^{{-}1}\left({\left[{\matrix{ {\bf I} \cr 0 \cr } } \right]} \right)^H\tilde{{\bf B}}_k^{{-}1} \left[{\matrix{ 0 \cr {\bf I} \cr } } \right]\cr =\;\; & {\bf D}_k^{cs} ( {{( {\bf D}_k^{cs} ) }^H{( {{\bf A}_k{\tilde{{\bf B}}}_k{\bf A}_k^H } ) }^ + {\bf D}_k^{cs} } ) ^{{-}1}\left({\left[{\matrix{ {\bf I} \cr 0 \cr } } \right]} \right)^H\quad\cr & \left[{\matrix{ {R_{k, D} + ( {\kappa_x\kappa_x^\ast{ + } \mu_x\mu_x^\ast } ) \Delta \sigma^2} & {( {\kappa_x\kappa_y^\ast{ + } \mu_x\mu_y^\ast } ) \Delta \sigma^2} \cr {( {\kappa_x^\ast \kappa_y + \mu_x^\ast \mu_y} ) \Delta \sigma^2} & {R_{k, T} + ( {\kappa_y\kappa_y^\ast{ + } \mu_y\mu_y^\ast } ) \Delta \sigma^2} \cr } } \right]\left[{\matrix{ 0 \cr {\bf I} \cr } } \right]\cr =\;\; & {\bf D}_k^{cs} \displaystyle{{( {\kappa_x\kappa_y^\ast{ + } \mu_x\mu_y^\ast } ) \Delta \sigma ^2} \over {R_{k, D} + ( {\kappa_x\kappa_x^\ast{ + } \mu_x\mu_x^\ast } ) \Delta \sigma ^2}} \cr & =\;\; \left[{\matrix{ {\cos \gamma_{_k }^c } \cr {\sin \gamma_{_k }^c \cdot e^{\,j\varphi_{_k }^c }} \cr } } \right]\displaystyle{{-\sin \gamma _{_k }^t \sin \gamma _{_k }^c e^{\,j( {\varphi_{_k }^t -\varphi_{_k }^c } ) }-\cos \gamma _k^t \cos \gamma _k^c } \over {MM^\ast ( R_D/\Delta \sigma ^2) + 1}}}} $$

where $M = \cos \gamma _k^t \sin \gamma _k^c e^{j\varphi _{_k }^c }-\cos \gamma _k^c \sin \gamma _k^t e^{j\varphi _k^t }$.

The suppression performance of target signal depends on the ratio between the composite clutter signal energy (Rk ,D) and the estimation error of noise variance (Δσ 2) (CENR). The simulation parameters are shown in Table 1.

Table 1. Simulation parameters

The relationship between target signal residual energy and CENR is shown in Fig. 3. With the increase of CENR in both H-polarized channel and V-polarized channel, the target residual energy gradually decreases.

Fig. 3. Relationship between target signal residual energy and CENR.

Experimental research

In the flow of signal processing of SDPPR, the extraction of reference signal is an important step. The reference signal can improve the SNR of target in the surveillance channel through matched filtering technology and time-domain clutter suppression method. In this section, we will verify the effectiveness of the proposed method from an experimental perspective.

Experimental scheme

The method proposed in this paper is an effective alternative when the SDPPR system RSR fails. In order to compare the performance of the proposed method, a reference antenna was added on the SDPPR system.

We will obtain the reference signal in two ways. (1) The reference signal is extracted from the reference channel using RSR method. The direct-path signal with high SNR can be obtained by using reference antenna point to the transmitting station. (2) Reference signal extraction in SDPPR based on multi-processing (RSMP) method. When the direct-path signal SNR of SDPPR system surveillance channel is low, the reference signal cannot be reconstructed, and the radar system fails. The combined processing method proposed in this paper can effectively extract the reference signal and restore the working state of the radar system.

Experimental setup

The experimental system consists of reference antenna, surveillance antenna, multi-channel receiver, and upper computer. The surveillance antenna [Reference Cao, Gong, Yi, Wang and Wan33] is shown in Fig. 4. The whole antenna system is composed of one pair of dual-polarization antenna and one pair of single polarization antenna. Dual-polarization antenna is called surveillance antenna and used to receive surveillance signal. Single polarization antenna is used to receive reference signal and called reference antenna. The polarization isolation of surveillance antenna in the frequency band (710–718 MHz) is >20 dB. The gain of surveillance antenna is ~10 dB.

Fig. 4. Antenna system.

As shown in Fig. 5, the green sector area represents surveillance antenna beam direction. The blue sector area represents reference antenna beam direction. The yellow elliptical area represents ground and buildings which bring clutter. The detection target was Boeing 737 plane.

Fig. 5. Experimental geometry.

Experimental analysis

The SDPPR system needs to reconstruct the reference signal from the surveillance channel. In order to verify the effectiveness of the method proposed in this paper, the surveillance antenna of the SDPPR system was pointed at an angle of 120° with the transmitter, so as to weaken the influence of direct-path signal on the surveillance signal. In this scenario, we can obtain the surveillance signal spectrum as shown in Fig. 6(a), and the constellation diagram on surveillance signal as shown in Fig. 6(b). When using the reference antenna to collect data, the received signal spectrum is shown in Fig. 7(a), and the signal constellation diagram is shown in Fig. 7(b).

Fig. 6. Surveillance channel data: (a) signal frequency spectrum and (b) constellation diagram.

Fig. 7. Reference channel data: (a) signal frequency spectrum and (b) constellation diagram.

In the effective bandwidth, the signal spectrum amplitude of the reference channel is at least 10 dB stronger than that of the surveillance channel, which also means that the SNR of the direct-path signal in the reference channel is much stronger than that in the surveillance channel. Therefore, the reference channel signal constellation is more focused than the surveillance channel constellation. In this scenario, the reference signal cannot be reconstructed from the surveillance channel.

The experimental data were collected from a Boeing 737 in the radar surveillance area. After obtaining the reference signal, we perform time-domain filtering and matched filtering on the V-polarized channel and H-polarized channel data, respectively. When the reference signal is obtained from the reference channel using signal reconstruction, the range Doppler (RD) map of the H-polarized channel and the V-polarized channel are shown in Fig. 8. When using the method in this paper to obtain the reference signal from the surveillance channel, the RD map of the H-polarized channel and the V-polarized channel are shown in Fig. 9.

Fig. 8. RD maps in the 80-th frame in RSR: (a) RD map of V channel and (b) RD map of H channel.

Fig. 9. RD maps in the 80-th frame in RSMP: (a) RD map in V channel and (b) RD map in H channel.

It can be seen from Figs 8 and 9, whether it is the H-polarized channel or the V-polarized channel, the reference signal obtained based on the reconstruction method has a better clutter suppression effect. This is because the SNR of the reference signal extracted from the surveillance channel is lower than that obtained from the reference channel reconstruction. However, the target peak is still evident in Fig. 9.

The RD map data of the H-polarized and V-polarized channels were fused [Reference Yi, Zhu and Cao34], and then the cell-average constant false alarm rate detection was performed, and the detection results are shown in Fig. 10. The red dots represent the target detection points. The green line represents the target real track obtained from the automatic-dependent surveillance broadcast system.

Fig. 10. Result of target detection: (a) detection result in RSR and (b) detection result in RSMP.

The detection result by using the RSR method as shown in Fig. 10(a). The target detection points are very continuous between the 50th and 100th seconds. The detection result by using the RSMP method as shown in Fig. 10(b). Although the reference signal cannot be reconstructed from the surveillance channel, the detection points obtained by using the RSMP method are still satisfactory.

In the scenario set in this paper, the statistics of the target detection results are shown in Table 2.

Table 2. Target detection statistical result

In Table 2, we can clearly see that the best result is obtained by using the reference signal reconstruction from reference channel (RSR-RC) method, followed by the reference signal extraction in multi-processing from surveillance channel (RSMP-SC) method. The reference signal reconstruction from surveillance channel (RSR-SC) method cannot obtain reference signal, and the target detection fails. For the SDPPR system, there are only one surveillance antenna and no reference antenna, so the conventional method to obtain reference signal is signal reconstruction. However, in urban environments, due to factors such as terrain and interference, signal reconstruction may fail. The method proposed in this paper can solve the above problems and is an effective alternative.

Conclusion

In this paper, a combined reference signal extraction method of SDPPR system based on sub-carrier processing, blind adaptive oblique projection technology and ECA is studied. This method is different from the signal reconstruction method. When the SNR of the direct-path signal is low and the reference signal cannot be reconstructed, the reference signal can still be extracted in the RSMP method. This can greatly improve the adaptability of the SDPPR system in complex urban environment scenarios. Therefore, the RSMP method can be used as an effective supplement to the RSR method.

Acknowledgements

The authors acknowledge Wuhan University for supporting the work of this paper. This work was supported by the National Nature Science Foundation of China (61967007, 61963016), the Key Research and Development Program of Jiangxi Province (20202BBEL53014, 20201BBF61012), Natural Science Foundation of Jiangxi Province (20212BAB202005), the Research Project of National Defence Science and Technology Key Laboratory Foundation (6142113180101), and the Jiangxi Province Department of Education Science and Technology Research Key Project (GJJ200667, GJJ170360).

Conflict of interest

The authors report no conflict of interest.

Yucheng Yi was born in Hunan, China, in 1985. He received his Ph.D. degree from the School of Electronic Information, Wuhan University. He is now a lecturer in the School of Information Engineering, East China Jiaotong University. His main research interests include radar signal processing and polarization array signal processing.

Lu Zhu was born in Jiangxi, China, in 1976. He received his Ph.D. degree from the Huazhong University of Science and Technology. He is now a professor and master supervisor in the School of Information Engineering, East China Jiaotong University. Recently, he has hosted and participated in more than 10 national research projects.

Long Liang was born in Shandong, China, in 1985. He received his Ph.D. degree from the School of Electronic Information, Wuhan University. His main research interests focus on airborne passive radar.

Yan Fu was born in Hubei, China, in 1989. He received his Ph.D. degree from the School of Electronic Information, Wuhan University. His main research interests include interference suppression, radar detection, and FM radar.

References

Colone, F, O'Hagan, DW, Lombardo, P and Baker, CJ (2009) A multistage processing algorithm for disturbance removal and target detection in passive bistatic radar. IEEE Transactions on Aerospace and Electronic Systems 45, 698722.CrossRefGoogle Scholar
Wan, XR, Yi, JX, Zhao, ZX and Ke, HY (2014) Experimental research for CMMB-based passive radar under a multipath environment. IEEE Transactions on Aerospace and Electronic Systems 50, 7085.CrossRefGoogle Scholar
Malanowski, M, Kulpa, K and Misiurewicz, J (2008) PaRaDe: passive radar demonstrator family development at Warsaw University of Technology. IEEE Microwaves Radar and Remote Sensing Symposium, Kiev, Ukraine, 22–24 September, pp. 75–78.CrossRefGoogle Scholar
Poullin, D (2005) Passive detection using digital broadcasters (DAB, DVB) with COFDM modulation. IEE Proceedings on Radar Sonar and Navigation 152, 143152.CrossRefGoogle Scholar
Coleman, C and Yardley, H (2008) Passive bistatic radar based on target illuminations by digital audio broadcasting. IET Proceedings on Radar Sonar and Navigation 2, 366375.CrossRefGoogle Scholar
Colone, F, Woodbridge, K, Guo, H, Mason, D and Baker, CJ (2011) Ambiguity function analysis of wireless LAN transmissions for passive radar. IEEE Transactions on Aerospace and Electronic Systems 47, 240264.CrossRefGoogle Scholar
Tan, DKP, Sun, H, Lu, Y, Lesturgle, M and Chan, HL (2005) Passive radar using global system for mobile communication signal: theory, implementation and measurements. IEE Proceedings on Radar Sonar and Navigation 152, 116123.CrossRefGoogle Scholar
Bartoletti, S, Conti, A and Win, MZ (2014) Passive radar via LTE signal of opportunity. IEEE International Conference on Communications Workshops, Sydney, NSW, Australia, 10–14 June, pp. 1–5.CrossRefGoogle Scholar
Zhao, JH, Ni, S, Yang, L, Zhang, Z, Gong, Y and Yu, X (2019) Multiband cooperation for 5G HetNets: a promising network paradigm. IEEE Vehicular Technology Magazine 14, 8593.CrossRefGoogle Scholar
Zhao, JH, Liu, J, Yang, L, Ai, B and Ni, S (2021) Future 5G-oriented system for urban rail transit: opportunities and challenges. China Communication 18, 112.CrossRefGoogle Scholar
Colone, F, Cardinali, R, Lombardo, P, Crognale, O, Cosmi, A, Lauri, A and Bucciarelli, T (2009) Space-time constant modulus algorithm for multipath removal on the reference signal exploited by passive bistatic radar. IET Radar Sonar and Navigation 3, 253264.CrossRefGoogle Scholar
Palmer, JE, Harms, HA, Searle, SJ and Davis, LM (2013) DVB-T passive radar signal processing. IEEE Transactions on Signal Processing 61, 21162126.CrossRefGoogle Scholar
Baczyk, MK and Malanowski, M (2010) Decoding and reconstruction of reference DVB-T signal in passive radar systems. IEEE International Radar Symposium (IRS). Vilnius, Lithuania, 16–18 June, pp. 1–4.Google Scholar
O'Hagan, DW, Kuschel, H, Heckenbach, J, Ummenhofer, M and Schell, J (2010) Signal reconstruction as an effective means of detecting targets in a DAB-based PBR. IEEE International Radar Symposium (IRS). Vilnius, Lithuania, 16–18 June, pp. 1–4.Google Scholar
Zemmari, R (2010) Reference signal extraction for GSM passive coherent location. IEEE International Radar Symposium (IRS). Vilnius, Lithuania, 16–18 June, pp. 1–4.Google Scholar
Wan, XR, Wang, JF, Hong, S and Tang, H (2011) Reconstruction of reference signal for DTMB-based passive radar systems. IEEE CIE International Conference on Radar. Chengdu, China, 24–27 October, pp. 165–168.Google Scholar
Searle, S, Davis, L and Palmer, J (2015) Signal processing considerations for passive radar with a single receiver. IEEE International Conference on Acoustics, Brisbane, Australia, 19–24 April, pp. 5560–5564.CrossRefGoogle Scholar
Yi, YC, Wan, XR, Yi, JX and Cao, XM (2019)Polarization diversity technology research in passive radar based on subcarrier processing. IEEE Sensors Journal 19, 17101719.CrossRefGoogle Scholar
Yi, YC, Wan, XR, Yi, JX and Cao, XM (2018) Polarization experimental research of passive radar based on digital television signal. Electronics Letters 54, 385387.CrossRefGoogle Scholar
Behrens, RT and Scharf, LL (1994) Signal processing applications of oblique projection operators. IEEE Transactions on Signal Processing 42, 14131424.CrossRefGoogle Scholar
Boyer, R and Bouleux, G (2008) Oblique projections for direction of arrival estimation with prior knowledge. IEEE Transactions on Signal Processing 56, 13741387.CrossRefGoogle Scholar
Mao, XP, Liu, AJ, Hou, HJ, Hong, H, Guo, R and Deng, WB (2012) Oblique projection polarisation filtering for interference suppression in high-frequency surface wave radar. IET Radar Sonar and Navigation 6, 7180.CrossRefGoogle Scholar
Mccloud, ML and Scharf, LL (2002) A new subspace identification algorithm for high-resolution DOA estimation. IEEE Transactions on Antennas and Propagation 50, 13821390.CrossRefGoogle Scholar
Cao, B, Zhang, QY, Liang, D, Wen, SM, Jin, L and Zhang, YQ (2010) Blind adaptive polarization filtering based on oblique projection. IEEE International Conference on Communications, Cape Town, South Africa, 23–27 May, pp. 1–5.CrossRefGoogle Scholar
Hou, HJ, Mao, XP and Li, SB (2012) A generalized oblique projection operator for interference suppression under colored noise. IEEE Radar Conference, Atlanta, USA, 7–11 May, pp. 687–692.Google Scholar
Yi, CL, Ji, ZY, Kirubarajan, T, Xie, JH and Hu, B (2016) An improved oblique projection method for sea clutter suppression in shipborne HFSWR. IEEE Geoscience and Remote Sensing Letters 13, 10891093.CrossRefGoogle Scholar
Tannaka, A and Imai, H (2014) Music-based DOA estimation by oblique projection along the signal subspace. IEEE Workshop on Statistical Signal Processing, Gold Coast, Australia, 2 July–29 June, pp. 300–303.CrossRefGoogle Scholar
Sun, ML, Xie, JH, Hao, ZW and Yi, CL (2012) Target detection and estimation for shipborne HFSWR based on oblique projection. IEEE International Conference on Signal Processing, Beijing, China, 21–25 October, pp. 386–389.CrossRefGoogle Scholar
Wang, J, Zhang, QY and Cao, B (2009) Multi-notch polarization filtering based on oblique projection. IEEE Global Mobile Congress, Shanghai, China, 12–14 October, pp. 1–6.CrossRefGoogle Scholar
Coon, J, Sandell, M, Beach, M and Mcgeehan, J (2006) Channel and noise variance estimation and tracking algorithms for unique-word based single-carrier systems. IEEE Transactions on Wireless Communication 5, 14881496.CrossRefGoogle Scholar
Boumard, S (2003) Novel noise variance and SNR estimation algorithm for wireless MIMO OFDM systems. IEEE Global Telecommunication Conference, San Francisco, USA, 1–5 December, pp. 1330–1334.CrossRefGoogle Scholar
Yi, JX, Wan, XR, Li, DS and Leung, H (2018) Robust clutter rejection in passive radar via generalized subband cancellation. IEEE Transactions on Aerospace and Electronic Systems 54, 19311946.CrossRefGoogle Scholar
Cao, XM, Gong, ZP, Yi, YC, Wang, BJ and Wan, XR (2016) Design of a dual-polarized Yagi-Uda antenna for the passive radar. IEEE Conference Publications. Guilin, China, 18–21 October, pp. 125–128.CrossRefGoogle Scholar
Yi, YC, Zhu, L and Cao, XM (2022) Signal fusion research in passive radar based on polarization diversity technology. International Journal of Microwave and Wireless Technologies, 114.Google Scholar
Figure 0

Fig. 1. Signal reconstruction process in passive radar.

Figure 1

Fig. 2. Signal processing flow of combined processing method.

Figure 2

Table 1. Simulation parameters

Figure 3

Fig. 3. Relationship between target signal residual energy and CENR.

Figure 4

Fig. 4. Antenna system.

Figure 5

Fig. 5. Experimental geometry.

Figure 6

Fig. 6. Surveillance channel data: (a) signal frequency spectrum and (b) constellation diagram.

Figure 7

Fig. 7. Reference channel data: (a) signal frequency spectrum and (b) constellation diagram.

Figure 8

Fig. 8. RD maps in the 80-th frame in RSR: (a) RD map of V channel and (b) RD map of H channel.

Figure 9

Fig. 9. RD maps in the 80-th frame in RSMP: (a) RD map in V channel and (b) RD map in H channel.

Figure 10

Fig. 10. Result of target detection: (a) detection result in RSR and (b) detection result in RSMP.

Figure 11

Table 2. Target detection statistical result