PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of sensorsMDPI Open Access JournalsMDPI Open Access JournalsThis articleThis JournalInstructions for authorssubscribe
 
Sensors (Basel). 2017 November; 17(11): 2554.
Published online 2017 November 6. doi:  10.3390/s17112554
PMCID: PMC5713469

A Robust High-Accuracy Ultrasound Indoor Positioning System Based on a Wireless Sensor Network

Abstract

This paper describes the development and implementation of a robust high-accuracy ultrasonic indoor positioning system (UIPS). The UIPS consists of several wireless ultrasonic beacons in the indoor environment. Each of them has a fixed and known position coordinate and can collect all the transmissions from the target node or emit ultrasonic signals. Every wireless sensor network (WSN) node has two communication modules: one is WiFi, that transmits the data to the server, and the other is the radio frequency (RF) module, which is only used for time synchronization between different nodes, with accuracy up to 1 μs. The distance between the beacon and the target node is calculated by measuring the time-of-flight (TOF) for the ultrasonic signal, and then the position of the target is computed by some distances and the coordinate of the beacons. TOF estimation is the most important technique in the UIPS. A new time domain method to extract the envelope of the ultrasonic signals is presented in order to estimate the TOF. This method, with the envelope detection filter, estimates the value with the sampled values on both sides based on the least squares method (LSM). The simulation results show that the method can achieve envelope detection with a good filtering effect by means of the LSM. The highest precision and variance can reach 0.61 mm and 0.23 mm, respectively, in pseudo-range measurements with UIPS. A maximum location error of 10.2 mm is achieved in the positioning experiments for a moving robot, when UIPS works on the line-of-sight (LOS) signal.

Keywords: ultrasonic indoor positioning system, wireless sensor network, time-of-flight for ultrasonic signal, least squares method

1. Introduction

Satellite signaling systems are universally well known. One example is the global positioning system (GPS), which has a high positioning accuracy at 0.15 m [1]. It is very convenient to use satellites to obtain locations. Unfortunately, satellite positioning signals normally cannot be received indoors due to obstructions from buildings, and people have to do most of their work (an estimated 87 percent overall) indoors [2].

Because of the reasons outlined above, indoor positioning systems (IPSs) have gained increasing attention from researchers across the globe in the last few decades, as IPSs are widely used in health-care as well as in home and building automation [3,4,5,6]. There are already many local positioning systems based on different technologies, and they each have their own advantages and weaknesses. Systems based on infrared beams make use of the scanning sweep method and are very inexpensive for users, but only have accuracy up to 57 cm [7]. Xiaohan Liu et al. [8] developed an IPS using a visible light communication platform, which has an accuracy of 10 cm, but it requires expensive infrastructure. The ultra-wide-band (UWB) is a great option for IPS and has a strong anti-interference capability, but it needs more infrastructure and is expensive for users [9]. IPS based on computer vision [10] has the highest precision of the all technology proposals, but it is too expensive and lacks wide applications. The other local positioning solutions consist of radio frequency identification (RFID) [11], WiFi [12], and ZigBee [13], and they have a low positioning accuracy.

Compared to these above technical solutions, a positioning system based on the ultrasonic signal and a wireless sensor network (WSN) has many advantages, including slow propagation speed (about 340 m/s), so that the hardware for the received ultrasonic signals can be simple, obtaining cost savings. The time-of-flight (TOF) of the signal from a transmitter device to a receiver is used to calculate the distance in the UWB and ultrasonic solutions, but a lower time synchronization accuracy between different nodes can be tolerated for the ultrasound’s speed, compared to the speed of electromagnetic wave propagation in the UWB system. The ultrasound signal is weakly influenced by the surroundings and has a negligible penetration of walls, and this can enhance its practical value in the proposed ultrasonic indoor positioning system. Besides, ultrasonic transducers can be integrated easily into the mobile phone system [14]. This is almost equivalent to adding a new human organ, and the location of the person is conveniently completed by it.

There are already many methods of ultrasonic signal processing for distance measurement, including non-destructive testing [15,16]. The envelope signal extracted by these methods has some harmonics that cannot be ignored, and this will affect the measurement accuracy if they are used in ultrasonic indoor positioning systems (UIPSs). The Cricket location system [17] detects the arrival of the ultrasound with a threshold value. It is simple, but has low precision due to noise interference. The random signal processing method is used to estimate ultrasonic TOF in some works in the literature [18,19], considering environmental noise, but these algorithms are too complex for UIPSs based on WSNs. The ultrasonic signals are encoded into various digital codes, such as binary frequency shift-keyed and code division multiple access codes [20,21,22,23]. These methods can indeed improve the stability of system, but are difficult to decode because of ultrasonic vibration inertia. The acoustical signals are modulated into four sets of orthogonal chirps with multiple access to solve signal interference in [24], but the chirp signals in conjunction with cross correlation are very complex, and the time resolution of the system for position will be reduced. The frequency-hopping spread spectrum (FHSS) has been introduced into the field of IPS to solve the multiple-access problem [25,26], but the ultrasonic transducer has a narrow frequency bandwidth, which will increase the difficulty and cost of realizing FHSS. Quadrature sampling is used to detect the ultrasonic envelope with a low frequency [27], which is appropriate for the WSN system and can save power, but it has low robustness. Power is not a problem. Considering that the positioning system works indoors, a new envelope detection method (EDM) is proposed in this paper. The EDM can obtain a high accuracy and robustness in UIPSs.

The paper is organized as follows. In Section 2, the UIPS architecture is described. A detailed account of the principle of the EDM is given in this section. Section 3 describes the algorithm to measure the TOF of the ultrasonic signal and the location of the target. The implementation of the UIPS is shown in Section 4, including the hardware and software. Section 5 gives the experiment results for the EDM in the simulation and practical applications. Finally, the main conclusions of the whole work are summarized in Section 6.

2. Global Description of the UIPS

2.1. UIPS Architecture

The UIPS has many WSN anchor nodes that are placed at known positions of the indoor space, as shown in Figure 1. The anchor nodes are on the interior ceiling, and there is a WSN node on the mobile robot system, which is the positioning target. The distance values from the target node to every anchor node and their address and time information will be transmitted to the server via wireless WiFi, when all distance measurements are finished. The coordinates for the mobile robot can be calculated by the server when the distance values to all anchors in the same moment are received. The server can forward these coordinates on to users (robots here).

Figure 1
The architecture of the ultrasonic indoor positioning system (UIPS).

There can be two work modes to transmit the ultrasonic signal: the first one is from the anchor nodes to the robot, shown in Figure 2a, and the second is vice versa, as shown in Figure 2b.

Figure 2
The two work modes for the ultrasonic signal: (a) and (b), Tx: transmitter, Rx: receiver.

We assume that the number of the anchors is N1 (N1>3), and the number of robots is N2 (N2>0) in a coverage area. Every anchor will send the ultrasonic signal separately at a different time, presented as t1, t2, (...) in Figure 2a. Thus, the signal interference between different anchor nodes is avoided when the time synchronization signal (TSS) arrives in the first work mode. There is some time drift between different anchors, for which there must be a compensation algorithm, but it is very complex, and the accuracy is not high. The advantage with this mode is that it needs only one TSS for positions of all robots, and saves power. If just one anchor transmits the signal in a TSS in order to void the compensation for time drift, there needs to be N1 TSSs for all robots, but the position will have large errors in both methods above when robots are moving. Because the measurements for the distances to all anchors are finished at different times, the position coordinates can only be calculated roughly through these distance values.

In the second work mode, the target nodes transmit the ultrasonic signal to the anchors, selected by the UIPS in this paper. Only one robot sends the signal that will be received by all anchors at the same time when the TSS for it arrives, and then the distances to all anchors are measured. The ultrasonic signal pulse is just 1 ms, which will be described later in this article. There is a reasonable assumption that the position for the moving robot will not change in so short a time, so this mode can be appropriate for mobile robot localization, although it needs more TSSs (Nr) for all robots. The information from the time synchronization node contains a robot address. The robot with this address will send the ultrasonic signal, and the remaining ones will not. The anchors distinguish the ultrasonic signals from each other with the address. The time synchronization node will cyclically transmit TSS to every robot in turn, so the resolution for UIPS will reduce when Nr grows larger.

2.2. Ultrasonic Envelope Detection

Ultrasound is a vibration wave, for which the amplitude will increase when it arrives at the peak. There will be a reducing process while the signal drops to the minimum amplitude value from the peak value. If the driving signal for the ultrasound is cut off immediately when the peak arrives, the signal will have a unique wave peak. The signal for this means can be described by Equation (1) [28].

x(t)=A(t)·sin(2πfct+θ(t))
(1)

where the parameters A(t), fc, and θ(t) are the signal amplitude, frequency, and phase, respectively. A(t) is also the envelope signal and is given by:

A(t)=βeα(tτ)2
(2)

where β is the amplitude parameter, and α is the bandwidth factor. We can see that A(t) will obtain its maximum value when t=τ from the equation.

The original ultrasonic signal x(t) from the received transducer is shown in Figure 3a, and the envelope A(t) is presented in Figure 3b. The ultrasonic envelope detection aims to extract A(t) from x(t). We consider that the time when the ultrasonic envelope arrives at its peak is the arrival time of the ultrasonic signal. That means the TOF for the ultrasonic envelope is τ in Equation (2), if t=0 is the start time for transmission of the signal. A rectangular window function w(t) with a length of N is defined as Equation (3). It will become w(tT), w(t2T), and w(t3T), (...), when the function is moving with a step of T, as shown in Figure 3c. T is the sampling period for the ultrasound signal. x(t) of the length of N·T is selected in every window, when x(t) multiplies the window function.

w(t)=10<tNT0tisothers.
(3)
Figure 3
(a) Ultrasonic signal; (b) Envelope; and (c) Rectangular window.

The digitalized values x(n) are provided by the analog to digital converter (ADC) from x(t) with a sampling rate of fs (1/T), given by:

x(n)=A(n)·sin(wcnT+θ(n))=A(n)·sin(wcnT)cos(θ(n))+A(n)·cos(wcnT)sin(θ(n))=sin(wcnT)cos(wcnT)A(n)cos(θ(n))A(n)sin(θ(n))
(4)

where A(n) and θ(n) stand for the amplitude sampling value of the envelope and the phase angle at nT time, respectively, and wc=2πfc.

If the window function w(t) is discretized into w(n) with the same frequency fs, and we can get Equation (5) in the window i,

x˜(i+j)=x(i+j)·w(ni+j)
(5)

where j=1,2,,N and i=0,1,2,. Because the value of N·T is small, at just 0.05 ms in this paper (where N=50 and T=1μs), ultrasonic signal is a mechanical vibration, with the characteristic of inertia. Its amplitude and phase have little change in such a short period of time. We can assume that the envelope value A(i+j) and the phase θ(i+j) are constant for all j values in every window i, presented as A˜(i) and θ˜(i), respectively, and then we can obtain Equation (6) from Equations (4) and (5), where x˜(i+1), x˜(i+2), ,x˜(i+N) can be obtained from ADC sampling.

x˜(i+1)=sin(wc(i+1)T)·A˜(i)cos(θ˜(i))+cos(wc(i+1)T)·A˜(i)sin(θ˜(i))x˜(i+2)=sin(wc(i+2)T)·A˜(i)cos(θ˜(i))+cos(wc(i+2)T)·A˜(i)sin(θ˜(i))x˜(i+3)=sin(wc(i+3)T)·A˜(i)cos(θ˜(i))+cos(wc(i+3)T)·A˜(i)sin(θ˜(i))=+x˜(i+N)=sin(wc(i+N)T)·A˜(i)cos(θ˜(i))+cos(wc(i+N)T)·A˜(i)sin(θ˜(i))
(6)

Equation (6) is written in the form of a matrix:

Xi=Ψi·Ci
(7)

where the transposes of Xi, Ci and Ψi are given by:

XiT=x˜(i+1)x˜(i+2)x˜(i+3)x˜(i+N)
(8)

CiT=A˜(i)cos(θ˜(i))A˜(i)sin(θ˜(i))
(9)

ΨiT=sin(wc(i+1)T)sin(wc(i+2)T)sin(wc(i+N)T)cos(wc(i+1)T)cos(wc(i+2)T)cos(wc(i+N)T)
(10)

Now, the problem is transformed into a method to estimate Ci of Equation (7). The estimation of Ci is C^i (assumed C^i=c^i1c^i2), and it can be determined using the least squares method for minimizing the following error function between Xi and X˜i (X˜i=ΨiCi^):

J(Ci^)=(XiΨiCi^)T(XiΨiCi^)
(11)

We can obtain Equation (12) in order to minimize J(Ci^):

JCiCi=Ci^=2ΨiT(XiΨiCi^)=0
(12)

and then, Ci^ is given by:

Ci^=(ΨiTΨi)1ΨiTXi
(13)

we can get A^(i) and θ^(i), which are the estimations of A˜(i) and θ˜(i), respectively, from Equations (9) and (13):

A^(i)=c^i12+c^i22
(14)

θ^(i)=arctanc^i2c^i1
(15)

3. Implementation of Positioning

The received nodes must know the time at which the transmitter node starts to generate the ultrasonic signal. Hence, there must be a time synchronization in these nodes, which is presented in the first part of this section, followed by TOF estimation and positioning.

3.1. Time Synchronization

The time protocol of the UIPS is based on reference broadcast synchronization (RBS) [29], which works in the data link layer. The architecture of synchronization for the UIPS is shown in Figure 4.

Figure 4
The architecture of time synchronization among different nodes.

The whole space is divided into several subspaces, and there is a reference node in every subspace. All anchor nodes of the subspace have a common ultrasonic signal coverage region, and the RBS node will be asked to broadcast the synchronization signal periodically to the anchor nodes and the target nodes by the server, when there are some positioning nodes that walk into the region. The transmitters will start to work, and the received nodes will begin sampling the signal as the synchronization broadcast is arriving.

The process of the transmission for TSS is shown in Figure 5. The central processing unit (CPU) of the reference node will start a timer count when it receives instructions from the server. The timer will generate an interrupt signal and count from zero again while the counter arrives at the TSS period, and then the CPU program will access the interrupt handing, where the synchronization frames are transmitted to the RF module. The frames will be encoded in the RF module and then transmitted to the receivers by the antenna. Receiving and decoding will be completed by the receivers with a reversed order of the sender when the TSS arrives at the receiving systems. Furthermore, the alignment for the received frames is finished and then an interrupt signal is generated for the CPU.

Figure 5
The process of the transmission for the time synchronization signal (TSS). RF: radio frequency.

The velocity of the electromagnetic wave is the approximate speed of light in a vacuum, so the propagation time differences from different nodes are quite minor and can be ignored in a subspace. The synchronization protocol based on RBS suffers from the uncertainties of the transmission and reception times, as well as the yield synchronization error, which is about 13 μs [30], and cannot meet the requirements of the UIPS. Several modifications are performed to reduce the error. First of all, the synchronization and other data transmissions have the same communication module as in the other system, but there are two different modules in the UIPS, and the modules will not be occupied by the data transmission process. The interrupt for transmitting and receiving have the highest priority, and will not be terminated by the other interrupt.

As described above, the rest factors for synchronization accuracy are in the RF module. CC2500 is selected as a synchronization module in UIPS, which can provide extensive hardware support for packet handling, data buffering, and burst transmissions [31]. Encoding/decoding, transmitting/receiving, and alignment are affected by the traffic rate of the synchronization module. A testing experiment is designed in order to verify this and calculate accuracy. The reference node transmits eight bytes of synchronization data to the other node. The structure for the data frame is shown in Figure 6.

Figure 6
The structure for the data frame, CRC: cyclic redundancy check.

The length of the address is two bytes in the figure, and a command is used to distinguish data frames into time synchronization or data transmission. This takes up one byte. The next three bytes are data bits. Because data transmission is completed by WiFi module, the data bits are reserved for future research. The last two bytes are cyclic redundancy check (CRC) bits, which indicates whether the received data has some errors. There are two nodes in the experiment, the data frames are sent once every 200 ms, the receiver will return the data to the reference immediately when the receiving is finished. The time from the beginning of transmission to the end of receiving for the reference node is recorded as ts in the sender. The statistics of 500 measurements for every rate are shown in Table 1.

Table 1
The statistics of the time synchronization experiment with different communication rates.

The difference between the maximum and minimum value of ts can indicate the precision of the time synchronization. This brings us to the conclusion that the higher the data communication rate is, the lower the differences are. The difference will be 1 μs when the rate is above 100 kb/s, so the precision of the time synchronization should be less than 1 μs because of one-sided communications in UIPS. The ranging measurement error, resulting from this, is less than 0.34 mm, assuming the ultrasonic velocity is 340 m/s.

3.2. Ultrasonic TOF Estimation

The analog-to-digital converter (ADC) of the received nodes starts to sample the ultrasound signal using a sampling frequency fs (fs=1T) when the synchronization data are arriving. The range of ultrasonic propagation is limited to 5 m in UIPS, so the time it takes to sample the signal is less than 15 ms, assuming the velocity is 340 m/s. The uncertain drift of the crystal oscillator of the received nodes can be ignored in this period. Therefore, the time instant at which the envelope reaches its maximum values (nmax) can be calculated with nmaxfs.

The parabolic interpolation method in which the adjacent samples around A^(nmax) can reduce the TOF estimation (ttof) error [32], is given by:

ttof=nmax+nxfs
(16)

where nx is the fractional correction term obtained from parabolic interpolation with three samples [27], given by:

nx=A^(nmax1)A^(nmax+1)2·A^(nmax1)2A^(nmax)+A^(nmax+1)
(17)

The phase information at nx can be obtained from Equation (15), and it can improve a greater resolution of ttof. The TOF calculated value by phase is shown in Equation (19) [20], where Int[] is the integer operation.

ttof=Intnmax+nxfs+θ^(nmax)2πfs
(18)

3.3. Estimation of Distances and Positioning

The speed of ultrasound v in air depends on many factors. The influence of temperature is only considered in the special indoor condition, and it can be estimated with Equation (19) [33], where ϕ is the temperature in degrees Celsius,

v=(331.3+0.6ϕ)m/s
(19)

Then, the distance di between the target node and the anchor node i can be estimated with ttof and v as:

di=v·ttof
(20)

The coordinate of the anchor node i is (xi,yi,zi), and it is assumed that the number of the anchor nodes is k. The coordinate of the target node is (x,y,z). Then, these equations are obtained as:

(xx1)2+(yy1)2+(zz1)2=d12(xx2)2+(yy2)2+(zz2)2=d22(xx3)2+(yy3)2+(zz3)2=d32++=(xxk)2+(yyk)2+(zzk)2=dk2
(21)

where the left and right equations are subtracted from the first and second equations, respectively, and the equation is expressed in matrix form as:

A·x=b
(22)

where A, x and b are referred to as:

A=2(x1x2)2(y1y2)2(z1z2)2(x1x3)2(y1y3)2(z1z3)2(x1x4)2(y1y4)2(z1z4)2(x1xk)2(y1yk)2(z1zk)
(23)

x=[xyz]T
(24)

b=x12x22+y12y22+z12z22+d22d12x12x32+y12y32+z12z32+d32d12x12x42+y12y42+z12z42+d42d12x12xk2+y12yk2+z12zk2+dk2d12
(25)

The close solution of Equation (22) for x can be obtained by using the least squares method from Equation (26) [27].

x=(AT·A)1·(AT·b)
(26)

4. Implementation of the UIPS

The implementation of UIPS is based on embedded system design. A 32-bit ARM Cortex™-M7 microprocessor is selected as the hardware, which features a single floating point unit (SFPU) and supports single precision data processing instructions and data types. The software is programmed in a bare metal machine in which the registers of hardware can be read or written directly. The real-time performance of the system is strengthened significantly.

4.1. Hardware

The node hardware of UIPS is designed for the core board with a size of 90 mm × 58 mm, shown in Figure 7. The core board is powered by a lithium battery (LB), which can be charged for reuse. It can be configured as a receiver anchor node and arranged in the interior ceiling. The specified extended board for the core board that is configured as the transmitter node can be used to control the mobile robot.

Figure 7
The photo of the UIPS node.

The microcontroller unit (MCU), the transmitter and receiver for ultrasound, the CC3200 and CC2500 modules, and the power system comprise the hardware system of UIPS, which is shown in Figure 8. UIPS has two STM32F722RCT6 32-bit MCUs, which have a running frequency up to 216 MHz and come from STMicroelectronics (http://www.st.com); one is the master, and the other is the slave with alternate sampling and calculation. This MCU integrates 512 Kbytes Flash memory with protection mechanisms and 256 Kbytes of static random access memory (SRAM) on chip. It has an abundance of peripheral resources and interfaces, such as a 12-bit ADC, a digital to analog converter (DAC), a timer that can produce the pulse width modulation (PWM) signal, universal synchronous or asynchronous receiver transmitters (USARTs) and a direct memory access (DMA) channel [34].

Figure 8
The hardware architecture of the UIPS node. ADC: analog-to-digital converter; LB: lithium battery; MCU: microcontroller unit; PWM: pulse width modulation; SPI: serial peripheral interface.

The signal strength of the ultrasonic transmitted transducer with a frequency of 40 kHz is too low to be sampled by the ADC directly, so it must be amplified firstly, and then processed by a band-pass filter to reduce the out-of-band noise. The ultrasonic received sensor can be driven by the PWM orsinusoidal signal, which can be generated by DAC. The former is adopted in UIPS. It can be obtained from the timer with high precision.

The CC3200 [35] module communicates with the master MCU by USART, which can realize the transparent data transmission between the master and the server via the WiFi wireless network. The CC2500 module is used for time synchronization and communicates with the master MCU through the SPI interface. Besides, the module is linked to the master MCU via an input/output (IO) port in hardware.

4.2. Software

The implementation of the algorithm for envelope detection in Section 2.2 requires a large number of matrix calculations using Equation (13), but some regularity can be found if Equation (10) is plugged into Equation (13). The elements of Ψi will be periodic and constant if the length N and sampling period T are constant, so all (ΨiTΨi)1ΨiT will also form a series of periodic coefficients. The coefficients in only one period are enough for the implementation of the algorithm, and are computed in advance using MATLAB and are written into the read-only memory (ROM) of the MCU in a table form. The computational complexity of the algorithm can be reduced significantly.

There are rich DMA channels and interrupt vectors on the MCU chip, with which the real-time and reliability performance of the software system can be enhanced. There is a joined IO port between the CC2500 module and the master MCU as described above. The module will produce a rising edge pulse signal when the time synchronization data come, which makes the software of the master go into an external interrupt vector with the highest priority of all. The system will have an alternating sampling and detection between the master MCU itself and the slave MCU, shown in Figure 6. The master will turn on the DMA channel and enable the timer to trigger the ADC for ultrasound signal acquisition with rate fs if it is time for the master itself. The software will run into the DMA interrupt service program when all M pieces of data are collected.

The envelope detection method presented in Section 2.2 is implemented in the DMA interrupt service program and then calculates the distance value from the target node to the anchor node; Hi=(ΨiTΨi)1ΨiT in the figure. The value will be transferred to the server through the WiFi network. The slave MCU will obtain the distance value like the master in the external interrupt program if it is its time, shown in the red dashed rectangle of the figure, and then transfer the value to the master through USART; finally, the master transfers it to the server. The whole software flow diagram is shown in Figure 9.

Figure 9
The software architecture and flow diagram of the UIPS receive node. DMA: direct memory access. IO: input/output.

5. Experiments

A simulation experiment using MATLAB is described to prove the validity of the envelope detection method in this section, meanwhile, the point-to-point ranging and the position for the mobile robot are presented to validate the UIPS’s performance.

5.1. Simulation

The analog ultrasonic signal and its envelope are generated from Equations (1) and (2) in Section 2.2, which only has a peak, as mentioned earlier. Noise is inevitable in any practicable system, so the white Gaussian series are added to the signals with a different signal-to-noise ratio (SNR) in order to test the filtering effect of the algorithm in Section 2.2. The Hilbert transform (HT) is widely used for envelope detection in engineering [36], by which the ultrasonic envelope is also extracted for comparison with the method in this paper. The results are shown in Figure 10.

Figure 10
The comparison for envelope detection based on the Hilbert transform (HT) and least squares method (LSM).

The envelope detection method based on LSM can detect the envelope of ultrasound with a high precision, comparing (a2) and (a4) in the figure, and there is no time delay for the peak between them. Filter effect analysis proves that the proposed algorithm can obtain a good filtering result, even if the SNR is very low (as shown from (b) to (d) of Figure 10). The method based on HT can detect the envelope too, but it has no filter effect from the third column on. This can only be achieved by means of some special digit filter. Besides, the HT is more complex.

5.2. Applications

The distance measurement experiment is completed to evaluate the accuracy of UIPS. A pair of nodes are at the same level as each other, which can guarantee that the nodes work on the LOS condition, avoiding interference for ultrasonic multi-path propagation. One is configured as the transmitter, and the other node as a receiver. The transmitter sent a PWM of 1 ms under the control of the time synchronization signal with a period of 200 ms. The ultrasonic signal completely disappears in such time, so that the influence of the last cycle will be avoided.

A commercial vision-based localization system (vicon) from Oxford Metrics Limited (www.vicon.com) is used to measure the separation between nodes, which has a precision of 1 mm. Firstly, some optics markers are planted around the ultrasound transducer symmetrically, shown in Figure 11, so that the position coordinates of two nodes can be obtained from the vicon system. They are (x1,y1,z1) and (x2,y2,z2), and the distance is:

d=(x1x2)2+(y1y2)2+(z1z2)2
(27)
Figure 11
An illustration of the distance measurement between two nodes.

Eight group distances are measured, which have some random bits. The quadrature sampling (QS) [27] is achieved by UIPS to compare the method based on LSM in this paper. The statistics are shown in Table 2, where X¯ is an average of 500 times of measurement for each distance, var is their variance, and |E|max is defined as:

|E|max=12(vmaxvmin)
(28)

where vmax and vmin represent respectively the maximum and minimum in all 500 measurements. |E|max can reflect a certain accuracy of UIPS.

Table 2
The statistics of distance measurements. QS: quadrature sampling.

We can see that X¯ error is less than 1.5 mm with the reference distance for two methods. The highest |E|max and var can obtain values of 0.61 and 0.23, but they will increase with distance. This is because of the influence of noise when the ultrasonic signal becomes weakened, and the SNR is smaller. |E|max and var in LSM are smaller than those in QS in all groups, so the method based on LSM has a higher accuracy and robustness.

The localization experiment for the mobile robot is performed in order to test the system function. Four anchor nodes are placed at a fixed position, which are configured as the receiving end, shown in Figure 12. Here, the coordinates of the anchor node are obtained from vicon as above. They are, respectively: (115.79,1004.64,1501.64), (1024.88,44.74,1685.52), (18.60,966.18,1577.05), and (968.89,44.40,1914.52). Each transducer of the anchor node has common coverage areas where the signal from the robot can be received by all. A and b can be calculated with Equations (23) and (25), given by

A=1818.182098.76367.76268.783941.64150.822169.362098.08825.76
(29)

b=[615.73+d22d12143.30+d32d121328.47+d42d12]T
(30)
Figure 12
An illustration of the positioning experiment for the mobile robot.

The motion with open-loop trajectory for the robot is a standard circle, which starts from different positions with various radii. There are also some optics markers around the ultrasound transducer of the node on the robot, so the trajectory can be obtained from vicon. The mobile robot will transmit the ultrasonic signal in each 200 ms cycle as with the distance measurement. The signal is transmitted with LOS in order to overcome the multi-path case, additionally, a short ultrasonic signal pulse, just 1 ms in this paper, can reduce interference for ultrasonic multi-path effect [37], so it is ignored here. The robot moves with a very low speed in the common coverage areas, so the Doppler effect of ultrasound is ignored in this experiment.

Four experiments are conducted with different radii in total, and the trajectory tracking results are shown in Figure 13.

Figure 13
The comparison of the robot motion with different radii circular trajectories from (ad). The position data sampling rate from vicon is much greater than that of UIPS.

We can see that the UIPS can follow the tracks of the robot with a maximum error of 10.24 mm. The error is 12 mm when the same experiments are completed though quadrature demodulation [38] with the same sampling frequency fs. |E|m is used to represent the maximum error in this paper, and is defined as:

|E|m=(xmxr)2+(ymyr)2+(zmzr)2
(31)

where (xm,ym,zm) is a position from UIPS, which deviates furthest from the vicon reference. (xr,yr,zr) is from vicon, because vicon is not in time synchronization with UIPS and has a much higher sampling frequency for positioning data than that of UIPS. (xr,yr,zr) is the nearest point to (xm,ym,zm) in all cases.

6. Conclusions

In this work, we present a new envelope detection method for the ultrasonic signal and develop an indoor position system based on this method, which can be also applied to other fields. Filter effect analysis from simulation proves that the method is effective. The system can fix a continuous position on the mobile robot with high and robust stability. We also present an improved time synchronization method based on RBS in the UIPS, which has an accuracy up to 1 μs. Network nodes have individual communication modules using the WiFi protocol, which is used widely in day-to-day life. Therefore, we will integrate the ultrasonic transducer with a mobile phone in the future and realize the precise positioning for personnel using the phone and the UIPS.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant Nos. 61333003, 61690212).

Author Contributions

Author Contributions

Jun Qi and Guo-Ping Liu conceived and designed the experiments; Jun Qi performed the experiments and analyzed the data; Guo-Ping Liu provided experimental platform and funds; Jun Qi wrote the paper; Guo-Ping Liu provided guidance and polishing for the paper.

Conflicts of Interest

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Nickitopoulou A., Protopsalti K., Stiros S. Monitoring dynamic and quasi-static deformations of large flexible engineering structures with GPS: Accuracy, limitations and promises. Eng. Struct. 2006;28:1471–1482. doi: 10.1016/j.engstruct.2006.02.001. [Cross Ref]
2. Klepeis N.E., Nelson W.C., Ott W.R., Robinson J.P., Tsang A.M., Switzer P., Behar J.V., Hern S.C., Engelmann W.H. The National Human Activity Pattern Survey (NHAPS): A resource for assessing exposure to environmental pollutants. J. Expo. Sci. Environ. Epidemiol. 2001;11:231–252. doi: 10.1038/sj.jea.7500165. [PubMed] [Cross Ref]
3. Kang-Wook K., Jihye K., Chang-Gun L., Junghee H. Accurate Indoor Location Tracking Exploiting Ultrasonic Reflections. IEEE Sens. J. 2016;16:9075–9088.
4. Liu H., Darabi H., Banerjee P., Liu J. Survey of wireless indoor positioning techniques and systems. IEEE Trans. Syst. Man Cybern. Part C Appl. Rev. 2007;37:1067–1080. doi: 10.1109/TSMCC.2007.905750. [Cross Ref]
5. Brena R.F., García-Vázquez J.P., Galván-Tejada C.E., Muñoz-Rodriguez D., Vargas-Rosales C., Fangmeyer J. Evolution of Indoor Positioning Technologies: A Survey. J. Sens. 2017 doi: 10.1155/2017/2630413. [Cross Ref]
6. Mainetti L., Patrono L., Sergi I. A survey on indoor positioning systems; Proceedings of the 2014 22nd International Conference on Software, Telecommunications and Computer Networks (SoftCOM); Split, Croatia. 17–19 September 2014; pp. 111–120.
7. Sayeef S., Madawala U.K., Handley P.G., Santoso D. Indoor personnel tracking using infrared beam scanning; Proceedings of the 2004 Position Location and Navigation Symposium; Monterey, CA, USA. 26–29 April 2004; pp. 698–705.
8. Xiaohan L., Makino H., Maeda Y. Basic study on indoor location estimation using visible light communication platform; Proceedings of the 2008 30th Annual International Conference of Engineering in Medicine and Biology Society; Vancouver, BC, Canada. 20–25 August 2008; pp. 2377–2380. [PubMed]
9. Ingram S.J., Harmer D., Quinlan M. Ultrawideband indoor positioning systems and their use in emergencies; Proceedings of the 2004 Position Location and Navigation Symposium; Monterey, CA, USA. 26–29 April 2004; pp. 706–715.
10. Mulloni A., Wagner D., Barakonyi I., Schmalstieg D. Indoor positioning and navigation with camera phones. IEEE Pervasive Comput. 2009;8 doi: 10.1109/MPRV.2009.30. [Cross Ref]
11. Sanpechuda T., Kovavisaruch L. A review of RFID localization: Applications and techniques; Proceedings of the 2008 Fifth International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology; Krabi, Thailand. 14–17 May 2008; pp. 769–772.
12. Rai A., Chintalapudi K.K., Padmanabhan V.N., Sen R. Zee: Zero-effort crowdsourcing for indoor localization; Proceedings of the 18th annual international conference on Mobile computing and networking; Istanbul, Turkey. 22–26 August 2012; pp. 293–304.
13. Goncalo G., Helena S. Indoor location system using ZigBee technology; Proceedings of the 2009 Third International Conference on Sensor Technologies and Applications; Athens, Greece. 18–23 June 2009; pp. 152–157.
14. Gualda D., Perez M.C., Urena J., Garcia J.C., Ruiz D., Garcia E., Lindo A. Ultrasonic LPS adaptation for smartphones; Proceedings of the 2013 International Conference on Indoor Positioning and Indoor Navigation; Montbéliard, France. 28–31 October 2013; pp. 1–8.
15. Gueuning F.E., Varlan M., Eugne C.E., Dupuis P. Accurate distance measurement by an autonomous ultrasonic system combining time-of-flight and phase-shift methods; Proceedings of the Instrumentation and Measurement Technology Conference; Brussels, Belgium. 4–6 June 1996; pp. 399–404.
16. Fritsch C., Ibanez A., Parrilla M. A digital envelope detection filter for real-time operation. IEEE Trans. Instrum. Meas. 1999;48:1287–1293. doi: 10.1109/19.816150. [Cross Ref]
17. Priyantha N.B., Chakraborty A., Balakrishnan H. The cricket location-support system; Proceedings of the 6th Annual International Conference on Mobile Computing and Networking; Boston, MA, USA. 6–11 August 2000; pp. 32–43.
18. Lu Z., Yang C., Qin D., Luo Y., Momayez M. Estimating ultrasonic time-of-flight through echo signal envelope and modified Gauss Newton method. Measurement. 2016;94:355–363. doi: 10.1016/j.measurement.2016.08.013. [Cross Ref]
19. Zhao Y., Fan X., Xu C.Z., Li X. ER-CRLB: An Extended Recursive Cramér–Rao Lower Bound Fundamental Analysis Method for Indoor Localization Systems. IEEE Trans. Veh. Technol. 2017;66:1605–1618. doi: 10.1109/TVT.2016.2553682. [Cross Ref]
20. Huang S.S., Huang C.F., Huang K.N., Young M.S. A high accuracy ultrasonic distance measurement system using binary frequency shift-keyed signal and phase detection. Rev. Sci. Instrum. 2002;73:3671–3677. doi: 10.1063/1.1502017. [Cross Ref]
21. Urena J., Hernandez A., Jimenez A., Villadangos J.M., Mazo M., Carcia J.C., Alvarez F.J., Marziani C.D., Perez M.C., Jiménez J.A., et al. Advanced sensorial system for an acoustic LPS. Microproces. Microsyst. 2007;31:393–401. doi: 10.1016/j.micpro.2007.02.003. [Cross Ref]
22. Seco F., Prieto J.C., Ruiz A.R.J., Guevara J. Compensation of multiple access interference effects in CDMA-based acoustic positioning systems. IEEE Trans. Instrum. Meas. 2014;63:2368–2378. doi: 10.1109/TIM.2014.2312511. [Cross Ref]
23. Álvarez F.J., Aguilera T., López-Valcarce R. CDMA-based acoustic local positioning system for portable devices with multipath cancellation. Dig. Signal Process. 2017;62:38–51. doi: 10.1016/j.dsp.2016.11.001. [Cross Ref]
24. Khyam M.O., Xinde L., Ge S.S., Pickering M.R. Multiple Access Chirp-Based Ultrasonic Positioning. IEEE Trans. Instrum. Meas. 2017;99:1–12. doi: 10.1109/TIM.2017.2737898. [Cross Ref]
25. Saad M.M., Bleakley C.J., Dobson S. Robust high-accuracy ultrasonic range measurement system. IEEE Trans. Instrum. Meas. 2011;60:3334–3341. doi: 10.1109/TIM.2011.2128950. [Cross Ref]
26. Gonzalez J.R., Bleakley C.J. High-precision robust broadband ultrasonic location and orientation estimation. IEEE J. Sel. Top. Signal Process. 2009;3:832–844. doi: 10.1109/JSTSP.2009.2027795. [Cross Ref]
27. Medina C., Segura J.C., De la Torre A. Ultrasound indoor positioning system based on a low-power wireless sensor network providing sub-centimeter accuracy. Sensors. 2013;13:3501–3526. doi: 10.3390/s130303501. [PMC free article] [PubMed] [Cross Ref]
28. Demirli R., Saniie J. Model-based estimation of ultrasonic echoes. Part I: Analysis and algorithms. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2001;48:787–802. doi: 10.1109/58.920713. [PubMed] [Cross Ref]
29. Elson J., Girod L., Estrin D. Fine-grained network time synchronization using reference broadcasts. ACM SIGOPS Oper. Syst. Rev. 2002;36:147–163. doi: 10.1145/844128.844143. [Cross Ref]
30. Cena G., Scanzio S., Valenzano A., Zunino C. Implementation and evaluation of the reference broadcast infrastructure synchronization protocol. IEEE Trans. Ind. Inform. 2015;11:801–811. doi: 10.1109/TII.2015.2396003. [Cross Ref]
31. Texas Instruments Incorporated CC2500. [(accessed on 2 October 2017)]; Available online: http://www.ti.com/lit/ds/symlink/cc2500.pdf.
32. Gueuning F., Varlan M., Eugene C., Dupuis P. Accurate distance measurement by an autonomous ultrasonic system combining time-of-flight and phase-shift methods. IEEE Trans. Instrum. Meas. 1997;46:1236–1240. doi: 10.1109/19.668260. [Cross Ref]
33. Khyam M.O., Ge S.S., Li X., Pickering M.R. Highly Accurate Time-of-Flight Measurement Technique Based on Phase-Correlation for Ultrasonic Ranging. IEEE Sens. J. 2017;17:434–443. doi: 10.1109/JSEN.2016.2631244. [Cross Ref]
34. ST Semiconductor Company STM32F722RC. [(accessed on 1 October 2017)]; Available online: http://www.st.com/zh/microcontrollers/stm32f722rc.html.
35. Texas Instruments Incorporated CC3200. [(accessed on 1 October 2017)]; Available online: http://www.ti.com/lit/ds/symlink/cc3220.pdf.
36. Luo H., Fang X., Ertas B. Hilbert transform and its engineering applications. AIAA J. 2009;47:923–927. doi: 10.2514/1.37649. [Cross Ref]
37. Ens A., Reindl L.M., Bordoy J., Wendeberg J., Schindelhauer C. Unsynchronized ultrasound system for TDOA localization; Proceedings of the 2014 International Conference on Indoor Positioning and Indoor Navigation (IPIN); Busan, Kore. 27–30 October 2014; pp. 601–610.
38. Angrisani L., Schiano Lo Moriello R. Estimating ultrasonic time-of-flight through quadrature demodulation. IEEE Trans. Instrum. Meas. 2006;55:54–62. doi: 10.1109/TIM.2005.861251. [Cross Ref]

Articles from Sensors (Basel, Switzerland) are provided here courtesy of Multidisciplinary Digital Publishing Institute (MDPI)