CN111947668A - Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector - Google Patents

Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector Download PDF

Info

Publication number
CN111947668A
CN111947668A CN202010871052.1A CN202010871052A CN111947668A CN 111947668 A CN111947668 A CN 111947668A CN 202010871052 A CN202010871052 A CN 202010871052A CN 111947668 A CN111947668 A CN 111947668A
Authority
CN
China
Prior art keywords
detector
measurement
representing
arrival time
relative
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.)
Granted
Application number
CN202010871052.1A
Other languages
Chinese (zh)
Other versions
CN111947668B (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.)
Central South University
Original Assignee
Central South University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN202010871052.1A priority Critical patent/CN111947668B/en
Publication of CN111947668A publication Critical patent/CN111947668A/en
Application granted granted Critical
Publication of CN111947668B publication Critical patent/CN111947668B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

The invention provides a method for measuring angle/ranging combined navigation of a Jupiter detector based on-line estimation, which comprises the steps of firstly, respectively taking the position and the speed of the detector and the position and the speed of a Jupiter as system state quantities, establishing a system state model according to the track dynamics, obtaining the star light angular distance measurement through an angle measuring sensor, obtaining the pulse arrival time measurement through an X-ray pulsar detector, respectively establishing a starlight angular distance measurement model and a pulse arrival time measurement model according to the starlight angular distance measurement and the pulse arrival time measurement, estimating and correcting the position and the speed of the wooden star on line by using unscented Kalman filtering, the method inhibits the influence of the ephemeris error of the Jupiter on the navigation precision, provides high-precision position and speed estimation information for the Jupiter detector, and has important practical significance on autonomous navigation of the Jupiter detector.

Description

Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector
Technical Field
The invention relates to the technical field of autonomous navigation of deep space probes, in particular to a method for measuring angles and distance of a wooden star probe by combining navigation based on-line estimation.
Background
China launches the first Mars exploration task-Tianqu I in China in 7, 23 months in 2020, develops the first real planet exploration in China, and plans to complete three tasks of winding, falling and patrolling once. With the development of carrier rockets and other deep space exploration technologies in China and the improvement of economic strength, detection tasks such as Venus and Jupiter and the like are developed gradually in the follow-up process.
For the planetary exploration task, the navigation accuracy has an important influence on the success or failure of the task. At present, navigation information is mainly provided for a detector through a ground measurement and control station, the method can meet the requirements of most near-earth space tasks, but when a deep space detection task at a longer distance is carried out, the ground radio measurement and control mainly has the problems of three aspects of extension during communication, navigation interruption possibly caused by interference of sunscals, celestial body shielding and the like, high operation cost and the like, and the requirement of the deep space detection task in the future on high-precision real-time navigation is difficult to meet. Therefore, it is necessary for the planetary probe to improve autonomous navigation capability of the probe.
The most mature autonomous astronomical navigation in the prior art is astronomical angle measurement navigation, and the method has the advantages of high instantaneous positioning precision and capability of providing direction information of the detector relative to a target celestial body. However, the farther the distance between the probe and the celestial body is, the lower the positioning accuracy of the angle measurement navigation is, and further, this method cannot directly provide information on the distance of the probe with respect to the target celestial body. The pulsar ranging navigation can provide high-precision positioning information for the detector, and the navigation precision is not influenced by the position between the detector and the celestial body. However, for the muxing detection task, because the current measurement level is limited, the muxing ephemeris contains a large error with the magnitude of about tens of kilometers, so that the position of the muxing relative to the sun cannot be accurately obtained, and the navigation precision is affected.
Disclosure of Invention
The invention provides a method for measuring angle/ranging combined navigation of a wooden star detector based on-line estimation, and aims to solve the problem of influence of a wooden star ephemeris error on navigation precision.
In order to achieve the above object, an embodiment of the present invention provides an online estimation-based method for integrated angle measurement/distance measurement navigation of a wooden star detector, including:
step 1, taking the position and speed of a detector and the position and speed of a wooden star as system state quantities, and establishing a system state model according to track dynamics;
step 2, obtaining star light angular distance measurement through an angle measuring sensor;
step 3, establishing a starlight angular distance measurement model according to the starlight angular distance measurement;
step 4, obtaining the pulse arrival time measurement through an X-ray pulsar detector;
step 5, establishing a pulse arrival time measurement model according to the pulse arrival time measurement;
and 6, estimating and correcting the position and the speed of the Jupiter on line through unscented Kalman filtering to obtain the position and speed estimation information of the detector.
Wherein, the step 1 specifically comprises:
taking the position and speed of the detector and the position and speed of the wooden star as system state quantities, as follows:
Figure BDA0002651120020000021
wherein, Xp=[rp vp]T,rpIndicating the position of the probe relative to the wooden star, vpRepresenting the velocity vector, X, of the detector relative to the Jupiterj=[rj vj]T,rjIndicating the position of the Jupiter relative to the sun, vjRepresenting the velocity vector of the muxing relative to the sun.
Wherein, the step 1 further comprises:
the system state equation, as follows:
Figure BDA0002651120020000022
wherein r ispIndicating the position of the probe relative to the wooden star, vpRepresenting the velocity vector of the probe relative to the Jupiter, rjIndicating the position of the Jupiter relative to the sun, vjRepresenting the velocity vector of the muxing relative to the sun,
Figure BDA0002651120020000023
are respectively rp、vp、rj、vjThe derivative of (1), g represents the 2 norm of the vector, g does not count as much as the vector3Represents the cube of g | |, mujDenotes the gravitational constant, μ, of JupitersDenotes the gravitational constant of the sun, rpsRepresenting the position vector of the detector relative to the sun, rsj=rp-rpsRepresenting the position vector of the sun relative to the Jupiter, wpRepresenting process noise, w, caused by various disturbances to the detectorjRepresenting process noise caused by various disturbances to which the wooden star is subjected;
formula (2) is represented as follows:
Figure BDA0002651120020000031
wherein,
Figure BDA0002651120020000032
the derivative of the state quantity X is represented,
Figure BDA0002651120020000033
representing time t
Figure BDA0002651120020000034
f (x (t), t) represents the system nonlinear state transfer function, w ═ 0 wp 0 wj]TRepresenting systematic process noise vectors, w (t) representing time tw。
Wherein, step 2 specifically includes:
obtaining the starlight angular distance between the detector and the wooden star and between the detector and the background fixed star by using the angular measurement sensor, and establishing a starlight angular distance measurement model by taking the starlight angular distance as the measurement, wherein the measurement model is as follows:
Figure BDA0002651120020000035
wherein alpha is1、α2And alpha3Respectively representing the angular distances r between the detector and the wooden star and three background starspIndicating the position of the probe relative to the wooden star, s1、s2And s3Respectively representing the direction vectors of three background stars under the inertial system.
Wherein, step 3 specifically includes:
measuring the angular separation of the stars as a quantity Z1=[α1 α2 α3]TEstablishing an expression of a starlight angular distance measurement model, which is as follows:
Z1=h1[X(t),t]+V1(t) (5)
wherein h is1[X(t),t]Non-linear continuous measurement function, V, representing angular separation of stars1(t) represents the measurement noise of the angular separation of the starlight at time t.
Wherein, step 4 specifically includes:
an X-ray pulsar detector is used to obtain a pulse arrival time measurement, and the pulse arrival time is used as a measurement to establish a pulse arrival time measurement model as follows:
Figure BDA0002651120020000041
wherein,
Figure BDA0002651120020000042
representing the time for the ith pulsar pulse to reach the solar system centroid,
Figure BDA0002651120020000043
represents the time of arrival of the ith pulsar pulse at the detector, b represents the position vector of the solar system centroid relative to the sun, c represents the speed of light, niRepresents the direction vector of the ith pulsar under the inertial system,
Figure BDA0002651120020000048
represents the distance, mu, from the ith pulsar to the center of mass of the solar systemsDenotes the gravitational constant of the sun, rpsRepresenting the position vector of the detector relative to the sun, rpIndicating the position of the probe relative to the wooden star, rjIndicating the position of the muxing relative to the sun.
Wherein, step 5 specifically includes:
measuring pulse arrival time as a quantity
Figure BDA0002651120020000044
An expression of a pulse arrival time measurement model is established as follows:
Figure BDA0002651120020000045
wherein,
Figure BDA0002651120020000046
representing the time for the ith pulsar pulse to reach the solar system centroid,
Figure BDA0002651120020000047
representing the time of arrival of the ith pulsar pulse at the detector, h2[X(t),t]Non-linear continuous measurement function, V, representing pulse arrival time2(t) represents the measurement noise of the pulse arrival time at time t.
Wherein, step 6 specifically includes:
when no pulse arrival time measurement is carried out at the filtering moment, a state estimation and an error covariance estimation of the detector are obtained through a system state model and a starlight angular distance measurement model in a fixed filtering period and by utilizing unscented Kalman filtering; and when the pulse arrival time at the filtering moment is measured, combining the system state model, the starlight angular distance measurement model and the pulse arrival time measurement model, and obtaining the state estimation and the error covariance estimation of the detector through unscented Kalman filtering to obtain navigation information.
The principle of the invention is as follows: due to the fact that the muxing ephemeris has large errors, the position and the speed of the muxing are used as states and expanded to state quantities, the position and the speed of the muxing are estimated on line through the star light angular distance measurement and the pulse arrival time measurement, the accuracy of the position and the speed of the muxing is improved, and the navigation precision is improved.
The scheme of the invention has the following beneficial effects:
according to the method for measuring the angle/distance of the wooden star detector based on the on-line estimation, the high-precision autonomous navigation of the deep space detector is achieved, the direction information of the detector relative to the wooden star is obtained through the star light angle distance measurement, meanwhile, the distance information of the detector relative to the sun is obtained through the pulse arrival time measurement, the influence of the ephemeris error of the wooden star on the accuracy of the measurement information is eliminated through the on-line estimation of the position and the speed of the wooden star, and the navigation precision is improved.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a flow chart of the inventive combined navigation method based on-line estimation for measuring angle/distance of the Mars probe;
FIG. 3 is a schematic diagram of the X-ray pulsar navigation principle of the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantages of the present invention more apparent, the following detailed description is given with reference to the accompanying drawings and specific embodiments.
The invention provides a method for measuring angles/ranging combined navigation of a Jupiter detector based on-line estimation, aiming at the problem of influence of ephemeris error of the conventional Jupiter on navigation precision.
As shown in fig. 1 to 3, an embodiment of the present invention provides a method for integrated angle/distance measurement and navigation of a wooden star finder based on online estimation, including: step 1, taking the position and speed of a detector and the position and speed of a wooden star as system state quantities, and establishing a system state model according to track dynamics; step 2, obtaining star light angular distance measurement through an angle measuring sensor; step 3, establishing a starlight angular distance measurement model according to the starlight angular distance measurement; step 4, obtaining the pulse arrival time measurement through an X-ray pulsar detector; step 5, establishing a pulse arrival time measurement model according to the pulse arrival time measurement; and 6, estimating and correcting the position and the speed of the Jupiter on line through unscented Kalman filtering to obtain the position and speed estimation information of the detector.
According to the method for measuring the angle/distance of the wooden star detector based on the online estimation, the position and the speed of the wooden star are used as the state to be expanded and maintained into the state quantity, the position and the speed of the wooden star are estimated online through the measurement of the star light angular distance and the measurement of the pulse arrival time, the accuracy of the position and the speed of the wooden star is improved, and the navigation precision is improved.
Wherein, the step 1 specifically comprises: taking the position and speed of the detector and the position and speed of the wooden star as system state quantities, as follows:
Figure BDA0002651120020000061
wherein, Xp=[rp vp]T,rpIndicating the position of the probe relative to the wooden star, vpRepresenting the velocity vector, X, of the detector relative to the Jupiterj=[rj vj]T,rjIndicating the position of the Jupiter relative to the sun, vjRepresenting the velocity vector of the muxing relative to the sun.
Wherein, the step 1 further comprises: the system state equation, as follows:
Figure BDA0002651120020000062
wherein r ispIndicating the position of the probe relative to the wooden star, vpRepresenting the velocity vector of the probe relative to the Jupiter, rjIndicating the position of the Jupiter relative to the sun, vjRepresenting the velocity vector of the muxing relative to the sun,
Figure BDA0002651120020000063
are respectively rp、vp、rj、vjThe derivative of (1), g represents the 2 norm of the vector, g does not count as much as the vector3Represents the cube of g | |, mujDenotes the gravitational constant, μ, of JupitersDenotes the gravitational constant of the sun, rpsRepresenting the position vector of the detector relative to the sun, rsj=rp-rpsRepresenting the position vector of the sun relative to the Jupiter, wpRepresenting process noise, w, caused by various disturbances to the detectorjRepresenting process noise caused by various disturbances to which the wooden star is subjected;
formula (2) is represented as follows:
Figure BDA0002651120020000064
wherein,
Figure BDA0002651120020000065
the derivative of the state quantity X is represented,
Figure BDA0002651120020000066
representing time t
Figure BDA0002651120020000067
f (x (t), t) represents the system nonlinear state transfer function, w ═ 0 wp 0 wj]TRepresenting the system process noise vector, w (t) represents w at time t.
Wherein, step 2 specifically includes: obtaining the starlight angular distance between the detector and the wooden star and between the detector and the background fixed star by using the angular measurement sensor, and establishing a starlight angular distance measurement model by taking the starlight angular distance as the measurement, wherein the measurement model is as follows:
Figure BDA0002651120020000068
wherein alpha is1、α2And alpha3Respectively representing the angular distances r between the detector and the wooden star and three background starspIndicating the position of the probe relative to the wooden star, s1、s2And s3Respectively representing the direction vectors of three background stars under the inertial system.
Wherein, step 3 specifically includes: measuring the angular separation of the stars as a quantity Z1=[α1 α2 α3]TEstablishing an expression of a starlight angular distance measurement model, which is as follows:
Z1=h1[X(t),t]+V1(t) (5)
wherein h is1[X(t),t]Non-linear continuous measurement function, V, representing angular separation of stars1(t) represents the measurement noise of the angular separation of the starlight at time t.
Wherein, step 4 specifically includes: an X-ray pulsar detector is used to obtain a pulse arrival time measurement, and the pulse arrival time is used as a measurement to establish a pulse arrival time measurement model as follows:
Figure BDA0002651120020000071
wherein,
Figure BDA0002651120020000072
representing the time for the ith pulsar pulse to reach the solar system centroid,
Figure BDA0002651120020000073
represents the time of arrival of the ith pulsar pulse at the detector, b represents the position vector of the solar system centroid relative to the sun, c represents the speed of light, niRepresents the direction vector of the ith pulsar under the inertial system,
Figure BDA0002651120020000074
represents the distance, mu, from the ith pulsar to the center of mass of the solar systemsDenotes the gravitational constant of the sun, rpsRepresenting the position vector of the detector relative to the sun, rpIndicating the position of the probe relative to the wooden star, rjIndicating the position of the muxing relative to the sun.
According to the method for measuring angle/distance by combining navigation based on online estimation and used for the wooden star detector, the accurate position estimation r of the detector relative to the wooden star can be obtained according to the formula (6)pAccurate knowledge of the location vector r of Jupiter relative to the sun is requiredjDirectly using the muxing ephemeris introduces muxing ephemeris error and affects the navigation precision, so r in the state quantity is usedjTo r tojAnd (6) online estimation and correction.
Wherein, step 5 specifically includes: measuring pulse arrival time as a quantity
Figure BDA0002651120020000075
An expression of a pulse arrival time measurement model is established as follows:
Figure BDA0002651120020000076
wherein,
Figure BDA0002651120020000077
representing the time for the ith pulsar pulse to reach the solar system centroid,
Figure BDA0002651120020000078
representing the time of arrival of the ith pulsar pulse at the detector, h2[X(t),t]Non-linear continuous measurement function, V, representing pulse arrival time2(t) represents the measurement noise of the pulse arrival time at time t.
Wherein, step 6 specifically includes: when no pulse arrival time measurement is carried out at the filtering moment, a state estimation and an error covariance estimation of the detector are obtained through a system state model and a starlight angular distance measurement model in a fixed filtering period and by utilizing unscented Kalman filtering; and when the pulse arrival time at the filtering moment is measured, combining the system state model, the starlight angular distance measurement model and the pulse arrival time measurement model, and obtaining the state estimation and the error covariance estimation of the detector through unscented Kalman filtering to obtain navigation information.
In the method for measuring angle/distance of a wooden star detector based on online estimation, the position and speed of the detector and the position and speed of a wooden star are respectively used as system state quantities, a system state model is established according to track dynamics, star angular distance measurement is obtained through an angle measuring sensor, pulse arrival time measurement is obtained through an X-ray pulsar detector, a star angular distance measurement model and a pulse arrival time measurement model are respectively established according to the star angular distance measurement and the pulse arrival time measurement, and when no pulse arrival time measurement is measured at the filtering moment, the state estimation and error covariance estimation of the detector are obtained through the system state model and the star angular distance measurement model in a fixed filtering period and through unscented Kalman filtering; when the pulse arrival time at the filtering moment is measured, the state estimation and error covariance estimation of the detector are obtained through the unscented Kalman filtering by combining a system state model, a starlight angular distance measurement model and a pulse arrival time measurement model, and the position and speed estimation information of the detector is obtained to obtain navigation information.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (8)

1. A method for measuring angle/distance of a wooden star detector and combining navigation based on-line estimation is characterized by comprising the following steps:
step 1, taking the position and speed of a detector and the position and speed of a wooden star as system state quantities, and establishing a system state model according to track dynamics;
step 2, obtaining star light angular distance measurement through an angle measuring sensor;
step 3, establishing a starlight angular distance measurement model according to the starlight angular distance measurement;
step 4, obtaining the pulse arrival time measurement through an X-ray pulsar detector;
step 5, establishing a pulse arrival time measurement model according to the pulse arrival time measurement;
and 6, estimating and correcting the position and the speed of the Jupiter on line through unscented Kalman filtering to obtain the position and speed estimation information of the detector.
2. The method for integrated muxing detector angle measurement/distance measurement navigation based on-line estimation as recited in claim 1, wherein the step 1 specifically comprises:
taking the position and speed of the detector and the position and speed of the wooden star as system state quantities, as follows:
Figure FDA0002651120010000011
wherein, Xp=[rp vp]T,rpIndicating the position of the probe relative to the wooden star, vpRepresenting the velocity vector, X, of the detector relative to the Jupiterj=[rj vj]T,rjIndicating the position of the Jupiter relative to the sun, vjRepresenting the velocity vector of the muxing relative to the sun.
3. The method for integrated muxing detector angle/distance measurement navigation based on-line estimation as recited in claim 2, wherein the step 1 further comprises:
the system state equation, as follows:
Figure FDA0002651120010000012
wherein r ispIndicating the position of the probe relative to the wooden star, vpRepresenting the velocity vector of the probe relative to the Jupiter, rjIndicating the position of the Jupiter relative to the sun, vjRepresenting the velocity vector of the muxing relative to the sun,
Figure FDA0002651120010000021
are respectively rp、vp、rj、vjThe derivative of (1), g represents the 2 norm of the vector, g does not count as much as the vector3Represents the cube of g | |, mujDenotes the gravitational constant, μ, of JupitersDenotes the gravitational constant of the sun, rpsRepresenting the position vector of the detector relative to the sun, rsj=rp-rpsRepresenting the position vector of the sun relative to the Jupiter, wpRepresenting process noise, w, caused by various disturbances to the detectorjRepresenting process noise caused by various disturbances to which the wooden star is subjected;
formula (2) is represented as follows:
Figure FDA0002651120010000022
wherein,
Figure FDA0002651120010000023
the derivative of the state quantity X is represented,
Figure FDA0002651120010000024
representing time t
Figure FDA0002651120010000025
f (x (t), t) represents the system nonlinear state transfer function, w ═ 0 wp 0 wj]TRepresenting the system process noise vector, w (t) represents w at time t.
4. The method for integrated muxing detector angle measurement/distance measurement navigation based on-line estimation as recited in claim 3, wherein the step 2 specifically comprises:
obtaining the starlight angular distance between the detector and the wooden star and between the detector and the background fixed star by using the angular measurement sensor, and establishing a starlight angular distance measurement model by taking the starlight angular distance as the measurement, wherein the measurement model is as follows:
Figure FDA0002651120010000026
wherein alpha is1、α2And alpha3Respectively representing the angular distances r between the detector and the wooden star and three background starspIndicating the position of the probe relative to the wooden star, s1、s2And s3Respectively representing the direction vectors of three background stars under the inertial system.
5. The method for integrated muxing detector angle measurement/distance measurement navigation based on-line estimation as recited in claim 4, wherein step 3 specifically comprises:
measuring the angular separation of the stars as a quantity Z1=[α1 α2 α3]TEstablishing an expression of a starlight angular distance measurement model, which is as follows:
Z1=h1[X(t),t]+V1(t) (5)
wherein h is1[X(t),t]Non-linear continuous measurement function, V, representing angular separation of stars1(t) represents the measurement noise of the angular separation of the starlight at time t.
6. The method for integrated muxing detector angle measurement/distance measurement navigation based on-line estimation as recited in claim 5, wherein the step 4 specifically comprises:
an X-ray pulsar detector is used to obtain a pulse arrival time measurement, and the pulse arrival time is used as a measurement to establish a pulse arrival time measurement model as follows:
Figure FDA0002651120010000031
wherein,
Figure FDA0002651120010000032
representing the time for the ith pulsar pulse to reach the solar system centroid,
Figure FDA0002651120010000033
represents the time of arrival of the ith pulsar pulse at the detector, b represents the position vector of the solar system centroid relative to the sun, c represents the speed of light, niRepresents the direction vector of the ith pulsar under the inertial system,
Figure FDA0002651120010000034
represents the distance, mu, from the ith pulsar to the center of mass of the solar systemsDenotes the gravitational constant of the sun, rpsRepresenting the position vector of the detector relative to the sun, rpIndicating the position of the probe relative to the wooden star, rjIndicating the position of the muxing relative to the sun.
7. The method for integrated muxing detector angle measurement/distance measurement navigation based on-line estimation as recited in claim 6, wherein the step 5 specifically comprises:
measuring pulse arrival time as a quantity
Figure FDA0002651120010000035
An expression of a pulse arrival time measurement model is established as follows:
Figure FDA0002651120010000036
wherein,
Figure FDA0002651120010000037
representing the time for the ith pulsar pulse to reach the solar system centroid,
Figure FDA0002651120010000038
representing the time of arrival of the ith pulsar pulse at the detector, h2[X(t),t]Non-linear continuous measurement function, V, representing pulse arrival time2(t) represents the measurement noise of the pulse arrival time at time t.
8. The method for integrated muxing detector angle measurement/distance measurement navigation based on-line estimation as recited in claim 1, wherein step 6 specifically comprises:
when no pulse arrival time measurement is carried out at the filtering moment, a state estimation and an error covariance estimation of the detector are obtained through a system state model and a starlight angular distance measurement model in a fixed filtering period and by utilizing unscented Kalman filtering; and when the pulse arrival time at the filtering moment is measured, combining the system state model, the starlight angular distance measurement model and the pulse arrival time measurement model, and obtaining the state estimation and the error covariance estimation of the detector through unscented Kalman filtering to obtain navigation information.
CN202010871052.1A 2020-08-26 2020-08-26 Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector Expired - Fee Related CN111947668B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010871052.1A CN111947668B (en) 2020-08-26 2020-08-26 Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010871052.1A CN111947668B (en) 2020-08-26 2020-08-26 Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector

Publications (2)

Publication Number Publication Date
CN111947668A true CN111947668A (en) 2020-11-17
CN111947668B CN111947668B (en) 2022-03-18

Family

ID=73366473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010871052.1A Expired - Fee Related CN111947668B (en) 2020-08-26 2020-08-26 Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector

Country Status (1)

Country Link
CN (1) CN111947668B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113375659A (en) * 2021-08-16 2021-09-10 中国人民解放军国防科技大学 Pulsar navigation method based on starlight angular distance measurement information

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070046530A1 (en) * 2005-08-25 2007-03-01 Honeywell International Inc. Methods and systems for satellite navigation
CN104006813A (en) * 2014-04-03 2014-08-27 中国人民解放军国防科学技术大学 Pulsar/starlight angle combination navigation method of high orbit satellite
CN106767765A (en) * 2016-11-17 2017-05-31 上海卫星工程研究所 Deep space based on position and orientation matrix intercept test the speed navigation fixed star method of selecting
CN107024211A (en) * 2017-06-22 2017-08-08 北京航空航天大学 A kind of deep space probe angle measurement/differential speed measuring/difference ranges Combinated navigation method
CN107144283A (en) * 2017-06-30 2017-09-08 上海航天控制技术研究所 A kind of high considerable degree optical pulsar hybrid navigation method for deep space probe
CN110940333A (en) * 2019-12-12 2020-03-31 中南大学 Deep space probe angle measurement and time delay combined navigation method based on online estimation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070046530A1 (en) * 2005-08-25 2007-03-01 Honeywell International Inc. Methods and systems for satellite navigation
CN104006813A (en) * 2014-04-03 2014-08-27 中国人民解放军国防科学技术大学 Pulsar/starlight angle combination navigation method of high orbit satellite
CN106767765A (en) * 2016-11-17 2017-05-31 上海卫星工程研究所 Deep space based on position and orientation matrix intercept test the speed navigation fixed star method of selecting
CN107024211A (en) * 2017-06-22 2017-08-08 北京航空航天大学 A kind of deep space probe angle measurement/differential speed measuring/difference ranges Combinated navigation method
CN107144283A (en) * 2017-06-30 2017-09-08 上海航天控制技术研究所 A kind of high considerable degree optical pulsar hybrid navigation method for deep space probe
CN110940333A (en) * 2019-12-12 2020-03-31 中南大学 Deep space probe angle measurement and time delay combined navigation method based on online estimation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
房建成 等: "深空探测器自主天文导航技术综述", 《飞控与探测》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113375659A (en) * 2021-08-16 2021-09-10 中国人民解放军国防科技大学 Pulsar navigation method based on starlight angular distance measurement information
CN113375659B (en) * 2021-08-16 2021-11-02 中国人民解放军国防科技大学 Pulsar navigation method based on starlight angular distance measurement information

Also Published As

Publication number Publication date
CN111947668B (en) 2022-03-18

Similar Documents

Publication Publication Date Title
Luo et al. A new Kalman filter-based in-motion initial alignment method for DVL-aided low-cost SINS
CN107655493B (en) SINS six-position system-level calibration method for fiber-optic gyroscope
CN103900576B (en) A kind of information fusion method of survey of deep space independent navigation
CN103017774B (en) Pulsar navigation method with single detector
CN107655485B (en) Cruise section autonomous navigation position deviation correction method
CN101692001B (en) Autonomous celestial navigation method for deep space explorer on swing-by trajectory
CN106052691B (en) Close ring error correction method in the mobile drawing of laser ranging
CN103954279B (en) Doppler's differential speed measuring model and the method with X-ray pulsar integrated navigation
CN109708663B (en) Star sensor online calibration method based on aerospace plane SINS assistance
CN104764449B (en) A kind of capture section deep space probe celestial self-navigation method based on ephemeris amendment
CN110940333B (en) Deep space probe angle measurement and time delay combined navigation method based on online estimation
CN112197765B (en) Method for realizing fine navigation of underwater robot
CN105737858A (en) Attitude parameter calibration method and attitude parameter calibration device of airborne inertial navigation system
CN111998855B (en) Geometric method and system for determining space target initial orbit through optical telescope common-view observation
CN106017480B (en) Depth Combinated navigation method towards deep space exploration capture section
CN111623785B (en) Deep space probe autonomous navigation method based on inter-satellite time delay measurement
CN106959457B (en) GLONASS almanac parameter estimation method for satellite navigation
CN111947668B (en) Online estimation-based angle measurement/distance measurement combined navigation method for wooden star detector
CN112525203A (en) Spacecraft autonomous astronomical navigation method based on angle constraint auxiliary measurement
CN105651287B (en) A kind of solar flare arrival time difference measurement and Combinated navigation method, system
CN114111804B (en) High-precision time zero point alignment method for carrier rocket multisource multiclass measurement data
CN109625331A (en) Satellite controller and satellite control method
CN113551669A (en) Short baseline-based combined navigation positioning method and device
CN103256932A (en) Replacement and extrapolation combined navigation method
CN116698048A (en) Combined navigation method based on pulsar/inter-satellite ranging/landmark

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220318

CF01 Termination of patent right due to non-payment of annual fee