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

Authors: Natee Thong-un; Minoru K. Kurosawa
Journal: Advances in Acoustics and Vibration
 (2014)

---

## Abstract

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.


---

## Body

## 1. 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 cross-correlation 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 idea is based on the expectation-maximization (EM) algorithm. The EM algorithm [7], although iterative in nature, is guaranteed under certain mild conditions to converge and, at convergence, to produce at least a local maximum. The EM algorithm has the desirable property of increasing the likelihood of finding a local maximum at each step. Moreover, the EM algorithm has been used to estimate the unknown parameters in the overlapping Gaussian signals [8]. Note that, for the overlapping LPM ultrasonic signals, the EM algorithm is interesting in that this method can estimate the unknown parameters of each signal. Simulation results indicate that the EM algorithm can provide the Doppler velocity estimation of the overlapping LPM ultrasonic signals.


## 2. Doppler Velocity Estimation for a Single LPM Echo

The discrete signal of the LPM ultrasonic signal can be modeled as [3](1)fn=Asin2πl0pbln⁡pbp1l0n+1
for   n=0,1,2,…,N-1 where A is the rectangular envelope pulse of the amplitude in a wide band frequency, pb is the sweep period band of the reference LPM ultrasonic signal, p1 is the starting period time, l0 is the length of the LPM ultrasonic signal, and N 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
(2)xn=sn+ηn,sn=Bsin2πl0pbln⁡αvdpbp1l0n-nTOF+1,
where nTOF is a sample when TOF happens; α(vd)=(v+vd)/v is the Doppler-shift function of the modulated LPM ultrasonic signal expressed in [4], where v is the ultrasonic velocity in air; and vd is the object velocity. α(vd) 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 B=A10-2.11vnTOF [1]. To verify the ability to estimate the Doppler velocity for a noisy situation, noise can be added to the model. The noise η(n) is represented as the additive white Gaussian noise (WGN). In this case, to estimate the Doppler velocity (vd), because no overlapping signals are present, it is not necessary to use the E-step of the EM algorithm. The estimation can be only performed in the maximization, M-step. In the statistical model, we assume that the x(n) of each sample is independent and has an identical distribution. Then, their probability density function (pdf) is the normal distribution illustrated as
(3)px;vd=∏n=0N-112πσ2exp⁡-12σ2xn-sn;vd2,
where σ2 is the variance of the Gaussian distribution and xn is the sampling data[x0,x1,…,xN-1]T. It can be expressed in a form of the noisy signal as
(4)ηn;0<n≤nTOF-1xn=sn-nTOF+ηn;nTOF≤n≤nTOF+Ld-1ηn;nTOF+Ld≤n≤N-1.
The pdf of (3) can be expanded to
(5)px;vd=∏n=0nTOF-112πσ2exp⁡-12σ2x2n·∏n=nTOFnTOF+Ld-112πσ2hhhhhhhhh×exp⁡-12σ2xn-sn-nTOF2·∏n=nTOF+LdN-112πσ2exp⁡-12σ2x2n.
The maximum likelihood estimation (MLE) function is simplified to
(6)p(x;vd)=12πσ2N/2exp⁡-12σ2∑n=0N-1x2nhhhhhhhh·∏n=nTOFnTOF+Ld-1exp⁡-12σ2-2xnsn-nTOF+s2n-nTOFhhhhhhhhhhhhhhhhhhhhhhhhhhh+s2n-nTOF-12σ2.
The MLE is found by only maximizing exp⁡⁡[-(1/2σ2)∑n=nTOFnTOF+Ld-1(-2x(n)s(n-nTOF) +s2(n-nTOF))] or, equivalently, by minimizing
(7)∑n=nTOFnTOF+Ld-1-2x(n)s(n-nTOF)+s2(n-nTOF).
Then Ld=(v/v+vd)L0 [4], where Ld and L0 are the lengths of Doppler shift and reference signals, respectively. When carefully scrutinizing in (7), we note that ∑n=nTOFnTOF+Ld-1s2(n-nTOF) does not deal with the estimation between the measured data x(n) and the model of the LPM ultrasonic signal s(n;vd); therefore, we can ignore it. Therefore, the MLE of nTOF is expressed by maximizing
(8)∑n=nTOFnTOF+Ld-1xnsn-nTOF.
Simply speaking, (8) is the correlation between the measured data x(n) and the model signal s(n). 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
(9)v^d=argmax⁡vd∈Θ∑n=nTOFnTOF+Ld-1x(n)s(n-nTOF),
where Θ is the region of interest and vd 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.Estimation searching and probability distribution of the Doppler velocity.
(a)
(b)
(c)
(d)


## 3. Doppler Velocity Estimation for Overlapping LPM Echoes

In the previous section, the Doppler velocity estimation of a LMP ultrasonic echo can be performed. However, 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(10)yn=∑m=1Msmn+ηn,
where sm(n) is the overlapping LPM ultrasonic echo model and y(n) 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 y(n) are independent and exhibit identical distributions, their probability density function (pdf) is the normal distribution according to
(11)py;vd1=∏n=0N-112πσ2exp⁡-12σ2yn-s1n2
or equivalently maximizing(12)v^d1=argmax⁡vd∈Θ∑n=nTOF1nTOF1+Ld1-1yns1n-nTOF1.
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(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 into the individual echoes by means of the expectation and maximization steps. First, we assume that ym is characterized as the unknown data for at mth LPM ultrasonic echo under the WGN process ηm(n) as
(13)ymn=smn+ηmn.
The unknown data of each individual echo can be formed as a linear transformation of the measured data y=∑m=1Mym. The E-step of the EM algorithm is able to compute ym in terms of the measured and estimated data as
(14)y^mkn=smkn+βmyn-∑i=1Msik,
where ∑m=1Mβm=1; for the example of simulation, if there are the two overlapping echoes, β1 and β2 are 0.5. k represents the number of iterations, which simply keeps a balance of the observation and the estimation of every echo. The expected datay^m is adjusted according to the difference between the observation and the estimation. Then, the M-step of the EM algorithm, which involves the maximization of the pdf, is expressed as
(15)v^dmk+1=argmax⁡vdm∑n=0N-1y^m·smn-nTOFm.
The definition can be concluded in the form of an algorithm for the Doppler velocity estimation of the M-number of the overlapping LPM ultrasonic echoes in Algorithm 1.(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)
(b)
(c)Algorithm 1.
(1) Use the Maximum likelihood estimation (M-step) to estimate the Doppler velocity of the first echo and set the starting iteration k=0.
(2) Make the initial guesses of the residual echoesvd20,vd30,…,vdM0.
(3) Form=2,3,…,M, compute the expected echoes (E-step):
(16)y^mkn=smkn+βmyn-∑i=1Msik.
(4) Form=2,3,…,M, compute the maximization (M-step):
(17)v^dmk+1=argmax⁡vdm∑n=0N-1y^m·smn-nTOFm.
(5) Adjustvd2k,vd3k,…,vdMk for the update:
(18)sm=2,3,…,Mk+1n=Bsin2πl0pbln⁡pbp1l0n-nTOFm+1αvdm=2,3,…,Mkhhhhhhhhhhhhhhhhhhhhhhhhh×pbp1l0n-nTOFm+1.
(6) Check the convergence criterion: if∥v^dmk+1-v^dmk∥≤tolerance, then stop.
(7) Setk→k+1 and back to step 3.
Note that, in this algorithm, the grid search is a main approach for each echo. Iterations for convergence of two and three overlapping LPM ultrasonic echoes are shown Figure3.Iterations for overlapping LPM ultrasonic echoes with a 0.05 m/s step of the grid-search method. (a) Two overlapping LPM ultrasonic echoes withvd1 and vd2 of 3.3 m/s and 2 m/s, respectively. (b) Three overlapping LPM ultrasonic echoes with vd1, vd2, and vd3 of 1 m/s, −2 m/s, and 5 m/s, respectively.
(a)
(b)


## 4. 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 Figures1(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.Table 1
Estimation results for the two overlapping LPM ultrasonic echoes with 10 dB SNR.
LPM ultrasonic signals
First echo
Second echo
Actual velocities (m/s)
−5
3.5
3.6
−1.6
3.3
2
Estimated velocities (m/s)
−4.95
3.55
3.65
−1.65
3.3
1.95Table 2
Estimation results for the three overlapping LPM ultrasonic echoes with 22 dB SNR.
LPM ultrasonic signals
First echo
Second echo
Third echo
Actual velocities (m/s)
4
3.5
2
1
−2
5
−4.5
3
−5
Estimated velocities (m/s)
4.05
3.6
2.05
1
−2.1
5.1
−4.45
3.1
−5.15(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.
(a)
(b)
(c)
(d)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.


## 5. 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 theM-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.



---

*Source: 921876-2014-11-12.xml*