Direct Position Determination with Single Sensor Based on Signal Periodicity

Due to practical limitations on size and cost, aerial vehicles generally cannot equip complicated sensors to form sensor array for target localization. In this paper, we investigate the direct position determination (DPD) of stationary source via single moving sensor. First, we analyze artificial signal structure and construct the DPD model with the frame periodicity of artificial signal. /e model incorporates Doppler information extracted from both transformation frames and adjacent samples into target localization. Secondly, we consider the effect of oscillator instability and present an iterative solution for joint estimation of target location and phase noise caused by oscillator imperfection. /e proposed technique fully exploits periodic structure of artificial wireless signal, which leads to significant enhancement in localization performance. Both theoretical analysis and simulations are presented to confirm its effectiveness.


Introduction
Passive localization of emitters is an important issue in many scientific and engineering applications. Traditional two-step localization methods extract intermediate parameter such as angle or time delay of incident signals to determine source positions [1]. In past decades, direct position determination (DPD) technique has become the focus of intensive research. DPD directly exploits sensor outputs to localize target without estimating intermediate parameter. It outperforms traditional two-step methods especially for low signal-tonoise ratio (SNR) and/or limited signal samples [2,3].
Among current DPD study, there has been considerable interest for geolocation based on single moving station. is is mainly due to the fact that both synchronization and data transmission between stations are not required when all processing is done within a single station [4]. Demissie first developed and analyzed a mathematical DPD model with moving station [5]. Based on array data, Capon iterative optimization technique and characteristic of noncircular source are introduced into DPD, respectively, for improving the positioning accuracy [6,7]. ese single-station methods are designed based on the angle information and array sensors. For moving station such as unmanned aerial vehicle, restrictions on payload dimensions and weight pose various challenges on the deployment of multichannel receiver or sensor array.
To address such issue, a variety of techniques such as reflector-aided geolocation have been proposed [4]. Among existing solutions, single-station DPD with Doppler shift is regarded as the most promising one. is approach utilizes Doppler shift caused by sensor motion to determine target location and enables single sensor geolocation [8,9]. ough effective and simple, single-station DPD based on Doppler would encounter two problems in practice. First, complexmodulated symbols would bring phase fluctuation into received signal, so it demands approaches that could separate unknown symbols and recover Doppler shift [10]. For instance, methods proposed in [11,12] utilize perfect temporal coherence of strictly noncircular signals to remove phase ambiguity caused by complex modulation. Single-channel technique used in [13,14] adopts phase curve fitting for joint estimation of both modulated symbols and Doppler information. However, these methods are limited to specific modulations, which makes them pale in practical application. Another inevitable challenge is the effect of oscillator instability. Slight phase noise of oscillator will lead to error in location result and must be considered [15].
In this paper, we consider exploiting signal periodicity to cope with above problem and investigate the DPD algorithm based on single moving sensor. It is worth noting that most artificial signals have periodic structure such as periodic radar pulse and communication signals that are coded in frame for transformation. is periodicity allows for time difference of arrival (TDOA) localization with single moving sensor [16,17]. Similarly, single-channel system presented in [18] virtually rotates receiving antenna with the period of target signal to provide angle information. Motivated by previous work, we analyze the signal structure and exploit signal frame periodicity to abstract Doppler information for target localization. Moreover, practical environment is considered, where phase noise induced by oscillator imperfections is not known exactly but obeys Gaussian distributions with known correlation coefficients. We establish the DPD model with oscillator phase noise. en, an iterative solution is provided for joint estimation of target location and phase noise. Simulation results show the proposed DPD algorithm can achieve better location accuracy compared with existing methods. e main contributions of this paper are listed as follows: (1) With the analysis of wireless signal structure, we divide each signal frame into periodic preamble and data block that contains random information symbols. en, a novel DPD model is provided for dealing with such partially periodic signal. It incorporates the Doppler information extracted from both transformation frames and adjacent samples into target localization.
(2) We combine the established signal model with the probability distribution of oscillator phase noise. To avoid multidimensional search, an iterative optimizing scheme is then presented for updating target position and phase noise alternately. e algorithm fully exploits the location information contained in Doppler shifts and minimizes the effect of oscillator imperfections.
(3) e Cramér-Rao bound (CRB) of the source position estimation based on the received signal model is derived.
e rest of the paper is organized as follows. Section 2 introduces the signal model and formulates the problem. Section 3 presents the DPD method with Doppler information of received signal. Performance analysis is given via CRB in Section 4. In Section 5, simulation experiments are provided for validating the proposed estimator. Finally, conclusion is drawn in Section 6.
Notation. roughout the paper, bold lowercase letters are used for vectors, and bold upper case letters are for matrices.
Superscripts T, H, and * represent the transpose, conjugate transpose, and complex conjugate, respectively. Small letter j denotes the imaginary number, i.e., �� � − 1 √ . symbolizes absolute value and ⊗ represents the Kronecker product. Re {·} and Im {·} signify the real and imaginary parts, respectively ⊙ represents the Hadamard product. represents the Euclidean norm. I m denotes m × m identity matrix. 1 n denotes n-dimensional unit vector, and 0 m is m-dimensional zero vector.

Data Model for Single Receiving Sensor.
Consider the scenario as pictured in Figure 1, where a flying sensor receives the direct path signal from stationary emitter. e sensor intercepts signals at K different points along the trajectory and collects each batch of data in a short interval. After downsampling, the baseband signal collected in the kth interception interval can be written as where p is the target coordinate vector; s(t) represents signal waveform; n k (t) is the additive noise; β k stands for complex fading factor; ϕ k (t) denotes random phase noise which is mainly introduced by oscillator drift [15]; and f k (p) denotes the Doppler shift caused by platform motion, that is, where f c represents signal frequency; c is the speed of light; and q k and v k denote coordinate vector and velocity vector of the receiving sensor at the kth interval, respectively. Equation (2) describes the relationship between f k (p) and source position, which implies that information of source position is involved in the residual frequency of r k (t). Generally, sensor position q k and velocity v k are considered to be known a priori. However, in contrast with active systems where transmitted signal is under control by the receiver, s(t) is not available for passive receivers. e complex-modulated symbols contained in s(t) would bring phase ambiguity into r k (t) and restrict the Doppler shift recovery [10]. Meanwhile, the existence of ϕ k (t) also poses challenges on the abstraction of Doppler shift.

Structure of Artificial Wireless Signal.
To get rid of complex-modulated sequence s(t), we first analyze the structure of artificial wireless signal in the sequel. As shown in Figure 2, it is well known that most artificial signals are characterized by frame transmission and each frame consists of preamble block and data block [19]. Data blocks of signals generally contain audio or video data while the preamble blocks carry advanced information like data coding parameters, service labels, current date, synchronization sequence, and so on. Note that both frame structure and preamble data remain almost unchanged once communication service begins [20]. Hence, we can draw the conclusion that there exists periodicity in artificial signals and frame length T F is exactly the length of period. Such frame periodicity phenomenon also occurs in signals of periodic radar pulse or the radar rotating with constant revolution speed.
Moreover, since pulse-shaping filters in digital modulation are real-valued, the phase ambiguity is mainly bought by information symbols of s(t). As such, there exists no phase change among samples within each symbol. As demonstrated in Figure 2, adjacent samples of an oversampled signal belong to the same symbol with a high probability, so phases of these samples generally remain the same. For narrow-band signal with low baud rate, this phenomenon is frequently encountered.
Both frame periodicity and the phase invariance between continuous samples has been successfully utilized in non-cooperative signal processing, such as direction of arrival estimation [19] and blind frequency recovery [21]. Based on above analysis, we list the assumptions in this study: (A1) Signal source s(t)can be divided into preamble sequence s a (t) and data sequence s b (t), that is, where T F denotes frame length and M is the number of involved frames. σ 2 a , σ 2 b are the power of s a (t), s b (t), respectively. In practice, T F can be obtained from signal standards or directly estimated from received data via TDOA measurements [15], so it is assumed to be known to the receiver. (A2) Narrow-band signal and a slow-fading environment are considered, that is, s(t) remains unchanged in a short period of signal samples and fades independently among observation intervals. Hence, s(t)satisfies where integer N is the number of involved samples and T s stands for sampling interval.
(A3) e oscillator phase noise ϕ k (t), k � 1, . . . , K is statistically independent from interval to interval. Within each observation interval, e jϕ k (t) is a Gaussian stationary process that satisfies where τ denotes time delay and correlation coefficients g(τ) and g(τ) are assumed to be known to the receiver.
Above assumptions describe the temporal coherence of artificial wireless signals. ough above discussions focus on Preamble block Data block Preamble block Mathematical Problems in Engineering the frame structure of artificial signals, the repetition of spreading code or cyclic prefix in orthogonal frequency division multiplexing system also meets Assumption A1. Interested readers may refer to [22] for details. Relying on above assumptions, the problem that we address now is to determine target position directly with r k (t) expressed in (1) and (2).

DPD Model with Signal Periodicity.
Based on Assumption A2, a new observation vector r(t) ∈ C N×1 can be constructed as follows: where scribes the effect of oscillator phase noise and diagonal matrix Γ 1,k � diag η 1,k . e N × 1 steering vector a k (p) has the following form: According to Assumption A1, r k ′ (t) ∈ C M×1 can be similarly constructed as where . , e jϕ k (t− (M− 1)T F ) ] T and Γ 2,k � diag η 2,k . e M × 1 steering vector b k (p) has the following form: With (6)-(9), we can construct array vector where noise vector Using properties of Hadamard product and Assumption A1, we find that Since the diagonal entries of Γ 2,k b k (p)b H k (p)Γ H 2,k are all 1, R b,k becomes the diagonal matrix with σ 2 b as its diagonal entry, i.e., R b,k � σ 2 b I M . en, with properties of Kronecker product and the independence between s a (t) and s b (t), covariance matrix of x k (t) leads to where is determined by temporal correlation of noise, so its elements are allowed to differ from each other. In the following, W k is supposed to satisfy W k � Σσ 2 k , where σ k denotes unknown noise power. MN × MN matrix Σ represents noise covariance matrix. Σ is generally known to the receiver, since it can be determined, for example, using sample statistics from a number of independent, identical experiments.
With above discussions, the Markov-like estimates of p with η k , k � 1, . . . , K can be obtained by fitting covariance matrix of received data to (12) in a weighted least squares sense. However, a multidimensional search is essential for estimating target position p and phase noise vector simultaneously, which has dominant computational complexity [23]. To solve the problem, an iteration optimizing scheme is adopted here. Specifically, the optimization of target position p is performed with fixed η k , k � 1, . . . , K in each iteration. en, the refined estimate of η k , k � 1, . . . , K could be obtained with respect to the latest estimated p. e updating procedure continues iteratively until converged. In the next subsection, the above iterative solution will be derived.

Updation of Target Position.
According to (12), R k can be expressed in vector form as where vec(·) is a vector obtained by stacking columns of the argument on top of each other.
T as the estimation of η k . en, the Markov-like estimate of p with η (b) k is obtained as where C k denotes weighting matrix, c k � vec(R k ), and R k is the covariance matrix constructed by data samples of x k (t).
Generally, C k is chosen to be the inverse of the asymptotic covariance of the residuals [23]. Substitution of (13) into (14) yields where C 0.5 k is a Hermitian square root factor of C k . According to (15), the estimation of ς k that minimizes the cost function yields Substituting (16) into (15) leads to p � arg min As K k�1 c H k C k c is constant during optimization, equation (17) becomes

Mathematical Problems in Engineering 5
Target position p can be directly updated via grid search with (18). Next section investigates the procedure for phase noise vector updating. (6) can be rewritten as

Updation of Phase Noise Vector. Equation
where A k (p) � diag a k (p) . en, the eigenvalue decomposition of its covariance matrix has the following form: where vector u k1 denotes the eigenvector related to the largest eigenvalue λ k1 . U k1 and Λ k1 contain noise eigenvectors and eigenvalues of R 1,k , respectively. Based on subspace theory, A k (p)η 1,k are orthogonal to the subspace spanned by U k1 [8].
Moreover, based on Assumption A3, the covariance matrix of η 1,k and η 2,k is assumed to satisfy where the m, nth entry of G 1 and G 2 is g((m − n)T s ) and g((m − n)T F ), respectively. Denote p as the vector containing estimated source location. With (20) and (21), η 1,k , k � 1, . . . , K can be updated by minimizing the following cost function: where real constant α is generally selected according to the noise level. Based on A3, the presumed oscillator phase noise is statistically independent from interval to interval. Hence, the cost function in (22) can be decoupled as Taking the complex gradient with respect to η 1,k and solving gives Similarly, equation (8) can be expressed as where B k (p) � diag b k (p) . Since R b,k � σ 2 b I M , the eigenvalue decomposition of covariance matrix of r k ′ (t) leads to where vector u k2 denotes the eigenvector related to the largest eigenvalue λ k2 . U k2 and Λ k2 contain noise eigenvectors and eigenvalues of R 2,k , respectively. Similar to (23), η 2,k , k � 1, . . . , K can be updated by minimizing the following cost function: Taking the complex gradient with respect to η 2,k and solving gives e steps of above target localization scheme are listed in Table 1, and related flowchart is shown in Figure 3. As shown in the table, step 3 is engaged to update target location based on received signal, while step 4 is applied for determining phase noise vector η 1,k , η 2,k , k � 1, . . . , K.
e iteration may be stopped by comparing estimation results between adjacent iteration. e main complexity of proposed method is involved in step 3. e computational complexity for estimating p via grid search is about o(8JKM 2 N 2 ), where J represents the number of grid points in the search. When there exists no oscillator phase noise, it merely requires steps 1∼3 to determine target location.

CRB on Position Estimation
In this section, we use the Cramér-Rao bound (CRB) on position estimation as a performance benchmark. By the properties of maximum-likelihood estimators, the asymptotic distribution of p is Gaussian with covariance matrix that equals the CRB.

Numerical Results
In this section, simulation results are provided to demonstrate the effectiveness of proposed algorithm. Without loss of generality, we assume a planar geometry and depict both sensor trajectory and target position in To enable DPD based on frame periodicity, target signal is assumed to have the transmission frame of 3.5 ms and each frame contains 12% periodic preamble symbols. Symbols are transmitted using QPSK modulation with 6 kBd. Signal is received with the sampling rate of 2.5 MHz and its carrier frequency f c � 1.5 GHz. Parameters related to target signal are summarized in Table 2. When adopting DPD processing, data related to 3 signal frames (M � 3) are collected for localization during each observation interval and signal waveform is assumed to remain unchanged within 3 adjacent samples (N � 3).
Fading factor β k , k � 1, . . . , K is independently and randomly drawn from a complex Gaussian random generator with zero mean and the variance of 1. e correlation coefficient is supposed to decrease linearly with τ, i.e., g(τ) � σ 2 g / � τ √ , where σ 2 g denotes the oscillator noise power. Unless otherwise stated, it is fixed at 0 dB. Root mean square error (RMSE) is adopted to evaluate the localization estimation, which is computed by where p(i) is the estimation result in the ith Monte Carlo trial. Firstly, we contrast proposed method with the two-step localization algorithm (designated as two-step) presented in [10] and the array-based DPD algorithm (designated as DPD-array) introduced in [9]. DPD-array determines target localization with both Doppler information abstracted from adjacent samples and the angle estimated with sensor array. Here DPD-array is applied with a uniform linear array with 3 antennas, each spacing half-wavelength apart. e twostep method utilizes nonlinear transformation to recover Doppler information from QPSK modulated signal samples. Figure 5 demonstrates position estimation results, and it shows that RMSE of the proposed method coincides with CRB when input SNR is larger than 3 dB. e result in Figure 5 indicates that the proposed method has the ability to achieve better estimation accuracy than both DPD-array and two-step methods. e improvement mainly results from the employment of signal frame periodicity. Since nonlinear transformations may amplify additive noise or produce interfering components, the achievable performance of two-step method cannot be satisfied at low SNR region.
Secondly, we evaluate the location RMSEs of tested methods with different source positions when SNR is 5 dB. As Figure 6 shows, the performance of each algorithm would degrade with the increase of X coordinate, since the transmitter and observer are spaced further apart with larger X. e result implies that the proposed method is able to attain better estimation accuracy with worse location geometry.
irdly, we discuss some important factors with respect to the proposed algorithm. Specifically, we change the number of preamble symbols to alter preamble percentage in Figure 7 and alter the involved signal frames M in Figure 8. Other parameters including frame length remain unchanged. Figure 7 shows the influence of percentage of preamble symbols on proposed estimator. Since nonperiodic data symbols cannot produce b k (p) defined in (9) for localization, the proposed method achieves lower RMSE with less preamble symbols involved in Figure 7. When SNR is 5 dB, the achievable RMSE with 5% preamble symbols approaches 0.5 km. RMSEs versus the number of involved signal frames M are shown in Figure 8. e results demonstrate that the performance of proposed method degrades with less signal frames. e phenomenon is due to the fact that we exploit Doppler information for target localization. e accuracy of Doppler shift estimation is mainly determined by signal collection time, which corresponds to MT F in the proposed method.
Moreover, we evaluate the proposed method when signal period is not exactly known to the receiver, i.e., T F � T F + T e , where T e is the error caused by oscillator or other imperfections. Figure 9 shows location estimates with distinct period errors. From the figure, we find that the performance of the proposed method degrades significantly with error larger than 1.5 μs.
is phenomenon can be predicated. When period error and symbol duration are of the level, the frame periodicity assumption given in (A1) would become invalid and the Doppler extracted from transformation frames cannot be utilized for DPD. Note that the proposed method still works well when T e < 1 μs, which indicates that it provides certain robustness to such imperfections. Figure 10 shows RMSEs versus the power of oscillator phase noise when SNR � 5 dB. As phase noise brings ambiguity for localization, it can be seen from Figure 10 that RMSEs would become larger with the increase of noise power. For phase noise smaller than 5 dB, its influence can be ignored.

Conclusion
In this paper, a novel DPD technique has been introduced to achieve single sensor geolocation. First, we analyze the typical structure of artificial wireless signals and establish the DPD model including oscillator instability. en, an iterative optimization scheme is presented for the joint estimation of target position and oscillator phase noise. e proposed method incorporates Doppler information caused by sensor movement into target localization and achieves good localization performance. ough analysis of the paper mainly considers single source, the proposed technique can also be applied to nonperiodic source with known sequence or multiple sources with distinct periods. DPD with frame periodicity has a number of associated research challenges, including low-complexity implementation, robust algorithm for wideband or multipath source, investigation of interference effects, and determination of the value of signal period. ese will be the focus of our future work.

Data Availability
e data used to support the findings of this study are available from the first author upon request (wang-cheng1988zz@gmail.com).