CN110426672B - Double-iteration positioning system based on interval analysis and application thereof - Google Patents

Double-iteration positioning system based on interval analysis and application thereof Download PDF

Info

Publication number
CN110426672B
CN110426672B CN201910664911.7A CN201910664911A CN110426672B CN 110426672 B CN110426672 B CN 110426672B CN 201910664911 A CN201910664911 A CN 201910664911A CN 110426672 B CN110426672 B CN 110426672B
Authority
CN
China
Prior art keywords
source
box
iteration
positioning
mid
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.)
Active
Application number
CN201910664911.7A
Other languages
Chinese (zh)
Other versions
CN110426672A (en
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.)
Jiangnan University
Original Assignee
Jiangnan University
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 Jiangnan University filed Critical Jiangnan University
Priority to CN201910664911.7A priority Critical patent/CN110426672B/en
Publication of CN110426672A publication Critical patent/CN110426672A/en
Application granted granted Critical
Publication of CN110426672B publication Critical patent/CN110426672B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/10—Position of receiver fixed by co-ordinating a plurality of position lines defined by path-difference measurements, e.g. omega or decca systems

Abstract

The invention discloses a double-iteration positioning system based on interval analysis and application thereof, belonging to the technical field of target positioning. The invention provides a double-iteration positioning system based on interval analysis for positioning a moving source on the basis of arrival time difference and arrival frequency difference. Unlike the conventional gaussian error model, the position of the radiation source and the TDOA-FDOA measurement are calculated only over the envelope interval. The algorithm adopts a double-iteration strategy to alternately obtain a position box and a speed box containing a positioned moving source, and takes a midpoint as a point estimation value. The effectiveness and the accuracy of the method are shown by a simulation experiment.

Description

Double-iteration positioning system based on interval analysis and application thereof
Technical Field
The invention relates to a double-iteration positioning system based on interval analysis and application thereof, belonging to the technical field of target positioning.
Background
Time difference of arrival (TDOA) and frequency difference of arrival (FDOA) measurements have been widely used in passive positioning systems such as radar, navigation, etc. The location of the mobile source based on TDOA-FDOA measurements becomes very complex due to the non-linear relationship of TDOA and FDOA.
For many years, various algorithms have been proposed to solve this problem. In particular, Ho and Xu convert a nonlinear TDOA and FDOA equation set into a quasi-linear equation set by introducing a set of pseudo-linearity parameters, and a two-step weighted least squares algorithm is proposed (k. c.ho and w.xu, "An Accurate algebratic Solution for moving source Location Using TDOA and FDOA measures," IEEE trans. signal process, vol.52, No.9, pp.2453-2463, sep.2004.). The solution of this method is closed and the lower cramer limit is reached at low noise. Ho and Lu in turn extend the TS-WLS method for use in positioning systems where sensor position errors exist (k.c. Ho, x.lu and l.kovavisearch, "Source Localization Using TDOA and FDOA measures in the Presence of the sensor Localization errors: Analysis and Solution," IEEE trans.signal process, vol.55, No.2, pp.684-696, feb.2007.). Later, Foy proposed a Taylor series method whose main idea was to solve for the first order Taylor expansion around the initial measurement and solve the objective function by iteration (W.H. Foy, "Position-location solution by Taylor-series estimation," IEEE trans. Aerosp. Electron. Syst., vol.AES-12, pp.187C194, Mar.1976.). In the article by Zhu, he proposed a dual-iteration localization algorithm based on TODA-FDOA measurements that brings estimates of position sources into solution for velocity values, thereby reducing computational cost (G.Zhu and D.Feng, "Bi-iterative method for moving source localization using TDOA and FDOA measures," Electron.Lett. vol., vol.51, No.1, pp.8-10, Jan.2015.).
However, the above methods all assume gaussian error models and can only provide point estimates of the source position and cannot provide confidence intervals.
Disclosure of Invention
The invention aims to provide a double-iteration positioning method based on interval analysis, which is characterized in that a positioning system is used for positioning a moving source, the positioning system comprises a movable external radiation source for emitting signals and an observation station fixed at an origin position, two pairs of antennas are simultaneously distributed on the observation station, one pair of antennas is used for receiving direct signals from the external radiation source, and the other pair of antennas is used for receiving reflected signals of a positioning target source, and the method comprises the following steps:
the first process is as follows:
3) construction of TDOA and FDOA model mo＝[do,fo]T∈[mo]Initial position box [ u ] brought into the source of movement to be positioned]And initial velocity box
Figure RE-GDA0002205297060000021
And is provided with QrIs an identity matrix;
4) calculating a position box of the positioned mobile source, wherein the weighted least square estimation value of the position box is as follows:
Figure RE-GDA0002205297060000022
the second process:
3) the position box [ u ] of the positioned moving source calculated by the first process is usedo]And initial velocityBox
Figure RE-GDA0002205297060000023
Model m substituting TDOA and FDOAo＝[do,fo]T∈[mo]In and is provided with QrIs an identity matrix.
4) Calculating the velocity box of the positioned moving source, wherein the weighted least square estimation value is as follows:
Figure RE-GDA0002205297060000024
updating the location Box [ u ] of the initial Source of movement]Speed box
Figure RE-GDA0002205297060000025
And recalculating the matrix Q by the updated position bin and velocity binrAnd
Figure RE-GDA0002205297060000026
repeating the two processes until the modulus of the estimation difference of the previous and subsequent two times is less than a given threshold value to obtain the position estimation and the speed estimation of the mobile source;
wherein c is the speed of light,
Figure RE-GDA0002205297060000027
and
Figure RE-GDA0002205297060000028
is an interval measure of TDOA and FDOA by vector mo]To represent
Figure RE-GDA0002205297060000029
And
Figure RE-GDA00022052970600000210
set of (a) um＝mid([u])，
Figure RE-GDA00022052970600000211
Position box [ u ]]Speed of mixingMeasuring box
Figure RE-GDA00022052970600000212
The initial value of the positioning target is obtained by the maximum coverage area of the external radiation source and the maximum speed of the positioned moving source respectively, the initial speed box is set as the maximum speed of the positioning target, and the initial position box is set as the maximum coverage area of the external radiation source.
In one embodiment of the invention, 10 of an initial iteration speed bin and a position bin are set-6Multiple a given threshold.
In one embodiment of the invention, the true measurement value moAt um＝mid([u]) A first order Taylor expansion of the vicinity of
Figure RE-GDA00022052970600000213
Wherein
Figure RE-GDA00022052970600000214
While
Figure RE-GDA00022052970600000215
In one embodiment of the invention, uoIs a least squares estimate of
Figure RE-GDA00022052970600000216
Wherein QrIs the covariance matrix of the TDOA measurements.
In one embodiment of the present invention, let moIs [ m ]o]Then, then
Figure RE-GDA0002205297060000031
In one embodiment of the invention, J is used1,m＝mid([J1]) As [ J1]Is estimated value ofThen, the position box for positioning the moving source is:
Figure RE-GDA0002205297060000032
in one embodiment of the invention, it is assumed that the position of the mobile source being located is approximately uo m＝mid([uo]) Then the true measured value moIn that
Figure RE-GDA0002205297060000033
A first order Taylor expansion of the vicinity of
Figure RE-GDA0002205297060000034
Wherein
Figure RE-GDA0002205297060000035
And
Figure RE-GDA0002205297060000036
having a similar form of equation, otherwise umBy uo mIn the alternative,
jacobian matrix J2Is shown as
Figure RE-GDA0002205297060000037
In one embodiment of the invention, the velocity of the moving source being positioned
Figure RE-GDA0002205297060000038
Is a least squares estimate of
Figure RE-GDA0002205297060000039
Wherein
Figure RE-GDA00022052970600000310
Is a covariance matrix of FDOA measurements.
In one embodiment of the present invention, let moIs [ m ]o]，
Figure RE-GDA00022052970600000311
In one embodiment of the invention, J is used2,m＝mid([J2]) Approximate estimate [ J2]Then, then
Figure RE-GDA00022052970600000312
The invention has the beneficial effects that:
the invention provides a double-iteration positioning system based on interval analysis for positioning a moving source on the basis of arrival time difference and arrival frequency difference, which comprises a movable external radiation source for transmitting signals and an observation station fixed at an origin position, wherein two pairs of antennas are simultaneously distributed on the positioning system, one pair is used for receiving direct signals from the external radiation source, and the other pair is used for receiving reflected signals of a positioning target source. The observation station at the origin can simultaneously detect the signal from the moving external radiation source and the reflected signal from the moving external radiation source by the positioned moving source. The time difference of arrival of the two signals at the observation station exists, and an equation of the distance difference between the moving external radiation source and the positioned moving source can be constructed through the principle. Similarly, according to the Kepler frequency shift effect, the observation station can obtain a frequency equation between the positioned mobile source and the external emitting source. Unlike the conventional Gaussian error model, the position of the radiation source and the TDOA-FDOA measurements are calculated only by a closed solution. In the system, a double iteration strategy is adopted to alternately obtain a position box and a speed box containing a positioned moving source, and a midpoint is used as a point estimation value. The simulation experiment shows the stability and accuracy of the method.
Drawings
FIG. 1: deviation of the algorithm provided by the present application, the two-step weighted least squares algorithm, from the maximum likelihood algorithm when the located mobile source is located at (280,325,275) meters.
FIG. 2: RMSE of the algorithm, two-step weighted least squares algorithm and maximum likelihood algorithm provided herein when the located mobile source is located at (280,325,275) meters.
FIG. 3: the proposed algorithm estimates the position boundaries.
FIG. 4: the proposed algorithm estimates the velocity boundary.
Detailed Description
Embodiment 1 double-iteration positioning method based on interval analysis
Introducing variables: assuming that there are M moving external radiation sources in a given three-dimensional space, their positions and velocities are respectively represented as
Figure RE-GDA0002205297060000041
And
Figure RE-GDA0002205297060000042
the 1 observation station is fixed at the origin position, and two pairs of antennas are simultaneously distributed on the observation station, wherein one pair of antennas is used for receiving direct signals from an external radiation source, and the other pair of antennas is used for receiving reflected signals of a positioning target source. The positioned target source is randomly deployed at any position u in the spaceo＝[xo,yo,zo]And at a speed
Figure RE-GDA0002205297060000043
And (6) moving. Based on TDOA and FDOA measurement signals acquired by the observation station, the time difference between the jth mobile radiation source signal and the positioning target source reflection signal reaching the observation station can be obtained as
Figure RE-GDA0002205297060000044
Doppler frequency difference of
Figure RE-GDA0002205297060000045
Figure RE-GDA0002205297060000046
Since the propagation speed of the signal is fixed, the distance difference and the speed difference,
Figure RE-GDA0002205297060000047
Figure RE-GDA0002205297060000048
wherein
Figure RE-GDA0002205297060000049
Indicating the distance of the jth external radiation source from the positioning target source,
Figure RE-GDA00022052970600000410
is differentiated with respect to time to obtain
Figure RE-GDA00022052970600000411
| | represents the euclidean distance, and ρ represents the unitized operation of the vector. Ro＝||uo| represents the distance from the observation station from which the positioning target originates,
Figure RE-GDA00022052970600000412
is composed of
Figure RE-GDA00022052970600000413
The result of the derivation of the time is obtained,
Figure RE-GDA00022052970600000414
to locate the range of the target from the observation station,
Figure RE-GDA00022052970600000415
the same can be obtained. c represents the signal propagation speed, approximately the speed of light.
The signals will receive various interferences during the transmission process, so the measured signals obtained by TDOA and FDOA have errors, and d in the measurement process is determined by assuming that the measurement errors follow a Gaussian distributionjoAnd fjoIs equal to
Figure RE-GDA0002205297060000051
Figure RE-GDA0002205297060000052
Wherein deltadAnd deltafObeying a gaussian distribution.
According to the 3 delta principle, the true measured value has more than 96% possibility of being positioned in the range of +/-3 delta of the average value, and the measured value can be approximately estimated
Figure RE-GDA0002205297060000053
And
Figure RE-GDA0002205297060000054
included within a range of ± 3 δ, averaged with actual measurements.
Figure RE-GDA0002205297060000055
Figure RE-GDA0002205297060000056
Wherein [ Delta d]And [ Delta d ]]Being the gaussian error boundary of TDOA and FDOA,
Figure RE-GDA0002205297060000057
and
Figure RE-GDA0002205297060000058
is an interval measurement value of TDOA and FDOA. For symbol simplification, use the vector mo]To represent
Figure RE-GDA0002205297060000059
And
Figure RE-GDA00022052970600000510
a set of [ m ]o]＝[[do]T,[fo]T]TWherein
Figure RE-GDA00022052970600000511
True TDOA and FDOA measurements mo＝[do,fo]T∈[mo]。
The main purpose of the following algorithm is to measure [ m ] by the interval of TDOA and FDOAo]And obtaining a position box and a speed box containing the real position and the speed value of the positioned target source.
(II) TDOA and FDOA dual-iteration positioning algorithm flow based on interval analysis:
the first step process:
1) construction of TDOA and FDOA models (m)o＝[do,fo]T∈[mo]) The initial position box brought into the positioned moving source is [ u ]]And initial velocity box
Figure RE-GDA00022052970600000512
And is provided with QrIs an identity matrix.
2) Calculating the velocity box of the positioned moving source, and then weighting the least square estimation value:
Figure RE-GDA00022052970600000513
the specific derivation process is as follows.
The second step is as follows:
5) the position box [ u ] of the positioned moving source calculated in the first step is usedo]And initial velocity box
Figure RE-GDA00022052970600000515
Model m substituting TDOA and FDOAo＝[do,fo]T∈[mo]In and is provided with
Figure RE-GDA00022052970600000514
Is a unit array.
6) Calculating the velocity box of the positioned moving source, wherein the weighted least square estimation value is as follows:
Figure RE-GDA0002205297060000061
the specific derivation process is referred to as derivation two in the following chapter.
Updating the location Box [ u ] of the initial Source of movement]Speed box
Figure RE-GDA0002205297060000062
And recalculates the matrix Q by the position bin and the velocity binrAnd
Figure RE-GDA0002205297060000063
repeating the above two processes until the modulus of the estimated difference between the previous and next two times is less than a given threshold or the iteration number is reached (the specific threshold and the iteration number are related to the environmental factors), and considering the environmental factors, the method selects 10 of the initial iteration speed box and the position box-6The multiple is a given threshold, and the iteration times are converted according to the given threshold. The application selects 10 with threshold values of initial iteration speed box and position box-6And stopping the iteration process and obtaining the position and the speed of the moving source.
The algorithm firstly adopts a strategy of a double iteration positioning algorithm, namely, estimation operation is carried out on the position and the speed of a positioned target source alternately. Velocity box [ u ]]Speed box
Figure RE-GDA0002205297060000064
Can be obtained by a priori knowledge, for example, the maximum coverage area of the external radiation source and the maximum velocity of the moving source being positioned. Finally obtained position box [ u ]o]Speed box
Figure RE-GDA0002205297060000065
Is included in the initial value u]And
Figure RE-GDA0002205297060000066
in (1).
Derivation one: from the formulas (2), (3), (5), (6), moAnd a mobile source position box [ u ]o]In a non-linear relationship. First assume [ u ]]Equivalent to um＝mid([u]) Then the true measured value moAt um＝mid([u]) A first order Taylor expansion of the vicinity of
Figure RE-GDA0002205297060000067
Wherein
Figure RE-GDA0002205297060000068
While
Figure RE-GDA0002205297060000069
Jacobian matrix J1Is shown as
Figure RE-GDA00022052970600000610
According to the formula (8), uoIs a least squares estimate of
Figure RE-GDA00022052970600000611
Wherein QrIs the covariance matrix of the TDOA measurements.
It can be derived from the above formula that u can be obtained only by acquiring the true position of the external radiation source and the accurate TDOA and FDOA measured valuesoThe value is obtained. Because m iso∈[mo]So that m in the formula (8)oIs [ m ]o]At this time, it can be
Figure RE-GDA00022052970600000612
When solving equation (10), it is necessary to calculate an interval matrix
Figure RE-GDA0002205297060000071
This would greatly increase the computational complexity. Therefore, this application adopts J1,m＝mid([J1]) As [ J1]The position box of the positioned moving source can be obtained by the following formula:
Figure RE-GDA0002205297060000072
derivation two: by equation (11), it can be assumed that the position of the positioned mobile source can be approximated as uo m＝mid([uo]) Then the true measured value moIn that
Figure RE-GDA0002205297060000073
A first order Taylor expansion of the vicinity of
Figure RE-GDA0002205297060000074
Wherein
Figure RE-GDA0002205297060000075
And
Figure RE-GDA0002205297060000076
having a similar form of equation, otherwise umBy uo mAnd (6) replacing.
Jacobian matrix J2Is shown as
Figure RE-GDA0002205297060000077
From equation (12), the velocity of the moving source to be located can be obtained
Figure RE-GDA0002205297060000078
Is a least squares estimate of
Figure RE-GDA0002205297060000079
Wherein
Figure RE-GDA00022052970600000710
Is a covariance matrix of FDOA measurements.
And solve for [ uo]The process of (a) is similar to (b), moIs replaced by [ mo]The following formula can be obtained
Figure RE-GDA00022052970600000711
Then further using J2,m＝mid([J2]) Approximate estimate [ J2]Then, then
Figure RE-GDA00022052970600000712
Through the derivation, the algorithm of the application adopts an initial position box u]Speed box
Figure RE-GDA00022052970600000713
As an initial input. First assume the velocity of the moving source being located
Figure RE-GDA00022052970600000714
It is known that the location box of the located mobile source can be updated using equation (11). The velocity box of the positioned moving source can then be updated by equation (16). Repeating the above two steps until the widths and the middle points of the position box and the speed box converge. The point estimate for the position and velocity of the moving source being located can be approximated by the midpoint of the bin.
Example 2 simulation experiment
The experiment assumes that the positioning system consists of 5 mobile external radiation sources, the radiation source positions being given in table 1. In the simulations that follow, we define the position of the positioned moving source in m and the velocity in m/s. For experimental convenience, the experiment assumed that the TDOA and FDOA signals were independent of each other, and the measured values of the signals were generated by adding true values to white Gaussian noise. Thus measuring TDOA andFDOA covariance matrix is δ2Q and 0.1. delta2Q, wherein Q is 0.5(I + I)TI) I is the identity matrix and δ is the noise variance. As the measurement error increases, the estimation bias and accuracy vary.
TABLE 1 radiation source location
Figure RE-GDA0002205297060000081
The present application performs 5000 Monte Carlo experiments, and for algorithm comparisons, defines the positional deviation of the positioned mobile source as
Figure RE-GDA0002205297060000082
Similarly positioned the moving source velocity offset is positioned as
Figure RE-GDA0002205297060000083
The position accuracy based on the root mean square error is defined as
Figure RE-GDA0002205297060000084
The position accuracy is defined as
Figure RE-GDA0002205297060000085
By varying the value of delta, variations in the measurement that cause noise can be achieved. With the increase of delta, the positioning deviation and the positioning precision of the algorithm and the two comparison algorithms are increased. However, as can be seen from fig. 1 and fig. 2, the algorithm provided by the present application can still maintain superior positioning deviation and positioning accuracy under high noise compared to the two-step weighted least square algorithm and the maximum likelihood algorithm.
The superiority of the algorithm relative to the point estimation result is verified by comparing with a two-step weighted least square algorithm and a maximum likelihood estimation algorithm. In addition, the point estimation and covariance matrix of the motion source can be obtained by using the least square method. The covariance matrix is then used to calculate the standard deviation of the position and velocity of the moving source being located. Finally, the position and velocity bins of the positioned moving source are obtained according to the 3-sigma principle. Fig. 3 and 4 show the positioning results of the position bin and the velocity bin obtained by the proposed algorithm. It can be seen that the upper and lower bounds of the position estimate and velocity estimate are very close to the true values of the located mobile source position and velocity before the noise power is 30 DB. When the measurement noise is greater than 30DB, the upper and lower bounds of the velocity and position estimates deviate from their true values, but their velocity and position estimation boxes still surround the true values. This shows that the algorithm can provide a range of position and velocity estimates, including real values, even under high noise conditions.
Although the present invention has been described with reference to the preferred embodiments, it should be understood that various changes and modifications can be made therein by those skilled in the art without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (10)

1. A double iteration positioning method based on interval analysis is characterized in that a positioning system is used for positioning a moving source, the positioning system comprises a movable external radiation source for emitting signals and an observation station fixed at an origin position, two pairs of antennas are simultaneously distributed on the observation station, one pair of antennas is used for receiving direct signals from the external radiation source, and the other pair of antennas is used for receiving reflected signals of a positioning target source, and the method comprises the following steps:
the first process is as follows:
1) construction of TDOA and FDOA model mo＝[do,fo]T∈[mo]Initial position box [ u ] brought into the source of movement to be positioned]And initial velocity box
Figure FDA0002858820340000011
And is provided with QrIs an identity matrix;
2) calculating a position box of the positioned mobile source, wherein the weighted least square estimation value of the position box is as follows:
Figure FDA0002858820340000012
wherein, J1,m＝mid([J1]) Is [ J1]Is determined by the estimated value of (c),
Figure FDA0002858820340000013
while
Figure FDA0002858820340000014
Figure FDA0002858820340000015
Is a Jacobian matrix;
the second process:
1) the position box [ u ] of the positioned moving source calculated by the first process is usedo]And initial velocity box
Figure FDA0002858820340000016
Model m substituting TDOA and FDOAo＝[do,fo]T∈[mo]In and is provided with
Figure FDA0002858820340000017
Is an identity matrix;
2) calculating the velocity box of the positioned moving source, wherein the weighted least square estimation value is as follows:
Figure FDA0002858820340000018
wherein, J2,m＝mid([J2]) Is [ J2]An estimated value of (d);
Figure FDA0002858820340000019
is a Jacobian matrix;
updating the location Box [ u ] of the initial Source of movement]Speed box
Figure FDA00028588203400000110
And recalculating the matrix Q by the updated position bin and velocity binrAnd
Figure FDA00028588203400000111
repeating the two processes until the modulus of the estimation difference of the previous and subsequent two times is less than a given threshold value to obtain the position estimation and the speed estimation of the mobile source;
wherein c is the speed of light,
Figure FDA00028588203400000112
and
Figure FDA00028588203400000113
is an interval measure of TDOA and FDOA by vector mo]To represent
Figure FDA00028588203400000114
And
Figure FDA00028588203400000115
set of (a) um＝mid([u])，
Figure FDA00028588203400000116
Position box [ u ]]Speed box
Figure FDA00028588203400000117
Is obtained by the maximum coverage area of the external radiation source and the maximum speed of the positioned moving source, respectively.
2. Double iterative localization method according to claim 1, characterized in that the true measurement value m isoAt um＝mid([u]) A first order Taylor expansion of the vicinity of
Figure FDA0002858820340000021
Wherein
Figure FDA0002858820340000022
While
Figure FDA0002858820340000023
Jacobi matrix
Figure FDA0002858820340000024
3. The dual-iteration localization method of claim 2, wherein u isoIs a least squares estimate of
Figure FDA0002858820340000025
Wherein QrIs the covariance matrix of the TDOA measurements.
4. A dual-iteration localization method according to claim 3, characterized by letting m beoIs replaced by [ mo]Then, then
Figure FDA0002858820340000026
5. The dual-iteration localization method of claim 4, wherein J is employed1,m＝mid([J1]) As [ J1]The position box to be located with the moving source is:
Figure FDA0002858820340000027
6. the dual-iteration localization method of claim 1, wherein the position of the localized mobile source is assumed to be approximately uo m＝mid([uo]) Then the true measured value moIn that
Figure FDA0002858820340000028
A first order Taylor expansion of the vicinity of
Figure FDA0002858820340000029
Wherein
Figure FDA00028588203400000210
And
Figure FDA00028588203400000211
having a similar form of equation, otherwise umBy uo mIn the alternative,
jacobian matrix J2Is shown as
Figure FDA00028588203400000212
7. The dual-iteration positioning method of claim 6, wherein the moving source velocity being positioned
Figure FDA00028588203400000213
Is a least squares estimate of
Figure FDA00028588203400000214
Wherein
Figure FDA00028588203400000215
Is a covariance matrix of FDOA measurements.
8. The dual-iteration positioning method of claim 7, wherein let moIs replaced by [ mo]，
Figure FDA0002858820340000031
9. The dual-iteration localization method of claim 8, wherein J is used2,m＝mid([J2]) Approximate estimate [ J2]Then, then
Figure FDA0002858820340000032
10. The method of any one of claims 1-9, setting 10 an initial iteration speed bin and a position bin-6Multiple a given threshold.
CN201910664911.7A 2019-07-23 2019-07-23 Double-iteration positioning system based on interval analysis and application thereof Active CN110426672B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910664911.7A CN110426672B (en) 2019-07-23 2019-07-23 Double-iteration positioning system based on interval analysis and application thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910664911.7A CN110426672B (en) 2019-07-23 2019-07-23 Double-iteration positioning system based on interval analysis and application thereof

Publications (2)

Publication Number Publication Date
CN110426672A CN110426672A (en) 2019-11-08
CN110426672B true CN110426672B (en) 2021-05-04

Family

ID=68411850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910664911.7A Active CN110426672B (en) 2019-07-23 2019-07-23 Double-iteration positioning system based on interval analysis and application thereof

Country Status (1)

Country Link
CN (1) CN110426672B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111239718B (en) * 2020-01-17 2022-11-15 电子科技大学 Multi-base-station target positioning method based on single-satellite radiation source
CN113835061B (en) * 2021-08-13 2022-07-05 中国人民解放军战略支援部队信息工程大学 Single-platform Doppler two-stage closed positioning method in presence of signal carrier frequency prior error
CN113945888B (en) * 2021-10-19 2022-05-03 江南大学 Interval passive positioning method and system based on TDOA

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102331581A (en) * 2011-05-27 2012-01-25 哈尔滨工业大学 Rapid positioning method of binary TDOA/FDOA satellite-to-earth integration positioning system
CN107703482A (en) * 2017-10-20 2018-02-16 电子科技大学 The AOA localization methods that a kind of closed solutions are combined with iterative algorithm
CN108322933A (en) * 2018-01-24 2018-07-24 江南大学 A kind of node accelerating gradient decline localization method
CN109633581A (en) * 2018-12-26 2019-04-16 杭州电子科技大学 Based on the localization method under external sort algorithm TDOA/FDOA error correction
CN109633591A (en) * 2019-01-18 2019-04-16 杭州电子科技大学 External illuminators-based radar is biradical away from localization method under a kind of observation station location error

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7688264B2 (en) * 2007-11-27 2010-03-30 The Boeing Company Method and apparatus for passive single platform geo-location

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102331581A (en) * 2011-05-27 2012-01-25 哈尔滨工业大学 Rapid positioning method of binary TDOA/FDOA satellite-to-earth integration positioning system
CN107703482A (en) * 2017-10-20 2018-02-16 电子科技大学 The AOA localization methods that a kind of closed solutions are combined with iterative algorithm
CN108322933A (en) * 2018-01-24 2018-07-24 江南大学 A kind of node accelerating gradient decline localization method
CN109633581A (en) * 2018-12-26 2019-04-16 杭州电子科技大学 Based on the localization method under external sort algorithm TDOA/FDOA error correction
CN109633591A (en) * 2019-01-18 2019-04-16 杭州电子科技大学 External illuminators-based radar is biradical away from localization method under a kind of observation station location error

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
An approximately efficient bi-iterative method for source position and velocity estimation using TDOA and FDOA measurements;Guo-Hui Zhu等;《Signal Processing》;20160101;全文 *
一种改进的双站时频差目标跟踪高斯和滤波算法;曹亚琴等;《计算机与数字工程》;20180531;第46卷(第5期);全文 *
利用TDOA 和FDOA 的单站多外辐射源目标定位算法;赵勇胜等;《四川大学学报( 工程科学版)》;20160430;第48卷(第1期);全文 *

Also Published As

Publication number Publication date
CN110426672A (en) 2019-11-08

Similar Documents

Publication Publication Date Title
CN110426672B (en) Double-iteration positioning system based on interval analysis and application thereof
CN110174643B (en) Positioning method based on arrival time difference without noise power information
CN107820206B (en) Non-line-of-sight positioning method based on signal intensity
CN107526070A (en) The multipath fusion multiple target tracking algorithm of sky-wave OTH radar
CN108761387B (en) Double-station time difference and frequency difference combined positioning method for fixed radiation source
KR20180091372A (en) Method for tracking target position of radar
CN110636436A (en) Three-dimensional UWB indoor positioning method based on improved CHAN algorithm
CN110471029B (en) Single-station passive positioning method and device based on extended Kalman filtering
Yang et al. Kalman filter based localization and tracking estimation for HIMR RFID systems
Hao et al. A novel Taylor series method for source and receiver localization using TDOA and FDOA measurements with uncertain receiver positions
Kazemi et al. Efficient closed-form solution for 3-D hybrid localization in multistatic radars
Koledoye et al. Improved MDS-based localization with non-line-of-sight RF links
CN111624549B (en) Passive filtering tracking method under non-common-view condition
Chen et al. TDOA/FDOA mobile target localization and tracking with adaptive extended Kalman filter
Molnár et al. Development of an UWB based indoor positioning system
Nicolalde-Rodríguez et al. Robust passive coherent location via nonlinearly constrained least squares
CN110536410B (en) Positioning method based on RSS and TDOA measurement in non-line-of-sight environment
Amiot et al. A hybrid positioning method based on hypothesis testing
Chen et al. IMM tracking of a 3D maneuvering target with passive TDOA system
CN112835020B (en) Rigid body positioning method for non-line-of-sight parameter estimation
US10156630B2 (en) Method for passively locating a non-movable transmitter
CN112068099B (en) Multi-radiation-source rapid positioning and speed measuring method and device based on error compensation
CN109884582B (en) Method for rapidly determining three-dimensional coordinates of target by utilizing one-dimensional direction finding
CN110401915B (en) SEKF and distance reconstruction combined moving target positioning method under NLOS condition
CN113608165A (en) Multi-station passive positioning method based on signal arrival time difference

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant