**2. System Model**

A distance measurement and positioning system using UWB technology is illustrated in Figure 1, where *di* and *ε<sup>i</sup>* are the thickness and the relative permittivity of the *i th* layer in the nondestructive environment, respectively. *s*(*t*) denotes the transmitted pulse signals, and *r*(*t*) is the received signal.

For the environment with three layers and one buried object, the *r*(*t*) is represented by

$$\begin{aligned} r(t) &= \sum r\_i(t) + r\_{ob}(t) + n(t), \\ r\_1(t) &= A\_1 s(t - \tau\_1), \\ r\_2(t) &= A\_1 A\_2 s(t - \tau\_1 - \tau\_2), \\ r\_3(t) &= A\_1 A\_2 A\_3 s(t - \tau\_1 - \tau\_2 - \tau\_3), \\ r\_{ob}(t) &= A\_{ob} s(t - \tau\_{ob}), \end{aligned} \tag{1}$$

where *Ai* denotes the amplitude factor which represents the reflection and transmission properties of the propagation environment, and its value depends on the properties of the *i*th layers, *τ<sup>i</sup>* is the traveling time of UWB pulse in the *i*th layer. *Aob* and *τob* are respectively the amplitude and traveling time of the signal which is reflected from the buried object. *n*(*t*) represents the additive white Gaussian noise.

At the receiver side, the procedure for determining the position of buried object is indicated in Figure 2. Accordingly, the position of buried object is determined by the distance from it to the transceiver. The LMFA method is applied to calculate the relative permittivity of environment and these distances.

**Figure 2.** The procedure for determining the position of buried object.

In the UWB system, *s*(*t*) is generated based on Gaussian [25] or Hermite [26] functions and their derivatives. Compared to Gaussian pulses, the Hermite pulse is useful for parallel data transmission with high data rates, but it is hard to achieve in the real world [27]. Therefore, the Gaussian pulses UWB signals are applied in this work. A typical Gaussian pulse usually takes the form [25]:

$$\log(t) = A\_p e^{-2\pi \left(\frac{t}{\mu\_p}\right)^2},\tag{2}$$

where *Ap* denotes the amplitude of pulse, *μ<sup>p</sup>* is a factor which influences the amplitude and the width of Gaussian pulse, also called time normalization factor, the width of a pulse becomes narrower when the *μ<sup>p</sup>* is reduced. The *nth* derivative of Gaussian pulse, named *nth*-order monocycle, is

$$\log\_n(t) = A\_{np} \frac{d^n}{dt^n} e^{-2\pi \left(\frac{t}{\mu\_p}\right)^2} \,, \tag{3}$$

where *Anp* is the amplitude of the *nth*-order monocycle. The shapes of different types of Gaussian monocycles are indicated in Figure 3.

**Figure 3.** The different types of Gaussian monocycles.

The selected Gaussian pulse shape must meet FCC requirements about power and frequency range used. According to Figure 4, it is observed that for the power spectral density (PSD), only the first-order Gaussian monocycle does not comply with the FCC requirements about the effective isotropic radiated power (EIRP) level. In this paper, we restrict our performance analysis of UWB system to the fourth-order Gaussian monocycle in (4). However, our method can be applied for any Gaussian pulse shape.

$$\log\_4(t) = A\_P \left[ -12\pi + 96\pi^2 \left( \frac{t}{\mu\_p} \right)^2 - 64\pi^3 \left( \frac{t}{\mu\_p} \right)^4 \right] e^{-2\pi \cdot \left( \frac{t}{\mu\_p} \right)^2}. \tag{4}$$

**Figure 4.** The power spectral density (PSD) of the different derivatives of Gaussian pulses.

It is observed that with impulse radio UWB (IR-UWB) signal, the correlation and therefore the propagation time can be calculated. However, the IR-UWB signal is not strong enough to be processed at the receiver in the case of multiple layers medium (the value of cross-correlation is lower than the noise lever). Hence, an orthogonal pseudo sequence *a*(*n*) is applied to the UWB signal in pulse position modulation (PPM) to ensure the cross-correlation is obtained when receiving many reflected signals from many surfaces and buried objects. The conventional UWB-PPM signal with the fourth-order Gaussian monocycle is given as:

$$s(t) = \sum\_{n=1}^{N} \mathbf{g}\_4(t - nT\_r - a\_n T\_{\text{FPM}}),\tag{5}$$

where *Tr* is pulse repetition period and *T*PPM is the time shift associated with binary PPM, *an* ∈ {0, 1}, and *N* is the length of UWB pulses. In this modulation technique, when *an* is 0, there is no additional time shift that modulates the pulse *g*4(*t*) and a time shift *T*PPM is added to *g*4(*t*) when *an* is 1. In our system model, *an* does not carry information, it only indicates whether the corresponding pulse is shifted or not.

The propagation distance *d* in the UWB-PPM system used for positioning the buried object in nondestructive medium is computed by:

$$d = \frac{c\tau}{2\sqrt{\varepsilon}},\tag{6}$$

where *c* is the velocity of light, and *τ*, *ε* are the traveling time of UWB pulses and the relative permittivity of the nondestructive medium, respectively.

The traveling time *τ* is calculated by the correlation between the received signal from the *i th* layer, *ri*(*t*), and the template waveform *ω*(*t*):

$$R\_i(\tau) = \int\_{-\infty}^{\infty} r\_i(t)\omega(t)dt. \tag{7}$$

With the conventional UWB-PPM systems, the template waveform at the receiver is (*g*4(*t*) − *g*4(*t* − *T*PPM)) and the correlation of this system is denoted by *R*0*<sup>i</sup>* and has the form:

$$R\_{0i}(\tau) = \int\_{-\infty}^{\infty} r\_i(t) [\mathcal{g}\_4(t) - \mathcal{g}\_4(t - T\_{\rm{PPM}})] dt. \tag{8}$$

Define the autocorrelation function of the fourth-order Gaussian monocycle in (4) as *RG*4(*τ*), we have

$$R\_{\rm G4}(\tau) = \int\_{-\infty}^{\infty} g\_4(t) g\_4(t - \tau) dt,\tag{9}$$

$$R\_{0i}(\tau) = R\_{G4}(\tau) - R\_{G4}(\tau - T\_{\rm PPV}).\tag{10}$$

The traveling time of UWB signal is the value that makes the correlation function hit its maximum; it corresponds to twice the transmission time from the transceiver to the reflective surface of the *i th* layer as shown in Figure 5 and given as

$$
\pi = \arg\, \pi(\max(\mathbb{R}\_i)),
\tag{11}
$$

where *max*(*Ri*) is the local extreme value corresponding to the reflected signals from the surface of the layers or object.

**Figure 5.** The local extreme values of correlation output corresponding to propagation distances in Figure 1.

According to (6), *τ* is a function of *ε* and the distance *d*, due to both *d* and *ε* are unknown; they should be estimated by large enough values of *τ* using LMFA. To increase the accuracy in the *τ* determination procedure, and in estimation of *d*, *ε*, the proposed UWB-PPM-ATS technique is applied.

#### **3. Proposal of Positioning Approach in Nondestructive Environments**

#### *3.1. Proposal of UWB-PPM-ATS*

From (1), (5), and (6), it can be seen that the time shift *T*PPM directly affects the quality of UWB-PPM systems in positioning applications. The performance of the system can be improved by selecting the appropriate value of the time shift. To increase the accuracy in the estimated distance and the position of buried objects by using the UWB-PPM system, the time shift in this system should be selected so that the ability to detect the received UWB pulse is the best. To accomplish this task, in UWB-PPM systems with the time shift *T*PPM is invariant, we recommend adjusting this time shift with a certain time constant to achieve its optimal value. The optimal value of *T*PPM defined in this paper is the value at which the UWB-PPM system gives the smallest error in distance estimation. So we propose an UWB pulse position modulation with an additional time shift (UWB-PPM-ATS). In the proposed technique, the pulse position will be changed with a time constant denoted by *ζ*. The signal of UWB-PPM-ATS is given as:

$$s(t) = \sum\_{n=1}^{N} g\_4(t - nT\_r - a\_n.(T\_{PPM} + \zeta)).\tag{12}$$

The "+" sign in (12) means that algebraic additions, so *ζ*, can take either positive or negative values. The effect of *ζ* on the quality of the system will be evaluated later.

An example of the conventional UWB, UWB-PPM, and UWB-PPM-ATS signal shapes are illustrated in Figure 6. In Figure 6, there is an example of a seven pulses sequence carrying seven bits 1010101. The pulse width is 0.28 ns, and it adopts the fourth derivative of the Gaussian function described by the blue line; the conventional PPM modulated pulses sequence is denoted by the red line. Whereby, when the bit value is 0, the transmitter will send out a pulse *g*40(*t*) as the original pulse (without any change from the unmodulated pulse) when the value bit is 1, the transmitter will send out a pulse *g*41(*t*) which is the original pulse with a time shift of *T*PPM. With the proposed UWB-PPM-ATS (black line), instead of the normal *T*PPM time shift, the new shift level is set to (*T*PPM + *ζ*).

**Figure 6.** The ultra-wide band (UWB) pulse shapes with pulse position modulation (PPM) and PPM-ATS.

The procedure for estimating distances in UWB-PPM-ATS system is the same as in the conventional system (UWB-PPM) and in (7); the template waveform is (*g*4(*t*) − *g*4(*t* − *T*PPM − *ζ*)). The correlation of the proposed system denoted by *R*1*<sup>i</sup>* and determined as:

$$R\_{1i}(\tau) = \int\_{-\infty}^{\infty} r\_i(t) [\mathcal{g}\_4(t) - \mathcal{g}\_4(t - T\_{\rm{PPM}} - \mathcal{\zeta})] dt. \tag{13}$$

So, we have:

$$R\_{1i}(\tau) = R\_{\rm G4}(\tau) - R\_{\rm G4}(\tau - T\_{\rm PPM} - \zeta) \big] dt. \tag{14}$$

The shapes of *RG*4(*τ*),*R*0*i*(*τ*) and *R*1*i*(*τ*) are shown in Figure 7.

Figure 7 indicates the autocorrelation function of the fourth-order Gaussian monocycle *RG*<sup>4</sup> and correlation functions of UWB-PPM *R*0*<sup>i</sup>* and UWB-PPM-ATS *R*1*<sup>i</sup>* systems (with additional time *ζ* = −0.08 ns and *ζ* = 0.08 ns); those functions are compared at the value of *τ* = 0. In the figure below inside Figure 7, the blue line with "+" sign denotes the value of *R*0*i*(*τ*), the red line denotes the value of *R*1*i*(*τ*) with *ζ* gets a negative value of −0.08 ns, and the dashed black line denotes the value of *R*1*i*(*τ*) with *ζ* gets a positive value of 0.08 ns. We can observe that the choice of a negative value of *ζ* makes *R*1*i*(*τ*) get the maximum value at *τ* = 0 and this value is greater than *R*1*i*(0)- with a positive value of *ζ* and also *R*0*i*(0) (both of these functions are not maximized at *τ* = 0). Thus *ζ* will be selected that makes the value of *RG*4(*τ* − *T*PPM − *ζ*) reaches to the minimum point of the *RG*<sup>4</sup> function which is denoted by *Rτopt* as shown in Figure 7. It is clear that, in the UWB-PPM-ATS scheme, the negative values of *ζ* give the better correlation function than positive values. The variation of the correlation function with the different of *ζ* values is illustrated in Figure 8. In Figure 8, the different values of *ζ* lead to the different shapes of the correlation function; the optimal value of *ζ* in this case is *ζopt* = −0.08 ns which makes *R*1*i*(0) have the maximum value. When the magnitude of the *ζ* increases close to the pulse width (*ζ* = −0.16 ns, −0.25 ns in Figure 8), the corresponding correlation functions will not reach the global maximum near the value of *τ* = 0 and its local maximum points have the values close to the global maximum. This leads to the higher error in determining the global maximum point of these functions. Therefore, depending on the specific parameter configuration of each UWB-PPM system, the value of *ζ* will be chosen such that *R*1*i*(0) in the (14) has the largest value.

**Figure 7.** Correlation functions of the conventional UWB-PPM and the proposed UWB-PPT-ATS systems with different time shift.

**Figure 8.** Correlation functions of proposed UWB-PPM-ATS with different additional time shift.

The procedure for estimating distances is based on the traveling time mentioned above and can be applied to a well-known propagation environment; however, it is unavailable for an unknown propagation environment. Consequently, we propose to apply LMFA [28] into the calculation of distances, relative permittivities, and determination of buried object location (depth and horizontal direction) in the case of the unknown propagation environment. The proposed method comes from the fact that the traveling time depends on the relative permittivities, buried object location, and transceiver position. So in order to find the value of these parameters, we change the position of the receiver horizontally and calculate the traveling time corresponding to each position of the transceiver. Based on the transceiver location and computed traveling time, the remaining parameters will be estimated.

#### *3.2. Estimation Algorithm*

Let us consider a specific system model of UWB-PPM systems illustrated in Figure 9, in which the environment is assumed to have a nondestructive structure that has two layers with relative permittivities of *ε*<sup>1</sup> and *ε*2.

**Figure 9.** Target system model in an environment with nondestructive structure of two layers.

To determine the environmental parameters and also the location of the target "T", the transceiver generates UWB pulses with rotation angles of the transceiver antenna are *α<sup>i</sup>* and get the reflected pulses. The traveling time *τ*1*<sup>i</sup>* corresponding to the *i th* time of emission is calculated according to (11), the geometric relationship between *τ*1*i*, *d*1, and *ε*<sup>1</sup> is indicated in (15).

$$\pi\_{1i} = \frac{2d\_1\sqrt{\varepsilon\_1}}{\varepsilon.\cos\alpha\_i},\tag{15}$$

where *α<sup>i</sup>* is the rotation angle of antenna and:

$$
\alpha\_{\bar{i}} = \bar{\imath}. \Delta \kappa;\tag{16}
$$

where Δ*α* is the minimum rotation angle of the transceiver antenna. Thus the values of *ε*<sup>1</sup> and *d*<sup>1</sup> in (15) are estimated based on known pairs (*αi*, *τ*1*i*). Then, to locate buried objects, the device is moved horizontally, and the transmitter emits pulses in the perpendicular direction after every movement step of Δ*x*. Denote *τ*2*<sup>i</sup>* is the traveling time of the transmitted signals, reflected from the buried object and returned to the transceiver at position of *xDi*, it is computed by:

$$\tau\_{2i} = \tau\_{10} + \frac{\sqrt{\varepsilon\_2 (d\_2^2 + (\mathbf{x}\_T - \mathbf{x}\_{Di})^2)}}{c},\tag{17}$$

where *xT* is the horizontal coordinates of the buried object and *τ*<sup>10</sup> is the traveling time of reflected signal at the interface between two layers of the environment in *d* direction. *τ*2*<sup>i</sup>* also computed according to (11), the values of *ε*2, *d*2, and *xT* in (17) are estimated based on known pairs (*τ*2*i*, *τ*10, *xDi*), and

$$
\mathbf{x}\_{Di} = \mathbf{i}.\Delta \mathbf{x}.\tag{18}
$$

The proposed algorithm to estimate the values of those parameters is indicated in Figure 10 and described in detail as below

**Figure 10.** Flowchart of estimated algorithms.

In Figure 10, the unknown parameter vectors are denoted by:

$$\mathbf{X}\_1 = (\hat{\varepsilon}\_1, \hat{d}\_1); \mathbf{X}\_2 = (\hat{\varepsilon}\_2, \hat{d}\_2, \hat{\mathbf{x}}\_T), \tag{19}$$

and the index *m* in Figure 10 is used to replace for 1, 2. With the known pairs (*τ*1*i*, *αi*) and (*τ*2*i*, *xDi*), the unknown parameter vectors **X**<sup>1</sup> and **X**<sup>2</sup> are determined such that the sum of the squares of the deviations *Si*(**X***i*) are minimized:

$$X\_1 \Longrightarrow \arg\min\_{\mathbf{X}\_1} S\_1(\mathbf{X}\_1) = \arg\min\_{\mathbf{X}\_1} \sum\_{i=1}^{M} \left[\tau\_{1i} - f\_1(a\_i, \mathbf{X}\_1)\right]^2,\tag{20}$$

$$X\_2 = \arg\min\_{\mathbf{X}\_2} S\_2(\mathbf{X}\_2) = \arg\min\_{\mathbf{X}\_2} \sum\_{i=1}^{M} [\tau\_{2i} - f\_2(\mathbf{x}\_{Di}, \mathbf{X}\_2)]^2,\tag{21}$$

where *f*1(*αi*,**X**1) and *f*2(*xDi*,**X**2) are:

$$f\_1(a\_i, \mathbf{X}\_1) = \frac{\widehat{d}\_1 \cdot \sqrt{\widehat{\varepsilon}\_1}}{c.cosa\_i},\tag{22}$$

$$f\_2(\mathbf{x}\_{Di'}, \mathbf{x}\_2) = \frac{\sqrt{\hat{\varepsilon}\_2(\hat{d}\_2^2 + (\hat{\mathbf{x}}\_T - \mathbf{x}\_{Di})^2)}}{c},\tag{23}$$

where. denotes the estimated values of the parameters; the vectors **X1** and **X2** were estimated according to the steps as shown in Figure 10 as follows:

+ Assign to **X**1, **X**<sup>2</sup> any arbitrary initialization values, denoted by **X**1−*int*, **X**2−*int*;

+ In each iteration step, the parameter vectors **X**1, **X**<sup>2</sup> are replaced by new estimates (**X**<sup>1</sup> + *δ*1), (**X**<sup>2</sup> + *δ*1) with *δ*1, *δ*<sup>2</sup> are updated step vectors and we have:

$$\begin{array}{llll} \mathbf{S}\_{1}(\mathbf{X}\_{1}+\delta\_{1}) \approx & [\boldsymbol{\tau}\_{1}-\mathbf{f}\_{1}(\mathbf{X}\_{1})]^{T}[\boldsymbol{\tau}\_{1}-\mathbf{f}\_{1}(\mathbf{X}\_{1})] \\ & -2[\boldsymbol{\tau}\_{1}-\mathbf{f}\_{1}(\mathbf{X}\_{1})]^{T}\mathbf{J}\_{1}\boldsymbol{\delta}\_{1}+\boldsymbol{\delta}\_{1}^{T}\mathbf{J}\_{1}^{T}\mathbf{J}\_{1}\boldsymbol{\delta}\_{1} \\ \mathbf{S}\_{2}(\mathbf{X}\_{2}+\delta\_{2}) \approx & [\mathbf{x}\_{D}-\mathbf{f}\_{2}(\mathbf{f}\_{2})]^{T}[\mathbf{x}\_{D}-\mathbf{f}\_{2}(\mathbf{X}\_{2})] \\ & -2[\mathbf{x}\_{D}-\mathbf{f}\_{2}(\mathbf{X}\_{2})]^{T}\mathbf{J}\_{2}\boldsymbol{\delta}\_{2}+\boldsymbol{\delta}\_{2}^{T}\mathbf{J}\_{2}^{T}\mathbf{J}\_{2}\boldsymbol{\delta}\_{2} \end{array}$$

with **J**1,2 is the Jacobian matrix, whose *i th* row equals **J**1*<sup>i</sup>* and **J**2*i*, respectively:

$$\mathbf{J}\_{1i} = \frac{\partial f\_1(\alpha\_i, \mathbf{X}\_1)}{\partial \mathbf{X}\_1},\tag{24}$$

$$\mathbf{J}\_{2i} = \frac{\partial f\_2(\mathbf{x}\_{Di}, \mathbf{X}\_2)}{\partial \mathbf{X}\_2}. \tag{25}$$

To get the minimum value of the sums **S**1, **S**2, the update step vectors *δ*1,2 is calculated such that the derivative of **S**1,2(**X**1,2 + *δ*1,2) with respect to *δ*1,2 has a result of zero, so, *δ*1,2 can be determined satisfying:

$$[\mathbf{J}\_1^T \mathbf{J}\_1 + \lambda \text{diag}(\mathbf{J}\_1^T \mathbf{J}\_1)] \delta\_1 = \mathbf{J}\_1^T [\pi\_1 - \mathbf{f}\_1(\mathbf{X}\_1)],\tag{26}$$

$$\left[\mathbf{J}\_2^T \mathbf{J}\_2 + \lambda \text{diag}(\mathbf{J}\_2^T \mathbf{J}\_2)\right] \delta\_2 = \mathbf{J}\_2^T [\mathbf{x}\_\mathcal{D} - \mathbf{f}\_2(\mathbf{X}\_2)],\tag{27}$$

where the damping factor *λ* (non-negative) is adjusted at each iteration. If **S**1,2 is reduced rapidly, a smaller value of *λ* can be used, whereas if in an iteration does not reduce the residual, *λ* can be increased. The update step vectors are computed as follows:

$$\delta\_1 = [\mathbf{J}\_1^T \mathbf{J}\_1 + \lambda \text{diag}(\mathbf{J}\_1^T \mathbf{J}\_1)]^{-1} \mathbf{J}\_1^T [\pi\_1 - \mathbf{f}\_1(\mathbf{X}\_1)] \tag{28}$$

$$\delta\_2 = [\mathbf{J}\_2^T \mathbf{J}\_2 + \lambda \text{diag}(\mathbf{J}\_2^T \mathbf{J}\_2)]^{-1} \mathbf{J}\_2^T [\mathbf{x}\_D - \mathbf{f}\_2(\mathbf{X}\_2)].\tag{29}$$

After a certain number of iterations, the output of LMFA is the final estimated values of system parameters that meet the constraint condition in (21), as shown in Figure 10.

#### **4. Simulation Results and Discussion**

#### *4.1. Simulation Parameters*

The accuracy of methods for determining the distance and characteristics of a multilayered reflective environments using UWB pulses is strongly dependent on UWB signal processing techniques. UWB-PPM is one of the candidates for positioning technology in a nondestructive environment with multiple reflective layers; with the proposed UWB-PPM-ATS, our trials indicated that the UWB-PPM pulses shifted with a certain time constant can be used to improve the precision in estimating the distance. All the numerical results in this paper were computed using Matlab.

First, we compare the exact estimation errors of the UWB-OOK, UWB-PPM and UWB-PPM-ATS systems with the results obtained from the simulations and the actual value of the parameters use (6), (8), (11), (13) and (30). Here OOK is one of the UWB pulse amplitude modulation techniques and has two level modulation with the bits are 0 and 1. When sending a bit 0, the transmitter will not send anything. When sending a bit 1, the transmitter will send a pulse [29], and the template waveform *ω*(*t*) at the receiver of UWB-OOK system is *g*4(*t*). In addition, the performance of UWB-OOK, UWB-PPM and UWB-PPM-ATS systems were evaluated for an environment with known characteristics (for example, here are three layers as indicated in Figure 1). Finally, the location determination technique in a unknown environment with multiple reflection layers using the LMFA based on the estimated parameters of those UWB systems is presented with the system model as shown in Figure 9. The comparisons between the above systems are evaluated in term of errors between the estimated values from the considered system and the true values. Based on the PSD plot of the Gaussian impulses shown in Figure 4, most of the numerical results presented in this section are based on analysis using the fourth-order Gaussian monocycle. The parameters of the example UWB systems are listed in Table 1 and follow [5].

To assess the performance of those UWB systems, the error of estimated distance is defined as follows:

$$
\delta\_d = \left| \hat{d}\_i - d\_{tri} \right|,
\tag{30}
$$

where ˆ *di* denotes the estimated distance and *dtri* is the true value.


**Table 1.** Simulation parameters.

In Figure 1, without generality, we assume that the propagation environment is heterogeneous with three layers: sand (dry), sandy soil (dry) and granite (dry) and the relative permittivities *ε*1,*ε*2,*ε*<sup>3</sup> are 4, 3, 5, respectively. Those layers are assumed to be dry to reduce the attenuation of environment, so propagation velocity in those layers are 15 cm/ns, 17.32 cm/ns, 13.42 cm/ns and the attenuation is 0.01 dB/m [30]. With a known transmission environment, those distance can be estimated using (6), where *τ<sup>i</sup>* is computed according to (8), (11) and (13). Then the errors of estimated distances are calculated according to (30) and illustrated in Figure 11 for OOK, PPM, ATS systems with *ζ* = −0.08 ns and *ζ* = −0.16 ns. In Figure 11, the black line denotes the average error of the OOK system, the blue line denote the error of the PPM system and the red lines denote the error of the ATS system, all with the same system parameters. Observe that the relative error of the UWB-OOK system is about 24%, of UWB-PPM is about 11%, of UWB-PPM-ATS with *ζ* = −0.08 ns is about 7%, with *ζ* = −0.16 ns is about 13%; so the UWB-PPM- ATS with *ζ* = −0.08 ns performs better for all distance values. This can be explained by comparing the correlation functions which are given in (8), (13) and Figures 7 and 8; we observe that *R*1*i*(0) is significantly greater than *R*0*i*(0) for *ζ* = −0.08 ns and *R*0*i*(0) is not the maximum value of *R*0*i*(*τ*), which leads to a smaller than average error in the PPM-ATS with *ζ* = −0.08 ns when compared to PPM scheme; but with *ζ* = −0.16 ns, *R*1*i*(0) is less than *R*0*i*(0) thus the PPM outperforms the PPM-ATS. Besides that, the UWB-OOK is one of the amplitude modulation techniques, thus it is affected by the transmission environment and more difficult to separate the received signal from noise interference in comparison with the PPM techniques which have a constant amplitude. The results in Figure 11 indicate that the time shift *ζ* directly affects the performance of PPM systems. The value of *ζ* is selected so that (*T*PPM + *ζ*) achieves the optimal value at which *RG*4(*τ* − *T*PPM − *ζ*) gets the minimum value. From (6) and (13), we also know that *ζ* is determined by the pulse shape and the pulse width normalization factor *μp*. Due to the influence of the UWB-PPM system to the choice of the time shift, a suitable value of *ζ* should be chosen according to the specific pulse shape and pulse width employed in a particular UWB application.

**Figure 11.** Comparison of distance estimation errors between OOK, PPM and the proposed UWB-PPM-ATS modulation techniques.

#### *4.2. Determine the Location of Buried Object by LMFA*

As mentioned above, the propagation distance can only be determined using the correlation function if the characteristic of the environment is well-known; however, this is impossible with unknown environment. Moreover, the position of buried object cannot be determined by only the propagation distances; it must be determined in both *x* and *d* directions (*xT*, *dT* in Figure 9). The parameters of each layer in the system model in Figure 9 are estimated sequentially using LMFA as presented in Section 3.2. The parameters of system model and the initialization vectors of the LMFA algorithm are listed in Table 2 in which two layers of environment are assumed to be dry sand and concrete with the corresponding propagation velocity of 16 cm/ns and 14.14 cm/s.


**Table 2.** Initialization parameters of the model.

First the parameters *ε*1, *d*<sup>1</sup> are estimated based on the propagation time *τ*1*<sup>i</sup>* and rotation angle *α<sup>i</sup>* with both UWB-PPM and UWB-PPM-ATS systems. Then, in a similar way, the parameters *d*2, *ε*2, and *xT* are estimated based on **ø2** and **xD**; therefore, the depth of buried object *dT* = *d*<sup>1</sup> + *d*<sup>2</sup> is also determined. The estimated results after 30 iterations and errors of UWB systems are listed in Table 3 and Figures 12–14. Figures 12 and 13 show the relationship between *τ*1*<sup>i</sup>* and *α<sup>i</sup>* in (15). In those figures, the dashed black line denote the curve of *τ*1*<sup>i</sup>* vs. *α<sup>i</sup>* with the true values of *d*<sup>1</sup> and *ε*1; the blue dots are *τ*1*<sup>i</sup>* values measured by simulated UWB systems and the red line denote the curve of *τ*1*<sup>i</sup>* vs. *α<sup>i</sup>* with the estimated values *d* <sup>1</sup> and *<sup>ε</sup>*1. Figure <sup>14</sup> shows the curve of (*τ*2*<sup>i</sup>* <sup>−</sup> *<sup>τ</sup>*10) depends on *xDi* according to (17). In those figures, the dashed black line denotes the curve of (*τ*2*<sup>i</sup>* − *τ*10) vs. *xDi* with the true values of *d*2, *ε*<sup>2</sup> and *xT*; the blue dots are (*τ*2*<sup>i</sup>* − *τ*10) values measured by simulated UWB systems and the red line denotes the curve of (*τ*2*<sup>i</sup>* − *τ*10) vs. *xDi* with the estimated values *d* 2, *<sup>ε</sup>*<sup>2</sup> and *<sup>x</sup>T*.

**Figure 12.** The results of estimating *d*<sup>1</sup> and *ε*<sup>1</sup> of the conventional UWB-PPM based on Levenberg–Marquardt Fletcher algorithm (LMFA).

**Figure 13.** The results of estimating *d*<sup>1</sup> and *ε*<sup>1</sup> of the proposed UWB-PPM-ATS based on LMFA.

**Figure 14.** The results of estimating *d*2, *ε*<sup>2</sup> and *xT* of the conventional UWB-PPM and the proposed UWB-PPM-ATS based on LMFA.

Table 3 and Figure 15 show the estimated results by the UWB systems. As seen, the modulation techniques for UWB pulses combine with nonlinear estimation method LMFA can be used to determine the thickness of different layers in the nondestructive environments, their relative permittivities and also the position of buried object. We observe that the UWB-PPM-ATS system outperforms the UWB-PPM system for all parameters of the model. This behavior results directly from the features of the correlation functions of different signals shown in Figures 7 and 8. When using UWB-PPM-ATS technique with the optimal value of *ζ*, the estimated traveling time values have a smaller error than using UWB-PPM technique, so the results from the LMFA of the PPM-ATS system give the higher accuracy than conventional PPM system.


**Table 3.** The results of the estimated parameters.

**Figure 15.** The true and estimated locations of buried object using the conventional UWB-PPM and the proposed UWB-PPM-ATS.

#### **5. Conclusions**

In this paper, the authors have proposed an UWB pulse position modulation with an additional time-shift technique for positioning applications in nondestructive environments of the UWB systems. The proposed method is utilized to determine the distances and the position of a buried object in heterogeneous environments. The proposed technique's quality is assessed based on the calculation errors compared to the actual values and comparison with the conventional modulation methods. Combining with the Levenberg–Marquardt Fletcher algorithm to determine the buried object location when the propagation environment is unknown, the efficiency of the proposed technique is better than the PPM-UWB system with the same model. Evaluating different UWB systems based on the calculation of errors provides reliable information for choosing modulation schemes with suitable parameters in the design of UWB systems for positioning applications, especially in nondestructive environments.

However, in this work, only one buried object location is determined; multiple buried object estimation algorithm will be researched in our future works. Additionally, the research on parameters of the buried object, such as shape, size, permittivity and so on, is left for the future work.

**Author Contributions:** Conceptualization, N.T.H. and P.T.H.; methodology, N.T.H. and N.L.C.; software, N.T.H.; validation, N.T.H. and P.T.H.; formal analysis, P.T.H. and N.L.C.; investigation, N.T.H.; resources, N.T.H.; data curation, P.T.H.; writing—original draft preparation, N.T.H.; writing—review and editing, N.T.H. and N.L.C.; visualization, N.L.C.; supervision, P.T.H.; project administration, P.T.H. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Conflicts of Interest:** The authors declare that there is no conflict of interest regarding the publication of this paper.

#### **Abbreviations**

The following abbreviations are used in this manuscript:

