Next Article in Journal
INS/GNSS Integration for Aerobatic Flight Applications and Aircraft Motion Surveying
Next Article in Special Issue
An Improved BLE Indoor Localization with Kalman-Based Fusion: An Experimental Study
Previous Article in Journal
Support Vector Data Description Model to Map Specific Land Cover with Optimal Parameters Determined from a Window-Based Validation Set
Previous Article in Special Issue
Using Dimensionality Reduction Techniques for Refining Passive Indoor Positioning Systems Based on Radio Fingerprinting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Platform Architecture for Decentralized Positioning Systems

Institute for Computing in Civil Engineering & Geo Information Systems, Rheinisch-Westfälische Technische Hochschule Aachen University, Mies-van-der-Rohe-Str. 1, Aachen 52074, Germany
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(5), 957; https://doi.org/10.3390/s17050957
Submission received: 13 February 2017 / Revised: 7 April 2017 / Accepted: 19 April 2017 / Published: 26 April 2017

Abstract

:
A platform architecture for positioning systems is essential for the realization of a flexible localization system, which interacts with other systems and supports various positioning technologies and algorithms. The decentralized processing of a position enables pushing the application-level knowledge into a mobile station and avoids the communication with a central unit such as a server or a base station. In addition, the calculation of the position on low-cost and resource-constrained devices presents a challenge due to the limited computing, storage capacity, as well as power supply. Therefore, we propose a platform architecture that enables the design of a system with the reusability of the components, extensibility (e.g., with other positioning technologies) and interoperability. Furthermore, the position is computed on a low-cost device such as a microcontroller, which simultaneously performs additional tasks such as data collecting or preprocessing based on an operating system. The platform architecture is designed, implemented and evaluated on the basis of two positioning systems: a field strength system and a time of arrival-based positioning system.

1. Introduction

Nowadays, the accurate localization of a user or an object is indispensable for Location-Based Services (LBS), such as asset tracking, inventory management or routing and navigation. Location sensing systems have become popular and emerged as a vital research area in recent years; therefore, research and commercial products have been developed in academic and industrial contexts. A wide range of technologies is used for positioning; they can be classified into outdoor, indoor, infrastructure-less and infrastructure-based systems. The efficiencies of indoor and outdoor positioning systems differ greatly from each other, due to the fact that indoor surrounding areas raise a challenge for position finding, especially for systems based on wireless technologies, because of factors, such as: signal scattering and attenuation due to the high density of obstacles, multipath reflections from walls and furniture or Non-Line-of-Sight (NLoS) [1]. Infrastructure-less positioning systems do not require any infrastructure (e.g., special hardware in a building). Examples of infrastructure-less methods are video-based (e.g., smart-phone camera [2], robot navigation [3,4]), speech source, inertial sensor or wireless sensor network-based methods. In contrast, infrastructure-based positioning systems need infrastructure such as permanently-installed hardware, electrical power, walls or tripods for mounting the reference points or Internet access [1]. Examples of infrastructure-based positioning systems are ultra-wideband- (UWB) or magnetic field-based localization systems [5,6]. Another example of infrastructure-based positioning solutions is RFID-based indoor localization systems, which can include RFID readers and battery-free tags [7,8]. Generally, the RFID systems enable the localization by using the proximity-based, distance-based [9,10] or the scene analysis [7] methods.
The last two classes of positioning systems can complement each other. Since, on the one hand, the infrastructure-based positioning system can be used as a complementary system in an infrastructure-less positioning system in order to provide a starting point or to support a long-term stability, on the other hand, the infrastructure-less based method can cover areas that cannot be reached from an infrastructure-based method [11].
Another class of positioning systems is the collaborative and device-free localization systems. The collaborative position localization technique is based on various units (sensor nodes) that collaborate to achieve their positions. This technique is also named cooperative or network position location [12]. The collaboration between the sensor nodes has been normally assumed to occur in Wireless Sensor Networks (WSNs) [13]. Device-free positioning systems, which can be used for perimeter security, enable the tracking of users without wearing any devices. Device-free localization systems can use technologies such as pressure or electric field sensors installed under the floor [14].
Therefore, a platform is required for the realization of the previously-mentioned localization systems. The architecture of a platform includes data processing and the system component interaction and design. The data processing can be computed on a server or in a Mobile Station (MS) in the case of the centralized, or decentralized and distributed approach, respectively. The system component can be an MS or an anchor.
To the best of our knowledge, there is neither a standard for the design, nor a detailed description of the architecture of a positioning system. The architecture of most localization systems is roughly described and divided into two parts: the sensor hardware and the positioning algorithm [15]. The first part relies on a variety of technologies, such as electromagnetic waves (e.g., UWB or Wireless LAN (WLAN)), or ultrasound. The positioning algorithm is based on various signal measurement methods, such as Time of Arrival (ToA), signal strength or angle [15].
In this paper, we focus on anchor-based indoor positioning systems. The main contributions of this paper are:
  • The proposition of a layer- and modular-based architecture, which enables the calculation of the localization on-the-fly on the MS. The position is computed on low-cost and resource-constrained devices such as microcontrollers. The proposed platform is non-proprietary and open, since it can be easily extended with other sensors, in order to enable the implementation of positioning systems based on other technologies. The proposed platform can interoperate with existing systems and protocols such as 6LoWPAN.
  • We propose a preprocessing method to remove outliers in measured data, which is convenient for resource-constrained devices, especially in terms of limited stack memory.
  • We introduce a method that calculates the 3D position based on the Singular Value Decomposition (SVD).
  • We propose a pre-processed method for the localization calculation, avoiding the execution of memory and computationally-expensive algorithms such as the SVD or Moore–Penrose, on resource-constrained devices.
  • We also demonstrate the feasibility to deploy the Moore–Penrose algorithm, which is based on SVD, on resource-constrained devices.
  • In order to improve the position estimate, we use the Gauss–Newton, as well as the Levenberg– Marquardt algorithms, which are derived in convenient form for resource-constrained devices.
The proposed localization platform architecture is designed, implemented and tested for two localization systems, a field strength system and a time of arrival-based system. Finally, we evaluate the energy consumption of both systems.
The remainder of the paper is organized as follows: Firstly, we give an architectural overview for localization systems in Section 2. Then, we review related works in Section 3. We introduce in Section 4 a general layered platform architecture for localization systems. In Section 5, we present two examples of the platform architecture. We give an experimental evaluation of both systems in Section 6. Finally, we conclude our paper in Section 7.

2. Architectural Overview of Localization Systems

Positioning systems can be classified into three categories: centralized, decentralized and distributed architecture. Figure 1 illustrates different architectures.
The centralized architecture is the most commonly-used one, in which the MSs exclusively communicate with the base station; whereas, the MS modules deliver data collected from sensors to the base station, such as distance or signal measurement. In the decentralized and distributed architecture, both the data processing and the localization calculation take place in the MS, and no data are sent to the base station. Optionally, the result of a position finding could be sent to the base station if necessary. Figure 1b illustrates the case where the position is computed on the MSs and no communication link is used. Furthermore, in a distributed MS network (cooperative evaluation) [16], the MSs can exchange data, to collectively achieve and enhance the accuracy and the efficiency of a localization process [17]. In this case, distributed localization algorithms [18] can be used, which favor cooperative sensing and positioning by involving multiple MSs. The decentralized and the distributed offer the following advantages: the scalability, robustness and energy-awareness.

3. Related Work

Firstly, we compare various Indoor Localization Systems (ILSs) from industry, as well as from the research area in terms of the used platforms and architecture. Finally, we give a brief overview of various standards for localization systems.

3.1. Centralized Indoor Localization Systems

By the centralized architecture, the MSs uniquely communicate with the base station, which usually has more processing and storage capacity, as well as energy resources than the MSs.
Ubisense has developed a Real-Time Locating System (RTLS) based on the UWB radio technology, which allows a simultaneous tracking of many tags by combining the Angle of Arrival (AoA) and the Time Difference of Arrival (TDoA) methods [19]. The position of the tags is calculated in a location server using the signals transmitted from the tags to the reference points, which are synchronized by timing cable.The Ekahau RTLS is a fingerprinting location system based on the Received Signal Strength Indicator (RSSI), which uses an existing WLAN network in order to enable the tracking of phones, bar code scanners or people wearing WLAN tags [20]. De Angelis et al. proposed a centralized indoor positioning system that is based on alternating current (AC) magnetic fields. The system performs in two phases: the calibration and the trilateration phase [21,22].
The advantages of a centralized architecture include the usage of lightweight MSs, since the whole complexity is shifted to the base station. This enables a high position resolution. The disadvantages comprise single point of failure (e.g., the base station failure), as well as poor performance with a large number of MSs. Therefore, the system can get into energy starvation in the case of continuous communication between the MSs and the base station. Furthermore, the power-saving techniques are hard to implement.

3.2. Decentralized Indoor Localization Systems

Most decentralized localization systems are based on smartphones, which use the on-board sensors to calculate the position of a user. Smartphones are significant information interfaces between the user and the environment, which process substantial computational capacity and communication capability [23]. The smartphones incorporate relatively low cost sensors, which do not facilitate an accurate localization. Furthermore, potential energy consumers are the screen, continuous use of the on-board sensors and the communication interfaces. The smartphone-based localization systems can be classified into signal- and inertial-based mobile ILS [24].
The signal-based smartphone ILSs enable the location of users by using existing infrastructure such as a WLAN [25,26]. Most of them are based on fingerprinting, which needs a calibration phase by manually collecting a huge set of training data [27]. In addition, a retraining process is indispensable if the deployment environment is altered. They do not enable continuous (smooth) localization, as well as show limited accuracy due the instability and unreliability of the RSS and the absence of a causal relationship between the Euclidean distance and the RSSI [28].
The smartphone inertial-based ILSs enable the localization of a user by using the on-board magnetometer and accelerometer [29,30]. Inertial-based ILSs allow a continuous localization only for a certain time due to the drift error of the inertial sensors. In addition, compass fluctuations due to a variable magnetic field induced by the indoor environment (e.g., ferromagnetic building materials) lead to inaccurate heading estimation. A start point is essential for the tracking, which requires the use of an external positioning system or the intervention of the user [24]. The accuracy of the localization depends on the smartphone orientation and position with respect to the user’s body [24].

3.3. Distributed and Cooperative Localization Systems

Schmid et al. present a proof-of-concept of an ad hoc localization system for persons in a Wireless Sensor Network (WSN) [31]. Yamagushi et al. present two approaches for collaborative indoor localization, which are the stop-and-go localization and the People-Centric Navigation (PCN) methods [13].
Although cooperative ILSs enable the improvement of the location coverage, as well as the location accuracy, especially in the case of poor geometric conditions [12], these systems have to overcome operational requirements, as well as technical challenges [13]; whereby, the operational requirements can be the protection of privacy: users do not allow sharing their positions, or the incentive of users to share their position for, e.g., getting service from the network [32]. Technical challenges include the selection of the reference nodes, the energy efficiency, which requests the optimization of space and time of localization, and the self-localization, particularly by mobile nodes [13]. Another technical challenge is the self-organization due the large numbers of MSs with random environmental characteristics. The error propagation is a serious problem in the distributed evaluation. In addition, most of the multihop localization techniques by the cooperative localization in WSNs are not implemented and are only treated at the theory level; or they were tested in simulated environments [33].

3.4. Standards for Localization

Although there are several positioning systems from the commercial or research fields, the ANSI 371.1 RTLS and the IEEE 802.15.4.x localization standards give only the specification of the physical and media access control (MAC) layer. Figure 2 illustrates this relationship.
The ANSI 371.1 RTLS standard specifies the physical layer for an RTLS, which is called WhereNet and developed by Zebra Technology Company [34]. The WhereNet supports both indoor and outdoor real-time positioning by using the ToA method and the 2.4-GHz Direct Sequence Spread Spectrum (DSSS) PHY [35].
The IEEE 802.15.4a is the first international standard specifying the wireless physical layer, in order to enable precision ranging [36]. The physical layer is based on UWB technology and supports Impulse Radio (IR) UWB and Chirp Spread Spectrum (CSS) [37]. The MAC layer is specified in the IEEE 802.15.4-2003 [38].

4. General Architecture of a Platform for Positioning Systems

We use the decentralized architecture due to the on-the-fly capability, robustness, scalability and energy-awareness of the MS.

4.1. System Interaction and Components

The proposed anchor-based ILS is composed of reference stations (anchors) and an MS, which includes sensors enabling a distance or signal strength measurement. As illustrated in Figure 3, the MS performs a ranging or signal strength measurement of the reference stations, whereby the data, as well as the position computation processing occur on the MS. The localization process consists of three distinct phases: measurement, preprocessing and position estimation. In the first phase, the MS collects data after performing a measurement on the anchors. The measurement data are preprocessed at the second phase, for example to eliminate outliers. Finally, the position is calculated at the third phase.
The suggested platform follows a modular-based architecture that ensures the portability and the extensibility of the system. As illustrated in Figure 4, the system architecture of the MS is divided into two layers: the System Layer (SL) and the Application Layer (AL).

4.2. System Layer of the Mobile Station

The SL consists of two sublayers: the hardware and the Operating System (OS). Both sublayers will be discussed in the next Section 4.2.1 and Section 4.2.2, respectively.

4.2.1. Hardware Layer

The hardware sublayer includes a computing unit (e.g., Microcontroller Unit (MCU)), sensory unit and electronic drivers. The sensory unit involves sensors such as UWB or magnetometer sensors to accomplish distance or magnetic field strength measurement, respectively. The electronic driver circuits enable, for example, interfacing with external driver units.

4.2.2. Operating System Layer

The OS sublayer includes a software interface driver module, which provides interfaces to communicate with sensors or other devices. Most digital sensors support standardized interfaces such as Serial Peripheral Interface (SPI), Inter Integrated Circuit (I2C) buses and Universal Asynchronous Receiver Transceiver (UART), which encourage the interoperability and the extensibility of the system with other technologies. A synchronization unit can serve to synchronize the received sensor data between a transmitter and receiver.
OSs for resource-constrained devices will be presented in the next Section 4.2.3.

4.2.3. Operating System for Microcontroller

A multithreaded Operating System (OS) is needed, since various tasks are running simultaneously inside an MS. These tasks may run quasi-parallel [39], performing for example a measurement and a processing of unknown position. The OS represents the upper sub-layer of the system layer, which abstracts the hardware sub-layer (cf. Figure 4). An OS can be essentially characterized by the following key design issues:
i
kernel structure, which can follow a monolithic model, layered approach or microkernel paradigm,
ii
the scheduler and
iii
the programming model [40].
Resource-constrained devices such as microcontrollers are characterized by limited computation and storage capacity, as well as battery power. The OSs for resource-constrained devices vary in the architecture, real-time support, scheduling, as well as the programming model and language [41]. Examples for these operating systems are: FreeRTOS [42], TinyOS [43], Contiki [44] and RIOT-OS [40]. Table 1 shows a comparison of various open-source OSs for resource-constrained devices. Although the FreeRTOS supports the most microcontrollers, FreeRTOS does not enable low power management features, like most real-time OSs, since the energy-savings of modern MCUs are platform specific [45,46]. Contiki and TinyOS lack support for real time, as well as some developer-friendly features, such as standard multi-threading and standard C and C++ programming language.

4.3. Application Layer

In order to ensure the portability and the extensibility of the system with various applications and positioning algorithms, the AL as the highest level of the proposed platform follows a modular-based architecture. The interoperability of the AL, as well as of the MS can be extended by using an open-standard format, such as the JavaScript Object Notation (JSON) [47], in order to enable data exchange with other devices: for example, data exchange between the MS and applications located on a PC. The AL is subdivided into two sublayers: the preprocessing and position computing sublayers (cf. Figure 4). Both sublayers of the AL will be discussed in the next Section 4.3.1 and Section 4.3.2, respectively.

4.3.1. Preprocessing

The preprocessing sublayer comprises data filtering, which, for instance, reduces the effect of statistical outliers from data delivered from the SL. The outliers can be filtered out by using the mean, the median or the Median Absolute Deviation (MAD) filters [48,49,50]. The preprocessing sublayer can also provide a calibration routine to correct uncertainties in the measured data.

4.3.2. Position Estimation

The top level sublayer provides the algorithmic core that computes the position. Commonly, a positioning algorithm is applied to estimate the unknown position of an MS based on measurements to reference points. The positioning algorithms must satisfy some practical requirements in order to be implemented in a practical system: for example, the algorithm should be robust against noisy measurements, otherwise the performance of the algorithm can decrease drastically. Before we describe some possible algorithms, we briefly introduce several measurement methods along with performance metrics of positioning algorithms.
(a)
Measurement methods:
The common measurement methods are based on signal strength, angular or distance observations, such as Received Signal Strength (RSS), Angle of Arrival (AoA) or Time of Arrival (ToA). The ToA is the most popular measurement technique, which can be estimated by using various ranging techniques, such as the one-way or two-way ToA, and Time Difference of Arrival (TDoA) [51]. Additionally, hybrid measurements can be used for the positioning, such as TDoA/AoA [52] or ToA/RSS [53].
(b)
Performance metrics:
The performance of a positioning algorithm can be evaluated by the following metrics: accuracy, precision, complexity, robustness, scalability, resilience to error and noise, coverage and cost [15].
  • Accuracy metrics: The localization accuracy metric shows how well the ground truth and estimated positions match. There are a number of accuracy metrics such as the Root Mean Square Error (RMSE), the Cumulative Distribution Function (CDF), the Probability Density Function (PDF) or the Frobenius metric (FROB) [54].
  • Cost metrics: A practical evaluation criterion is the cost of an algorithm, which is often a trade-off against accuracy. Common cost metrics are: algorithm complexity, convergence time, power consumption, reference to node ratio [55].

4.3.3. Positioning Algorithms

Positioning algorithms can be classified into two groups: deterministic and probabilistic methods [12]. Deterministic methods determine directly the position based on the measurements by applying the lateration or the least squares method. The lateration is a popular location algorithm, which computes the position of an unknown MS by measuring its distance or angle from multiple reference positions. The algorithm is called trilateration if the number of the reference points is three, otherwise, it is called multilateration. On the other hand, probabilistic methods [56,57], which are known as Bayesian methods, enable the position finding by considering the uncertainty of the measurements [12]. The Bayesian methods proceed in two steps: the prediction and correction steps.

5. Examples of a Platform for Positioning Systems

Based on the general discussion of an open platform, we now present two exemplary platforms for positioning systems. The first positioning system is ToA based and makes use of the UWB technology, while the second is RSS based and uses Direct Current (DC)-pulsed magnetic signals. We apply the decentralized architecture for the two positioning systems. The UWB- and the magnetic signal-based positioning systems will be described in Section 5.1 and Section 5.2, respectively.

5.1. UWB-Based Positioning System

The proposed system is composed of various UWB transceivers with known positions as reference stations (RSs) and one MS, which also incorporates a UWB transceiver (see Figure 5). The UWB transceiver computes a peer-to-peer distance to the reference stations using the two-way ToA measurement technique [15]. The position is calculated on the MS based on the measured distances to the RSs.
The system, as well as the application layer of the UWB-based positioning system will be described in the next Section 5.1.1 and Section 5.1.2, respectively.

5.1.1. System Layer of the UWB-Based Mobile Station

As illustrated in Figure 6, the SL is based on the previously discussed template architecture (see Figure 4) and is composed of a hardware and an operating system sublayer, which will be described in this section.
(a)
Hardware:
The hardware comprises four subsystems: the power unit supplying the sensor board with energy, the Microcontroller Unit (MCU), the sensory unit and the driver circuits (see Figure 7). The hardware layer is implemented based on the STM32F4, ARM Cortex-M4 core operating at 168 MHz and the UWB module P440 ranging sensor from TIME DOMAIN®, which enables ranging measurements with an accuracy of a few centimeters. The properties of the used MCU and the UWB ranging sensor are summarized in the Table 2 and Table 3, respectively.
(b)
Operating system:
We use the RIOT-OS [61], which is an open-source IoT operating system developed at the “Freie Universität Berlin”. RIOT-OS is based on a microkernel architecture, which was deployed for a rescue scenario to track and monitor fire fighters. In order to fulfill severe real-time requirements for hard industrial or emergency scenarios, the micro-kernel provides a zero-latency interrupt handling and prioritized threads with a minimum context-switching time. RIOT-OS implements a tickless scheduler, in order to achieve a maximum energy savings and to support deep-sleep mode by all resource-constrained MCUs. Furthermore, the RIOT-OS supports the 6LoWPAN and the IPv6 protocols, which enables the interoperability with existing systems and protocols [62]. Based on the architecture of the RIOT-OS, we developed and integrated device drivers for the P440 ranging sensor and the UART controller. These software driver components build the driver module and are part of the system layer (cf. Figure 6).

5.1.2. Application Layer of the UWB-Based Mobile Station

As mentioned in Section 4.3.1, the preprocessing sublayer serves to remove the outliers in the data delivered from the SL. Therefore, we use the median filter to remove noise from the measured distances captured from the UWB sensor [63]. We apply the shell sort algorithm to implement the median filter, which does not require recursion, such as the quick sort algorithm [64]. Although the iterative shell sort algorithm is slower than the quick sort algorithm, it is suitable for resource-constrained devices such as MCUs with a limited stack size. The AL incorporates a command shell for the interaction with a user or an application by using the serial interface. Furthermore, the data are exchanged by using the JSON between the MS and other systems, such as applications located on a PC. This can enhance the interoperability and is achieved by using a minimalistic JSON parser at the MCU.
In this section, we firstly present an algebraic multilateration approach for the position estimation. Then, we describe two methods to calculate the linear square method for the ranging-based positioning system. Finally, we derive the equations to compute a non-linear least squares (NLS) method in a convenient form for MCUs.
(a)
Algebraic multilateration method:
Assume ( x , y , z ) and ( x i , y i , z i ) for i = 1 , 2 , . . . , n are the coordinates of the MS and of n reference points, respectively. In addition, the measured distances between the reference points and the MS are d i . The unknown location of the MS is the intersection of the spheres, whose equations are:
( x x i ) 2 + ( y y i ) 2 + ( z z i ) 2 = d i 2 i = 1 , . . . , n .
The system of nonlinear equations in (1) can be solved by different methods [48,65,66,67]. We solved it by transforming the system of equations into a matrix form [68]. The algorithm used is not related to a specific anchor, since most algorithms subtract the coordinates of a specific anchor for the linearization of the equation system. Additionally, the algorithm gives a measure of the solvability of the multilateration problem and provides a recursive least square approach to update the position [68]. The solution of the linearized system is completely determined if the distances from four reference points are known. The problem requires the estimation of the unknown position x = ( x , y , z ) such that:
A x = b ,
where the matrix A and the vector b have the following forms [68]:
A = 1 2 x 1 2 y 1 2 z 1 1 2 x 2 2 y 2 2 z 2 1 2 x 3 2 y 3 2 z 3 1 2 x n 2 y n 2 z n ,
b = d 1 2 x 1 2 y 1 2 z 1 2 d 2 2 x 2 2 y 2 2 z 2 2 d 3 2 x 3 2 y 3 2 z 3 2 d n 2 x n 2 y n 2 z n 2 .
Equation (2) can be solved based on the QR decomposition in the microcontroller. An efficient method to implement the QR decomposition is the Householder transformation [69]. If A is ill-conditioned or singular, the x position can be computed by using the Moore–Penrose pseudo-inverse algorithm [68,69]. The pseudo-inverse is the best linear reconstruction operator in the mean square sense [70], which is more robust and reliable than the Householder transformation, but substantially more computationally expensive:
x = A + b ,
whereby A + is the pseudo-inverse of the matrix A [68,69]. The pseudo-inverse matrix can be computed based on the SVD of the matrix A [71]. SVD enables the calculation of the underdetermined and overdetermined systems of linear equations. Furthermore, the SVD is more robust to numerical errors [72], but it is computationally expensive.
(b)
Preprocessed pseudo-inverse matrix:
Since the matrix A in (3) depends only on the coordinates of the RSs, the constant matrix ( A + ) can be computed externally, e.g., in a PC. This method enables saving the resources of the MS, which can be initialized with the matrix A + for, e.g., with the help of a serial communication. In this case, the computation of a new x position in Equation (2) is reduced to a matrix multiplication: A + b .
(c)
Non-linear least squares method: Gauss–Newton method:
The algebraic multilateration method does not always provide a good estimation due to the measurement uncertainties [68]. In this case, the NLS method can be used to improve the position calculated from the algebraic multilateration method. This method is based on the minimization of the squares of the errors:
F ( x , y , z ) = i = 1 n f i 2 ( x , y , z ) ,
whereby f i ( x , y , z ) is the error function:
f i ( x , y , z ) = ( x x i ) 2 + ( y y i ) 2 + ( z z i ) 2 d i i = 1 , . . . , n .
Minimizing the sum of the square errors is a common problem in the area of applied mathematics, which can be solved, for instance, with the Gauss–Newton or the Levenberg–Marquardt algorithms [73]. We use the Gauss–Newton method to improve the estimated position, which is calculated by using the pseudo-inverse algorithm.
Since the Gauss–Newton method requires the first derivatives, we define the following Jacobian matrix:
J f = f 1 x f 1 y f 1 z f 2 x f 2 y f 2 z f 3 x f 3 y f 3 z f n x f n y f n z .
We introduce the error function vector f :
f = f 1 , f 2 , f 3 , , f n T .
Starting by an initial position guess x ( 1 ) = ( x ˜ , y ˜ , z ˜ ) T calculated from the algebraic multilateration method, the Gauss–Newton method proceeds by the following iterations:
x ( k + 1 ) = x ( k ) + s ( k )
s ( k ) = ( J f T ( k ) J f ( k ) ) 1 J f T ( k ) f ( k ) ,
where x ( k ) is the k-th approximation of the position and s ( k ) is the k-th error correction vector. We calculate s by using the Moore–Penrose pseudo-inverse algorithm, since the QR-decomposition, such as the QR-Householder algorithm, can fail due to the singularity or bad conditioning of the matrices.
Using Equations (7) and (8) leads to:
J f T J f = ( i = 1 n ( x x i ) 2 ( f i + d i ) 2 i = 1 n ( x x i ) ( y y i ) ( f i + d i ) 2 i = 1 n ( x x i ) ( z z i ) ( f i + d i ) 2 i = 1 n ( x x i ) ( y y i ) ( f i + d i ) 2 i = 1 n ( y y i ) 2 ( f i + d i ) 2 i = 1 n ( y y i ) ( z z i ) ( f i + d i ) 2 i = 1 n ( x x i ) ( z z i ) ( f i + d i ) 2 i = 1 n ( y y i ) ( z z i ) ( f i + d i ) 2 i = 1 n ( z z i ) 2 ( f i + d i ) 2 ) ,
and:
J f T f = i = 1 n ( x x i ) f i ( f i + d i ) , i = 1 n ( y y i ) f i ( f i + d i ) , i = 1 n ( z z i ) f i ( f i + d i ) T .
Equations (12) and (13) are composed of sum terms, which can be implemented in the microcontroller, for example by the use of a for-loop; whereby, the upper bound of the loop is n, which is equal to the reference point’s number. The matrix in Equation (12) is symmetrical; this property can be used to reduce the computational burden by computing only the upper or lower part of the matrix. Finally, the terms ( f i + d i ) 2 and f i f i + d i , which appear in each sum term of the matrix elements, can be computed only once by each iteration in Equations (12) and (13), respectively.

5.2. Magnetic Field-Based Positioning System

The Magnetic Indoor Local Positioning System (MILPS) is a magnetic field based positioning system, which is a representative for an RSS-based localization system [74]. The decentralized MILPS is based on DC-pulsed magnetic signals that show no special multipath effects and have good characteristics for penetrating various obstacle [6]. The MILPS enables a decentralized control of the individual coils (reference stations), as well as the decentralized synchronization of the entire system without the need of communication technology. In our previous work [6], we showed the possibility to synchronize and to control the coils and MS based on low-cost real-time clocks (RTCs). Furthermore, we calculated the distances to the anchors on the MS in a two-dimensional scenario. Based on the previous work, we developed this exemplary platform for MILPS, whereas the MS and the RSs are equipped with real-time clocks, and the MS additionally incorporates a magnetic sensor (see Figure 8). The developed MILPS in this work enables the calculation of an optimized three-dimensional position on the MS based on the measurement of the magnetic field, as well as the elevation angle to the anchors.
The system and the application layer of the MILPS, which is illustrated in Figure 9, will be described in the next Section 5.2.1 and Section 5.2.2, respectively.

5.2.1. System Layer of MILPS

Such as the platform template in Figure 4, the MILPS system layer is composed of a hardware and an operating system sublayer, which will be described in this section.
(a)
Hardware:
The MILPS has a similar hardware architecture as the UWB-based localization system described before; therefore, the hardware of the MS comprises four subsystems: the power unit, the Microcontroller Unit (MCU), the sensory unit and the driver circuits. The hardware layer is implemented with the LPC2387 ARM7 core (see Table 4). The sensory unit of the MS includes the 3D-magnetic sensor HMR2300, which offers a range of ± 2 G [75]. The properties of the used HMR2300 magnetometer are summarized in Table 5. Furthermore, the hardware sublayer of the MS incorporates the DS3234 RTC in order to synchronize the gathered magnetic data, which are generated from the coils. The hardware of the MS is illustrated in Figure 10.
Each coil is driven via a Control Driver Unit (CDU), which includes a LPC2387-MCU, an RTC and a driver circuit. The driver circuit enables the CDU to interface with an H-bridge, in order to control the voltage polarity. Figure 11 illustrates the CDU, as well as the control of the coils.
(b)
Operating system:
Based on the architecture of the RIOT-OS, we developed and integrated the device drivers for the DS3234 real-time clock and the HMR2300 magnetometer. These software driver components are part of the SL (cf. Figure 9). A decentralized synchronization mechanism enables periodic control of the coils and the MS using the Time Division Multiple Access (TDMA) scheme. At the initialization phase, the RTCs of the MS and the coils are set to the same time. In the operating mode, the coils are activated in fixed duration slots, which are cyclically organized. Simultaneously, the MS gathers the magnetic data from the magnetometer, which can be assigned to the source coils based on the predefined time slots [6].

5.2.2. Application Layer of MILPS

The AL is the highest level of the MS, which follows a modular-based architecture. It is subdivided into two sublayers: the preprocessing and position computing layers (cf. Figure 7b). The first sublayer includes the data filtering module, which uses the the median filter to remove the outliers from the gathered magnetic data, which are delivered from the system layer. The top level sublayer represents the algorithmic core that computes the position on the MS. Both ALs incorporate a command shell for the interaction with a user or an application by using the serial interface, similar to the UWB-based system. Furthermore, the AL of the MS includes a minimal JSON parser. In the following section, we describe the algorithms for the position estimation.
Theoretically, the magnetic field B i generated from the coil i, is given by the following equation:
B i = K r i 3 1 + 3 sin 2 θ i i = 1 , 2 , . . . , n .
In this context, K = μ 0 N I F 4 π , where N describes the number of turns of the wire, I is the current running through the coil, F expresses the base area of the coil, μ 0 is the permeability of free space, r i is the distance between the MS and coil i and θ i is the MS elevation angle relative to the coil plane [74,77].
In the following subsections, we describe three methods to calculate the position of the MS: the first one is the algebraic multilateration method; the second and third methods are the NLS-based methods, which use the estimated position from the first method as a start value. We used the NLS Gauss–Newton, as well as the Levenberg–Marquardt algorithm, since the algebraic multilateration does not always deliver an optimized position and due to the nonlinear measurement model. Furthermore, we derive the equations to compute these NLS methods in a convenient form for MCUs.
(a)
Algebraic multilateration method:
In the two-dimensional case (2D), when the coil i and the magnetometer lay on the same horizontal plane, θ i is equal to zero. Thus, (14) is reduced to the following equation:
B i = K r i 3 i = 1 , 2 , . . . , n .
The position of the MS is computed by using the algebraic multilateration method [68], based on the distances to n coils, calculated according to:
r i = K B i 3 3 i = 1 , 2 , . . . , n .
In the general three-dimensional case (3D), the unknown elevation angles θ i can be estimated by using a three-axis accelerometer, which enables the measurement of the pitch angle β and the roll angle κ of the MS [74]. Based on the measured pitch and roll angles, the elevation angles between the MS and reference stations can be calculated as follows:
θ i = arctan 3 4 tan I i ± 3 4 tan I i 2 + 1 2 ,
whereby I i is the inclination of the magnetic field from coil i, which is calculated using:
I i = arcsin B x , i sin β + B y , i cos β sin κ + B z , i cos β cos κ B i ,
where B x , i , B y , i and B z , i are the magnetic field components in the coordinate system of the sensor, which is integrated in the MS, and B i = B x , i 2 + B y , i 2 + B z , i 2 is the magnetic field magnitude. Hence, the distances r i between the MS and the reference stations can be calculated based on the estimated elevation angles θ i :
r i = k 1 + 3 sin 2 θ i B i 3 i = 1 , 2 , . . . , n .
In this case, the position of the MS can also be calculated by using the algebraic multilateration algorithm. Therefore, similar to the UWB-based localization system, the pseudo-inverse matrix can be processed in a computing unit such as a PC or a laptop, in order to initialize the MCU with the preprocessed result. In this way, the MCU has only to compute a matrix multiplication, in order to estimate an MS’s position.
(b)
Gauss–Newton algorithm:
Based on (16), the ( x , y , z ) coordinates of the MS can be computed by solving the following nonlinear system of equations:
f i ( x , y , z ) = k r i 3 1 + 3 sin 2 θ i B i = 0 i = 1 , 2 , . . . , n
whereby, r i = ( x x i ) 2 + ( y y i ) 2 + ( z z i ) 2 , sin θ i = z z i r i ; ( x i , y i , z i ) and ( x , y , z ) are the coordinates of the i-th coil and the MS, respectively. B i is the measured magnetic field strength.
For simplicity, we set u = ( x , y , z ) ; the Gauss–Newton algorithm iteratively finds the best estimate u ^ , which minimizes the sum of squares:
u ^ = a r g min u i = 1 n ( f i ( u ) ) 2 .
The Gauss–Newton method starts with an initial guess x ( 1 ) calculated by the direct method and proceeds iteratively (see Equation (10)) [73,78]; where J f is the Jacobian matrix of the function f = ( f 1 , f 2 , . . . , f n ) at u k . The Jacobian matrix J f is calculated based on Equation (20) to:
J f = f 1 x f 1 y f 1 z f 2 x f 2 y f 2 z f i x f i y f i z f n x f n y f n z ,
whereby, f i x , f i y and f i z are respectively equal to:
f i x = 3 K Δ x i Δ x i 2 + Δ y i 2 + 5 Δ z i 2 d m i ,
f i y = 3 K Δ y i Δ x i 2 + Δ y i 2 + 5 Δ z i 2 d m i ,
and:
f i z = 12 K Δ z i 3 d m i ,
where Δ x i = x x i , Δ y i = y y i , Δ z i = z z i and d m i = Δ x i 2 + Δ y i 2 + 4 Δ z i 2 1 2 Δ x i 2 + Δ y i 2 + Δ z i 2 3 .
The iteration process stops when the updates become sufficiently small. Furthermore, the initial guess x ( 1 ) is calculated by using the algebraic multilateration method and the distances r i as described in the first method. Similar to the UWB-based ILS, the Gauss–Newton algorithm uses the Moore–Penrose to calculate the error correction in (10), as well as to choose the x ( k ) with the minimum error value.
If the elevation angle θ i is unknown, it can be approximated by replacing the term 1 + 3 sin 2 θ i by “1.5” in (19) [77] or by using an iterative algorithm [79].
  • Levenberg–Marquardt method:
    The Levenberg–Marquardt method is also an algorithm for solving the NLS problems that is based on the trust-region approach [73]. The advantage of the trust-region strategy is the stability against the rank-deficiency of the Jacobian matrix J f , which is one of the weaknesses of the Gauss–Newton method [73]. Like the Gauss–Newton method, the Levenberg–Marquardt method proceeds iteratively (see Equation 10), whereby the error correction vector is equal to:
    s ( k ) = ( J f T ( k ) J f ( k ) + μ 2 ( k ) I ) 1 J f T ( k ) f ( k ) ,
    where I is the identity matrix and μ is the damping parameter. The error vector s is calculated by using the QR-Householder instead of the Moore–Penrose pseudo-inverse algorithm, since it is less computing time consuming, and the Levenberg–Marquardt method is robust against the rank-deficiency of the Jacobian matrix J f .
    The initial damping-parameter μ ( 1 ) can be calculated based on the matrix A ( 1 )
    A ( 1 ) = J f T ( 1 ) J f ( 1 ) ,
    as follows:
    μ ( 1 ) = τ · m a x i { a i i ( 1 ) } ,
    where a i i ( 1 ) are the diagonal elements of the matrix A ( 1 ) and τ is chosen by the user. As a rule of thumb, a small value of τ should be chosen (e.g., τ = 10 6 ), if the initial guess x ( 1 ) is believed to be a good approximation; otherwise, τ = 10 3 or τ = 1 can be used. Furthermore, the value of μ ( k ) can be updated based on the gain ratio ϱ [80]:
    ϱ = G ( k ) G ( k + 1 ) G ( k ) g T ( k ) g ( k ) ,
    G ( k ) = f T ( k ) f ( k ) ,
    G ( k + 1 ) = f T ( k + 1 ) f ( k + 1 ) ,
    g ( k ) = f ( k ) + J f ( k ) s ( k ) .
    μ ( k + 1 ) = 2 · μ ( k ) if ϱ β 0 μ ( k ) 2 if ϱ β 1 ,
    where by 0 < β 0 < β 1 < 1 (e.g., β 0 = 0.2 , β 1 = 0.8 ).

6. Experimental Evaluation

In this section, we present the results of the experimental evaluation of the two presented positioning systems in Subsection 6.2 and Subsection 6.3, respectively. The aim of the this evaluation is to demonstrate the feasibility to implement the proposed platform; therefore, we will not address issues, such as the impact of the MS or the anchors’ placement and selection on the localization accuracy. These issues are treated in [12,81,82]. For the evaluation, we give a brief summary of the complexity of the used algorithms, as well as the results of the accuracy measurements of the UWB-based system and MILPS. Furthermore, we evaluated the computing time of the algorithms on the STM32F407 and the LPC2387 MCUs, which are running at 168 MHz and 72 MHz, respectively. Finally, we evaluated the energy consumption of the algorithms on both systems as well as of the UWB-based ILS and the MILPS in Section 6.4.

6.1. Complexity of the Algorithms Used

The algorithms used in this paper are based on the matrix multiplication, the Moore–Penrose and the QR-Householder algorithm, and their complicity is summarized in Table 6; whereby, m and n are the number of rows and columns of the matrix A , respectively, while, n and p are the number of rows and columns of the matrix B , respectively.

6.2. UWB-Based Localization System Evaluation

A static measurement setup was performed, in order to inspect the positioning performance of the UWB system, whereby four UWB reference transceivers were located in the corners of a 6 m × 7 m room and the UWB mobile stations are placed in several points and at three different heights in the room. Hence, the location of the MS is measured at 27 different locations, whereby the measurement is repeated fifty times at each location. The MSs lie in one meter grid (see Figure 12).

6.2.1. Accuracy Evaluation

Figure 13a illustrates the three-dimensional position errors by using the algebraic multilateration and the Gauss–Newton methods. Figure 13b shows the positioning error, which is defined as the Euclidean distance between the estimated and true position. Figure 13c contains the empirical CDF of the position error of all locations by using the algebraic multilateration method; the error in the x and y coordinates is less than 3.5 cm, while the error in the z coordinate is less than 25.3 cm. The error in the z coordinate results from the unfavorable geometrical configuration, since the reference stations are located approximately at the same height. By further applying the Gauss–Newton method, the positioning error is reduced to 2.2 cm in the x and y coordinate, as well as to 11.2 cm in the z coordinate (see Figure 13d).

6.2.2. Computing Time Measurement

The A + matrix is calculated based on the Moore–Penrose method, and only once, in the initial phase at the start of the MCU or the positioning application. Based on the computed A + matrix, the localization of the MS is determined by the algebraic multilateration method. The position of the MS can be improved by using the Gauss–Newton algorithm, which uses the position delivered from the algebraic multilateration method as a starting point. The Gauss–Newton proceeds iteratively up to the desired accuracy or the maximal iteration number is reached. The average iteration number in this experiment is five. Despite the constrained computing resources of the STM32F407 MCU, the computing time of each positioning is in the order of 0.032 ms without using the Gauss–Newton algorithm. In contrast, the mean estimated position time increases approximately up to 7.9 ms by using the Gauss–Newton algorithm. The evaluation of the described computing steps is summarized in Table 7.

6.3. MILPS Evaluation

For the computing time and the accuracy evaluation of the MILPS, we placed four coils inside two rectangular rooms. The two rooms are separated by a wall, which is 0.58 m thick (see Figure 14). Two coils are placed in a room of a surface of 5.3 × 6.7 m 2 , while the other coils are placed in a room of a surface of 5.79 × 6.7 m 2 . The MSs are placed in 27 various positions, since each MS represented in Figure 14 is located at three different heights of 0.655 m, 1.423 m and 2.308 m. The true positions of the reference and mobile stations are determined by geodetic methods with millimeter accuracy using a tachymeter. We choose this configuration to demonstrate that the MILPS can measure the position even if the coils and the MSs are separated by walls.

6.3.1. Accuracy Measurement

Figure 15 presents a comparison between the algebraic multilateration, Gauss–Newton and Levenberg–Marquardt methods in terms of errors, which are illustrated in Figure 15a. Figure 15b shows the positioning error, which is defined as the Euclidean distance between the estimated and true position. As illustrated by this figure, the position errors of Point Numbers 7 and 22 are out of the bound by the Gauss–Newton algorithm, since it diverges. Figure 15c shows the experimental results of the positioning error obtained from the algebraic multilateration method, which are represented by the CDF. In this example, the error in the x and y components is lower than 30 cm. However, the z component of the MS coordinates shows the worst performance, since three coils were placed at nearly equal heights.
The application of the Gauss–Newton method reduced the errors in the x and y coordinates to 10 cm, but impairs the z coordinate values (see Figure 15d). In contrast, by using the Levenberg–Marquardt method, the errors in the z coordinates are limited to 0.6 m for 93 % of the measured points and to 1.1 m for other points (see Figure 15e). Furthermore, Figure 15a shows that the Levenberg–Marquardt method has generally lower deviations in all coordinate components compared with other methods.

6.3.2. Computing Time Measurement

As explained in Section 5.2.2a, we used the algebraic multilateration method to estimate the position of the MSs illustrated in Figure 14. Therefore, the A + matrix is also calculated based on the Moore–Penrose method, only once, at the start of the MCU. Similar to the UWB-based system, the position estimation of the MSs can be optimized by using the Gauss–Newton or the Levenberg–Marquardt algorithms. Both algorithms utilize the position carried out by the algebraic multilateration as a starting point and proceed iteratively up to the desired accuracy or when the maximal iteration number is reached. The average iteration number in this experiment is seven and six for the Gauss–Newton and the Levenberg–Marquardt method, respectively. The mean time of calculating a starting position and the computing resource-constrained LPC2387-MCU is approximately 0.1 ms. The computing time of an estimated position increases approximately up to 33 ms or 21 ms by using the Gauss–Newton or the Levenberg–Marquardt algorithm, respectively. The evaluation of the described computing steps and the used NLS methods is summarized in Table 8.

6.4. Energy Consumption

The energy consumption of the algorithms is measured based on the measurement of the drain-source current in the supply line, which is powered by a reference voltage supply V c c ( V c c = 5 V). Hence, the energy used for each localization processing task can be calculated by integrating the electric power over the times, which are summarized in Table 7 and Table 8. We measured a current consumption of about 75 mA and 69 mA at the ambient temperature of 26 C for the STM32F4 and the LPC2387 MCU in the active mode, respectively. The measured energies for the localization algorithms by the UWB-based ILS and the MILPS are summarized in Table 9.
We also measured the energy consumption of the MS performed for the UWB-based ILS, as well as for the MILPS by using the aforementioned method for the energy consumption of the algorithms. We measured a current consumption of about 410 mA and 27 mA by the UWB transceiver and the magnetometer, respectively; whereas, the measurement times by the UWB-ILPS, as well as by the MILPS are 120 ms and 1 s, respectively. The total energies, which are required for a position estimation, are calculated based on the current drain of the sensors, as well as the energy consumption of the MCU (see Table 9). The energy usages of the UWB-based ILS and the MILPS are summarized in Table 10, whereby the energy consumptions of the UWB transceiver and the magnetometer are 246 mWs and 405 mWs, respectively. The UWB solution is more energy-efficient than MILPS, since the UWB sensor measures much faster than the magnetic sensor (approximately 9 x ).

7. Conclusions

In this article, we present a platform for indoor location systems that is designed and implemented for a decentralized architecture, as well as tested for two positioning technologies. The investigated systems use the UWB and magnetic technologies with different measurement methods, which are the time of arrival and the field strength techniques, respectively. The suggested platform is modular and layer based, in order to enable a better reusability of the software components and extensibility with various localization technologies and algorithms. The interoperation with other systems can be reached by using standardized interfaces and data format. Furthermore, the use of an IoT-capable OS enables the MS to be a part of the IoT network. The presented platform can be also a basis for a distributed localization approach.

Acknowledgments

The authors thank Naouar Guerchali for the support in the implementation and evaluation of the Levenberg–Marquardt algorithm in MATLAB.

Author Contributions

Zakaria Kasmi conceived the research, designed the platform architecture and implemented the hardware and software components. Furthermore, he integrated the software components in RIOT-OS, performed the experiments and data evaluation, and wrote all parts of the article. Abdelmoumen Norrdine offered valuable advices to migrate the positioning algorithms into RIOT-OS, he helped by the experiments and data evaluation and intensively reviewed the paper. Jörg Blankenbach gave suggestions to the experimental part and reviewed the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mautz, R. Indoor Positioning Technologies; ETH Zurich, Department of Civil, Environmental and Geomatic Engineering, Institute of Geodesy and Photogrammetry: Zurich, Switzerland, 2012. [Google Scholar]
  2. Werner, M.; Kessel, M.; Marouane, C. Indoor positioning using smartphone camera. In Proceedings of the 2011 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Centro Cultural Vila Flor, Guimaraes, Portugal, 21–23 September 2011; pp. 1–6. [Google Scholar]
  3. Zhao, Y.G.; Cheng, W.; Jia, L.; le Ma, S. The Obstacle Avoidance and Navigation Based on Stereo Vision for Mobile Robot. In Proceedings of the 2010 International Conference on Optoelectronics and Image Processing (ICOIP), Haiko, China, 11–12 November 2010; Volume 2, pp. 565–568. [Google Scholar]
  4. Desouza, G.N.; Kak, A.C. Vision for mobile robot navigation: A survey. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 24, 237–267. [Google Scholar] [CrossRef]
  5. Mahfouz, M.R.; Fathy, A.E.; Kuhn, M.J.; Wang, Y. Recent trends and advances in UWB positioning. In Proceedings of the IEEE MTT-S International Microwave Workshop on Wireless Sensing, Local Positioning (IMWS 2009) and RFID 2009, Cavtat, Croatia, 24–25 September 2009; pp. 1–4. [Google Scholar]
  6. Kasmi, Z.; Norrdine, A.; Blankenbach, J. Towards a Decentralized Magnetic Indoor Positioning System. Sensors 2015, 15, 29799. [Google Scholar] [CrossRef] [PubMed]
  7. Ni, L.M.; Liu, Y.; Lau, Y.C.; Patil, A.P. LANDMARC: Indoor location sensing using active RFID. In Proceedings of the First IEEE International Conference on Pervasive Computing and Communications (PerCom 2003), Fort Worth, TX, USA, 23–26 March 2003; pp. 407–415. [Google Scholar]
  8. Zhao, Y.; Smith, J.R. A battery-free RFID-based indoor acoustic localization platform. In Proceedings of the 2013 IEEE International Conference on RFID (RFID), Johor Bahru, Malaysia, 4–5 September 2013; pp. 110–117. [Google Scholar]
  9. Bouet, M.; dos Santos, A.L. RFID tags: Positioning principles and localization techniques. In Proceedings of the 2008 1st IFIP Wireless Days, Dubai, UAE, 24–27 November 2008; pp. 1–5. [Google Scholar]
  10. Zhang, Y.; Li, X.; Amin, M. Principles and Techniques of RFID Positioning. In RFID Systems; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2010; pp. 389–415. [Google Scholar]
  11. Hol, J.D.; Dijkstra, F.; Luinge, H.; Schon, T.B. Tightly coupled UWB/IMU pose estimation. In Proceedings of the 2009 IEEE International Conference on Ultra-Wideband, Vancouver, BC, Canada, 9–11 September 2009; pp. 688–692. [Google Scholar]
  12. Zekavat, R.; Buehrer, R.M. Handbook of Position Location: Theory, Practice and Advances, 1st ed.; Wiley-IEEE Press: Hoboken, NJ, USA, 2011. [Google Scholar]
  13. Yamaguchiy, H.; Higuchi, T.; Higashino, T. Collaborative Indoor Localization of Mobile Nodes. In Proceedings of the International Conference on Mobile Computing and Ubiquitous (ICMU2012), Okinawa, Japan, 23–25 May 2012; pp. 156–163. [Google Scholar]
  14. Kivimäki, T.; Vuorela, T.; Peltola, P.; Vanhala, J. A review on device-free passive indoor positioning methods. Int. J. Smart Home 2014, 8, 71–94. [Google Scholar] [CrossRef]
  15. Liu, H.; Darabi, H.; Banerjee, P.; Liu, J. Survey of Wireless Indoor Positioning Techniques and Systems. IEEE Trans. Syst. Man Cybern. C Appl. Rev. 2007, 37, 1067–1080. [Google Scholar] [CrossRef]
  16. Yao, J.; Balaei, A.; Hassan, M.; Alam, N.; Dempster, A. Improving Cooperative Positioning for Vehicular Networks. IEEE Trans. Veh. Technol. 2011, 60, 2810–2823. [Google Scholar] [CrossRef]
  17. Zhou, Z.; Cui, J.H.; Zhou, S. Localization for Large-Scale Underwater Sensor Networks. In Networking 2007. Ad Hoc and Sensor Networks, Wireless Networks, Next Generation Internet; Lecture Notes in Computer Science; Akyildiz, I., Sivakumar, R., Ekici, E., Oliveira, J., McNair, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; Volume 4479, pp. 108–119. [Google Scholar]
  18. Han, G.; Xu, H.; Duong, T.; Jiang, J.; Hara, T. Localization algorithms of Wireless Sensor Networks: A survey. Telecommun. Syst. 2013, 52, 2419–2436. [Google Scholar] [CrossRef]
  19. Ubisense Ltd. The Ubisense Precise Real-time Location System. 2016. Available online: http://www.ubisense.net (accessed on 22 April 2017).
  20. Ekahau Inc. The Ekahau Real Time Location System. 2016. Available online: https://www.airistaflow.com (accessed on 22 April 2017).
  21. De Angelis, G.; Pasku, V.; De Angelis, A.; Dionigi, M.; Mongiardo, M.; Moschitta, A.; Carbone, P. An Indoor AC Magnetic Positioning System. IEEE Trans. Instrum. Meas. 2015, 64, 1275–1283. [Google Scholar] [CrossRef]
  22. De Angelis, G.; De Angelis, A.; Pasku, V.; Moschitta, A.; Carbone, P. A hybrid outdoor/indoor Positioning System for IoT applications. In Proceedings of the 2015 IEEE International Symposium on Systems Engineering (ISSE), Rome, Italy, 29–30 September 2015; pp. 1–6. [Google Scholar]
  23. Yang, Z.; Wu, C.; Liu, Y. Locating in Fingerprint Space: Wireless Indoor Localization with Little Human Intervention. In Proceedings of the 18th Annual International Conference on Mobile Computing and Networking (Mobicom ’12), Istanbul, Turkey, 22–26 August 2012; ACM: New York, NY, USA, 2012; pp. 269–280. [Google Scholar]
  24. Galván-Tejada, C.E.; García-Vázquez, J.P.; Galván-Tejada, J.I.; Delgado-Contreras, J.R.; Brena, R.F. Infrastructure-Less Indoor Localization Using the Microphone, Magnetometer and Light Sensor of a Smartphone. Sensors 2015, 15, 20355. [Google Scholar] [CrossRef] [PubMed]
  25. Zhuang, P.; Wang, D.; Shang, Y. SMART: Simultaneous indoor localization and map construction using smartphones. In Proceedings of the The 2010 International Joint Conference on Neural Networks (IJCNN), Barcelona, Spain, 18–23 July 2010; pp. 1–8. [Google Scholar]
  26. Jiang, Y.; Pan, X.; Li, K.; Lv, Q.; Dick, R.P.; Hannigan, M.; Shang, L. ARIEL: Automatic Wi-fi Based Room Fingerprinting for Indoor Localization. In Proceedings of the 2012 ACM Conference on Ubiquitous Computing (UbiComp ’12), Pittsburgh, PA, USA, 5–8 September 2012; ACM: New York, NY, USA, 2012; pp. 441–450. [Google Scholar]
  27. Martin, E.; Vinyals, O.; Friedland, G.; Bajcsy, R. Precise Indoor Localization Using Smart Phones. In Proceedings of the 18th ACM International Conference on Multimedia (MM ’10), Firenze, Italy, 25–29 October 2010; ACM: New York, NY, USA, 2010; pp. 787–790. [Google Scholar]
  28. Wang, Y.; Zhao, B.; Jiang, Z. Rssi-based smooth localization for indoor environment. Sci. World J. 2014, 2014, 639142. [Google Scholar] [CrossRef] [PubMed]
  29. Park, K.; Shin, H.; Cha, H. Smartphone-based pedestrian tracking in indoor corridor environments. Pers. Ubiquitous Comput. 2013, 17, 359–370. [Google Scholar] [CrossRef]
  30. Kang, W.; Han, Y. SmartPDR: Smartphone-Based Pedestrian Dead Reckoning for Indoor Localization. IEEE Sens. J. 2015, 15, 2906–2916. [Google Scholar] [CrossRef]
  31. Schmid, J.; Völker, M.; Gädeke, T.; Weber, P.; Stork, W.; Müller-Glaser, K.D. An approach to infrastructure-independent person localization with an IEEE 802.15.4 WSN. In Proceedings of the 2010 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Zurich, Switzerland, 15–17 September 2010; pp. 1–9. [Google Scholar]
  32. Pourabdollah, A.; Meng, X.; Jackson, M. Towards low-cost collaborative mobile positioning. In Proceedings of the Ubiquitous Positioning Indoor Navigation and Location Based Service (UPINLBS), Kirkkonummi, Finland, 14–15 October 2010; pp. 1–5. [Google Scholar]
  33. Bal, M.; Liu, M.; Shen, W.; Ghenniwa, H. Localization in cooperative Wireless Sensor Networks: A review. In Proceedings of the 13th International Conference on Computer Supported Cooperative Work in Design (CSCWD 2009), Santiago, Chile, 22–24 April 2009; pp. 438–443. [Google Scholar]
  34. Mao, G.; Fidan, B.; Mao, G.; Fidan, B. Localization Algorithms and Strategies for Wireless Sensor Networks; Information Science Reference–Imprint of IGI Publishing: Hershey, PA, USA, 2009. [Google Scholar]
  35. Gu, Y.; Lo, A.; Niemegeers, I. A survey of indoor positioning systems for wireless personal networks. IEEE Commun. Surv. Tutor. 2009, 11, 13–32. [Google Scholar] [CrossRef]
  36. Karapistoli, E.; Pavlidou, F.N.; Gragopoulos, I.; Tsetsinas, I. An overview of the IEEE 802.15.4a Standard. IEEE Commun. Mag. 2010, 48, 47–53. [Google Scholar] [CrossRef]
  37. Sahinoglu, Z.; Gezici, S. Ranging in the IEEE 802.15.4a Standard. In Proceedings of the 2006 IEEE Annual Wireless and Microwave Technology Conference, Clearwater Beach, FL, USA, 4–5 December 2006; pp. 1–5. [Google Scholar]
  38. Salman, N.; Rasool, I.; Kemp, A.H. Overview of the IEEE 802.15.4 standards family for Low Rate Wireless Personal Area Networks. In Proceedings of the 2010 7th International Symposium on Wireless Communication Systems (ISWCS), York, UK, 19–22 September 2010; pp. 701–705. [Google Scholar]
  39. Tripathi, A.; Karnik, N. Trends in multiprocessor and distributed operating systems designs. J. Supercomput. 1995, 9, 23–49. [Google Scholar] [CrossRef]
  40. Baccelli, E.; Hahm, O.; Günes, M.; Wählisch, M.; Schmidt, T.C. RIOT OS: Towards an OS for the Internet of Things. In Proceedings of the 32nd IEEE International Conference on Computer Communications (IEEE INFOCOM), Turin, Italy, 14–19 April 2013. [Google Scholar]
  41. Farooq, M.O.; Kunz, T. Operating Systems for Wireless Sensor Networks: A Survey. Sensors 2011, 11, 5900–5930. [Google Scholar] [CrossRef] [PubMed]
  42. FreeRTOS Operating System Organisation. Available online: http://www.freertos.org (accessed on 22 April 2017).
  43. Levis, P.; Madden, S.; Polastre, J.; Szewczyk, R.; Whitehouse, K.; Woo, A.; Gay, D.; Hill, J.; Welsh, M.; Brewer, E.; et al. TinyOS: An operating system for sensor networks. In Ambient Intelligence; Springer: Berlin, Germany, 2004. [Google Scholar]
  44. Dunkels, A.; Gronvall, B.; Voigt, T. Contiki—A Lightweight and Flexible Operating System for Tiny Networked Sensors. In Proceedings of the First IEEE Workshop on Embedded Networked Sensors (Emnets-I), Tampa, FL, USA, November 2004. [Google Scholar]
  45. Will, H.; Schleiser, K.; Schiller, J. A real-time kernel for wireless sensor networks employed in rescue scenarios. In Proceedings of the IEEE 34th Conference on Local Computer Networks ( LCN 2009), Zürich, Switzerland, 20–23 October 2009; pp. 834–841. [Google Scholar]
  46. Simonović, M.; Saranovac, L. Power management implementation in FreeRTOS on LM3S3748. Serbian J. Electr. Eng. 2013, 10, 199–208. [Google Scholar] [CrossRef]
  47. Friesen, J. Java XML and JSON, 1st ed.; Apress: New York, NY, USA, 2016. [Google Scholar]
  48. Blankenbach, J.; Willert, V. Robuster räumlicher Bogenschnitt–Ein Ansatz zur robusten Positionsberechnung in Indoor-Szenarien. Allgemeine Vermessungs-Nachrichten (AVN) 2009, 9, 320–327. [Google Scholar]
  49. Leys, C.; Ley, C.; Klein, O.; Bernard, P.; Licata, L. Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. J. Exp. Soc.Psychol. 2013, 49, 764–766. [Google Scholar] [CrossRef]
  50. Rousseeuw, P.J.; Croux, C. Alternatives to the median absolute deviation. J. Am. Stat. Assoc. 1993, 88, 1273–1283. [Google Scholar] [CrossRef]
  51. Al-Ammar, M.A.; Alhadhrami, S.; Al-Salman, A.; Alarifi, A.; Al-Khalifa, H.S.; Alnafessah, A.; Alsaleh, M. Comparative Survey of Indoor Positioning Technologies, Techniques, and Algorithms. In Proceedings of the 2014 International Conference on Cyberworlds (CW), Santander, Spain, 6–8 October 2014; pp. 245–252. [Google Scholar]
  52. Cong, L.; Zhuang, W. Hybrid TDOA/AOA mobile user location for wideband CDMA cellular systems. IEEE Trans. Wirel. Commun. 2002, 1, 439–447. [Google Scholar] [CrossRef]
  53. Huang, J.Y.; Wan, Q. Comments on “The Cramer-Rao Bounds of Hybrid TOA/RSS and TDOA/RSS Location Estimation Schemes”. IEEE Commun. Lett. 2007, 11, 848–849. [Google Scholar] [CrossRef]
  54. Efrat, A.; Forrester, D.; Iyer, A.; Kobourov, S.G.; Erten, C.; Kilic, O. Force-directed Approaches to Sensor Localization. ACM Trans. Sen. Netw. 2010, 7, 27. [Google Scholar] [CrossRef]
  55. Yedavalli, K.; Krishnamachari, B. Sequence-Based Localization in Wireless Sensor Networks. IEEE Trans. Mob. Comput. 2008, 7, 81–94. [Google Scholar] [CrossRef]
  56. Baggio, A.; Langendoen, K. Monte Carlo localization for mobile wireless sensor networks. Ad Hoc Netw. 2008, 6, 718–733. [Google Scholar] [CrossRef]
  57. Peng, R.; Sichitiu, M.L. probabilistic, constraint-based localization for wireless sensor networks. In Proceedings of the Second Annual IEEE Sensor and Ad Hoc Communications and Networks, Santa Clara, CA, USA, 26–29 September 2005; pp. 541–550. [Google Scholar]
  58. STMicroelectronics. Discovery Kit with STM32F407VG MCU, Rev. 5; STMicroelectronics: Geneva, Switzerland, 2016.
  59. Time Domain®. Data Sheet/User Guide PulsON® 440, Rev. 1; Time Domain: Huntsville, AL, USA, 2016.
  60. STMicroelectronics. STM32F405xx STM32F407xx data sheet, Rev. 7; STMicroelectronics: Geneva, Switzerland, 2016.
  61. Hahm, O.; Baccelli, E.; Petersen, H.; Wählisch, M.; Schmidt, T.C. Demonstration Abstract: Simply RIOT: Teaching and Experimental Research in the Internet of Things. In Proceedings of the 13th International Symposium on Information Processing in Sensor Networks (IPSN ’14), Berlin, Germany, 15–17 April 2014; IEEE Press: Piscataway, NJ, USA, 2014; pp. 329–330. [Google Scholar]
  62. Vasseur, J.P.; Dunkels, A. Interconnecting Smart Objects with IP: The Next Internet; Morgan Kaufmann Publishers Inc.: San Francisco, CA, USA, 2010. [Google Scholar]
  63. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes in C: The Art of Scientific Computing, 2nd ed.; Cambridge University Press: New York, NY, USA, 1992. [Google Scholar]
  64. Gupta, D.P. Design and Analysis of Algorithms; Phi Learning: Delhi, India, 2012. [Google Scholar]
  65. Pradhan, S.; Hwang, S.S. Mathematical analysis of line intersection algorithm for TOA trilateration method. In Proceedings of the 2014 Joint 7th International Conference on and Advanced Intelligent Systems (ISIS), 15th International Symposium on Soft Computing and Intelligent Systems (SCIS), Kita-Kyushu, Japan, 3–6 December 2014; pp. 1219–1223. [Google Scholar]
  66. Silva, B.; Pang, Z.; Åkerberg, J.; Neander, J.; Hancke, G. Experimental study of UWB-based high precision localization for industrial applications. In Proceedings of the 2014 IEEE International Conference on Ultra-WideBand (ICUWB), Paris, France, 1–3 September 2014; pp. 280–285. [Google Scholar]
  67. Beutel, J. Geolocation in a picoradio environment. Diploma Thesis, ETH Zurich, Electronics Laboratory, Zurich, Switzerland, 1999. [Google Scholar]
  68. Norrdine, A. An Algebraic Solution to the Multilateration Problem. In Proceedings of the 15th International Conference on Indoor Positioning and Indoor Navigation, Sydney, Australia, 13–15 November 2012; Volume 13. [Google Scholar]
  69. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes 3rd Edition: The Art of Scientific Computing, 3rd ed.; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
  70. Daubechies, I. Ten Lectures on Wavelets; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1992. [Google Scholar]
  71. Noble, B.; Daniel, J. Applied Linear Algebra; Prentice-Hall: Upper Saddle River, NJ, USA, 1988. [Google Scholar]
  72. Leach, S. Singular Value Decomposition—A Primer; Department of Computer Science, Brown University: Providence, RI, USA, 1995; Unpublished Manuscript. [Google Scholar]
  73. Nocedal, J.; Wright, S. Numerical Optimization; Springer Series in Operations Research and Financial Engineering; Springer: New York, NY, USA, 2000. [Google Scholar]
  74. Norrdine, A. Präzise Positionierung und Orientierung innerhalb von Gebäuden; Technische Univ. Darmstadt, Geodätisches Institut: Darmstadt, Germany, 2009. [Google Scholar]
  75. Honeywell International Inc. Smart Digital Magnetometer HMR2300, Rev. J.; Honeywell International Inc.: Morris Plains, NJ, USA, 2012.
  76. NXP Semiconductors. LPC2387 Product Data Sheet, Rev. 5.1; NXP Semiconductors: Eindhoven, The Netherlands, 2013.
  77. Prigge, E. A Positioning System with No Line-of-Sight Restrictions for Cluttered Environments; Stanford University: Stanford, CA, USA, 2004. [Google Scholar]
  78. Björck, A. Numerical Methods for Least Squares Problems; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1996. [Google Scholar]
  79. Blankenbach, J.; Norrdine, A. Magnetic Indoor Local Positioning System. In Indoor Wayfinding and Navigation; CRC Press: Boca Raton, FL, USA, 2015; pp. 53–80. [Google Scholar]
  80. Dahmen, W.; Reusken, A. Numerik für Ingenieure und Naturwissenschaftler; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  81. Chen, Y.; Kobayashi, H. Signal strength based indoor geolocation. In Proceedings of the IEEE International Conference on Communications (ICC 2002), New York, NY, USA, 28 April–2 May 2002; Volume 1, pp. 436–439. [Google Scholar]
  82. Krishnakumar, A.S.; Krishnan, P. On the accuracy of signal strength-based estimation techniques. In Proceedings of the IEEE 24th Annual Joint Conference of the IEEE Computer and Communications Societies, Miami, FL, USA, 13–17 March 2005; Volume 1, pp. 642–650. [Google Scholar]
  83. Neapolitan, R. Foundations of Algorithms; Jones & Bartlett Learning: Burlington, MA, USA, 2015. [Google Scholar]
  84. Martinez, D. High Performance Embedded Computing Handbook: A Systems Perspective; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  85. Smoktunowicz, A.; Wróbel, I. Numerical aspects of computing the Moore–Penrose inverse of full column rank matrices. BIT Numer. Math. 2012, 52, 503–524. [Google Scholar] [CrossRef]
Figure 1. Position evaluation: (a) central, (b) decentral and (c) distributed.
Figure 1. Position evaluation: (a) central, (b) decentral and (c) distributed.
Sensors 17 00957 g001
Figure 2. Localization system architectures. (a) Industrial and academic area, (b) ANSI371.1 RTLS and IEEE 802.15.4.x.
Figure 2. Localization system architectures. (a) Industrial and academic area, (b) ANSI371.1 RTLS and IEEE 802.15.4.x.
Sensors 17 00957 g002
Figure 3. Components’ interaction.
Figure 3. Components’ interaction.
Sensors 17 00957 g003
Figure 4. The architecture of the mobile station.
Figure 4. The architecture of the mobile station.
Sensors 17 00957 g004
Figure 5. UWB-based Indoor Localization System (ILS) principle. M S , Mobile Station; R S i , Reference Station i.
Figure 5. UWB-based Indoor Localization System (ILS) principle. M S , Mobile Station; R S i , Reference Station i.
Sensors 17 00957 g005
Figure 6. System architecture of a time of arrival-based MS.
Figure 6. System architecture of a time of arrival-based MS.
Sensors 17 00957 g006
Figure 7. UWB-based localization system hardware. (a) STM32F407 Discovery Board [58], (b) UWB P440 transceiver [59], (c) simplified hardware block diagram.
Figure 7. UWB-based localization system hardware. (a) STM32F407 Discovery Board [58], (b) UWB P440 transceiver [59], (c) simplified hardware block diagram.
Sensors 17 00957 g007
Figure 8. Magnetic indoor local positioning system (MILPS) principle.
Figure 8. Magnetic indoor local positioning system (MILPS) principle.
Sensors 17 00957 g008
Figure 9. MILPS: System architecture of the Mobile Station (MS).
Figure 9. MILPS: System architecture of the Mobile Station (MS).
Sensors 17 00957 g009
Figure 10. MILPS: Mobile station hardware. (a) Mobile station hardware, (b) simplified hardware block diagram.
Figure 10. MILPS: Mobile station hardware. (a) Mobile station hardware, (b) simplified hardware block diagram.
Sensors 17 00957 g010
Figure 11. MILPS reference stations hardware. (a) Control Driver Unit (CDU), (b) simplified hardware block diagram.
Figure 11. MILPS reference stations hardware. (a) Control Driver Unit (CDU), (b) simplified hardware block diagram.
Sensors 17 00957 g011
Figure 12. Experimental setup for position measures between various MSs and four reference stations.
Figure 12. Experimental setup for position measures between various MSs and four reference stations.
Sensors 17 00957 g012
Figure 13. UWB-based system: scatter plots, position errors, and empirical CDFs. ML, Multilateration; GNM, Gauss–Newton Method. (a) Scatter plot of the ML and the GNM, (b) position error of the ML and the GNM, (c) CDF of the points estimated by the ML, (d) CDF of the estimated positions after GNM.
Figure 13. UWB-based system: scatter plots, position errors, and empirical CDFs. ML, Multilateration; GNM, Gauss–Newton Method. (a) Scatter plot of the ML and the GNM, (b) position error of the ML and the GNM, (c) CDF of the points estimated by the ML, (d) CDF of the estimated positions after GNM.
Sensors 17 00957 g013
Figure 14. Experimental setup for position measures between various MSs and four coils.
Figure 14. Experimental setup for position measures between various MSs and four coils.
Sensors 17 00957 g014
Figure 15. MILPS: scatter plots, position errors, and empirical CDFs. GNM, Gauss–Newton Method; ML, Multilateration; LVM, Levenberg–Marquardt Method. (a) Scatter plot of the ML, the GNM and the LVM, (b) position error of the ML, the GNM and the LVM, (c) CDF of the points estimated by the multilateration, (d) CDF of the estimated positions after GNM, (e) CDF of the estimated positions after LVM.
Figure 15. MILPS: scatter plots, position errors, and empirical CDFs. GNM, Gauss–Newton Method; ML, Multilateration; LVM, Levenberg–Marquardt Method. (a) Scatter plot of the ML, the GNM and the LVM, (b) position error of the ML, the GNM and the LVM, (c) CDF of the points estimated by the multilateration, (d) CDF of the estimated positions after GNM, (e) CDF of the estimated positions after LVM.
Sensors 17 00957 g015
Table 1. Comparison of operating systems. FIFO, First-In First-Out; TOS, TinyOS.
Table 1. Comparison of operating systems. FIFO, First-In First-Out; TOS, TinyOS.
OSArchitectureReal TimeSchedulingProgramming ModelProgramming Language
FreeRTOSMonolithicFull SupportRound-robin preemptive and cooperativeThreadsC
TinyOSMonolithicNo SupportFIFOPrimarily event driven, support for TOS threadsnesC
ContikiModularPartial SupportEvent basedProtothreads and eventsC with some constraints
RIOT-OSMicrokernelFull SupportTickless, preemptive scheduling with prioritiesThreadsC and C++
Table 2. The properties of the STM32F407 MCU. MCU, Microcontroller Unit [60].
Table 2. The properties of the STM32F407 MCU. MCU, Microcontroller Unit [60].
MCUFamilyVendorFrequencyRAMFlash
STM32F407ARM Cortex-M4ST Microelectronics168 MHz192 KB1024 KB
Table 3. The properties of the deployed P440 ranging sensor [59].
Table 3. The properties of the deployed P440 ranging sensor [59].
AccuracyMax. Operating RangeMax. Ranging RateFrequency RangeTransmission Power
2.1 cm300 m–1100 m125 Hz 3.1 GHz– 4.8 GHz50 μ W
Table 4. The properties of the deployed LPC2387 MCU. [76].
Table 4. The properties of the deployed LPC2387 MCU. [76].
MCUFamilyVendorFrequencyRAMFlash
LPC2387ARM7NXP72 MHz96 KB512 KB
Table 5. The properties of the HMR2300 magnetometer. FS, Full Scale [75].
Table 5. The properties of the HMR2300 magnetometer. FS, Full Scale [75].
RangeSample RateResolutionAccuracy
± 2 Gauss (G)up to 154 Hzup to 70 μ G 0.5 % FS (over ± 1 G)
Table 6. Complexity of various algorithms.
Table 6. Complexity of various algorithms.
AlgorithmComplexity (Flops)
Matrix multiplication: A m , n × B n , p m p ( 2 n 1 ) [83]
QR-Householder 2 m n 2 2 3 n 3 [84]
Moore–Penrose pseudoinverse 14 m n 2 + 16 3 n 3 [85]
Table 7. Mean computing times of the algorithms used by the UWB-based System. Computing times measured on an STM32F407 running at 168 MHz.
Table 7. Mean computing times of the algorithms used by the UWB-based System. Computing times measured on an STM32F407 running at 168 MHz.
AlgorithmComputing Time [μs]
A + for the multilateration method (at the start)2115
Multilateration method32
Gauss–Newton per iteration1561
Table 8. Mean computing times of the used algorithms by MILPS. Computing times measured on an LPC2387 running at 72 MHz.
Table 8. Mean computing times of the used algorithms by MILPS. Computing times measured on an LPC2387 running at 72 MHz.
AlgorithmComputing Time (μs)
A + for the multilateration method (at the start)4563
Multilateration method92
Gauss–Newton method per iteration4645
Levenberg–Marquardt method per iteration3467
Table 9. List of measured energy consumption values of the algorithms by the UWB-based ILS and MILPS.
Table 9. List of measured energy consumption values of the algorithms by the UWB-based ILS and MILPS.
AlgorithmUWB-Based ILS Energy (μWs)MILPS Energy (μWs)
A + by the multilateration method (at the start) 793.13 1574.24
Multilateration method12 31.74
Gauss–Newton per iteration 585.38 1602.53
Levenberg–Marquardt method per iteration 1196.12
Table 10. List of measured energy consumption of the UWB-based ILS and the MILPS for a position estimation.
Table 10. List of measured energy consumption of the UWB-based ILS and the MILPS for a position estimation.
Localization SystemEnergy (mWs)
UWB-based ILS 246 + 0.793 + 0.012 + 5 × 0.585 249.73
MILPS (Gauss–Newton) 405 + 1.6 + 0.032 + 7 × 1.63 418.04
MILPS (Levenberg–Marquardt) 405 + 1.6 + 0.032 + 6 × 1.21 413.89

Share and Cite

MDPI and ACS Style

Kasmi, Z.; Norrdine, A.; Blankenbach, J. Platform Architecture for Decentralized Positioning Systems. Sensors 2017, 17, 957. https://doi.org/10.3390/s17050957

AMA Style

Kasmi Z, Norrdine A, Blankenbach J. Platform Architecture for Decentralized Positioning Systems. Sensors. 2017; 17(5):957. https://doi.org/10.3390/s17050957

Chicago/Turabian Style

Kasmi, Zakaria, Abdelmoumen Norrdine, and Jörg Blankenbach. 2017. "Platform Architecture for Decentralized Positioning Systems" Sensors 17, no. 5: 957. https://doi.org/10.3390/s17050957

APA Style

Kasmi, Z., Norrdine, A., & Blankenbach, J. (2017). Platform Architecture for Decentralized Positioning Systems. Sensors, 17(5), 957. https://doi.org/10.3390/s17050957

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop