CN109633724B - Passive target positioning method based on single-satellite and multi-ground-station combined measurement - Google Patents
Passive target positioning method based on single-satellite and multi-ground-station combined measurement Download PDFInfo
- Publication number
- CN109633724B CN109633724B CN201910038540.1A CN201910038540A CN109633724B CN 109633724 B CN109633724 B CN 109633724B CN 201910038540 A CN201910038540 A CN 201910038540A CN 109633724 B CN109633724 B CN 109633724B
- Authority
- CN
- China
- Prior art keywords
- satellite
- radiation source
- target
- positioning
- ground
- 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
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000005855 radiation Effects 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 2
- 230000004927 fusion Effects 0.000 claims 1
- 238000013459 approach Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/46—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being of a radio-wave signal type
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
The invention belongs to the technical field of electronic countermeasure, and particularly relates to a passive target positioning method based on single-satellite and multi-ground-station combined measurement. Aiming at a positioning scene of an aerial target radiation source, the invention respectively measures the cosine angles of the directions from which target radiation source signals reach the target radiation source signals to the respective ground stations through a single satellite and a plurality of ground stations, and the time difference and the frequency difference from which the target radiation source signals reach the satellite and a plurality of ground observation stations, and gives a weighted least square analytic solution of the target position through pseudo-linear processing of the cosine angles of the directions and fusing a time difference measurement equation and a frequency difference measurement equation. The method can realize single instantaneous high-precision positioning, has no positioning ambiguity problem, and the mean square error of the positioning error can approach the Clarmet-Roche lower limit.
Description
Technical Field
The invention belongs to the technical field of electronic countermeasure, and particularly relates to a passive target positioning method based on single-satellite and multi-ground-station combined measurement.
Background
On one hand, the key points of the modern electronic information war are the competition of electromagnetic space between enemy and my and the detection and monitoring of key targets of the opponent, so as to obtain information of strategic deployment, platform type and the like of enemy units. Therefore, the reconnaissance and positioning technology for the target plays a very important role in the field of electronic countermeasure. Compared with an active positioning method, passive positioning has the advantages of long acting distance, hidden reception, difficulty in being discovered by the opposite side and the like, and has important significance for improving the survival capability and the operational capability of the weapon system in the electronic warfare environment.
On the other hand, with the development of the aerospace technology and the breakthrough of the satellite-borne detection technology, passive detection and positioning of a ground radiation source by using satellite-borne electronic detection equipment has become one of the leading-edge issues of the space military technology which is disputed and developed in various countries at present. Although the single-satellite passive positioning system has the advantages of low cost, flexible positioning system and the like, the single-satellite passive positioning system realizes the positioning of the radiation source based on one observation satellite, so the obtained target information amount is less, and the positioning accuracy is not high enough. Therefore, the single star is considered to be combined with a plurality of ground observation stations to realize the space-ground combined target positioning.
Disclosure of Invention
The invention aims to solve the problems and provides a passive target positioning method based on single-satellite and multi-ground-station direction finding and combined time difference and frequency difference. Aiming at a positioning scene of an aerial target radiation source, a single satellite and a plurality of ground stations are used for respectively measuring the cosine angles of the directions from which target radiation source signals reach to the single satellite and the plurality of ground stations, and the time difference and the frequency difference from which the target radiation source signals reach to the satellite and the plurality of ground observation stations, and a weighted least square analytic solution of a target position is given by carrying out pseudo-linear processing on the cosine angles of the directions and fusing a time difference measurement equation and a frequency difference measurement equation. The method can realize single instantaneous high-precision positioning, has no positioning ambiguity problem, and the mean square error of the positioning error can approach the Cramer-Rao lower limit (CRLB).
The technical scheme adopted by the invention is as follows:
according to the method, cosine angles of the directions of the target radiation source signals reaching the single satellite and the multiple ground stations and time difference and frequency difference of the target radiation source signals reaching the satellite and the multiple ground observation stations are measured respectively, and a weighted least square analytic calculation formula of a target position is given through algorithm derivation. The localization model is shown in fig. 1. In the figure, the coordinate system of ground fixation is recorded in the { e system (x) e ,y e ,z e ) At x, the position of the satellite is S,e =[x s,e ,y s,e ,z s,e ] T The position of the target radiation source is x T,e =[x t,e ,y t,e ,z t,e ] T And the positions of the three ground observation stations are respectively marked as x r1,e =[x r1,e ,y r1,e ,z r1,e ] T ,x r2,e =[x r1,e ,y r1,e ,z r1,e ] T And x r3,e =[x r3,e ,y r3,e ,z r3,e ] T . In the system b of the star coordinate system, the position of the target radiation source is x T,b =[x t,b ,y t,b ,z t,b ] T . In a satellite-borne direction finding positioning system, the following coordinate transformation relationship is provided: x is a radical of a fluorine atom T,b =M(x S,e -x T,e ) Detailed description of the invention
Which comprises the following steps: satellite x S,e Corresponding subsatellite point latitude B in geodetic coordinate system S And longitude L S And by satellite x S,e The satellite-borne attitude sensor outputs several angle information: yaw angle psi, pitch angle theta, roll angle phi. Simultaneously, in the formula:
the invention mainly comprises the following steps:
a. respectively measuring the cosine angles of the directions from which the target radiation source signals reach to the single satellite and the multiple ground stations, and the time difference and the frequency difference between the target radiation source signals respectively reaching the satellite and the multiple ground observation stations;
b. establishing a target positioning model by performing pseudo-linear processing on the direction cosine angles and fusing a time difference measurement equation and a frequency difference measurement equation;
c. and (c) solving the model established in the step (b) by adopting a weighted least square method to obtain the target position.
Specifically, in the step a, the invention is based on the following principle:
the various observations of the positioning system can be expressed as follows:
A. single-star measurement target radiation source signal arrival direction cosine angle alpha s And beta s The expression of the measured value is:
in the formula: d x =[1,0,0],D y =[0,1,0],D z =[0,0,1],Respectively represent the direction cosine angle alpha s And beta s The true value of (c) is,andrespectively, representing the measurement error thereof. Meanwhile, the following steps are included:
in the formula:true value, γ, representing the cosine angle of the direction s Representing the cosine angle alpha of direction s And beta s The resulting direction cosine angle is calculated.
B. Ith ground observation station x ri,e Measuring the cosine angle alpha of the arrival direction of the target radiation source signal ri And beta ri The expression of the measured value is:
similarly, in the formula:respectively represent the direction cosine angle alpha ri And beta ri The true value of (a) is,andrespectively represent the measurement error, and simultaneously, the following are:
in the formula:representing true values of direction cosine angles, gamma ri Representing the cosine angle alpha of direction ri And beta ri The resulting direction cosine angle is calculated.
C. Converting the time difference measurement value into a distance difference, and then expressing the measurement value as follows:
in the formula:representative frequency difference ρ s,ri The true value of (a) is,representing its measurement error.
D. The frequency difference measurement can be converted to a rate of change of the range difference over time, and the expression:
in the formula:representative frequency difference xi s,ri The true value of (a) is,representing its measurement error.
Based on various measured values of the positioning system, the positioning and the speed fixing of the hollow target radiation source are realized at the same time, and the derivation of a weighted least square positioning algorithm is carried out.
In the step b, the method specifically adopted by the invention is as follows:
a) From the above formula A:
namely:
and according to the first-order Taylor series expansion principle, the following can be obtained:
the finishing is carried out to obtain:
b) From the above formula A:
the following can be obtained:
according to the first-order Taylor series expansion principle, the following can be obtained:
the carrying-in finishing can be carried out as follows:
c) From the above formula B:
the same can be obtained:
d) From the above formula B:
the following can be obtained by the same way:
e) From the above formula D:
further, according to the above definition of the cosine angle of direction a, it can be obtained:
namely:
similarly, according to the first-order Taylor series expansion principle:
the finishing can be carried out as follows:
wherein let H = M (v) S,e -v T,e ) Then, there are:
then on the basis of the above algorithm derivation:
wherein, the specific expressions of b, A, N and N are detailed as follows:
wherein u is s And u ri The following are explicit:
wherein:
wherein C is s And C ri The following are explicitly stated:
thus, a weighted least squares solution of the position and velocity vector Ω of the target radiation source can be obtained as:
since the matrix N contains the target position to be obtained, the method is firstly to solve the target by the least square solutionGiving out the target position and speed needed in the matrix N to obtain the matrix N, and further obtaining the weighted least square solution of the target radiation source
By combining the scheme, the method provided by the invention is subjected to the following error analysis:
the influence factors of the positioning accuracy of the system for positioning the hollow target radiation source mainly include the direction-finding error of the satelliteDirection finding error of ground stationTime difference measurement errorAnd frequency difference measurement errorThe Clarmet-Rao lower limit (CRLB) of the positioning error in the presence of these four systematic measurement errors is analyzed and calculated below, and the geometric distribution GDOP of the positioning error is given. The specific expression of the measurement equation of the system is as follows:
if the parameter vector is set to be omega, the Jacobian matrix of the system observation equation relative to the target position vector omegaThe method comprises the following specific steps:
the theoretical precision bound of the target positioning error without the constraint of the WGS-84 earth ellipsoid model is obtained as follows:
CRLB(Ω)=G T W -1 G
the passive positioning system has the advantages that high-precision positioning and constant speed in a wide area range are realized by jointly measuring the time difference, the frequency difference and the direction cosine angles by the single satellite and the plurality of ground stations, the conventional satellite-based passive positioning system is subjected to auxiliary optimization by combining auxiliary observation information of the ground stations, the space-ground combined passive positioning is realized, the passive positioning system has the advantages of passive positioning and space electronic satellite reconnaissance, and the precision and the real-time performance of target positioning and tracking are improved.
Drawings
FIG. 1 is a diagram of a WGS-84 model-based positioning model for two-star direction-finding positioning;
FIG. 2 is a GDOP plot of positioning error;
FIG. 3 is a simulation diagram of the variation of the position positioning resolving error with the direction finding error of the satellite;
FIG. 4 is a simulation diagram of the variation of the position positioning resolving error with the direction finding error of the ground station;
FIG. 5 is a diagram showing a simulation of the variation of the position-locating resolving error with time-difference measurement error;
FIG. 6 is a simulation diagram of the variation of velocity positioning calculation error with frequency difference measurement error.
Detailed Description
The following verification and explanation of the above positioning method with reference to the drawings first make the following reasonable assumptions on the system model:
1. assuming that the satellite is a low orbit satellite, the orbit altitude is relatively low, typically 500km to 1000km;
2. unifying satellite attitude measurement errors existing in engineering practice into satellite direction finding errors;
3. the measurement errors are assumed to follow a gaussian distribution with a mean value of zero, and the errors are independent of each other.
(1) GDOP map of positioning error:
as shown in FIG. 2, assume that the orbit height of a single star is H S =800km, and the longitude and latitude of the corresponding satellite sub-satellite point are respectively (L) S ,B S ) = (103 °,37 °), in the figure the star points are marked, the velocity vector is v S,e =[6,4,2] T km/s. The longitude and latitude of the three ground observation stations are respectively (L) r1 ,B r1 )=(100°,34°),(L r2 ,B r2 ) = (101 °,36 °) and (L) r1 ,B r1 ) = (102 °,34 °), see triangle dot notation in the figure, simulations all assuming the ground observation station to be stationary. Set the velocity vector of the target to v T,e =[0.2,0.15,0.02] T km/s. Direction finding error of satelliteThe root mean square error of (2) is set to be 0.1 degrees, and the direction-finding errors of three ground stationsThe root mean square error is set to be 0.5 DEG, and the time difference measurement error is set to beThe size of the root mean square error of (2) is set to be 10us, and the frequency difference measurement error isThe size of the root mean square error is set to be 10Hz, and a GDOP diagram for positioning and resolving the target position of single-satellite and multi-ground-station direction finding and combined time difference and frequency difference is obtained through simulation.
As can be seen from the upper graph, the GDOP graphs of the positioning errors are approximately symmetrically distributed about the single star and the intersatellite points of the three ground stations, the root mean square error of the positioning errors of the target positions near the ground stations is small, a certain deflection exists near the local intersatellite point of the single star, and the root mean square error of the positioning errors in the wide-area longitude and latitude range is stably distributed.
Assuming a single star orbit height ofH S =800km, the longitude and latitude of the satellite are (L) S ,B S ) = (103 °,37 °), in the figure the star points are marked, the velocity vector is v S,e =[6,4,2] T km/s. The longitude and latitude of three ground observation stations are respectively (L) r1 ,B r1 )=(100°,34°),(L r2 ,B r2 ) = (101 °,36 °) and (L) r1 ,B r1 ) = (102 °,34 °), see the triangular dot marks in the figure, and the ground observation station is assumed to be stationary in the simulation. Setting the longitude, latitude and elevation of the aerial target radiation source as (L) t ,B t ,H t ) = 102 °,36 °,15km, the velocity vector of the target radiation source is v T,e =[0.2,0.15,0.02] T km/s. The satellite, ground station and target positions in the following simulations are all set as described above.
(2) Influence of satellite direction finding errors:
as shown in fig. 3, the direction-finding error of the ground stationThe root mean square error is set to 0.5 DEG, and the time difference measurement errorThe root mean square error of (1) is set to be 10us, and the frequency difference measurement error isAll set as 10Hz, and the direction-finding error of the satellite The size of the root mean square error is changed from 0.1 degree to 1 degree, positioning simulation calculation is carried out, and the change of the position positioning calculation error along with the direction-finding error of the satellite is obtained.
(3) Influence of ground station direction finding error:
as shown in FIG. 4, the direction finding error of the satelliteThe root mean square error is set to be 0.3 DEG, and the time difference measurement error is set to beThe size of the root mean square error of (2) is set to be 10us, and the frequency difference measurement error isAll set as 10Hz, and the direction-finding error of the ground stationThe size of the root mean square error is changed from 0.2 degrees to 2 degrees, positioning simulation calculation is carried out, and the change of the positioning calculation error along with the direction-finding error of the ground station is obtained.
(4) Influence of time difference measurement error:
direction finding error of satelliteThe root mean square error is set to be 0.3 degrees, and the direction-finding error of the ground stationThe root mean square error of (2) is set to 0.5 DEG, and the frequency difference measurement errorAre all set to 10Hz, and the time difference measurement error n is ρ The size of the root mean square error is changed from 5us to 50us, positioning simulation calculation is carried out, and the change of the positioning calculation error along with the difference measurement error is obtained, as shown in figure 5.
(5) Influence of frequency difference measurement error:
direction finding error of satelliteThe root mean square error is set to be 0.3 degrees, and the direction-finding error of the ground stationThe root mean square error is set to 0.5 DEG, and the time difference measurement error n ρ The root mean square error of (2) is set to 10us, and the frequency difference measurement error is setThe size of the root mean square error is changed from 2Hz to 20Hz, positioning simulation calculation is carried out, and the change of the speed calculation error along with the time difference measurement error is obtained, as shown in figure 6.
As can be seen from fig. 3, 4, 5, 6, the solution proposed herein can well approximate CRLB, only slightly deviating from CRLB by about 0.1km after measurement errors are large.
Claims (1)
1. The passive target positioning method based on the single-satellite and multi-ground-station combined measurement is characterized by comprising the following steps of:
a. respectively measuring the cosine angles of the directions from which target radiation source signals arrive to the single satellite and the multiple ground stations, and the time difference and the frequency difference between the target radiation source signals respectively arrive at the satellite and the multiple ground observation stations; the specific method comprises the following steps:
setting the coordinate system of ground fixation { e system (x) e ,y e ,z e ) At x, the position of the satellite is S,e =[x s,e ,y s,e ,z s,e ] T The number of the ground observation stations is 3, and the positions are x respectively r1,e =[x r1,e ,y r1,e ,z r1,e ] T ,x r2,e =[x r1,e ,y r1,e ,z r1,e ] T And x r3,e =[x r3,e ,y r3,e ,z r3,e ] T The position of the target radiation source is x T,e =[x t,e ,y t,e ,z t,e ] T In the system b of the star coordinate system, the position of the target radiation source is x T,b =[x t,b ,y t,b ,z t,b ] T In the satellite-borne direction finding positioning system, the following coordinate transformation relationship is provided: x is the number of T,b =M(x S,e -x T,e ) WhereinIncluding satellite x S,e Corresponding subsatellite point latitude B in geodetic coordinate system S And longitude L S And by satellite x S,e The satellite-borne attitude sensor outputs several angle information: yaw angle ψ, pitch angle θ, roll angle Φ:
a1, measuring the cosine angle alpha of the arrival direction of the target radiation source signal by a single star s And beta s :
In the formula: d x =[1,0,0],D y =[0,1,0],D z =[0,0,1],Respectively representing the cosine angle alpha of the direction s And beta s The true value of (c) is,andrespectively representing the measurement errors;
order to
In the formula:true value, γ, representing the cosine angle of the direction s Is represented by a direction cosine angle α s And beta s Calculating the direction cosine angle obtained by the measured value;
a2, i-th ground observation station x ri,e Measuring the cosine angle alpha of the arrival direction of the target radiation source signal ri And beta ri :
In the formula:respectively represent the direction cosine angle alpha ri And beta ri The true value of (c) is,andrespectively representing the measurement errors;
order to
In the formula:representing true values of direction cosine angles, gamma ri Is represented by a direction cosine angle α ri And beta ri Calculating the direction cosine angle obtained by the measured value;
a3, obtaining time differences of target radiation source signals respectively reaching a satellite and a ground observation station, and converting time difference measurement values into distance differences:
in the formula:representative time difference p s,ri The true value of (a) is,representing the measurement error thereof;
a4, obtaining frequency differences of target radiation source signals respectively reaching the satellite and the ground observation station, and converting frequency difference measurement values into the change rate of distance differences along with time:
in the formula:representative frequency difference xi s,ri The true value of (a) is,representing the measurement error thereof;
b. establishing a target positioning model by performing pseudo-linear processing on the direction cosine angles and fusing a time difference measurement equation and a frequency difference measurement equation; the specific method comprises the following steps:
b1, performing pseudo linearization on the direction cosine angle obtained in the step a1 to obtain:
b2, performing pseudo linearization on the direction cosine angle obtained in the step a2 to obtain:
b3, fusing the time difference measurement equation and the frequency difference measurement equation, and establishing a positioning model, which specifically comprises the following steps:
from the above pseudo-linearization equation
Wherein, let H = M (v) S,e -v T,e ) Then, there are:
then the fusion yields:
wherein,
c. b, solving the model established in the step b by adopting a weighted least square method to obtain a target position; the specific method comprises the following steps:
and solving by adopting a weighted least square method to obtain a weighted least square solution of the target radiation source as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910038540.1A CN109633724B (en) | 2019-01-16 | 2019-01-16 | Passive target positioning method based on single-satellite and multi-ground-station combined measurement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910038540.1A CN109633724B (en) | 2019-01-16 | 2019-01-16 | Passive target positioning method based on single-satellite and multi-ground-station combined measurement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633724A CN109633724A (en) | 2019-04-16 |
CN109633724B true CN109633724B (en) | 2023-03-03 |
Family
ID=66062126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910038540.1A Active CN109633724B (en) | 2019-01-16 | 2019-01-16 | Passive target positioning method based on single-satellite and multi-ground-station combined measurement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633724B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110988851B (en) * | 2019-12-19 | 2023-11-21 | 南京邮电大学 | Different-orbit single-star time-sharing frequency measurement positioning method based on star optimization |
CN111736187B (en) * | 2020-06-22 | 2022-05-24 | 北京理工大学 | High-precision high-sensitivity single-satellite GNSS positioning method based on passive synthetic aperture |
CN112114341B (en) * | 2020-08-13 | 2022-01-25 | 中国人民解放军军事科学院国防科技创新研究院 | Low-earth-orbit satellite cooperative frequency measurement passive positioning method |
CN112179355B (en) * | 2020-09-02 | 2023-05-26 | 西安交通大学 | Attitude estimation method aiming at typical characteristics of luminosity curve |
CN112285645B (en) * | 2020-12-25 | 2021-03-26 | 中国人民解放军国防科技大学 | Positioning method, storage medium and system based on high-orbit satellite observation platform |
CN112904325A (en) * | 2021-01-20 | 2021-06-04 | 西安电子科技大学 | Double-star offshore target positioning method based on microwave forwarding |
CN114518585B (en) * | 2022-01-28 | 2022-10-21 | 深圳市北斗天宇科技有限公司 | Target positioning method, device, terminal and storage medium based on single satellite |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5099245A (en) * | 1987-10-23 | 1992-03-24 | Hughes Aircraft Company | Vehicle location system accuracy enhancement for airborne vehicles |
US7436359B1 (en) * | 2007-05-30 | 2008-10-14 | Northrop Grumman Systems Corporation | Method for single satellite geolocation of emitters using an ambiguous interferometer array |
CN103969622A (en) * | 2014-04-25 | 2014-08-06 | 西安电子科技大学 | Time difference positioning method based on multiple motion receiving stations |
CN104535964A (en) * | 2014-12-24 | 2015-04-22 | 大连大学 | Helmet type microphone array sound source positioning method based on low-frequency diffraction delay inequalities |
CN104880723A (en) * | 2015-04-29 | 2015-09-02 | 中国电子科技集团公司第三十六研究所 | Method and device for positioning of ground radiation sources through single-satellite direction finding |
WO2016202748A1 (en) * | 2015-06-17 | 2016-12-22 | Thales | Method and device for locating an electromagnetic emission source and system implementing such a method |
CN106842118A (en) * | 2016-12-23 | 2017-06-13 | 安徽四创电子股份有限公司 | For the time difference positioning method of Distributed Multi positioning monitoring system |
CN106908819A (en) * | 2017-03-14 | 2017-06-30 | 西安电子科技大学 | Height rail double star time-varying high receives the when frequency difference estimation method of signal |
CN106950533A (en) * | 2016-12-05 | 2017-07-14 | 南京理工大学 | A kind of sane beam forming combined positioning-method of self adaptation for suppressing interference |
CN107015199A (en) * | 2017-05-09 | 2017-08-04 | 南京航空航天大学 | A kind of double unmanned plane direction finding time difference positioning methods for considering UAV Attitude angle |
CN108226978A (en) * | 2018-01-15 | 2018-06-29 | 电子科技大学 | A kind of Double-Star Positioning System method based on WGS-84 models |
CN108445446A (en) * | 2018-02-06 | 2018-08-24 | 中国人民解放军战略支援部队信息工程大学 | A kind of passive method for locating speed measurement and device |
CN108562295A (en) * | 2018-07-23 | 2018-09-21 | 电子科技大学 | A kind of three station time difference orbit determination methods based on two body Model of synchronous satellite |
CN108761387A (en) * | 2018-06-01 | 2018-11-06 | 西安邮电大学 | A kind of dual station TDOA-FDOA joint location method of Fixed Emitter |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8188919B2 (en) * | 2009-06-16 | 2012-05-29 | L3 Communications Integrated Systems, L.P. | Globally-convergent geo-location algorithm |
TWI498581B (en) * | 2014-01-03 | 2015-09-01 | Ind Tech Res Inst | Satellite positioning method, satellite positioning apparatus, and computer-readble medium |
-
2019
- 2019-01-16 CN CN201910038540.1A patent/CN109633724B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5099245A (en) * | 1987-10-23 | 1992-03-24 | Hughes Aircraft Company | Vehicle location system accuracy enhancement for airborne vehicles |
US7436359B1 (en) * | 2007-05-30 | 2008-10-14 | Northrop Grumman Systems Corporation | Method for single satellite geolocation of emitters using an ambiguous interferometer array |
CN103969622A (en) * | 2014-04-25 | 2014-08-06 | 西安电子科技大学 | Time difference positioning method based on multiple motion receiving stations |
CN104535964A (en) * | 2014-12-24 | 2015-04-22 | 大连大学 | Helmet type microphone array sound source positioning method based on low-frequency diffraction delay inequalities |
CN104880723A (en) * | 2015-04-29 | 2015-09-02 | 中国电子科技集团公司第三十六研究所 | Method and device for positioning of ground radiation sources through single-satellite direction finding |
WO2016202748A1 (en) * | 2015-06-17 | 2016-12-22 | Thales | Method and device for locating an electromagnetic emission source and system implementing such a method |
CN106950533A (en) * | 2016-12-05 | 2017-07-14 | 南京理工大学 | A kind of sane beam forming combined positioning-method of self adaptation for suppressing interference |
CN106842118A (en) * | 2016-12-23 | 2017-06-13 | 安徽四创电子股份有限公司 | For the time difference positioning method of Distributed Multi positioning monitoring system |
CN106908819A (en) * | 2017-03-14 | 2017-06-30 | 西安电子科技大学 | Height rail double star time-varying high receives the when frequency difference estimation method of signal |
CN107015199A (en) * | 2017-05-09 | 2017-08-04 | 南京航空航天大学 | A kind of double unmanned plane direction finding time difference positioning methods for considering UAV Attitude angle |
CN108226978A (en) * | 2018-01-15 | 2018-06-29 | 电子科技大学 | A kind of Double-Star Positioning System method based on WGS-84 models |
CN108445446A (en) * | 2018-02-06 | 2018-08-24 | 中国人民解放军战略支援部队信息工程大学 | A kind of passive method for locating speed measurement and device |
CN108761387A (en) * | 2018-06-01 | 2018-11-06 | 西安邮电大学 | A kind of dual station TDOA-FDOA joint location method of Fixed Emitter |
CN108562295A (en) * | 2018-07-23 | 2018-09-21 | 电子科技大学 | A kind of three station time difference orbit determination methods based on two body Model of synchronous satellite |
Non-Patent Citations (4)
Title |
---|
"A Novel Method for Determining DOA from Far-Field TDOA or FDOA";Karleigh J. Cameron等;《Electrical Engineering and Systems Science》;20181231;1-7 * |
"WGS-84模型下时差频差半定规划定位算法";李万春,等;《航空学报》;20170731;第38卷(第7期);288-295 * |
"利用TDOA和FDOA的单站多外辐射源目标定位算法";赵勇胜,等;《四川大学学报(工程科学版)》;20161231;第48卷(第S1期);170-177 * |
"单星测频无源定位法";陆安南,等;《通信学报》;20040930(第09期);160-168 * |
Also Published As
Publication number | Publication date |
---|---|
CN109633724A (en) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633724B (en) | Passive target positioning method based on single-satellite and multi-ground-station combined measurement | |
CN109613583B (en) | Passive target positioning method based on single star and ground station direction finding and combined time difference | |
CN101270993B (en) | Remote high-precision independent combined navigation locating method | |
CN104165640B (en) | Near-space missile-borne strap-down inertial navigation system transfer alignment method based on star sensor | |
CN103674034B (en) | Multi-beam test the speed range finding revise robust navigation method | |
CN101339244B (en) | On-board SAR image automatic target positioning method | |
CN102998690B (en) | Attitude angle direct resolving method based on global position system (GPS) carrier wave double-difference equation | |
CN105242285B (en) | A kind of Navigation of Pilotless Aircraft data based on satellite communication are spoofed recognition methods | |
CN103900576B (en) | A kind of information fusion method of survey of deep space independent navigation | |
CN103744052B (en) | A kind of double star being applied to aerial target location surveys time-of-arrival direction finding method and device | |
CN107132542B (en) | A kind of small feature loss soft landing autonomic air navigation aid based on optics and Doppler radar | |
CN111380518B (en) | SINS/USBL tight combination navigation positioning method introducing radial velocity | |
CN108917764B (en) | Distance-only relative navigation method for double-star formation | |
CN102323571B (en) | Distribution method of satellite-borne dual-antenna SAR (Synthetic Aperture Radar) interferometric calibrator with comprehensive overall parameter | |
CN110081905B (en) | Light wave atmospheric refraction error calculation method based on single-station photoelectric theodolite | |
CN103674032A (en) | Satellite autonomous navigation system and method integrating pulsar radiation vector and timing observation | |
CN102116628A (en) | High-precision navigation method for landed or attached deep sky celestial body detector | |
CN104049269A (en) | Target navigation mapping method based on laser ranging and MEMS/GPS integrated navigation system | |
Liu et al. | Accurate mapping method for UAV photogrammetry without ground control points in the map projection frame | |
CN111505608B (en) | Laser pointing on-orbit calibration method based on satellite-borne laser single-chip footprint image | |
CN117192578A (en) | Shipborne measurement and control antenna shafting parameter calibration method for tracking unmanned aerial vehicle | |
CN116105730A (en) | Angle measurement-only optical combination navigation method based on cooperative target satellite very short arc observation | |
CN102607563B (en) | System for performing relative navigation on spacecraft based on background astronomical information | |
CN114063054A (en) | Method for backtracking and calculating rocket launching point parameters | |
RU2723199C1 (en) | Method and system for determining orientation of spacecraft in space with autonomous correction of light aberration effect |
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 |