5G Positioning: An Analysis of Early Datasets

Global Navigation Satellite Systems (GNSSs) are nowadays the prevailing technology for positioning and navigation. However, with the roll-out of 5G technology, there is a shift towards ‘hybrid positioning’: indeed, 5G time-of-arrival (ToA) measurements can provide additional ranging for positioning, especially in environments where few GNSS satellites are visible. This work reports a preliminary analysis, the processing, and the results of field measurements collected as part of the GINTO5G project funded by ESA’s EGEP programme. The data used in this project were shared by the European Space Agency (ESA) with the DICA of Politecnico di Milano as part of a collaboration within the ESALab@PoliMi research framework established in 2022 between the two organizations. The ToA data were collected during a real-world measurement campaign and they cover a wide range of user environments, such as indoor areas, outdoor open sky, and outdoor obstructed scenarios. Within the test area, eleven self-made replica 5G base stations were set up. A trolley, carrying a self-made 5G receiver and a data storage unit, was moved along predefined trajectories; the trolley’s accurate trajectories were determined by a total station, which provided benchmark positions. In the present work, the 5G data are processed using the least squares method, testing and comparing different strategies. Therefore, the primary goal is to evaluate algorithms for position determination of a user based on 5G observations, and to empirically assess their accuracy. The results obtained are promising, with positional accuracy ranging from decimeters to a few meters in the worst cases.


Introduction
Global Navigation Satellite Systems (GNSSs) have been the dominant technology for positioning, navigation, and timing for decades [1,2].Our smartphones seamlessly fuse GNSS data with readings from other sensors to improve position estimation accuracy; the deployment of 5G technology across the world is also expected to open new opportunities in the domain of positioning, not just communication.Firstly, 5G networks can broadcast GNSS corrections for improved accuracy.To continue, 5G can augment GNSSs by providing additional range measurements in environments where only a few satellites are visible.The accurate determination of location information is a strategic byproduct of 5G cellular networks [3][4][5], as it has the potential to benefit numerous commercial applications, from individual to public services.This covers application domains such as transportation, public safety, retail, and healthcare.The integration with GNSSs [6][7][8] becomes relevant in urban areas, where reliance on GNSSs alone could be challenging [9].Compared to previous mobile generations such as LTE, 5G technology features a new radio-access technology called new radio (NR), which offers several advantages for precise positioning.Designed by the 3rd Generation Partnership Project (3GPP), 5G NR is engineered to meet a variety of performance metrics, serving enhanced mobile broadband (eMBB), ultrareliable low-latency communication (URLLC), and massive machine-type communication (mMTC) applications [10].Next, we present the most important key enablers for accurate positioning.It is important to realize that while this list is not directly pertinent to the following discussion of our experiment, it could be useful as a summarized introduction to the technological aspects.
New positioning reference signals: to enable more accurate positioning than LTE with release 16 [11], the 3GPP updated the positioning reference signal (PRS) and sounding reference signal (SRS), providing the downlink (DL) and uplink (UL) signals, respectively [12].Traditional signals such as channel state information reference signals (CSI-RSs) and synchronization signals (SSs) have limitations that make them less suitable for positioning: interference caused by adjacent cells reduces their ability to detect a sufficient number of neighboring 5G base stations; furthermore, the signals of neighboring cells overshadow those from distant ones (near-field effect), hindering the detection of the latter.The new signals are specifically designed to improve audibility thanks to the muting concept [13]: multiple base stations can transmit PRS in a coordinated manner by literally 'muting' the less relevant PRS transmissions to avoid interference.
New positioning methods: 5G can provide different measurements that are related to the position of the user equipment (UE).In summary, these measurements can be based either on angles or on distances [14,15].To elaborate, angular measurements include the UL's angle of arrival (AoA) and DL's angle of departure (AoD), distance-based measurements include time of arrival (ToA), DL and UL time difference of arrival (TDoA), round trip time (RTT) and multi-cell round trip time (MC-RTT).The positioning methods (or estimation algorithms) may be implemented in one or more of the following modes: UE-assisted mode, UE-based mode, stand-alone mode, and network-based mode [14].
New spectrum: Frequency Range 2 (FR2) is an ultrahigh-frequency band allocated for 5G in the millimetre-wave (mmWave) region, spanning from 24.3 GHz to 52.6 GHz.This portion of the spectrum complements the existing Frequency Range 1 (FR1), which covers the bands below 6 GHz and recently expanded to 7 GHz.It addresses the issue of spectrum scarcity in wireless communication systems and enables high data rates, capacity, and bandwidth with low latency, while also providing superior positioning accuracy [3,16].In this frequency range, radio signals experience penetration and diffraction losses, resulting in a dominant line-of-sight (LOS) component and minimal multipath effects [17].However, the use of millimeter-wave wireless signals also poses challenges, such as high path loss, which can be mitigated through the adoption of new specialized compensation techniques, such as beamforming and highly directional antennas.
Wideband carriers: NR provides a significant bandwidth improvement over LTE; while LTE provides a maximum of 20 MHz, NR provides up to 100 MHz in Frequency Range 1 (FR1: 450 MHz to 7 GHz) and 400 MHz in Frequency Range 2 [18].The variance in the delay estimation is inversely proportional to the bandwidth of the signals; this implies that as the signal bandwidth widens, the uncertainty in the delay estimation decreases as the main lobe of the correlation function becomes narrower.Narrower lobes are more easily discriminated by the receiver and, therefore, the distinction between direct and reflected paths improves [19].
Massive multiple-input multiple-output (MIMO) and beamforming: massive MIMO wireless communication refers to the idea of equipping cellular base stations with very large quantities of antennas.Such a massive number of antennas causes interference problems, which can be mitigated by deploying the beamforming technique: a process designed to produce the radiated beam patterns of the antennas by building up the processed signals in the direction of the desired terminals and canceling beams of interfering signals [20].The benefits of using beamforming in massive MIMO systems include enhanced energy efficiency, improved spectral efficiency, increased system security, and applicability for mmWave bands.
The main purpose of this paper is to implement the positioning using 5G time-ofarrival (ToA) observations and to assess the accuracy and the reliability of the estimated coordinates from experimental data.These real-world 5G ranging measurements have been provided by the ESA and were produced during the GINTO5G project, described in the next section.Different strategies to process these data will be compared.In general, when time series of observations are processed to estimate trajectories, a proper filtering process like Kalman filtering [21,22], allows the smoothing of the noise and blunders in individual epochs.On the contrary, we process the epochs independently, using the Least Squares method [23], because we are interested in the accuracy and reliability of the single-epoch solutions, not in the filtering results.In the rapidly evolving field of 5G positioning, the present research introduces an innovative approach in the use of the sounding reference signal.Distinct from the 3GPP specifications for 5G positioning [24], our methodology makes use of the SRS in the downlink direction and proposes signal design changes that would allow it (see Section 2.4 for details).Furthermore, most of the previous research studies on 5G positioning, and in particular positioning using the SRS signal [12,25], are based on simulations rather than real-world experiments carried out in a controlled measurement setup.This setup generates at least two major benefits: repeatability of measurement conditions, and more tangible results than purely theoretical studies.The paper is structured as follows: Section 2 describes the GINTO5G experiment and data collection; Section 3 describes in detail the strategies applied to process the data; Section 4 discusses the results; and finally, Section 5 presents the conclusions.

The GINTO5G Experiment
As part of ESA's GINTO5G project, a number of field tests were carried out in 2021 with the aim of collecting and analyzing 5G signals with the specific aim of extracting pseudoranges to estimate the location of the device being tested.For the 5G experiment, an adaptation of the sounding reference signals (SRSs) [24], configured on the downlink instead of the uplink transmission, has been used.This signal was transmitted on 3.7 GHz, spanned 100 MHz in bandwidth, and consisted of a custom frame structure (further details are in Section 2.4).To check the trajectories estimated by 5G the benchmark trajectory was estimated by a total station.It should be specified that although the survey was performed in 2021, this paper summarizes [26] and reports the first experimental results.Technical details of the GINTO5G experiment are described in [27][28][29][30][31] and are just briefly summarized here.

5G Testbed
The GINTO5G experiment, a collaborative initiative of the ESA with several European companies and universities, represents an important step in validating the potential of 5G for positioning and navigation.A designated testing area was selected to host the specialized infrastructure needed for the experimental 5G signal transmission and reception.This setup is unique in that it is not a conventional operational setup.Instead, it is a selfmade engineering construct of a compact 5G private network, tailored specifically for the scope of the GINTO5G project, namely, positioning with 5G signals.The testbed is located at the Fraunhofer Institute in Germany, close to Nürnberg.The measurement setup consists of a transmitting block and a receiving block.The transmitting block is composed of eleven transmitters, paired with an IT ELITE Antenna SEC3710 DP, that emulate a mini 5G private network able to transmit 5G positioning signals.The receiving block consists of a USRP X300 with an omnidirectional antenna programmed to receive 5G signals from the transmitters.Moreover, a Leica MS50 total station is used to survey the ground truth of each trajectory to be estimated.

Experimental Campaign
All measurements were conducted at the L.I.N.K. Test and Evaluation Center of the Fraunhofer IIS campus in Nürnberg (Germany).
As shown in Figure 1, eleven transmitters emulating 5G transmitting reference points (TRPs) were distributed across the indoor and outdoor area (loading zone + driveway) and their position was surveyed accurately using a total station prior to the beginning of the 5G data acquisition tests.A trolley was used to host the receiving unit, able to acquire 5G signals, during its movement on a-priori defined trajectories to different areas of the campus.The trolley, shown in Figure 1, was also equipped with a battery backup system and a PC to log the measurements during the field trials.

Execution
The total number of trajectories taken into consideration for this experiment can be clustered into three sub-groups depending on the location of the measurement: loading zone, driveway, or indoor area.

Data Recording
This section describes the necessary steps needed to be performed for processing the raw data and, therefore, obtaining ToA measurements [27,32,33].The receiving unit recorded IQ samples generated using the Matlab 5G toolbox and transmitted by the transmitting points.
The computation of 5G NR ToA is performed based on the SRSs sent on the downlink channels.This represents an intentional deviation from the 3GPP specifications, which define the use of the SRS on the uplink channel.The choice to adapt the transmission of the SRS to the downlink stems from the fact that existing tools and signal post-processing algorithms were developed for the SRS signal.The adaptation to the SRS signal consisted in prefixing each SRS symbol with a unique secondary synchronization signal (SSS) sequence; this SSS sequence helps to match a specific SRS burst to a transmitting antenna.The SRS signal is transmitted as 10 ms frames, as described in the 3GPP specifications, and its configuration is characterized by a 1 OFDM (orthogonal frequency-division multiplexing) symbol and comb2 transmission pattern with the following bandwidth configuration parameters (Table 1): Table 1.Bandwidth configuration parameters, as specified in [24].

Parameter
Value This received signal was processed offline as follows: first, the SSS correlation step looked for the maximum SSS peak to allow for different transmitting antennas to be distinguished; followed by SRS correlation based on a fast Fourier transform approach.The maximum SRS correlation peak represents the ToA for the signal with respect to the sampling rate.For each 10 ms frame, one ToA value was generated.
Although 5G NR transmitters are in general synchronized by a common 10 MHz clock source, the synchronization is not perfect.There is a delay caused by a number of factors, including different cable lengths and the connection between the antenna and ADC (analog to digital converter), and it needs to be compensated for.Furthermore, additional, nonconstant latency can be introduced by each USRP on start-up, necessitating the performance of latency compensation at least once each time the setup is turned on.The following equation gives the effect of latency on the ToA measurement: where τ is the time of flight of the signal from the base station, and dt R and dt BS are, respectively, the clock offset of the receiver and of the base station with respect to a given reference time scale.In the GINTO5G experiment, all the base stations were electronically synchronized and the residual synchronization error resulted in post-set-up calibration with a magnitude of a few meters, which was used to correct the measured ToA.This instrumental calibration is not discussed here because it was preliminary to this work [31].
As a result, the ToA contains just the geometry and the clock offset of the receiver.As further detailed in Section 3, the receiver clock offset was removed by time differencing pairs of ToA measurements, leading, therefore, to TDoA measurements.The calibrated observations contain errors that can be ascribed to three different sources: • clock jitter; • radio channel effects (e.g., multipath); • the accuracy of the reference measurements of the distances between antennas.
Although clock jitter has a general impact on the ToA quality, it cannot be directly compensated for, and in general, low-jitter clocks (picosecond jitter) have to be taken into account during system design.The equipment and conditions used in this work were experimental and well controlled, a situation difficult to reproduce on a standard massmarket device: picosecond-jitter clocks indeed suggest that such results cannot currently be reproduced with mass-market devices/smartphones.This clearly will be an object of future investigations.The coordinates of the base stations and the reference trajectories are given by a terrestrial survey with a Leica MS50 total station: considering the used instrumentation and the surveying technique, we can assume that the accuracy of the benchmark results is at the centimeter level, at least a magnitude better than the FR1 accuracy.As a result, the main error sources in the measurements are channel effects like multipath, and accordingly measuring line-of-sight signals is crucial for a low-error calibration.

Methods
The section focuses on a description of the algorithms implemented to process the ToA observations.The available GINTO5G datasets include 5G time of arrival (ToA) observations (Figure 5, left), their signal-to-noise ratio (SNR) values, local Cartesian coordinates of the base stations, and the ground truth output from the total station.In the processing, the time difference of arrival (TDoA) values will be utilized as the observations.As specified in Section 2, the underlying assumption is that all the base stations are synchronized; hence their clock offsets are set to zero.The clock term of the receiver can be eliminated by differentiating two simultaneous ToA observations.In fact, given the receiver R and the two ToA to two base stations choosing BSi as the reference station, the relevant TDoA is given by Multiplying by the signal velocity c, the observed ToA can be converted to metric pseudoranges: where ρ R is the geometric distance between the base station and the user receiver.Note that the differences of the metric observations, c • TDoA i,j , are directly differences of distances.From a geometrical point of view, the pre-elimination of the clock in the TDoA values does not present any specific advantage with respect to the processing of undifferenced observations (see, for example, [34]); indeed, it is counter balanced by the reduction in the input observations.However, we decided to process the TDoA to investigate this specific approach, and also because it is very popular in the technical literature.Given these premises, our objectives in this study are: 1.
Investigate possible strategies for the choice of the reference station in the TDoA (Section 3.1); 3.
Estimate the positions using the least squares method in single epochs (Section 3.2); 4.
Analyze the results and the statistics of the estimated trajectories with respect to the ground truth (Section 4).
All the base stations of the experiment are at the same height, which is also approximately the height of the trolley.As was preliminarily discussed in the technical documentation of GINTO5G, with such a configuration, the estimation of the vertical coordinate causes an ill conditioning of the resulting system; therefore, only horizontal coordinates will be estimated.

TDoA Analysis
In order to verify the accuracy and reliability of the 5G observations, the measured TDoA are compared with the reference TDoA obtained from the known distances between the base stations (BSs) and the receiver, as determined by the coordinates provided by the total station benchmark TDoA).
In each trajectory, different approaches are tested for choosing the reference station for TDoA computation: • use one BS as reference for all the epochs; • use the pivot method, e.g., following a scheme BS1-BS2, BS2-BS3, and so on, for each epoch; • choose as reference station for each epoch the station with the best SNR: in this case, the configuration can change between epochs.
It is worth noting that, in pivoting, the selection of a specific scheme is not important.Indeed, by a proper propagation of the covariance matrix (as in Equation ( 10)), the solution of the least squares method for different pivoting strategies should be identical up to numerical rounding [34].The results of this analysis are discussed in Section 4.1.

Least Squares Algorithm
In the single-epoch solution, the unknowns are the horizontal coordinates of the receiver [x R , y R ], while the height is kept constant.
Since the least squares method is applicable only to a set of linear equations, the processing involves the linearization of the observation equations with respect to the receiver's coordinates.The starting point of the process is the general formula for the distance between BS and R: where ρBS R and ẽBS R are, respectively, the approximated distance and the unitary vector from BS to R; dx R is the difference between the true and the approximated coordinates of R. Considering two base stations (BS1 and BS2), and remembering that the TDoA will be the input of the algorithm, the difference in the distances can be written as where ρ BS1,BS2 R is the difference in the two distances, from the two BS to R. ∆e BS1,BS2 R is the difference of the two unitary vectors, from BS1 and BS2 to R.
Given all the available n ToA in one epoch, n − 1 TDoA can be derived.A system can be written introducing the classical notation of the least squares method: where y is the vector of the observables, b is the vector containing the known terms, and A is the design matrix, whose n − 1 rows are simply the transpose of ∆e BSi,BSj R .The final LS solution is given by where Q is the cofactor matrix of TDoA, y 0 are the available observations, and N = A T Q −1 A.
In addition to the unknowns, their covariance matrix can be estimated: where σ2 is the a posteriori variance.Even the observations' residuals and normalized residuals can be computed according to [23].
Q is obtained as follows, given the cofactor matrix C of the undifferenced ToA: where ∆ is the matrix of the linear transformation from ToA to TDoA.In our specific case, arranging the stations in order to have the reference one in the first position: We assume that the available ToAs are not correlated: therefore, when weights are applied, C is defined as where the weights can be, for example, given by In case weights are not applied, each η = 1.The weighting of the ToA will be discussed in the following sections.Figure 5 (right) shows the LS processing scheme.
Linearized least squares requires iterations: in our application, the initial approximate coordinates are set in the barycenter of the base stations, at each following iteration the estimated values at the previous iteration are used.The algorithm iterates up to a maximum number of 20 iterations, unless the convergence (set to one millimeter) is reached.
To assess the accuracy of the solutions, for each epoch the estimated position is compared with the ground truth: where t is the epoch, [ x5G , ŷ5G ] is the estimated position, [x TS , y TS ] is the position from the total station.The average and standard deviation of the errors for x and y as well as the 2D plane are then calculated for each trajectory: where N is the number of epochs.This analysis helps to identify any patterns in the errors and provides a better understanding of the accuracy of the 5G trajectory.Least squares convergence and the final results are discussed in Section 4.2.For our data, no correlation exists between the SNR values and the distances between the stations and the receiver.This is due to the small range of distances and covered area.The objective of the preliminary SNR investigation is to search for a criterion to select the optimal reference station for the TDoA.In the following, we assess the influence of the SNR of the reference station on the quality of the data.For each take, the following scenarios are considered for the TDoA computation: • one reference station for the whole trajectory (this choice is repeated for each one station); • the pivot method, with the schema: BS7-BS8, BS8-BS9, BS9-BS10, BS10-BS11; • at each epoch the reference station is the station with the best SNR.
For all the above cases we compute the TDoA individual errors, as the differences between the observed TDoA and the true single differences of the distances between the stations and the receiver.Table 3 shows the statistics on the TDoA errors for each case.In the analysis of the errors, a few observations of some stations are clearly affected by isolated blunders of more than 100 m.These observations, less than ten, are excluded by the following analysis and least squares solutions.Table 3.Take 01: errors (measured minus ground truth) in TDoA for different choices of reference station.Ref. i: base station i as reference for all the epochs; Pivot: pivot scheme; Best SNR: for each epoch the best SNR at each epoch.Mean: mean; St. dev.: standard deviation.Cum.err.: cumulative error, square root of the sum of the squares of the errors divided by the population.
[10] Ref. [ In our data, no correlation exists between the SNR of the reference station and the resulting TDoA errors; therefore, the SNR sorting is not an optimal criterion to select the best reference station in TDoA computation.The above conclusions, discussed in detail, are relevant to Take 01, but all the other takes provide similar results.

Least Squares Solution
In this section, we will present the positions and trajectories estimated using the least squares method in a single epoch utilizing the TDoAs as inputs.The presentation of the results will start with the trajectories in the loading zone, followed by the indoor trajectories, and conclude with the trajectories in the driveway zone.Taking into account the results of the previous section, no SNR weighting is applied to the observations.Moreover, each trajectory is estimated by using one reference station for all the epochs.In the least squares iterations, we anticipate that in all the processing the convergence is reached in less than five iterations.

Loading Zone Trajectories: Takes 01, 02, and 03
In the loading zone, only the outdoor stations (7,8,9,10,11) are used.In this area, a good line of sight (LOS) is expected.The trajectory of Take 01 is analyzed first.The trajectory is repeatedly estimated by using as the TDoA reference, firstly, station 7, then 8, and so on, up to 10.By a visual check of the results and the residuals, the observations of station 7 are clearly affected by the largest errors, both locally correlated and sparse outliers.This is probably due to the presence of a reflecting surface near the station, and the resulting multipath.By excluding station 7 from the dataset, the results clearly improve, due to a significant reduction in spurious patterns and isolated outliers; compare, for example, Figures 9 and 10.Considering that, in general, station 7 showed the best SNR, this confirms that in our data the SNR does not provide a useful index of quality.The final statistics of Take 01 excluding station 7 are shown in Table 4 for different choices of reference stations; the results are almost homogeneous.

2D Error
Ref.Note that the exclusion of station 7 locally improves the quality of the results but does not significantly change the general statistics.After removing station 7, the redundancy of the the least squares method estimates in single epochs is very low (one, equal to three TDoA observations minus two unknowns).4.2.2.Indoor Trajectories: Takes 06, 07, and 08 For indoors, we consider the six indoor base stations (1,2,3,4,5,6).In this case, different choices of the reference station were also tested; however, none of them significantly outperforms the others in term of the error statistics.Therefore, to avoid repetition, we present just the results obtained using station 1 as the reference (Figures 13-15 and Table 6).Take 08 exhibits excellent results, while Take 06 shows a slightly worse performance, and Take 07 provides the worst results.In any case, the indoor results are generally good considering all the disturbances, mainly interferences and multipath, that can occur in such scenarios.Blunders are present in single epochs; as premised in the introduction, they could be smoothed using time-series filtering, but that is outside the scope of the present paper.Table 7 presents an analysis of the residuals obtained using the least squares method for all the trajectories.The observations of TDoAs BS1-BS2 provide the worst results: this is particularly significant in Take 06, while in the other takes the residuals are more homogeneous.As for the loading zone, we remove the data of the worst station, in this case station 2. Tables 8 and 9 report the statistics of Take 06, which improve significantly; the other takes do not significantly improve.In the driveway processing, only the outdoor base stations were used.The driveway trajectories start at the loading zone, where LOS is available for all the stations, then the trajectories move to NLOS conditions.
The positioning (Figures 16 and 17) along the north-west part of the trajectories suffers from shadowing; once the trolley moves closer to the building and enters a non-line-ofsight (NLOS) condition the solutions are completely blinded and meaningless.As the usual statistics of errors (mean, standard deviation, . . . ) are useless, Figure 18 displays the percentages of the errors' magnitudes.The results are very bad and not usable in practice, but they were expected because of the NLOS conditions of this scenario.

395
This experiment presented the opportunity to accomplish the following aims:

Conclusions
This experiment provided the opportunity to accomplish the following aims: • to implement algorithms for positioning using 5G observations that were applied to 5G data from the GINTO5G experiment; • to experimentally assess the accuracy of the positioning in an environment where the deployment of the base stations was carefully controlled and optimized; • to conduct experimental research on techniques aimed at identifying and reducing measurement errors.
The processing focused on TDoA observations, processed in a single epoch using the least squares method.The accuracies of the estimated positions reach the decimeter and the meter levels, respectively, in outdoor and indoor scenarios.Within this specific experiment, the SNRs of the base stations do not exhibit a significant correlation with distance with respect to the receiver; moreover, the use of the SNR as a criterion to choose the reference station in TDoA does not improve the quality of the differenced observations.This can be explained by the fact that the GINTO5G experimental area was relatively small, spanning 40 m × 40 m; moreover, the data of the experiment were calibrated in order to remove the clock biases of the base stations, and this explains the good quality of the results.In our scenario, the redundancies of the single-epoch solutions are small.Some outliers exist when obstructions and reflecting surfaces are present, especially indoors.Clearly, with larger redundancies in single epochs, or by time-series filtering, the estimates could improve.In any case, the obtained results in LOS have accuracies that range from decimeters to meters.
Given these experimental premises, the achieved results are satisfying and promising for further research about mass-market applications.In conclusion, the hybridization of 5G and GNSS signals for positioning holds great promise for a wide range of applications, from autonomous vehicles and augmented reality to smart cities and beyond.
been obtained for Takes 02 and 03; the exclusion of station 7 improves the statistics in these cases as well.The results of using least squares on Takes 02 and 03 are presented in Figures11 and 12and Table5.Note that in Take 02 the worst results are in the same area as those of Take 01.

Figure 17 .
Figure 17.Take 05, reference station 4: Least squares results.Positioning (Figures16 and 17) along the north-west part of the trajectories suffers 389 from shadowing: once the trolley moves closer to the building and enters in a Non-Line of 390 Sight (NLOS) condition the solutions are completely blindered and meaningless.As the 391 usual statistics of errors (mean, standard deviation,...) are useless, Figure18displays the 392 percentages of error's magnitudes: the results are really bad, and not usable in practice, but 393 they were expected because of NLOS conditions of this scenario.

Figure 18 .
Figure 18.Pie charts representing the error magnitude in Takes 04 and 05. 396

Figure 18 .
Figure 18.Pie charts representing the error magnitudes in Takes 04 and 05.

Table 2 .
Take 01: statistics of SNRs (dB) of the base stations.

Table 4 .
Take 01: statistics of 2D errors for each reference station after excluding station 7.