CN110412504A - It is associated with based on angle with the passive track-corelation of time difference information and localization method - Google Patents

It is associated with based on angle with the passive track-corelation of time difference information and localization method Download PDF

Info

Publication number
CN110412504A
CN110412504A CN201910738048.5A CN201910738048A CN110412504A CN 110412504 A CN110412504 A CN 110412504A CN 201910738048 A CN201910738048 A CN 201910738048A CN 110412504 A CN110412504 A CN 110412504A
Authority
CN
China
Prior art keywords
target
observation station
angle
error
cost matrix
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
CN201910738048.5A
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.)
University of Electronic Science and Technology of China
Original Assignee
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201910738048.5A priority Critical patent/CN110412504A/en
Publication of CN110412504A publication Critical patent/CN110412504A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/0257Hybrid positioning
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/10Position of receiver fixed by co-ordinating a plurality of position lines defined by path-difference measurements, e.g. omega or decca systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/14Determining absolute distances from a plurality of spaced points of known location

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

The invention belongs to electronic countermeasure technology fields, and in particular to a kind of to be associated with based on angle with the passive track-corelation of time difference information and localization method.Method of the invention mainly uses multi radiation sources data association algorithm, is associated to obtain cost matrix C by the angular observation data of multiple observation stations and time difference observation data simultaneously respectively1(p, q) and C2(p, q), then it is utilized respectively angle measurement and the standard deviation of range difference measurement value error obtains total cost matrix C (p, q), C'=1./C is enabled, the position by solving the local maximum that C' corresponds to number of targets positions target.Beneficial effects of the present invention are that the present invention accurately can complete to be associated with to multistation data, and finally accurately estimate the position of target, and method is simple, works well.

Description

It is associated with based on angle with the passive track-corelation of time difference information and localization method
Technical field
The invention belongs to electronic countermeasure technology fields, and in particular to a kind of passive multistation based on angle and time difference information is more Target association and localization method.
Background technique
In electronic countermeasure field, due to the development of perturbation technique, active location technology poor anti jamming capability lacks hidden Property weakness start gradually to become obvious, and passive location technology has very strong superiority relative to active location technology. It has many advantages, such as remote operating distance, lobe-on-receive, strong antijamming capability and is not easy to be realized by other side, for improving system in electricity Survival ability and fighting efficiency under son war environment have a very important role.Passive location mainly passes through radiation emission Electromagnetic wave signal reaches the time difference of observation station obtaining the azimuth of radiation source, pitch angle or electromagnetic wave and realizes.Therefore Passive location generally requires multiple base station collaborations and completes.But when there are multiple targets, passive multistation location is easy to appear void Decoy point causes to position mistake.So carrying out data correlation for Multiple targets observation data, false target point is rejected for standard Determine that position is of great significance.
In passive location system, direction cross positioning is presented as that a plurality of straight line crosses in space, be with most extensively, A kind of longest localization method of search time.But meanwhile it there is also some problems, such as in a practical situation it measurement miss Difference is larger, can have biggish position error particularly with long distance positioning;If it is in the case where carrying out Multi-target position, Multiple observation stations can generate a large amount of false crosspoint.And positioning using TDOA is presented as that multiple hyperboloids cross in space, with it Anti-interference, the advantages that Yi Shixian and precision are high, occupies more and more important position in passive location.Positioning using TDOA is to pass through Measurement emitter Signals reach the time difference of different observation stations to determine target position, need multiple base stations to cooperate with, therefore also hold Easily there is false target point.So angle information and time difference information joint are carried out data correlation, for realizing that target is accurately fixed Position has great research significance.
Summary of the invention
In view of the above-mentioned problems, being closed the invention proposes a kind of based on the passive track-corelation of angle information and time difference information Connection and localization method.
The technical solution adopted by the present invention is that:
It is a kind of that the positioning of track-corelation is realized based on being associated with localization method for angle information and time difference information, it is described Be in the positioning that crossed using angle and the straight line of the progress of time difference information and hyperbola, to false target reject to Realize the data correlation of multi-source data.The track-corelation location model of this method is as shown in Figure 1.
It has been known that there is 3 observation stations, and wherein observation station 1 is main website, and observation station 2 and observation station 3 are slave station, their coordinate point Not are as follows: (x1,y1)、(x2,y2)、(x3,y3), the coordinate of target 1 and target 2 is unknown.Assuming that 3 observation stations observe target respectively 1 azimuth is θ1,1、θ2,1、θ3,1If the coordinate of radiation source 1 is (xp,yq), then by formula:
When only existing single target, the coordinate of radiation source 1 can be found out by azimuth information, wherein niIt is i-th The angle measurement error of observation station.
If being positioned using time difference information to single target, it may be assumed that when known to 1 position of target, target 1 reaches The distance between main website 1 and slave station 2,3 difference are respectively △ d21、△d31, according to hyp definition as can be seen that target 1 is located at It is fixed point with main website and slave station, is that fixed length is formed by hyperbola with the distance between the fixed point of target 1 to two difference.When target position When setting unknown, it is assumed that the coordinate of unknown object 2 is (x in Fig. 1p,yq), with main website 1, the distance of slave station 2 and slave station 3 is distinguished For Indicate target 2 reach the distance between main website and slave station k it is poor, then just like Lower equation:
When being related to single goal positioning, above equation can directly calculate the position (x of unknown object 2p,yq), target position Set be two-by-two observation station the intersection point on hyperbola is formed by using range difference as fixed length.
Above method can only be positioned with single target, when there are multiple targets, each observation The problem of azimuth information or time difference information of multiple targets can all be observed by standing, this relates to data correlation.For example sentence Which angle information is the corresponding same target in disconnected each observation station out, could more accurately be positioned in this way, otherwise The influence of false target may be will receive, to cause positioning mistake.Since angle measurement error is opposite in actual measurement It is larger, if will cause biggish position error in the farther away situation in target range observation station, therefore, in order to ensure that target positions Accuracy, the present invention use angle information and the united method of time difference information, association algorithm flow chart as shown in Fig. 2, tool Body method is as follows:
Assuming that observation station and target are all located on X/Y plane, if main website position coordinates (x1,y1), the coordinate of slave station is (xk, yk), wherein k=2,3 ..., K, K are observation station sum.Then there is j-th of target to the distance of k-th observation station and arrives main website Apart from its difference are as follows:Wherein k=2,3 ..., K, j=1,2 ..., J, J be target sum,Table Show measurement error.It enablesFor the azimuth information of j-th of target of k-th of observation station measurement, wherein k=1, 2 ..., K, j=1,2 ..., J, θk,jIndicate k-th of observation station to j-th of target azimuth information true value, nk,jTable Show measurement noise.
Objective plane is divided into the grid of P × Q range, each mesh point represents a position coordinates in objective plane (xp,yq), wherein p=1,2 ..., P, q=1,2 ... Q traverses each mesh point in grid plan:
1, point (x is calculatedp,yq) azimuth relative to each observation station:Meanwhile Calculate this for the distance of main website and slave station and they the distance between it is poor:Its Middle k=2,3 ..., K.
2, each Searching point (x is calculatedp,yq) azimuth relative to observation stationThe side observed with observation station Parallactic angleBetween errorWherein k=1,2 ..., K, j=1,2 ..., J.
Range difference of this relative to observation station is calculated simultaneouslyThe range difference observed with observation stationBetween ErrorWherein k=2,3 ..., K, j=1,2 ..., J.
3, the cost matrix being made of the angle overall error searched for every time is calculatedWherein K=1,2 ..., K, p=1,2 ..., P, q=1,2 ..., Q and the cost matrix being made of range difference overall errorWherein k=2,3 ..., K, p=1,2 ..., P, q=1,2 ..., Q.
4, combine the cost matrix C based on angle1With the cost matrix C based on range difference2Obtain total cost matrix C (p, Q)=1/ σ1*C1+1/σ2*C2, wherein σ1And σ2Respectively represent the standard of observation station angle measurement and range difference measurement value error Difference.
5, C'=1./C is enabled, the position by solving the local maximum that C' corresponds to number of targets positions target.
Beneficial effects of the present invention are that can accurately be associated to multistation data, and finally accurately estimate target Position, method is simple, works well.
Detailed description of the invention
Fig. 1 is multistation angle and time difference Model for joint positioning figure;
Fig. 2 is multi radiation sources data association algorithm flow chart;
Fig. 3 is the Multi-target position puppet spectral peak figure based on angle;
Fig. 4 is the Multi-target position puppet spectral peak figure based on angle and the time difference;
Fig. 5 is positioning result figure;
Fig. 6 is correct localization with threshold variation curve graph.
Specific embodiment
The present invention is described in detail below with reference to specific example:
It is fixed that this example intersects the above-mentioned passive track-corelation direction finding based on angle information and time difference information using matlab Position method is verified, and for simplicity, is made the following assumptions to algorithm model:
1. main website, slave station and target are all in X/Y plane;
2. main website and slave station range accuracy having the same;
3. all Engineering Errors are all added in range error;
4. assuming that target is static or movement velocity is extremely low.
Assuming that target area is the rectangular region of 100km × 80km, 500 Monte Carlo emulation, 4 mesh are carried out Target coordinate is respectively (15,70), (35,75), (58,75), (85,70), unit km.Use 1 main website and 3 slave stations pair Above-mentioned target is positioned, and the coordinate of main website and slave station is respectively (5,30), (30,15), (70,15), (95,30), and unit is Km, the range error at all stations obey the Gaussian Profile that mean value is zero, and mutually indepedent between the range error at all stations.500 Multistation data correlation effect and locating effect are obtained after secondary Monte Carlo emulation.
Track-corelation data correlation effect:
To compare direction cross positioning method based on angle information and being closed based on the target of angle information and time difference information Connection localization method, in addition sets (10,50) for target, (13,50), (50,30), (85,60), wherein the first two target position Relatively.It is to be obtained after the angle information of each target observed to multiple observation stations is associated as shown in Figure 3 Target puppet spectral peak figure, as can be seen from the figure since two of them target range closely causes very much corresponding two spectrums cutting edge of a knife or a sword to become one Spectrum cutting edge of a knife or a sword to can not be accurately positioned, Fig. 4 be in situation identical with the target position Fig. 3 to the angle information of each target and when What poor information obtained after being associated has the pseudo- spectral peak figure of 4 peak values, Fig. 3 is compared, in Fig. 4 corresponding to 4 target positions It is high-visible to compose cutting edge of a knife or a sword, it follows that the object localization method based on angle information and time difference information is excellent for the resolution ratio of distance In the direction cross positioning method based on angle information.
For the locating effect for more intuitively showing algorithm, Fig. 5 gives the 100th locating effect.In figure, border circular areas Correct localization region is represented, estimated location is located in border circular areas then identification and correct positioning has been carried out, can from Fig. 5 To find out, estimated location and target actual position very close to, it is believed that accurate target position estimation may be implemented in the algorithm.
Error minimum value, max value of error and average error representated by three curves respectively indicate multiple in Fig. 6 The minimum value of the distance between target state estimator position and actual position, maximum value and to multiple target state estimator positions and true position The distance between set the average value that summation is averaged.Thresholding is the judgment criteria of target positioning correctness, when target is estimated It indicates that positioning is correct when counting position and actual position distance lower than thresholding, and then indicates to position mistake beyond thresholding.Accuracy is Refer to and positions the percentage that correct number accounts for total degree in 500 positioning carried out in the case where some determines thresholding.It can be with by Fig. 6 Find out, correct localization is gradually promoted with the increase of thresholding, the energy when actual position and Estimated Position Error are in 3km or so Reach very high correct localization, it can be seen that the passive track-corelation direction finding based on angle information and time difference information is handed over Fork localization method works well.

Claims (1)

1. being associated with based on angle with the passive track-corelation of time difference information and localization method, observation station and target is enabled to be all located at XY In plane, main website position coordinates (x1,y1), the coordinate of slave station is (xk,yk), wherein k=2,3 ..., K, K are observation station sum; There is distance of j-th of target to k-th of observation stationWith the distance for arriving main websiteDifference are as follows:Wherein J=1,2 ..., J, J are target sums,It indicates measurement error, enablesFor the jth of k-th of observation station measurement The azimuth information of a target, wherein k=1,2 ..., K, j=1,2 ..., J, θk,jIndicate k-th of observation station to j-th of mesh The true value of target azimuth information, nk,jIndicate measurement noise.Objective plane is divided into the grid of P × Q range, each net Lattice point represents a position coordinates (x in objective planep,yq), wherein p=1,2 ..., P, q=1,2 ... Q, traversal grid are flat Each mesh point in face;It is characterized in that, the described method comprises the following steps:
S1, Searching point (x is calculatedp,yq) azimuth relative to each observation station:Meanwhile it calculating This for the distance of main website and slave station and they the distance between it is poor:
S2, each Searching point (x is calculatedp,yq) azimuth relative to observation stationThe azimuth observed with observation stationBetween errorRange difference of this relative to observation station is calculated simultaneouslyWith The range difference that observation station is observedBetween error
The cost matrix that S3, calculating are made of the angle overall error searched for every timeAnd by away from The cost matrix of deviation overall error composition
S4, combine the cost matrix C based on angle1With the cost matrix C based on range difference2Obtain total cost matrix C (p, q) =1/ σ1*C1+1/σ2*C2, wherein σ1And σ2Respectively represent the standard deviation of observation station angle measurement and range difference measurement value error;
S5, C'=1./C is enabled, the position by solving the local maximum that C' corresponds to number of targets positions target.
CN201910738048.5A 2019-08-12 2019-08-12 It is associated with based on angle with the passive track-corelation of time difference information and localization method Pending CN110412504A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910738048.5A CN110412504A (en) 2019-08-12 2019-08-12 It is associated with based on angle with the passive track-corelation of time difference information and localization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910738048.5A CN110412504A (en) 2019-08-12 2019-08-12 It is associated with based on angle with the passive track-corelation of time difference information and localization method

Publications (1)

Publication Number Publication Date
CN110412504A true CN110412504A (en) 2019-11-05

Family

ID=68367061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910738048.5A Pending CN110412504A (en) 2019-08-12 2019-08-12 It is associated with based on angle with the passive track-corelation of time difference information and localization method

Country Status (1)

Country Link
CN (1) CN110412504A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308418A (en) * 2020-03-10 2020-06-19 慧众行知科技(北京)有限公司 Steady method for two-dimensional positioning of target with unknown height
CN111337875A (en) * 2020-04-13 2020-06-26 北京航天长征飞行器研究所 Space multi-station multi-system dimension reduction positioning method
CN113390406A (en) * 2021-06-16 2021-09-14 电子科技大学 Multi-target data association and positioning method based on passive multi-sensor system
CN115524662A (en) * 2022-10-27 2022-12-27 中国电子科技集团公司信息科学研究院 Direction finding time difference combined positioning method and system, electronic equipment and storage medium
RU2790347C1 (en) * 2022-04-14 2023-02-16 Общество с ограниченной ответственностью "Специальный Технологический Центр" Difference-range method for determining the location of objects

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108061877A (en) * 2017-12-14 2018-05-22 电子科技大学 A kind of passive track-corelation direction cross positioning method based on angle information
CN108089148A (en) * 2017-12-14 2018-05-29 电子科技大学 A kind of passive track-corelation direction cross positioning method based on time difference information
CN108152789A (en) * 2018-01-03 2018-06-12 电子科技大学 Utilize the passive track-corelation data correlation and localization method of RSS information
CN108802674A (en) * 2018-07-19 2018-11-13 中国人民解放军战略支援部队信息工程大学 It is a kind of for the combined method for searching and device that directly position
CN109001670A (en) * 2018-06-01 2018-12-14 中国人民解放军战略支援部队信息工程大学 A kind of distributed passive location method and device for combining the time difference and angle
CN109917333A (en) * 2019-01-31 2019-06-21 西南电子技术研究所(中国电子科技集团公司第十研究所) Merge the passive location method of AOA observed quantity and TDOA observed quantity

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108061877A (en) * 2017-12-14 2018-05-22 电子科技大学 A kind of passive track-corelation direction cross positioning method based on angle information
CN108089148A (en) * 2017-12-14 2018-05-29 电子科技大学 A kind of passive track-corelation direction cross positioning method based on time difference information
CN108152789A (en) * 2018-01-03 2018-06-12 电子科技大学 Utilize the passive track-corelation data correlation and localization method of RSS information
CN109001670A (en) * 2018-06-01 2018-12-14 中国人民解放军战略支援部队信息工程大学 A kind of distributed passive location method and device for combining the time difference and angle
CN108802674A (en) * 2018-07-19 2018-11-13 中国人民解放军战略支援部队信息工程大学 It is a kind of for the combined method for searching and device that directly position
CN109917333A (en) * 2019-01-31 2019-06-21 西南电子技术研究所(中国电子科技集团公司第十研究所) Merge the passive location method of AOA observed quantity and TDOA observed quantity

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIHAO YIN: "A simple and accurate TDOA- AOA localization method using two stations", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING LETTERS》 *
熊健: "一种改进的多站AOA/TDOA联合无源定位闭合形式算法", 《电讯技术》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308418A (en) * 2020-03-10 2020-06-19 慧众行知科技(北京)有限公司 Steady method for two-dimensional positioning of target with unknown height
CN111308418B (en) * 2020-03-10 2021-11-23 慧众行知科技(北京)有限公司 Steady method for two-dimensional positioning of target with unknown height
CN111337875A (en) * 2020-04-13 2020-06-26 北京航天长征飞行器研究所 Space multi-station multi-system dimension reduction positioning method
CN113390406A (en) * 2021-06-16 2021-09-14 电子科技大学 Multi-target data association and positioning method based on passive multi-sensor system
CN113390406B (en) * 2021-06-16 2022-05-24 电子科技大学 Multi-target data association and positioning method based on passive multi-sensor system
RU2790347C1 (en) * 2022-04-14 2023-02-16 Общество с ограниченной ответственностью "Специальный Технологический Центр" Difference-range method for determining the location of objects
CN115524662A (en) * 2022-10-27 2022-12-27 中国电子科技集团公司信息科学研究院 Direction finding time difference combined positioning method and system, electronic equipment and storage medium
CN115524662B (en) * 2022-10-27 2023-09-19 中国电子科技集团公司信息科学研究院 Direction finding time difference joint positioning method, system, electronic equipment and storage medium

Similar Documents

Publication Publication Date Title
CN110412504A (en) It is associated with based on angle with the passive track-corelation of time difference information and localization method
CN108089148B (en) A kind of passive track-corelation direction cross positioning method based on time difference information
CN108061877B (en) A kind of passive track-corelation direction cross positioning method based on angle information
CN110988923B (en) Satellite navigation interference source passive positioning method based on multi-platform cooperation
CN109581317A (en) One kind being based on the matched corner object localization method of echo-peak
CN106526549A (en) False target identification method with combination of two-coordinate radar and three-coordinate radar statistics
CN108152789B (en) Utilize the passive track-corelation data correlation and localization method of RSS information
CN107229033A (en) Multiple target reaching time-difference localization method based on height dimension sectioning search
CN106872942A (en) For the positioning precision calculation method of Distributed Multi positioning monitoring system
CN105188082A (en) Evaluation method for RSS (Received Signal Strength)/AOA (Angle of Arrival)/TDOA (Time Difference of Arrival) positioning performance under indoor WLAN (Wireless Local Area Network) environment
CN108614268A (en) The acoustics tracking of low altitude high speed airbound target
CN103376447A (en) Method for achieving three-dimension positioning of non-cooperative double multistatic radar target
CN104330788A (en) Radio location method based on reach path reverse tracking
Lei et al. Optimal distribution for four-station TDOA location system
CN104021292A (en) Dim target detection and tracking method based on formation active networking
CN106680783B (en) Method based on station site error's blending algorithm confrontation false targets
CN104793177A (en) Microphone array direction finding method based on least square methods
CN105787081A (en) Radiation platform correlation method based on radiation source space positions
Giacometti et al. Accuracy studies for TDOA-AOA localization of emitters with a single sensor
Kazemi et al. Efficient closed-form solution for 3-D hybrid localization in multistatic radars
CN104569915A (en) Positioning method used in multiple photoelectric detection systems and based on target movement model
CN108562893A (en) A kind of external illuminators-based radar multistation combined tracking method
CN103728608A (en) Antenna arrangement method for improving MIMO-OTH radar detecting performance in ionized layer double-Gaussian model
Ding et al. Multistatic ultra-wideband localization for NLOS environments
Wang et al. The research of indoor three-dimensional positioning algorithm based on ultra-wideband technology

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20191105

RJ01 Rejection of invention patent application after publication