Next Article in Journal
Optimization and Analysis of Surface Roughness, Flank Wear and 5 Different Sensorial Data via Tool Condition Monitoring System in Turning of AISI 5140
Previous Article in Journal
Testing Multi-Frequency Low-Cost GNSS Receivers for Geodetic Monitoring Purposes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

2-D Joint Sparse Reconstruction and Micro-Motion Parameter Estimation for Ballistic Target Based on Compressive Sensing

1
National Lab of Radar Signal Processing, Xidian University, Xi’an 710071, China
2
School of Electronics and Communication Engineering, Sun Yat-Sen University, Guangzhou 510275, China
*
Author to whom correspondence should be addressed.
Submission received: 29 June 2020 / Revised: 4 August 2020 / Accepted: 4 August 2020 / Published: 5 August 2020
(This article belongs to the Section Physical Sensors)

Abstract

:
The sparse frequency band (SFB) signal presents a serious challenge to traditional wideband micro-motion curve extraction algorithms. This paper proposes a novel two-dimension (2-D) joint sparse reconstruction and micro-motion parameter estimation (2D-JSR-MPE) algorithm based on compressive sensing (CS) theory. In this technique, the 2D-JSR signal model and the micro-motion parameter dictionary are established based on the segmented SFB echo signal, in which the idea of piecewise effectively reduces the model complexity of ballistic target. With the accommodation of the CS theory, the 2D-JSR-MPE of the echo signal is transformed into solving a sparsity-driven optimization problem. Via an improved orthogonal matching pursuit (OMP) algorithm, the high-resolution range profiles (HRRP) can be reconstructed accurately, and the precise micro-motion curves can be simultaneously extracted on phase accuracy. The employment of 2-D joint processing can effectively avoid the interference of the sparse reconstruction error caused by cascaded operation in the subsequent micro-motion parameter estimation. The proposed algorithm benefits from the anti-jamming characteristic of the SFB signal and 2-D joint processing, thus remarkably enhancing its accuracy, robustness and practicality. Extensive experimental results are provided to verify the effectiveness and robustness of the proposed algorithm.

1. Introduction

The ballistic missile has become the leading weapon in modern wars. Missiles are often accompanied by decoys in the middle of the flight trajectory, and the micro-motion characteristics of warheads and decoys are obviously variant depending on their mass distribution. This makes estimating the target’s micro-motion parameters an important technique for target recognition in the missile defense system [1,2,3,4]. At present, methods such as phase derived ranging (PDR) technique are commonly used to extract micro-motion features of the ballistic target by processing wideband echo signals. By recourse to the phase information of range profile of the wideband echo signal, this type of approach can accurately depict the motion state of the scattering center and extract the micro-motion feature, thereby making the ranging accuracy reach the level of half a wavelength [5,6]. However, the traditional PDR algorithm above is based on a high-resolution range profile (HRRP) and does not use the 2-D coherent accumulation gain of echo signal, so it is susceptible to noise. In addition, the wideband radar is prone to environmental interference, which undermines its stability. The wideband signal requires the receiver to have the corresponding instantaneous bandwidth, posing a severe challenge to the design of the radar system and limiting its practicality in reality. In order to solve these problems, a sparse frequency band (SFB) signal has been developed, including a sparse orthogonal frequency division multiplexing (OFDM) signal and a sparse stepped-frequency signal [7,8]. Only transmitting part of the full frequency band (FFB), the SFB signal can effectively avoid certain interference frequency bands and reduce the data sampling amount, thereby improving the anti-jamming ability of the radar system and decreasing the operation burden of the radar hardware. However, the SFB signal severely challenges the traditional micro-motion curve extraction algorithms. Hence, it is of great significance to study the micro-motion extraction in SFB.
For the SFB signal, the reconstruction of HRRP is essential. The existing sparse reconstruction methods can be classified into four main categories. The first is nonlinear filtering represented by sequence CLEAN (S-CLEAN) [9,10]. It suppresses the high grating and side lobes of the pulse response function by remaining the peak value of the main lobe. The second is a bandwidth extrapolation (BWE) algorithm, such as the Burg algorithm [11,12]. This kind of method usually uses the linear prediction model to fit observation data, then adopts the spectrum estimation method to estimate the model parameters, and improves resolution through aperture extrapolation. The third is the parametric spectrum estimation technique, including the RELAX algorithm [13,14]. This type of algorithm establishes a parametric model for the echo signal, and precisely estimates the positions and amplitudes of the scattering centers through spectrum estimation so as to achieve high-resolution imaging. However, the above algorithms are all easily affected by noise and model errors, so the sparse reconstruction algorithm based on compressive sensing (CS) theory has been proposed in recent years. This kind of algorithm fulfills sparse data reconstruction by solving an l 1 -norm optimization problem, and thus is highly robust against noise [15,16]. When it comes to the micro-motion target, the 2-D coupling problem exists in the echo signal duo to micro-motion. Thus, if the 2-D joint micro-motion parameter estimation can be achieved simultaneously with sparse reconstruction, higher accuracy in estimation and stronger robustness can be guaranteed. However, by implementing sparse reconstruction in a single dimension, the above four types of methods fail to make the most of the 2-D coupling information and the 2-D coherent accumulation gain of the echo signal, which adversely affects the accuracy of signal reconstruction and parameter estimation as well as the robustness against noise.
Inspired by the above problems, this paper presents a 2-D joint sparse reconstruction and micro-motion parameter estimation (2D-JSR-MPE) algorithm aimed at the SFB signal. In this technique, a novel 2-D joint sparse reconstruction signal model combined with a segmented micro-motion characteristic parameter dictionary is established, in which the piecewise processing effectively reduces the model complexity of ballistic target. Based on this, a 2D-JSR-MPE optimization function is developed by means of the CS sparse representation, transforming the sparse reconstruction and parameter estimation into solving a sparsity-driven optimization problem with l 1 -norm. With the accommodation of the improved orthogonal matching pursuit (OMP) algorithm [17], micro-motion parameter estimation and sparse reconstruction can be simultaneously realized, thereby the high-precision HRRP and micro-motion curve can be obtained. Two-dimensional joint processing makes it possible to prevent the interference of the sparse reconstruction error in micro-motion parameter estimation, thus effectively avoiding the adverse effect of error transmission. The proposed algorithm estimates the micro-motion parameters by processing the phase information, so the estimation accuracy can reach the level of half a wavelength. Moreover, since the 2-D joint processing makes the most of the 2-D coupling information and 2-D coherent accumulation gain of the echo signal, and the SFB signal performs well in anti-jamming with low computational complexity, the proposed algorithm has the advantages of high accuracy, strong robustness, good practicality and low operational burden. Abundant experimental results corroborate the superiority of the proposed algorithm.
This paper is organized as follows. Section 2 introduces the 2D-JSR-MPE signal model of the ballistic target. Section 3 explicates the 2D-JSR-MPE algorithm. Section 4 provides the experimental results to illustrate the effectiveness of the proposed algorithm. The paper ends with a brief conclusion in Section 5.

2. Signal Model

2.1. The Micro-Motion Model of Ballistic Target

The micro-motion model of ballistic target is shown in Figure 1a. ( x , y , z ) is the radar coordinate system and ( X , Y , Z ) are the target body coordinate system. RLOS is the radar line of sight, and the elevation angle and the azimuth angle of RLOS are α and β , respectively. The geometry model of ballistic target is shown in Figure 1b. H denotes the height of the cone, r refers to the radius of the bottom, and d represents the distance between the mass center and the top of the cone. For the ballistic cone target model established in this paper, there are three equivalent scattering centers on the target [2], namely p 1 , p 2 and p 3 in Figure 1b. When moving in precession, the target makes a spin motion on its axis of symmetry at an angular velocity of ω s , and makes a conic motion on the conic axis OM at an angular velocity of ω c . The precession angle is θ . When the target moves in nutation, on the basis of precession, the cone axis of the target oscillates in accordance with the sine wave within a certain range. The oscillation angle is κ ( t m ) = δ sin ( ω v t m ) ; δ and ω v stand for the maximum amplitude and the angular velocity of oscillation, respectively. Given that the ballistic cone target studied in this paper is axisymmetric, its spin motion does not affect the radar echo of the target and is therefore not taken into consideration.
In this paper, it is assumed that the translational motion of the target has been compensated, then the instantaneous radial distance between the k th scattering center and the radar can be calculated by the projection from the k th scattering center to RLOS
r k ( t m ) = Y k sin γ ( t m ) + Z k cos γ ( t m ) + R 0
where R 0 represents the distance from the mass center of target to the radar; γ ( t m ) denotes the angle between the symmetry axis of the target and RLOS at t m . When the target moves in precession, the cosine of the angle between the symmetry axis of the target and RLOS can be expressed as (2)
cos γ ( t m ) = cos θ cos ( θ + α ) + sin θ sin ( θ + α ) cos ( ω c t m + φ )
where φ denotes initial phase. When the target moves in nutation, the cosine of the angle between the symmetry axis of the target and RLOS can be expressed as (3)
cos γ ( t m ) = cos θ [ cos ( θ + α ) cos [ κ ( t m ) ] sin ( α + 3 θ ) sin [ κ ( t m ) ] ] + sin θ cos ( ω c t m ) [ cos ( α + 3 θ ) sin [ κ ( t m ) ] sin ( θ + α ) cos [ κ ( t m ) ] ]
Assigning (2) or (3) to (1), the theoretical expression of the instantaneous distance between each scattering center and the radar can be given by (4)
{ r k ( t m ) = R 0 + R k ( t m ) , k = 1 , 2 , 3 R 1 ( t m ) = d cos γ ( t m ) R 2 ( t m ) = ( H d ) cos γ ( t m ) r sin γ ( t m ) R 3 ( t m ) = ( H d ) cos γ ( t m ) + r sin γ ( t m )
where R k ( t m ) represents the instantaneous micro distance of each scattering center. It can be seen that the instantaneous micro distance of each scattering center changes with time in a form approximating to the sine wave and includes the target’s motion and geometry parameters.

2.2. D-JSR-MPE Signal Model

It is assumed that the radar transmits a linear frequency modulated (LFM) signal. After preprocessing such as range compression and de-modulation to the baseband, the received signal of the k th scattering center can be given by (5) [18]
S k ( f r , t m ) = σ k · exp [ j 4 π ( f c + f r ) c · R k ( t m ) ]
where σ k represents the scattering coefficient of the k th scattering center; c denotes the light velocity. From (4), it can be seen that the instantaneous micro distance of each scattering center can be expressed as a sum of the multi-order Sine functions, but it is complex to cope with such forms in practice. Thus, underpinned by the idea of piecewise [19], in a short period of time, using a second-order polynomial model to approximate the instantaneous micro distance of scattering points can reduce the complexity of the model effectively. Moreover, because of the short dwell time, the second-order polynomial is accurate enough to fit the instantaneous micro-motion trajectory of the ballistic target, which is conducive to the subsequent reconstruction of the micro-motion curve. In this paper, the echo signal is segmented evenly, and the approximate expression of the instantaneous micro distance of the j th segment of the k th scattering point can be written as (6)
R j k ( t m ) = a j k + b j k t m + c j k t m 2
where a j k , b j k and c j k denote micro-motion parameters. By estimating them accurately, the micro-motion curve can be reconstructed. Assigning (6) to (5), the j th segmented echo signal can be given by (7)
S j k ( f r , t m ) = σ k · exp [ j 4 π ( f c + f r ) c · ( a j k + b j k t m + c j k t m 2 ) ]
After sampling processing, the discrete form of (7) can be expressed as (8)
S j k ( n , m ) = σ k · exp [ j 4 π ( f c + n Δ f ) c · ( a j k + b j k ( m Δ t ) + c j k ( m Δ t ) 2 ) ]
where n refers to the range frequency bin index, and m the azimuth slow time index. The ranges of n and m are as follows: n [ N / 2   + 1 : N / 2   ] and m [ M / 2   + 1 : M / 2   ] , where N refers to the total number of the range frequency bin of FFB signal, and M stands for the number of pulses contained in the j th segmented echo signal. Δ f corresponds to the range frequency sampling interval, and Δ t the slow time sampling interval.
Equation (8) is the expression corresponding to the FFB signal. When the range frequency band is sparse, the geometry model of the SFB signal is shown in Figure 2. It is assumed that the SFB signal consists of Q sparse subbands, in which the q th subband is L q in length, running from f q to f q + L q 1 , and f q denotes the initial range frequency bin index of the q th subband. It is assumed that the q th subband signal can be represented by the matrix S j , q k which is L q × M in size. Then, (8) can be expressed in the form of matrix, as shown in (9)
S j k = [ S j , 1 k S j , q k S j , Q k ] L ¯ × M
where L ¯ = q = 1 Q L q . Assuming that there are a total of K scattering centers in the target, the SBF echo signal of the target can be expressed as (10)
S j = [ S j , 1 S j , q S j , Q ] L ¯ × M
where S j , q = k = 1 K S j , q k . Considering the inevitable noise and combined with (8) and (10), S j can be re-expressed as the following matrix form
S j = D j F r · Y j T · F a T + N j
where represents the Hadamard product and D j L ¯ × M denotes the micro-motion parameter matrix as (12)
D j = [ exp ( j 4 π ( Δ f + f c ) c P j ( Δ t ) ) exp ( j 4 π ( Δ f + f c ) c P j ( 2 Δ t ) ) exp ( j 4 π ( Δ f + f c ) c P j ( M Δ t ) ) exp ( j 4 π ( L 1 Δ f + f c ) c P j ( Δ t ) ) exp ( j 4 π ( L 1 Δ f + f c ) c P j ( 2 Δ t ) ) exp ( j 4 π ( L 1 Δ f + f c ) c P j ( M Δ t ) ) exp ( j 4 π ( L ¯ Δ f + f c ) c P j ( Δ t ) ) exp ( j 4 π ( L ¯ Δ f + f c ) c P j ( 2 Δ t ) ) exp ( j 4 π ( L ¯ Δ f + f c ) c P j ( M Δ t ) ) ] L ¯ × M
where P j ( m Δ t ) = b j ( m Δ t ) + c j ( m Δ t ) 2 , m = 1 , 2 , M , the ranges of b j and c j are as follows: b j [ η / 2   : η / 2   ] and c j [ μ / 2   : μ / 2   ] , in which η and μ should be large enough to ensure that b j and c j can include the center frequencies and chirp rates of all sub-segments echoes. F r = [ F r 0 ; ; F r q ; ; F r ( Q 1 ) ] L ¯ × N represents the partial Fourier dictionary of range dimension, F r q = [ f r q f q ; ; f r q f q + L q 1 ] L q × N , where f r q x = [ 1 , ω r x , , ω r ( N 1 ) x ] , x = f q , f q + 1 , , f q + L q 1 , ω r = exp [ j 2 π / N   ] . F a M × M denotes the standard Fourier dictionary of azimuth dimension, whose detailed form is not given here. Y j M × N corresponds to the 2-D full-resolution ISAR image and N j L ¯ × M the additive complex Gaussian white noise matrix.
Equation (11) is the 2D-JSR-MPE signal model proposed in this paper. In the next section, we will explicate the HRRP reconstruction and micro-motion curve estimation of scattering centers based on the CS theory.

3. D-JSR-MPE Algorithm

3.1. The Principle of 2D-JSR-MPE Based on CS

The CS theory points out that by analyzing the sparsity of the echo signal, constructing the measurement matrix, and using the existing optimization algorithms, the sparse signal can be reconstructed with a little data sampling amount. Combined with the derivation in the previous section, in order to facilitate the construction and solution of the optimization function, the vectorization operation is carried out for (11)
v e c [ S j ] = v e c [ D j F r · Y j T · F a T ] + v e c [ N j ]
where v e c [ · ] refers to the vectorization operation. For clarity, we mark s j = v e c [ S j ] ( L ¯ · M ) × 1 , y j = v e c [ Y j T ] ( N · M ) × 1 , and n j = v e c [ N j ] ( L ¯ · M ) × 1 . Based on this, (13) can be expressed as (14) according to the relationship between vectorization operation and Kronecker product [20,21,22]
s j = G j · ( F a F r ) · y j + n j = G j · W · y j + n j
where represents Kronecker product; W ( L ¯ · M ) × ( N · M ) denotes the 2-D joint sparse reconstruction dictionary, and W = F a F r . In particular, G j ( L ¯ · M ) × ( L ¯ · M ) is a diagonal matrix composed of all elements in D j as the diagonal elements. Based on the CS theory and (14), we can establish a 2D-JSR-MPE optimization function as (15)
min y j 1 , s t . s j G j · W · y j 2 < ε j
where ε j represents the noise threshold; · i denotes the i -norm of vector. According to the CS theory, 2D-JSR-MPE can be realized by solving the l 1 -norm optimization problem shown in (15). Therefore, in the next subsection, we will elucidate an algorithm which can quickly and accurately address the optimization problem shown in (15).

3.2. 2D-JSR-MPE Based on the Improved OMP Algorithm

Due to the large size of G j · W , it is difficult for common algorithms to solve (15). In this subsection, via fast Fourier transform (FFT), the improved OMP algorithm is employed to jointly realize HRRP reconstruction and micro-motion curve estimation. The improved OMP algorithm mainly includes the following four steps:
(1)
Micro-motion curve parameters estimation. Firstly, the micro-motion curve parameters of the first scattering center in s j are estimated by (16)
a j 1 , b j 1 , c j 1 = max | ( H j · W ) H s j |
where H j = A j G j , and A j ( L ¯ · M ) × ( L ¯ · M ) is a diagonal matrix composed of all elements in a j as diagonal elements as (17)
a j = [ exp ( j 4 π ( Δ f + f c ) c a j ) exp ( j 4 π ( Δ f + f c ) c a j ) exp ( j 4 π ( Δ f + f c ) c a j ) exp ( j 4 π ( L 1 Δ f + f c ) c a j ) exp ( j 4 π ( L 1 Δ f + f c ) c a j ) exp ( j 4 π ( L 1 Δ f + f c ) c a j ) exp ( j 4 π ( L ¯ Δ f + f c ) c a j ) exp ( j 4 π ( L ¯ Δ f + f c ) c a j ) exp ( j 4 π ( L ¯ Δ f + f c ) c a j ) ] L ¯ × M
For solving (16), the traditional method needs to assign all values of a j 1 , b j 1 , c j 1 to H j , respectively, the computational complexity of each iteration is O ( N M 2 L ¯ ) , which greatly increases the computational amount. In this paper, FFT is introduced into the process of solving a j 1 , b j 1 , c j 1
U = F F T 2 { S j D j }
where F F T 2 { · } represents the 2-D fast Fourier transform (2D-FFT). Based on (18), a ^ j 1 , b ^ j 1 , c ^ j 1 , the optimal estimation of a j 1 , b j 1 , c j 1 , can be determined by searching the maximum value of U . Then, the basis matrix corresponding to SFB Φ j = [ H j , a ^ j 1 , b ^ j 1 , c ^ j 1 ] can be constructed by a ^ j 1 , b ^ j 1 , c ^ j 1 , and the homologous basis matrix with respect to FFB established as Φ j F F B .
(2)
Scattering center estimation. After obtaining Φ j , the weighted least square estimation (WLSE) method is adopted to estimate the principal component of the scattering centers in y j , and the corresponding scattering coefficient is Ψ ^ j = [ Ψ ^ j 1 ] = ( Φ j H Φ j ) 1 Φ j H S j . On this basis, we can reconstruct the corresponding SFB observation data S ^ j = Φ j Ψ ^ j .
(3)
Residual signal estimation. First, the residual SFB signal S ¯ j is calculated by subtracting S ^ j from S j , i.e., S ¯ j = S j S ^ j . Then, according to step 1 and step 2, the parameter estimation and scattering center extraction of the principle component signal in S ¯ j are carried out. The estimated micro-motion curve parameters are a ^ j 2 , b ^ j 2 , c ^ j 2 . In this case, a ^ j 2 , b ^ j 2 , c ^ j 2 are used to update the basis matrix to Φ j = [ H j , a ^ j 1 , b ^ j 1 , c ^ j 1 , H j , a ^ j 2 , b ^ j 2 , c ^ j 2 ] , and the scattering coefficients are estimated as Ψ ^ j = [ Ψ ^ j 1 ; Ψ ^ j 2 ] = ( Φ j H Φ j ) 1 Φ j H S j . It should be noted that the scattering coefficient of the principal component signal estimated in the first time and Φ j F F B are also updated here.
(4)
Iterative operation. Repeat Step 2 and Step 3 until S ¯ j 2 < ε j is met.
Through the above processing, the micro-motion curve of the k th scattering center of the j th segment can be reconstructed by (19)
R ^ j k ( t m ) = a ^ j k + b ^ j k t m + c ^ j k t m 2
Moreover, the reconstructed HRRP of the j th segment can be obtained using range inverse FFT (IFFT) for the reconstructed FFB signal S ^ j F F B = Φ j F F B Ψ ^ j . The final reconstructed HRRPs and micro-motion curves can be obtained by splicing and integrating the processing results of each subsegment. Through 2-D joint processing, sparse reconstruction and micro-motion parameter estimation are carried out simultaneously, which avoids the cascaded error transfer resulting from “reconstruction first and then estimation”, thus improving the accuracy and robustness of parameter estimation. Since the proposed algorithm taps the phase information to estimate the micro-motion parameters, the estimation accuracy can reach the level of half a wavelength, rendering the proposed algorithm highly accurate in terms of parameter estimation.
The computational complexity of the proposed algorithm mainly arises from the FFT operation in Step 1 and the matrix inversion operation in Step 2. The computational complexity of the matrix inversion operation in Step 2 is O ( k = 1 K k 3 ) , which can be neglected since there are only a few equivalent scattering centers for the ballistic target studied in this paper. Therefore, it is Step 1 that mainly determines the complexity of the proposed algorithm, which can be expressed by the number of complex multiplications [23]. It can be seen from (18) that 2D-FFT is employed to estimate the micro-motion parameters of each subsegment. The computational complexity of 2D-FFT is O ( M · N log 2 ( M · N ) ) , and the search ranges of b j and c j are η and μ , respectively, so the computational complexity corresponding to each subsegment is O ( η · μ · M · N log 2 ( M · N ) ) . Given that the echo signal is divided into J segments, the total computational complexity of the proposed algorithm is O ( J · η · μ · M · N log 2 ( M · N ) ) . Through FFT and Hadamard product, the computational complexity of each iteration of the proposed algorithm O ( M · N log 2 ( M · N ) ) is much smaller than that of direct matrix operation O ( N M 2 L ¯ ) , so the efficiency of the 2-D joint processing algorithm can be effectively improved by fast operation.
In order to describe the proposed algorithm clearly, a flowchart is given in Figure 3.

4. Experiments and Analyses

The geometry model of the ballistic target shown in Figure 1b is adopted in this simulation experiment, and the main simulated target and radar system parameters are shown in Table 1 and Table 2, respectively. In this experiment, three comparison algorithms are employed to illustrate the advantages of the proposed algorithm. The ‘zero-padding FFT (ZP-FFT)’ method [8] carries out a zero padding operation on the part of the vacant frequency band, and obtains the range profiles of the scattering centers through range IFFT. For the ‘S-CLEAN’ method, the HRRPs of the scattering centers are reconstructed by the S-CLEAN algorithm [10], and then the micro-motion curve is extracted by the traditional PDR algorithm. The ‘2-D cascaded processing algorithm (2D-CPA)’ uses the OMP algorithm [17] to reconstruct HRRP in the range dimension first, and then superimposes the target support area along the range bin. For the superimposed one-dimension (1-D) vector, the modified discrete chirp Fourier transform (MDCFT) [24] is utilized to estimate the parameters and reconstruct the micro-motion curve in the azimuth dimension. With respect to the proposed algorithm, for clarity, we name it as ‘2D-JSR-MPE’. It is worth noting that the HRRPs obtained by ZP-FFT have high gating and side lobes, so ZP-FFT is not applied to the subsequent estimation of micro-motion parameters in this experiment. In what follows, the effectiveness and superiority of the proposed algorithm are demonstrated by the experimental results obtained by the three comparison algorithms and the proposed algorithm at different signal to noise ratios (SNRs) and sparse rates (SRs).
First, the precession motion of the target is studied, and the experiment is carried out at SNR = 10 dB. Figure 4a shows the FFB echo waveform in range-frequency and azimuth-time domains at SNR = 10 dB with its corresponding HRRPs in Figure 4b. It can be seen that due to shielding effect, only one micro-motion curve corresponding to scattering centers p 2 or p 3 can be observed [25]. It is widely acknowledged that the sparse sampling pattern can be classified into two categories: gap missing sampling (GMS) and random missing sampling (RMS). They are only different in the sampling method, while identical in the subsequent experimental processing. The SFB signal (such as the sparse OFDM signal) is usually transmitted in the form of GMS. Therefore, in this experiment, we use GMS to illustrate the performance of the algorithms. In addition, two kinds of SRs are set, SR1 = 0.75 and SR2 = 0.5. The range dimension sampling points corresponding to SR1 and SR2 are three fourths and half as many as the FFB sampling points in Figure 4a, respectively. Figure 5 shows the SFB waveforms at SR1 and SR2 as well as the HRRPs obtained by ZP-FFT. As shown in Figure 5b,d, the HRRPs obtained by ZP-FFT have high gating and side lobes with low resolution, which seriously hinders the follow-up work.
In the following experiment, we study the influence of the SR on the algorithms. When the target is in precession, the reconstructed results of HRRP obtained by 2D-JSR-MPE, 2D-CPA and S-CLEAN in the case of SR1 and SR2 are shown in Figure 6. It can be seen that in the case of SR1, all three algorithms can reconstruct good HRRPs. However, at SR2, due to the influence of high grating and side lobes, severe errors occur in the HRRPs reconstructed by S-CLEAN. Figure 7 shows the micro-motion curves reconstructed by the three algorithms corresponding to precession at the two SRs. Since the S-CLEAN method is used to extract the maximum peak point to achieve the suppression of high grating and side lobes, with the decrease in SR, the error arising from extracting peak point leads to the micro-motion curve reconstruction failure. 2D-CPA fails to make full use of 2-D coupling information and 2-D coherent accumulation gain of echo signal due to cascaded processing, and is prone to the estimation error transmission, so its reconstruction accuracy is insufficient. In contrast, 2D-JSR-MPE can accurately reconstruct the micro-motion curve in the same condition.
In order to study the influence of SNR on the algorithms, HRRPs are reconstructed by 2D-JSR-MPE, 2D-CPA and S-CLEAN at SNR = −10 dB, 0 dB and 10 dB. The homologous results are shown in Figure 8. It can be seen that there are a lot of false points in the HRRPs reconstructed by 2D-CPA and S-CLEAN at low SNR, but 2D-JSR-MPE succeed in reconstructing HRRPs accurately. Figure 9 reveals the micro-motion curves reconstructed by the three algorithms under the three SNR conditions. As SNR decreases, 2D-CPA and S-CLEAN are seriously disturbed by noise, and the error of the reconstructed micro-motion curves increases gradually. However, compared with the above two comparison algorithms, 2D-JSR-MPE can still reconstruct the accurate micro-motion curves at different SNRs. The experiments above corroborate the good robustness of the proposed algorithm against SR and SNR.
In addition, based on the precession motion of the target, the oscillating motion of the symmetry axis of the target is introduced into the experiment to estimate the micro-motion curves of the target with nutation motion. Figure 10 and Figure 11 are the corresponding micro-motion curves reconstructed by the three algorithms at the two kinds of SRs and the three types of SNRs when the target is in nutation. They suggest that the experimental conclusions drawn in the complex nutation are the same as those acquired in the precession, which further verify the effectiveness of the proposed algorithm.
For the purpose of analyzing the impact of SR and SNR on the performances of the proposed algorithm and the two comparison algorithms quantitatively, the estimated correct rate (ECR) of the reconstructed micro-motion curve is defined by (20)
E C R = ( 1 t m | F r ( t m ) F e ( t m ) | t m | F r ( t m ) | ) 100 %
where F e ( t m ) represents the estimated value of the micro-motion curve and F r ( t m ) the true value of the micro-motion curve. Firstly, the influence of SR on the performance of the algorithms is analyzed. In the case of precession motion, with other conditions unchanged, SR ranges from 0.25 to 1 in step of 0.125 at SNR = 10 dB. Without loss of generality, 40 groups of echo data are collected at each SR, followed by the reconstruction of their micro-motion curves. After calculating the ECR of each curve, the average value of 40 groups of ECRs can be obtained, which is shown in Figure 12a. Compared with the other two algorithms, the proposed one has higher ECRs at all SRs.
Finally comes the analysis of the influence of SNR on the algorithms’ performance. Assuming other conditions remain unchanged, the robustness of the algorithms against noise is analyzed with SR = 0.75 and SNR ranging from −15 to 20 dB in step of 5 dB. The average ECR is calculated by the same method, and the result is shown in Figure 12b. It can be seen that the ECRs of the proposed algorithm are above 90% at different SNRs, higher than those of the other two algorithms. This serves as convincing proof of the good robustness of the proposed algorithm against noise.

5. Conclusions

This paper proposes a novel 2D-JSR-MPE algorithm. Through piecewise processing, we establish the 2-D joint sparse reconstruction signal model and the target’s micro-motion characteristic parameter dictionary, in which the idea of piecewise effectively reduces the model complexity of ballistic target. Based on the CS theory, an improved OMP algorithm is employed to solve a sparse optimization problem with l 1 -norm. With the help of 2-D joint processing, the precise reconstruction of HRRP and micro-motion curve can be realized simultaneously, thereby avoiding the error transmission between HRRP reconstruction and micro-motion parameters estimation. Making full use of the 2-D coupling information and 2-D coherent accumulation gain of echo signal, the proposed algorithm can accurately reconstruct the micro-motion curve of the ballistic target on phase accuracy. Extensive experimental results demonstrate that the proposed algorithm can be implemented with high accuracy and strong robustness at different SNRs and SRs. Further efforts to improve the performance of our work in this paper in the case of sparse frequency band and sparse aperture (SFB-SA) echo signal of ballistic target are underway.

Author Contributions

J.W. is responsible for all the theoretical work, the implementation of the experiments and the writing of the manuscript. L.Z., S.S. and H.L. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China with grant numbers 61771372 and 61771367, the National Science Foundation for Distinguished Young Scholars with grant number 61525105, and the Fund for Foreign Scholars in University Research and Teaching Programs (the 111 Project) with grant number B18039.

Acknowledgments

The authors thank the anonymous reviewers for their valuable comments to improve the paper quality.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, V.C.; Li, F.; Ho, S.S.; Wechsler, H. Micro-Doppler effect in radar: Phenomenon, model, and simulation study. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 2–21. [Google Scholar] [CrossRef]
  2. Gao, H.; Xie, L.; Wen, S.; Kuang, Y. Micro-Doppler Signature Extraction from Ballistic Target with Micro-Motions. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1969–1982. [Google Scholar] [CrossRef]
  3. Xing, Y.; You, P.; Yong, S. Parameter Estimation of Micro-Motion Targets for High-Resolution-Range Radar Using Online Measured Reference. Sensors 2018, 18, 2773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wu, Y.; Lu, H.; Zhao, F.; Zhang, Z. Estimating Shape and Micro-Motion Parameter of Rotationally Symmetric Space Objects from the Infrared Signature. Sensors 2016, 16, 1722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Camp, W.W.; Mayhan, J.T.; O’Donnell, R.M. Wideband radar for ballistic missile defense and range-doppler imaging for satellites. Linc. Lab. J. 2000, 12, 267–280. [Google Scholar]
  6. Liu, Y.; Zhu, D.; Li, X.; Zhuang, Z. Micromotion characteristic acquisition based on wideband radar phase. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3650–3657. [Google Scholar] [CrossRef]
  7. Su, X.; Liu, Z.; Chen, X.; Wei, X. Closed-Form Algorithm for 3-D Near-Field OFDM Signal Localization under Uniform Circular Array. Sensors 2018, 18, 226. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, L.; Qiao, Z.; Xing, M.; Li, Y.; Bao, Z. High-resolution ISAR imaging with sparse stepped-frequency waveforms. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4630–4651. [Google Scholar] [CrossRef]
  9. Bose, R.; Freeman, A.; Steinberg, B.D. Sequence CLEAN: A modified deconvolution technique for microwave images of contiguous targets. IEEE Trans. Aerosp. Electron. Syst. 2002, 38, 89–97. [Google Scholar] [CrossRef]
  10. Martorella, M.; Acito, N.; Berizzi, F. Statistical CLEAN Technique for ISAR Imaging. IEEE Trans. Geosci. Remote Sens. 2007, 45, 3552–3560. [Google Scholar] [CrossRef]
  11. Li, H.J.; Farhat, N.H.; Shen, Y. A new iterative algorithm for extrapolation of data available in multiple restricted regions with applications to radar imaging. IEEE Trans. Antennas. Propag. 1987, 35, 581–588. [Google Scholar]
  12. Gupta, I.; Beals, M.; Moghaddar, A. Data extrapolation for high resolution radar imaging. IEEE Trans. Antennas. Propag. 1994, 42, 1540–1545. [Google Scholar] [CrossRef]
  13. Huan, S.; Zhang, M.; Dai, G.; Gan, H. Low Elevation Angle Estimation with Range Super-Resolution in Wideband Radar. Sensors 2020, 20, 3104. [Google Scholar] [CrossRef] [PubMed]
  14. Lane, R.O.; Copsey, K.D.; Webb, A.R. A Bayesian approach to simultaneous autofocus and super-resolution. Proc. Spie Int. Soc. Opt. Eng. 2004, 5427, 133–142. [Google Scholar]
  15. Rauhut, H.; Schnass, K.; Vandergheynst, P. Compressed sensing and redundant dictionaries. IEEE Trans. Inf. Theory 2008, 54, 2210–2219. [Google Scholar] [CrossRef] [Green Version]
  16. Chung, H.; Park, Y.M.; Kim, S. Wideband DOA Estimation on Co-prime Array via Atomic Norm Minimization. Sensors 2020, 13, 3235. [Google Scholar] [CrossRef]
  17. Zhang, L.; Duan, J.; Qiao, Z.J.; Xing, M.; Bao, Z. Phase adjustment and isar imaging of maneuvering targets with sparse apertures. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 1955–1973. [Google Scholar] [CrossRef]
  18. Shao, S.; Zhang, L.; Liu, H.; Zhou, Y. Accelerated translational motion compensation with contrast maximisation optimisation algorithm for inverse synthetic aperture radar imaging. IET Radar Sonar Navig. 2019, 13, 316–325. [Google Scholar] [CrossRef]
  19. Voros, J. Modeling and parameter identification of systems with multisegment piecewise-linear characteristics. IEEE Trans. Autom. Control 2002, 47, 184–188. [Google Scholar] [CrossRef]
  20. Jahromi, M.J.; Kahaei, M.H. Two-dimensional iterative adaptive approach for sparse matrix solution. Electron. Lett. 2014, 50, 45–47. [Google Scholar] [CrossRef]
  21. Duarte, M.F.; Baraniuk, R.G. Kronecker Compressive Sensing. IEEE Trans. Image Process. 2012, 21, 494–504. [Google Scholar] [CrossRef] [PubMed]
  22. Shao, S.; Zhang, L.; Liu, H. High-resolution ISAR imaging and motion compensation with 2-D joint sparse reconstruction. IEEE Trans. Geosci. Remote Sens. 2020, 1–21. [Google Scholar] [CrossRef]
  23. Shao, S.; Zhang, L.; Liu, H.; Zhou, Y. Spatial-variant contrast maximization auto focus algorithm for ISAR imaging of maneuvering targets. Sci. China Inf. Sci. 2019, 62, 37–39. [Google Scholar] [CrossRef] [Green Version]
  24. Xu, G.; Xing, M.; Zhang, L.; Duan, J.; Chen, Q.; Bao, Z. Sparse Apertures ISAR Imaging and Scaling for Maneuvering Targets. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 2942–2956. [Google Scholar] [CrossRef]
  25. Niu, J.; Zhu, Y.; Liu, Y.; Li, X.; Kuang, G. Radar micro-motion feature extraction of cone target based on Hilbert-Huang transform. In Proceedings of the IEEE International Conference on Information Science and Technology, Hubei, China, 23–25 March 2012; pp. 401–404. [Google Scholar]
Figure 1. Ballistic target model: (a) Micro-motion model; (b) Geometry model.
Figure 1. Ballistic target model: (a) Micro-motion model; (b) Geometry model.
Sensors 20 04382 g001
Figure 2. The geometry model of SFB signal.
Figure 2. The geometry model of SFB signal.
Sensors 20 04382 g002
Figure 3. Flowchart of 2D-JSR-MPE.
Figure 3. Flowchart of 2D-JSR-MPE.
Sensors 20 04382 g003
Figure 4. FFB situation: (a) Echo waveform; (b) HRRPs.
Figure 4. FFB situation: (a) Echo waveform; (b) HRRPs.
Sensors 20 04382 g004
Figure 5. SFB situation. SR1: (a) Echo waveform; (b) HRRPs obtained by ZP-FFT. SR2: (c) Echo waveform; (d) HRRPs obtained by ZP-FFT.
Figure 5. SFB situation. SR1: (a) Echo waveform; (b) HRRPs obtained by ZP-FFT. SR2: (c) Echo waveform; (d) HRRPs obtained by ZP-FFT.
Sensors 20 04382 g005
Figure 6. The reconstructed results of HRRP of the target with precession motion obtained by different algorithms with different SRs at SNR = 10 dB: (a) FFB situation; (b) 2D-JSR-MPE; (c) 2D-CPA; (d) S-CLEAN.
Figure 6. The reconstructed results of HRRP of the target with precession motion obtained by different algorithms with different SRs at SNR = 10 dB: (a) FFB situation; (b) 2D-JSR-MPE; (c) 2D-CPA; (d) S-CLEAN.
Sensors 20 04382 g006
Figure 7. The reconstructed results of micro-motion curves of the target with precession motion obtained by different algorithms with different SRs at SNR = 10 dB: (a) 2D-JSR-MPE; (b) 2D-CPA; (c) S-CLEAN.
Figure 7. The reconstructed results of micro-motion curves of the target with precession motion obtained by different algorithms with different SRs at SNR = 10 dB: (a) 2D-JSR-MPE; (b) 2D-CPA; (c) S-CLEAN.
Sensors 20 04382 g007
Figure 8. The reconstructed results of HRRP of the target with precession motion obtained by different algorithms with different SNRs at SR = 0.75: (a) FFB situation; (b) 2D-JSR-MPE; (c) 2D-CPA; (d) S-CLEAN.
Figure 8. The reconstructed results of HRRP of the target with precession motion obtained by different algorithms with different SNRs at SR = 0.75: (a) FFB situation; (b) 2D-JSR-MPE; (c) 2D-CPA; (d) S-CLEAN.
Sensors 20 04382 g008
Figure 9. The reconstructed results of micro-motion curves of the target with precession motion obtained by different algorithms with different SNRs at SR = 0.75: (a) 2D-JSR-MPE; (b) 2D-CPA; (c) S-CLEAN.
Figure 9. The reconstructed results of micro-motion curves of the target with precession motion obtained by different algorithms with different SNRs at SR = 0.75: (a) 2D-JSR-MPE; (b) 2D-CPA; (c) S-CLEAN.
Sensors 20 04382 g009
Figure 10. The reconstructed results of micro-motion curves of the target with nutation motion obtained by different algorithms with different SRs at SNR = 10 dB: (a) 2D-JSR-MPE; (b) 2D-CPA; (c) S-CLEAN.
Figure 10. The reconstructed results of micro-motion curves of the target with nutation motion obtained by different algorithms with different SRs at SNR = 10 dB: (a) 2D-JSR-MPE; (b) 2D-CPA; (c) S-CLEAN.
Sensors 20 04382 g010
Figure 11. The reconstructed results of micro-motion curves of the target with nutation motion obtained by different algorithms with different SNRs at SR = 0.75: (a) 2D-JSR-MPE; (b) 2D-CPA;(c) S-CLEAN.
Figure 11. The reconstructed results of micro-motion curves of the target with nutation motion obtained by different algorithms with different SNRs at SR = 0.75: (a) 2D-JSR-MPE; (b) 2D-CPA;(c) S-CLEAN.
Sensors 20 04382 g011
Figure 12. ECRs obtained by different algorithms: (a) Different algorithms with different SRs at SNR = 10 dB; (b) Different algorithms with different SNRs at SR = 0.75.
Figure 12. ECRs obtained by different algorithms: (a) Different algorithms with different SRs at SNR = 10 dB; (b) Different algorithms with different SNRs at SR = 0.75.
Sensors 20 04382 g012
Table 1. Main simulated parameters of the ballistic target.
Table 1. Main simulated parameters of the ballistic target.
ParametersValues
Height of the target4.0 m
Distance between the mass center and the top of the target2.7 m
Radius of the bottom0.3 m
Spinning frequency2 Hz
Coning frequency4 Hz
Oscillating frequency1 Hz
Precession angle10°
Table 2. Main simulated parameters of the radar system.
Table 2. Main simulated parameters of the radar system.
ParametersValues
Carrier frequency10 GHz
Bandwidth1 GHz
Pulse repetition frequency4 kHz
Coherent processing interval1 s

Share and Cite

MDPI and ACS Style

Wei, J.; Shao, S.; Zhang, L.; Liu, H. 2-D Joint Sparse Reconstruction and Micro-Motion Parameter Estimation for Ballistic Target Based on Compressive Sensing. Sensors 2020, 20, 4382. https://doi.org/10.3390/s20164382

AMA Style

Wei J, Shao S, Zhang L, Liu H. 2-D Joint Sparse Reconstruction and Micro-Motion Parameter Estimation for Ballistic Target Based on Compressive Sensing. Sensors. 2020; 20(16):4382. https://doi.org/10.3390/s20164382

Chicago/Turabian Style

Wei, Jiaqi, Shuai Shao, Lei Zhang, and Hongwei Liu. 2020. "2-D Joint Sparse Reconstruction and Micro-Motion Parameter Estimation for Ballistic Target Based on Compressive Sensing" Sensors 20, no. 16: 4382. https://doi.org/10.3390/s20164382

APA Style

Wei, J., Shao, S., Zhang, L., & Liu, H. (2020). 2-D Joint Sparse Reconstruction and Micro-Motion Parameter Estimation for Ballistic Target Based on Compressive Sensing. Sensors, 20(16), 4382. https://doi.org/10.3390/s20164382

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