Micromotion Feature Extraction of Space Target Based on Track-Before-Detect

1 Institute of Information and Navigation, Air Force Engineering University, Xi’an 710077, China 2Collaborative Innovation Center of Information Sensing and Understanding, Xi’an 710077, China 3Key Laboratory for Information Science of ElectromagneticWaves, Ministry of Education, FudanUniversity, Shanghai 200433, China 4Department of Electrical and Computer Engineering, National University of Singapore, Singapore 11576


Introduction
The monitoring and recognition of space targets have become more and more important since these space targets such as fragments of satellites, rocket bodies, and other space debris are hazardous to aerospace activities with the increasing outer space explorations [1].Generally, a space target has complex micromotion such as spinning, precession, and rolling, in addition to the body translation [2].These micromotions will induce the micro-Doppler (m-D) effect, which can be regarded as a unique signature, providing additional information for classification, recognition, and identification of the target [3].Therefore, since the concept of micro-Doppler was introduced to the radar signal processing field [4,5], the micromotion feature extraction of space target has drawn extensive attention of scholars [6][7][8][9][10][11][12][13][14][15][16].
Based on the fixed scatterer model [7], the micromotion feature of space target was usually analyzed with the timefrequency distribution technique, and the micromotion feature of target was extracted with parameter transformation method [7].However, the time-frequency methods have some limitations due to the fact that the change of micro-Doppler frequency with time is highly nonlinear [8].A feature extraction method based on empirical-mode decomposition (EMD) was proposed to estimate the oscillation frequency of the truck's surface in [9].To get the higher potential payoffs from the exploitation of m-D effect in wideband radar system, the micromotion feature extraction for wideband radar based on complex image orthogonal matching pursuit decomposition is proposed in [10].
However, many space targets are of conical and conecylinder shapes, these rotationally symmetric targets cannot be described exactly with the fixed scatterer model.Their scattering centers of the edge of the conic bottom slide with the radar line of sight (LOS) that means the scattering centers local at the two intersections of the conic bottom edge and the radar wave incident plane [11].Based on the slide scatterer model, the mathematical expression of micro-Doppler is derived, which is validated by the darkroom data in [12,13].However, these works did not provide effective micromotion feature extraction methods for the slide scatterer model.Unfortunately, most of the existing micromotion feature extraction methods for the fixed scatterer model will not function properly because the micro-Doppler curve of a slide scatterer is much different from that of a fixed scatterer.Although some methods based on the complex generalized Radon transform and extended Hough transform are proposed for micromotion feature extraction based on the slide scatterer model [14,15], the computation load of these methods is heavy, and the methods need high pulse repetition frequency (PRF) to avoid the frequency domain aliasing phenomenon.In this case, a lot of radar resources needs to be allocated for micromotion feature extraction.
In fact, the existing micromotion feature extraction methods of space target are implemented after target detection and tracking.It is necessary to allocate the limited radar resources for target detecting, tracking, and feature extraction, successively.In the case of multitarget monitoring, the allocation contradiction of radar resources will be serious.To overcome this problem, we intend to seek a novel approach which can implement target detecting and tracking simultaneously when extracting the micromotion features of a space target, which can be available for both the fixed scatterer model and the slide scatterer model.
In recent years, the track-before-detect (TBD) technology has shown good performance in weak target detecting and tracking.Unlike traditional techniques that declare the presence of a target at each scan, the TBD technology processes more consecutive scans jointly, and then it declares the presence of a target and its track [16][17][18].Through the interscan accumulation, the TBD technology can improve the probability of target detection.The idea of establishing a feedback loop between feature extraction and TBD of target has been provided in our previous preliminary work [19].In this paper, a micromotion feature extraction method of space target based on TBD technology is further proposed.In the method, the tracks of target scatterers can be obtained during the TBD of target, and the tracks are the range-slow-time curves of target scatterers indeed.We call the range-slowtime curve as "range trajectory" in this paper.On this basis, we attempt to add the micromotion feature extraction into the process of target detecting and tracking by establishing a feedback loop, namely, that the micromotion feature parameters are extracted via fitting the obtained range trajectory according to the mathematical expression of m-D effect, and in turn the extracted micromotion feature parameters are utilized to update the parameters of TBD adaptively.As a result, the micromotion feature extraction, detecting, and tracking of space target can be implemented simultaneously with the information feedback, which can provide real-time and effective information for target recognition.By changing the mathematical expressions of different kinds of m-D effects for range trajectory fitting, the fixed scatterer model, the slide scatterer model, and any kinds of micromotion forms can be processed with the proposed method.
This paper is organized as follows.Taking the conical target contains a fixed scatterer and two slide scatterers as an example, the micromotion feature of the target precession is analyzed in Section 2. Combining with the TBD technology, the micromotion feature extraction method is presented in detail in Section 3. Simulations are presented in Section 4 and some conclusions are made in the last section.

Micromotion Feature of Conical Target with Precession
When the space target is flying outside the atmosphere, it needs to spin around its own axis of symmetry to maintain stability.Meanwhile, it needs to cone around a spatial directional axis due to the lateral force caused by projectile separation and bait release.These two types of rotation formed the precession of target.For a conical target, it usually consists of three dominant scatterers (i.e., one tipcone scatterer and two cone-base scatterers), where the tipcone scatterer can be treated as a fixed tip-cone scatterer, and the cone-base scatterers are slide cone-base scatterers.
The geometry of conical target with precession is shown in Figure 1, where  is the target center, and the target is rotating around -axis with precession angle  and precession frequency  rad/s.The distance between the radar and the target center is denoted as   .The radar LOS is denoted as  ⇀ n , -axis is in the plane determined by  ⇀ n and -axis, and -axis is established according to the right-hand rule.The LOS can be represented as  ⇀ n = sin   → y + cos   → z , where  is the angle between LOS and -axis.The target consists of three scatterers denoted as , , and , respectively.The radius of the base circle is  0 , the distance between the target centroid and scatterer  is denoted as ||, and the vertical distance between the target centroid and the bottom circle is represented as ||.According to the rigid body dynamics, the velocity of target should be along -axis and the velocity is given as V.
Assume that the angle between -axis and the projection of vector  → oa on the  plane is  0 at time  = 0, then the angle between LOS and  → oa is denoted as , and it can be calculated by cos  = sin  sin  sin ( +  0 ) + cos  cos .
The projection of  → oa,  → oc, and  → od in the line-of-sight direction at time  can be represented as   (),   (), and   (), respectively; that is, Reference [1] points out that, in the wideband radar system, the micromotion feature of target can be described by range-slow-time image, where the peak of range profile appears to be a range-slow-time curve (range trajectory) which is determined by the scatterer range ().And the range trajectory reflects the micromotion feature of target scatterer.Therefore, from (2)-( 4), the micromotion feature of cone-tip scatterer  is sinusoid, while that of the cone-base scatterers  and  is quasisinusoid, which deviates from the standard sinusoid.
The impact of the shielded effect on the proposed micromotion feature extraction method will be discussed in Section 3.
Next, a simple simulation is given to validate the correctness of theoretical analysis.The pulse duration is   = 1 s, carrier frequency is   = 10 GHz, signal bandwidth is  = 3 GHz, pulse repetition frequencies is PRF = 60 Hz, and the coherent processing time is   = 1 s.The radius of the base circle is  0 = 1 m, the distance between the target center and cone-tip scatterer  is || = 3 m, and the vertical distance between the target center and the bottom circle is || = 0.3 m.
To observe the micromotion feature of target more intuitively, we assume that   = 0 m and V = 0 m/s.The precession angle is  = 15 ∘ , precession frequency is  = 8 rad/s, the angle between LOS and -axis is  = 145 ∘ , and the initial angle is  0 = /100 rad.We can obtain the range-slow-time image as shown in Figure 2.
From Figure 2, we can see that the range profile peak of scatterer  changes with slow-time following the sinusoidal form, while those of the sliding scatterers  and  are quasisinusoids (deviating from the standard sinusoid form), which confirms the theoretical analysis.

Micromotion Feature Extraction Based on TBD
Due to the requirement of additional observation pulses after target detection and tracking, the existing micromotion feature extraction methods usually occupy much radar resources.Also, the real-time performance and radar efficiency are not satisfied.So we intend to combine the micromotion feature extraction with target detecting and tracking in the way of information feedback and then implement the micromotion feature extraction, detecting and tracking of space target simultaneously.It will save the radar resources and improve the radar efficiency and real-time performance of micromotion feature extraction.
In this section, we establish a feedback loop between micromotion feature extraction and TBD of target.Firstly, based on TBD technology, the target scatterer range trajectory information is backtracked along with the energy accumulation process, and then the micromotion feature parameters of target can be fitted with these obtained trajectory information.In return, the extracted micromotion feature parameters are used to update the state transition set of TBD adaptively, and the results of micromotion feature extraction are considered into the declaration of the target presence.As a result, micromotion feature extraction, target detecting, and tracking can be completed simultaneously.

Range Trajectory Backtracking and Micromotion Feature
Extraction.For simpleness, the range trajectory backtracking for single scatterer based on TBD is described firstly.The monitoring area is divided into   ×   grids according to the range and azimuthal angle, and each grid is denoted as a state (, ),  = 1, . . .,   ,  = 1, . . .,   , which represents the position (( −   /2) ⋅ Δ +  0 , ( −   /2)Δ +  0 ), where Δ and Δ are the stepped increasement of the range and azimuthal angle, and ( 0 ,  0 ) is the center of the monitoring area.Assume the radar transmits wideband signal ( t ).At each scan, the beams towards all the azimuthal angles (  angles) are formed by the radar.The beam width is denoted as   , each state (, ) will be hit by  + 1 successive beams, where Obviously, the sequence number of these +1 beams should be ( − /2,  − /2 + 1, . . .,  + /2).At the th scan, the echo signal of the state x  = (  ,   ) of the th beam can be represented as where   (⋅) represents the beampattern of the th beam and   (  ) represents the obtained transmit gain at the   th azimuthal angle with the th beam and  (  ,  ) is the backscattered amplitude of the state (  ,   ) at the th scan.
After performing range compression, the high-resolution range profile (HRRP) can be obtained as Thus, at the th scan, the measured value of each state x  = (  ,   ) can be defined as At the th scan, the cumulative energy of state x  is denoted as (x  ), and it can be calculated as where Γ(x  ) is the state transition set.Γ(x  ) contains all the possible state x −1 which can transit to state x  .Set an appropriate detection threshold for the cumulative energy function (x  ) after  scans accumulation, and the state sequence whose cumulative energy is larger than the threshold can be backtracked according to where   (x  ) is the backtracking function, which is used to record the state corresponding to the maximum cumulative energy at each scan.Assume that the recorded state of the th scan is x  = (  ,   ), the estimated scatterer range trajectory can be denoted as () = (  −  /2)⋅Δ+ 0 ,  = 1, . . ., , and the angle trajectory is Θ() = (  −  /2)Δ+ 0 ,  = 1, . . ., .
For the cone-tip scatterer , the estimated range trajectory should be equal to (2), where Δ is the time interval between the two adjacent scans.Without loss of generality, it holds Δ = 1/PRF.Similarly, for the cone-base scatterers  and , the estimated range trajectories () should be equal to ( 3) and ( 4), respectively.However, the TBD method described above can obtain only one scatterer range trajectory, which cannot meet the requirement of getting range trajectories of each scatterers.Although some TBD methods for multitarget detecting and tracking have been proposed [20,21], they require that the state of different targets can not be the same.However, in this paper, the different scatterers are usually at the same azimuthal angle and may be of the same range at some scans (i.e., the range trajectories may be intersected).That means the different scatterers may have the same state, which leads to the existing method which can not be used.To resolve this problem, we improve the TBD method as follows.
Assume the target consists of  observable scatterers.At the th scan, for each azimuthal angle   , select any  states to form an expanded state (( ,1 ,   ), ( ,2 ,   ), . . ., ( , ,   )), where the  states can be the same; that is,  ,1 can be equal to  ,2 when  1 ̸ =  2 .The measured value of the expanded state y  = {( ,1 ,   ), ( ,2 ,   ), . . ., ( , ,   )} is defined as On the basis, the cumulative energy shown as ( 9) is conducted in terms of the expanded state, where the state correlation is necessary.For example, there are 6 kinds of states correlation approaches when  = 3, shown as Figure 3.
Among these 6 kinds of states correlation approaches, only one is consistent with the actual scatterers trajectories.In [22], we have proposed a state correlation approach which can be used for states correlation.
After  scans, the cumulative energy (y  ) should be compared with a detection threshold   : where  is a constant which will affect the target detecting performance, in this paper we call  as "detection threshold coefficient."The state sequences whose cumulative energy is larger than the threshold can be backtracked according to (10).
It should be pointed out that the states correlation leads to the dependence of the energy accumulation of each scatterer.A specific expression of the detection threshold coefficient is difficult to be derived with a given false alarm probability.Just as [20], the detection threshold coefficient can be chosen from Monte-Carlo experiments.
Assume that the number of the cumulative energy (y  ) which is larger than the threshold is .Thus,  state sequences will be obtained, and each state sequence contains  range trajectories   () = (  −   /2) ⋅ Δ +  0 ( = 1, . . ., ,  = 1, . . ., ) and  angle trajectories Θ  () = (  −   /2)Δ +  0 ( = 1, . . ., ,  = 1, . . ., ).Due to the fact that intersections number of any two range trajectories will be small, a condition for selecting the reasonable state sequence from the  state sequences is defined as where ( 1 () =  2 ()) represents the intersections number of  1 () and  2 () and  is a constant used to control the ratio of intersections number to scans number .The state sequences which satisfy the condition shown as ( 14) are selected, and the corresponding range trajectories and angle trajectories are obtained from the target scatterers.
For each possible value of , the TBD procedure proposed above can be conducted to obtain the range trajectories of  scatterers.The shielded effect has been discussed in Section 2: when  falls with in Area 1 or Area 3, three scatterers are observable (corresponding to  = 3); when  falls within Area 2 or Area 4, two scatterers are observable (corresponding to  = 2).Therefore, the possible values of  are  = 3 and  = 2. Thus, the TBD procedure of  = 3 and  = 2 should be conducted.
In the case of  = 3, if there is state sequence which satisfies (14) which can be selected out, three range trajectories (i.e.,  1 (),  2 (), and  3 ()) can be obtained.Firstly, we assume that  1 () is the range trajectory of the cone-tip scatterer  and  2 () and  3 () are the range trajectories of the cone-base scatterers  and , respectively. 1 (),  2 (), and  3 () can be represented as where () represents the error induced from the noise and TBD procedure.We can estimate the micromotion feature parameter vector P A = [||,  0 , , , , V, ||,  0 ,   ] by fitting  1 (),  2 (), and  3 () according to the curve form of standard sinusoid shown as (2) and quasisinusoid shown as (3)-( 4) with the least squares method: min Equation ( 16) can be solved with the Levenberg-Marquardt method [23], which is sensitive to the initial values.Therefore, how to set the appropriate initial values is proposed as follows.The EMD method [24] can be used to separate  1 () into a set of intrinsic mode functions (IMF) which is descended by frequency.Thus, we can obtain where  is the error from () and EMD method.According to (3) and ( 4), we can get Similarly, separating  + () with the EMD method, we can obtain The initial values || ini ,  0ini ,  ini ,  ini ,  ini , V ini , || ini ,  0ini , and  ini for (19) can be obtained according the equation set: where  sin(⋅) is the inverse function of sin (⋅).In ( 20 All (  ) is compared with each other, and the estimated micromotion feature parameter vector is defined as The fitting error is denoted as It should be pointed that PA and Ê are obtained under the assumption that  1 () is the range trajectory of the cone-tip scatterer .Therefore, we rewrite them as PA ( 1 ) and Ê( 1 ).
Next, we assume  2 () and  3 () are the range trajectory of the cone-tip scatterer, respectively.The micromotion feature parameter vectors PA ( 2 ) and PA ( 3 ) can be obtained, and the corresponding fitting errors are denoted as Ê( 2 ) and Ê( 3 ).The micromotion feature parameter vector corresponding to the minimum fitting error is selected out as the final micromotion feature parameter vector The fitting error is denoted as Similarly, in the case of  = 2, the micromotion feature parameter vector can be obtained.What should be pointed out is that if the two range trajectories are corresponding to the cone-base scatterers  and , respectively, the parameter || can not be obtained.

Adaptive Update of State Transition
Set.Just as mentioned in (18), the state transition set contains all the possible state y −1 which can transit to state y  , and the construction of state transition set will significantly affect the efficiency and performance of TBD and micromotion feature extraction.In this paper, the state transition set is updated adaptively, by forecasting the state of scatterers in next scan according to the extracted micromotion feature parameters.
Assume that, in the th scan, the cumulative energy (y  ) is larger than the threshold; the extracted micromotion feature parameter vector can be obtained with the method proposed in Section 3.
where mean (⋅) represents the mean value of a vector.Obviously, the smaller  ,y  is, the higher consistency is.Obviously, if  ,y  is relatively small, it shows that the consistency of the extracted micromotion feature parameter vectors is high, we can consider that the precision of the micromotion feature extraction is well, and the forecast of state information in the next scan is accurate.In this case, the searching range gate can be decreased to reduce the computation load.On the contrary, if  ,y  is relatively large, the forecasted precision is low and the searching range gate should be increased appropriately.Therefore,   (⋅) should be an increasing function, which is defined as where  is a constant coefficient,   is the minimum of searching range gate, which is set to be 5 in this paper, and the maximum value of   (⋅) is limited to 20.

Adaptive Starting and Ending of Precession Feature Extraction.
In TBD algorithms, the detection threshold affects the target detecting performance directly.However, how to set an appropriate detection threshold is difficult, especially in the multiscatterers environment.
As the micromotion feature extraction has been included into the process of target detecting and tracking by establishing a feedback loop, the pertinent issue now is how to declare the presence of a target and when to start and end the micromotion feature extraction algorithm.
In the proposed method, the energy accumulation value and the result of micromotion feature extraction are both taken into full consideration to address the problem above, which is different from the traditional signal processing method.Just as we know, in traditional method, only when the energy accumulation value is larger than the threshold, it declares the presence of a target, and then the extraction of its micromotion feature will be implemented by transmitting additional observation pulses.In comparison, in the proposed method, by adding the micromotion feature extraction into the process of target detecting and tracking, on condition that the precision of micromotion feature extraction is high enough, it can declare the presence of a target and the micromotion feature is extracted successfully at the same time, although the energy accumulation value dose not achieve the required level.
Assume the minimum and maximum total number of scans that are jointly processed in TBD are   and   , respectively.Firstly, set two detection threshold coefficients: the lower detection threshold coefficient  1 and the higher detection threshold coefficient  2 .In the th ( ≥   ) scan if the cumulative energy (y  ) is larger than   2 =  2 ⋅  ⋅ , we declare the presence of a target and the micromotion feature parameters can be obtained with the proposed method described in Section 3.1.On the other hand, if the cumulative energy (y  ) is larger than   1 =  1 ⋅  ⋅  and smaller than   2 , backtracking the target scatterers trajectories according to (10) and (14).On this basis, the extracted micromotion feature parameter vector at th scan can be obtained.Go on to cumulate energy for the data of  + 1th scan, and extract the micromotion feature parameters at  + 1th scan.The consistency of the extracted micromotion feature parameter vectors  ,y  can be calculated according to (29).If  ,y  and fitting error Ê,y  are both relatively small (satisfies  ,y  <   and Ê,y  <   ), we can declare the presence of a target and get the micromotion feature parameters.The energy accumulation is no longer needed.Otherwise, update the state transition set according to (27), and the energy accumulation of the data of  + 2th scan is needed.Repeat the steps described above until it satisfies  ,y  <   and Ê,y  <   or (y  ) >   2 , or it reached the   th scan.In conclusion, the flowchart of micromotion feature extraction of space target based on TBD is shown in Figure 4.
There are three points should be noticed for the proposed method.
(1) The values of  1 and  2 are chosen from Monte-Carlo experiments with the given false alarm probability  FA, 1 and  FA, 2 in the case the presence of a target is declared when the energy accumulation value is larger than the threshold.The corresponding detection probabilities are denoted as  , 1 and  , 2 .However, the result of micromotion feature extraction is taken into consideration to declare the presence of a target in this paper; thus the final false alarm probability  FA and detection probability   obtained from the proposed method are hard to be calculated from a specific expression with the variables of  FA, 1 ,  , 1 ,  FA, 2 , and  , 2 .Numerous experiments have shown that  FA will be a little higher than  FA, 2 .Therefore the value of  FA, 2 can be chosen according to the desired false alarm probability, which is set as  FA, 2 = 0.005 in this paper.The value of  FA, 1 will affect the detection probability and computation load, and it is set as  FA, 1 = 0.5 in this paper from numerous experiments with the consideration of maximizing the detection probability and minimizing the computation load.
(2) When the value of  is larger than the observable scatterers number, the cumulative energy and the trajectories which are corresponding to the target will not be different from that corresponding to the noise.Thus, the target will not be detected.On the contrary, if the value of  is smaller than the observable scatterers number, only  trajectories can be obtained.Thus, some scatterers trajectories will be lost and the corresponding micromotion feature parameters can not be extracted.Therefore, the value of  should be set as the maximum possible value firstly, which can be preset by experience.Then the value of  decreases gradually until a target is detected and the micromotion feature is extracted or  = 1.
(3) The computational load of the proposed method is proportional to the search range.Therefore, to reduce the computational load, we assume that target detecting and tracking with narrow-band radar is conducted firstly, and a lower threshold is used to detect the possible targets; then the obtained coarse position and velocity information of targets are utilized to determine the search range of the proposed method.In the following text, for the state of y  , the cumulative energy, range trajectories, angle trajectories, estimated micromotion feature parameter vector, fitting error, consistency of the extracted micromotion feature parameter vector, state transition set of target, and noise in next scan are denoted as   (y  ),   (y  ),  y  , ,  y  , , Θ y  , , Θ y  , , PA,y  , , PA,y  , , Ê,y  , , Ê,y  , ,  ,y  , ,  ,y  , , Γ  (y +1 ), and Γ  (y +1 ), respectively.
Using any two values of the range trajectory, the probable slope can be obtained, which can be utilized for the coarse slope compensation of trajectories, as shown in Figure 8.
Based on the range trajectories of target and noise, respectively, the estimated micromotion feature parameter vectors PA,y  , and PA,y  , can be obtained according to (13)- (22).At this time, the energy accumulation values   (y  ) and   (y  ) are both smaller than threshold   2 , and the consistencies  ,y  , and  ,y  , can not be obtained due to the inexistence of the extracted micromotion feature parameter vector at the former scan.Therefore, it is necessary to continue to carry out energy accumulation.For   (y  ) and   (y  ), respectively, update state transition sets Γ  (y +1 )  and Γ  (y +1 ) according to (27) by utilizing the obtained precession feature parameter vectors PA,y  , and PA,y  , .After 24 scans, for   (y  ), we can obtain the consistency  ,y  , = 3.92 and fitting error Ê,y  , = 0.26, which does not satisfy  ,y  <   and Ê,y  <   .On the contrary, for   (y  ), the consistency  ,y  , = 0.08 and fitting error Ê,y  , = 0.05, which satisfies  ,y  <   and Ê,y  <   .Therefore, we declare the presence of a target, and the micromotion feature parameter vector PA,y  , can be obtained, as shown in Table 1.And the update process of the consistency, fitting error, and the size of searching range gate are shown in Figure 9.
From Figure 8, we can see that the consistency of the extracted micromotion feature parameter vector is increasing (i.e.,  ,y  is decreasing) with energy accumulation, and  ,y  , can reduce to the value smaller than   = 0.1, while  ,y  , is always higher than   = 0.1.Similarly, the size  of searching range gate has the same change trend with the consistency.At the same time, the fitting error increases with energy accumulation, and Ê,y  , can stabilize at the value which is smaller than   , while Ê,y  , will be higher than   .Therefore, the state sequence corresponding to the target can be selected out, and the other state sequences induced from noise can be eliminated.Furthermore, the obtained micromotion feature parameter vector shown in Table 1 is closed to the theoretical value, which illustrates the effectiveness of the proposed method.Next, the detection performance and micromotion feature extraction performance of the proposed method are compared with the traditional TBD method and the extended Hough transform method, respectively.With the proposed method, the false alarm probability and the detection probability are  FA = 0.0052 and   = 0.8823.With the traditional TBD method, fixing the false alarm probability as 0.0052, the detection probability is 0.4827.Obviously, the detection performance of the proposed method is better than that of the traditional TBD method.Further, the micromotion feature extraction performance is defined as the estimated error of each micromotion feature parameter.With the extended Hough transform method, the micromotion feature parameter vector is obtained shown in Table 2. Due to the low SNR, the micromotion feature parameter vector can not be extracted effectively with the extended Hough transform method, while the high micromotion feature extraction performance can be obtained with the proposed method, shown as Table 1.
Finally, the detection performance and micromotion feature extraction performance with different SNRs are considered.Fixing the false alarm probability as 0.0052, the detection probability of the proposed method is better than that of the traditional TBD method, shown in Figure 10.Without loss of generality, we think the micromotion feature extraction is successful when the estimated error of each micromotion feature parameter is less than 10%.The success of micromotion feature extraction with different SNRs is shown in Figure 11.We can see that the success of micromotion feature extraction with the proposed method is much higher than that with the extended Hough transform method when SNR is lower than 11 dB.What is more important is that the micromotion feature parameters can be extracted without transmitting additional pulses, which can save the radar resources and provide real-time information for target recognition.

Conclusions
A micromotion feature extraction method for space target based on TBD is proposed in this paper.This method establishes a feedback loop between micromotion feature extraction and TBD of target.As a result, the micromotion feature can be extracted concurrent with target detecting and tracking, and the detection performance can be improved.The steps of the method are described in detail, and some simulations are given to illustrate its effectiveness.It is noted that the cone-shaped target is taken as an example

Figure 1 :
Figure 1: The geometry of conical target with precession.

Figure 4 :Figure 5 :
Figure 4: Micromotion feature extraction of space target based on TBD.

Figure 8 :
Figure 8: Backtracked range trajectories with coarse slope compensation; (a) induced by target; (b) induced by noise.

Figure 9 :
Figure 9: Update process.(a) Consistency of the extracted micromotion feature parameters.(b) Fitting error.(c) Size of searching range gate.

Table 2 :
Extracted micromotion feature parameter vector in with extended Hough transform method.