Micro-Doppler Effect and Sparse Representation Analysis of Underwater Targets

At present, the micro-Doppler effects of underwater targets is a challenging new research problem. This paper studies the micro-Doppler effect of underwater targets, analyzes the moving characteristics of underwater micro-motion components, establishes echo models of harmonic vibration points and plane and rotating propellers, and reveals the complex modulation laws of the micro-Doppler effect. In addition, since an echo is a multi-component signal superposed by multiple modulated signals, this paper provides a sparse reconstruction method combined with time–frequency distributions and realizes signal separation and time–frequency analysis. A MicroDopplerlet time–frequency atomic dictionary, matching the complex modulated form of echoes, is designed, which effectively realizes the concise representation of echoes and a micro-Doppler effect analysis. Meanwhile, the needed micro-motion parameter information for underwater signal detection and recognition is extracted.


Instruction
Victor C. Chen introduced the concept of micro-Doppler effects in radar [1].The micro-Doppler effect is a physical phenomenon caused by micro-motions, such as vibrations and rotations of target components, which are inherent features of moving targets.Extracting micro-Doppler features is of great significance for improving the performance of target detection and recognition.At present, there are many studies on this topic in the radar and air sonar fields [2][3][4].However, in the field of underwater acoustic signal processing, it is only in the beginning stages [5,6].When a large underwater target, such as a submarine, unmanned underwater vehicle, or submersible, is underway, its components will vibrate, rotate, and perform other complex periodic micro-motions due to the action of forces from the engine, propeller, and other mechanical equipment, resulting in a micro-Doppler effect.The micro-Doppler effects of underwater targets are more complicated, due to the influence of an inhomogeneous water medium and turbulence on the vibration of the target shell and the rotation of the propeller.
The magnitude of micro-Doppler frequency shifts is directly proportional to the product of the vibration rate and displacement offset and inversely proportional to the emission signal wavelength [1].The Doppler modulation of vibration can be observed if the product is large enough and the signal wavelength is small enough.
Micro-motions produce complex modulation of sonar emission signals, and the target echo is a nonstationary time-varying signal.The main theories of nonstationary signal analysis are joint time-frequency distribution [7] and instantaneous frequency analysis [8].Different kernel functions can be used to analyze the local joint time-frequency information of multi-component signals successfully, but the inherent cross-term problem of quadratic Sensors 2023, 23, 8066 2 of 15 time-frequency distribution has not been solved.Instantaneous frequency analysis can obtain micro-Doppler frequency estimation, but the analytical signal obtained by the Hilbert transform is only suitable for analyzing single-component signals.For multi-component signals superposed by the reflected waves from target components, determining how to analyze and interpret the modulation components contained in the echo and extracting the micro-Doppler features is still a challenge.
In recent years, the sparse representation theory of signals has attracted broad attention [9,10], and has achieved good application effects in many aspects of signal processing and deep learning [11][12][13][14].Mallat and Zhang proposed a new time-frequency distribution method based on the Wigner Ville distribution [15].The method constructed a time-frequency distribution without interatomic cross-term interference and with good time-frequency aggregation.This time-frequency distribution can adaptively adjust the characteristic parameters of time-frequency atoms with the local characteristics of the analyzed signal.Therefore, the time-frequency distribution characteristics of multicomponent time-varying signals can be revealed, which provides a method for analyzing micro-Doppler effects.
In this paper, the micro-Doppler effects of underwater targets are studied.According to the micro-motion analysis results of underwater targets, the echo models of some typical micro-motion components are established.A time-frequency atomic dictionary, the MicroDopplerlet dictionary, is designed.The modulation law and micro-Doppler effects of echoes are analyzed by sparse reconstruction combined with time-frequency distributions.

Micro-Motion Analysis of Underwater Targets
The vibrating shell and rotating propeller are the main generators of the micro-motions of underwater targets.The magnitude of the micro-motion frequency and displacement are the basis of micro-Doppler effects.The propeller blade is long and has a very high tip velocity, causing a large amount of velocity and displacement [16].Table 1 calculates the micro-Doppler frequency shift of 5-blade and 7-blade propellers and shells under typical conditions.The results show that propellers have obvious micro-Doppler effects.
The shell of underwater targets can be classified as a single-or double-layer cylindrical shell.Aron correctly described the vibration problem of cylindrical shells for the first time [17].Flugge [18] and Love [19] proposed methods to calculate the vibration of a thin shell with an infinite length and fixed length, establishing the famous Love first approximation theory and obtaining the calculation formula for the free vibration frequency of an infinite-length cylindrical shell.Subsequently, based on the Love theory, many scholars have developed various theories for approximate cylindrical shells [20][21][22].The results show that low frequencies correspond to the large radial displacement, middle frequencies correspond to the axial displacement component, and high frequencies correspond to the circumferential displacement.Yamada [23] used the transfer matrix method to calculate the free vibration of a double-layer cylindrical shell.Raheb and Wagner [24] established the mechanical equilibrium equation of a double-layer stiffened shell.The force of a stiffened double-layer cylindrical shell is expressed in the form of shell displacement, which provides a theory and method for the calculation of shell displacement.
On the other hand, Tang et al. conducted theoretical research on the vibration characteristics of ring-ribbed cylindrical shells under different boundary conditions and discussed the influence of different parameter changes on the natural frequency of shell structures [25].Shang et al. used ANSYS software to obtain the normal displacement, normal vibration velocity, surface sound pressure, normal sound intensity, and other parameters of the outer surface of double-layer cylindrical shells [26].The submarine is usually simplified to the model shown in Figure 1.Wang [27] carried out experiments on the vibration characteristics of different parts of a simplified model of a single-layer shell.When a 500 N external force was applied to the shell, the vibration displacement of the hemisphere shell A was 3.98 microns, the vibration displacement of the middle cylindrical shell B was 8.91 microns, and the vibration displacement of the conical shell C was 12.59 microns.Huang [28] con-ducted a numerical simulation experiment on the vibration characteristics of double-layer shells.When a 1000 N external force was applied to the shell, the vibration displacement of the hemisphere shell A was 5.62 microns, the vibration displacement of the middle cylindrical shell B was 7.94 microns, and the vibration displacement of the conical shell C was 3.16 microns.The shell of underwater targets can be classified as a single-or double-layer cylindrical shell.Aron correctly described the vibration problem of cylindrical shells for the first time [17].Flugge [18] and Love [19] proposed methods to calculate the vibration of a thin shell with an infinite length and fixed length, establishing the famous Love first approximation theory and obtaining the calculation formula for the free vibration frequency of an infinite-length cylindrical shell.Subsequently, based on the Love theory, many scholars have developed various theories for approximate cylindrical shells [20][21][22].The results show that low frequencies correspond to the large radial displacement, middle frequencies correspond to the axial displacement component, and high frequencies correspond to the circumferential displacement.Yamada [23] used the transfer matrix method to calculate the free vibration of a double-layer cylindrical shell.Raheb and Wagner [24] established the mechanical equilibrium equation of a double-layer stiffened shell.The force of a stiffened double-layer cylindrical shell is expressed in the form of shell displacement, which provides a theory and method for the calculation of shell displacement.
On the other hand, Tang et al. conducted theoretical research on the vibration characteristics of ring-ribbed cylindrical shells under different boundary conditions and discussed the influence of different parameter changes on the natural frequency of shell structures [25].Shang et al. used ANSYS software to obtain the normal displacement, normal vibration velocity, surface sound pressure, normal sound intensity, and other parameters of the outer surface of double-layer cylindrical shells [26].The submarine is usually simplified to the model shown in Figure 1.Wang [27] carried out experiments on the vibration characteristics of different parts of a simplified model of a single-layer shell.When a 500 N external force was applied to the shell, the vibration displacement of the hemisphere shell A was 3.98 microns, the vibration displacement of the middle cylindrical shell B was 8.91 microns, and the vibration displacement of the conical shell C was 12.59 microns.Huang [28] conducted a numerical simulation experiment on the vibration characteristics of double-layer shells.When a 1000 N external force was applied to the shell, the vibration displacement of the hemisphere shell A was 5.62 microns, the vibration displacement of the middle cylindrical shell B was 7.94 microns, and the vibration displacement of the conical shell C was 3.16 microns.The propeller and shafting-shell coupled system have an important influence on the shell vibration.The unsteady external excitation force of the propeller during operation and the friction force of the shafting and bearing caused by the propeller rotation are important excitation sources for the vibration of the shell structure.At present, scholars have paid attention to the exciting force of the propeller on the shell [29][30][31][32][33]. Hassan et al. conducted an experimental test on the hydrodynamic performance of a 5-blade high-skew propeller [16].The results showed that the propeller's excitation force could reach 158,750 N when the advance coefficient was 0.3, and the rotational velocity was 600 r/m.This thrust force is amplified by shafting and will cause the shell vibration greatly.The research results of reference [34] showed that the unsteady excitation force of the propeller can be amplified at the natural vibration frequency, and the excitation force transmitted to the shell at the first natural frequency can be amplified nearly 20 times.Zhang [35] established the coupled shafting-hull structure model and obtained an unsteady excitation force from a 7-blade propeller with a rotating velocity of 110 r/m using a semi-empirical formula, among which the longitudinal unsteady excitation force was about 700 N.
Based on the above research results, if the propeller excitation force is known, the excitation force transmitted to the shell can be obtained through the shafting amplification characteristics.Then, according to the displacement generated by the unit excitation force, the vibration displacement of the shell can be calculated.
The estimation results of the micro-motion amplitude, micro-motion frequency, and micro-Doppler frequency of propellers and different components of shells are given in Table 1.When the wavelength is 50 mm, according to Equations ( 11) and ( 16), the calculated micro-Doppler frequency is above 1 kHz.It can be seen that the propeller has a high velocity, large displacement, and significant micro-motion characteristics, and the micro-Doppler frequency shift is large.The vibration amplitude of the submarine shell can reach a centimeter level, and the micro-Doppler frequency shift can be hundreds of hertz.

Micro-Doppler Effect and Echo Models of Typical Components
Generally speaking, micro-Doppler effects are not only related to the wavelength and the product of micro-motion frequency and displacement but also to the observation direction of sonar.When studying micro-Doppler effects, the included angle between the vibration direction and the incident wave direction should be considered to establish the echo model.However, for simple harmonic vibration points and planes and rotation propellers, the simplified models based on ignoring the included angle are still universal.Simplified models can describe the complex modulation law of micro-Doppler effects.The simplified echo models of these three types are established below.
Based on the comprehensive research results, the calculation formula for vibration amplitude is: where F is the applied propeller excitation force, ρ is the vibration displacement generated by a unit Newtonian force, and K is the amplification factor of the excitation force transmitted from the first natural frequency of the propeller to the shell.

Scattering Point of Simple Harmonic Vibration
Smaller target components can be regarded as scattering points.When a scattering point vibrates, the echo phase is modulated, and the frequency changes periodically.The sideband frequency near the carrier frequency is generated, the micro-Doppler effect.The geometric relationship between vibration point and sonar is established as shown in Figure 2. The target coordinate system x y z is the translation of the sonar coordinate system xyz along the z axis.Supposing that the scattering point t P vibrates along the z axis, and the equilib- A is the amplitude, v f is the vibration frequency, and is the angular frequency.The simple harmonic vibration of scattering points can be regarded as a spring oscillator system, and the displacement law is obtained from kinematics: sin( ) where ϕ is the initial phase.The displacement is a sine function varying with time t .
The distance from the scattering point to the sonar is: be the sonar emission signal.Where 0 f is the carrier frequency and u  is the complex envelope.Under ideal conditions, the echo is the delay of the emis- sion signal, that is: where σ is the reflection coefficients.
is the phase function, and the time derivative of ( ) t Φ is instantaneous frequency, which is: The phase modulation causes micro-Doppler effects.The micro-Doppler frequency is:

Scattering Plane of Simple Harmonic Vibrations
The large area shell vibration can be regarded as a vibration scattering plane.The amplitude and phase of echoes are all modulated.The amplitude flash and the frequency change periodically.The geometrical relationship between the vibration plane and the observation sonar is shown in Figure 3.The center of the vibration plane is located at 0 O (0, 0, ) z ′ , perpendicular to the z axis.The scattering points on the plane vibrate along the z axis.The vibration angular frequency v ω is the same, but the amplitude is differ- ent.The amplitude distribution function is ( , )  Supposing that the scattering point P t vibrates along the z axis, and the equilibrium position is O (0, 0, z 0 ), A v is the amplitude, f v is the vibration frequency, and ω v = 2π f v is the angular frequency.The simple harmonic vibration of scattering points can be regarded as a spring oscillator system, and the displacement law is obtained from kinematics: where φ is the initial phase.The displacement is a sine function varying with time t.The distance from the scattering point to the sonar is: Let s(t) = ue j2π f 0 t be the sonar emission signal.Where f 0 is the carrier frequency and u is the complex envelope.Under ideal conditions, the echo is the delay of the emission signal, that is: where σ is the reflection coefficients.Φ(t) = 2π f 0 t − 4π λ A m sin(ω v t + ϕ) is the phase function, and the time derivative of Φ(t) is instantaneous frequency, which is: The phase modulation causes micro-Doppler effects.The micro-Doppler frequency is:

Scattering Plane of Simple Harmonic Vibrations
The large area shell vibration can be regarded as a vibration scattering plane.The amplitude and phase of echoes are all modulated.The amplitude flash and the frequency change periodically.The geometrical relationship between the vibration plane and the observation sonar is shown in Figure 3.The center of the vibration plane is located at O (0, 0, z 0 ), perpendicular to the z axis.The scattering points on the plane vibrate along the z axis.The vibration angular frequency ω v is the same, but the amplitude is different.The amplitude distribution function is A v (x , y ).The echo model can be established by analyzing the displacement law of simple harmonic points.At time t , the displacement of any scattering point ( , ,0) t P x y ′ ′ on the vibration plane away from the equilibrium position is ( , )sin The distance from the sonar to the scattering point is: Substituting Equation ( 6) into t z′ in the case of t z z ′ >> 0 , the above formula becomes approximately: The reflected wave of t P can be expressed as: A double integral of the reflected wave is carried out over the whole plane.The echo is: where S is the integral region and ds is the area element.For a circular plane of radius R , let the amplitude distribution function be . After integration, the echo is: Due to the existence of the sinc function, the amplitude presents the characteristics of a time domain flash.The phase is modulated by a sinusoidal function.The micro-Doppler frequency is:

Blades of Rotation Propeller
Similar to the vibration plane, when the propeller rotates, the blade echo is modulated in both amplitude and phase.The geometric relationship between blades and sonar is shown in Figure 4. Let the blade's length be L , the center is located at , and the blades rotate at an angular velocity of v ω .At time t, the displacement of any scattering point P t (x , y , 0) on the vibration plane away from the equilibrium position is z t = A v (x , y ) sin ω v t.The distance from the sonar to the scattering point is: Substituting Equation ( 6) into z t in the case of z 0 >> z t , the above formula becomes approximately: The reflected wave of P t can be expressed as: A double integral of the reflected wave is carried out over the whole plane.The echo is: where S is the integral region and ds is the area element.For a circular plane of radius R, let the amplitude distribution function be After integration, the echo is: Due to the existence of the sinc function, the amplitude presents the characteristics of a time domain flash.The phase is modulated by a sinusoidal function.The micro-Doppler frequency is:

Blades of Rotation Propeller
Similar to the vibration plane, when the propeller rotates, the blade echo is modulated in both amplitude and phase.The geometric relationship between blades and sonar is shown in Figure 4. Let the blade's length be L, the center is located at O (0, 0, z 0 ), and the blades rotate at an angular velocity of ω v .
Sensors 2023, 23 Since 0 p z l >> , the distance is approximately: The reflected wave of scattering point t P is: ( ) The reflected wave is integrated into the blade length L , and the echo of the blade is: The micro-Doppler frequency is: The size of the micro-Doppler frequency is related to the blade length, rotation angular velocity, and wavelength.The echo amplitude flash and micro-Doppler shift are unique characteristics of underwater targets in the time and frequency domains.Micro-Doppler features are considered important features of underwater target recognition.The flash characteristics can be observed in the time domain, while the micro-Doppler frequency shift requires joint time-frequency analysis.
For a propeller with N identical blades, N blades are equally spaced.The initial rotation angle of each blade is: According to the superposition principle, the echo is: Let the distance from any scattering point P t to O on the blade be l p .If the position is P 0 (0, y 0 , z 0 ) at t = 0 and the position is P t (0, y t , z t ) at t = t, then y t = −l p sin(ψ 0 + ω v t) and z t = l p cos(ψ 0 + ω v t), where ψ 0 is the initial rotation angle.In the sonar coordinate system, the coordinate of P t is (y t , z 0 + z t ), that is, P t [−l p sin(ψ 0 + ω v t), z 0 + l p cos(ψ 0 + ω v t)].The distance from P t to O is: Since z 0 >> l p , the distance is approximately: The reflected wave of scattering point P t is: The reflected wave is integrated into the blade length L, and the echo of the blade is: The micro-Doppler frequency is: The size of the micro-Doppler frequency is related to the blade length, rotation angular velocity, and wavelength.The echo amplitude flash and micro-Doppler shift are unique characteristics of underwater targets in the time and frequency domains.Micro-Doppler features are considered important features of underwater target recognition.The flash characteristics can be observed in the time domain, while the micro-Doppler frequency shift requires joint time-frequency analysis.
For a propeller with N identical blades, N blades are equally spaced.The initial rotation angle of each blade is: According to the superposition principle, the echo is: Sensors 2023, 23, 8066 8 of 15 The micro-Doppler frequency of the kth blade is:

Sparse Reconstruction Combined with Time-Frequency Distribution
The echo of a moving target is a nonstationary signal superimposed by many complex modulation components.

Sparse Representation of Echo Data
Supposing that an existing overcomplete dictionary D = {g r (t), r = 1, 2, ..., K}, where g k (t) represent atoms.In the Hilbert space H = R N where K > N, for any echo data x(t) ∈ H, if M atoms are selected adaptively in dictionary D, the nonlinear approximation of the echo is: where the set of atoms corresponding to the integer subscript set I M is B = {g r (t), r ∈ I M }.
Let Γ = {I m } be the set of all subscripts and c = {c r , r ∈ I M } be the nonlinear representation coefficients.If the potential of the subscript set is Card(I M ) < N, then c is the sparse representation of the signal in dictionary D. In particular, if Card(I M ) = min Card(I m ), c is called an optimal sparse representation.Equation ( 20) is written in matrix form x = Dc, where x = {x(t), t = 0, 1, ..., N − 1}.The solution of coefficient c can be expressed as a non-convex optimization model: where ||•|| 0 is the L0 norm.The sparse representation results of MP after M iterations are: In the formula, c m = R m x(t), g r m is the sparse decomposition coefficient, where • represents the inner product, g r m is the optimal atom chosen in each iteration, R M x(t) is the residual signal obtained after the Mth iteration, and R 0 x(t) = x(t).Ignoring the residual signal, sparse representation of the echo is solved by the Wigner Ville distribution, namely: The first term, on the right of the above equation, is the sum of the Wigner Ville distribution of the self-term between atoms, and the second term is the cross-term between atoms.The first term can correctly represent the time-frequency distribution structure of signal x(t), and can be obtained using the optimal atom g r m selected in each iteration and the decomposition coefficient c m .Therefore, a time-frequency distribution without interatomic cross-term interference can be obtained by retaining the first term of the above equation: Clearly, the properties of this time-frequency distribution are determined by the Wigner Ville distribution, and it has a good time-frequency aggregation and time-frequency resolution.

Design of Time-Frequency Atomic Dictionary
The dictionary design is an important aspect of sparse decomposition and the key to analyzing micro-Doppler effects.The design of time-frequency atoms should benefit both time-frequency representation of nonstationary signals and the concise representation of signals.
The micro-Doppler effect is a complex modulation form.The time-frequency characteristic of echoes is an oscillation curve in the time-frequency plane.Although Gabor, Chirplet, Dopplerlet, and other traditional analytical dictionaries have certain versatility, they cannot easily and succinctly represent such complex modulation signals.Based on the time-frequency characteristics of the signal, a kind of time-frequency atom can be designed to match the echo modulation form to represent the echo succinctly.To match the modulation form of the echo, ignoring the known carrier frequency, the time-frequency atoms can be constructed as follows: where w is a window function to fit the envelope of the signal, u is a displacement factor, and A is the oscillation amplitude factor.ω = 2π f r is the oscillation frequency factor.The amplitude and phase of an echo are all modulated.According to the modulation form of an echo, the constructed time-frequency atoms are: The time-frequency atoms represented in Equations ( 25) and ( 26) are called Mi-croDopplerlet atoms.The time-frequency distribution of an echo can be obtained by substituting the optimal atoms selected by MP into Equation (24).The designed timefrequency atoms match the echo modulation form.Therefore, the atoms can describe the complex modulation law and time-frequency distribution of the echo, which can represent the echo succinctly.When the echo is sparsely decomposed, the instantaneous frequency of the selected optimal time-frequency atom equals the micro-Doppler frequency of each echo component.Therefore, the micro-Doppler frequency is: where Φ(t) = arg{g r } and arg{•} means extracting the phase information.The value of C is determined by amplitude scintillation and amplitude flash.For a simple harmonic vibration point, C = 2, and for a simple harmonic vibration plane and rotation propeller, C = 1.
When the envelope and wavelength parameters of the emission signal are known, the time-frequency atoms have only three parameters: A, ω, and u.Discretizing these three parameters, the set obtained is the time-frequency atomic dictionary, the MicroDopplerlet dictionary.The discretizing method is given below.
Constrained by equal interval sampling, u is the equal interval value, that is, u = n/ f s , where f s is the signal sampling frequency.The micro-Doppler frequency f md is a cosine signal, so the sampling frequency should satisfy f s > 2 f md .A and ∆ have a linear relationship, namely A = m∆, where ∆ is a discrete step.The relationship between ω and f md is complex, affected not only by the linear relationship, but also by nonlinear sinusoidal relationships.For simplicity, the equal interval discretization approximation can be used.ω is represented by the corresponding digital angular frequency Ω, that is, ω = Ω f s .Because of f s > 2 f md , Ω = 2π f md / f s < π.Therefore, Ω can be discretized evenly in the range of [0, π).
If the amplitude A µ and frequency f µ of the micro-Doppler frequency are obtained by FFT, the vibration frequency of a vibration point and plane is ω v = 2π f µ .The vibration amplitude of a vibration point is A v = λA µ /2ω v .The vibration amplitude of a vibration plane is A v = λA µ /ω v .The rotation velocity of the blade is Ω = 2π f µ and the blade length is L = λA µ /Ω.

Computer Simulation
Simulation experiment 1: Micro-Doppler effect of a typical single component Research on the micro-Doppler effect of the reflected waves of simple harmonic vibration points and planes and a rotation propeller by sparse reconstruction combined with time-frequency distributions can reveal the changing law of time-frequency and amplitude characteristics and extract the micro-Doppler frequency of the reflected waves.Suppose the emission signal frequency is 30 kHz, the signal length is 300 ms, and the sampling frequency is 300 kHz.When the sound velocity is 1500 m/s, the signal wavelength is 5 cm.The amplitude of the vibration point is 2 cm, and the vibration frequency is 100 Hz.In the case of a radius of the vibration circular plane of 1 m, vibration amplitude of 4 cm, vibration frequency of 60 Hz, blade length of a 5-blade propeller of 1 m, and rotational velocity of 10 Hz, the reflected waves, time-frequency distribution, and micro-Doppler frequency of three micro-motion components are shown in Figures 5-7, respectively.The micro-Doppler frequency of a vibration point varies with a sine law.The micro-Doppler effect of a vibration plane and propeller is more complex.Besides the micro-Doppler frequency presenting sinusoidal variation, the reflected wave amplitude also appears as a flash effect.
Compared with the vibration point and plane, blades have a larger size, so the micro-Doppler effect is more obvious.The five micro-Doppler frequencies of a 5-blade propeller are shown in Figure 7c.Since each blade with equal length is installed at equal intervals, the amplitude and frequency of the five micro-Doppler frequencies are the same, and the phase difference between blades is 2π/5.
Simulation experiment 2: Sparse representation analysis of multi-component superimposed wave The vibration amplitude of the vibration point is set to 1 cm and the vibration frequency to 100 Hz.The vibration amplitude of the vibrating plane is set to 1 cm and the vibration frequency to 60 Hz, with a 5-blade propeller with a length of 1 m and a speed of 10 Hz.The experimental results are simulated as follows.
Let the target echo be the superposition of the reflected waves of a simple harmonic vibration point and plane and rotation blades.The echo is sparsely decomposed into the designed MicroDopplerlet time-frequency atomic dictionary by MP, and the timefrequency analysis is performed on the reflected wave of separated single components to extract the micro-Doppler frequency.The parameters are extracted, including the vibration amplitude and frequency of the point and plane, the length of the blade, and the rotational velocity.
In simulation experiment 2, the micro-motion parameters of the three components' emission signal parameters are the same as those in simulation 1.The signal-to-noise ratio is 15 dB.The sparse decomposition grids are divided as follows: amplitude interval of micro-motion is 1 mm, frequency interval of micro-motion is 1 Hz, and time delay is 1 ms.The number of MP iterations is 10.The simulation results are shown in Figures 8-11.
quency to 100 Hz.The vibration amplitude of the vibrating plane is set to 1 cm and the vibration frequency to 60 Hz, with a 5-blade propeller with a length of 1 m and a speed of 10 Hz.The experimental results are simulated as follows.
Let the target echo be the superposition of the reflected waves of a simple harmonic vibration point and plane and rotation blades.The echo is sparsely decomposed into the designed MicroDopplerlet time-frequency atomic dictionary by MP, and the time-frequency analysis is performed on the reflected wave of separated single components to extract the micro-Doppler frequency.The parameters are extracted, including the vibration amplitude and frequency of the point and plane, the length of the blade, and the rotational velocity.
In simulation experiment 2, the micro-motion parameters of the three components' emission signal parameters are the same as those in simulation 1.The signal-to-noise ratio is 15 dB.The sparse decomposition grids are divided as follows: amplitude interval of micro-motion is 1 mm, frequency interval of micro-motion is 1 Hz, and time delay is 1 ms.The number of MP iterations is 10.The simulation results are shown in Figures 8-11             Figure 8 shows the MicroDopplerlet time-frequency atoms of the vibration point, plane, and blades.These time-frequency atoms match the respective modulation forms of the reflected waves.Figure 9 shows that the designed time-frequency atoms and MP are used to decompose and separate the echo of the three components correctly.Figure 10a shows the time-frequency distribution.In Figure 10b, the five sine waves with a large amplitude and low frequency are the micro-Doppler frequencies of the five blades.The micro-Doppler frequency amplitude can reach 1.25 kHz.The two sinusoidal waves with a smaller amplitude and higher frequency are the micro-Doppler frequencies of the point and plane, respectively.The micro-Doppler frequency amplitude of the vibration point is about 510 Hz, and that of the vibration plane is about 300 Hz.

Conclusions
Based on the analysis of the micro-motions of an underwater target, the echo models of a simple harmonic vibration point and plane and rotation propeller were established in this paper, and the complex modulation characteristics of micro-Doppler effects are revealed.A method of sparse reconstruction combined with time-frequency distributions was presented, and a MicroDopplerlet time-frequency atomic dictionary matching with echo modulation was designed to realize the sparse decomposition of echoes and the analysis of micro-Doppler effects.The conclusions are as follows: (1) The vibration amplitude of the underwater target shell is at the centimeter level, and the vibration frequency ranges from tens of hertz to hundreds of hertz.The micro-Doppler frequency shift can reach tens of hertz to hundreds of hertz.The propeller has a large blade length and rotation rate, and its micro-Doppler frequency shift can reach hundreds to thousands of hertz; (2) The computer simulation shows that using the designed MicroDopplerlet timefrequency atomic dictionary not only separates vibration points, vibration surfaces, and propellers but can also distinguish the superimposed modulation components of the same type of components.In addition, the micro-Doppler frequency and vibration frequency, vibration amplitude, and the propeller rotation speed and length parameters of each micro-motion component are extracted.The extracted parameters are consistent with the given results, proving the correctness of the designed time-frequency atomic dictionary and the established model.

Figure 1 .Figure 1 .
Figure 1.Simplified model of submarine.The propeller and shafting-shell coupled system have an important influence on the shell vibration.The unsteady external excitation force of the propeller during operation and the friction force of the shafting and bearing caused by the propeller rotation are important excitation sources for the vibration of the shell structure.At present, scholars have paid attention to the exciting force of the propeller on the shell [29-33].Hassan et al. con-Figure 1. Simplified model of submarine.

Figure 2 .
Figure 2. Geometric relationship between vibration point and sonar.
v A x y ′ ′ .The echo model can be established by analyzing the displacement law of simple harmonic points.

Figure 2 .
Figure 2. Geometric relationship between vibration point and sonar.

Figure 3 .
Figure 3. Geometric relationship between vibration plane and sonar.

Figure 3 .
Figure 3. Geometric relationship between vibration plane and sonar.

Figure 4 .
Figure 4. Geometric relationship between blades and sonar.
It is difficult to use joint time-frequency distribution to analyze micro-Doppler effects and local joint time-frequency information.Sparse decomposition can better separate multi-component signals and effectively reveal the time-varying characteristics of nonstationary signals.It provides a method for analyzing micro-Doppler effects and for feature extraction of underwater targets.This paper presents a method of sparse reconstruction combined with time-frequency distributions.A time-frequency atomic dictionary is designed according to the modulation form of micro-motion components' echo.Multi-component signals are separated into singlecomponent signals by sparse decomposition of the echo.The Wigner Ville distribution is used to analyze the separated single-component signals.This method constructs a time-frequency distribution without interatomic cross-term interference and with good time-frequency aggregation.It can analyze the micro-Doppler effects of echoes well.
Figure8shows the MicroDopplerlet time-frequency atoms of the vibration point, plane, and blades.These time-frequency atoms match the respective modulation forms of the reflected waves.Figure9shows that the designed time-frequency atoms and MP are used to decompose and separate the echo of the three components correctly.Figure10ashows the time-frequency distribution.In Figure10b, the five sine waves with a large amplitude and low frequency are the micro-Doppler frequencies of the five blades.The micro-Doppler frequency amplitude can reach 1.25 kHz.The two sinusoidal waves with a smaller amplitude and higher frequency are the micro-Doppler frequencies of the point and plane, respectively.The micro-Doppler frequency amplitude of the vibration point is about 510 Hz, and that of the vibration plane is about 300 Hz.Figure11compares the true and estimated values of the micro-motion parameters of the point, plane, and blades.(a), (b), and (c) are the true values, and (d), (e), and (f) are the estimated values; the extracted parameters are consistent with the true values.The simulation results show that the micro-Doppler frequency processed by FFT can correctly estimate the parameters, such as the velocity and length of blades, the vibration amplitude, and the frequency of the vibration point and plane.

Table 1 .
Estimation of micro-motion amplitude, micro-motion frequency, and micro-Doppler frequency of propellers and different shells.