Doppler Velocity Estimation of Overlapping Linear-Period-Modulated Ultrasonic Waves Based on an Expectation-Maximization Algorithm

The occurrence of an overlapping signal is a significant problem in performing multiple objects localization. Doppler velocity is sensitive to the echo shape and is also able to be connected to the physical properties of moving objects, especially for a pulse compression ultrasonic signal. The expectation-maximization (EM) algorithm has the ability to achieve signal separation. Thus, applying the EM algorithm to the overlapping pulse compression signals is of interest. This paper describes a proposed method, based on the EM algorithm, of Doppler velocity estimation for overlapping linear-period-modulated (LPM) ultrasonic signals. Simulations are used to validate the proposed method.


Introduction
Recently, acoustic systems have been used in many industrial and educational applications because such systems have the advantages of a variety of acoustic transducers, low cost, small size, and simple hardware.For example, ultrasonic distance measurement is a modern and relatively flexible approach for environment recognition.The ultrasonic distance measurement determination uses the time-of-flight (TOF) approach, which is the time for sound waves to travel between the sound source and the objects.In general, TOF can simply be computed by performing a cross-correlation between the transmitted and received signals [1].As a result, TOF can be measured by the maximum peak during a storage time.
A linear-frequency-modulated (LFM) signal, which has a linear variation in an instantaneous frequency with time, is an example of a pulse compression signal used for ultrasonic distance measurement [1].However, for the case of the distance measurement of the moving object, the LFM signal is unsuitable because of the shifted frequency of the Doppler effect.The problem is that the cross-correlation cannot completely be achieved between the transmitted and received signals.Therefore, a linear-period-modulated (LPM) signal, the period of which is linearly swept with time, has been demonstrated to solve the Doppler-shift problem of the LFM signal [2,3].A LPM signal can be used in the crosscorrelation by providing the TOF with the Doppler-shift compensation via Doppler velocity estimation [4].
However, the Doppler-shift problem is more challenging for the multiple-echoes case.The overlapping signals due to many objects will exhibit many Doppler velocities.If each Doppler velocity can be completely separated, then the Doppler velocity can be estimated for each individual echo.In the case of a pulse compression signal, the LFM echoes can be discriminated by the use of the Fractional Fourier Transform [5,6].The Fractional Fourier Transform method can identify the time delay and the length of each echo.The length of the echo can be used to determine the Doppler-shifted frequency of the reflectors, which is directly related to the Doppler velocity.However, the Fractional Fourier Transform cannot be adapted for use with a nonlinear frequency modulated signal, such as the LPM signal.Therefore, this paper describes the proposed method of estimating the Doppler velocities embedded in the overlapping LPM signals.This proposed

Doppler Velocity Estimation for a Single LPM Echo
The discrete signal of the LPM ultrasonic signal can be modeled as [3] for  = 0, 1, 2, . . .,  − 1 where  is the rectangular envelope pulse of the amplitude in a wide band frequency,   is the sweep period band of the reference LPM ultrasonic signal,  1 is the starting period time,  0 is the length of the LPM ultrasonic signal, and  is the total number of samples.When we consider the LPM ultrasonic echo produced by reflection from the moving object, it can be expressed in a form of the noisy signal as where  TOF is a sample when TOF happens; (V  ) = (V + V  )/V is the Doppler-shift function of the modulated LPM ultrasonic signal expressed in [4], where V is the ultrasonic velocity in air; and V  is the object velocity.(V  ) is directly related to the length of the received LPM ultrasonic echo.The amplitude of the received LPM ultrasonic echo accounting for attenuation along the distance travelled is  = 10 −2.11VTOF [1].To verify the ability to estimate the Doppler velocity for a noisy situation, noise can be added to the model.The noise () is represented as the additive white Gaussian noise (WGN).In this case, to estimate the Doppler velocity (V  ), because no overlapping signals are present, it is not necessary to use the -step of the EM algorithm.The estimation can be only performed in the maximization, -step.In the statistical model, we assume that the () of each sample is independent and has an identical distribution.Then, their probability density function (pdf) is the normal distribution illustrated as where  2 is the variance of the Gaussian distribution and () is the sampling data[(0), (1), . . ., ( − 1)]  .It can be expressed in a form of the noisy signal as The pdf of (3) can be expanded to The maximum likelihood estimation (MLE) function is simplified to The MLE is found by only maximizing exp or, equivalently, by minimizing Then   = (V/V+V  ) 0 [4], where   and  0 are the lengths of Doppler shift and reference signals, respectively.When carefully scrutinizing in (7), we note that ∑  TOF +  −1 = TOF  2 ( −  TOF ) does not deal with the estimation between the measured data () and the model of the LPM ultrasonic signal (; V  ); therefore, we can ignore it.Therefore, the MLE of  TOF is expressed by maximizing Advances in Acoustics and Vibration Simply speaking, (8) is the correlation between the measured data () and the model signal ().Therefore, (8) provides the maximum likelihood function with the exact Doppler velocity under the peak value of the cross-correlation.It can be defined as where Θ is the region of interest and V  is the Doppler velocity estimator.This idea was evaluated by performing a MATLAB computer simulation.The period of the transmitted LPM ultrasonic signal was swept from 20 s to 50 s.The sampling frequency rate was 12.5 MHz.A pulse length was equal to 3.247 ms.In the environment setup, the propagation velocity of an ultrasonic wave in air was 343.6426 m/s at temperature 20 ∘ C, and the attenuation was set at 2.11 dB/m along the travelling distance.In the simulation, the Doppler velocity of the interesting region was varied from −5 m/s to 5 m/s, and the instantaneous position of the moving object far away from the speaker was 1 m.The grid-search method was used for determining the Doppler velocity via the maximum likelihood function.The search had a step resolution equivalent to 0.01 m/s. Figure 1 shows the maximum Doppler velocity estimation, which was computed from (9), when the actual values of which were assumed that an object was moving in and out to a sound source with 4.8 and −2.6 m/s, respectively.For each Doppler velocity, 100 trials were performed to determine the probability distribution.The SNR of the setup with randomness was 0 dB.

Doppler Velocity Estimation for Overlapping LPM Echoes
In the previous section, the Doppler velocity estimation of a LMP ultrasonic echo can be performed.However, Advances in Acoustics and Vibration the estimation will be more complicated if M-number of overlapping echoes are present.To proceed further, considering the overlapping signals, the LPM discrete signal is where   () is the overlapping LPM ultrasonic echo model and () is the measured echo.In general, the amplitudes under the overlapping echoes are directly related to their distance; therefore, the other echoes except for the first echo can be treated as noise.When the observations () are independent and exhibit identical distributions, their probability density function (pdf) is the normal distribution according to or equivalently maximizing Equations ( 11) and ( 12) represent the correlation of the measured data with the first echo.Assume that two overlapping signals are at distances of 0.7 m and 0.75 m from the sound source, which is pictured in Figure 2  into the individual echoes by means of the expectation and maximization steps.First, we assume that   is characterized as the unknown data for at mth LPM ultrasonic echo under the WGN process   () as The unknown data of each individual echo can be formed as a linear transformation of the measured data  = ∑  =1   .The -step of the EM algorithm is able to compute   in terms of the measured and estimated data as where ∑  =1   = 1; for the example of simulation, if there are the two overlapping echoes,  1 and  2 are 0.5. represents the number of iterations, which simply keeps a balance of the observation and the estimation of every echo.The expected data ŷ is adjusted according to the difference between the observation and the estimation.Then, the -step of the EM algorithm, which involves the maximization of the pdf, is expressed as The definition can be concluded in the form of an algorithm for the Doppler velocity estimation of the -number of the overlapping LPM ultrasonic echoes in Algorithm 1.

Algorithm 1.
(1) Use the Maximum likelihood estimation (step) to estimate the Doppler velocity of the first echo and set the starting iteration  = 0.

Simulation Results and Discussion
The algorithm was tested by performing simulations of the LPM ultrasonic echoes.In the case of one Doppler velocity estimation, the Doppler velocity and distance to the moving object were determined under 100 times according to support statistics.The probability distribution of the estimated Doppler velocity is illustrated in Figures 1(c) and 1(d).Then, we considered examples of two overlapping LPM ultrasonic echoes, the distances of which, with a sound source, were 0.65 m and 0.7 m, and three overlapping echoes, the distances of which were 0.5 m, 0.65 m, and 0.75 m, as shown in Figure 4, under the additive WGN signal.The Doppler velocities in each object were assumed that the maximum velocity was 5 m/s and the minimum velocity was −5 m/s.In addition, the maximum distance of the multiple objects is 1 m.The Doppler velocity estimation results of the overlapping LPM ultrasonic echoes are tabulated in Table 1.
The SNR under consideration was 10 dB.Table 2 describes the estimation results of the three closed overlapping LPM ultrasonic echoes.In the table, the actual parameters in the simulation are listed.The SNR used to support the Doppler velocity estimation of three overlapping echoes was 22 dB.The step size in both cases for the grid search was 0.05 m/s.
Ultrasonic velocity measurement by high-resolution Doppler velocity estimation was evaluated with computer simulations using MATLAB.In simulation, the environment was assumed that SNR of the reflected echo was approximately 20 dB by adding normally distributed random noise.As for the Doppler velocity estimation referred in [4], which requires two cycles of LPM signal, estimation is based on time duration of two peaks in the cross-correlation function.It does not have any problem with only one moving object but when the number of moving objects is more than one object, it is difficult to estimate the Doppler velocities due to many peaks.Hence, the proposed method is designed to estimate the Doppler velocity of two and three objects using only one cycle of LPM signal.However, this is based on EM algorithm and requires convergence to produce at least a local maximum.It is under the iterative procedure in nature.

Conclusion
In this paper, we first presented the Doppler velocity estimation of the overlapping LPM ultrasonic echoes using the EM algorithm.The single echo results were achieved by the M-step of algorithm, that is, the maximum likelihood estimation.The cases of three and two overlapping LPM ultrasonic echoes were addressed by first performing the Doppler velocity estimation of the first echo, followed by the estimation of the Doppler velocities of the residual echoes using the EM algorithm.The grid search was used for the searching engine.Our proposed method was evaluated by performing MATLAB computer simulations.The simulation results confirmed the validity of the proposed method.

Figure 1 :
Figure 1: Estimation searching and probability distribution of the Doppler velocity.

Figure 2 :
Figure 2: (a) Two overlapping LPM ultrasonic echoes at distances of 0.7 m and 0.75 m from the sound source.(b) Maximizing function of two overlapping echoes: the first and second echoes with 3.3 m/s and 2 m/s, respectively.(c) Percentage of overlap.
(a).The first echo has a Doppler velocity of 3.3 m/s, and the second echo has a Doppler velocity of 2 m/s.The maximum likelihood function of the Doppler velocity of the first echo is shown in Figure 2(b).In addition, Figure 2(c) shows the ability of the model to estimate the Doppler velocity of the first echo, in which the second echo is embedded, when the amount of overlap is changed.With the estimate of the Doppler velocity of the first echo, the Doppler velocity of the residual echo now must be estimated.To accomplish this estimation, the EM algorithm allows us to separate the overlapping echoes

Figure 3 :
Figure 3: Iterations for overlapping LPM ultrasonic echoes with a 0.05 m/s step of the grid-search method.(a) Two overlapping LPM ultrasonic echoes with V 1 and V 2 of 3.3 m/s and 2 m/s, respectively.(b) Three overlapping LPM ultrasonic echoes with V 1 , V 2 , and V 3 of 1 m/s, −2 m/s, and 5 m/s, respectively.

Figure 4 :
Figure 4: (a) Two overlapping LPM ultrasonic echoes at distances of 0.65 m and 0.7 m from a sound source.(b) Estimation of two overlapping LPM ultrasonic echoes.(c) Three overlapping LPM ultrasonic echoes at distances of 0.5 m, 0.65 m, and 0.75 m from a sound source.(d) Estimation of the three overlapping LPM ultrasonic echoes.

Table 1 :
Estimation results for the two overlapping LPM ultrasonic echoes with 10 dB SNR.

Table 2 :
Estimation results for the three overlapping LPM ultrasonic echoes with 22 dB SNR.