CN109613583B - Passive target positioning method based on single star and ground station direction finding and combined time difference - Google Patents

Passive target positioning method based on single star and ground station direction finding and combined time difference Download PDF

Info

Publication number
CN109613583B
CN109613583B CN201910002322.2A CN201910002322A CN109613583B CN 109613583 B CN109613583 B CN 109613583B CN 201910002322 A CN201910002322 A CN 201910002322A CN 109613583 B CN109613583 B CN 109613583B
Authority
CN
China
Prior art keywords
satellite
radiation source
target
positioning
time difference
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
CN201910002322.2A
Other languages
Chinese (zh)
Other versions
CN109613583A (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.)
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 CN201910002322.2A priority Critical patent/CN109613583B/en
Publication of CN109613583A publication Critical patent/CN109613583A/en
Application granted granted Critical
Publication of CN109613583B publication Critical patent/CN109613583B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/46Determining 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention belongs to the technical field of electronic countermeasure, and particularly relates to a passive target positioning method based on single star and ground station direction finding and combined time difference measuring. According to the method, aiming at a positioning scene of an aerial target radiation source, the single star and the ground station are used for respectively measuring the directions of the target radiation source signals reaching the respective residual angles, and the time difference of the target radiation source signals reaching the satellite and the ground observation station is used for giving a weighted least square analysis solution of the target position through pseudo-linearization processing of the directions of the residual angles and fusion of a time difference measurement equation. The method can realize single instantaneous high-precision positioning, does not have the problem of positioning and de-blurring, and the mean square error of positioning errors can approach the lower limit of the Keramen-Row.

Description

Passive target positioning method based on single star and ground station direction finding and combined time difference
Technical Field
The invention belongs to the technical field of electronic countermeasure, and particularly relates to a passive target positioning method based on single star and ground station direction finding and combined time difference measuring.
Background
On the one hand, in a modern electronic information battle system, by utilizing an electronic reconnaissance means, electromagnetic wave signals from enemy battle units are measured, analyzed, sorted and identified to obtain strategic deployment of the enemy units, platform type and other information, the electronic reconnaissance system has important reference values for strategic decision making of the higher layers of the enemy, protection of the enemy battle units and strategic striking of the enemy. The passive reconnaissance positioning technology has the advantages of good concealment, long acting distance, high anti-interference capability, high system safety protection and the like, occupies an increasingly important position in an electronic reconnaissance system, and has been valued by the military of various countries.
On the other hand, along with development of the aerospace technology and breakthrough of the satellite-borne detection technology, the passive reconnaissance positioning is more and more dependent on the electronic reconnaissance satellite, and the single-satellite passive positioning system only utilizes one observation satellite to realize positioning of the ground radiation source, so that the single-satellite passive positioning system has the advantages of low cost, flexible positioning system and the like, but simultaneously, the single-satellite passive positioning system only utilizes one observation satellite to realize positioning of the radiation source, the target information amount obtained by a single platform is less, the positioning precision is not high enough, and under the normal condition, the space position and the motion parameter of the target radiation source can be estimated only by continuously observing the ground target radiation source for a plurality of times. At present, a single-star passive positioning method based on a kinematics principle obtains relatively high positioning precision through information such as the change rate of the arrival angle of a target radiation source signal, the change rate of Doppler frequency and the like, but the method needs multiple observations and cannot realize instantaneous single-time high-precision positioning.
Disclosure of Invention
The invention aims to solve the problems and provides a passive target positioning method based on single-satellite and ground station direction finding and combined time difference. Aiming at a positioning scene of an aerial target radiation source, measuring the directions and the residual angles of the target radiation source signals reaching the respective directions respectively through a single satellite and a ground station, and obtaining the time difference of the target radiation source signals reaching the satellite and the ground observation station respectively, and obtaining a weighted least square analysis solution of the target position through pseudo-linearization processing of the directions and the residual angles and fusing a time difference measurement equation. The method can realize single instantaneous high-precision positioning, does not have the problem of positioning and de-blurring, and the mean square error of positioning errors can approach the Cramer-Row lower limit (CRLB).
And (3) giving a weighted least square analytical solution of the target radiation source under the condition constraint by pseudo-linearization of a direction finding equation and fusion of WGS-84 ellipsoid constraint. The method can uniformly realize the positioning solution of the target when the two satellites are in the same orbit and different orbits, the satellite configuration is flexible under the positioning method, the iterative solution is not needed, the calculated amount is small, and the problem of positioning solution ambiguity does not exist.
The technical scheme adopted by the invention is as follows:
according to the method, the single star and the ground station are used for respectively measuring the complementary angles of directions of target radiation source signals reaching the single star and the ground station, and the time difference of the target radiation source signals respectively reaching the satellite and the ground observation station is deduced through an algorithm, so that a weighted least square analytic calculation formula of the target position is given. The positioning model is shown in fig. 1. In the figure, the satellite position is x under the ECEF of the ground fixed coordinate system S,e =[x s,e ,y s,e ,z s,e ] T The position of the ground observation station is x r,e =[x r,e ,y r,e ,z r,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 star coordinate system b, the position of the target radiation source is x T,b =[x t,b ,y t,b ,z t,b ] T From document [1]As can be seen, in the satellite-borne direction-finding positioning system, there is the following coordinate conversion relationship: x is x T,b =M(x S,e -x T,e ) The specific expression is
Figure GDA0003987147420000021
The method comprises the following steps: satellite x S,e Corresponding latitude B of point under satellite in geodetic coordinate system S And longitude L S By satellite x S,e Several angle information output by the satellite-borne attitude sensor: yaw angle ψ, pitch angle θ, roll angle φ. And simultaneously, in the formula:
Figure GDA0003987147420000022
Figure GDA0003987147420000023
the invention mainly comprises the following steps:
a. measuring the residual angles of directions of the target radiation source signals reaching the satellite and the ground observation station respectively through a single satellite and a ground station, and the time difference of the target radiation source signals reaching the satellite and the ground observation station respectively;
b. establishing a target positioning model by pseudo-linearization processing of the directional residual angle and fusing a time 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 present invention is based on the following principle:
the observations of the positioning system can be expressed as follows:
A. single star measuring target radiation source signal arrival direction residual angle alpha 1 and β1 The expression of the measured value is:
Figure GDA0003987147420000031
in the formula :Dx =[1,0,0],D y =[0,1,0],D z =[0,0,1],
Figure GDA0003987147420000032
Respectively represent the direction residual angle alpha 1 and β1 Is true of>
Figure GDA0003987147420000033
and />
Figure GDA0003987147420000034
Representing their measurement errors, respectively. Meanwhile, there are:
Figure GDA0003987147420000035
in the formula :
Figure GDA0003987147420000036
true value, gamma, representing the angle of the directional residual 1 Representing the angle alpha of the tailline from direction 1 and β1 The calculated direction residual angle is calculated from the measured value of (2).
B. Ground observation station measures the direction residual angle alpha of arrival of target radiation source signal 2 and β2 The expression of the measured value is:
Figure GDA0003987147420000037
similarly, in the formula:
Figure GDA0003987147420000038
respectively represent the direction residual angle alpha 2 and β2 Is true of>
Figure GDA0003987147420000039
and />
Figure GDA00039871474200000310
Representing the measurement errors respectively, and simultaneously:
Figure GDA00039871474200000311
in the formula :
Figure GDA00039871474200000312
true value, gamma, representing the angle of the directional residual 2 Represented by azimuth angle alpha 2 And pitch angle beta 2 The calculated direction residual angle is calculated from the measured value of (2).
C. The time difference measurement value is converted into a distance difference, and the expression of the measurement value is as follows:
ρ=ρ 0 +n ρ
=d t,s -d t,r +n ρ
=||(x T,e -x S,e )||-||(x T,e -x r,e )||+n ρ
in the formula :nρ Is the distance difference measurement error converted from the time difference measurement error.
The positioning of the target radiation source is converted into the problem of solving the position of the target radiation source by utilizing the plurality of groups of nonlinear equations, and from the basic principle of geometric positioning, the weighted least square solution of the target can be determined by measuring a spatial positioning line formed by direction finding, a spatial positioning hyperboloid formed by time difference measurement, intersecting the line and intersecting the line and the plane.
In the step b, the method specifically comprises the following steps:
a) From the above formula A:
Figure GDA0003987147420000041
namely:
Figure GDA0003987147420000042
and according to the first-order Taylor series expansion principle, the method can be as follows:
Figure GDA0003987147420000043
the carrying-in arrangement can be obtained:
Figure GDA0003987147420000044
b) From the above formula A:
Figure GDA0003987147420000045
/>
the method can obtain:
Figure GDA0003987147420000046
and according to the first-order Taylor series expansion principle, the method can be as follows:
Figure GDA0003987147420000047
the carrying-in arrangement can be obtained:
Figure GDA0003987147420000051
c) From the above formula B:
Figure GDA0003987147420000052
the same principle can be obtained:
Figure GDA0003987147420000053
d) From the above formula B:
Figure GDA0003987147420000054
the same principle can be obtained:
Figure GDA0003987147420000055
e) From the above formula A, B:
Figure GDA0003987147420000056
in the formula ,
Figure GDA0003987147420000057
the arrangement can be:
Figure GDA0003987147420000058
in the formula ui ,C i ,n i (i=1, 2) is represented by the following expression:
Figure GDA0003987147420000061
then, in summary, it is available:
b-Ax T,e =Nn
wherein ,
Figure GDA0003987147420000062
Figure GDA0003987147420000063
Figure GDA0003987147420000064
/>
Figure GDA0003987147420000071
Figure GDA0003987147420000072
thus, a weighted least squares solution for the target radiation source can be obtained as:
Figure GDA0003987147420000073
wherein q=cov (Nn) =nwn T And (3) recording:
Figure GDA0003987147420000074
is the covariance matrix of the noise. Since the matrix N contains the position of the object to be solved, it is done here by first solving for +.>
Figure GDA0003987147420000075
Giving the required target position in the matrix N to obtain the matrix N, and further obtaining the weighted least square solution of the target radiation source>
Figure GDA0003987147420000076
In combination with the scheme, the method provided by the invention is used for carrying out the following error analysis:
positioning the air target radiation source, wherein the influence factors of the system positioning accuracy mainly comprise the direction finding error of the satellite
Figure GDA0003987147420000077
Direction finding error of ground station>
Figure GDA0003987147420000078
Time difference measurement error n ρ Analytical calculations are performed below on the lower keramen-lo limit (CRLB) of the positioning error in the presence of these three measurement errors, and a Geometric Distribution (GDOP) diagram of the positioning error is thus given. The specific expression of the measurement equation of the system is as follows: />
Figure GDA0003987147420000081
Let the parameter vector be x T,e Then the system observation equation is related to the target position vector x T,e Jacobian matrix of (a)
Figure GDA0003987147420000082
The method comprises the following steps:
Figure GDA0003987147420000083
the covariance matrix of the measured noise is recorded as:
Figure GDA0003987147420000084
the theoretical accuracy bound for the target positioning error is:
CRLB(x T,e )=G T W -1 G
the invention has the beneficial effects that the auxiliary observation information of the ground station is introduced to carry out auxiliary optimization on the existing positioning system based on the satellite passive positioning, so that the combined passive positioning of the earth and the sky is realized, the advantages of passive positioning and space electronic satellite reconnaissance are realized, and the accuracy and the instantaneity of target positioning and tracking are improved.
Drawings
FIG. 1 is a positioning model diagram of a two-star direction finding positioning based on a WGS-84 model;
FIG. 2 is a GDOP map of positioning errors;
FIG. 3 is a simulation plot of the position solution error as a function of the direction finding error of the satellite;
FIG. 4 is a simulation plot of the position solution error as a function of ground station direction error;
fig. 5 is a simulation diagram of the variation of the positioning solution error with time difference measurement error.
Detailed Description
The above positioning method is verified and described below with reference to the accompanying drawings, and the following reasonable assumption is made 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 zero mean and are independent of each other.
(1) GDOP map of positioning error:
as shown in FIG. 2, assume that the orbital height of a single star is H S =800 km, the corresponding satellite points below the satellite have a longitude and latitude of (L S ,B S ) = (103 °,37 °), see star point labels in the figure. The longitude and latitude of the ground observation station are respectively (L) r ,B r ) = (100 °,34 °), see triangle point labels in the figure. Error of direction finding of satellite
Figure GDA0003987147420000091
The root mean square error of (2) is set to 0.1 DEG, the direction finding error of the ground station is +.>
Figure GDA0003987147420000092
The root mean square error of (2) is set to 0.5 DEG, and the time difference measurement error n ρ The root mean square error of (2) is set to be 10us, and the GDOP map for the positioning calculation of the direction finding and the combined time difference of the single star and the ground station is obtained through simulation.
From the above figures, it can be seen that the GDOP map of positioning errors is approximately symmetrically distributed about the undersea point of a single star, has certain deformation and deflection locally, and is stably distributed in the wide-area longitude and latitude range.
Setting the orbit height of a single star as H S =800 km, the corresponding satellite points below the satellite have a longitude and latitude of (L S ,B S ) = (103 °,37 °), the longitude and latitude of the ground observation station are (L r ,B r ) = (100 °,34 °), the longitude and latitude of the earth surface target radiation source are (L t ,B t ) = (102 °,36 °) when the target radiation source is about 300km from the ground observation station. In the following simulations the satellite, ground station and target positions were set as described above.
(2) Influence of satellite direction finding errors:
as shown in fig. 3, the direction finding error of the ground station
Figure GDA0003987147420000101
The root mean square error of (2) is set to 1 DEG, and the time difference measurement error n ρ The root mean square error of (2) is set to 10us, and the direction finding error of the satellite is set to +.>
Figure GDA0003987147420000102
The root mean square error of the satellite is changed from 0.1 degree to 2.5 degrees, and positioning simulation solution is carried out to obtain the change of the positioning solution error along with the direction finding error of the satellite.
(3) Effect of ground station direction finding error:
error of direction finding of satellite
Figure GDA0003987147420000103
The root mean square error of (2) is set to 0.3 DEG, and the time difference measurement error n ρ The root mean square error of (2) is set to 10us, and the direction finding error of the ground station is +.>
Figure GDA0003987147420000104
The root mean square error of the ground station is changed from 0.1 degrees to 3.1 degrees, and positioning simulation solution is carried out to obtain the change of the positioning solution error along with the direction finding error of the ground station.
(4) Effect of time difference measurement errors:
error of direction finding of satellite
Figure GDA0003987147420000105
The root mean square error of (2) is set to 0.3 DEG, and the direction finding error of the ground station
Figure GDA0003987147420000106
The root mean square error of (2) is set to 1 DEG, and the time difference measurement error n ρ The root mean square error of (1) is changed from 1us to 31us, and positioning simulation solution is carried out to obtain the change of the positioning solution error with time difference measurement error.
As can be seen from fig. 3, 4 and 5, the solution proposed herein approximates the CRLB well, and deviates slightly from the CRLB by about 0.1km only after the measurement error is large.

Claims (1)

1. The passive target positioning method based on single star and ground station direction finding and joint time difference is characterized by comprising the following steps:
a. measuring the residual angles of directions of the target radiation source signals reaching the satellite and the ground observation station respectively through a single satellite and a ground station, and the time difference of the target radiation source signals reaching the satellite and the ground observation station respectively; the specific method comprises the following steps:
the satellite position is x under the ECEF fixed coordinate system S,e =[x s,e ,y s,e ,z s,e ] T The position of the ground observation station is x r,e =[x r,e ,y r,e ,z r,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 star coordinate system b, 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 conversion relation exists: x is x T,b =M(x S,e -x T,e), wherein
Figure FDA0003987147410000011
Including satellite x S,e Corresponding latitude B of point under satellite in geodetic coordinate system S And longitude L S By satellite x S,e Is based on satellite-borne attitude sensingSeveral angle information output by the device: yaw angle ψ, pitch angle θ, roll angle φ:
Figure FDA0003987147410000012
Figure FDA0003987147410000013
a1, measuring the direction residual angle alpha of arrival of the target radiation source signal by a single star 1 and β1
Figure FDA0003987147410000014
in the formula :Dx =[1,0,0],D y =[0,1,0],D z =[0,0,1],
Figure FDA0003987147410000015
Respectively represent the direction residual angle alpha 1 and β1 Is true of>
Figure FDA0003987147410000016
and />
Figure FDA0003987147410000017
Representing the measurement errors thereof respectively;
order the
Figure FDA0003987147410000018
in the formula :
Figure FDA0003987147410000019
true value, gamma, representing the angle of the directional residual 1 Representing the angle alpha of the tailline from direction 1 and β1 Measuring meter of (2)Calculating the obtained direction residual chord angle;
a2, measuring the direction residual angle alpha of arrival of the target radiation source signal by the observation station 2 and β2
Figure FDA0003987147410000021
in the formula :
Figure FDA0003987147410000022
respectively represent the direction residual angle alpha 2 and β2 Is true of>
Figure FDA0003987147410000023
and />
Figure FDA0003987147410000024
Representing the measurement errors thereof respectively; />
Order the
Figure FDA0003987147410000025
in the formula :
Figure FDA0003987147410000026
true value, gamma, representing the angle of the directional residual 2 Represented by azimuth angle alpha 2 And pitch angle beta 2 A direction residual angle obtained by calculation of the measured value of (2);
a3, obtaining time differences of the target radiation source signals respectively reaching the satellite and the ground observation station, and converting the time difference measured value into a distance difference:
ρ=||(x T,e -x S,e )||-||(x T,e -x r,e )||+n ρ
in the formula :nρ A distance difference measurement error converted from the time difference measurement error;
b. establishing a target positioning model by pseudo-linearization processing of the directional residual angle and fusing a time difference measurement equation; the specific method comprises the following steps:
b1, pseudo-linearizing the direction residual angle obtained in the step a1 to obtain:
Figure FDA0003987147410000027
b2, pseudo-linearizing the direction residual angle obtained in the step a2 to obtain:
x r,e cosβ 2 -y r,e cosα 2 -(cosβ 2 D x -cosα 2 D y )x T,e =sinα 2 (y r,e -y t,e )n α2 +sinβ 2 (x t,e -x r,e )n β2
Figure FDA0003987147410000031
b3, fusing a time difference measurement equation, and establishing a positioning model, wherein the method specifically comprises the following steps of:
from the pseudo-linearization equation described above
Figure FDA0003987147410000032
in the formula ui ,C i ,n i See the following expression:
Figure FDA0003987147410000033
fusion results in:
b-Ax T,e =Nn
wherein ,
Figure FDA0003987147410000034
Figure FDA0003987147410000041
Figure FDA0003987147410000042
/>
Figure FDA0003987147410000043
Figure FDA0003987147410000044
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:
solving by adopting a weighted least square method to obtain a weighted least square solution of the target radiation source, wherein the weighted least square solution is as follows:
Figure FDA0003987147410000045
in which q=cov (N) =nwn T And (3) recording:
Figure FDA0003987147410000046
is the covariance matrix of the noise; the matrix N contains the position of the target to be solved, and the least square solution of the target is firstly adopted
Figure FDA0003987147410000047
Giving the required target position in the matrix N to obtain the matrix N, and further obtaining the weighted least square solution of the target radiation source>
Figure FDA0003987147410000051
/>
CN201910002322.2A 2019-01-02 2019-01-02 Passive target positioning method based on single star and ground station direction finding and combined time difference Active CN109613583B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910002322.2A CN109613583B (en) 2019-01-02 2019-01-02 Passive target positioning method based on single star and ground station direction finding and combined time difference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910002322.2A CN109613583B (en) 2019-01-02 2019-01-02 Passive target positioning method based on single star and ground station direction finding and combined time difference

Publications (2)

Publication Number Publication Date
CN109613583A CN109613583A (en) 2019-04-12
CN109613583B true CN109613583B (en) 2023-05-26

Family

ID=66016077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910002322.2A Active CN109613583B (en) 2019-01-02 2019-01-02 Passive target positioning method based on single star and ground station direction finding and combined time difference

Country Status (1)

Country Link
CN (1) CN109613583B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991571B (en) * 2019-04-17 2022-11-18 电子科技大学 Long-distance distributed passive direction finding positioning array arrangement method
CN110061801B (en) * 2019-04-25 2020-11-17 电子科技大学 Abnormal target discrimination system and method in aviation monitoring system
CN110988851B (en) * 2019-12-19 2023-11-21 南京邮电大学 Different-orbit single-star time-sharing frequency measurement positioning method based on star optimization
CN111273223A (en) * 2020-01-20 2020-06-12 西安交通大学 Passive positioning countermeasure method based on double-antenna time delay
CN111308530B (en) * 2020-02-17 2021-12-03 中国人民解放军战略支援部队信息工程大学 Short wave multi-station and single-satellite cooperative direct positioning method based on two-dimensional direction of arrival
CN112051543B (en) * 2020-08-30 2023-08-08 西南电子技术研究所(中国电子科技集团公司第十研究所) Passive positioning system for ground terminal of low-orbit satellite mobile communication
CN112462328B (en) * 2020-11-17 2024-01-26 安徽理工大学 Weighted orientation method based on sensor array arrival time difference measurement subset
CN113189625A (en) * 2021-04-29 2021-07-30 中国电子科技集团公司第三十六研究所 Error correction method and device based on single-star interferometer direction finding system and satellite
CN113642252B (en) * 2021-08-31 2022-05-17 中国人民解放军国防科技大学 Target positioning method and device based on single satellite, electronic equipment and medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772228A (en) * 2016-11-23 2017-05-31 山西奥克斯电子系统工程中心 Aerial target radiation source localization method based on arriving signal intensity
CN106886034A (en) * 2017-01-12 2017-06-23 清华大学 A kind of relay type based on single user multiple spot cheats the localization method in source
CN108226978A (en) * 2018-01-15 2018-06-29 电子科技大学 A kind of Double-Star Positioning System method based on WGS-84 models

Family Cites Families (1)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772228A (en) * 2016-11-23 2017-05-31 山西奥克斯电子系统工程中心 Aerial target radiation source localization method based on arriving signal intensity
CN106886034A (en) * 2017-01-12 2017-06-23 清华大学 A kind of relay type based on single user multiple spot cheats the localization method in source
CN108226978A (en) * 2018-01-15 2018-06-29 电子科技大学 A kind of Double-Star Positioning System method based on WGS-84 models

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
单星测频无源定位法;陆安南,等;《通信学报》;20040930(第09期);160-168 *

Also Published As

Publication number Publication date
CN109613583A (en) 2019-04-12

Similar Documents

Publication Publication Date Title
CN109613583B (en) Passive target positioning method based on single star and ground station direction finding and combined time difference
CN109633724B (en) Passive target positioning method based on single-satellite and multi-ground-station combined measurement
CN101270993B (en) Remote high-precision independent combined navigation locating method
CN101893440B (en) Celestial autonomous navigation method based on star sensors
CN104165640B (en) Near-space missile-borne strap-down inertial navigation system transfer alignment method based on star sensor
CN102998690B (en) Attitude angle direct resolving method based on global position system (GPS) carrier wave double-difference equation
CN107132542B (en) A kind of small feature loss soft landing autonomic air navigation aid based on optics and Doppler radar
CN108226978B (en) Double-satellite positioning method based on WGS-84 model
CN105548976A (en) Shipborne radar offshore precision identification method
CN111380518B (en) SINS/USBL tight combination navigation positioning method introducing radial velocity
CN102323571B (en) Distribution method of satellite-borne dual-antenna SAR (Synthetic Aperture Radar) interferometric calibrator with comprehensive overall parameter
CN103674032A (en) Satellite autonomous navigation system and method integrating pulsar radiation vector and timing observation
CN103913180A (en) Mounting angle calibration method for onboard large-view-field high-precision star sensor
CN105184002A (en) Simulation analysis method for pointing angle of data transmission antenna
CN103900576A (en) Information fusion method for autonomous navigation of deep space detection
CN111505608B (en) Laser pointing on-orbit calibration method based on satellite-borne laser single-chip footprint image
CN107300700B (en) Agile synthetic aperture radar satellite bunching mode attitude maneuver demand calculation method
CN103322969B (en) A kind of boat-carrying radio survey device antenna deformation measurement method
CN113218577A (en) Outfield measurement method for star point centroid position precision of star sensor
CN106996779B (en) Ultraviolet sensors systematic error on-orbit calibration method based on GNSS
CN102607563B (en) System for performing relative navigation on spacecraft based on background astronomical information
RU2723199C1 (en) Method and system for determining orientation of spacecraft in space with autonomous correction of light aberration effect
Giorgi Attitude determination
CN116105730A (en) Angle measurement-only optical combination navigation method based on cooperative target satellite very short arc observation
Wang et al. Single-satellite positioning algorithm based on direction-finding

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