Proposal of UWB-PPM with Additional Time Shift for Positioning Technique in Nondestructive Environments

: The ultra-wide band (UWB) technology has many advantages in positioning and measuring systems; however, powers of UWB signals rapidly reduce while traveling in propagation environments, hence detecting UWB signals are difficult. Various modulation techniques are applied for UWB signals to increase the ability for detecting the reflected signal from transmission mediums, such as pulse amplitude modulation (PAM), pulse position modulation (PPM), and so on. In this paper, we propose an ultra-wide band pulse position modulation technique with optimized additional time shift (UWB-PPM-ATS) to enhance the accuracy in locating buried object in nondestructive environments. Moreover, the Levenberg–Marquardt Fletcher algorithm (LMFA) is applied to determine the medium parameters and buried object location simultaneously. The inﬂuences of proposed modulation technique on determining system’s parameters, such as a propagation time, distance, and properties of the medium are analyzed. Calculation results indicate that the proposed UWB-PPM-ATS gives higher accuracy than the conventional one such as UWB-OOK and UWB-PPM in both homogeneous and heterogeneous environments. Furthermore, the LMFA with the proposed UWB-PPM-ATS outperforms the LMFA with the traditional modulation method, especially for unknown propagation environment. OOK, PPM and the proposed UWB-PPM-ATS modulation techniques.


Introduction
With ultra-wide bandwidth, the ultra-wide band (UWB) signal is considered as an ideal locating technique in a short-range with high spatial resolution. As defined by the Federal Communications Commission (FCC), UWB technology has a center frequency higher than 2.5 GHz, or if less than 2.5 GHz, there must be a minimum bandwidth ratio of 0.2 [1], or the minimum bandwidth must reach 500 MHz [2]. To avoid affecting other narrow band systems, the rules of FCC allow the effective isotropic radiated power (EIRP) level of UWB devices to be lower than −41.3 dBm/MHz in the frequencies range of 3.1 to 10.6 GHz [3], so the UWB devices can work for more extended periods than narrow band systems with the same battery power, and due to the use of very narrow pulses, UWB signals are better able to penetrate in the nondestructive environments. General modulation techniques are used for UWB signals such as pulse amplitude modulations (PAMs), On-Off Keying (OOK), and pulse position modulation (PPM) [4]. One can use time-hopping (TH) in UWB systems to create TH-PPM, TH-BPSK signal types [5], or design a generator circuit which generates the 4-th and 5-th order derivative of Gaussian pulses in TH-QPSK system applied to multipath channels [6]. Each modulation technique has a different application range. The choice of the right modulation configuration not only increases the efficiency of system implementation but also maximizes the • Firstly, the UWB-PPM signal which is used in determining the distance under the homogeneous and heterogeneous propagation conditions is mathematically analyzed.
• Secondly, based on the analysis of received signal, the delay time and then the propagation distances are estimated. • Thirdly, to increase the accuracy of the estimated distance, an enhanced UWB-PPM modulation technique, called ultra-wide band pulse position modulation with additional time shift (UWB-PPM-ATS) is proposed and compared with other techniques such as conventional UWB-PPM, UWB-OOK with the properties of transmission environment is known. • Finally, the Levenberg-Marquardt nonlinear estimation algorithm is applied to estimate the system parameter and the target when the propagation environment is unknown.
The remainder of this paper is organized as follows. Section 2 describes the system model for estimating distances in nondestructive environments. The proposed system model and the parameters are presented in Section 3. The simulation results are provided in Section 4, and finally, conclusions and further work are discussed in Section 5.

System Model
A distance measurement and positioning system using UWB technology is illustrated in Figure 1, where d i and ε i 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 where A i denotes the amplitude factor which represents the reflection and transmission properties of the propagation environment, and its value depends on the properties of the ith layers, τ i is the traveling time of UWB pulse in the ith layer. A ob 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.  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]: where A p denotes the amplitude of pulse, µ p 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 µ p is reduced. The n th derivative of Gaussian pulse, named n th -order monocycle, is where A np is the amplitude of the n th -order monocycle. The shapes of different types of Gaussian monocycles are indicated in Figure 3. 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.  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: where T r is pulse repetition period and T PPM is the time shift associated with binary PPM, a n ∈ {0, 1}, and N is the length of UWB pulses. In this modulation technique, when a n 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 a n is 1. In our system model, a n 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: 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, r i (t), and the template waveform ω(t): 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 0i and has the form: Define the autocorrelation function of the fourth-order Gaussian monocycle in (4) as R G4 (τ), we have 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 where max(R i ) is the local extreme value corresponding to the reflected signals from the surface of the layers or object. 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.

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: 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 + ζ). 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 1i and determined as: So, we have: The shapes of R G4 (τ),R 0i (τ) and R 1i (τ) are shown in Figure 7. Figure 7 indicates the autocorrelation function of the fourth-order Gaussian monocycle R G4 and correlation functions of UWB-PPM R 0i and UWB-PPM-ATS R 1i 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 0i (τ), the red line denotes the value of R 1i (τ) with ζ gets a negative value of −0.08 ns, and the dashed black line denotes the value of R 1i (τ) with ζ gets a positive value of 0.08 ns. We can observe that the choice of a negative value of ζ makes R 1i (τ) get the maximum value at τ = 0 and this value is greater than R 1i (0)with a positive value of ζ and also R 0i (0) (both of these functions are not maximized at τ = 0). Thus ζ will be selected that makes the value of R G4 (τ − T PPM − ζ) reaches to the minimum point of the R G4 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 1i (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 1i (0) in the (14) has the largest value.  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.

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 ε 1 and ε 2 . 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 α i and get the reflected pulses.
The traveling time τ 1i corresponding to the i th time of emission is calculated according to (11), the geometric relationship between τ 1i , d 1 , and ε 1 is indicated in (15).
where α i is the rotation angle of antenna and: where ∆α is the minimum rotation angle of the transceiver antenna. Thus the values of ε 1 and d 1 in (15) are estimated based on known pairs (α i , τ 1i ). 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 τ 2i is the traveling time of the transmitted signals, reflected from the buried object and returned to the transceiver at position of x Di , it is computed by: where x T is the horizontal coordinates of the buried object and τ 10 is the traveling time of reflected signal at the interface between two layers of the environment in d direction. τ 2i also computed according to (11), the values of ε 2 , d 2 , and x T in (17) are estimated based on known pairs (τ 2i , τ 10 , x Di ), and The proposed algorithm to estimate the values of those parameters is indicated in Figure 10 and described in detail as below In Figure 10, the unknown parameter vectors are denoted by: and the index m in Figure 10 is used to replace for 1, 2. With the known pairs (τ 1i , α i ) and (τ 2i , x Di ), the unknown parameter vectors X 1 and X 2 are determined such that the sum of the squares of the deviations S i (X i ) are minimized: where f 1 (α i , X 1 ) and f 2 (x Di , X 2 ) are: where . denotes the estimated values of the parameters; the vectors X 1 and X 2 were estimated according to the steps as shown in Figure 10 as follows: + Assign to X 1 , X 2 any arbitrary initialization values, denoted by X 1−int , X 2−int ; + In each iteration step, the parameter vectors X 1 , X 2 are replaced by new estimates (X 1 + δ 1 ), (X 2 + δ 1 ) with δ 1 , δ 2 are updated step vectors and we have: with J 1,2 is the Jacobian matrix, whose i th row equals J 1i and J 2i , respectively: 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: 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: 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.

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: whered i denotes the estimated distance and d tri is the true value. 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 , ε 3 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 τ i 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 1i (0) is significantly greater than R 0i (0) for ζ = −0.08 ns and R 0i (0) is not the maximum value of R 0i (τ), 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 1i (0) is less than R 0i (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 R G4 (τ − 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.

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 (x T , d T 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. First the parameters ε 1 , d 1 are estimated based on the propagation time τ 1i and rotation angle α i with both UWB-PPM and UWB-PPM-ATS systems. Then, in a similar way, the parameters d 2 , ε 2 , and x T are estimated based on ø 2 and x D ; therefore, the depth of buried object d T = d 1 + d 2 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 τ 1i and α i in (15). In those figures, the dashed black line denote the curve of τ 1i vs. α i with the true values of d 1 and ε 1 ; the blue dots are τ 1i values measured by simulated UWB systems and the red line denote the curve of τ 1i vs. α i with the estimated values d 1 and ε 1 . Figure 14 shows the curve of (τ 2i − τ 10 ) depends on x Di according to (17).
In those figures, the dashed black line denotes the curve of (τ 2i − τ 10 ) vs.
x Di with the true values of d 2 , ε 2 and x T ; the blue dots are (τ 2i − τ 10 ) values measured by simulated UWB systems and the red line denotes the curve of (τ 2i − τ 10 ) vs.
x Di with the estimated values d 2 , ε 2 and x T .   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.

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.