Connected Objects Geo-Localization Based on SS-RSRP of 5G Networks

: The Global Positioning System (GPS) is not the only way to solve connected objects’ geo-localization problems; it is also possible to use the mobile network infrastructure to geo-locate objects connected to the network, using antennas and signals designed for voice and data transfer, such as the 5th generation network. 5G is considered as a least expensive solution because there is no speciﬁc equipment to set up. As long as the object is in an area covered by the network, it connects to the nearest 5G Micro-Cell (MC). Through exchange of signals with the MC node we can locate the object. Currently, this location is very fast with less than 5 s but not very precise because it depends on the number of MC antennas of the operator in question and their distance. This paper presents a novel technique to geo-locate connected object in a covered 5G area. We exploit the 5G SS-RSRP used for signal quality measurement, to estimate the distance between two Connected Objects (COs) in move and in a dense urban area. The overall goal is to present a new concept laying on the 5G SS-RSRP signalling. The proposed solution takes into consideration the Deterministic and the Stochastic effect of the received signals which is not treated by the previous works. The accuracy is optimum even after approaching to the distance of one meter which is not reached in previous works too. Our method can also be deployed in the upcoming 5G network since it relies on 5G signals itself. This work and that of Wang are both based on RSRP and give comparable theoretical complexities therefore comparable theoretical execution times as well. However, to obtain a reliable learning Wang requires a huge amount of data which makes it difﬁcult to get a real time aspect of this algorithm. The use of RSRP and the elimination of the learning phase will give more chance to our work to achieve desired performances. Numerical results show the appropriateness of the proposed algorithms and good location accuracy of around one meter. The Cramer Rao Lower Bound derivations shows the robustness of the proposed estimator and consolidate the work.


Introduction
Global Positioning System (GPS) has for decades been the unique practical way of global positioning. This long period of pioneering work has allowed GPS to acquire major technical improvements in terms of equipment and theory. But the most important benefit of GPS remains linked to the change of human habits by making localization an unavoidable feature in our lives and in different fields such as transports, military, health and games. In spite of this remarkable accomplishment, GPS still has limitations that slow down and keep on slowing down progress in the field of localization. These limitations are mainly related to the accuracy of the satellite signals used or their availability.
Indeed, for security reasons, GPS limits the accuracy of signals used in the civil domains giving a geo-localization of the order of a few meters. The security and the strategic nature of this limitation confirms that this will not be lifted in the coming years.
Location-based applications such as vehicle collision detection, real-time Assisted Living for healthcare providers or visually impaired people localization will continue to be subject of scientific papers until we have improved the accuracy to a few centimeters. However, the next generation of mobile network, 5G, promises a leap in performance technologically compared to 4G. They can offer a flow rate of 10 times higher than the current 4G. The latency is also reduced tremendously since it is divided by 10 compared to 4G and LTE-A. The 5G network ensures robust reliability and radio link availability. The network is also able to simultaneously connect a large number of objects with a better energy efficiency [1]. All those criteria make the 5G a strong candidate to assist GPS features.
Localization solutions based on radio signals have evolved theoretically and fundamentally over the last few decades for the 4G and older generations. Works in 2G and 3G such as [2][3][4] have accuracy over 100 m. Ref. [5] use trilateration techniques in 3G and 4G networks and have achieved accuracies of 50 m. Refs. [6,7] use delay techniques and obtain accuracies of 10 m.
However, the next generation of mobile network, 5G, promises a leap in technology performance compared to 4G and previous generations. They can offer a flow rate of 10 times higher than the current 4G. The latency is also reduced tremendously since it is divided by 10 compared to 4G and LTE-A. The 5G network ensure robust reliability and radio link availability. The network is also capable to simultaneously connect a large number of objects with a better energy efficiency. All those criteria make the 5G a strong candidate to assist GPS features.
The performance promised by the arrival of 5G is estimated in the latest 3GPP standard update to be less than 4 m. This will unlock new opportunities for location-based applications such as vehicle collision detection or precise topographic mapping. Recent research works are trying to improve the accuracy to reach the 3GPP objectives and specially to satisfy the requirements of new applications [8].
The works in [9][10][11] use mmWave and a massive number of antennas, confirming achievable positioning accuracy below one meter. In [12] authors use RSS measurements and distributed massive MIMO and show an accuracy of 35 m from a single base station. The works in [13][14][15] use multipath assisted techniques with advanced tracking algorithms for positioning. Other works such as [16][17][18][19] propose collaborative positioning solutions and localization solutions in ultra-dense networks where the simulation of certain works shows the ability to reach accuracies of less than one meter.
Khan et al. [20] provided an algorithm to predict the location of a person in 5G networks using RSS measurements. For this purpose, the authors use in the first step the relative coordinates by using Isomap and in a second step these coordinates are transformed by a Procrustes analysis. The performance of their algorithm is evaluated using the Cramer-Rao lower bound. Talvitie et al. [21] address high-speed object positioning and tracking in a 5G network. The authors obtain a positioning accuracy of less than 3 m for real data collected on a high-speed train route between Shanghai and Beijing. Authors use an extended Kalman filter to predict the positions.
The objective of this paper is to present an innovative method to exploit the 5G radio network to assist GPS in Geo-localization. The evaluation was carried out on a simulated signal on air and compared to recently published works. Parameters, such as Synchronization Signal-Reference Signal Received Power (SS-RSRP), distance between two COs were estimated. The method was validated together with the suitable distribution law to inform about the accuracy of the precision. In fact, distance is the only exact metric that we can use as a basis for experimental validation of our method. Whether it is the distance between two COs or the distance between two positions of a single COs, this measurement can be obtained easily and with the desired level of accuracy. Knowing D ij this will be a key solution to many problems in industry 4.0, e.g., (inter-car distance detection, inter-robot distance detection and even inter-drone etc ...). The distance D ij will become fully known without the deployment of any sensors or GPS equipment. In 5G concept any single device will be covered immediately by the nearest 5G antenna (MC). Hence the coverage is guaranteed everywhere, certainly in outdoors or even in indoors. Another important application, knowing D ij , will give full knowledge of a realtime constellation of Connected Objects (COs) in move. In this work, following are the four global achievements: • We demonstrated analytically and by simulation that the proposed algorithm gives good accuracy when approaching the 0.55 m barrier. This barrier has never been reached in previous works. • We have full knowledge of the error distribution, as it follows the Rayleigh law. This knowledge will help to tune the solution to other positioning scenarios and applications. • The CRLB is derived to showcase the robustness of the proposed estimator at low and high SNR. Theoretical analysis and simulation both confirm the validity of the proposed model. • A comparative study of complexity is also performed with two recent works Wang in [22] and Kim in [23].

Path Loss Model
In urban dense environment, the propagated radio signals are subject to many electromagnetic phenomena such as, reflection, refraction, scattering and diffraction. Hence the received signal will be guided through multiple paths, and this will create the multi-paths fading subject to destructive or constructive interference. In large scale fading cases, obstacles trigger shadowing effects and create multi-paths and in small scale fading interference will occur and create multi-paths as well.
The behavior of the Radio Signal Strength (RSS) in time and in space, is two folds: deterministic and stochastic process. The degradation of power signal strength is defined as Path Loss PL. Many factors could be used for the estimation of PL, such as multi-paths phenomenon, distance from transmitter, penetration losses through walls, etc. The RSS expression is the combination of deterministic effects and stochastic effects using Gaussian model as where d represents the l 2 -norm distance between the transmitter at (x t , y t ) and the receiver (x r , y r ). PL d 0 is the loss at a reference distance d 0 . The constant n is generally greater than 2 and specific to propagation environment. χ σ is the main stochastic behavior of PL and it could be modeled as zero mean Gaussian random variable with variance σ 2 . However, when the channel is time variant the wide-sense stationary process assumption is lost and the shadowing or the multi-paths fading effects will dominate χ σ [24]. The RSS could be modeled as described by the equation below: From (2) the path loss is considered as large scale fading (very slow process) and the remaining part is considered as small scale fading (fast process) due to multi-paths and shadowing. The main part responsible for shadowing effects is Φ dBm = log 10(Φ mW ) and Φ ∼ N(0, σ 2 Φ ) is a zero mean Gaussian distribution with variance σ 2 Φ . Next to shadowing, multi-paths fading could be also modeled in function of space and time as Ξ dBm (d, t) = 10 log 10(Ξ mW ) and Ξ dBm (d, t) ∼ N(0, σ 2 Ξ ). In many commercial wireless devices the path loss (RSS) can be recoded using the Received Signal Strength Indicator (RSSI) metric. For the rest of this paper, the RSS will be used as RSSI recoded values from wireless devices.

SS-RSRP Model
Synchronization Signal-Reference Signal Received Power (SS-RSRP) is useful for 5G system analysis. It can also shed light on emerging 5G Public Networks (PN) operation and thus help develop local licensing regulations and interference control mechanisms. Generally SS-RSRP is measured through the use of pilots symbols, specific to the deployed cell commonly known as Cell specific Reference Signal (CRS) in 4G networks and altered in 5G by the Synchronization Signal Block (SSB). The accurate measurement of SS-RSRP mainly in micro cells and dense urban area and where interferences become a potential source of signal degradation, makes the accuracy of received SS-RSRP a challenging task. Henceforth, this subsection, will derive analytically the received SS-RSRP based on the full knowledge of RSS signal.
Another approximation of PL is used in literature as where P E is the transmission power of T x antenna, P ca is the cable losses and G R x is the R x antenna gain.
The proposed model will be elaborated based on Figure 1, along with the following assumptions: 1 We will assume in this paper, that there is no cable losses to simplify the calculation process.
The antenna gain depends on the direction of arrival to the Macro-Cell (MC) Currently OFDM is the potential waveform for 5G networks, hence we consider the k th -OFDM symbol in baseband model. The transmitter and receiver are both equipped with one antenna (n r = n t ). The channel type deployed in this paper is a selective Rayleigh channel. The analytical function of the received k th -OFDM symbol is expressed as where u k is the transmitted symbol, w k is zero mean white Gaussian complex noise of variance N 0 2 and L is the number of multipath. Let us define the equivalent channel matrix C (n c +n g )×(n c +n g ) , where n c is the number of subcarrier and n g is the number of subcarrier reserved for the Cycle Prefix (CP) interval. Thus one can rewrite the k th -OFDM received signal in matrix notation as y k = G k u k + w k , where y k ∈ C (n c +n g )×1 and w k ∈ C (n c +n g )×1 represents additive white Gaussian noise channel at time k. When removing CP, the frequency domain signal is generated again by a DFT process and the whole mathematical expression is defined as is the cycle prefix removing matrix for the first n g elements of y k . This will lead to another expression of the DFT-output as z k = (F ⊗ I n t )ξ 3 (F −1 ⊗ I n t ) x k + w k = H k x k + w k , ξ 2 = [0 n c n g I n g ] ⊗ I n g . The matrix ξ 3 ∈ C n c ×n c is the circulant matrix, and ξ 3 = ξ 2 G k ξ 1 where w k is the noise in the frequency domain with zero mean and variance σ 2 w and H k is C n c ×n c frequency domain matrix defined as H k = (F ⊗ I n c )ξ 3 (F −1 ⊗ I n c ) . Now we can express analytically the channel matrix in frequency domain H k as [26,27]: where nth block H k (n) represents the OFDM channel gain at the nth subcarrier and can be written as The received frequency-domain signal z n,k at the DFT output can be expressed as In 5G NR the Synchronization Signal Block SSB is used for signal measurement, therefore N special sub-carriers dedicated to SSB measurement are deployed. The SS-RSRP signal is defined as average of channel power during one measurement period of (SSB). Hence The frequency domain channel H n,k is a diagonal , with the following expression where P h = E{h p k (l) h m k−k (l ) * } denotes the mean power to each channel coefficient in time domain. On the other side, the mathematical expression of RSSI at the receiver could be interpreted as: Henceforth, the SS-RSRP in dBm could be re-written as SS RSRP [dBm] = RSSI[dBm] + 10 log 10(n c ) − 10 log(N) When combining (2) and (12), we can derive the expression of the angleθ i and the frontal distance with the MC aŝ Finally, when we refer to Figure 1 the distance between two COs in 5G cell will be estimated as

Distance Estimation
The initial phase of the proposed solution is to simulate a realistic RSS environment using Matlab. The model setup deploys the log-normal shadowing model as presented in (2). All parameters are set using empirical observations and Gaussian noises Φ dBm (d) and Ξ dBm (d, t). The MC position is supposed well known fixed at (0,0) and any estimated distance will be calculated based on (l 2 ) norm. Each environment has its own empirical path loss constants n as in (2). Hence it will be crucial to extrapolate correct n values as described. The n value is obtained by setting the derivative of the following mean square error estimate J(n) = Σ k i=1 (−PL d 0 − χ σ − 10n log 10( d i d 0 )) 2 to zero, ∂J(n) ∂n = 0. The Algorithm 1 describes the initial phase for a realistic RSS environment, measured at the receive antenna, where P t denotes the transmit power (dBm) of the transmitter, d 0 is the reference distance d 0 (m), Pd 0 is the received power at the distance d 0 and PL d 0 = P t − Pd 0 is the path loss at the reference distance. The Algorithm 1 will serve to generate a realistic RSS value to be used in the next estimation phase as described in Algorithm 2. The Algorithm 2 summarizes the needed steps to estimate the distance between two COs. The COs are supposed to be in move while deploying a 5G NR radio access network. rss(x, y) =P t − PL d0 − 10n log 10(  (12) 4 Angle estimation as in (14) 5 Frontal distance estimation as in (14) (12) 9:θ j = ( 70π Angle estimation as in (14) 10 Frontal distance estimation as in (14) 11: end for 12: for i in 1 to q do 13: for j in 1 to p do 14: Mutual distance estimation as in (16) 15: end for 16: end for

Numerical Results and Validation
We consider a 5G NR transmission system whose parameters, taken by 5G Release 15 [28], are listed in Table 1. To obtain all the numerical results shown in the next of this section, 10,000 independent realizations were simulated. Perfect knowledge of channel state has been supposed in all our performance evaluation trials.

MSE Method
We used Matlab software Release 2019 Ubuntu version to compile Algorithms 1 and 2. Initially a realistic radio signal strength environment is generated using Algorithm 1. The two COs are supposed to be mobile within one meter square with a full coverage ensured by an MC. Afterward Algorithm 2 is enabled to generate the estimated distance between the two COs. We illustrated the whole process as in Figure 2.

Distribution law Observation
Step : 0 Step : 1 Step : 2 Step In order to scale the accuracy of the proposed Algorithm, we used the Minimum Mean Square Error (MSE). As it is illustrated in Figure 3 the MSE tends to zeros at the very first rounds of Algorithm 2. We notice that the convergence time for Algorithm 2 is not immediate, this is due to the needed calculation time required for RSSI generation in Algorithm 1 and also the processing time needed for Algorithm 2. Henceforth the MSE values become zeros after a short transitional phase.

Accuracy Quantification
We can quantify the accuracy of our method as the deviation in horizontal and vertical coordinates (x, y) compared to the estimated coordinates (x,ỹ). Q = (x −x, y −ỹ) = (X, Y). Normal distribution appears to be the most convenient choice, especially in the light of Figure 4. One can consider the deviations at x and y are independent random variables and follow Gaussian laws of zero mean and variance σ 2 = 0.1758 m. The intervals next to the parameter are the 95% confidence intervals for the distribution parameters. These values are collected after Matlab processing of Figure 4. The density of joint probability of deviations (X, Y) is given by: The characterization of the precision in x and in y does not give information about the probability of the distance between the desired point and the point obtained. To obtain this probability, we change Cartesian coordinates to polar coordinates according to the transformation: Thus the Jacobian of transformation is r: This will lead to : Hence to find the marginal, it is necessary to integrate on the other random variable, that is: The angle is therefore uniformly distributed over [0, 2π]. On the other hand, we deduce: The probability density followed by the distance R is that of a Rayleigh variable as illustrated in Figure 5. As it is shown in Figure 5, the precision of the proposed method is robust when we are approaching the barrier of one meter. The precision becomes weak while approaching the origin. At distance greater than 0.55 m the error probability tends to the zero.

Complexity and Performance Analysis
Our solution uses 2 algorithms. The theoretical complexity depends only on the second one since the first one is executed off-line and the complexity of our solution will be O(n 2 ). The Algorithm described in [22] is based on SS-RSRP measurements and a step using Random forest classifier with the machine learning Library Sklearn. We are interested here in the Algorithm without considering the learning part which is also executed off-line. The theoretical complexity of this Algorithm is O(n 2 ), Kim et al. in [23] proposed a multiplemodel PHD filter and map fusion for cooperative positioning and mapping in vehicular networks with 5G mmWave communication links. The proposed framework is based on 3 algorithms. The Map correction Algorithm uses repetitively and successively the two other algorithms. The theoretical complexity of these algorithms is O(n 3 ). Presented works in Table 2 are all compliant with 5G standards but with different complexities. The theoretical complexity of Kim et al. is higher than the other two algorithms and therefore the theoretical execution time is higher. Our work and that of Wang are both based on SS-RSRP and give comparable theoretical complexities therefore comparable theoretical execution times. However, to obtain a reliable learning Wang's method requires a huge amount of data which makes it difficult to be used in real time. The use of the SS-RSRP and the absence of a learning phase are the main advantages of our work that allowed us to achieve the desired performance.

CRLB Analysis
The Cramer-Rao Bound (CRB) is a second performance criterion that gives a lower bound to the mean square error of estimation in the set of unbiased estimates. It is used here as a performance measure for the proposed estimator. The observed signal is described by (16). It can be written as where w(n) is a white Gaussian noise following N(0, σ 2 ). σ 2 denotes the variance of the white Gaussian noise. In this equation we have four unknown variables To derive the CRLB, first, we write the likelihood function as Since we have The Fisher Information matrix I(Ψ) doses exist and the CRLB = 1 I(Ψ) . In fact After some mathematical manipulations we find We assume that the sum of sines or cosines becomes negligible for large N. The CRLB gives the lowest variance of our proposed estimator D ij . Figure 6 shows that the proposed estimator follows at all SNR ranges the CRLB with almost the same gap. At one meter we notice a gap of 9 dB, this is due to the mathematical approximation done to derive the CRLB. When the SNR values increase asymptotically the variance of the estimator becomes closer to the CRLB limit.

Conclusions
In this paper we proposed an innovative method to estimate the distance between two mobile COs from the measurements of SS-RSRP when 5G is deployed. We have provided a sophisticated analysis, naturally assuming that the SS-RSRP values depend on the distance from the COs to the MC. In particular, it has been shown that the use of these 5G signals makes it possible to provide good location accuracy without resorting to the GPS system. Using simulations and a detailed analytical study, it has been shown that the proposed algorithms have achieve good location accuracy of distance of around one meter.  Data Availability Statement: The study did not report any data.

Conflicts of Interest:
The authors declare no conflict of interest.