Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Spatiotemporal Variation and Stability of Rice Planting Using Landsat–MODIS Fusion Images from 1990 to 2020
Next Article in Special Issue
Analysis of Electromagnetic Wave and Multipath Suppression from Overhead Perspective
Previous Article in Journal
A Case Study on the Convection Initiation Mechanisms over the Northern Edge of Tarim Basin, Xinjiang, Northwest China
Previous Article in Special Issue
HCM-LMB Filter: Pedestrian Number Estimation with Millimeter-Wave Radar in Closed Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Range Sidelobe Iterative Suppression Algorithm for Extended Target with Non-Grid Multiple Scattering Points

1
The School of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210000, China
2
The 723 Institute of CSSC, Zhenjiang 212000, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(19), 4811; https://doi.org/10.3390/rs15194811
Submission received: 26 June 2023 / Revised: 17 September 2023 / Accepted: 27 September 2023 / Published: 3 October 2023
(This article belongs to the Special Issue Advances in Radar Systems for Target Detection and Tracking)

Abstract

:
This paper concentrates on extended targets to suppress high sidelobes in non-grid multiple scattering points. Notice that the problem essentially resolves the high sidelobes of the non-grid multiple scatterers after locating positions; we proposed a novel iterative algorithm based on an adaptive pulse compression algorithm on the basis of an off-grid position estimation. First, the grid offsets of the multiple scattering points through decoherence and super-resolution techniques are estimated by employing an enhanced MUSIC algorithm. Then, an optimization model is constructed for the non-grid multiple scattering points of the target. Finally, a sidelobe suppression algorithm based on the minimum mean square error (MMSE) criterion is presented. Numerical results reveal that the proposed algorithm can achieve better estimation performance on the offsets of the non-grid multiple scattering points, while it suppresses range sidelobes effectively.

1. Introduction

Matched filter (MF) is widely applied in pulse compression to solve the contradiction between target detection performance and high-range resolution. However, MF can produce high-range sidelobes. Since a weak target is easily covered by the sidelobes of a strong target, it results in the degradation of target detection performance [1]. Therefore, extensive research has been focused on the sidelobe suppression algorithm in recent years.
It is common to employ the windowing method in sidelobe suppression. However, the windowing algorithm brings about the broadening width of the main lobe and the loss of signal-to-noise ratio (SNR) [2,3,4]. Mismatch filter (MMF) can also improve the detection ability of weak targets and the ratio of main lobe and sidelobe, but MMF will also increase the main lobe width and introduce SNR loss [5,6,7]. The CLEAN algorithm proposed in [8] needs to search the locations of the strong targets and then eliminate the influence of the range sidelobes. It can suppress the range sidelobes of the strong targets, while the extraction effect of weak targets is not satisfactory. The least square (LS) method proposed in [9] can acquire better sidelobe suppression performance. However, the processing of scattering points outside the window causes a serious impact on the processing of targets in the window [10,11].
In [12], an adaptive pulse compression (APC) algorithm, which takes the matched filter output as prior information, is proposed by Shannon D. Blunt and Karl Gerlach according to the minimum mean square error (MMSE) criterion. In the APC algorithm, adaptive filters are iteratively calculated with the power value of the adjacent range unit and can effectively suppress range sidelobes. In [13,14], the authors compare the APC algorithm with the mismatched filter and LS method, and prove that the APC algorithm has remarkable sidelobes suppression effect. In addition, considering the Doppler mismatch existing in the target echo, Blunt proposed an APC method based on Doppler compensation to avoid the performance degradation due to Doppler mismatch [15].The adaptive pulse compression repair (PCR) method is proposed in [16], and it performs APC on the result of MF and has stronger robustness against Doppler mismatch and systematic errors. In order to further reduce the calculation amount of the APC algorithm, Blunt combined a dimension reduction algorithm and proposed the fast adaptive pulse compression (FAPC) method, which includes DFAPC and CFAPC [17,18]. FAPC preserves most performance advantages of APC while reducing the computational effort, but reduces the degrees of freedom of the algorithm and can suppress weak targets in dense targets’ environments.
The above APC algorithms assume that the targets are located on the sampling points. Since the targets are randomly distributed with different ranges, the targets will not strictly fall on the sampling points, resulting in sampling mismatch. In APC processing, when targets are not located at the sampling points, the APC result will have great performance loss. In [19,20], Henke proposed to oversample the echo signal to suppress the influence of distance sampling mismatch, but this approach will lead to an increase in computing costs. In [21,22,23], an APC algorithm based on linear constraint minimum variance (LCMV) criterion is proposed to solve the distance sampling mismatch problem by setting the main lobe’s width and interference zero-point’s constraint conditions. However, it is difficult to classify relative values of strong and weak targets and conduct quantitative operations.
Based on the above analysis, an iterative suppression algorithm for range sidelobe of the targets with non-grid multiple scattering points is considered in this paper. In our method, an improved MUSIC algorithm is used to estimate the location of multiple scattering points by decoherence and super-resolution [24,25,26,27], and an MMSE filter for multiple scattering points is derived to realize the suppression of high sidelobe. The innovation lies in that the MUSIC algorithm originally used for angle measurement can estimate the accurate position of the non-grid targets by constructing the subdivided steering vector after adding a step of dechirp [28,29,30]. On the basis of obtaining the accurate position, the original grid APC algorithm is improved to the non-grid APC algorithm, which can restrain the high-distance sidelobe. In this way, it avoids covering the small targets and completes the demand of sidelobe suppression. Numerical results verify the effectiveness of the proposed algorithm. The simulation parameters are set to refer to the parameters of the actual maritime detection radar. Parameters of some maritime detection radars are tested data. These echo signals have dense targets and wide sidelobes, and the strong target covers the weak targets, so the algorithm can completely realize the function of detecting the signal.
The structure of this paper is organized as follows. In Section 2, the echo signal model of the target with non-grid multiple scattering points is constructed. Section 3 introduces the range sidelobe iterative suppression algorithm based on the non-grid multiple scattering points model, including the improved MUSIC algorithm and the target sidelobe iterative suppression algorithm based on the non-grid scattering points model. Section 3 also conducts an analysis and comparison of the complexity of the range dimension sidelobe suppression algorithms based on grid targets’ echo signals and non-grid targets’ echo signals. In Section 4, the proposed algorithm is verified by the simulated and practical measured data, and the algorithm results are provided and analyzed. Section 5 provides conclusions.

2. Model

Target Echo Signal Model of Non-Grid Multiple Scattering Points

The non-grid multiple scattering points model established in this chapter is shown below. The echo signal with unknown targets is constructed by a main scattering point and multiple weak scattering points. Simultaneously, the scattering points are positioned on the non-grid range points. In this model, the strong point sidelobe will mask the weak point signal.
Let the radar transmit pulse be s 0 t and it can be rewritten as follows:
s 0 ( t ) = rect t T p exp j π K t 2
where K is the frequency modulation slope, T p is the pulse width, rect t = 1 ,   t 1 2 0 ,   else   .
The echo signal y ( t ) of the multiple scattering points can be expressed as follows:
y ( t ) = q = 0 Q 1 α q rect t τ q T p exp j π K t τ q 2 + n t = q = 0 Q 1 α q s 0 t τ q + n t
where the number of strong and weak scattering points is set as Q , in which there are Q 1 non-grid and Q 2 on-grid scattering points, τ 0 is the time delay of the strongest point, τ q q = 0 , 1 , , Q 1 is the time delay of the rest points, n t is the additive white noise.
In radar signal processing process, the echo signal y ( t ) is processed by MF, which takes the convolution of the echo signal with complex conjugation of the transmitted signal. The pulse compression filter output y MF t is as follows:
y MF t = y t s 0 * t = l = 0 Q 1 B T p sin π B t τ q π B t τ q
where * represents the convolution operation, * represents the complex conjugate operation, B = K T is the bandwidth of the transmitted signal and determines the range resolution of the radar system.
In this model, it is assumed that the strong scattering point cannot be distinguished by the range resolution with other weak scattering points. At the same time, strong scattering point sidelobes can also mask weak scattering points.
Note that the scattering points are usually positioned on the non-grid range points, and sidelobes are difficult to suppress by the classical APC algorithms.

3. Methods

When dealing with dense targets, sidelobes of the multiple targets are a challenge for effectively mitigating pulse compression. It becomes imperative to conduct an analysis of the echo and accurately ascertain the locations of non-grid scattering points of the target. Consequently, the suppression of high sidelobes can be achieved based on these estimated positions.
In this section, the target model with non-grid multiple scattering points is constructed. Section 3.1 proposes the modified MUSIC algorithm to estimate the time-delay offsets of non-grid multiple scattering points in echo signal. Section 3.2 constructs the non-grid multiple scattering points model according to the delay offset of the targets. The targets’ sidelobes are suppressed by employing the distance sidelobe suppression algorithm based on MMSE criterion.

3.1. Scattering Points Offset Estimation Algorithm Based on Modified MUSIC

In the non-grid multiple scattering points model proposed in this paper, the scattering points are located on the non-grid sampling points. Traditional estimation algorithms, such as the maximum likelihood (ML) and window estimation algorithm, are inadequate for accurately estimating and locating non-grid scattering points. Compared with the above two algorithms, the MUSIC algorithm can measure the offsets of multiple scattering points simultaneously with higher accuracy and resolution and has stable and excellent performance.
Classical MUSIC algorithm is based on the principle that echo signal subspace and noise signal subspace are orthogonal when the covariance matrix is full rank. However, the echo signals in this model are coherent. Hence, it is necessary to preprocess the coherent signal using the spatial smoothing technique to modify the MUSIC algorithm for obtaining the correct estimation of the coherent model. The modified MUSIC algorithm ensures the accuracy of the estimation.
In practical applications, the value of τ l cannot be calculated by direct measurement. Since the linear frequency modulation signal is used in this paper, multiple scattering points superimpose echo signals with different offsets, which makes direct use of MUSIC algorithm unfeasible. Therefore, Dechirp processing is required for echo signals. The single frequency signal corresponding to the multiple scattering points is obtained and then processed by the MUSIC algorithm. The Dechirp reference signal is constructed based on the transmitted signal, s ref t = exp j π K t τ ref 2 , where τ ref is the reference offset.
To obtain the signal after Dechirp processing, the echo signal is multiplied by the reference signal as follows:
s rd t = y t × s ref t = q = 0 Q 1 s q t , τ q = q = 0 Q 1 rect t τ q T p exp j 2 π K τ ref τ q t exp j π K τ q 2 τ ref 2
We can find that the signal in Equation (4) is processed by multiplying a rectangular window. The length and position of the window are determined based on the length and position of the common part of the multiple scattering points echo signal. The result of multiplying a window is equivalent to truncating the signal containing the q scattering points and the peak point of the signal s ˜ q t , τ q . This can be expressed as follows:
s ˜ t = rect t t 0 T w s rd t = rect t t 0 T w q = 0 Q 1 s q t , τ q = q = 0 Q 1 s ˜ q t , τ q s ˜ q t , τ q = α q rect t t 0 T w rect t τ q T p exp j 2 π K τ ref τ q t exp j π K τ q 2 τ ref 2
where T w = L w / f s , L w is the signal’s length. The exp j π K τ q 2 τ ref 2 in Equation (5) is a constant term. Its single frequency is f q = K τ ref τ q , where τ q can be expressed as follows:
τ q = τ ref f q / K , q = 1 , 2 , , Q
Since the echo signal has Q scattering points, the Dechirp echo matrix can be constructed according to the τ q obtained by Equation (6), as shown in Equation (7):
G = s ˜ 0 t , τ 0 s ˜ 1 t , τ 1 s ˜ Q 1 t , τ Q 1 T
Therefore, the signal model can be constructed as follows:
X = G S + N
where X represents the echo signal vector after Dechirp processing, G corresponds to a matrix of size L × L w , which is obtained by intercepting multiple scattering points’ Dechirp echo matrix. S represents the amplitude of the scattering points, and N corresponds to the noise signal vector. According to Equation (8), the super-resolution of the MUSIC algorithm can be used to determine the specific grid offsets of the scattering points.
Therefore, the covariance matrix of the signal after Dechirp can be written as follows:
R s = E X X H
The echo signals of the multiple scattering points are coherent. In order to obtain correct estimation of the coherent model, the spatial smoothing technique is used to preprocess the coherent signals. The matrix R xx is as follows:
R xx = 1 d j = 1 d R s j
where d represents the number of submatrices, R s j corresponds to the submatrix intercepted by R s , and R xx represents the covariance matrix obtained after spatial smoothing and decoherence operation. Singular value decomposition (SVD) is then performed on the covariance matrix R xx , which is as follows:
R xx = U S V T
where the decomposition of R xx results in an eigenvalue matrix, represented by S = diag λ 1 , λ 2 , λ M , and an eigenvector matrix, represented by U = u 1 , u 2 , u M .
R xx is a semi-positive definite matrix with M eigenvalues, among which Q eigenvalues are positive and M Q eigenvalues are close to zero. This enables R xx to be separated into two subspaces: one subspace with Q dimension representing the signal space, and the other subspace with M Q dimension representing the noise signal space.
To obtain the precise non-grid offsets, the time-delay grid is divided near the peak point, and each grid is subdivided Z times to construct a subdivided Dechirp matrix, represented by G .
To obtain the precise non-grid offsets, intercept part of the signal with length 10 before and after the peak point, refine the time axis of this signal to Z times of the sampling rate, then construct a subdivided Dechirp matrix G .
The spatial spectral function P MUSIC can be written as follows:
P MUSIC p q = 1 a H p q U S V H a p q
The non-grid offsets subdivision point corresponding to the maximum value of P MUSIC is used to determine the position p = p 0 , p 1 , , p Q 1 . The position information of scattering points can be determined with an accuracy of up to two decimal places. The specific positions of the non-grid scattering points can be obtained through the following calculation: the offsets between the grid point positions and their non-grid point, which gives rise to specific positions P and the non-grid delay offsets Δ . These are given by the following:
Δ = Z L Z L P q = l + Δ q
where Q represents the number of scattering points. represents the rounding function, P = P 0 , P 1 , , P Q 1 represents the specific positions of scattering points, while P q represents the specific position of the non-grid scattering point after adding the offset on the grid point, which is accurate to two decimal places. The non-grid offset of the signal is denoted by Δ = Δ 0 , Δ 1 , , Δ Q 1 T , Δ q 0.5 , 0.5 . The flow chart of MUSIC algorithm is shown in Figure 1.

3.2. Target Sidelobe Iterative Suppression Algorithm Based on Non-Grid Model

The non-grid offsets Δ q ( q = 1 , 2 , , Q ) obtained by the MUSIC algorithm can be used to determine the delay τ q = l T s + Δ q of the target point at the range unit l , where T s represents the sampling time interval, let t = n T s , T s = 1 / f s , and f s represents the sampling frequency. Under this sampling condition, the signal model can be converted into a vector form to obtain the frequency modulation signal sequence s , which is as follows:
s = s 0 , s 1 , , s N 1 T
where N represents the number of sampling points within a pulse width.
The radar echo signal model of the scattering point q at the range unit l after being sampled is given by the following:
y ˜ l , Δ q = x T l , Δ q s e Δ q + v ˜ l = x T l , Δ q s ˜ Δ q + v ˜ l
where x l , Δ q = x l , Δ q , x l 1 , Δ q , , x l N , Δ q T represents N -point continuous sampling when the q th scattering point is non-grid, and v ˜ l is additive Gaussian white noise vector. represents the Hadamard product, e Δ q = e j g Δ q , e j g T s + Δ q , , e j g N 1 T s + Δ q T represents the phase mismatch between the linear frequency modulation signal and the grid point caused by Δ q , and g represents the frequency modulation function of the signal.
The pulse compression output of the matched filter after sampling at the q th scattering point is obtained by convolving the complex conjugate of the transmitted signal with the echo, and is expressed in discrete form as follows:
x MF l , Δ q = s ˜ H Δ q y ˜ l , Δ q = s ˜ H Δ q A T l , Δ q s ˜ Δ q + s ˜ H Δ q v ˜ l
where y ˜ l , Δ q = y ˜ l , Δ q , y ˜ l + 1 , Δ q , , y ˜ l + N 1 , Δ q T represents N -point continuous sampling of the echo pulse corresponding to the range unit l when sampling mismatch occurs. Its specific expression is given by the following:
y ˜ l , Δ q = A T l , Δ q s ˜ Δ q + v ˜ l
where A l , Δ q = x l , Δ q , x l + 1 , Δ q , , x l + N 1 , Δ q is the N -point continuous sampling matrix of distance dimension when the q th scattering point is non-grid, and v ˜ l = v ˜ l , v ˜ l + 1 , , v ˜ l + N 1 is the noise vector.
The previous item of the matched filter output x MF l , Δ q can be expressed as follows:
s H Δ q A T l , Δ q s Δ q = s H Δ q x T l , Δ q x T l + 1 , Δ q x T l + N 1 , Δ q s Δ q = s H Δ q s T x l , Δ q s T x l + 1 , Δ q s T x l + N 1 , Δ q = s H Δ q 0 0 s 0 Δ q s 1 Δ q 0 s 0 Δ q s 0 Δ q s 1 Δ q s N 1 Δ q s 1 Δ q 0 s N 1 Δ q s N 1 Δ q 0 T x ^ l , Δ q = s H Δ q S ˜ T Δ q x ^ l , Δ q
where x ^ l , Δ q = x l + N 1 , Δ q x l , Δ q x l N + 1 , Δ q T , let S ˜ Δ q s Δ q = r ^ Δ q , r ^ Δ q = r N + 1 Δ q r 0 Δ q r N 1 Δ q T be the autocorrelation function of the reference transmitted signal s Δ q .
Then, the range unit l matched by the q th scattering point after filtering can be expressed as follows:
x MF l , Δ q = x ^ T l , Δ q r ^ Δ q + u l
Using the pulse compression technique, most of the energy of the target concentrated in several range units allows for the application of processing windows in the construction of filter coefficients and reduces computation. The x ¯ MF l , Δ q calculated through the processing window can be expressed as follows:
x ¯ MF l , Δ q = x ˜ T l , Δ q r ˜ Δ q + u l
where x ˜ l , Δ q = x l + K d , Δ q x l , Δ q x l K r , Δ q T , r ˜ = r K r Δ q r 0 Δ q r K d Δ q T .
Then, the sampling value of K r + K d + 1 point matched by the q th scattering point can be expressed as follows:
x MF l , Δ q = B T l , Δ q r ˜ Δ q + u ˜ l
where x MF l , Δ q = x ¯ MF l K r , Δ q , x ¯ MF l K r + 1 , Δ q , x ¯ MF l + K d , Δ q , u ˜ l = u l K r u l u l + K d T .
B l , Δ q = x ˜ l K r , Δ q x ˜ l K r + 1 , Δ q Δ q x ˜ l + K d , Δ q = x l K r + K d , Δ q x l + 2 K d , Δ q x l 2 K r , Δ q x l K r + K d , Δ q
Echo signal model y ˜ l , Δ q and its pulse compression signal model x MF l , Δ q were constructed according to the non-grid offset Δ q of range unit l calculated by the above algorithm.
In the process of non-grid APC processing on the pulse compression signal, when the non-grid scattering points estimated by the MUSIC algorithm are processed, in order to develop MMSE filter on the signal model of the non-grid points, the cost function is designed by using the MF signal model.
Based on RMMSE criterion, the cost function is constructed by using the output result of matching filter, the cost function is as follows:
J l = q = 0 Q 1 E x l , Δ q ω H l x MF l , Δ q 2
where E represents the statistical expectation, and ω l is a K r + K d + 1 × 1 size of MMSE filtering coefficient vector, which is the unique calculation of signal amplitude x l for each individual non-grid scattering points cell. K = K r + K d + 1 is the length of MMSE filter.
If J l take the gradient with respect to ω H l and set it to zero, that is the following:
ω H J ω , ω H = 0
which can be calculated as follows:
ω H J ω , ω H = q = 0 Q 1 E x MF l , Δ q x l , Δ q x MF l , Δ q x MF H l , Δ q ω l = 0
it is equivalent to the following:
ω l = q = 0 Q 1 E x MF l , Δ q x MF H l , Δ q 1 E x MF l , Δ q x l , Δ q
If the impulse response of each range unit is not correlated, that is when E x m x n = 0 ,   m n , and independent from noise statistics, that is E x m v n = 0 , then we can obtain the following:
E x MF l x MF H l = q = 0 Q 1 E B T l , Δ q r ˜ Δ q r ˜ T Δ q B l , Δ q + R ˜ = q = 0 Q 1 n = 2 K r 2 K d ρ ˜ l + n r ˜ n Δ q r ˜ n T Δ q + R ˜ = q = 0 Q 1 1 n = 2 K r 2 K d ρ l + n , Δ q r ˜ n Δ q r ˜ n T Δ q + q = 0 Q 2 1 n = 2 K r 2 K d ρ l + n r ˜ n r ˜ n T + R ˜
E x MF l , Δ q x l , Δ q = q = 0 Q 1 E B T l , Δ q x l , Δ q r ˜ Δ q = q = 0 Q 1 ρ ˜ l I K r + K d + 1 r ˜ Δ q = q = 0 Q 1 1 ρ l , Δ q I K r + K d + 1 r ˜ Δ q + q = 0 Q 2 1 ρ l I K r + K d + 1 r ˜
where ρ l , Δ q = x l , Δ q 2 represents the expected power of x l , Δ q , ρ l = x l 2 represents the expected power of x l , I K r + K d + 1 is the identity matrix of K r + K d + 1 × K r + K d + 1 dimension, R ˜ = E u ˜ l u ˜ H l represents the noise covariance matrix of size K r + K d + 1 × K r + K d + 1 , so we can obtain the ω l as follows:
ω l = ρ ˜ l C l + R ˜ 1 r ˜ Δ q
C l = q = 0 Q 1 1 n = 2 K r 2 K d ρ l + n , Δ q r ˜ n Δ q r ˜ n T Δ q + q = 0 Q 2 1 n = 2 K r 2 K d ρ l + n r ˜ n r ˜ n T
where r ˜ n Δ q means to shift the signal r ˜ Δ q by n bits, leaving a blank bit to fill in the zero. For example, when n = 2 , r ˜ 2 Δ q = 0 0 r K r Δ q r K d 2 Δ q T , when n = 2 , r ˜ 2 Δ q = r K r + 2 Δ q r K d Δ q 0 0 T .
To calculate the filter coefficient ω l of the range unit l according to Equations (27) and (28), prior information of the output power of the unit before and after the grid point is required. After obtaining the filter coefficient ω l , it is conjugate transposed in order to suppress the sidelobe at the range unit l . The flow chart of the target sidelobes iterative suppression algorithm based on the non-grid model is shown in Figure 2.

3.3. Computational Complexity Analysis

The range dimension sidelobes suppression algorithm based on non-grid scattering points echo signal differs from that based on grid point echo signal by accounting for non-grid issues in practical applications. The distance dimension sidelobes suppression algorithm based on non-grid scattering points echo signal is introduced, and the reference template is adjusted according to the offsets to reduce the computational complexity of the algorithm using a small processing window.
The detailed implementation steps of these two algorithms are provided in Table 1 and Table 2, based on the discussion of the small window non-grid APC processing flow using both echo and matched filter output.
The computational complexity of these two algorithms is analyzed below in Table 3:
Table 3 displays the computational complexity of the APC algorithm based on echo and the non-grid APC algorithm based on matched filter output with small window, as analyzed below. The complexity of the APC algorithm based on echo is O N 3 , while that of the non-grid APC algorithm based on matched filter output with small window is O K 3 , where K w = K r + K d + 1 . Figure 3 compares the computational complexity of the adaptive pulse compression algorithm based on echo and matched filter. Figure 3 indicates that for N / K > 5 , the small window treatment can reduce computation by at least two orders of magnitude.
After comparing the echo-based and MF-based APC algorithms, we show that the complexity of the new algorithm is two orders of magnitude lower than that of the echo-based APC algorithm. The comparison with other algorithms is provided in the following Table 4.
From Table 4, it is evident that the proposed algorithm exhibits a lower computational complexity compared to the traditional APC algorithm, as well as FAPC and LCMV-APC. This highlights the computational superiority of the proposed algorithm.

4. Results

In both simulation and field tests, the LFM signal is used as the radar transmitting signal. In Section 3.1, multiple scattering points are set and the MUSIC algorithm is used to modify the echo signal to estimate and verify the number and positions of scattering points. In Section 3.2, the scattering points in the echo signal with an unclear number of targets are estimated and their positions are computed. The signal model is constructed to process the suppressed sidelobes of the MF signal, and the effectiveness of the algorithm proposed is verified in this section.

4.1. Simulation Result

For the simulation test, the signal parameters were set as follows: bandwidth B = 4   MHz , pulse width T p = 20   μ s , sampling frequency f s = 5   MHz , and pulse repetition period PRI which is T r = 744.8   μ s . Four points were set, and their location information is provided in Table 5.
The positions of the subdivided peak point and scattering points are calculated by inputting the echo signal into the MUSIC algorithm, which obtained the super-resolution results. Figure 4 shows the resulting super-resolution image, where each extreme point indicates the specific position of a peak point or scattering points on the non-grid. The four extreme points in Figure 4 correspond to the position offsets of P 1 , P 2 , P 3 , and P 4 , based on the setting results.
The echo signal model is reconstructed using the position information of the scattering points, and the pulse compression signal is calculated as a template. The template and the original pulse compression signal were then input into the adaptive pulse compression algorithm for computation. The results of the algorithm are presented in Figure 5, and the resulting scattering point positions are provided in Table 6.
From Figure 5 and Table 6, it can be observed that the sidelobes of the signal after applying the non-grid APC algorithm is 60 dB lower than that of the matched filter output result.
In addition, we analyze the processing advantages of the non-grid APC algorithm in the case of non-grid points. The proposed method can probe small targets and distinguish four targets clearly. The corresponding performance is shown in Figure 6 and Figure 7.
Figure 6 shows that the windowing algorithm cannot solve the non-grid multiple scatterers problem well. Moreover, it broadens the sidelobes and causes the results to deteriorate. It can be seen from Figure 7 that the adaptive sidelobe suppression technique based on the grid model cannot obtain good sidelobe suppression performance for non-grid multi-scattering targets, while the adaptive sidelobe suppression technique based on the non-grid model can suppress the processing sidelobe to the noise level. The two very near targets cannot be distinguished in the output of MF. However, they can be distinguished by the adaptive sidelobe processing based on the non-mesh model. The weak targets which are masked by the high sidelobes of the MF are also shown by the adaptive sidelobe processing technique based on the non-mesh model. It can be seen from the figure that the width of the main lobe of the pulse compression signal processed by the adaptive sidelobe suppression algorithm based on the non-grid model is significantly narrower. At the same time, the signal-to-noise ratio loss is very small.
Moreover, the close proximity of the four points causes the sidelobes of the high signal amplitude point to spread over the low signal amplitude points. However, after applying the proposed non-grid sidelobe suppression algorithm, the sidelobes of all points are effectively suppressed. Compared with counterparts, the proposed algorithm is more efficient, with low computational amount and fast processing speed.

4.2. Experimental Results of Tested Data

The tested data, which is measured by radar, have a bandwidth B = 4   MHz , a pulse width T p = 20   μ s , a sampling frequency f s = 5   MHz , and a pulse repetition period PRI which is T r = 744.8   μ s . Substituting the echo signal into the MUSIC algorithm, the resulting super-resolution image is shown in Figure 8, which was used to calculate the positions of the subdivided peak point and scattering points.
Figure 8 presents the super-resolution image, with each extreme point indicating the precise position of the peak point and scattering points on the non-grid, which are located at 938.83, 937.45, 935.11, 933.41, and 930.55. Additionally, the echo signal model is reconstructed based on the targets’ positions, and a pulse compression signal was calculated as a template. After substituting the resultant signal and the original pulse compression signal into the non-grid APC algorithm, a comparison between the non-grid APC and the original matched filter result is presented in Figure 9.
Figure 9 shows that the sidelobes are reduced by almost 20 dB with the implementation of the non-grid APC algorithm in comparison to the traditional MF method. The proposed non-grid APC algorithm enables effective suppression of the sidelobes when scattering points’ positions and peak signal coordinates are estimated well.
We apply the non-grid APC algorithm to conduct batch processing experiments on 190 groups of tested signals, and the processing results are shown in Figure 10. The figure shows that the proposed algorithm effectively suppresses sidelobe waves for two-dimensional signals.
In Figure 10, the abscissa represents the sampling points in the distance dimension, and the ordinate represents the pulse index; the sea ship targets are vertically connected into an area. As shown in Figure 10a, the highlighted yellow area represents the peak point of the scattering points, and the left and right sub-bright light yellow areas represent the weak scattering points and sidelobes. The proposed method suppresses the sidelobes effectively in the sub-bright region. Therefore, the number and characteristics of ship targets are distinct as shown in Figure 10b. As shown in Figure 10, these results show that the proposed algorithm improves the overall performance and accuracy of the radar system.

5. Conclusions

In this paper, we propose a novel range sidelobe suppression technique based on the non-grid multiple scattering points model. The improved MUSIC algorithm is utilized for time-delay super-resolution and accurate estimation of the scattering points’ position offset information. By constructing the non-grid scattering signal model using precise non-grid offsets information, an iterative algorithm is employed to design optimal filters for each range unit adaptively. This approach enhances the detection performance of dense targets on non-grid significantly, while ensuring detection accuracy and performance. Since the sidelobes of strong points are easily confounded by noise and interference, traditional pulse compression methods may not be effective in suppressing the sidelobes. Through the adaptive pulse compression method, the proposed algorithm can suppress the sidelobe heavily and improve the anti-jamming ability in complex environments. In this way, it can accurately estimate the number and position of the targets and increase the sensitivity of the radar system. Therefore, the radar system detects more distant targets. In addition, traditional pulse compression methods are prone to false positives. However, the proposed algorithm reduces the probability of false alarm significantly and improves the accuracy and reliability of target detection by suppressing the sidelobes heavily. Therefore, this algorithm solves the challenge of detecting non-grid dense points effectively, and improves the overall detection performance and accuracy.
We achieve positive results based on this study while there is still room for improvement. In complex environments, such as the presence of strong interference or signal attenuation, the suppression effect of the algorithm may be affected, resulting in inaccurate positioning. We shall delve deeper into this algorithm, aspiring to attain greater breakthroughs in the times to come.

Author Contributions

Conceptualization, J.Z.; methodology, J.Z.; software, J.X.; validation, J.X.; investigation, Y.Z.; resources, J.C.; writing—original draft preparation, J.X.; writing—review and editing, J.Z.; supervision, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation (62171220).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiangbo, L.; Hongsheng, Y.; Dongbo, H.; Minh, N.D. Patch Match Filter: Efficient Edge-Aware Filtering Meets Randomized Search for Fast Correspondence Field Estimation. In Proceedings of the 2013 IEEE Conference on Computer Vision and Pattern Recognition, Portland, OR, USA, 23–28 June 2013; pp. 1854–1861. [Google Scholar]
  2. Gerlach, K.; Blunt, S.D.; Picciolo, M.L. Robust Adaptive Matched Filtering Using the FRACTA Algorithm. IEEE Trans. Aerosp. Electron. Syst. 2004, 40, 929–945. [Google Scholar] [CrossRef]
  3. Li, K.; Yuan, L. The FIR Window Function Design Based on Evolutionary Algorithm. In Proceedings of the 2011 International Conference on Mechatronic Science, Electric Engineering and Computer (MEC), Jilin, China, 19–22 August 2011; pp. 1797–1800. [Google Scholar]
  4. Wu, X.; Liu, K.; Luo, M.; Shen, Q. LMMSE Channel Estimation Algorithm Based on Delay Time-Domain Windowing in AF Cooperative Systems. In Proceedings of the 2011 6th International ICST Conference on Communications and Networking in China (CHINACOM), Harbin, China, 17–19 August 2011; pp. 760–764. [Google Scholar]
  5. Blinchikoff, H.J. Range Sidelobe Reduction for the Quadriphase Codes. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 668–675. [Google Scholar] [CrossRef]
  6. Fang, L.; Dong, D.; Boroyevich, D.; Mattavelli, P.; Wang, S. Improving High-Frequency Performance of an Input Common Mode EMI Filter Using an Impedance-Mismatching Filter. IEEE Trans. Power Electron. 2014, 29, 5111–5115. [Google Scholar]
  7. Rao, M.N. Design of Mismatched Filters for Long Binary Codes. In Proceedings of the 2008 IEEE Radar Conference, Rome, Italy, 26–30 May 2008; pp. 1–6. [Google Scholar]
  8. Bose, R. Lean CLEAN: Deconvolution Algorithm for Radar Imaging of Contiguous Targets. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 2190–2199. [Google Scholar] [CrossRef]
  9. Zrnic, B.; Zejak, A.; Petrovic, A.; Simic, I. Range Sidelobe Suppression for Pulse Compression Radars Utilizing Modified RLS Algorithm. In Proceedings of the 1998 IEEE International Symposium on Spread Spectrum Techniques and Applications, Sun City, South Africa, 2–4 September 1998; pp. 1008–1011. [Google Scholar]
  10. Yanming, L. Study on the Technology of Quadric Surface Extracting Base on Least Square Method. In Proceedings of the 2011 Second International Conference on Mechanic Automation and Control Engineering, Hohhot, China, 15–17 July 2011; pp. 5328–5331. [Google Scholar]
  11. Jin, E.S.; Liu, L.L.; Bo, Z.Q.; Klimek, A. Parameter Identification of the Transformer Winding Based on Least-Squares Method. In Proceedings of the 2008 IEEE Power and Energy Society General Meeting-Conversion and Delivery of Electrical Energy in the 21st Century, Pittsburgh, PA, USA, 20–24 July 2008; pp. 1–6. [Google Scholar]
  12. Blunt, S.D.; Gerlach, K. Adaptive Pulse Compression. In Proceedings of the 2004 IEEE International Radar Conference, Philadelphia, PA, USA, 18–21 December 2004; pp. 271–276. [Google Scholar]
  13. Blunt, S.D.; Gerlach, K. Adaptive Pulse Compression via MMSE Estimation. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 572–584. [Google Scholar] [CrossRef]
  14. Blunt, S.D.; Gerlach, K. A Novel Pulse Compression Scheme Based on Minimum Mean-Square Error Reiteration. In Proceedings of the 2003 Proceedings of the International Conference on Radar, Adelaide, SA, Australia, 3–5 September 2003; pp. 349–353. [Google Scholar]
  15. Blunt, S.D.; Smith, K.J.; Gerlach, K. Doppler-Compensated Adaptive Pulse Compression. In Proceedings of the 2006 IEEE Conference on Radar, Verona, NY, USA, 24–27 April 2006; pp. 1–6. [Google Scholar]
  16. Gerlach, K.; Blunt, S.D. Radar Pulse Compression Repair. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1188–1195. [Google Scholar] [CrossRef]
  17. Blunt, S.D.; Higgins, T. Achieving Real-Time Efficiency for Adaptive Radar Pulse Compression. In Proceedings of the 2007 IEEE Radar Conference, Waltham, MA, USA, 17–20 April 2007; pp. 116–121. [Google Scholar]
  18. Blunt, S.D.; Higgins, T. Dimensionality Reduction Techniques for Efficient Adaptive Pulse Compression. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 349–362. [Google Scholar] [CrossRef]
  19. Henke, D.; Mccormick, P.; Blunt, S.D.; Higgins, T. Practical Aspects of Optimal Mismatch Filtering and Adaptive Pulse Compression for FM Waveforms. In Proceedings of the 2015 IEEE Radar Conference, Arlington, VA, USA, 10–15 May 2015; pp. 1149–1155. [Google Scholar]
  20. Henke, D. Robust Optimal and Adaptive Pulse Compression for FM Waveforms; University of Kansas: Kansas, MS, USA, 2015. [Google Scholar]
  21. VanTrees, H.L. Optimum Array Processing Part IV of Detection, Estimation and Modulation Theory; John Wiley and Sons Inc.: New York, NY, USA, 2002; pp. 513–518. [Google Scholar]
  22. Geng, Z.; Deng, H.; Himed, B. Adaptive Radar Beamforming for Interference Mitigation in Radar Wireless Spectrum Sharing. IEEE Signal Process. Lett. 2015, 22, 484–488. [Google Scholar] [CrossRef]
  23. Su, H.; Liu, H.; Shui, P.; Bao, Z. Adaptive Beamforming for Nonstationary HF Interference Cancellation in Skywave Over-the-Horizon Radar. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 312–324. [Google Scholar] [CrossRef]
  24. Blunt, S.D.; Gerlach, K. Adaptive Pulse Compression Repair Processing. In Proceedings of the 2005 IEEE International Radar Conference, Arlington, VA, USA, 29 July 2005; pp. 519–523. [Google Scholar]
  25. Shackelford, A.K.; De Graaf, J.; Talapatra, S.; Blunt, S.D.; Gerlach, K. Adaptive Pulse Compression: Preliminary Experimental Measurements. In Proceedings of the 2007 IEEE Radar Conference, Waltham, MA, USA, 17–20 April 2007; pp. 234–237. [Google Scholar]
  26. Baghel, V.; Panda, A.; Panda, G. An Efficient Hybrid Adaptive Pulse Compression Approach to Radar Detection. In Proceedings of the 2013 International Conference on Signal Processing and Communication, Noida, India, 12–14 December 2013; pp. 335–339. [Google Scholar]
  27. Qianrong, L.; Xiji, W.; Zhicheng, W.; Xiaochao, M.; Na, Y. Robust Adaptive Pulse Compression Under Encode-Decoder Network for Synthetic Aperture Radar. In Proceedings of the 2021 Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Bali, Indonesia, 1–5 November 2021; pp. 1–5. [Google Scholar]
  28. Qianrong, L.; Wenming, T.; Qingqing, L.; Bingqi, Z.; Ke, D.; Xiangzhen, Y. Robust Adaptive Pulse Compression for Synthetic Aperture Radar. In Proceedings of the 2021 China International SAR Symposium (CISS), Shanghai, China, 3–5 November 2021; pp. 1–6. [Google Scholar]
  29. Higgins, T.; Blunt, S.D.; Gerlach, K. Gain-Constrained Adaptive Pulse Compression via An MVDR Framework. In Proceedings of the 2009 IEEE Radar Conference, Pasadena, CA, USA, 12–16 October 2009; pp. 1–6. [Google Scholar]
  30. Akdemir, Ş.B.; Candan, Ç. Study of Angle of Arrival Estimation Performance of MVDR and Capon Methods under Intermittent Interference. In Proceedings of the 2018 Signal Processing and Communications Applications Conference (SIU), Izmir, Turkey, 2–5 May 2018; pp. 1–4. [Google Scholar]
Figure 1. Flow chart of target scattering points offset estimation algorithm based on modified MUSIC.
Figure 1. Flow chart of target scattering points offset estimation algorithm based on modified MUSIC.
Remotesensing 15 04811 g001
Figure 2. Flow chart of target sidelobe iterative suppression algorithm based on non-grid model.
Figure 2. Flow chart of target sidelobe iterative suppression algorithm based on non-grid model.
Remotesensing 15 04811 g002
Figure 3. Comparison of computational complexity of algorithm based on echo and MF.
Figure 3. Comparison of computational complexity of algorithm based on echo and MF.
Remotesensing 15 04811 g003
Figure 4. Results of location estimation by multiple points simulation MUSIC algorithm.
Figure 4. Results of location estimation by multiple points simulation MUSIC algorithm.
Remotesensing 15 04811 g004
Figure 5. Shows the results of the multi-objective simulation of the non-grid APC algorithm.
Figure 5. Shows the results of the multi-objective simulation of the non-grid APC algorithm.
Remotesensing 15 04811 g005
Figure 6. Comparison of windowing algorithm results with non-grid APC algorithm.
Figure 6. Comparison of windowing algorithm results with non-grid APC algorithm.
Remotesensing 15 04811 g006
Figure 7. Comparison between traditional grid APC and non-grid APC algorithm.
Figure 7. Comparison between traditional grid APC and non-grid APC algorithm.
Remotesensing 15 04811 g007
Figure 8. Location estimation results of MUSIC algorithm for tested data.
Figure 8. Location estimation results of MUSIC algorithm for tested data.
Remotesensing 15 04811 g008
Figure 9. Results of sidelobe suppression by non-grid APC algorithm based on tested data.
Figure 9. Results of sidelobe suppression by non-grid APC algorithm based on tested data.
Remotesensing 15 04811 g009
Figure 10. Comparison of pulse-by-pulse sidelobe suppression processing in sea surface multiple targets radar beam scanning: (a) matching processing; (b) sidelobe suppression processing.
Figure 10. Comparison of pulse-by-pulse sidelobe suppression processing in sea surface multiple targets radar beam scanning: (a) matching processing; (b) sidelobe suppression processing.
Remotesensing 15 04811 g010
Table 1. Implementation steps of APC algorithm based on echo.
Table 1. Implementation steps of APC algorithm based on echo.
  • Initialization: ω 0 = s , the matched filter output x MF was calculated by x ^ M F l = s H y ˜ l = s H A T l s + s H v l , let i = 1 ,   x ^ 0 = x MF ;
  • Calculate   power   estimates   ρ i l = x ^ i 1 l 2 , and C i l matrix is constructed by C l = n = N + 1 N 1 ρ l + n s n s n H ;
  • Calculate   ω i l   matrix   by   ω l = ρ l C l + R 1 s ;
  • Update   x ^ i l   by   x ^ l = ω H l y ˜ l ;
  • Repeat Step 2 until the desired accuracy is achieved.
Table 2. Implementation steps of small-window non-grid APC algorithm based on matched filter output.
Table 2. Implementation steps of small-window non-grid APC algorithm based on matched filter output.
  • Initialization: Calculate the matched filter output x MF , let i = 1 ,   x ^ 0 = x MF ;
  • The Δ is calculated using Equation (13), calculate power estimates ρ i l = x ^ i 1 l 2 ,   ρ i l , Δ l = x ^ i 1 l , Δ l 2 and C i l matrix is constructed by Equation (30);
  • Calculate   ω i l matrix by Equation (29);
  • Update   x ^ i l   by   x ^ l = ω ˜ l x M F l ;
  • Repeat Step 2 until the desired accuracy is achieved.
Table 3. The computational complexity of each range unit in each iteration.
Table 3. The computational complexity of each range unit in each iteration.
StepsEcho APCMatch the Filtered Non-Grid APC
C l matrix construction 2 N 1 N 2 2 K r + K d + 1 3 K r + K d + 1 2
Filter weight vector N 3 + N 2 K r + K d + 1 3 + K r + K d + 1 2
Distance   cell   x l estimation N K r + K d + 1
Table 4. Comparison of the complexity of each algorithm.
Table 4. Comparison of the complexity of each algorithm.
AlgorithmComputational Complexity
APC algorithm based on echo O N 3
FAPC O N log N
LCMV-APC O N 2 N p
Non-grid APC O K 3
Table 5. Table of simulation scenarios.
Table 5. Table of simulation scenarios.
TargetsThe Index of Sampling Points of the Target Distance DimensionTarget SNR/dB
1931.0920
2934.1840
3947.4280
4939.1175.56
Table 6. Simulation results.
Table 6. Simulation results.
TargetsThe Index of Sampling Points of the Target Distance DimensionGround Truth Amplitude/dBNon-Grid APC Amplitude/dB
1931.09−20−19.65
2934.18−40−39.83
3947.4200
4939.11−4.43−4.31
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, J.; Zhang, Y.; Zhang, J.; Chen, J. Range Sidelobe Iterative Suppression Algorithm for Extended Target with Non-Grid Multiple Scattering Points. Remote Sens. 2023, 15, 4811. https://doi.org/10.3390/rs15194811

AMA Style

Xu J, Zhang Y, Zhang J, Chen J. Range Sidelobe Iterative Suppression Algorithm for Extended Target with Non-Grid Multiple Scattering Points. Remote Sensing. 2023; 15(19):4811. https://doi.org/10.3390/rs15194811

Chicago/Turabian Style

Xu, Jing, Yanan Zhang, Jindong Zhang, and Jiarui Chen. 2023. "Range Sidelobe Iterative Suppression Algorithm for Extended Target with Non-Grid Multiple Scattering Points" Remote Sensing 15, no. 19: 4811. https://doi.org/10.3390/rs15194811

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