CN102914764A - Time difference positioning method for resisting sampling synchronous error of receiver - Google Patents

Time difference positioning method for resisting sampling synchronous error of receiver Download PDF

Info

Publication number
CN102914764A
CN102914764A CN2012103800688A CN201210380068A CN102914764A CN 102914764 A CN102914764 A CN 102914764A CN 2012103800688 A CN2012103800688 A CN 2012103800688A CN 201210380068 A CN201210380068 A CN 201210380068A CN 102914764 A CN102914764 A CN 102914764A
Authority
CN
China
Prior art keywords
error
synchronous error
receiver
collection
location parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN2012103800688A
Other languages
Chinese (zh)
Inventor
罗来源
万群
文飞
熊瑾煜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
FIFTY SEVENTH RESEARCH INSTITUTE OF CHINESE PEOPLE'S LIBERATION ARMY GENERAL STAFF
University of Electronic Science and Technology of China
Original Assignee
FIFTY SEVENTH RESEARCH INSTITUTE OF CHINESE PEOPLE'S LIBERATION ARMY GENERAL STAFF
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by FIFTY SEVENTH RESEARCH INSTITUTE OF CHINESE PEOPLE'S LIBERATION ARMY GENERAL STAFF, University of Electronic Science and Technology of China filed Critical FIFTY SEVENTH RESEARCH INSTITUTE OF CHINESE PEOPLE'S LIBERATION ARMY GENERAL STAFF
Priority to CN2012103800688A priority Critical patent/CN102914764A/en
Publication of CN102914764A publication Critical patent/CN102914764A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The patent of the invention discloses a high-precision range difference-of-arrival (RDOA) positioning method for synchronous error existing in a sampling clock of a receiver. For the sampling synchronization of sampling clocks of a plurality of receivers, a clock reference is provided by generally using a global position system (GPS) and the high-precision synchronous sampling is realized. However, the leading edge of a pulse signal output by the sampling clock of each receiver has errors, so that sampling of the receiver is excited by a GPS clock pulse to generate a random error of a sampling period, and further the positioning precision of the RDOA is deteriorated. According to the positioning method disclosed by the invention, the influence on the positioning of the RDOA due to the synchronous error of the sampling clock is overcome by using the random change rule of the synchronous error of the sampling clock, so that the positioning precision of the RDOA is remarkably improved.

Description

A kind of time difference positioning method of anti-receiver sample-synchronous error
Technical field
The invention belongs to Passive Location, the high precision when relating in particular to a plurality of receiver sampling clocks and having synchronous error, arrive the location from range-difference measurements method fast.
Background technology
The high precision passive positioning has important using value in radar, sonar, satellite and systems such as communicating by letter.Arriving range difference (Range Difference Of Arrival, RDOA) location is a kind of important high-precision locating method, obtains a wide range of applications in sonar, radar and other radio-positioning navigation field.Existing RDOA localization method is broadly divided into geometry bracketing method and statistical decision method two large classes.Wherein, how much bracketing methods mainly utilize signal source that the RDOA measured value comprises and the geometry site between the receiver, the position of directly calculating target; Statistical decision method can be utilized unnecessary RDOA measured value, according to the statistical property of measuring error, and the optimum position result.Statistical decision method has a lot, comprises that mainly KF (Kalman Filtering), TSE (Taylor Series Expansion), DAC (Divide And Conquer), SI (Spheral Interpolation), SX (Spheral Intersection), PX (Plane intersection) and some improve one's methods.These Classic Geiger location approachs have obtained widespread use in passive positioning.
Patent of the present invention is mainly in the situation that there is the high precision RDOA localization method of synchronous error in the receiver sampling clock.For realizing the synchronous of a plurality of receiver sampling clocks, general using GPS (GPS) provides clock reference in passive location system, realizes the high-precise synchronization sampling.Its ultimate principle is that each receiver sampling clock excites forward position and the gps time of the pulse signal of realizing each receiver sampling clock output synchronous through the time clock of GPS.But because there is error in the forward position of the pulse signal of each receiver sampling clock output, so that there is the stochastic error in a sampling period in the sampling of receiver after the gps clock pulse excitation, cause the RDOA bearing accuracy to worsen.Therefore, be necessary to develop a kind of high-precision locating method of anti-receiver sampling clock synchronous error.
Summary of the invention
The objective of the invention is to exist the situation of synchronous error that a kind of high-precision RDOA localization method is provided for the receiver sampling clock.
According to the method, utilize the characteristics of receiver sampling clock synchronous error in the actual RDOA positioning system, model corresponding synchronous error model, then when the estimating target position, estimate simultaneously the sample-synchronous error, finally obtain hi-Fix by eliminating the sample-synchronous error.
The inventive method comprises the steps:
1) according to receiver sample-synchronous error model and systematic parameter (sampling rate), the discrete synchronous error collection of structure.
2) get concentrated a bit in measuring distance is poor, the offseting of discrete synchronous error, then use least square method to estimate corresponding location parameter, and calculate and estimate residual error.Traversal error collection obtains the location parameter corresponding with discrete synchronous error collection and estimates collection and set of residuals.
3) the search set of residuals is estimated as optimal estimation with the location parameter that wherein least residual is corresponding, also can obtain corresponding sample-synchronous estimation of error simultaneously.
Description of drawings
Fig. 1 has shown that in the range observation variance be σ lThe positioning error accumulated probability curve of the inventive method and SI method in the time of=1.5 meters.
Fig. 2 has shown that in the range observation variance be σ 1The positioning error accumulated probability curve of the inventive method and SI method in the time of=3.0 meters.
Fig. 3 has shown that in the range observation variance be σ 1The positioning error accumulated probability curve of the inventive method and SI method in the time of=6.0 meters.
Embodiment
The RDOA location relates to a plurality of research stations and observes simultaneously same emitter Signals.Take two-dimensional localization as example, emissive source and research station are positioned at same plane, and each research station receives the same signal of emissive source radiation, calculate the RDOA with respect to reference station.Here do not consider the measuring process of RDOA.
Suppose in rectangular coordinate system, the position of N research station is S i(x i, y i), i=1,2 ..., N, radiation source are positioned at S (x, y), are defined as follows parameter:
S iDistance to initial point
r i = x i 2 + y i 2 , i = 1,2 . . . , N - - - ( 1 )
S is to the distance of initial point
r = x 2 + y 2 - - - ( 2 )
S is to S iDistance
d i = ( x - x i ) 2 + ( y - y i ) 2 , i = 1,2 . . . , N - - - ( 3 )
S iAnd S jBetween the emissive source signal propagation distance that records poor
d ij=d i-d j i,j=1,2,...,N,i≠j (4)
N research station is total
Figure BSA00000786574100024
Individual different distance is measured combination, namely produces Individual range difference, wherein only having N-1 is independently, all the other range differences are fully definite by this N-1 range difference.Might as well select S 1Be reference station, calculate one group of independently range difference set, be designated as d=(d 1, d 2..., d N) TIf, there is measuring error, the range difference measured value is d I1=d i-d 1+ n I1, i=2,3 ..., N, wherein n I1Expression S iAnd S 1Between the range difference measuring error, write as vector form N=(n 21, n 31..., n N1) TFor the distributed model of range difference measuring error, suppose n I1Obey the distribution of zero-mean, and n I1And n J1Separate, i here, j=2 ..., N; I ≠ j.
The sampling of a plurality of receivers have error synchronously the time, range difference measured value model can be expressed as d I1=d i-d 1+ n I1+ v I1, i=2,3 ..., N, wherein v I1Expression S iAnd S 1Between range difference measure by the error that causes synchronously and v I1In discrete set { θ, 0, θ }, evenly distribute v I1~U{-θ, 0, θ }, θ is the signal propagation distance in the sampling period.
In location Calculation, with the source to reference station apart from d 1Separate parameter as one and estimate separately, but with its as and the separate redundant parameter of radiation source coordinate directly do least-squares estimation, obtaining d 1After the initial estimation of radiation source coordinate, recycle the restriction relation that exists between them, construct another independently linear least-squares problem, to obtain comparatively accurate location estimation.In little error territory, if the statistical property of known measuring error can also be weighted processing.
The implementation step of localization method of the present invention is as follows:
1) the discrete synchronous error collection of structure is by error v I1Model, the order
v = v 21 v 31 . . . v N 1
Then v has M=3 N-1Plant combination, the set that makes v is I.
2) to any v ∈ I, calculate d I1=d I1-v I1, i=2,3 ..., N.Introduce redundant variables d 1, from (1) to (4) formula conversion obtains the linear orientation equation
1 2 · ( d i 1 2 - r i 2 + r 1 2 ) = - x i 1 x - y i 1 y - d i 1 d 1 , i = 2 , . . . , N - - - ( 5 )
Only consider little error condition, write as matrix form
δ=Ap+ε (6)
Wherein
A = - x 21 - y 21 - d 21 . . . . . . . . . - x N 1 - y N 1 - d N 1 , p = x y d 1 ,
x i1=x i-x 1,y i1=y i-y 1,i=2,...,N;
Suppose Q ε=cov{ ε }=T 1Q 0T 1, the weighted least-squares of p is estimated
p ^ WLS = ( A T Q ϵ - 1 A ) - 1 A T Q ϵ - 1 δ - - - ( 8 )
The variance of estimating is
Q p = cov { p ^ WLS } = ( A T T 1 - 1 Q 0 - 1 T 1 - 1 A ) - 1 - - - ( 9 )
Because d 1 2 = ( x - x 1 ) 2 + ( y - y 1 ) 2 , The structure linear equation
V=Hu+ε 1 (10)
V = ( p ( 1 ) - x 1 ) 2 ( p ( 2 ) - x 2 ) 2 p 2 ( 3 ) , H = 1 0 0 1 1 1 , u = ( x - x 1 ) 2 ( y - y 1 ) 2 , ϵ 1 ≈ T 2 · δp = 2 ( x - x 1 ) 0 0 0 2 ( y - y 1 ) 0 0 0 2 d 1 · Δ p 1 Δ p 2 Δ p 3
Suppose Then the weighted least-squares of u is estimated as
u ^ WLS = ( H T Q ϵ 1 - 1 H ) - 1 H T Q ϵ 1 - 1 V - - - ( 11 )
The variance of estimating is
Q u = cov { Δ u ^ WLS } = ( H T Q ϵ 1 - 1 H ) - 1 = ( H T T 2 - 1 Q p - 1 T 2 - 1 H ) - 1 - - - ( 12 )
Then utilize u and x=[x sy s] TBetween relation, calculate location estimation and the Positioning estimation error covariance thereof of emissive source.
Calculate e=d-v-b, wherein b=[b 21... b M1] T, b m 1 = ( x m - x ^ s ) 2 + ( y m - y ^ s ) 2 - ( x 1 - x ^ s ) 2 + ( y 1 - y ^ s ) 2 ,
Figure BSA000007865741000317
With
Figure BSA000007865741000318
By
Figure BSA000007865741000319
Obtain.The set E that the set U that the estimation that allows v traversal set I obtain M u forms and the estimation of M e form.
3) to e ∈ E, search makes the residual error cost function
Figure BSA000007865741000320
Minimum e is corresponding
Figure BSA000007865741000321
Be final estimation, and obtain corresponding location estimation
Figure BSA00000786574100041
With
Figure BSA00000786574100042
The positioning performance of two kinds of methods is assessed in experiment by the positioning error accumulated probability curve of contrast the inventive method and Chan ' s SI method.Use 4 base stations in the experiment, the base station location parameter is respectively S 1: (35.9,1.4) km, S 2: (20.0,0.0) km, S 3: (5.9,2.3) km, S 4: (32.0,7.0) km.Locating area be 1 dead ahead, base station, 200~400km, about-rectangular area of 100~100km.Locating area is divided into 10km * 10km grid, and target is placed on the summit of grid (totally 201 * 201 impact points) and tests.Sample frequency is 40MHz/s, and the sampling period is 25ns.Measuring error n lStandard deviation sigma iGetting respectively 1.5 meters, 3.0 meters and 6.0 meters tests.
Fig. 1~Fig. 3 illustrates respectively the method for patent of the present invention and Chan ' s SI method at σ iGet respectively the positioning error accumulated probability curve in 1.5 meters, 3.0 meters and the 6.0 meters situations.Transverse axis is relative positioning error
( p ^ ( 1 ) - x ) 2 + ( p ^ ( 2 ) - y ) 2 d 1 × 100 %
Can find out from Fig. 1~Fig. 3, when having the sample-synchronous error, the method for patent of the present invention can significantly improve bearing accuracy.For example, when the standard deviation of range observation is 1.5 meters, the relative positioning error of the method for patent of the present invention is respectively 70%, 96% and 100% less than 0.5%, 1% and 1.8% probability, and the relative positioning error of Chan ' s SI method is respectively 43%, 73% and 92% less than 0.5%, 1% and 1.8% probability.When the standard deviation of estimating in the time difference was larger, the method for patent of the present invention is compared with Chan ' s SI method had clear superiority (referring to Fig. 2 and Fig. 3) equally.
Although the method for with reference to the accompanying drawings this patent having been invented is described with way of example, the invention is not restricted to above-mentioned these details, and the present invention contains various modification or the change of covering within the claim scope.
Industrial applicability
The localization method that the present invention can be proposed be applied to arrive the location from range-difference measurements system, satisfies the hi-Fix requirement in the situation that there is synchronous error in the receiver sampling clock.

Claims (5)

1. the time difference positioning method of an anti-receiver sample-synchronous error, it is characterized in that: utilize the actual characteristics that arrive receiver sampling clock synchronous error in range difference (RDOA) positioning system, set up corresponding synchronous error model, when the estimating target position, estimate simultaneously the sampling clock synchronous error, finally obtain hi-Fix by eliminating the sampling clock synchronous error.Its concrete steps are:
According to receiver sampling clock synchronous error model and systematic parameter (sampling rate), the discrete synchronous error collection of structure.
Get concentrated a bit in arriving the range difference measurement, the offseting of discrete synchronous error, then use least square method to estimate corresponding location parameter, and residual error is estimated in calculating.The discrete synchronous error collection of traversal obtains the location parameter corresponding with discrete synchronous error collection and estimates collection and set of residuals.
The search set of residuals is estimated as optimal estimation with the location parameter that wherein least residual is corresponding, also can obtain corresponding sampling clock synchronous error simultaneously and estimate.
2. the time difference positioning method of anti-receiver sample-synchronous error as claimed in claim 1 is characterized in that: by sampling clock synchronous error v I1Determine N-1 dimensional vector v=[v 21... v N1], wherein N represents the receiver number, again by sampling clock synchronous error v I1Discrete synchronous error collection { θ, 0, θ }, the set of determining all possible N-1 dimensional vector v is error collection I: I = { v 1 , v 2 , . . . , v 3 v - 1 }
v 1=[-θ,0,...,0],v 2=[0,0,...,0],v 3=[θ,0,...,0],v 4=[0,-θ,...,0],...,
Figure FSA00000786574000012
Wherein θ is the signal propagation distance in the sampling period.
3. the time difference positioning method of anti-receiver sample-synchronous error as claimed in claim 1, it is characterized in that: get a vector among the error collection I, in arriving the range difference measurement, offset, then use least square method to estimate corresponding location parameter, and calculate corresponding residual error; Traversal error collection I obtains and error collection I one to one location parameter estimation collection and set of residuals.
4. the time difference positioning method of anti-receiver sample-synchronous error as claimed in claim 1 is characterized in that: determine that least square solution is
p=(A TA) -1A Tδ
Wherein () TThe transposition of representing matrix or vector,
A = - x 21 - y 21 - d ^ 21 . . . . . . . . . - x N 1 - y N 1 - d ^ N 1 , δ = 1 2 d ^ 21 2 - r 2 2 + r 1 2 d ^ 13 2 - r 3 2 + r 1 2 . . . d ^ N 1 2 - r N 2 + r 1 2
x I1=x i-x 1, y I1=y i-y 1,
Figure FSA00000786574000015
v I1The sampling clock synchronous error, i.e. i-1 element of a vector v among the error collection I, d I1Measure for arriving range difference, i=2 ..., N;
Figure FSA00000786574000016
(x i, y i) be the position coordinates of N receiver, i=1 ..., N;
Determine that location parameter is the first two element of p
Figure FSA00000786574000021
Determine location parameter
Figure FSA00000786574000022
Corresponding residual error is ε=(δ-Ap) T(δ-Ap), wherein () TThe transposition of representing matrix or vector;
Each vector v among the traversal error collection I is determined respectively corresponding location parameter
Figure FSA00000786574000023
And residual epsilon, the location parameter that all are corresponding
Figure FSA00000786574000024
Form the location parameter collection, all corresponding residual epsilon form set of residuals.
5. the time difference positioning method of synchronous error when anti-receiver as claimed in claim 1 is sampled is characterized in that: the search set of residuals is estimated as final source location with the location parameter that wherein least residual is corresponding and estimates.
CN2012103800688A 2012-11-24 2012-11-24 Time difference positioning method for resisting sampling synchronous error of receiver Pending CN102914764A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012103800688A CN102914764A (en) 2012-11-24 2012-11-24 Time difference positioning method for resisting sampling synchronous error of receiver

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012103800688A CN102914764A (en) 2012-11-24 2012-11-24 Time difference positioning method for resisting sampling synchronous error of receiver

Publications (1)

Publication Number Publication Date
CN102914764A true CN102914764A (en) 2013-02-06

Family

ID=47613209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012103800688A Pending CN102914764A (en) 2012-11-24 2012-11-24 Time difference positioning method for resisting sampling synchronous error of receiver

Country Status (1)

Country Link
CN (1) CN102914764A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001670A (en) * 2018-06-01 2018-12-14 中国人民解放军战略支援部队信息工程大学 A kind of distributed passive location method and device for combining the time difference and angle

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030017832A1 (en) * 2001-07-18 2003-01-23 Anderson Robert J. Method for estimating TDOA and FDOA in a wireless location system
CN102628931A (en) * 2011-07-14 2012-08-08 刘怡光 Linear relation-based time difference positioning algorithm with high precision

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030017832A1 (en) * 2001-07-18 2003-01-23 Anderson Robert J. Method for estimating TDOA and FDOA in a wireless location system
CN102628931A (en) * 2011-07-14 2012-08-08 刘怡光 Linear relation-based time difference positioning algorithm with high precision

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《传感技术学报》 20120731 姚山峰 等 "时差定位系统的缩比实验原理分析" 第973-979页 1-5 第25卷, 第7期 *
《信号处理》 20051031 周非 等 "一种新的稳定的无线定位算法" 第498-501页 1-5 第21卷, 第5期 *
周非 等: ""一种新的稳定的无线定位算法"", 《信号处理》 *
姚山峰 等: ""时差定位系统的缩比实验原理分析"", 《传感技术学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001670A (en) * 2018-06-01 2018-12-14 中国人民解放军战略支援部队信息工程大学 A kind of distributed passive location method and device for combining the time difference and angle

Similar Documents

Publication Publication Date Title
JP5352422B2 (en) Positioning device and program
KR101224512B1 (en) Positioning method of unknown signal generator based on tdoa method
CN101887128B (en) Method for determining inter-frequency deviation of navigation satellite of global satellite navigation system
CN104714244A (en) Multi-system dynamic PPP resolving method based on robust self-adaption Kalman smoothing
JP5078082B2 (en) POSITIONING DEVICE, POSITIONING SYSTEM, COMPUTER PROGRAM, AND POSITIONING METHOD
CN102508245B (en) Method for determining equivalence of elevation estimation precisions of satellite-borne multi-frequency and multi-baseline InSARs (interferometric synthetic aperture radars)
Ostroumov et al. Compatibility analysis of multi signal processing in APNT with current navigation infrastructure
CN109782289B (en) Underwater vehicle positioning method based on baseline geometric structure constraint
CN104020480A (en) Satellite navigation method for interactive multi-model UKF with self-adapting factors
CN101014874A (en) Independent positioning device and independent positioning method
CN105372691A (en) Long baseline satellite formation GNSS relative positioning method based on ambiguity fixing
CN101682394A (en) Use satellite-signal to determine that the geographical position of tellurian reflector separates
CN107607032B (en) GNSS deformation monitoring system
CN101971047A (en) Device and method for the real-time monitoring of the integrity of a satellite navigation system
CN101680940A (en) Determine that based on weighted least-squares estimation the geographic position of tellurian transmitter separates
CN102749637A (en) Method for realizing accurate positioning of vehicle-mounted GPS (Globe Positioning System)
CN105510945A (en) PPP positioning method applied to satellite navigation landing outfield detection
CN113671546A (en) High-precision relative motion vector algorithm based on double differences between carrier observed value epochs
US20240159529A1 (en) Systems and methods for extending the spatial coverage of a reference pressure network
CN105182367A (en) Novel method for extracting BDS satellite ionosphere puncture point electron concentration
CN106597479A (en) Method for estimating movement speed of ionospheric disturbance by utilizing Beidou base station array data
O'Keefe et al. An investigation of tightly-coupled UWB/low-cost GPS for vehicle-to-infrastructure relative positioning
WO2018056859A1 (en) Collision avoidance method based on relative radio navigation
US20220244407A1 (en) Method for Generating a Three-Dimensional Environment Model Using GNSS Measurements
CN102023301A (en) Satellite selection method for medium earth orbit satellite search and rescue system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130206