Next Article in Journal
Dynamic Synchronous Capture Algorithm for an Electromagnetic Flowmeter
Next Article in Special Issue
A Novel Auto-Sorting System for Chinese Cabbage Seeds
Previous Article in Journal
Hazards in Motion: Development of Mobile Geofences for Use in Logging Safety
Previous Article in Special Issue
A Novel Method of Identifying Paddy Seed Varieties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Quadrilateral Geometry Classification Method and Device for Femtocell Positioning Networks

1
Department of Communications Engineering, Yuan-Ze University, Taoyuan 320, Taiwan
2
Communication Research Center, Yuan-Ze University, Taoyuan 320, Taiwan
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(4), 817; https://doi.org/10.3390/s17040817
Submission received: 13 December 2016 / Revised: 28 March 2017 / Accepted: 6 April 2017 / Published: 10 April 2017
(This article belongs to the Special Issue Innovative Sensing Control Scheme for Advanced Materials)

Abstract

:
This article proposes a normalization multi-layer perception (NMLP) geometry classifier to autonomously determine the optimal four femtocell evolved Node Bs (FeNBs), which can use time difference of arrival (TDOA) to measure the location of the macrocell user equipment (MUE) with the lowest GDOP value. The iterative geometry training (IGT) algorithm is designed to obtain the training data for the NMLP geometry classifier. The architecture of the proposed NMLP geometry classifier is realized in the server of the cloud computing platform, to identify the optimal geometry disposition of four FeNBs for positioning the MUE located between two buildings. Six by six neurons are chosen for two hidden layers, in order to shorten the convergent time. The feasibility of the proposed method is demonstrated by means of numerical simulations. In addition, the simulation results also show that the proposed method is particularly suitable for the application of the MUE positioning with a huge number of FeNBs. Finally, three quadrilateral optimum geometry disposition decision criteria are analyzed for the validation of the simulation results.

1. Introduction

The latest Long Term Evolution-Advanced (LTE-A) standard anticipates the increasing use of small cells (known as femtocells) to provide the geo-location information required to meet the emerging communication and networking needs of future smart city applications [1,2]. An overview of the LTE positioning methods is provided in [3], which includes the analysis of previous indoor localization methods. The aim of the paper [3] is focused on the investigation of the floor detection techniques in indoor environments. Based on the massive multiple-input multiple-output (MIMO) antennas and the millimeter wave communication technologies, the small cell concept will appear on the fifth generation (5G) cellular networks [4]. To realize the seamless coverage, a larger number of small cells have to be densely deployed in Heterogeneous network (HetNet) scenarios. This trend results in an attractive technological study of LTE outdoor localization. In providing location-specific services, the location of the moving macrocell user equipment (MUE) must be known with a high degree of accuracy. Then, the FeNBs can steer the digital beam-forming (DBF) to suppress the interfering signal generated from the uplink transmission of the outdoor MUE. Additionally, it could allow the call of the MUE to be handed to the indoor FeNBs, thereby improving the system efficiency.
MUE localization is conventionally performed using a standalone GNSS/GPS system [5]. However, such systems not only incur a high battery consumption, but may also lack the precision required to provide true location-specific services. In some cases (e.g., indoor environments or built-up urban areas), the GPS signal may be intermittent or too weak to perform localization. Notably, the GPS measurements obtained using the MUE cannot be directly used for localization purposes in femtocell networks, since they cannot be transmitted in real time to the femtocell evolved Node Bs (FeNBs) [6]. Various methods have been suggested for improving the MUE localization performance in wireless networks [7,8]. Thus, in a more recent 3GPP specification (Release 11) [9], it is proposed that MUE positioning be performed by the radio access net-work (RAN) itself, rather than the MUE.
Existing MUE positioning methods, including those based on time difference of arrival (TDOA), direction of arrival (DOA), or hybrid TDOA / DOA measurements, generally utilize a Gauss-Newton Interpolation (GNI) algorithm to estimate the location of the MUE in two-dimensional (2D) space [8,10]. Multi-station TDOA positioning systems use optical communications to send the received signal of each FeNB to a reference FeNB, which then uses this information to calculate the TDOA between itself and each transmitting FeNB. Such systems have the advantages of a high accuracy and low system complexity. Furthermore, the common channel error among the different FeNBs is easily eliminated, and thus, the synchronization problem becomes trivial.
TDOA positioning systems require at least four FeNBs to locate the moving MUE in three-dimensional (3D) space. However, the accuracy of the localization results depends on the geometry disposition of the selected FeNBs. The geometric dilution of precision (GDOP) contour map of four stations is simulated in [11]. It shows that the square distribution shape provides the best coverage area, the lozenge shape is the second best, the other irregular quadrilateral (IQ) shapes are the third best, and the straight line has the worst coverage area. Still, to the best of the author’s knowledge, there is no other optimization approach provided to autonomously determine the optimum geometry disposition of four FeNBs. Accordingly, the present study proposes a normalization multi layer perception (NMLP) geometry classifier to autonomously determine the optimal four FeNBs with a lower GDOP value and an iterative geometry training (IGT) algorithm to obtain the training data for the NMLP geometry classifier. The most well-known neural network is the MLP architecture [12], which is widely used for solving problems related to the classification of the different patterns. Normalizing the inputs can make training faster and reduce the chances of getting stuck in local optima with the exact same outputs that the MLP had before [13]. Having obtained the optimal set of FeNBs, the position of the MUE is estimated using the TDOA method and updated using a tracking filter (such as an extended Kalman filter or an adaptive α-β-γ filter [14], which will not be included in this article), as soon as the MUE transmission cannot be continuously received.
In addition, three optimum geometry disposition decision criteria are analyzed for the square, lozenge, and other quadrilateral shape of the four FeNBs for positioning the moving MUE when two or more of the same quadrilateral shapes are determined. The derived analytical expressions are generally applicable to geometries where the MUE is surrounded by the selected FeNBs. A comparison of the analytical results with simulations using the typical geometries of outdoor positioning systems shows good agreement.
The structure of the paper is described as follows. Section 2 presents the femtocell positioning network model. The architecture of the NMLP geometry classifier and the IGT procedure for generating the training data are described in Section 3. In Section 4, three optimum geometry disposition decision criteria are also considered, to determine the optimum shape of four FeNBs in the square, lozenge, and irregular quadrilateral, respectively. In Section 5, the performance of the proposed NMLP geometry classifier is simulated in the outdoor geo-location scenario. It consists of 16 FeNBs located in two adjacent multi-floor buildings. Finally, Section 6 concludes the paper.

2. Femtocell Positioning Network Model

The scenarios under study include the combination of a macro cell with femtocells in outdoor environments, which are operated at the same frequency band. It is assumed that the interference coordination and smart digital beam forming techniques may be needed to solve the problems of co-channel interference and signal attenuations. Figure 1 shows an illustrative femtocell network consisting of multiple FeNBs deployed throughout two adjacent multi-floor buildings. An assumption is made that the FeNBs are connected to a cloud-computing platform through the Internet. Moreover, the coordinates of each FeNB and the training data are pre-stored in the database of the cloud-computing platform. The training data of the NMLP geometry classifier are determined using an IGT algorithm based on the TDOA measurements obtained by each FeNB and the GDOP metric. The architecture of the NMLP geometry classifier is realized in the server of the cloud-computing platform to perform the searching of the optimal set of FeNBs for positioning purposes.
The uplink transmissions from the MUE are conventionally handled using the Single Carrier-frequency division multiple access (SC-FDMA) scheme [15]. Moreover, the signals transmitted from the MUE are detected by digital time delay estimation (TDE) receivers [16], utilizing a digital cross-correlation technique in the frequency domain. In the multi-FeNB TDOA positioning system proposed in this study, the received signal of each FeNB is sent through the Internet to a reference FeNB, which calculates the TDOA between itself and each FeNB such that the position of the MUE can be further derived.
Assume that the moving MUE emits an SC-FDMA signal s[n] in 3D space. The received signal at the reference FeNB (FeNB0) thus has the form:
x 0 [ n ] = s 0 [ n ] + w 0 [ n ] , n [ 0 , M 1 ]
where w 0 [ n ] is the discrete white Gaussian noise at FeNB0; and the variable n is defined as the sampling instant. The received signals at three other arbitrarily-chosen FeNBs are given by:
x i [ n ] = s i [ n D i 0 ] + w i [ n ] n [ 0 , M 1 ] ,   i = 1 , 2 , 3
where w i [ n ] is the discrete white Gaussian noise at FeNBi. The discrete time delay between FeNBi and FeNB0 is estimated as:
l ^ i = D i 0 ,   i = 1 , 2 , 3
Four-FeNB TDOA localization schemes calculate the possible MUE location based on the TDOA measurements of three FeNB receiver pairs lying on separate hyperbola [10]. In other words, the MUE location is obtained by solving three hyperbolic equations, i.e.:
r i , 0 = ( x ^ x 0 ) 2 + ( y ^ y 0 ) 2 + ( z ^ z 0 ) 2 ( x ^ x i ) 2 + ( y ^ y i ) 2 + ( z ^ z i ) 2 ;   i = 1 , 2 , 3  
where ( x i , y i , z i ) and ( x 0 , y 0 , z 0 ) are the coordinates of the i th and reference FeNBs, respectively; and ( x ^ , y ^ , z ^ ) is the estimated MUE location. In addition, r i , 0 is the differential distance between the MUE and the master FeNB, and the MUE and FeNBi, respectively, and is obtained from the measured time delay between FeNBi and FeNB0 as:
r i , 0 = c D i 0 M f s + e r i , 0 = c τ ^ i , 0 + e r i , 0 = c ( τ ^ i τ ^ 0 ) + e r i , 0      i = 1 , 2 , 3
where M is the size of the FFT, τ ^ i , 0 is the analog time delay between FeNBi and FeNB0, and f s is the sampling frequency of the analog-to-digital converter. Finally, e r i , 0 is the measurement error with standard deviation σ e i .
The accuracy of range-based positioning location (PL) systems depends on the geometric relationship between the locations of the FeNBs and the location of the MUE, since errors arising from channel or hardware perturbations in the TDOA estimates are transformed into geographic position errors in the hyperbolic solution [10]. The positioning accuracy of different geometric FeNB configurations can be evaluated using the GDOP parameter, which is defined as the ratio of the root mean square (RMS) position error to the RMS ranging error. The GDOP for an unbiased estimator and a 3D hyperbolic system are given by [10]:
GDOP = 1 N n = 1 N ( x x ^ n ) 2 + ( y y ^ n ) 2 + ( z z ^ n ) 2 σ e i
where N is the number of measurements; σ e i is the standard deviation of the measurement error e r i , 0 ; and ( x , y , z ) and ( x ^ n , y ^ n , z ^ n ) are the true location and the nth estimated location of the MUE, respectively. To calculate the GDOP, it is necessary to know the positions of the FeNBs around the moving MUE. A symbol table is shown in Table 1, which is used to keep track of the symbols and notations used in the paper.

3. Principle of NMLP Geometry Classifier

The architecture of the NMLP geometry classifier is shown in Figure 2. It consists of input data normalization, data scaling, an input layer, two hidden layers, and an output layer. The input vectors include:
s i 0 s i 1 = ( x i 1 x i 0 , y i 1 y i 0 , z i 1 z i 0 ) s i 0 s i 2 = ( x i 2 x i 0 , y i 2 y i 0 , z i 2 z i 0 ) s i 0 s i 3 = ( x i 3 x i 0 , y i 3 y i 0 , z i 3 z i 0 ) s i 1 s i 3 = ( x i 3 x i 1 , y i 3 y i 1 , z i 3 z i 1 )
where the coordinates of P set of the four FeNBs are S i = s i 0 ( x 0 , y 0 , z 0 ) s i 1 ( x 1 , y 1 , z 1 ) s i 2 ( x 2 , y 2 , z 2 ) s i 3 ( x 3 , y 3 , z 3 ) ;   i   =   1 , 2 , ... , P . , where s i j is the jth component of the ith input vector and P is the total number of four FeNBs in the coverage range of the target MUE. Note that each component in a vector of a quadrilateral is constituted by four FeNBs. The data normalization is obtained by:
S i 0 S i 1 = S i 0 S i 1 L i , max = ( ( x i 1 x i 0 ) , ( y i 1 y i 0 ) , ( z i 1 z i 0 ) ) = ( d i 1 , d i 2 , d i 3 ) S i 0 S i 2 = S i 0 S i 2 L i , max = ( ( x i 2 x i 0 ) , ( y i 2 y i 0 ) , ( z i 2 z i 0 ) ) = ( d i 4 , d i 5 , d i 6 ) S i 0 S i 3 = S i 0 S i 3 L i , max = ( ( x i 3 x i 0 ) , ( y i 3 y i 0 ) , ( z i 3 z i 0 ) ) = ( d i 7 , d i 8 , d i 9 ) S i 1 S i 3 = S i 1 S i 3 L i , max = ( ( x i 3 x i 1 ) , ( y i 3 y i 1 ) , ( z i 3 z i 1 ) ) = ( d i 10 , d i 11 , d i 12 )
where L max = maximum   of { s i 0 s i 1 , s i 1 s i 2 , s i 2 s i 3 , s i 3 s i 0 } .
The scaled input and output data are required for the MLP, in order to provide the correct value for an activation function. The linear scaling of data is the transformation from the original data variability within the interval [ d q , min , d q , max ] to the interval [0.1, 0.99], and can be performed using the Formula [10]:
η i q = ( d i q d q , min ) ( b a ) d q , max d q , min ,   q = 1 , ... , Q ; i = 1 , ... , P
where [ a , b ] = [ 0.1 , 0.99 ] . d i q is the original length of the ith input component at the qth neuron of the input layer. The MLP network of the geometry classifier is composed of 12 (Q) inputs, which represent 12 features of the NMLP geometry classifier.
The sigmoid activation function [12] used is non-linear, differentiable, and is defined by:
y ( a n ( k ) ) = 1 / ( 1 + e α a n ( k ) )
the neuron outputs of last layer where m l = 5 is the total number of neurons in the last layer.
Given the scaled training data: D = {( S i 1 ( k ) , D i 1 ( k ) ) , ..., ( S i m l ( k ) , D i m l ( k ) ) } where S i q ( k ) , D i q ( k ) are the given input and the desired output at the kth iteration.
The weight matrices are updated by:
W n m ( s ) ( q + 1 ) = W n m ( s ) ( q ) + μ ( s ) ( q ) δ n ( s ) ( q ) x o u t , n ( s 1 ) ( q )
where the learning rate is given by:
μ ( v ) = μ 0 / ( 1 + ( k / t ) )   ( 100 < t < 500 )
For the hidden layer:
δ j ( s ) ( q ) = ( Σ h = 1 n s + 1 δ h ( s + 1 ) ( q ) w h j ( s + 1 ) ( q ) ) g ( v j ( s ) ( q ) )
For the output layer:
δ j ( s ) ( q ) = ( d o u t p u t ( q ) x o u t , j ( s ) ( q ) ) g ( v j ( s ) ( q ) )
During the training, the network weights are adjusted in order to reduce the mean square error (MSE) obtained by [12]
MSE = Σ q = 1 m l ( s i q d i q ) 2 m l
The performance of the NMLP geometry classifier will be evaluated in terms of the percentages for correct classification. It is defined on the difference between the desired output and the estimated output of the NMLP geometry classifier. The classification accuracy is defined as [12]:
classification   accuracy = N c P × 100 %
where N c is the total number of correct classifications and P is the total number of four FeNBs within the sphere.
The flowchart of the IGT algorithm is shown in Figure 3, which is used to generate the training data for the NMLP geometry classifier. As shown in Figure 3, the IGT algorithm comprises five steps. In the first step, four FeNBs are chosen such that the distance between each pair of FeNBs is less than the pre-specified diameter of a sphere (expressed in meters). In other words, four FeNBs are selected such that the following equations are satisfied:
( x i x 0 ) 2 + ( y i y 0 ) 2 + ( z i z 0 ) 2 2 R   ,   i = 1 , 2 , 3 ( x i x i + 1 ) 2 + ( y i y i + 1 ) 2 + ( z i z i + 1 ) 2 2 R   ,   i = 1 , 2 ( x i x i + 2 ) 2 + ( y i y i + 2 ) 2 + ( z i z i + 2 ) 2 2 R   ,   i = 1
where ( x i , y i , z i ) T i = 0 , 1 , 2 , 3 are the coordinates of the four FeNBs; R is the sphere radius; and T denotes the transpose of the matrix. In the second step, the selected set of FeNBs is updated as required, to ensure that all four FeNBs are located on the same plane. Note that the sufficient condition for such a case occurs when the parallelepiped volume constructed by the three vectors s 0 s 1 , s 0 s 2 , s 0 s 3 is equal to zero. In other words:
x 1 x 0 y 1 y 0 z 1 z 0 x 2 x 0 y 2 y 0 z 2 z 0 x 3 x 0 y 3 y 0 z 3 z 0 = 0
where the FeNB located at coordinates ( x 0 , y 0 , z 0 ) T is chosen as the reference station. In the third step, the geometry disposition of the four FeNBs is evaluated to ensure that the FeNBs form either a square or a lozenge shape. In this case, the sufficient condition is taken as two of the vectors having a perpendicular orientation to one another. In other words, the inner product of the two vectors is equal to zero, i.e.:
s 0 s 2 s 1 s 3 = 0
In the fourth step, a check is made to confirm that the four FeNBs are not located on a straight-line.
If four FeNBs are located on a straight-line, this increases the radius, and if not, then the shape is an IQ. In general, two vectors are parallel if their outer product is equal to zero. Thus, the four FeNBs are not located on a straight-line if A B C is valid, where A, B, and C are defined, respectively, as follows:
A :   s 0 s 1 × s 0 s 2 0 B :   s 0 s 1 × s 0 s 3 0 C :   s 0 s 2 × s 0 s 3 0
If the conditions associated with Steps 1, 2, and 3 are satisfied, the optimal four-FeNB disposition has a low GDOP. Otherwise, if the conditions associated with Steps 1, 2, and 4 are satisfied, the acceptable four-FeNB configuration is determined with a higher GDOP. A lower GDOP represents a better localization result. As soon as either case is satisfied, the coordinates of the moving MUE are estimated using the TDOA method. Finally, in the fifth step, a check is made to confirm that the four FeNBs and MUE are located within the designated sphere, i.e.:
( x i x ^ ) 2 + ( y i y ^ ) 2 + ( z i z ^ ) 2 [ 200 + ( r 1 ) × 20 ] 2 = R 2     i = 0 , 1 , 2 , 3   ,   r = 1 , 2 , 3 ,
where the coverage radius R of the MUE is initiated with 200 m and increased with the increment of 20 m.

4. Optimum Geometry Disposition Decision Criteria

It was noted that the optimum quadrilateral shapes, consisting of four FeNBs, can be selected to measure the location of the moving MUE. The square is the best, the lozenge is next, and the IQ is the worst. When two or more of the same quadrilateral shapes are selected, the following three optimum geometry disposition decision criteria are considered, to determine the best accuracy of the TDOA positioning measures. The GDOP formulas of three optimum geometry disposition decision criteria are derived. Since the location of four FeNBs is fixed, the quadrilateral shape of the FeNBs does not change due to the movement of the MUE. Therefore, in order to simplify the proof of selecting the optimum shapes of four FeNBs, three different criteria are proved in the two dimensional space, and the results are still valid in the three dimensional space.
Criterion 1: Square optimum geometry disposition decision criterion depends on the distance between the MUE and the square center. When the distance is shorter, the positioning accuracy is better.
Let the distance from the MUE to the ith FeNB be given by [11]:
r i , T D O A = ( x x i ) 2 + ( y y i ) 2 ,   i = 0 , 1 , 2 , 3
The differential distance is obtained from subtracting two TDOA measurements between FeNBi and FeNB0.
Δ r i , T D O A = r i , T D O A r 0 , T D O A ,   i = 1 , 2 , 3
The TDOA measurements for three pairs of FeNBs are given as:
t i t 0 = Δ r i , T D O A c + ( ε i ε 0 ) c i = 1 , 2 , 3
where c is the light velocity and the measurement error is expressed as:
ε = [ ε 0 ε 1 ε 2 ε 3 ]
It is assumed that the errors in all actual measurements are random, independent, and have an identical root mean square (RMS) value σ r 2 .
E [ ε i ] = ε ¯ E [ ε i ε j ] = 0 E [ ε i 2 ] = σ r 2 + ε ¯ 2 , i = 0 , 1 , 2 , 3 σ r 2 = E [ ( ε i E [ ε i ] ) 2 ] = E [ ε i 2 ] ε ¯ 2
The TDOA measurement errors for three pairs of FeNBs can be expressed as the matrix form.
[ ( ε 1 ε 0 ) ( ε 2 ε 0 ) ( ε 3 ε 0 ) ] = [ 1 1 0 0 1 0 1 0 1 0 0 1 ] [ ε 0 ε 1 ε 2 ε 3 ] = A ε
The error covariance matrix is:
Q = E [ ( A ε ) ( A ε ) T ] = A E [ ε ε T ] A T = ( σ r 2 + ε ¯ 2 ) A A T = ( σ r 2 + ε ¯ 2 ) [ 2 1 1 1 2 1 1 1 2 ]
The partial derivatives of the noise free measurement Equations (23) with respect to the unknown MUE coordinates (x, y) are given as:
d Δ r i = ( x x i r i , T D O A x x 0 r 0 , T D O A ) d x + ( y y i r i , T D O A y y 0 r 0 , T D O A ) d y ;   i = 1 , 2 , 3
The matrix is expressed as:
[ d Δ r 1 , T D O A d Δ r 2 , T D O A d Δ r 3 , T D O A ] = [ x x 1 r 1 , T D O A x x 0 r o , T D O A + y y 1 r 1 , T D O A y y 0 r o , T D O A x x 2 r 2 , T D O A x x 0 r o , T D O A + y y 2 r 2 , T D O A y y 0 r o , T D O A x x 3 r 3 , T D O A x x 0 r o , T D O A + y y 3 r 3 , T D O A y y 0 r o , T D O A ] [ d x d y ]
where:
H = [ x x 1 r 1 , T D O A x x 0 r o , T D O A x x 2 r 2 , T D O A x x 0 r o , T D O A x x 3 r 3 , T D O A x x 0 r o , T D O A y y 1 r 1 , T D O A y y 0 r o , T D O A y y 2 r 2 , T D O A y y 0 r o , T D O A y y 3 r 3 , T D O A y y 0 r o , T D O A ]
When the error covariance matrix is not diagonal, as in Equation (28), the GDOP will be given by:
GDOP = t r a c e ( G )
where the matrix G is the variance of the linear optimum unbiased estimator divided by the variance of measurement noise [17].
G = 1 σ r 2 ( H T Q 1 H ) 1
For the simplification of proving the square optimum geometry disposition decision criterion, the coordinates of FeNB0 (1, 0), FeNB1 (0, 1), FeNB2 (−1, 0), FeNB3 (0, −1), and MUE (x, y) are substituted into Equation (29), to yield:
H = [ x x 2 + ( y 1 ) 2 x 1 ( x 1 ) 2 + y 2 x + 1 ( x + 1 ) 2 + y 2 x 1 ( x 1 ) 2 + y 2 x x 2 + ( y + 1 ) 2 x 1 ( x 1 ) 2 + y 2 y 1 x 2 + ( y 1 ) 2 y ( x 1 ) 2 + y 2 y ( x + 1 ) 2 + y 2 y ( x 1 ) 2 + y 2 y + 1 x 2 + ( y + 1 ) 2 y ( x 1 ) 2 + y 2 ]
When the MUE moves from (1, 0) to (0.1, 0), the GDOPs for different average TDOA position measurement errors (0, 1, 2, 3, 4 m) are computed and these are shown in Figure 4. It shows that the GDOP decreases with the distance between the MUE and the square center. The GDOP also decreases with the average TDOA position measurement error.
Criterion 2: Lozenge optimum geometry disposition decision criterion depends on the complementary angle difference (CAD), which is the difference between two complementary angles of a lozenge shape. When the CAD of the lozenge optimum geometry disposition is smaller, the positioning accuracy is better.
For simplification, the GDOP formula is derived by substituting the coordinates of FeNB0 (1, 0), FeNB2 (−1, 0), FeNB1 (0, N), FeNB3(0, −N), and MUE (0.1, 0) into Equation (31), to yield:
H = [ 1 10 N 2 + 0.01 + 1 2 1 10 N 2 + 0.01 + 1 N N 2 + 0.01 0 N N 2 + 0.01 ]
Substituting Equations (28) and (35) into Equation (33) yields:
G = [ ( 5 + ε ¯ 2 ) ( N 2 + 0.01 ) 200 N 2 + 3 0 0 ( 5 + ε ¯ 2 ) ( N 2 + 0.01 ) 10 N 2 ]
Then the GDOP of lozenge shape is expressed as:
G D O P = G 1 , 1 + G 2 , 2 = ( ε ¯ 2 + 5 ) ( N 2 + 0.01 ) 200 N 2 + 3 + ( ε ¯ 2 + 5 ) ( N 2 + 0.01 ) 10 N 2
The CAD formula is derived as:
CAD = 2 ( arccos ( 1 1 + N 2 ) arccos ( N 1 + N 2 ) )
Figure 5 is generated from Equations (37) and (38), where N varies from 0 to 1 and σ r 2 = 5. It shows that the GDOP decreases with the CAD. The GDOP also decreases with the average TDOA position measurement error. When N = 1, CAD = 0°. This is the square shape, which has a smaller GDOP value than the lozenge shape. When N = 0, CAD = 180°. This is the straight line, which has the largest GDOP value.
Criterion 3: Irregular-quadrilateral optimum geometry disposition decision criterion depends on the distance between the diagonal intersection and the center of the IQ. When the distance is shorter, the positioning accuracy of the IQ optimum geometry disposition is better.
For simplification, the GDOP formula is derived by substituting the coordinates of FeNB0 (0.5, 0), FeNB1 (0, y 1 ), FeNB2 (−0.5, 0), FeNB3 (0, −1), and MUE (0.1, 0) into Equation (31), to yield:
H = [ 1 10 y 1 2 + 0.01 + 9 4 4.083 2.3495 y 1 y 1 2 + 0.01 0 0.995 ] ,
Figure 6 is generated from Equations (28), (32), (33), and (39), where y 1 varies from 0 to 1. It shows that the GDOP decreases with the distance between the diagonal intersection and the center of the quadrilateral and the average TDOA position measurement error.

5. Simulations

In performing the simulations, it was assumed that the time delay between FeNBi and FeNB0 was successfully measured with a detection probability of Pd = 0.9, while the false alarm probability was equal to PFA = 10−4. Moreover, the signal-to-noise ratio (SNR) of the TDE receivers with a linear optimum filter was assumed to be −13 dB over the Rayleigh channels, with Pd = 0.9 and PFA = 10−4. The simulations considered the outdoor geo-location scenario shown in Figure 1, consisting of 16 FeNBs located in adjacent multi-floor buildings. The total number of four FeNBs in the coverage range of the target MUE is P = 1820. Moreover, the moving MUE was assumed to be located between two buildings and positioned at the boundary of a macrocell network. The coordinates of the MUE trajectory are (60, 0, 0), (60, 10, 0), (60, 20, 0), (60, 30, 0), (60, 40, 0), (60, 50, 0), (60, 60, 0), (60, 70, 0), (60, 80, 0), and (60, 90, 0). In simulations, the SC-FDMA signal was considered as it is used by LTE MUE as an uplink access scheme [7,15]. The SC-FDMA waveforms are specified with bandwidth 10 MHz, subcarrier spacing 15 KHz, FFT/IFFT size 1024, and preamble length 320. The sixteen FeNBs in accordance with their coordinate information are shown in Table 2. (Note that the coordinates are expressed in units of meters.) Any four FeNBs constituted a quadrilateral group. Moreover, the MUE is assumed to be moving between two adjacent buildings. The simulation tool is Matlab. After the simulations, the performance of the proposed NMLP geometry classifier using the multi-station TDOA MUE localization scheme was numerically evaluated.
The simulations are performed as follows. First, the IGT algorithm is used to calculate the training data of all quadrilateral FeNBs. Then, the NMLP geometry classifier is trained to determine the optimal shape by searching all of the quadrilateral sets of FeNBs within the detection range of the MUE signal. The weight matrices of MLP networks for the first, second, and third layers are W n m ( 2 ) , W n m ( 3 ) , and W n m ( 4 ) , respectively. The matrix dimensions are 12 × 6, 6 × 6 and 6 × 5, respectively. The desired outputs obtained from the training results are classified as five different shapes, which include [0.99, 0.01, 0.01, 0.01, 0.01], [0.01, 0.99, 0.01, 0.01, 0.01], [0.01, 0.01, 0.99, 0.01, 0.01], [0.01, 0.01, 0.01, 0.99, 0.01], and [0.01, 0.01, 0.01, 0.01, 0. 99] for the square, lozenge, IQ, straight-line, and non-planar shapes, respectively. Finally, the optimum geometry disposition decision criteria are used to determine the best accuracy of the TDOA positioning measures.
Figure 7 shows that the MLP network using six by six neurons for the hidden layers has the shortest convergence time. It achieved a classification accuracy of 98.0659%. Figure 8 shows that the NMLP geometry classifier can improve the convergent time, and the classification accuracy achieved was 98.0330%. Table 3 shows the estimated NMLP network outputs of five neurons for five different shapes. The estimated NMLP geometry classifier outputs are [0.9418, 0.0186, 0.0587, 0.0145, 0.0173], [0.0504, 0.8951, 0.0751, 0.0291, 0.0912], [0.0124, 0.0186, 0.9427, 0.0624, 0.0174], [0.0107, 0.0181, 0.0114, 0.9755, 0.0130], and [0.0178, 0.0136, 0.0584, 0.0175, 0.9581] which are classified as square, lozenge, IQ, straight-line, and non-planar shapes, respectively.
The classification results after training are shown in Table 4, which shows the number of geometry classifications for three cases. For all FeNBs, a unique square geometry is selected by the proposed NMLP geometry classifier to estimate the current MUE position obtained using the TDOA scheme. Here, the TDOA measurement error of each FeNB is assumed to be Gaussian distributed with a mean of zero and a standard deviation of σ e i = 5   m . Moreover, the true coordinates of the MUE are assumed to be (60, 50, 0). The GDOP value is calculated as 1.533 m. For the case not including FeNB4 in the adjacent multi-floor buildings, P reduces to 1465, and two lozenge geometry groups are selected by the proposed NMLP geometry classifier. The CAD of the first lozenge geometry group is 0.8° and the CAD of the second lozenge group is 3.1°. The GDOP value of the former group is 2.012 m and the GDOP value of the latter group is 2.432 m. Then, four FeNBs of a lozenge geometry with a GDOP value of 2.012 m are determined by the proposed NMLP geometry classifier. The decision criterion 2 is verified. The positioning accuracy is better when the CAD of the lozenge optimum geometry disposition is smaller. The TDOA measurement error of the lozenge geometry group is larger than the square geometry group if the FeNB4 is not deployed in the adjacent multi-floor buildings. For the case not including FeNB4 and FeNB8 in the adjacent multi-floor buildings, P reduces to 1001, and all of the square and lozenge geometry groups will be removed. In accordance with the IQ optimum geometry disposition decision criterion, the IQ group with the minimum GDOP value of 4.354 will be selected from 263 IQ groups.
When the MUE moves from the position (60, 0, 0) to (60, 90, 0), the trajectory of the moving path is shown in Figure 1. The proposed NMLP geometry classifier determines the first three best shapes for the three cases. The GDOP values of the first three best quadrilateral shapes for all FeNBs are listed in Table 5, which shows that the GDOP of the square geometry group is smaller when the distance between the MUE and the square center is shorter. The GDOP of the square geometry group is smaller than the lozenge geometry. When two lozenge geometry groups are selected by the proposed NMLP geometry classifier, the GDOP for the CAD of 0.8° is smaller than the CAD of 3.1°. In this case, the former is determined. Both decision criteria 1 and 2 are verified by the simulation results of Table 5. Table 6 shows the GDOP values of the first three best quadrilateral shapes for the case not including FeNB4 in the adjacent multi-floor buildings. It shows that the GDOP of the lozenge geometry group is smaller than the IQ geometry. The GDOP values of the first three best quadrilateral shapes for the case not including FeNB4 and FeNB8 in the adjacent multi-floor buildings are shown in Table 7, which demonstrates that the GDOP values are smaller when the distance between the diagonal intersection and the center of the IQ is shorter. The decision criterion 3 is verified. When four FeNBs with coordinates (0, 0, 50), (0, 180, 50), (180, 0, 50), and (180, 180, 50) are added into the buildings, the total number of four FenBs in the coverage range of the target MUE increases to P = 4845. The simulation results show that the GDOP value is calculated as 1.053 m and the classification accuracy achieved is 97.196%.

6. Conclusions

This paper proposes a novel NMLP geometry classifier and IGT algorithm for autonomously reducing the positioning error of four-FeNB TDOA measurements. Three optimum geometry disposition decision criteria are analyzed for the square, lozenge, and other IQ shape of the four FeNBs, for positioning the moving MUE. The derived analytical expressions are generally applicable to geometries where the MUE is surrounded by the selected FeNBs. The simulation results have confirmed that the proposed NMLP geometry classifier and optimum geometry disposition decision criteria can provide accurate outdoor geo-location information on the MUE for indoor femtocells to support its quality of services in HetNet. A comparison of the analytical results with simulations using the typical geometries of outdoor positioning systems shows good agreement.
The computational burden is mainly caused by the offline IGT algorithm simulations pertaining to the construction of the training data. Once the initial learning of the NMLP is finished using the collected training data to optimize network weights based on the NMLP network outputs of five neurons for five different shapes, the additional computation overhead in NMLP is light. Since the operations of the NMLP geometry classification are executed in the server, it is assumed that the computation ability of the server is powerful enough to determine the optimal geometry disposition for femtocell positioning networks. The information regarding the training data that governs the NMLP geometry classification operation of femtocell positioning networks is stored in a database. As soon as a new operation condition occurs, the new training data is generated from the simulation and included in the database to rerun the IGT training procedure of Figure 3. The simulation results show that the proposed method is particularly suitable for the application of MUE positioning with a huge number of FeNBs.
How to optimally determine four FeNBs to position the UE in indoor scenarios is another hot research topic in current telecommunication industries and academics. The proposed NMLP geometry classifier and IGT algorithm will be applied to the study of indoor femtocell positioning networks.

Acknowledgments

The research work was supported by Ministry of Science and Technology, Taiwan, R. O. C. (MOST 105-2218-E-155-001).

Author Contributions

Jeich Mar conceived, designed and wrote the paper; Tsung Yu Chang and Yu Jie Wang performed the experiments; Jeich Mar, Tsung Yu Chang and Yu Jie Wang analyzed the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cimmino, A.; Pecorella, T.; Fantacci, R.; Granelli, F.; Rahman, T.F.; Sacchi, C.; Carlini, C.; Harsh, P. The Role of Small Cell Technology in Future Smart City Applications. Trans. Emerg. Telecommun. Technol. 2013, 25, 11–20. [Google Scholar] [CrossRef]
  2. Boudreau, G.; Panicker, J.; Guo, N.; Wang, N.; Vrzic, S. Interference Coordination and Cancellation for 4G Networks. IEEE Commun. Mag. 2009, 47, 74–81. [Google Scholar] [CrossRef]
  3. Del Peral-Rosado, J.A.; Bavaro, M.; Lopez-Salcedo, J.A.; Seco-Granados, G.; Chawdhry, P.; Fortuny-Guasch, J.; Crosta, P.; Zanier, F.; Crisci, M. Floor Detection with Indoor Vertical Positioning in LTE Femtocell Networks. In Proceedings of the 2015 IEEE Globlecom Workshops, San Diego, CA, USA, 6–10 December 2015; pp. 1–6. [Google Scholar]
  4. Wang, C.-X.; Haider, F.; Gao, X.; Yang, Y.; Yuan, D.; Aggoune, H.; Haas, H.; Fletcher, S.; Hepsaydir, E. Cellular Architecture and Key Technologies for 5G Wireless Communication Networks. IEEE Commun. Mag. 2014, 52, 122–130. [Google Scholar] [CrossRef]
  5. A 3GPP Description Document-Overview of 3GPP Release 9. 4 September 2014, Volume 3. Available online: www.3gpp.org (accessed on 28 January 2017).
  6. Zhang, J.; de la Roche, G. Femtocells: Technologies and Deployment; John Wiley: New York, NY, USA, 2010. [Google Scholar]
  7. Sesia, S.; Toufik, I.; Baker, M. LTE: The UMTS Long Term Evolution: From Theory to Practice, 2nd ed.; John Wiley and Sons: New Jersey, NJ, USA, 2009; pp. 423–436. [Google Scholar]
  8. 3GPP TS 36.305. Stage 2 Functional Specification of UE Positioning in E-UTRAN, Release 9; 3GPP: Valbonne, France, 2010. [Google Scholar]
  9. Johansson, J.; Hapsari, W.A.; Kelley, S.; Bodog, G. Minimization of Drive Tests in 3GPP Release 11. IEEE Commun. Mag. 2012, 50. [Google Scholar] [CrossRef]
  10. Levanon, N. Lowest GDOP in 2-D scenarios. IEE Proc. Radar 2000, 147, 149–155. [Google Scholar] [CrossRef]
  11. Isa, I.S.; Saad, Z.; Omar, S.; Osman, M.K.; Ahmad, K.A.; Mat Sakim, H.A. Suitable MLP Network Activation Functions for Breast Cancer and Thyroid Disease Detection. In Proceedings of the 2010 IEEE Second International Conference Computational Intelligence, Modeling and Simulation, Bali, Indonesia, 28–30 September 2010. [Google Scholar]
  12. Kasevski, M.; Tiwsin, V.; Pozdukhy, A. Machine Learning for Spatial Environment Data, Theory, Application, and Sofware, 1st ed.; EPFL Press: Lausanne, Switzerland, 2009. [Google Scholar]
  13. Liberti J.C., Jr.; Rappaport, T.S. Smart Antennas for Wireless Communications: IS-95 and Third Generation CDMA Applications; Prentice Hall: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  14. Manolakis, D.G.; Ingle, V.K.; Kogon, S.M. Statistical and Adaptive Signal Processing; McGraw-Hill Book Co.: New York, NY, USA, 2000. [Google Scholar]
  15. Myung, H.G.; Goodman, D.J. Single Carrier FDMA, a New Air Interface for Long Term Evolution; Wiley: New York, NY, USA, 2008. [Google Scholar]
  16. Chen, J.; Huang, Y.; Benesty, J. Audio Signal Processing for Next-Generation Multimedia Communication System; Kluwer Academic Publisher: Boston, MA, USA, 2004; pp. 197–223. [Google Scholar]
  17. Kay, S.M. Fundamentals of Statistical Signal Processing: Estimation Theory; Prentice Hall: Englewood Cliffs, NJ, USA, 1993. [Google Scholar]
Figure 1. Outdoor geo-location scenario between two multi-floor buildings.
Figure 1. Outdoor geo-location scenario between two multi-floor buildings.
Sensors 17 00817 g001
Figure 2. Architecture of the NMLP geometry classifier.
Figure 2. Architecture of the NMLP geometry classifier.
Sensors 17 00817 g002
Figure 3. IGT procedure for generating the training data.
Figure 3. IGT procedure for generating the training data.
Sensors 17 00817 g003
Figure 4. GDOP for different distances.
Figure 4. GDOP for different distances.
Sensors 17 00817 g004
Figure 5. GDOP for different CADs.
Figure 5. GDOP for different CADs.
Sensors 17 00817 g005
Figure 6. GDOP for different distances.
Figure 6. GDOP for different distances.
Sensors 17 00817 g006
Figure 7. MSE of the MLP network using a different number of neurons for the first and second layers.
Figure 7. MSE of the MLP network using a different number of neurons for the first and second layers.
Sensors 17 00817 g007
Figure 8. MSE comparison between NMLP and MLP networks for six by six neurons in the first and second layers.
Figure 8. MSE comparison between NMLP and MLP networks for six by six neurons in the first and second layers.
Sensors 17 00817 g008
Table 1. Symbol table.
Table 1. Symbol table.
e r i , 0 Measurement error
σ e i Standard deviation of the measurement error
w i [ n ] Discrete white Gaussian noise
r i , 0 The differential distance between the MUE and the master FeNB, and the MUE and FeNBi
τ ^ i , 0 Analog time delay between FeNBi and FeNB0
f s Sampling frequency of the analog-to-digital converter
W n m ( s ) ( q ) Weight
μ 0 Learning rate
cLight velocity
GThe variance of the linear optimum unbiased estimator divided by the variance of measurement noise
ε Measurement error
RSphere radius
TTranspose of the matrix
N c Total number of correct classifications
PTotal number of 4 FeNBs within the sphere
Table 2. Coordinates of sixteen FeNBs.
Table 2. Coordinates of sixteen FeNBs.
FeNB1: (0, 0, 30)FeNB5: (0, 0, 40)
FeNB2: (0, 120, 30)FeNB6: (0, 120, 60)
FeNB3: (120, 0, 30)FeNB7: (120, 0, 20)
FeNB4: (120, 120, 30)FeNB8: (120, 120, 40)
FeNB9: (0, 60, 5)FeNB13: (0, 0, 20)
FeNB10: (0, 60, 60)FeNB14: (0, 15, 30)
FeNB11: (120, 0, 0)FeNB15: (0, 30, 40)
FeNB12: (120, 100, 60)FeNB16: (0, 45, 50)
Table 3. Estimated NMLP output for different geometries.
Table 3. Estimated NMLP output for different geometries.
Neuron 1 OutputNeuron 2 OutputNeuron 3 OutputNeuron 4 OutputNeuron 5 Output
square0.94180.01860.05870.01450.0173
lozenge0.05040.89510.07510.02910.0912
IQ0.01240.01860.94270.06240.0174
straight line0.01070.01810.01140.97550.0130
not co-plane0.01780.01360.05840.01750.9581
Table 4. Number of geometry classifications.
Table 4. Number of geometry classifications.
SquareLozengeIQStraight LineNot Co-Plane
all FeNBs1621741593
FeNB 4 not included0213941220
FeNB4,8 not included002635733
Table 5. GDOP of the first three best shapes for all FeNBs.
Table 5. GDOP of the first three best shapes for all FeNBs.
CADSquare (0° CAD) (Distance of MUE to Square Center)Lozenge (0.8° CAD)Lozenge (3.1° CAD)
MUE Trajectory
(60, 0, 0)8.114 (67.082 m)10.90112.648
(60, 10, 0)6.359 (58.31 m)7.5069.796
(60, 20, 0)4.781 (50 m)6.9057.203
(60, 30, 0)3.203 (42.426 m)4.6315.049
(60, 40, 0)1.806 (36.056 m)2.5643.593
(60, 50, 0)1.533 (31.623 m)2.0122.432
(60, 60, 0)1.053 (30 m)1.4351.915
(60, 70, 0)1.546 (31.623 m)2.2162.609
(60, 80, 0)2.371 (36.056 m)2.6513.140
(60, 90, 0)2.689 (42.426 m)3.2063.710
Table 6. GDOP of the first three best shapes not including FeNB4.
Table 6. GDOP of the first three best shapes not including FeNB4.
CADLozenge (0.8° CAD)Lozenge (3.1° CAD)IQ (0 m)
MUE Trajectory
(60, 0, 0)10.90112.64813.703
(60, 10, 0)7.5069.79611.27
(60, 20, 0)6.9057.2039.656
(60, 30, 0)4.6315.0496.157
(60, 40, 0)2.5643.5934.658
(60, 50, 0)2.0122.4323.604
(60, 60, 0)1.4351.9152.161
(60, 70, 0)2.2162.6092.958
(60, 80, 0)2.6513.1403.411
(60, 90, 0)3.2063.7103.907
Table 7. GDOP of the first three best shapes not including FeNB4, FeNB8.
Table 7. GDOP of the first three best shapes not including FeNB4, FeNB8.
DistanceIQ (0 m)IQ (8.75 m)IQ (15.052 m)
MUE Trajectory
(60, 0, 0)15.47416.46617.476
(60, 10, 0)12.86214.91116.143
(60, 20, 0)10.78111.16812.569
(60, 30, 0)6.6567.1147.572
(60, 40, 0)5.0056.4366.632
(60, 50, 0)4.3545.4446.733
(60, 60, 0)3.0154.0646.256
(60, 70, 0)3.4823.9244.719
(60, 80, 0)3.8154.3276.191
(60, 90, 0)3.2525.2355.592

Share and Cite

MDPI and ACS Style

Mar, J.; Chang, T.Y.; Wang, Y.J. A Quadrilateral Geometry Classification Method and Device for Femtocell Positioning Networks. Sensors 2017, 17, 817. https://doi.org/10.3390/s17040817

AMA Style

Mar J, Chang TY, Wang YJ. A Quadrilateral Geometry Classification Method and Device for Femtocell Positioning Networks. Sensors. 2017; 17(4):817. https://doi.org/10.3390/s17040817

Chicago/Turabian Style

Mar, Jeich, Tsung Yu Chang, and Yu Jie Wang. 2017. "A Quadrilateral Geometry Classification Method and Device for Femtocell Positioning Networks" Sensors 17, no. 4: 817. https://doi.org/10.3390/s17040817

APA Style

Mar, J., Chang, T. Y., & Wang, Y. J. (2017). A Quadrilateral Geometry Classification Method and Device for Femtocell Positioning Networks. Sensors, 17(4), 817. https://doi.org/10.3390/s17040817

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