A Novel TDoA-Based Method for 3D Combined Localization Techniques Using an Ultra-Wideband Phase Wrapping-Impaired Switched Beam Antenna

: This paper presents a novel Time Difference of Arrival-based approach suitable for single-anchor positioning systems, implemented by phase wrapping-impaired array antenna. With the latter being a typical occurrence in large Switched Beam Antenna (SBA) operating in the low microwave range. The proposed method takes advantage from the large bandwidth of radio link, established between the anchor and the positioning target, by generating an unambiguous equivalent phase relationship between antenna array elements. The technique is validated by adopting a relatively large SBA antenna operating in the 4.75–6.25 GHz bandwidth, and capable to position a target in a 3D domain. Experimental data, carried out in the 4–7 GHz frequency bandwidth, show that by dealing properly with the inherent constraint of phase wrapping issues, it is possible to get a signiﬁcant improvement on the elevation angle with respect to methods not capable to deal with phase reconstruction and thus operating in a phase-less context. Combining range and angle errors, the associated cumulative distribution function error in 90% of cases shows an error of 0.13 m. Author Contributions: Conceptualization, M.P. and A.C.; methodology, M.P.; software, M.P. and E.C.; validation, A.C. and G.C.; formal M.P.; investigation, M.P.; resources, G.C.; data processing, M.P.; writing—original M.P. and E.C.; and


Introduction
Wireless positioning technologies enable a wide range of application scenarios where positional information provides key points for the development of unattended and noninvasive context-aware services [1,2]. In recent years, the emerging of paradigms that can benefit from localization in many applications, such as Internet of Things and Smart Environment [3], led to an increasing interest in both the scientific community and industry. As a consequence, such new scenarios have fostered different localization techniques and the investigation of related technologies [4].
Several methods are available in the literature, due to the very different cases of studies. They can be classified according to the specific signal features involved in the position estimation, such as Direction of Arrival (DoA), Time of Arrival (ToA), Time Difference of Arrival (TDoA), and range evaluation by acquiring the received signal strength [5]. These techniques rely on the established radio links between the target to be positioned and two or more spacial references, also called anchors. In this respect, different antenna structures and methodologies were proposed during the last years. In particular, the Switched Beam Antenna (SBA) technology got notoriety due to the optimal trade-off between simplicity and effectiveness [6][7][8][9]. Such architecture allows to achieve the necessary spatial diversity to properly evaluate received signal features and perform estimation methods.
The main drawback for SBAs is due to their structural dimensions, because the high directivity requirement for each antenna element leads to significantly large overall dimensions. In turn, it is common that high directivity SBAs have begun to face phase impairment problems caused by increased distance between antenna center of phases [10].
Such issue is well known in Multiple Input Multiple Output (MIMO) antenna arrays design [11], and an example of its effect is a reduction of field of view in radar imaging systems [12,13].
As well as the mentioned wireless systems, wireless positioning estimation accuracy is also affected by the phase impairing issue, in particular when enabled by maximumlikelihood (ML)-based algorithms [14] such as the ones in [15,16]. The trivial workaround to overcome the phase impairing problem is to implement DoA estimation relying only on scalar measurements and ignoring phase information [5,17,18]. Nevertheless, the reduction of measurement information certainly leads to a worsening of overall accuracy [19,20], so a methodology to recover corrupted phase information becomes of fundamental importance for a phase-conscious SBA-based positioning system.
In this paper, we propose a novel method based on TDoA measurements able to reconstruct a coherent relative phase information using a natively phase-impaired SBA architecture. The phase reconstruction is done exploiting wideband signals, which provides high time-resolution [17], making them well suited for accurate TDoA estimations. Exploiting the TDoA information together with the signal pulse shaping, a relative phase information can be reconstructed given that a time synchronization between different antenna readings is given. Note that dealing with time differences, an absolute time reference between transmitter and receiver is not needed, thus the proposed method does not require any signal coding scheme needed to embed propagation delay information within the transmitted data packet.
To validate the proposed technique, in this work we provide a massive campaign of localization experiments carried out in an anechoic chamber, and laboratory instrumentation instead of dedicated electronic systems. This permits the provision of a pure experimental proof of the method based on wideband signals, without any kind of communication protocol and embedded information; furthermore, the use of instrumentation equipments permits reducing significantly measurement uncertainty, thus leading to more reliable conclusions. To emulate the lack of an absolute time synchronization between transmitters and receivers, we introduced a random time delay between different localization experiments. This paper is organized as it follows. Section 2 describes the proposed TDoA method, introducing the analytical model. Section 3 provides the experimental validation of the TDoA method. In particular it describes the design criteria of the positioning system architecture, and discusses the independent experimental validation results.

Application Scenario and Analytical Model
The key point of the method proposed in this work is the capability of extracting coherent phase-displacement information given an array system that is intrinsically phaseimpaired in the working frequency band. For a better understanding of phase-impairing problem, consider an ideal Uniform Circular Array (UCA) in a tridimensional scenario [21] as in Figure 1.
In Figure 1, an RF-continuous wave is impinging towards the UCA, propagating on a trajectory defined by angle (ϑ, ϕ) with respect to array vertical direction. Considering a simple wave propagation model, the propagating plain wave at frequency f 0 impinging on each i-th antenna can be written as where x i , y i , z i are equal to i-th antenna center of phase coordinates respect reference system O as indicated in Figure 1, and k is the wave vector which portrays the wave propagation direction and speed. Applying (1), the phase difference between signals received at two different i-j antennas results and it is time unvarying but dependent on physical antenna positions and signal trajectory. The phase difference is directly proportional to the projection of i-j antennas distance on wave propagation direction ( Figure 1). As described in [16,22], the tridimensional DoA estimation can be placed as finding the solution of the following problem: where C(ϑ, ϕ) is a ML estimator (i.e., MUSIC [21,23,24], ESPRIT [25], minVAR [16,18]), S is the steering vector (or rather the vector comprising amplitude and phase signal readings for each antenna), and M(ϑ, ϕ) is the reference map as shown in [16]. The reference map contains the ideal steering vectors given by model [26,27] for every (ϑ, ϕ) DoA, and the ideal goal for the ML estimator is to get its minimum value only when θ ,φ exactly corresponds to the real DoA. Different ML estimators can overcome different kinds of inequalities between measured and ideal steering vectors, and as shown in [16,22,28], the minVAR estimator resulted to be the most feasible one in case of phase-less scalar power antenna or received signal strength indicator (RSSI) readings.
Considering an array with N antenna elements, the generic steering vector S 0 is defined as where each s 0i is the amplitude-phase measurement of sinusoidal signal received by i-th antenna element, having directional power gain G i (ϑ, ϕ), obtained using an ideal coherent receiver chain. All phase measurements shall be consistent between themselves, or rather they have to keep their sorting order considering the signal propagation path ( Figure 1). By this the receiver chains for each antenna must share the same phase reference α 0 . As it will be shown in Section 3.1, a SBA architecture is a feasible solution to enable such measurement conditions . As proven in [16,21] and indicated in [23,24], the information related to signal DoA is kept within the steering vector through the ratio between its own terms, or rather such information is strictly related to steering vector direction on its N-dimensional vectorial space. Considering the implementation of an ML estimator able to compare steering vectors directions more than single components absolute values, as, for instance, applying a MUSIC [21,23,24] or a minVAR estimator [16] in place of a typical LSE estimator [18], every term of any measured steering vector can be divided for a constant factor. By this, it is possible to choose a generic R-th reference antenna between the available ones and put where G R = G R (ϑ, ϕ) but (ϑ, ϕ) is fixed within a specific measurement, so G R is a constant value for all i-th elements.
Dealing with an incoherent sensing chain, as, for instance, the one for RSSI-based DoA positioning system [29][30][31][32], only the magnitude of each s i term is taken in account, so the information related to phase displacement between different antennas is disregarded reducing the amount of available information and so the potential DoA estimation accuracy.

Phase-Wrapping Issues
Now, let us consider to implement a coherent sensing chain able to get complete amplitude-phase measurements for each term s i of (5). Dealing with a simple scenario as the one in Figure 1, it is trivial that it must result where the sorting order is strictly dependent on signal direction of arrival (ϑ, ϕ). Note that such a sequence is preserved also applying a linear transform on each measured phase, so normalizing the steering vector respect to a generic R-th reference antenna as in (5) it and the ∆α iR terms are the ones depicted in (3).
Considering an ideal sensing chain the phase difference terms as defined in (3) is always monotone with respect to antenna distance. Unfortunately, a real phase detection is always impaired by phase wrapping issue, so measured ∆α iR terms must be placed aŝ where the introduction of modulus operator (mod) has dramatic effects on phase information coherency, disrupting the phase terms sorting order as given in (7). To avoid such phenomena, a condition can be placed imposing that the maximum ∆α iR term should be lower than 2π. Because phase displacement terms are dependent on signal DoA, this would require to dynamically change the reference antenna considering an a priori knowledge of expected direction of arrival, which is nonsensical.
A more affordable condition could be the bounding of the maximum distance between center of phases of different antennas, so that given any reference antenna all terms of (7) will be surely lower than 2π.
Considering (3), such condition can be written as and given an array structure it allows to identify a sub-domain for signal DoAs where the phase ambiguity is avoided (or rather, a "Field of View" similar to one defined by Nyquist sampling conditions for MIMO arrays [12,33,34]). By array geometry, surely it can be put that so it results max and by (11) it is clear how much the dimensions of the entire array structure must be constrained in order to achieve an effective signal DoA estimation on the overall (ϑ, ϕ) domain.
As an example an array working in f 0 = 2.45 GHz, the IEEE 802.11bg WiFi frequency band should have the maximum distance between antennas center of phases equal to λ 0 = 12.24 cm, while shifting to the upper f 0 = 5.2 GHz band for 802.11n such distance lowers to λ 0 = 5.7 cm. Dealing with a SBA system, as, for instance, the one described in Section 3, each antenna of the array will have physical dimensions comparable to λ 0 , thus fulfilling such dimensional constraints will lead to serious device engineering issues.

Phase Reconstruction Method
The constraint on array structure dimensions is given by condition (11), which is a direct consequence of (3). A straight way to relax such condition is to increase the signal wavelength λ 0 , but this would require to shift working frequency band.
This work aims to introduce a method able to unwrap actual phase detections by synthetically lowering the equivalent signal wavelength as placed in (3). For propagating single tone signals, the wavelength is given by [35] where v phase is the phase speed related to single tone propagation on propagation direction, c 0 is the speed of light, and ε R , µ R are the relative dielectric and diamagnetic constants for propagation medium. Checking (12), it is clear that there is not any degree of freedom left for varying the equivalent λ 0 parameter. Dealing with non-monochromatic signals, or rather modulated RF signals, the propagating speed of the entire signal envelope on propagating direction is no more equal to phase velocity but it equals the group velocity (13) [35].
The group speed definition as given in (13) is too complex to be useful for inferring relation between modulating signal parameters and equivalent λ G obtained. Nevertheless, by (13) is clear how further degrees of freedom can be added by introducing a modulated signal in place of a single tone.
In place of a single-tone signal, consider using a multi-tone composed signal such as where B is the pass-band signal bandwidth, M is the number of tones (even), and A i is the signal amplitude for the i-th tone. To simplify the scenario, (14) can be simplified considering a constant signal amplitude A i = A 0 leading to signal in Figure 2a.   (14). Because the signal is composed by a discrete number of harmonics, its envelope is periodic with period and pulse width equal to It is possible to exploit signal periodicity trying to fit a simpler signal model on transmitted signal. For instance in Figure 2b a monochromatic signal coarsely fits on actual signal, preserving the matching simply the signal maximums to envelope maximums. The fit signal model can be written as where f FIT and β FIT are chosen to achieve the hard matching between real signal and fit signal maximums. Note that f FIT is known by reference signal periodicity and it results while β FIT provides the exact alignment of maximums by shifting the fit signal. Until now, the proposed signal model has been fit considering a specific starting time for the signal envelope, thus providing a fixed value for β FIT . Nevertheless, we are dealing with a propagating signal, so changing the observation point we expect to observe a time shift on signal observations in Figure 2b.
Similarly to (1), in the fit signal model the starting phase shift due to changing of observation point can be put as where r 0 is the initial observation point, r the new one, and k FIT the wave vector describing the equivalent fit signal propagation. Note that the latter can be defined as where the v group parameter is still unknown. Despite this, Equation (1) can be rewritten as a typical single-tone propagating wave impinging on the i-th antennas of an array where, exactly as in (1), r i is the i-th antenna center of phase position vector with respect to reference system O, as indicated in Figure 1.
Because the model in (20) exactly resembles the fundamental propagation model, to avoid the phase-wrapping issue the condition in (10), it can be put again that but observing Figure 2b it is clear that λ FIT >> λ 0 . Following (19), the equivalent wavelength for signal fit can be written as so theoretically the condition in (21) can be relaxed considering a multi-tone transmitted signal as defined in (14), optimizing the number of tones M and the bandwidth B. Relaxing condition (22) straightly enables to recover phase-information on phase-impaired scenario as in Figure 1. A pictorial proof of such concept is shown in Figure 3. In Figure 3, three hypothetical impinging signal traces, with transmitted signal defined as in (1), are depicted. As in Figure 2b, the wavelength for f 0 frequency is too small respect to the signal period; in fact, it is not possible to get a meaningful phase difference value considering the RF signal carrier. Conversely, the maximum distance between signal peaks is large, making possible to define a unambiguous phase delay ∆α 21 , ∆α 31 also for higher time delays, or higher distances between different antennas. As a drawback dealing with periodic signals still leaves a risk of ambiguity, and in next Section such point will be discussed.
By a practical standpoint it means that the receiver chain should sample the received signal envelope depicted in Figure 2, and a postprocessing unit has to extract signal features in order to get peak amplitude A i , proportional to each antenna gain G i (ϑ, ϕ), and relative phase β FIT parameters related in some way to received signal relative time shiftings (Figure 3). The right time shifting information is preserved only if the antenna switching mechanism is synchronized to the transmitted signal period, to ensure the coherency of time domains between different acquisitions.
Considering Figure 3, it is conceptually clear that doing feature estimations on different antenna signals and embedding these information in a steering vector, the latter shall contain an overall phase displacement information not impaired by wrapping. By this, a generic DoA estimation algorithm can be involved without any kind of issue.
Note that in (22) the v group is still undetermined, and despite a generic definition of such term is given in (13), achieving an analytical formula for λ FIT could be a challenging task. In the following, a computational way to get the λ FIT and the k FIT terms is given.
Consider that to have an ideal coherent IQ receiver connected through a switch to each different antenna of the array, as in Figure 4. For either I and Q demodulator output channels, the signal is sampled fulfilling the Nyquist sampling condition for signal given in (14) providing f SAMPLE ≥ B, with bandwidth B defined as in Figure 2. The antenna switching mechanism is synchronized to signal period T period to ensure time domain coherency between different antenna readings. Given a transmitted signal s TX (t) = s(t), an ideal radio-link of distance R would provide a received signal equal to where τ corresponds to the propagation delay, G R the receiver antenna gain, and H 0 the coefficient modeling transmission gain and radio link path-loss, as the effects of the radio channel. Applying concepts of dynamic system modeling, it can also be written that where H( f ) is the radio channel transfer function with related impulsive response h(t).
Considering a wide-band signal, the path-loss cannot be defined by a single value, so the channel transfer function must be written as in the case of broadband transmitting and receiving antennas with constant gains G T and G R along the entire signal bandwidth, and a generic phase misalignment γ 0 between f 0 receiver local oscillator and transmitter one [36]. Factorizing the variation term on signal components frequencies with respect to the transform function reference value in the carrier frequency f 0 , the transfer function can be rewritten as The postprocessing unit in Figure 4 is able to get a complete time observation of complex envelope signal as in Figure 2, because IQ demodulation provides the frequency shifting of bandpass signal centered in f 0 to baseband keeping spectrum asymmetries and relative tone phases. The base-band converted IQ digital signal equals to where the f 0 frequency related term is removed due to base-band conversion. The resulting Fourier transform for base-band signal is so its spectrum is exactly the one depicted in Figure 2, shifted to a center frequency f 0 = 0 Hz. To deal with variation in path-loss, a slight modification in transmitted signal must be taken in account to equalize signal alteration due to propagation. Practically, the transmitted signal spectrum is changed to and this allows to get the received signal Let consider again a multi-antenna structure as the one in Figure 1, with the objective of identifying a steering vector as in (5). The Antenna Switch state is changed by the uC after getting enough samples for each antenna. This leads to get different signal streams s BBRX (t, i), one for each antenna, as where the delay τ i and the G i (ϑ, ϕ) antenna gain coefficients take in account the radio-link variations due to the swapping of receiving antenna. The antenna array is not distributed (Figure 1) thus for every antenna the radio-link distance R is considered to be the same thus H 0 is not dependent on antenna index. Note that the Local Oscillator phase misalignment term γ 0 is supposed to be the same between all the antennas, and this condition can be satisfied implementing a PLL to get the phase-lock with equal phase offset on the f 0 carrier signal.
In order to define a steering vector as in (5) a single amplitude-phase coefficient for the equivalent wavelength signal must be produced by each different antenna signal observation. Consider applying the condition on monochromatic signal fit as proposed in Figure 2, imposing equivalence between equivalent and received signals peaks, given that all the different antenna readings share the same time reference as specified in Figure 4.
Assuming a monochromatic signal propagation the i-th steering vector term is and taking in account that the equivalent signal must have the same peak as the real one it must be The normalized term respect the R-th reference antenna results in Because the steering vector term is the complex coefficient which models phase shift and amplitude gain given by the i-th antenna reception on signal, considering the baseband signal observations it results that thus the imposition of peaks equivalence is placed as where τ PKi is known by signal observations. Merging (35) with (36) it results that τ PKi T period (37) so the s i terms can be extracted by features measurements of received s BBRX (t, i) signals, each one evaluated in the peak instant time τ PKi .
As indicated in Section 2.1, all steering vector terms can be normalized respect the same coefficient. Considering a generic R-th reference antenna, each term s i can be rewritten as T period (38) and the latter shows the measurement formula for each antenna of the array. Note that the phase term could lead to phase-wrapping issues and its presence is related to ambiguity given by signal periodicity (Figure 3). To avoid such drawback a further condition can be put on exponential argument in (38). In fact, if it is placed max i,j∈ [1,N] the exponential term can be neglected avoiding any phase-wrapping issue. In addition to this it must be assured that the sample time T sample is lower than minimum delay between different antennas to be sure that every antenna will measure a different peak arrival time.
Considering the maximum signal propagation speed of c 0 , the equation above can be placed with respect to distance between antenna center of phases and when the latter equations are satisfied the steering vector obtained as in (38) will present consistent phase information. Note that while the conditions above enable to get a phase consistent steering vector, the analytical form for k FIT (ϑ, ϕ) as in (19) is still undetermined due to lack of knowledge about the v group term.

Definition of a Valid Observation Model
As summarized in Section 2.1, a valid observation model is needed to calculate the ideal expected steering vectors belonging to the reference map [16,22,24]. Such a model must be entirely defined, as in (4) the reference map M(ϑ, ϕ) shall be defined as a function of DoA domain only. The easiest way to define the latter is to implement the definition for the s i terms as given in (32).
Let consider the expected physical steering-vector related to RF signal propagation. Considering the relationship between the baseband signal and the RF signal it results that whereτ PKi is the real peak delay respect the beginning of signal transmission. Until now, the τ PKi terms have been defined as the time instant associated to the i-th received signal peak within the baseband observation window.
The real received peak timeτ PKi is related to measured τ PKi as where t 0i is the time offset between the transmission and the beginning of signal period observation. Despite such offset is unknown if the period synchronization between different antenna time windows is ensured ( Figure 4) we can suppose that t 0i = t 0 is constant for the entire set of readings. By this the difference between measured time peaks result Considering the transmitted signal spectrum as in (29) and the transfer function as in (26) it results Considering that for each antenna the shape of transmitted pulse is the same it results thatτ so finally we get and the steering vector term defined as this has exactly the same form as given in (5) for the propagation of a monochromatic signal. The signal fit is defined as monochromatic, so let force the signal propagation speed to c 0 , putting Forcing the equivalence between (47) and the signal fit steering vector in (34) results in thus the phase equivalence condition is The second term of (49) surely shows phase wrapping if the array is not rightly scaled considering the f 0 wavelength, as described in Section 2.2.
The first term argument is the phase term related to baseband fit signal as put in Section 2.3. If conditions on transmitting signal in (40) are fulfilled, it is given that such term will not show any phase-wrapping issue thus (50) becomes and merging with (19) and consequently Note that following this, the v group (ϑ, ϕ) and λ FIT (ϑ, ϕ) functions have no more actual physical meaning, even resulting different between different antennas. They become terms needed to build an equivalent model able to describe the expected phase displacement among the array elements. Anyway, now the equivalent wave vector is defined through a deterministic function of the DoA domain which leads to expected evaluations consistent with measurements; thus, the latter can be used to build an effective reference map M(ϑ, ϕ) as required in (4).

TDoA Method Experimental Validation
This section describes a possible hardware suitable for the implementation for the methodology introduced above. To prove the method effectiveness, a massive campaign of DoA estimations was carried out by means of a phase-impaired anchor device based on a SBA array. Thus, Cumulative Density Functions (CDFs) of DoAs estimation errors, with and without the use of phase reconstruction, are compared to assess the improvements given by the proposed solution.

Anchor Switched Beam Antenna Design
The SBA-based anchor device is capable to enable the Space Division Multiplexing Access (SDMA) approach, whose structure is shown in Figure 5. The structure is engineered to be well suited for DoA estimations for phase-less observations, ensuring a complete coverage of hemispherical domain surrounding to the device, according to [37,38]. The six antennas that compose the anchor system are designed to operate in the nominal bandwidth 4-7 GHz. The individual antenna element that composes the SBA, is a 2 × 2 array of patch antennas, whose design was proposed in [39,40]. They are and arranged to cover a semi-spherical range, following the design technique proposed in [41].
Each antenna element works in circular polarization (CP), increasing the robustness of the radio link and thus contrasting the multi-path impairment, regardless of the relative orientation of the target. In Figure 6, the radiation pattern of the six patch antennas implemented in the SBA at the center frequency of 5.5 GHz is reported; for the purpose of the localization purpose, these can be considered to be identical. The Half-Power Beamwidth is reasonably maintained constant across the entire working frequency band. A detailed antenna performance is provided in Table 1, where a distinction between nominal and operative frequency bandwidth is made, with the latter being the frequency range where the realized antenna gain is 8 ± 2 dB.  The implemented positioning procedures take advantage of the capability to steer the radiation pattern to different direction and thus areas in the space domain. The beam control is achieved through the sequential activation of the six antenna elements. Indeed, the SBA rotation is performed by means of a microwave six-pole single-through digitally controlled by a microcontroller, similarly to the system definition of Figure 4.
To provide the needed directivity throughout the frequency band, the antenna element itself is based on a 2 × 2 array, so increasing the overall element aperture, it reaches a dimension of 77 mm × 77 mm. Antenna elements are arranged side by side and are tilted 34 deg with respect to the vertical direction, and each one covers an angular aperture of 72 deg on ϕ angle. Table 2 lists the center points of each antenna element, considering the reference system as in Figure 7.  Observing Table 2 and Figure 7, the directions leading to the worse phase impairments should be the ones with ϑ = 90 deg. In this regard, Figure 8 (top) shows the ideal expected phase differences respect to the reference antenna 1 for center frequency f 0 = 5.5 GHz, while on the bottom the same are represented considering the phase wrapping phenomena. In the latter, the phase ambiguity is manifested for several ranges of ϑ, for example the most evident in the plots are those around ϑ = 90 deg.  Figure 9 shows the expected phase differences considering an equivalent signal at frequency f FIT . The overall smaller values of phase differences related to signal fit could lead to consider about phase measurement sensitivity problems. Nevertheless, actual measurements are not carried out in phase, but, in time, where the time sampling is given according to phase measurement resolution requirements as described in Section 2.3.

Experimental Setup
To prove the effectiveness of the method, a massive campaign of measurement was carried out in a semi-anechoic chamber scenario with the SBA anchor installed on a rotating platform controlled by a computer, thus spanning the entire azimuthal domain. Measurements were acquired through the utilization of a calibrated vector network analyzer (VNA), connected to a computer to get measurements synchronized with the SBA antenna switching, as depicted in Figure 10. The VNA equipment is able to provide a measurement of the radio-link transfer function as shown in Section 2.3. Synchronized measurement routines are used to fulfill the conditions on sampling periodicity as given in Section 2.3, Figure 4. Due to high communication times between computer and VNA, the synchronization condition is satisfied acquiring more than one signal period for each antenna measurement. The adoption of the laboratory instrumentation is due to the need to remove all the possible cause of measurement uncertainty, and thus permits us to focus on the methodological aspects of this work. In particular, the validation of calibrated VNA as equivalent to time domain signal sampling was demonstrated by the authors of [15]. Figure 11 shows a picture of the laboratory set-up schematically represented in the Figure 10, which includes absorbing panels and the VNA hidden by the same panels; the SBA is installed in a rotating platform, while the target in positioned on a teflon pole and moved in different position and getting different heights. The target is emulated by only the transmitting antenna, which is identical to the SBA individual element; the reader can refer to the Table 1 for a summary of the transmitting antenna features.
Nevertheless, to prove method robustness towards a lack of synchronization between the signal transmission and signal reception in the experiment, a random time shift has been introduced in each set of time observations throughout the SBA antennas. Each time observation is given as the Inverse Fourier Transform of an S 21 measurement trace considering a bandwidth B = 3 GHz centered in f 0 = 5.5 GHz and a number of points M = 801.

Experimental Measurements
In this section, we report the results of a massive campaign of measurements carried out putting the target emulated by the transmitting antenna on different heights and different distances; Figure 12 reports the entire set of target positions, with respect to the SBA anchor. An entire set of time observation was taken, for each target position and for a convenient number of SBA rotation angles, spanning the azimuthal domain ϕ with a step of 1 deg.  In Figure 13 are shown the typical measured calibrated data composing a steering vector, acquired by the instrumentation in terms of the magnitude of the transmission parameter S21. By observing the plots we can easily observe from the data associated to the SBA element 2, that the effective antenna bandwidth is between 4.75 GHz to 6.25 GHz, thus justifying its use as reduced bandwidth defined herein above. Note that in this example, antenna 2 is the one aligned with the transmitting antenna, while the others result directed elsewhere according to the SBA geometry shown in Figure 6; still according with the SBA geometry antenna 1 points upward, that is θ = 0 degree and thus provides a low trace in the data plots. The second set of plots shown in Figure 13, provides instead the steering vector data at center frequency, as a function of the azimuthal angle. As expected, in this set of data the shapes of the curve resemble the one of the radiation diagram shown in Figure 6. For each rotation the signal DoA and the related error estimations were estimated following the method introduced in Section 2.1, adopting the MUSIC ML estimator in its scalar implementation [21], dealing with phase-less steering vectors, and in its standard implementation [23,24], with phase-reconstructed data. To assess the validity of the phase reconstruction process, the estimations were also carried out using steering vectors completed with pure phase information related to the center frequency f 0 . The results about all these cases are compared in the Figure 14 where are reported the CDFs for the obtained DoA estimation errors.
In particular, in Figure 14a,b the distribution functions are calculated with respect to each angle subject to estimation, while Figure 14c shows the combined estimation error considering the euclidean distance between the estimated DoA and the real one as withθ,φ being the estimated DoA values and ϑ, ϕ the actual ones. By observing Figure 14 it is clear how much the proposed phase reconstruction method is able to recover a coherent phase information. This is especially noticeable when looking at the accuracy improvement on azimuthal domain ϕ (Figure 14b), where the introduction of center frequency phase information actually worsened overall results. This is coherent with results in Figure 8 where severe phase wrapping phenomena were expected along the azimuthal domain, especially for lower azimuth angles.
In the azimuthal direction the implementation of phase recovery did not lead to accuracy improvements but it prevented the disruption of steering vectors information. Conversely on elevation direction the re-introduction of phase allowed to overcome the lack of antenna elements due to SBA array geometry, giving an huge improvement on estimation accuracy, as shown in Figure 14a.
The experiment was repeated by considering the reduced bandwidth, 4.75-6.25 GHz, where the SBA antenna provides a maximally gain according to Table 1. The DoA estimation error in the reduced bandwidth is reported in the Figure 14c, in comparison with those carried out in the nominal bandwidth. From the data we can recognize that reducing the bandwidth, to that where the antenna performs better, the noise bandwidth reduces and the overall SNR increases thus leading to slightly better DoA accuracy. In the Figure 14c we can estimate that the CDF improves of 1.5% at 20 degree, while it remains constant at 10 degrees.   In addition to pure DoA estimations, the involvement of wide-band signals enables to calculate range estimations. Applying the method already proposed in [15], measuring the peak delay with respect to the transmission time range estimation is possible by simply calculating where d is the index of the antenna showing the maximum absolute peak value for its related time observation, τ PKd is the time of arrival of peak, and t 0 is the signal transmission beginning time. Note that here the knowledge of t 0 parameter is mandatory, thus range estimations require an absolute synchronization between transmitter and receiver achievable through the VNA or applying specific signal coding schemes [42][43][44].
The range resolution is related to timing resolution, and the latter is approximable to inverse of signal bandwidth following signal condition shown in (40), so where with a bandwidth B = 3 GHz it results ∆R ≈ 10 cm. Furthermore dealing with a one-way transmission there is not any issue related to range ambiguity. Position estimation is calculated as where θ ,φ is the estimated DoA for each transmitting position respect the SBA reference system. The positioning error is calculated as the euclidean distance between estimation and real position. Figure 15 shows the CDF for both the range estimation error and the positioning error, being the latter a combination of the error due to range and angle estimations. As shown by plots, a correct time offset compensation let the range estimation error to be comparable with range resolution, i.e., ≈10 cm. The CDFs of final position estimations follow the trends of Figure 14c, upholding the efficacy of the proposed phase reconstruction method.
Furthermore, in this case we report the positioning estimation carried out by considering the reduced operation bandwidth of 4.75-6.25 GHz. We see from the curves, that the range estimation exhibits a slight reduction of the accuracy, nominally at 0.1 m the CDF reduces of 3.7%, as result of the analysis bandwidth reduction. Nevertheless, the combined range-DoA positioning error maintains the same slight advantage from reduced noise bandwidth.

Simulated DoA Accuracy versus SNR
For the sake of a better comprehension of the experimental results, we report briefly the results of a set of simulations carried out by implementing the analytical model of Section 2.1 and applying the technique described in Section 2.3. During this analysis, we considered both the SBA and the transmitting antenna as implemented of radiators characterized by the features described in Table 1 and with the radiation pattern of Figure 6. We adopted the awgn{.} function available in Matlab, which can estimate the incoming signal power and to add a proper level of white noise, such that the signal-to-noise ratio (SNR) reach a given value. The results of this analysis were conducted at different SNR levels, ranging from −15 dB to 10 dB, and the corresponding CDF curves are reported in Figure 16. From the CDF curve associated to SNR = −15 dB, we can recognize the shape already shown in Figure 14c; under a closer inspection, we see that the CDF values at errors 10 degree and 20 degree are in agreement for less than 3% with those experimentally derived. We can thus conclude that the overall SNR affecting the measurement is in the range of −15 dB, and that the analytical model is capable to reproduce with accuracy of less then 3% the DoA experiments. We should not be surprised about the capability of the method to recover the DoA estimate with good accuracy, because in spite the large bandwidth adopted, the method introduces a process gain due to inverse FFT involved in the phase reconstruction.

Conclusions
This work presented a new methodology for recovering phase information corrupted by phase wrapping issues on phase-impaired antenna arrays intended for DoA estimations. The proposed method resulted effective in fixing the loss of accuracy in DoA estimations, leading to an huge accuracy improvement as expected by the proper inclusion of the phase information.
The cost for this benefit is the involvement of wideband signals, but dealing with modern data communication protocols such drawback could be less cumbersome than expected. Furthermore, modern ultra-wideband (i.e., fraction bandwidth > 20%)-based systems often provide a time synchronization protocol which also enables range estimation capability thus allowing to implement complete and cost-effective positioning solutions.
An example of implementation based on calibrated VNA measurements has been shown, leaving to further investigation the development of a system capable of data communication, implementing the functionalities of systems shown in [16,22]. Nevertheless, this work has introduced also some design principles and algorithmic concepts useful for actual system development.