CN104965191A - Two-site time difference positioning method - Google Patents
Two-site time difference positioning method Download PDFInfo
- Publication number
- CN104965191A CN104965191A CN201510349211.0A CN201510349211A CN104965191A CN 104965191 A CN104965191 A CN 104965191A CN 201510349211 A CN201510349211 A CN 201510349211A CN 104965191 A CN104965191 A CN 104965191A
- Authority
- CN
- China
- Prior art keywords
- delta
- theta
- path difference
- value
- sin
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/04—Position of source determined by a plurality of spaced direction-finders
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
In view of the state of the art that the current plane time difference positioning needs to use three sites, the invention provides a method of only using two sites to realize target positioning. Firstly, through carrying out simplified processing on a one-dimensional two-base positioning equation, a single-base direction finding formula with a measurement basis in the midpoint of the two sites is acquired; then, a least square method is used for correcting accuracy of the single-base direction finding formula, on the basis of features that a compensation function obtained through curve fitting is not sensitive to target distance changes, a method of using a target bearing ideal value, a single-base path difference measurement value and the actually-used base line length value for compensating direction finding deviation is given. On the basis, according to a same difference relationship between the single-base path difference and the two-base path difference proved by a two-base three-site path difference equation, and an actually-measured single-base path difference value and an arrival angle through correction, two path difference values corresponding to a one-dimensional two-base array are obtained recursively, and thus two independent positioning equations meeting requirements of two-dimensional plane positioning can be equivalently acquired.
Description
Technical field
The present invention relates to radio-location technology, be specifically related to a kind of dual station time difference positioning method.
Background technology
Positioning using TDOA is that a kind of time of arrival of being measured the same pulse of same radiation source by multiple receiver is poor, is then obtained the method for radiation source positions according to the position of each platform and the time difference by central processing station.Under modern measure technology condition, due to the raising of time difference measurements precision, time difference position technolot has become a kind of more accurate localization method, and is widely used.
The mistiming being arrived two stations by target emanation signal can determine the hyperboloid that a pair is focus with two stations, the hyperbolic position therefore time difference locating technology is otherwise known as.According to existing hyperbolic navigation positioning principle, two dimensional surface at least must use three stations, obtain two independently path differences thus, in order to construct two independently positioning equations.Three dimensions then needs to arrange four stations, and in order to produce three pairs of hyperboloids, intersect to obtain line by face, face, segment-Facet intersection invocation point realizes the location Calculation to target location.
Summary of the invention
The state of the art of three websites must be used for existing plane positioning using TDOA, the object of the invention is to utilize through single base direction finding value of Curve fitting compensation, and the equal difference relation between single double-basis path difference give a kind of dual station time difference measurement that only uses can the method for realize target location.
The present invention is achieved through the following technical solutions.
First be single base direction finding formula at single baseline point midway place by obtaining a witness mark benchmark to the simplify processes of one dimension double-basis eikonal equation.
Subsequently, least square method is adopted to give the compensation method that effectively can be improved single base direction finding accuracy.First according to the target range that location systems art index is given, or by believable information estimating target distance; The theoretical value of the angle of arrival, the actual base length node-by-node algorithm path difference used is utilized subsequently on the detection interval of whole target azimuth; Further for the consideration to substantial measurement errors, and unusual in order to avoid calculating appearance, artificial correction is carried out to path difference theoretical value; Then pre-estimating the deviation of the target angle of arrival by the path difference value through revising, the correction to single base direction finding value can be realized by least square method on this basis.
Although utilize the method for curve can obtain the direction finding Exact Solutions of restraining, curve-fitting method inherently needs the Exact Solutions using direction finding, and single base direction finding only can provide approximate solution comparatively accurately.The reason finally breaking away from this predicament analyzes to find that the penalty function obtained by curve is mainly relevant to the change of base length and path length difference, and insensitive to the change of target range, especially when base length is less than 100 kilometers.Just based on this characteristic, just can adopt the approximate evaluation value of target radial distance, by the ideal value of target azimuth, base length and actual measurement path difference, calculate the offset of direction finding deviation.
In fact, based on this feature, another kind of compensation way is according to system works distance, and precalculate and store the direction finding bias compensation value corresponding to certain specific base length, the actual measurement approximate value then when real work according to the target angle of arrival carries out real-time calling.
On the basis by compensating the accurate direction finding solution of acquisition, can demonstrate,prove have a kind of equal difference relation between single base path difference and double-basis path difference by double-basis three station eikonal equation, based on this equal difference relation, utilize single base path difference value of actual measurement and the angle of arrival through revising, just recursion can obtain two the path difference values corresponding to one dimension double-basis array, the acquisition that just energy is of equal value thus meets two location-independent equations required for two dimensional surface location.
Specifically comprise the following steps:
The time difference measurement of step 1, the baseline mid point angle of arrival.If form a positioning using TDOA system by two stations, utilize the time difference of measuring and obtaining, be calculated as follows the target angle of arrival θ in two station baseline midpoint:
In formula: Δ R
13it is the path difference between two stations; v
cfor the light velocity; Δ t is the time difference, D
0for base length.
Step 2, on whole object position-sensitive detecting interval, calculate the theoretical value of path difference, and revised.According to the target range that location systems art index is given, or by believable information estimating target distance, utilize the theoretical value of the angle of arrival, the actual base length used, the distance by cosine law node-by-node algorithm baseline two ends:
Obtain path difference thus:
ΔR
i=r
1i-r
3i(4)
In formula: r
2it is baseline midpoint target range.Different from formula (1), for succinctly, the subscript of the path difference calculated value when the different angle of arrival no longer comprises site information.
For the consideration to measuring error, and unusual in order to avoid occurring when calculating, artificial correction is carried out to path difference theoretical value;
ΔR
ai=a
0ΔR
i(5)
In formula: a
0< 1, span is about about 0.5.
Step 3, utilize least square method to the correction of target angle of arrival measured value.By the theoretical value of the angle of arrival, actual base length and corrected path difference value calculate the deviation data (θ of the angle of arrival
i, Δ ε
i):
Δε
i=2D
0sinθ
i-ΔR
ai(6)
And the linear combination function be constructed as follows:
Δε=c
1f
1+c
2f
2+c
3f
3+c
4f
4+c
5f
5+c
6f
6+c
7f
7+c
8f
8(7)
Wherein: c
ifor undetermined constant.F
iknown function for following:
f
3=f
1+f
2
f
5=(f
1+f
2)
3
f
6=(f
1+f
2)
7
f
7=(f
1-f
2)
3
f
8=(f
1-f
2)
8
Following linear equation is set up by least square method:
Ac=y (8)
Wherein:.
Solve coefficient thus:
c=A\y (11)
The deviation delta ε that so far just can try to achieve by curve fitting method revises single base direction finding value:
In formula: Δ R
a=a
0Δ R
13for the path difference measured value through revising.Absolute value in formula is to eliminate in the imaginary part phenomenon that likely there will be close to deviation delta ε during baseline axis direction.
Practical application also can precalculate according to operating distance and store the direction finding bias compensation value corresponding to certain specific base length, and the actual measurement approximate value then when real work according to the target angle of arrival carries out real-time calling.
Step 4, demonstrate,prove to obtain by double-basis three station eikonal equation the equal difference relation that has between single base path difference and double-basis path difference, based on this equal difference relation, utilizing single base path difference value of actual measurement and the angle of arrival through revising, just recursion can obtain two the path difference values corresponding to one dimension double-basis array.
From geometric relationship, single base path difference of corresponding array total length is the path difference sum of adjacent two baselines:
ΔR
13=r
1-r
3=(r
1-r
2)+(r
2-r
3)=ΔR
1+ΔR
2(13)
And the length sum of adjacent two baselines equals array total length.As by single base path difference by one dimension double-basis path difference direction finding formula:
In path difference Δ R
1or Δ R
2displace, then have:
Therefrom can obtain the result being similar to arithmetic series:
Wherein:
Step 5, try to achieve the distance and bearing in array end points target by one dimension double-basis positioning equation, at the distance and bearing of baseline right-hand member place target be:
Wherein:
ΔR
13=v
cΔt (22)
Characteristic of the present invention:
1. decrease the quantity of detection site, not only reduce arrangement and the use cost of positioning system, but also reduce the difficulty of addressing.
2., because only comprising a time difference measurement, so for same base length, the overall measurement error of the time difference locating technology system being more than or equal to three than website quantity is little by the overall measurement error of bistatic location.This is because total observational error will be the component sum of measuring error of multiple different independently metering system, if independently observed quantity is more, then overall measurement error will more be difficult to be lowered.
3., under positioning using TDOA system, the locating and tracking precision of target is main relevant to the relative geometrical relation between target and observer when measuring error etc. remains unchanged.Owing to only there being two stations, the impact of arrangement manner on the locating and tracking precision of target will effectively be reduced, and thereby simplify the analysis of positioning system.
4. not only can be used for early detection, more can be used for in-plant navigator fix.
Accompanying drawing explanation
Fig. 1: one dimension Bistatic
Fig. 2: compensate curve during different radial distances
Fig. 3: relative range error during different stop spacing
Fig. 4: angle measurement error
Embodiment
Further illustrate the present invention below in conjunction with accompanying drawing 1-Fig. 4 how to realize.
Embodiment
A kind of dual station time difference positioning method.Accompanying drawing 1 is one dimension Bistatic; Figure 2 shows compensate curve during different radial distances; Relative range error when accompanying drawing 3 is different stop spacings; Accompanying drawing 4 is angle measurement error.
This patent proposes a kind of dual station time difference positioning method, this cardinal principle of an independently path difference value realize target location that only utilizes is: obtaining on the basis of accurate direction finding value compensating based on single baseline time difference measurement and matching, the equal difference relation had between single base path difference and double-basis path difference is utilized to solve path difference value corresponding to double-basis line, two independent equations be met required for plane positioning that thus just can be of equal value.
This method solving double-basis path difference by the single base path difference of actual measurement must rely on direction finding solution accurately, the main cause being finally able to engineer applied is the impact of change on compensate value of target range is insensitive, just can improve the accuracy of single base direction finding thus by least square method.
One, three path difference location, station
1, fundamental equation
For the one dimension double-basis array shown in accompanying drawing 1, following two geometry subsidiary equations can be listed by the cosine law:
In formula: θ
2it is the target angle of arrival of central site; D
ithe length of baseline; r
ifor radial distance.
Because having:
x=r
2sinθ
2(3)
In formula: x is the horizontal ordinate in rectangular coordinate system.
Therefore geometry subsidiary equation can be rewritten as:
2, interior point location
According to corresponding to baseline D
1and D
2path difference solve positioning equation, then the middle end points 2 of receiving array will become the public connecting end point of adjacent two baselines, claim this kind of locator meams be based on put in baseline path difference location.
Corresponding to baseline D
1and D
2eikonal equation respectively:
ΔR
1=r
1-r
2(6)
ΔR
2=r
2-r
3(7)
In formula: Δ R
iit is path difference.
Above-mentioned eikonal equation is substituted into geometry assist type (4) and (5), after transposition arranges, has following binary once linear system of equations:
Therefrom directly can solve the target range measured based on path difference:
And the target angle of arrival:
As equal in adjacent two baselines: D
1=D
2=D
0, then range finding and direction finding expression are:
3, outer point location
If one of selected path difference is the total length of array, then some outer point will become the points of common connection of two baselines, claim this kind of locator meams to be locate based on the path difference of exterior point.
Corresponding to baseline total length D
13eikonal equation be:
ΔR
13=r
1-r
3(14)
By above formula with corresponding to baseline D
2path difference formula (7) substitute into geometry assist type (4) and (5), transposition arrangement after have following binary once linear system of equations:
Therefrom directly can solve the target range at website 3 place:
Can solve again:
The target angle of arrival of trying to achieve thus on the right side of one dimension double-basis array is:
As equal in adjacent two baselines: D
1=D
2=D
0, then range finding and direction finding solution are:
Two, the path difference relation between single double-basis
From geometric relationship, the path difference of corresponding array total length equals the path difference sum of adjacent two baselines:
ΔR
13=r
1-r
3=(r
1-r
2)+(r
2-r
3)=ΔR
1+ΔR
2(22)
And the length sum of adjacent two baselines equals array total length.As by this relational expression respectively by the path difference Δ R in one dimension double-basis path difference direction finding formula (13)
1or Δ R
2displace, then have:
Therefrom can solve two path difference values corresponding to one dimension Bistatic respectively:
Wherein:
Demonstrate,prove thus, in one dimension Bistatic, the path difference value corresponding to each single baseline is 1/2nd of the path difference value corresponding to total base length just, then increases and decreases an identical departure.Analog computation shows, small departure Δ a only could be solved by correct when the angle of arrival itself accurately solves, and as adopted approximate direction finding value, then the result obtained will be disabled.
Three, single base direction finding and curve compensation
1, single base direction finding
According to existing analysis, two adjacent high-order path differences
can be thought equal by approximate, so equidistant direction finding formula (13) can be reduced to:
Cause:
ΔR
13=r
1-r
3=(r
1-r
2)+(r
2-r
3)=ΔR
1+ΔR
2
There is single base direction finding formula:
As can be seen here: the witness mark reference point of single baseline direction finding is the point midway place at single baseline, instead of at the end points place, left and right of single baseline.And analog computation result shows, be not only correct to Short baseline to the single base direction finding formula (29) after reference mark shift-corrected, but also can be applicable to comparatively Long baselines.
2, the compensation of deviation
Analog computation shows, simplifies because of approximate, has a small deviation, single base direction finding formula (29) can be rewritten as this reason between the calculated value and theoretical exact value of single base direction finding formula:
In formula: Δ ε is departure.And in order to simplify subsequent analysis, start the subscript omitting the angle of arrival from here.
By transposition process, formula (30) can be represented as:
Δε=2D
0sinθ-ΔR
13(31)
Deviation can be regarded as the difference between approximate path difference and actual measurement path difference thus.On the basis of formula (31), following linear combination function can be constructed:
Δε=c
1f
1+c
2f
2+c
3f
3+c
4f
4+c
5f
5+c
6f
6+c
7f
7+c
8f
8(32)
In base: c
ifor undetermined constant; f
ifor known function.
By sunykatuib analysis, utilize only known parameter D
0with Δ R
13, each function item can be expressed as:
f
3=f
1+f
2
f
5=(f
1+f
2)
3
f
6=(f
1+f
2)
7
f
7=(f
1-f
2)
3
f
8=(f
1-f
2)
8
Some function items are wherein that the mode calculated by hand simulation is obtained, the most basic function f
1and f
2in fact be exactly the sine and cosine functions of the approximate target angle of arrival.
By direction finding value and the path difference value acquisition deviation data (θ of whole search coverage
i, Δ ε
i), and set up following linear equation:
Ac=y (33)
Wherein:
Least square solution be c=A y.After trying to achieve deviation delta ε by matching, when in the end calculating the exact value of the angle of arrival, attention must add absolute value, namely has:
This is to eliminate in curve fitting process, and when close to baseline axis direction, departure Δ ε likely there will be imaginary number phenomenon.
3, to the consideration of path difference measuring error
In fact, when actual measurement, path difference value itself, inevitably with measuring error, to this, can adopt a coefficient to carry out approximate representation path difference measured value departing from exact value:
ΔR
13w=(1±w)ΔR
13(37)
In formula: ± w simulates measured value to the deviation of exact value.
Simulation subsequently shows, when supposing that measured value is less than ideal value, utilize curve fitting method still can obtain original or even better Compensate approximate result, but when measured value is greater than exact value, by appearance because path difference value is greater than single base length, thus cause the function containing extraction of square root:
Abnormal situation is there is when calculating.Be artificial reduction path difference value to this solution, such as, be multiplied by the coefficient (such as getting 0.5) that is less than 1 on the right of expression:
ΔR
13w=0.5(1±w)ΔR
13
After process like this, still can ensure the convergence of compensation calculation when w gets+0.5, and the calculating of direction finding accuracy is not affected, and still can keep original level.But should be noted that the path difference now in single base direction finding formula when calculating, and when carrying out least square method and calculating, all should adopt the path difference value through artificially simulating amendment.Finally from the needs of actual computation, and in order to be concise in expression, the measuring error coefficient of reflection path difference and artificial correction factor are combined, given path difference amendment type is:
ΔR
13a=a
0ΔR
13(38)
In formula: a
0be one and considered the measuring error of path difference and the coefficient of artificial modifying factor.
4, the variation characteristic of deviation
Based on curve fitting method, need to know direction finding value accurately to the acquisition of deviation data, but measuring the obtainable angle of arrival with the path difference between dual station is only an approximate value, and as do not obtained direction finding solution accurately, then aforesaid whole analysis seems to be only a theoretic proof of pure mathematics.If still can obtain deviation data when the angle of arrival is not also determined accurately, then the penalty method based on curve could be used for engineering reality.
Further sunykatuib analysis proves, as done larger depreciation correction to path difference artificially in advance, the calculated value such as getting path difference is the half of actual measured value, upper trifle has illustrated to do so can't affect the correct of corrected Calculation and accuracy, then when base length is determined, the impact of change on compensate value of target range is very little.If base length is only in the scope below 100 kilometers, then the compensation of deviation will present substantially not with the characteristic of radial distance change.
Figure 2 shows when base length is 50 kilometers, to the deviation compensation curve of different radial distances target.Based on compensating substantially not with the characteristic of radial distance change, we just for specific base length, can obtain by precalculating and preserve deviation data.Single base direction finding formula (29) can be first adopted to know the angle of arrival of target roughly comparatively accurately during practical engineering application, then from the deviation data of storage, select corresponding deviation data in order to revise to direction finding formula according to measured angle value, thus just can obtain direction finding value more accurately.
Analog computation shows, when baseline is larger, such as, is greater than hundreds of kilometer, and the polymerization property of bias compensation value will be deteriorated, and now according to different radial distances, can adopt the mode of segmented compensation.
Four, only based on the positioning error of single base path difference
1, basic skills
Use the geometric relationship shown in accompanying drawing 1, be located at position 1 and 3 place and be provided with two detection sites, and adopt time difference measurement to obtain corresponding to single Long baselines D
13path difference Δ R
13, compensated direction finding deviation by curve subsequently, make to be positioned at Long baselines mid point, the direction finding accuracy at position 2 place effectively improves, and can be used for solving the deviation delta a corresponding to double-basis path difference.Now, recycle relational expression (25) or (26) that channel syndrome obtains, just can accurately obtain each path difference value corresponding to one dimension double-basis array.
On this basis, just the present entity site location of the fructufy obtained by exterior point positioning analysis can be utilized the range finding of target and direction finding:
Wherein:
ΔR
13=v
cΔt
2, relative range error
For ease of error analysis, if transition function:
q=2(ΔR
13-2ΔR
2) (40)
Namely have:
When not considering the measuring error of path difference, the differential of distance to time difference Δ t is:
Wherein:
And have:
For path difference Δ R
2, in current analysis, temporarily wherein the comprised deviation delta ε being used for compensating path difference is regarded as a small constant, partial differential obtained thus is:
According to error theory, ignore the error effect of spacing between station, the relative range error only produced by time difference measurement is:
In formula: σ
Δ tfor the root-mean-square error of time difference measurement.
Fig. 3 gives relative range error curve during different stop spacing, as seen from the figure, be greater than after 3 kilometers at baseline, in larger angle of arrival region, substantially can meet the technical requirement of 0.1%R.
3, angle measurement error
Have direction finding formula (21) direct differentiation:
The angle measurement error only produced by time difference measurement is:
Fig. 4 shows angle measurement error curve during different stop spacing.
4, conclude the speech
Analysis result shows, the distance accuracy only obtained by single base time difference measurement than double-basis time difference measurement result as well, this mathematical form should be rational, because, for double-basis positioning using TDOA, general orientation error will be two error component sums that are different, independently time difference measurement, and for single base positioning using TDOA, then only comprises the error component of a time difference measurement.
Claims (6)
1. a dual station time difference positioning method, is characterized in that utilizing through single base direction finding value of Curve fitting compensation, and the equal difference relation between single double-basis path difference obtains and meets required for two dimensional surface location two independently positioning equations.
2. method according to claim 1, is characterized in that by obtaining a witness mark benchmark to the simplify processes of one dimension double-basis positioning equation be single base direction finding formula at dual station mid point:
3. method according to claim 1, is characterized in that calculating path difference theoretical value on whole object position-sensitive detecting interval, and is revised.According to the target range that location systems art index is given, or by believable information estimating target distance, utilize the theoretical value of the angle of arrival, the actual base length used, the distance by cosine law node-by-node algorithm baseline two ends:
Obtain path difference thus:
ΔR
i=r
1i-r
3i(4)
In formula: v
2it is the target range of baseline midpoint.Different from formula (1), for succinctly, the subscript of the path difference calculated value when the different angle of arrival no longer comprises site information.
Further for the consideration to substantial measurement errors, and unusual in order to avoid calculating appearance, artificial correction is carried out to path difference theoretical value;
ΔR
ai=a
0ΔR
i(5)
In formula: a
0< 1, span is about about 0.5.
4. the method according to claim 1,2 and 3, is characterized in that utilizing the measured value of least square method to the target angle of arrival to revise.By the theoretical value of the angle of arrival, actual base length and corrected path difference value calculate the deviation data (θ of the angle of arrival
i, Δ ε
i):
Δε
i=2D
0sinθ
i-ΔR
ai(6)
And the linear combination function be constructed as follows:
Δε=c
1f
1+c
2f
2+c
3f
3+c
4f
4+c
5f
5+c
6f
6+c
7f
7+c
8f
8(7)
Wherein: c
ifor undetermined constant.F
iknown function for following:
f
3=f
1+f
2
f
5=(f
1+f
2)
3
f
6=(f
1+f
2)
7
f
7=(f
1-f
2)
3
f
8=(f
1-f
2)
8
Following linear equation is set up by least square method:
Ac=y (8)
Wherein:
Solve coefficient thus:
c=A\y (11)
The deviation delta ε that so far just can try to achieve by curve fitting method revises single base direction finding value:
In formula: Δ R
a=a
0Δ R
13for the path difference measured value through revising.Absolute value in formula is to eliminate in the imaginary part phenomenon that likely there will be close to deviation delta ε during baseline axis direction.
Practical application also can according to system works distance, and precalculate and store the direction finding bias compensation value corresponding to certain specific base length, the actual measurement approximate value then when real work according to the target angle of arrival carries out real-time calling.
5. the method according to claim 1 and 4, it is characterized in that demonstrate,proving to obtain by double-basis three station eikonal equation the equal difference relation that has between single base path difference and double-basis path difference, based on this equal difference relation, utilizing single base path difference value of actual measurement and the angle of arrival through revising, just recursion can obtain two the path difference values corresponding to one dimension double-basis array.
From geometric relationship, single base path difference of corresponding array total length is the path difference sum of adjacent two baselines:
ΔR
13=r
1-r
3=(r
1-r
2)+(r
2-r
3)=ΔR
1+ΔR
2(13)
And the length sum of adjacent two baselines equals array total length.As by single base path difference by one dimension double-basis path difference direction finding formula:
In path difference Δ R
1or Δ R
2displace, then have:
Therefrom can obtain the result being similar to arithmetic series:
Wherein:
6. the method according to claim 1 and 5, is characterized in that trying to achieve distance and bearing in array end points target by one dimension double-basis positioning equation:
Wherein:
ΔR
13=v
cΔt (22)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510349211.0A CN104965191A (en) | 2015-06-23 | 2015-06-23 | Two-site time difference positioning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510349211.0A CN104965191A (en) | 2015-06-23 | 2015-06-23 | Two-site time difference positioning method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104965191A true CN104965191A (en) | 2015-10-07 |
Family
ID=54219232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510349211.0A Withdrawn CN104965191A (en) | 2015-06-23 | 2015-06-23 | Two-site time difference positioning method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104965191A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054134A (en) * | 2016-05-20 | 2016-10-26 | 东南大学 | Rapid positioning method based on TDOA |
CN111323746A (en) * | 2020-03-19 | 2020-06-23 | 哈尔滨工程大学 | Double-circular-array azimuth-equivalent delay inequality passive positioning method |
CN113238217A (en) * | 2021-06-03 | 2021-08-10 | 哈尔滨工业大学 | Distributed high-frequency ground wave radar combined positioning method based on interior point method |
US11397256B2 (en) * | 2017-03-15 | 2022-07-26 | Cypress Semiconductor Corporation | Estimating angle measurements for source tracking using a phased array system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3103386A1 (en) * | 1981-02-02 | 1982-08-12 | C. Plath Gmbh Nautisch-Elektronische Technik, 2000 Hamburg | Method for locating a modulated transmitter |
CN103091660A (en) * | 2012-09-24 | 2013-05-08 | 郁涛 | Virtual base line direction-finding method based on phase jump correction |
CN104049237A (en) * | 2013-03-15 | 2014-09-17 | 郁涛 | Passive positioning technology based on difference measurement |
-
2015
- 2015-06-23 CN CN201510349211.0A patent/CN104965191A/en not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3103386A1 (en) * | 1981-02-02 | 1982-08-12 | C. Plath Gmbh Nautisch-Elektronische Technik, 2000 Hamburg | Method for locating a modulated transmitter |
CN103091660A (en) * | 2012-09-24 | 2013-05-08 | 郁涛 | Virtual base line direction-finding method based on phase jump correction |
CN104049237A (en) * | 2013-03-15 | 2014-09-17 | 郁涛 | Passive positioning technology based on difference measurement |
Non-Patent Citations (2)
Title |
---|
郁涛: "大跨距双站交会角的单站观测", 《航空电子技术》 * |
郁涛: "对相似递推程差的一种修正方法", 《现代导航》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054134A (en) * | 2016-05-20 | 2016-10-26 | 东南大学 | Rapid positioning method based on TDOA |
US11397256B2 (en) * | 2017-03-15 | 2022-07-26 | Cypress Semiconductor Corporation | Estimating angle measurements for source tracking using a phased array system |
CN111323746A (en) * | 2020-03-19 | 2020-06-23 | 哈尔滨工程大学 | Double-circular-array azimuth-equivalent delay inequality passive positioning method |
CN111323746B (en) * | 2020-03-19 | 2023-05-05 | 哈尔滨工程大学 | Direction-equivalent time delay difference passive positioning method for double circular arrays |
CN113238217A (en) * | 2021-06-03 | 2021-08-10 | 哈尔滨工业大学 | Distributed high-frequency ground wave radar combined positioning method based on interior point method |
CN113238217B (en) * | 2021-06-03 | 2024-03-08 | 哈尔滨工业大学 | Distributed high-frequency ground wave radar joint positioning method based on interior point method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103487034B (en) | Method for measuring distance and height by vehicle-mounted monocular camera based on vertical type target | |
CN109807911B (en) | Outdoor patrol robot multi-environment combined positioning method based on GNSS, UWB, IMU, laser radar and code disc | |
CN103091660A (en) | Virtual base line direction-finding method based on phase jump correction | |
CN105547297A (en) | Indoor positioning method based on UWB positioning system | |
CN104965191A (en) | Two-site time difference positioning method | |
CN103376447A (en) | Method for achieving three-dimension positioning of non-cooperative double multistatic radar target | |
CN106443578A (en) | Indoor positioning base station coordinate calibration method | |
CN105044669A (en) | Three-station time-difference-measuring stereoscopic positioning method | |
CN101413785A (en) | Error compensation method of positioning system based on double-rotating laser plane transmitter network | |
CN105738915A (en) | Three-dimensional radar measuring method and three-dimensional radar measuring device | |
CN108594215A (en) | A kind of multi-point cooperative high precision wireless distance-measuring and positioning method | |
CN104363649A (en) | UKF (unscented Kalman filter)-based WSN (wireless sensor network) node location method with constraint conditions | |
CN102207380B (en) | High-precision horizontal axis tilt error compensation method | |
CN110516350A (en) | A kind of ERS point tolerance modification method based on anisotropy weighting | |
CN101847262A (en) | Fast three-dimensional point cloud searching and matching method | |
CN105277919A (en) | Single-point difference positioning method | |
Pop et al. | Verification of a triangulation network in Cluj-Napoca for future topographic surveys | |
CN102155956B (en) | High-precision method for compensating horizontal axis tilt error of vertical angle | |
CN110365420A (en) | A kind of underwater sound sensing network node position error modification method of combination velocity of sound inverting | |
CN111912431B (en) | Method for testing positioning accuracy of mobile robot navigation system | |
CN101266153A (en) | Mapping engineering top total station accuracy assessment method | |
CN102620745A (en) | Airborne inertial measurement unite (IMU) collimation axis error calibration method | |
CN109001670B (en) | Distributed passive positioning method and device combining time difference and angle | |
CN113702928B (en) | Intersection radar calibration method and equipment | |
CN104330078B (en) | Combined measuring method based on three-point resection model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20151007 |