Next Article in Journal
Person Re-Identification with RGB–D and RGB–IR Sensors: A Comprehensive Survey
Next Article in Special Issue
A Novel Swarm Intelligence Algorithm with a Parasitism-Relation-Based Structure for Mobile Robot Path Planning
Previous Article in Journal
Biometric Recognition: A Systematic Review on Electrocardiogram Data Acquisition Methods
Previous Article in Special Issue
Research on Smart Tourism Oriented Sensor Network Construction and Information Service Mode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Indoor Positioning Design for Mobile Phones via Integrating a Single Microphone Sensor and an H2 Estimator

1
Department of Mechanical Engineering, National Pingtung University of Science and Technology, Pingtung 912301, Taiwan
2
National Chung-Shan Institute of Science and Technology, Taoyuan 32546, Taiwan
3
Department of Systems and Naval Mechatronic Engineering, National Cheng Kung University, Tainan 701401, Taiwan
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(3), 1508; https://doi.org/10.3390/s23031508
Submission received: 28 December 2022 / Revised: 26 January 2023 / Accepted: 27 January 2023 / Published: 29 January 2023
(This article belongs to the Special Issue Smart Mobile and Sensing Applications)

Abstract

:
An indoor positioning design developed for mobile phones by integrating a single microphone sensor, an H2 estimator, and tagged sound sources, all with distinct frequencies, is proposed in this investigation. From existing practical experiments, the results summarize a key point for achieving a satisfactory indoor positioning: The estimation accuracy of the instantaneous sound pressure level (SPL) that is inevitably affected by random variations of environmental corruptions dominates the indoor positioning performance. Following this guideline, the proposed H2 estimation design, accompanied by a sound pressure level model, is developed for effectively mitigating the influences of received signal strength (RSS) variations caused by reverberation, reflection, refraction, etc. From the simulation results and practical tests, the proposed design delivers a highly promising indoor positioning performance: an average positioning RMS error of 0.75 m can be obtained, even under the effects of heavy environmental corruptions.

1. Introduction

Indoor positioning designs for mobile communication systems have attracted increasing attention recently due to emergency and security concerns. With the increase in time spent by mobile phone users in indoor environments, a mandatory bill, called Enhanced 911 (E911) was passed by the USA, providing the specific requirements for the indoor positioning accuracy of mobile communication systems. The indoor positioning accuracy requirements of E911 are less than 100 m (67% calls) and 300 m (95% calls) [1]. However, the well-known outdoor positioning system, GPS, is not useful for guaranteeing the accuracy of indoor positioning because of the shielding effect of the positioning signal transmission of satellites, i.e., the positioning ability of GPS is weak in regards to indoor environments such as shopping malls, hypermarkets, office building, etc. A variety of theoretical and available technologies have been proposed in recent decades for achieving the requirements of the indoor positioning [2,3,4,5], based on methods using infrared rays (IR), ultrasound, radio-frequency identification (RFID), wireless local area networks (WLAN), Bluetooth, audible sound, and other technologies. Nevertheless, not all of the above-mentioned methods can be applied to mobile communication systems due to complexity of implementation and the costs of the hardware and software; hence, a new design which is suitable for performing indoor positioning, based on the integration of a single microphone of the mobile communication system and one positioning estimation design. is proposed in this investigation.
Many indoor positioning methods have been proposed in recent years, and an accurate positioning estimation is the common goal of these designs. The earliest design was based on an infrared ray (IR) system [5,6,7,8,9], and until now, it was the most simple and common design for the purpose of indoor positioning. It is possible that the indoor positioning performance of these types of IR-based designs is acceptable in indoor environments which are well-arranged. However, practically, environmental light sources, such as florescent light and sunlight, are always tricky problems which strongly reduce position accuracy [8]. To compensate for the effect of this environment disturbance, several filters have been developed [8,10]. Applications of IR-based indoor positioning designs are constrained due to these interferences, and constructing such a light-based system comes at a high cost. Radio frequency (RF)-based technologies [11,12] are the other designs used for achieving the indoor positioning goal. Categories of RF-based positioning methods are mainly divided into two groups: 1. the radio frequency identification (RFID) method, and 2. the Wireless Local Area Network (WLAN) method. WhereNet is a real-time indoor and outdoor positioning design developed by Zebra Technologies via the RFID method for users who are in intricate environments, such as libraries, offices, etc. [4,13]. Key parts comprising this positioning system include tags, positioning antennas, processors, servers, where ports, and a software algorithm using the differential time of arrival method (TDOA) for calculating the locations of moving tags.
The major disadvantage of this kind of positioning method is that is necessitates an enormous cost for building numerous infrastructures in the working area. As to the WLAN positioning design, most of WLAN-based algorithms are developed via adopting the received signal strength of the WLAN signals. Generally speaking, WLAN-based designs possess a low-cost feature due to the popularity of WLAN infrastructures in indoor environments. By building on this advantage, a RADAR positioning system has been developed by a Microsoft research group, combining received signal strength detection with the triangulation positioning method. Unfortunately, the received signal strength of WLAN is naturally and inevitably affected by various environmental uncertainties, such as the multipath effect, the no line-of-sight effect [14,15,16,17], etc.; hence, some auxiliary designs are proposed to reduce the effect of environmental uncertainties [15,18,19], e.g., a radio map using the fingerprinting method, or the approximation design of a specific environment using the fuzzy logic approach. Fingerprinting techniques work well when the stored information of WLAN access points (APs) increase significantly, i.e., more APs are needed [20]. Algorithms based on sound detection provide another potential method for indoor positioning designs. The Active Bat system was developed by AT&T Cambridge by mimicking the navigation behavior of bats [21].
For improving the accuracy of the Active Bat system, a new design called the Cricket system combines the overall design of the Active Bat system with an extra RF method [22,23]. For the above developments, multi-sensors and ultrasound designs are adopted. Daredevil, developed by Microsoft, uses audible sound sources, providing an indoor positioning ability for mobile phones with at least two embedded microphones [4]. In the work in [14], the mobile phone is used as an emitter, and it collocates with the Wi-Fi network to achieve higher accuracy. A contrasting design, in which the handheld devices are arranged as receivers for some predefined emitters, is another popular design because mobile phone users always need the real-time display of the mobile phones’ monitors [24,25]. The positioning algorithm based on TOA or TDOA methods requires that multi-sensors be used, i.e., the total costs will be higher than those using the single sensor design; besides, the TOA or TDOA methods are degraded by four main factors: 1. background noise, 2. the multi-path effect [26,27], 3. non line-of-sight propagation [28,29], and 4. mobile synchronization recovery [30,31].
Therefore, an indoor positioning algorithm which can be easily implemented in a mobile phone using just one microphone is proposed. Unlike in the above design, in which mobile phones were set up as emitters, in this study, the mobile phone’s microphone is arranged as the receiver of tagged sound sources using distinct frequencies. Four tagged sound sources with distinct frequencies are broadcasted from speakers placed in the corners of an indoor space, and for reducing the total cost of the hardware, only one mobile phone microphone is used to collect messages of different tagged sound sources. Regarding the development of the methodology, one novel indoor positioning method combining the received signal strength (RSS) method, the fast Fourier transform (FFT) method, H2 estimation design, and the intersection of circles method is proposed. In this proposed method, RSS is used to calculate the strength of each tagged sound source, FFT is adopted to analyze the spectrum of the collected data of tagged sound sources, the H2 estimation design is used to purify the noisy tagged sound source, and output denoised sound pressure level (SPL), and an accurately estimated position of the user in an indoor environment can then be solved by the intersection of circles method.

2. The Proposed Indoor Positioning Algorithm

The overall schematic of the proposed indoor positioning design is shown as Figure 1. In this investigation, this procedure is separated into three stages. A RSS analysis of the first stage is used to verify the pressure of the received signals measured by the single microphone of the mobile phones; additionally, in this stage, one system identification method—recursive least square (RLS)—and the famous power spectrum analysis tool—FFT—are utilized to simulate the walk behavior of users and make SPL selections regarding which two tagged sound sources will be adopted. In the second stage, one novel estimation design possesses an effective reduction ability regarding environmental corruptions. A set of purified SPLs can be obtained by using this proposed estimation design, and the first two of these four SPLs, which exhibit the strongest intensities, will be used as inputs of geometry equations in the third stages. Based on the purified SPLs, four geometry equations, which are functions of the users’ indoor position ( x r ( k ) , y r ( k ) ) , can be easily determined from the relationship between tagged sound sources and the microphone of the user’s mobile phone. The users’ indoor position ( x r ( k ) , y r ( k ) ) can be further solved by the intersection of circles method in the third stage.

2.1. Data Acquisition and RSS Analysis

In the following, a detailed mathematical expression for these three stages will be derived. In the first stage, the received raw data of the four tagged sound sources are as follows: S n , for n = 1 to 4, with distinct tagged frequencies, should be identified. As to the choice of these two tagged sound sources, there are two steps arranged previously: calculations of magnitudes and frequencies. For identifying the magnitude and frequency, FFT and RSS analysis are used. The tagged frequency of each sound source cannot be higher than 20 kHz due to the physical limitation of the standard microphone used, and the sampling frequency of the mobile phone is 8 kHz; hence, tag frequencies of the sound sources are selected as 15.5 kHz, 16.5 kHz, 17.5 kHz, and 18.5 kHz, respectively. For separating the four tagged sound sources, a specific frequency tag f S n is assigned for each of the four tagged sound sources, as follows:
O n : ( 14000 + 1000 n ) H z < f S n < ( 15000 + 1000 n ) H z ,   for   n = 1   to   4
Denote the received signals of the tagged sound sources S n as z S n , for n = 1 to 4, which contain all surrounding sounds. In this stage, two tagged sound candidates will be selected based on the messages of magnitudes and frequencies of these four tagged sound sources. For the purpose of calculating magnitudes (dB) and frequencies (Hz) of the tagged sound sources, the famous power spectrum method FFT is used. According to the definition of FFT, the magnitude for these four distinct frequencies can be expressed as:
F S n ( k ) = FFT { z S n ( k ) } ,   for   n = 1   to   4
A mean value magnitude adopted for calculation of the intensities of the tagged sound sources is defined as the following:
C S n = 1 B κ = 0 B F S n ( k ) ,   for   n = 1   to   4
where B is the number of bins. Equations (2) and (3) will be utilized to identify the received tagged sound sources.

2.2. Indoor Positioning Algorithm

The geometry relationship between the carried receiver and distinct tagged sound sources is illustrated in Figure 2. In Figure 2, the carried receiver is set up as the mobile phone and is initialized as u r ( 0 ) = { x r ( 0 ) , y r ( 0 ) } . Suppose at least four tagged sound sources, with distinct frequencies, are placed in four corners of an indoor space, and their coordinates are represented as ( x 1 , y 1 ) , ( x 2 , y 2 ) , ( x 3 , y 3 ) , and ( x 4 , y 4 ) .
Figure 2 is an irregular-shaped indoor plane, and there are four tagged sound sources placed in four corners. Tagged sound sources are denoted as S n , for n = 1 to 4, and tag frequencies of these sound sources are set up as f S n . As to the initial locations of the tagged sound sources, they are placed in ( x n , y n ) of the test indoor plane. The moving data of the carried receiver in every time instant is defined as:
u r = { ( x r ( k ) , y r ( k ) ) R 2 } ,   for   1 k N
where N is number of indoor positioning iterations.
In the following, the transformation of SPLs and the relative distances will be detailed.
Data collection of SPLs: The carried receiver collects four different SPL values from the tagged sound sources at each sampling time, which can be expressed as a set LSn:
L S n = { l S n ( k ) R 1 } ,   for   n = 1   to   4 ,   and   1 k N
where l S n is the measured SPL of S n with the tag frequency f S n .
Relative distances between each tagged sound source and the receiver can be expressed as the following four sets.
R S n = { r S n ( k ) R 1 } ,   for   n = 1   to   4 ,   and   1 k N
where r S n ( k ) is the relative distance from the positions of tagged sound sources S n to the receiver. Based on Equations (5) and (6), the instantaneous relative distance between each tagged sound source and the receiver can be calculated from the SPL difference l S n ( k ) and l S n ( k 1 ) as
r S n ( k ) = r S n ( k 1 ) / 10 l S n ( k ) l S n ( k 1 ) 20 + ( R C + 16 π r 2 2 R C + 16 π r 1 2 ) 1 2 ,   for   n = 1   to   4 ,   and   1 k N
where R c is the uncertainty room constant. Due to the fact that R c cannot be measured in prior, it is regarded as a modeling uncertainty. Based on this, the corrupted relative distance r S n ( k ) can be further expressed as below:
r S n ( k ) = r S n ( k 1 ) / 10 l S n ( k ) l S n ( k 1 ) 20 + w ( k ) ,   for   n = 1   to   4 ,   and   1 k N
where initial values of r S n ( k ) and l S n ( k ) are r S n ( 0 ) and l S n ( 0 ) , respectively, and w ( k ) = ( R C + 16 π r 2 2 R C + 16 π r 1 2 ) 1 2 .
The next section will introduce the method of extracting the noiseless SPLs with distinct tagged frequencies and derive the geometry mathematical formulation for finding the instantons positions of the mobile phone users.

2.2.1. Geometric Mathematical Formulation

In Figure 3, there are four right angle triangles inside a quadrilateral, and each of them describes the relationship between the receiver and each tagged sound sources.
Based on Figure 3, four equations are derived, and these equations express the geometric relationships between the receiver and the tagged sound sources in an indoor plane.
( x 1 x r ( k ) ) 2 + ( y 1 y r ( k ) ) 2 = r S 1 ( k ) 2
( x 2 x r ( k ) ) 2 + ( y 2 y r ( k ) ) 2 = r S 2 ( k ) 2
( x 3 x r ( k ) ) 2 + ( y 3 y r ( k ) ) 2 = r S 3 ( k ) 2
( x 4 x r ( k ) ) 2 + ( y 4 y r ( k ) ) 2 = r S 4 ( k ) 2
where ( x 1 , y 1 ) , ( x 2 , y 2 ) , ( x 3 , y 3 ) , and ( x 4 , y 4 ) are the coordinates of the tagged sound sources, and these positions are fixed and known. ( x r ( k ) , y r ( k ) ) is the coordinate of the receiver at time instant k . r S n ( k ) , n = 1   to   4 are the relative distances of the receiver to the tagged sound sources S n at time instant k and can be calculated from Equation (8).
From the above mathematical formulation for the geometric relationship between the tagged sound sources and the receiver, four equations are obtained. However, only two of these four equations will be utilized to solve the two unknowns: ( x r ( k ) , y r ( k ) ) , and two of these four equations—Equations (9)–(12)—will be selected according to the measured SPLs by using Equations (2) and (3). The criterion for selecting the two equations is l M a x = M a x   L S n and l Sec = M a x { L S n M a x   L S n } , e.g., if the intensity sequence of the measured SPLs is in sequence: l S 2 > l S 3 > l S 1 > l S 4 based on Equations (3), (10) and (11), the following will then be selected and combined as a set of binary quadratic equations as:
( x 2 x r ( k ) ) 2 + ( y 2 y r ( k ) ) 2 = r S 2 ( k ) 2
( x 3 x r ( k ) ) 2 + ( y 3 y r ( k ) ) 2 = r S 3 ( k ) 2

2.2.2. Solution by Intersection of Two Circles Method

For solving the solution ( x r ( k ) , y r ( k ) ) of the binary quadratic equations, the intersection of two circles method is applied. The reason for using this method is that it offers calculation convenience and a low computational burden for calculators of mobile phones.
In Figure 4, S 1 and S 2 are tagged sound sources which are selected from Figure 3. The coordinates are ( x 1 , y 1 ) and ( x 2 , y 2 ) , respectively. Two overlapped circles can be illustrated when a radius r S 1 ( k ) and r S 2 ( k ) , which are the relative distance between the receiver and the tagged sound source, are assigned for tagged sound sources S 1 and S 2 , respectively. Suppose r S 1 ( k ) and r S 2 ( k ) are large enough to intersect with each other. Two intersect point, u r 1 ( k ) and u r 2 ( k ) , can be obtained, as shown in Figure 4. Theoretically, one of these two points, ( u r 1 ( k ) and u r 2 ( k ) ), would be the solution of the corresponding binary quadratic equations.
Denoting the distance S 1 S 2 ¯ as d , it can be expressed as:
d = a + b = ( x 2 x 1 ) 2 + ( y 2 y 1 ) 2
where a is S 1 u 0 ¯ and b is S 2 u 0 ¯ .
There are two equations obtained from triangles Δ S 1 u r 1 u 0 and Δ S 2 u r 1 u 0 as
a 2 + h 2 = r S 1 ( k ) 2 b 2 + h 2 = r S 2 ( k ) 2
Using d = a + b , the variable a can be solved as:
a = ( r S 1 ( k ) 2 r S 2 ( k ) 2 + d 2 ) / ( 2 d )
The solution of the point u 0 ( x u 0 , y u 0 ) can be obtained
( x u 0 , y u 0 ) = ( x 1 , y 1 ) + a d ( x 2 x 1 , y 2 y 1 )
Substituting a into Equation (16) to solve h , the coordinate of the intersection point u r ( k ) = ( x r ( k ) , y r ( k ) ) can be obtained as below:
x r ( k ) = x u 0 ± h d ( y 2 y 1 ) y r ( k ) = y u 0 ± h d ( x 2 x 1 )
x r ( k ) = x 1 + a d ( x 2 x 1 ) ± r S 1 ( k ) 2 a 2 d ( y 2 y 1 ) y r ( k ) = y 1 + a d ( y 2 y 1 ) ± r S 1 ( k ) 2 a 2 d ( x 2 x 1 )
In Equation (20), two position solutions of the receiver are obtained simultaneously, but only one of them is the correct answer. There is a simple method to determine whether it is the correct one: the correct solution u r ( k ) must be inside the indoor plane.
However, there are some special cases in which two solutions of Equation (20) are both inside the quadrilateral S 1 S 2 S 3 S 4 . One judgment is proposed for the reasonable selection of these special cases by considering the moving velocity of the users. Normally, a moving velocity of 10 m/s is the maximum limitation for most of users. To calculate the moving velocity by the difference of current positioning solution u r ( k ) and the previous solution u r ( k 1 ) , one checking condition can found for the selection of the correct solution as:
Δ u r ( k ) = u r ( k ) u r ( k 1 ) t 10 m s
where t is the sampling time of the system.
The solution of Equation (20) based on the accurate measurement of SPLs l S n ( k ) , is n = 1   to   4 . Naturally, environmental corruptions, such as reverberation, interference, etc., are always contained in the measurement process of SPLs, hence the SPLs l S n ( k ) should be denoised before calculating Equations (8) and (20). For treating this noise reduction problem regarding the corrupted SPLs l S n ( k ) , an estimation design for effectively removing the environmental corruptions and purifying the SPLs l S n ( k ) is proposed. In the following, a systematic estimation design, combining RLS system modeling and a steady-state optimal estimator, is developed for denoising the corrupted SPLs.

2.3. System Modeling

Before estimating the correct SPLs via adopting an optimal estimation method, each of the received SPLs should first be mathematically modeled. Suppose a set of the measured SPLs can be described as:
L a l l = [ l S 1 , , l S 4 ]
where L a l l is the set of all measured SPL data, and l S n , for n = 1 to 4 is the set of measured SPL data from the tagged sound sources.
The mathematical models of the received SPLs can be expressed as the following white-noise driven difference equation:
l S n ( k ) = i = 1 m a S n i l S n ( k i ) + w S n ( k ) ,   for   n = 1 , , 4
Z S n ( k ) = l S n ( k ) + v S n ( k )
where Z S n ( k ) is the noisy measurement output of the tagged sound source n . Moreover, w S n ( k ) and v S n ( k ) are Gaussian white noise, with a zero mean and which are uncorrelated with l S n ( k ) . a S n i , for i = 1, …, m are identifiable system parameters of the tagged sound source n , and m is the system order.
A regressive form is used to express the difference equation in Equation (23):
l S n ( k ) = Ψ S n ( k ) T λ ^ S n ( k ) + w S n ( k )
where Ψ S n ( k ) = [ l S n ( k 1 )     l S n ( k m ) ] T is the regression vector which contains the measured SPL data, and λ ^ S n ( k ) = [ a S n 1     a S n m ] is the parameter vector.
Equation (25) can be further formulated as a state space form, as follows:
L S n ( k + 1 ) = Ψ S n L S n ( k ) + Λ w S n ( k )
Z S n ( k ) = Ω L S n ( k ) + Π v S n ( k )
where Ψ S n = [ a S n 1 a S n 2 a S n m 1 a S n m 1 0 0 0 0 0 1 0 ] m × m , Λ = [ 1 0 0 0 ] m × 1 , Ω = [ 1 0 0 0 ] T 1 × m , Π = 1 and L S n ( k ) = [ l S n ( k )   l S n ( k 1 )     l S n ( k m + 2 )   l S n ( k m + 1 ) ] T is the state vector, and w S n ( k ) and v S n ( k ) are white noises.
Remark 1.
The system order m is a user decided parameter. Theoretically, a more accurate model can be found by increasing this system order. Surely, a high order system needs more calculation time and storage memory. There exists a tradeoff when selecting the system parameter m. System parameters a S n i , for i = 1 to m can be optimally calculated by several identification methods, such as the recursive least squares method (RLS), the stochastic subspace identification method (SSI), and the system realization using information matrix (SRIM). In this investigation, the system parameters  a S n i in Equation (23) will be optimally determined by the RLS method.
The RLS algorithm is often used for searching the optimal parameters of systems with a set of unknown parameters by using the input and output measured raw data. The standard identification process of this algorithm can be expressed as follows:
ε S n ( k ) = Z S n ( k ) Ψ S n ( k ) T λ ^ S n ( k 1 )
Q ¯ S n ( k ) = 1 γ [ Q ¯ S n ( k 1 ) Q ¯ S n ( k 1 ) Ψ S n ( k ) λ ^ S n ( k ) T Q ¯ S n ( k 1 ) γ + Ψ S n ( k ) T Q ¯ S n ( k 1 ) λ ^ S n ( k ) ]
λ ^ S n ( k ) = λ ^ S n ( k 1 ) + Q ¯ S n ( k ) Ψ S n ( k ) ε S n ( k )
where Q ¯ S n ( k ) is the estimation of the coefficient covariance at time instant k , λ ^ S n ( k ) is the identified parameter, Ψ S n ( k ) is the input data, ε S n ( k ) is the prediction error, Z S n ( k ) is the measurement output, and γ is the forgetting factor. The range of forgetting factor γ is usually given within 0.95 to 1.

2.4. H2 Estimation Design

Based on the identified system parameters λ ^ S n ( k ) in Equation (30), an H2 estimator is proposed for eliminating the influence of the environmental corruptions.
Equations (26) and (27) represent the state-space system of measured SPLs, and the purified SPLs h S n ( k ) can be reconstructed as:
h S n ( k ) = J L S n ( k )
where J is a constant matrix that is set up to draw out the purified SPL h S n ( k ) from state vector L S n ( k ) . The designed target is to hunt for the estimation L ^ S n ( k ) from the measured SPL Z S n ( k ) , which is corrupted by environmental noises; hence, the state estimator for purifying the corrupted SPL is formulated as the following:
L ^ S n ( k + 1 ) = Ψ S n L ^ S n ( k ) + G S n [ Z S n ( k ) Ω L ^ S n ( k ) ] h ^ S n ( k ) = J L ^ S n ( k )
where G S n m × 1 is the designed estimation gain in a steady state.
Define the estimation error between purified SPL signal and estimation signal as follows:
e ˜ S n ( k ) = h S n ( k ) h ^ S n ( k ) = J L S n ( k ) J L ^ S n ( k ) = J L ˜ S n ( k )
where L ˜ S n ( k ) = L S n ( k ) L ^ S n ( k )
The performance index of the H2 estimation design of the indoor positioning problem can be expressed by using the mean square error of estimation error e ˜ S n ( k ) as:
X S n = E { e ˜ S n ( k + 1 ) e ˜ S n ( k + 1 ) T } = E { J L ˜ S n ( k + 1 ) L ˜ S n ( k + 1 ) T J T }
where e ˜ S n ( k + 1 ) = J L ˜ S n ( k + 1 )
Furthermore, Equation (34) can be presented as:
X S n = E { t r ( J L ˜ S n ( k + 1 ) L ˜ S n ( k + 1 ) T J T ) } = t r ( J E { L ˜ S n ( k + 1 ) L ˜ S n ( k + 1 ) T } J T )
From Equation (26), the estimation error L ˜ S n ( k + 1 ) at a steady state can be described as:
L ˜ S n ( k + 1 ) = L S n ( k + 1 ) L ^ S n ( k + 1 ) = Ψ S n L S n ( k ) + Λ w S n ( k ) { Ψ S n L ^ S n ( k ) + G S n [ Z S n ( k ) Ω L ^ S n ( k ) ] } = Ψ S n ( L S n ( k ) L ^ S n ( k ) ) + Λ w S n ( k ) G S n [ Ω L S n ( k ) + Π v S n ( k ) Ω L ^ S n ( k ) ] = Ψ S n L ˜ S n ( k ) + Λ w S n ( k ) G S n [ Ω L ˜ S n ( k ) + Π v S n ( k ) ] = ( Ψ S n G S n Ω ) L ˜ S n ( k ) + Λ w S n ( k ) G S n Π v S n ( k )
After some mathematical derivations, the H2 estimation design for background noise reduction of the indoor positioning design could be summarized as the following Theorem 1.
Theorem 1.
An H2 steady state estimator for the indoor positioning problem can be constructed if a positive-definite matrix  E S n = E S n T can be found such that the following LMIs hold
[ E S n E S n Λ D S n Π ( E S n Ψ S n D S n Ω ) Λ T E S n I 0 0 Π T D S n T 0 I 0 ( E S n Ψ S n D S n Ω ) T 0 0 E S n ] > 0  
where D S n = E S n G S n , and the steady state covariance of the estimation error is bounded by
E { e ˜ S n ( k + 1 ) e ˜ S n ( k + 1 ) T } < t r ( J E S n 1 J T )
Remark 2.
Proof of Theorem 1 is given inAppendix A.
From Equation (37), D S n and E S n can be computationally calculated. Based on these two parameters, the estimation gain of the H2 estimation design can be obtained by using G S n = E S n 1 D S n . Additionally, by substituting G S n into Equation (32), the H2 estimator can be constructed.
The process of constructing the proposed H2 estimation design is summarized in the following steps:
  • Step 1. Given the R S n , T S n as identity matrices, and J as a constant matrix based on the extraction of desired state variables.
  • Step 2. Solve the LMI form of Equation (37) for obtaining the positive matrix D S n and E 2 S n .
  • Step 3. Calculate the estimation gain G S n based on solutions of E S n and D S n in Step 2.
  • Step 4. Substituting the estimation gain G S n into Equation (32), the H 2 estimation design can be constructed as Equation (39).
L ^ S n ( k + 1 ) = Ψ S n L ^ S n ( k ) + G S n [ Z S n ( k ) Ω L ^ S n ( k ) ] h ^ S n ( k ) = J L ^ S n ( k )

3. Verification of Indoor Positioning Performance

To verify the proposed indoor positioning system, simulation results and practical tests will be discussed and compared. Before discussing the positioning performance of this proposed design, the environmental settings of the hardware and software used will be detailed first. Next, the simulation results, which contain one scenario, will be discussed and analyzed. The practical implementation and testing of this proposed design will be verified after the simulation process. Finally, the comparisons of the simulation results and the practical tests will be discussed.

3.1. Testing Environment Setting

3.1.1. Arrangements of Hardware

Hardware adopted for the indoor positioning verification of this proposed design comprises four speakers: Philips AT10, for broadcasting tagged sound sources, one iPhone6 microphone used as a receiver, and a dual core CPU of an iPhone 6 adopted as the calculator of the proposed indoor positioning algorithm. Specifications of the Philips AT10 and the iPhone 6 are displayed in Figure 5 and are listed in Table 1 and Table 2.

3.1.2. Software Design

After choosing the hardware above, the distinct tag frequencies can be calculated by the FFT method using the measured raw data of the receiver. The analog to digital resolution of iPhone 6 is 16 bits; hence, the measured analog sound pressure of these four tagged sound sources can be presented with 2048 bits.
Distinct tagged frequencies are selected from the following frequency spans:
15   k H z < f S 1 < 16   k H z 16   k H z < f S 2 < 17   k H z 17   k H z < f S 3 < 18   k H z 18   k H z < f S 4 < 19   k H z
The setting parameters for the distinct tagged sound source are shown in Table 3, including tag frequencies, relative distances, and SPL values, respectively. These values are defined as the initial values and are fixed.
For analyzing the effects of the inevitable environmental corruptions which infiltrate to the located tagged sound sources. Two background sounds: 1. one recorded frame of a really noisy pub and 2. a song named: Free loops are used as the environmental corruptions for verifying the robust property of this proposed indoor positioning design.

3.2. Test Results

3.2.1. The Configuration for Practical Results

For simulations, a flow chart with three stages each playing specific roles in the indoor positioning process, is illustrated in Figure 6. The corruptions of environments, such as reverberant effects and environmental uncertainties, are set up as random noises, with the magnitude 15 dB for mimicking the practical situations.
The parameters for the proposed H2 estimator in Equation (39) are given in Table 4, and for conserving computational power, the system order m for the SPL model in Equation (23) is selected as 15. The steady state values λ ^ S n ( k ) , G S n , and E S n , for n=1 to 4, of SPL models with respect to tagged sound sources S1 to S4 are listed in Table 4. For saving space, only λ ^ S 1 ( k ) , G S 1 , and E S 1 of the tagged sound source S1 are listed in Table 4.

3.2.2. Testing Scenario

As shown in Figure 7, a prearranged circle path (simulation data) within an indoor area (length: 45 m × width: 40 m), which has 4 tagged sound sources (+) allocated in the four corners of this area, is used. The radius of this prearranged circle path is 10 m, and the initial point of the receiver is at the point (xr = 25 me, yr = 20 m). In Figure 7, the mobile phone user (red dot) walks counterclockwise along the circle path (arrow).
Figure 8 shows the indoor positioning result, which only uses measured SPLs and Equations (6) and (20) to calculate the user’s position (xr, yr), without the help of an estimation design. The root mean square (RMS) error of this simulation is 1.59 m. As for the positioning result of utilizing the H2 estimation design, it is plotted in Figure 9, and the RMS error is 0.77 m.
From Figure 9 and Figure 10, it is obvious that the indoor positioning accuracy can be effectively improved by the proposed estimation design, and the indoor positioning performance of the proposed method is superior to that without any estimation design.

3.3. Practical Experiments

Figure 11 shows the overall testing process of our proposed indoor positioning design. Similarly, this process for practical tests is also divided into three stages. In the first stage, the real time measured SPL data is used to build the mathematical walking model of the mobile phone user, and an RSS analysis using the FFT is adopted for identifying the first two tagged sound sources of the allocated four tagged sound sources. In stage 2, purified SPLs can be delivered, and in the third stage, real time positions of the mobile phone user can be solved by using the intersection of circles method.
The practical experiment of this proposed design is conducted in the building of the Department of Systems and Naval Mechatronic Engineering, National Cheng Kung University, No.1, University Road, Tainan City, Taiwan. The tested building and layouts of the floors are shown as Figure 12. The iPhone6 is chosen as a receiver, and four tagged sound sources are allocated in appropriated positions, according to the requirements of the experiments. The sampling frequency of the received data is 44,100 Hz. Four distinct frequencies: 15.5 kHz, 16.5 kHz, 17.5 kHz, and 18.5 kHz, are selected for the purpose of identifying the detected sound sources. After the above arrangements in the test environment and apparatus, the same test pattern—a circle path—is assigned beforehand to verify the indoor positioning abilities of this proposed estimation design and the scheme without any estimation design.
This real test is conducted in an area with dimensions: (40 m × 45 m) and covered with four tagged sound sources. The relative distance between the tagged sound sources (S1 to S4) is 18 m, as shown in Figure 13, and the cross stars represent the locations of the tagged sound sources. In this scenario, the test path has a circle trajectory, with a radius of 5 m. A total of 7998 SPL data sampled within 40 s are acquired in this test, and each real-time indoor positioning point is drawn per 0.01 Hz. The mobile phone user will walk counterclockwise along this circle path, as illustrated in Figure 14.
Figure 15 shows the total indoor positioning result, without using any estimation design, and the RMS error with respect to the circle reference path is 2.49 m. Figure 16 shows the positioning result of the proposed H2 estimation design, and the RMS error, which corresponds to the circle reference path, is 0.74 m. Figure 17 shows the indoor positioning error histories of the proposed estimation design with respect to the indoor positioning scheme without any estimation design.
From Figure 16 and Figure 17, it can be seen that the proposed indoor positioning design highly improves the positioning accuracy, even under the effect of the worst background noises, and it outperforms the design that uses only measured SPLs, without any estimation design.

3.4. Comparisons of Simulation Results and Practical Results

As shown in Table 5, the RMS error of the indoor positioning result, without the help of estimation designs, is 1.59 m within an area with a length of 45 m and a width of 40 m.
As shown in Table 6, the RMS error of the indoor positioning performance of the H2 estimation design is 0.77 m, based on the same test environment and condition.
In real experiments, the covered areas are similar to those of simulations. Table 7 shows the RMS error of the indoor positioning, without using the estimation design, within an area with a length of 45 m and a width of 40 m.
As for the H2 estimation design, the positioning performance is listed in Table 8. Under the same test condition, the proposed H2 estimation design delivers an indoor positioning performance, with an RMS error of 0.77 m.
From comparisons of Table 5, Table 6, Table 7 and Table 8, it can be seen that the indoor positioning performances of this proposed method in simulation and practical experiments are similar because the noise levels of the environmental corruptions simulated by Equation (8) have been tuned beforehand, according to off-line measured environmental corruptions for the purpose of approximating the noise level of the true corruption in the real environment. Similar results can be obtained for the indoor positioning scheme without any estimation design. As a whole, the proposed indoor positioning design yields better positioning performance than that without any estimation design, and it is robust in regards to the environmental background noises.

4. Conclusions

An effective and accurate indoor positioning function is necessary for users who carry mobile phones because they spend a lot of time in shielding spaces. For achieving this design target, an indoor positioning scheme with the H2 estimation design, which can purify the corrupted SPLs of measured sound sources, is successfully developed for mobile phone users in this investigation. Although the test conditions vary randomly, and indoor environments are challenging due to the unknown background noises, simulation results and practical tests obviously show the promising indoor positioning performance of this proposed method: an average positioning RMS error of 0.75 m can be obtained. Two main contributions can be summarized for this investigation: 1. A compact indoor positioning system, with a high indoor positioning accuracy and capable of execution on mobile phones, is developed because the proposed H2 estimation design possesses the following natural potentials: a low power consumption for computation and an easy-to-implement filter structure. 2. This investigation provides a positioning possibility other than GPS for use in commonly unshielded spaces, as well as indoor environments.

Author Contributions

Conceptualization, Y.-H.C., Y.-Y.C. and P.-Y.C.; methodology, Y.-H.C. and Y.-Y.C.; software, Y.-H.C. and Y.-Y.C.; validation, Y.-H.C., Y.-Y.C. and P.-Y.C.; formal analysis, Y.-H.C., Y.-Y.C. and P.-Y.C.; investigation, Y.-H.C., Y.-Y.C. and P.-Y.C.; resources, Y.-H.C. and Y.-Y.C.; data curation, Y.-H.C., Y.-Y.C. and P.-Y.C.; writing—original draft preparation, Y.-H.C., Y.-Y.C. and P.-Y.C.; writing—review and editing, Y.-H.C. and Y.-Y.C.; visualization, Y.-Y.C. and P.-Y.C.; supervision, Y.-H.C. and Y.-Y.C.; funding acquisition, Y.-H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of Taiwan, grant number MOST 111-2221-E-020-023.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

From Equations (35) and (36), the estimation covariance can be further derived as
E { L ˜ S n ( k + 1 ) L ˜ S n ( k + 1 ) T } = E { [ ( Ψ S n G S n Ω ) L S n ( k ) + Λ w S n ( k ) G S n Π v S n ( k ) ] [ L ˜ S n ( k ) T ( Ψ S n G S n Ω ) T + w S n ( k ) T Λ T v S n ( k ) T Π T G S n T ] } = E { ( Ψ S n G S n Ω ) L ˜ S n ( k ) L ˜ S n ( k ) T ( Ψ S n G S n Ω ) T + Λ w S n ( k ) w S n ( k ) T Λ T + G S n Π v S n ( k ) v S n ( k ) T Π T G S n T } = ( Ψ S n G S n Ω ) E { L ˜ S n ( k ) L ˜ S n ( k ) T } ( Ψ S n G S n Ω ) T + Λ Λ T + G S n Π Π T G S n T
where L S n ( k ) , w S n ( k ) , and v S n ( k ) are mutually orthogonal, and the covariance matrix E { w S n ( k ) w S n ( k ) T } = I m × m and E { v S n ( k ) v S n ( k ) T } = I m × m are assumed as the identity matrix. The other covariance matrix E { L ˜ S n ( k ) L ˜ S n ( k ) T } , at a steady state in practical design, is constant (i.e., k ) and is represented as E { L ˜ S n ( k ) L ˜ S n ( k ) T } = Θ S n .
Combining Equation (35) and Equation (A1), the mean-square error can be rewritten as:
X S n = t r ( J E { L ˜ S n ( k + 1 ) L ˜ S n ( k + 1 ) T } J T ) = t r ( J [ ( Ψ S n G S n Ω ) Θ S n ( Ψ S n G S n Ω ) T + Λ Λ T + G S n Π Π T G S n T ] J T ) = t r ( J [ ( Ψ S n G S n Ω ) Θ S n ( Ψ S n G S n Ω ) T Θ S n + Λ Λ T + G S n Π Π T G S n T ] J T ) + t r ( J Θ S n J T )
Equation (A2) shows that the mean-square error X S n must have an upper bound as below:
X S n t r ( J Θ S n J T )
under the following inequality holds:
( Ψ S n G S n Ω ) Θ S n ( Ψ S n G S n Ω ) T Θ S n + Λ Λ T + G S n Π Π T G S n T < 0
Given that E S n = Θ S n 1 , multiplying both sides of Equation (A4) by E S n and selecting D S n = E S n G S n , the Equation (A4) can be rewritten as:
( E S n Ψ S n D S n Ω ) E 2 S n 1 ( E S n Ψ S n D S n Ω ) 1 E S n + E S n Λ Λ T E S n + D S n Π Π T D S n T < 0
For acquiring solution E S n in Equation (A5), the Schur complement is employing to Equation (A5) for finding out a positive solution E S n , yielding the following equivalent LMI form:
[ E S n E S n Λ D S n Π ( E S n Ψ S n D S n Ω ) Λ T E S n I 0 0 Π T D S n T 0 I 0 ( E S n Ψ S n D S n Ω ) T 0 0 E S n ] > 0
This is Equation (37), and the proof is completed.

References

  1. Commission, P.C. FCC Amended Report to Congress on the Deployment of E-911 Phase II Services Tier III Service Providers, 1st ed.; Wireless Telecommunications; FCC: Washington, DC, USA, 2005. [Google Scholar]
  2. Atia, M.M.; Liu, S.; Nematallah, H.; Karamat, T.B.; Noureldin, A. Integrated Indoor Navigation System for Ground Vehicles With Automatic 3-D Alignment and Position Initialization. IEEE Tran. on Veh. Technol. 2015, 64, 1279–1292. [Google Scholar] [CrossRef]
  3. Constandache, I.; Agarwal, S.; Tashev, I.; Choudhury, R.R. Daredevil: Indoor location using sound. SIGMOBILE Mob. Comput. Commun. Rev. 2014, 18, 9–19. [Google Scholar] [CrossRef]
  4. WHERENET Real-Time Locating System. Available online: http://www.wherenet.com/ (accessed on 1 August 2008).
  5. Casas, R.; Cuartielles, D.; Marco, A.; Gracia, H.J.; Falco, J.L. Hidden Issues in Deploying an Indoor Location System. IEEE Perv. Comput. 2007, 6, 62–69. [Google Scholar] [CrossRef]
  6. NDI Measurement You Can Trust. Available online: http://www.ndigital.com/ (accessed on 1 August 2008).
  7. States, R.A.; Pappas, E. Precision and repeatability of the Optotrak 3020 motion measurement system. J. Med. Eng. Technol. 2006, 30, 11–16. [Google Scholar] [CrossRef] [PubMed]
  8. Fernando, X.N.; Krishnan, S.; Sun, H.; Kazemi-Moud, K. Adaptive denoising at infrared wireless receivers. Prceeding SPIE 2003, 5074, 199–207. [Google Scholar]
  9. Harter, A.; Hopper, A.; Steggles, P.; Ward, A.; Webster, P. The anatomy of a context-aware application. In Proceedings of the 5th Annual ACM/IEEE International Conference on Mobile Computing and Networking, Seattle, WA, USA, 15 August 1999. [Google Scholar]
  10. Chunhan, L.; Yushin, C.; Gunhong, P.; Jaeheon, R.; Seung-Gweon, J.; Seokhyun, P. Indoor positioning system based on incident angles of infrared emitters. Industrial Electronics Society, 2004. ECON 2004. In Proceedings of the 30th Annual Conference of IEEE, Busan, Republic of Korea, 2 November 2004. [Google Scholar]
  11. Kaemarungsi, K.; Krishnamurthy, P. Properties of indoor received signal strength for WLAN location fingerprinting. Mobile and Ubiquitous Systems: Networking and Services, 2004. MOBIQUITOUS 2004. In Proceedings of the First Annual International Conference on, Boston, MA, USA, 26 August 2004. [Google Scholar]
  12. Tsung-Nan, L.; Po-Chiang, L. Performance comparison of indoor positioning techniques based on location fingerprinting in wireless networks. In Proceedings of the 2005 International Conference on Wireless Networks, Communications and Mobile Computing, Maui, HI, USA, 13 June 2005. [Google Scholar]
  13. Chon, H.D.; Jun, S.; Jung, H.; An, S.W. Using RFID for Accurate Positioning. J. Glob. Post. Sys. 2004, 3, 32–39. [Google Scholar] [CrossRef] [Green Version]
  14. Wi-Fi Design, WLAN Planning and Site Survey Tool, Wi-Fi Spectrum Analysis. Available online: http://www.ekahau.com/ (accessed on 1 August 2008).
  15. King, T.; Kopf, S.; Haenselmann, T.; Lubberger, C.; Effelsberg, W. COMPASS: A probabilistic indoor positioning system based on 802.11 and digital compasses. In Proceedings of the 1st International Workshop on Wireless Network Testbeds, Experimental Evaluation & Characterization, Los Angeles, CA, USA, 29 September 2006. [Google Scholar]
  16. Bahl, P.; Padmanabhan, V.N. RADAR: An in-building RF-based user location and tracking system. INFOCOM 2000. In Proceedings of the Nineteenth Annual Joint Conference of the IEEE Computer and Communications Societies, Tel Aviv, Israel, 6 August 2000. [Google Scholar]
  17. Wang, Y.; Jia, X.; Lee, H.K. An indoors wireless positioning system based on wireless local area network infrastructure. In Proceedings of the 6th International Symposium on Satellite Navigation Technology Including Mobile Positioning and Location Services, Melbourne, Australia, 22 July 2003. [Google Scholar]
  18. Kim, Y.; Chon, Y.; Cha, H. Smartphone-Based Collaborative and Autonomous Radio Fingerprinting. IEEE Trans. Syst. Man Cybern. Part C 2012, 42, 112–122. [Google Scholar] [CrossRef]
  19. Saeed, A.; Kosba, A.E.; Youssef, M. Ichnaea: A Low-Overhead Robust WLAN Device-Free Passive Localization System. IEEE J. Select. Top. Sig. Proc. 2014, 8, 5–15. [Google Scholar] [CrossRef]
  20. Gu, Y.; Lo, A.; Niemegeers, I. A survey of indoor positioning systems for wireless personal networks. IEEE Comm. Surv. Tutor. 2009, 11, 13–32. [Google Scholar] [CrossRef] [Green Version]
  21. The Bat Ultrasonic Location System. Available online: http://www.cl.cam.ac.uk/research/dtg/attarchive/bat/ (accessed on 1 August 2008).
  22. Priyantha, N.B. The Cricket Indoor Location System; Massachusetts Institute of Technology: Cambridge, MA, USA, 2005. [Google Scholar]
  23. Priyantha, N.B.; Chakraborty, A.; Balakrishnan, H. The Cricket location-support system. In Proceedings of the 6th annual international conference on Mobile computing and networking, Boston, MA, USA, 1 August 2000. [Google Scholar]
  24. Hoflinger, F.; Hoppe, J.; Zhang, R.; Ens, A.; Reindl, L.; Wendeberg, J.; Schindelhauer, C. Acoustic indoor-localization system for smart phones. In Proceedings of the 2014 11th International Multi-Conference on Systems, Signals and Devices (SSD14), Barcelona, Spain, 11 February 2014. [Google Scholar]
  25. Liu, K.; Liu, X.; Li, X. Guoguo: Enabling Fine-Grained Smartphone Localization via Acoustic Anchors. IEEE Trans. Mob. Comput. 2016, 15, 1144–1156. [Google Scholar] [CrossRef]
  26. Polotti, P.; Sampietro, M.; Sarti, A.; Tubaro, S.; Crevoisier, A. Acoustic localization of tactile interactions for the development of novel tangible interfaces. In Proceedings of the 8th Int. Conference on Digital Audio Effects (DAFX-05), Madrid, Spain, 20 September 2005. [Google Scholar]
  27. Crevoisier, A.; Polotti, P. Tangible acoustic interfaces and their applications for the design of new musical instruments. In Proceedings of the 2005 Conference on New Interfaces for Musical Expression, Vancouver, BC, Canada, 26 May 2005. [Google Scholar]
  28. Qi, Y.; Kobayashi, H.; Suda, H. On time-of-arrival positioning in a multipath environment. IEEE Trans. Veh. Technol. 2006, 55, 1516–1526. [Google Scholar] [CrossRef]
  29. Jin, Y.; O’Donoughue, N.; Moura, J.M. Position location by time reversal in communication networks. In Proceedings of the 2008 IEEE International Conference on Acoustics, Speech and Signal Processing, Las Vegas, NV, USA, 31 March 2008. [Google Scholar]
  30. Carotenuto, R.; Merenda, M.; Iero, D.; Corte, F.G.D. Mobile Synchronization Recovery for Ultrasonic Indoor Positioning. Sensors 2020, 20, 702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Gualda, D.; Rubio, M.C.; Urena, J.; Bachiller, S.; Villadangos. J., M.; Hernandez, A.; Garcia, J.J.; Jimenez, A. LOCATE-US:Indoor Positionong for Mobile Device Using Encoded Ultrasonic Signals, Inertial Sensors and Graph-Matching. Sensors 2021, 21, 1950. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flow chart of the proposed indoor positioning algorithm.
Figure 1. Flow chart of the proposed indoor positioning algorithm.
Sensors 23 01508 g001
Figure 2. Geometric relationship between the carried receivers with respect to tagged sound sources.
Figure 2. Geometric relationship between the carried receivers with respect to tagged sound sources.
Sensors 23 01508 g002
Figure 3. The geometric relationship between the tagged sound sources and the receiver in an indoor plane.
Figure 3. The geometric relationship between the tagged sound sources and the receiver in an indoor plane.
Sensors 23 01508 g003
Figure 4. Intersection of two circles.
Figure 4. Intersection of two circles.
Sensors 23 01508 g004
Figure 5. The used devices: (a) Wireless speaker (Philips AT10). (b) iPhone6.
Figure 5. The used devices: (a) Wireless speaker (Philips AT10). (b) iPhone6.
Sensors 23 01508 g005
Figure 6. The overall flow chart of these proposed indoor positioning designs.
Figure 6. The overall flow chart of these proposed indoor positioning designs.
Sensors 23 01508 g006
Figure 7. A circle path for the indoor positioning test.
Figure 7. A circle path for the indoor positioning test.
Sensors 23 01508 g007
Figure 8. The positioning result without using any estimation design.
Figure 8. The positioning result without using any estimation design.
Sensors 23 01508 g008
Figure 9. The indoor positioning result of the proposed H2 estimation design (simulation).
Figure 9. The indoor positioning result of the proposed H2 estimation design (simulation).
Sensors 23 01508 g009
Figure 10. Histories of positioning errors of the proposed estimation design with respect to the positioning scheme without using any estimation design (simulation).
Figure 10. Histories of positioning errors of the proposed estimation design with respect to the positioning scheme without using any estimation design (simulation).
Sensors 23 01508 g010
Figure 11. The flow chart of the proposed indoor positioning algorithm for practical tests.
Figure 11. The flow chart of the proposed indoor positioning algorithm for practical tests.
Sensors 23 01508 g011
Figure 12. The test building and the layouts of floors for this proposed estimation design.
Figure 12. The test building and the layouts of floors for this proposed estimation design.
Sensors 23 01508 g012
Figure 13. The test environment in the first floor of the building.
Figure 13. The test environment in the first floor of the building.
Sensors 23 01508 g013
Figure 14. The circle path for the indoor positioning test in an indoor plane.
Figure 14. The circle path for the indoor positioning test in an indoor plane.
Sensors 23 01508 g014
Figure 15. The indoor positioning result without any estimation design.
Figure 15. The indoor positioning result without any estimation design.
Sensors 23 01508 g015
Figure 16. The indoor positioning result of the proposed H2 estimation design (real test).
Figure 16. The indoor positioning result of the proposed H2 estimation design (real test).
Sensors 23 01508 g016
Figure 17. Histories of positioning errors of the proposed estimation design with respect to the positioning scheme without using any estimation design (real test).
Figure 17. Histories of positioning errors of the proposed estimation design with respect to the positioning scheme without using any estimation design (real test).
Sensors 23 01508 g017
Table 1. Specifications of the Philips AT10.
Table 1. Specifications of the Philips AT10.
Total output power3 W
Frequency response60–20,000 Hz, ± 3 dB
Signal to noise ratio>70 dBA
MP3 bit rate32–320 kbps
Dimensions239 × 104 × 127 mm
Speaker impedance6 ohm
Table 2. Specifications of iPhone 6.
Table 2. Specifications of iPhone 6.
ChipA8 chip with 64 bits
CPUDual-core 1.4 GHz Typhoon
(ARM v8-based)
Number of microphonesTriple microphones, bottom, front, rear
(only one could be arbitrarily used by the developer.)
Dimensions138.1 × 67 × 6.9 mm
Table 3. Setting values of each tagged sound source.
Table 3. Setting values of each tagged sound source.
Number
n
Sound Source
S n
Frequency
f Sn
Relative Distance
r S n ( 0 )
SPL Value
l S n ( 0 )
1S115.5 kHz0.1 m111.5 dB
2S216.5 kHz0.1 m117.5 dB
3S317.5 kHz0.1 m114.9 dB
4S418.5 kHz0.1 m116.1 dB
Table 4. Initial values of the proposed H2 estimation design with respect to the tagged sound source S1.
Table 4. Initial values of the proposed H2 estimation design with respect to the tagged sound source S1.
VariablesDefinitionsValues
L ^ S 1 Initialization of estimation states 0 m × 1
Q S 1 + Error covariance I m × m
T S 1 System disturbance covariance 0.5
R S 1 Measurement noise covariance 0.01
J Constant matrix I m × m
ε S 1 Positive value 20
F v Positive definite weighting matrix 80
F w Positive definite weighting matrix 80
U Positive definite weighting matrix I m × m
The system steady state parameters λ ^ S 1 and the weighting matrices of H2 estimation design under steady-state conditions.
λ ^ S 1 = [ 0.2132 0.1922 0.1712 0.1503 0.1294 0.1084 0.0875 0.0666 0.0457 0.0248 0.0039 0.017 0.0379 0.0587 0.0796 ] , G S 1 = [ 0.2106 0.9824 2.274 e 4 2.026 e 4 1.778 e 4 1.531 e 4 1.283 e 4 1.035 e 4 7.879 e 5 5.405 e 5 2.932 e 5 4.609 e 6 2.01 e 5 4.479 e 5 6.948 e 5 ]
E S 1 = [ 0.01682 0 0 0 0.01682 0.01682 0.01682 0 0 0 0.01682 ] 15 × 15
Table 5. The positioning results without estimation designs in the simulation.
Table 5. The positioning results without estimation designs in the simulation.
Reference PathCovered AreaRMS Error
Circle45 m × 40 m1.59 m
Table 6. The positioning results of the proposed H2 estimation design in the simulation.
Table 6. The positioning results of the proposed H2 estimation design in the simulation.
Reference PathCovered AreaRMS Error
Circle45 m × 40 m0.77 m
Table 7. The positioning results without any estimation design in the practical experiments.
Table 7. The positioning results without any estimation design in the practical experiments.
Reference PathCovered AreaRMS Error
Circle45 m × 40 m2.49 m
Table 8. The indoor positioning of the proposed H2 estimation design in the practical experiments.
Table 8. The indoor positioning of the proposed H2 estimation design in the practical experiments.
Reference PathCovered AreaRMS Error
Circle40 m × 45 m0.74 m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Y.-H.; Chang, P.-Y.; Chen, Y.-Y. Indoor Positioning Design for Mobile Phones via Integrating a Single Microphone Sensor and an H2 Estimator. Sensors 2023, 23, 1508. https://doi.org/10.3390/s23031508

AMA Style

Chen Y-H, Chang P-Y, Chen Y-Y. Indoor Positioning Design for Mobile Phones via Integrating a Single Microphone Sensor and an H2 Estimator. Sensors. 2023; 23(3):1508. https://doi.org/10.3390/s23031508

Chicago/Turabian Style

Chen, Yung-Hsiang, Pei-Yu Chang, and Yung-Yue Chen. 2023. "Indoor Positioning Design for Mobile Phones via Integrating a Single Microphone Sensor and an H2 Estimator" Sensors 23, no. 3: 1508. https://doi.org/10.3390/s23031508

APA Style

Chen, Y. -H., Chang, P. -Y., & Chen, Y. -Y. (2023). Indoor Positioning Design for Mobile Phones via Integrating a Single Microphone Sensor and an H2 Estimator. Sensors, 23(3), 1508. https://doi.org/10.3390/s23031508

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