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.
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 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, , where is the reference offset.
To obtain the signal after Dechirp processing, the echo signal is multiplied by the reference signal as follows:
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
scattering points and the peak point of the signal
. This can be expressed as follows:
where
,
is the signal’s length. The
in Equation (5) is a constant term. Its single frequency is
, where
can be expressed as follows:
Since the echo signal has
scattering points, the Dechirp echo matrix can be constructed according to the
obtained by Equation (6), as shown in Equation (7):
Therefore, the signal model can be constructed as follows:
where
represents the echo signal vector after Dechirp processing,
corresponds to a matrix of size
, which is obtained by intercepting multiple scattering points’ Dechirp echo matrix.
represents the amplitude of the scattering points, and
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:
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
is as follows:
where
represents the number of submatrices,
corresponds to the submatrix intercepted by
, and
represents the covariance matrix obtained after spatial smoothing and decoherence operation. Singular value decomposition (SVD) is then performed on the covariance matrix
, which is as follows:
where the decomposition of
results in an eigenvalue matrix, represented by
, and an eigenvector matrix, represented by
.
is a semi-positive definite matrix with eigenvalues, among which eigenvalues are positive and eigenvalues are close to zero. This enables to be separated into two subspaces: one subspace with dimension representing the signal space, and the other subspace with 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 times to construct a subdivided Dechirp matrix, represented by .
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 .
The spatial spectral function
can be written as follows:
The non-grid offsets subdivision point corresponding to the maximum value of
is used to determine the position
. 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
and the non-grid delay offsets
. These are given by the following:
where
represents the number of scattering points.
represents the rounding function,
represents the specific positions of scattering points, while
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
. 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
obtained by the MUSIC algorithm can be used to determine the delay
of the target point at the range unit
, where
represents the sampling time interval, let
,
, and
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
, which is as follows:
where
represents the number of sampling points within a pulse width.
The radar echo signal model of the scattering point
at the range unit
after being sampled is given by the following:
where
represents
-point continuous sampling when the
th scattering point is non-grid, and
is additive Gaussian white noise vector.
represents the Hadamard product,
represents the phase mismatch between the linear frequency modulation signal and the grid point caused by
, and
represents the frequency modulation function of the signal.
The pulse compression output of the matched filter after sampling at the
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:
where
represents
-point continuous sampling of the echo pulse corresponding to the range unit
when sampling mismatch occurs. Its specific expression is given by the following:
where
is the
-point continuous sampling matrix of distance dimension when the
th scattering point is non-grid, and
is the noise vector.
The previous item of the matched filter output
can be expressed as follows:
where
, let
,
be the autocorrelation function of the reference transmitted signal
.
Then, the range unit
matched by the
th scattering point after filtering can be expressed as follows:
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
calculated through the processing window can be expressed as follows:
where
,
.
Then, the sampling value of
point matched by the
th scattering point can be expressed as follows:
where
,
.
Echo signal model and its pulse compression signal model were constructed according to the non-grid offset of range unit 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:
where
represents the statistical expectation, and
is a
size of MMSE filtering coefficient vector, which is the unique calculation of signal amplitude
for each individual non-grid scattering points cell.
is the length of MMSE filter.
If
take the gradient with respect to
and set it to zero, that is the following:
which can be calculated as follows:
it is equivalent to the following:
If the impulse response of each range unit is not correlated, that is when
, and independent from noise statistics, that is
, then we can obtain the following:
where
represents the expected power of
,
represents the expected power of
,
is the identity matrix of
dimension,
represents the noise covariance matrix of size
, so we can obtain the
as follows:
where
means to shift the signal
by
bits, leaving a blank bit to fill in the zero. For example, when
,
, when
,
.
To calculate the filter coefficient
of the range unit
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
, it is conjugate transposed in order to suppress the sidelobe at the range unit
. 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
, while that of the non-grid APC algorithm based on matched filter output with small window is
, where
.
Figure 3 compares the computational complexity of the adaptive pulse compression algorithm based on echo and matched filter.
Figure 3 indicates that for
, 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.