CN104316048A - Method for building universal pulsar-based autonomous navigation measurement model - Google Patents
Method for building universal pulsar-based autonomous navigation measurement model Download PDFInfo
- Publication number
- CN104316048A CN104316048A CN201410542106.4A CN201410542106A CN104316048A CN 104316048 A CN104316048 A CN 104316048A CN 201410542106 A CN201410542106 A CN 201410542106A CN 104316048 A CN104316048 A CN 104316048A
- Authority
- CN
- China
- Prior art keywords
- rightarrow
- pulsar
- centerdot
- delta
- solar system
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a method for building a universal pulsar-based autonomous navigation measurement model. The method comprises calculating the geometric delay of a pulsar pulse signal, the physical delay of the signal, time coordinate transformation correction and double-star orbital motion delay, respectively, and then building the universal pulsar-based autonomous navigation measurement model and a pulsar-based autonomous navigation observation equation. According to the method for building the universal pulsar-based autonomous navigation measurement model, a brand-new pulsar-based navigation measurement model is provided, which is capable of theoretically reaching the calculation accuracy of 1ns; and meanwhile, the space motion velocity of the pulsar and the time delay effect of the impulsion double-star orbital motion are taken into account.
Description
Technical field
The present invention relates to a kind of pulsar independent navigation measurement model, be applicable to all kinds of spacecraft independent navigation fields based on pulsar.
Background technology
The pulsar rotation period is highly stable, and its Spin parameters and uranometry parameter can high-precision measurings.One component is distributed in space all directions, and its locus coordinate and rotation model form inertia space-time frame of reference through the pulsar accurately measured, and can be used for all kinds of aircraft independent navigation.The basis of pulsar independent navigation builds high-precision pulsar navigation measurement model and observations of pulsar equation.At present, in the associated pulsation star independent navigation document delivered, about the measurement model of pulsar navigation, substantially the model of doctor Sheikh has all been continued to use (see Sheikh S I.the Use of Variable Celestial X Ray Sources for Spacecraft Navigation [D], Ph.D Dissertation.Dept.of Aerospace Engineering, University of Marylan, UMI, 2005).This model provides, and pulsar pulse arrives aircraft moment and the difference expression formula arriving the solar system barycenter moment, and this expression formula strictly can not describe the relation between pulsar pulse arrival aircraft moment and pulsar impulse ejection moment.Based on the pulsar navigation observation equation that Sheikh model is set up, introduce solar system barycenter inadequately relative to heliocentric position vector, do not eliminate solar system Shapiro delayed impact completely, do not comprise the single order Doppler delay that the relative solar system centroid velocity of pulsar causes, the orbital motion not relating to binary pulsar postpones correct problems.In addition, Sheikh model observes the Pulsar timing model adopted inconsistent with the radio timing of setting up pulsar ephemeris (rotation model parameter and uranometry parameter), when actual navigation application, and meeting drawing-in system error artificially.Application number be 200710191527.7 Chinese patent document disclose a kind of Combinated navigation method of Multi-information acquisition.Application number be 200910063267.4 Chinese patent document disclose a kind of pulsar/CNS Combinated navigation method.These documents, discuss the filtering algorithm problem of pulsar navigation, in principle, still based on measurement model and the observation equation of Sheikh.
Summary of the invention
In order to overcome the deficiencies in the prior art, the invention provides a kind of strict pulsar navigation measurement model of universality, this model strictly describes the theory relation between pulsar pulse arrival aircraft moment and pulsar impulse ejection moment, and is consistent with Pulsar timing observation model.Meanwhile, the navigation model provided take into account the impact of pulsar spatial movement speed and the impact of binary pulsar orbital motion, is a kind of more applicable pulsar navigation model.
The technical solution adopted for the present invention to solve the technical problems comprises the following steps:
1. build pulsar navigation measurement model
(1) geometric delay of pulsar pulse signal is calculated
Wherein, at ephemeris with reference in epoch, solar system center-of-mass system, pulsar unit vector is
relative solar system centroid distance is R
0, spatial movement speed is
observation moment satellite relative solar system barycenter vector is
observation moment and ephemeris are Δ t with reference to difference epoch, and the light velocity is c;
(2) physical delay of signal is calculated
Wherein, PB
sSbe the main celestial body number of the solar system, G is Newtonian gravitational constant, M
ka kth celestial body in solar system quality,
the vector of the relative kth of a satellite celestial body in solar system, r
kbe
mould, m
Θsolar mass, r
Θbe satellite to sun distance, ψ is the subtended angle of the sun and pulsar relative satellite;
(3) computing time, coordinate conversion corrected Δ t
c
Wherein t
0corresponding JD=2443144.5003725, t
tCGthe TCG moment,
w
0terrestrial gravitation gesture and the rotation gesture sum at Geoid place,
the position vector of observation moment satellite relative to the earth's core,
be the earth's core relative to solar system velocity vector of enter of mass, t
tCBbe the TCB moment, celestial body in solar system is in the gravitational potential at the earth's core place
m
jthe quality removing extratellurian J celestial body in solar system, r
eJthe distance that the earth arrives this celestial body,
the impact of relativity higher order term,
represent all asteroid impacts of the solar system, t
bBarrive the double star barycenter moment, t with the pulsar pulse that double star barycentric coordinate time represents
posthat pulsar ephemeris is with reference to epoch;
(4) calculate double star orbital motion and postpone Δ t
ob, Δ t
obcomprise Roemer delay, Einstein postpones, Shapiro postpones and aberration postpones;
(5) the pulsar independent navigation measurement model of universality
The clock correction supposing satellite clock is δ t, and it is t that the X-ray detector be arranged on satellite observes the pulsar pulse obtained arrive the satellite moment within integral time
obs, the impulse ejection moment of pulsar rotation model prediction is expressed as t
p, closing therebetween is t
p-t
obs=-Δ
g-Δ
sh+ Δ t
c-Δ t
ob+ δ t, obtains the pulsar navigation measurement model of universality:
2. build pulsar independent navigation observation equation
The invention has the beneficial effects as follows: give brand-new pulsar navigation measurement model, model comprises pulsar pulse geometric delay, physical delay and time coordinate conversion and calculates, and in theory, model can reach 1ns computational accuracy.Meanwhile, model take into account pulsar spatial movement speed and binary pulsar orbital motion time delay effect, is a kind of more applicable pulsar navigation model.
Embodiment
Below in conjunction with embodiment, the present invention is further described, the present invention includes but be not limited only to following embodiment.
The pulsar independent navigation measurement model of universality, refers to the measurement model of all kinds of spacecraft independent navigations based on X-ray pulse Solo and bi-orbiters.Pulsar navigation measurement model describes to measure the pulsar pulse that obtains and arrive the spacecraft moment (TOA) and the theoretical model of relation between the pulsar impulse ejection moment.When building pulsar navigation measurement model, we are decomposed into the vacuum geometric delay (for X ray wave band, Dispersion of Media postpone can ignore) of pulsar pulse signal, physical delay and reference frame conversion it and calculate.To the calculating of this three, sets forth analytical expression and algorithm, all can reach the precision being better than 1ns in theory.For binary pulsar, navigation measurement model must calculate the delay that pulsar orbital motion causes.Pulsar orbital motion postpones to comprise geometric delay, companion star's gravitation postpones, time coordinate conversion and the aberration that causes relative to double star barycenter transverse velocity due to pulsar postpone.The calculating that binary pulsar orbital motion postpones is based on single order post newton method orbital motion.
On the pulsar independent navigation measurement model basis building universality, give the observation equation of pulsar navigation.Pulsar navigation observation equation is relevant with pulsar pulse TOA measuring method.Consider aircraft relative to pulsar motion to the impact building high s/n ratio pulse profile, we utilize the position vector adopted value of pulsar navigation measurement model and observation moment aircraft, each x-ray photon TOA that X-ray detector measurement is obtained, convert the photo emissions moment of pulsar to, thus pulsar position (not being solar system barycenter), build high s/n ratio integrated pulse profile.Again by comparing with the full sized pules profile of pulsar position, obtain the pulsar impulse ejection moment observed.In this case, observations of pulsar equation describes, the relation between the same impulse ejection moment that impulse ejection moment of pulsar rotation model prediction and observing on board the aircraft obtains.The difference of the two is the error of position of aircraft vector adopted value and the function of aircraft clock clock correction.
Below for the navigation of the earth satellite of apply pulse star, the construction method of pulsar navigation measurement model is described:
1. build pulsar navigation measurement model
(1) geometric delay of signal
Signal Geometrical propagation postpones mainly accurate Calculation and observes moment pulsar relative to the position vector of spacecraft, and for reaching 1ns precision, computing formula must comprise single order item and the second order term of pulsar spatial movement speed.Suppose that in solar system center-of-mass system, pulsar (or binary pulsar barycenter) unit vector is at ephemeris (comprising pulsar Spin parameters and uranometry parameter) with reference to epoch
relative solar system centroid distance is R
0, spatial movement speed is
(can be analyzed to proper motion in right ascension, proper motion in declination and radial velocity component), observation moment satellite relative solar system barycenter vector is
observation moment and ephemeris are Δ t with reference to difference epoch, and the light velocity is c, are ignoring signal propagation R
0in constant range delay situation, the geometric delay of pulsar pulse signal can calculate with following formula.
Formula (1) is also the precise calculation that solar system Roemer postpones, and includes speed
doppler delay, even if for the larger pulsar of movement velocity, formula (1) can reach the computational accuracy of 1ns.Formula (1) have ignored
and above item, for pulsar navigation application, ignore the error of item much smaller than 1ns.
(2) physical delay of signal
Signal physical delay refers to that solar system Shapiro postpones, and the expression formula only comprising single order item is as follows:
On the right of above formula, Section 1 is constant term, is ignored, and consider the second order term (only calculating sun impact) that Shapiro postpones, the actual Shapiro Delay computing formula adopted is
In formula (2), PB
sSbe the main celestial body number of the solar system, G is Newtonian gravitational constant, M
ka kth celestial body in solar system quality,
the vector of the relative kth of a satellite celestial body in solar system, r
kbe
mould.M
Θsolar mass, r
Θbe satellite to sun distance, ψ is the subtended angle of the sun and pulsar relative satellite.Formula (2) is consistent with the correlation formula of the radio timing model determining pulsar ephemeris.
(3) reference frame conversion
The fundamental reference system of pulsar navigation is solar system barycenter inertial reference system, and for the purpose of convenience of calculation, pulsar navigation theoretical model needs the conversion carrying out reference frame to calculate.In reference frame conversion calculates, relativity transformer effect of space length is less than 1ns, can ignore.The most complicated is the calculating that between different reference frame, time coordinate is changed, and corrects also referred to as Einstein.Suppose that earth satellite clock is synchronous with terrestrial time TT, need calculating observation moment TT to geocentric coordinate time (TCG) conversion, TCG to solar system barycentric coordinate time (TCB) conversion and TCB to the transfer problem of binary pulsar barycentric coordinate time.TT and TCG closes
TCG-TT=L
G(t
TCG-t
0) (3)
Wherein t
0corresponding JD=2443144.5003725, t
tCGthe TCG moment,
w
0being terrestrial gravitation gesture and the rotation gesture sum at Geoid place, is known quantity.TCG and TCB closes
In formula (4)
the position vector of observation moment satellite relative to the earth's core,
be the earth's core relative to solar system velocity vector of enter of mass, provided by position of the earth ephemeris, t
tCBbe the TCB moment, celestial body in solar system is in the gravitational potential at the earth's core place
m
jthe quality of extratellurian J the celestial body in solar system of removing (not containing asteroid), r
eJit is the distance that the earth arrives this celestial body.In integrand
be the impact of relativity higher order term, the impact of relativity higher order term can be analyzed to linear term and periodic term, and its periodic term is negligible.
represent all asteroid impacts of the solar system, asteroid impact can be analyzed to linear term and periodic term equally, and its periodic term is very little, can ignore completely.Usually
with
utilize U.S. JPL celestial body in solar system ephemeris to analyze and research separately to determine, their value is respectively
with
formula (4) integration item must utilize JPL DE ephemeris to adopt numerical integrating to calculate, integration starting point t
0corresponding to JD=2443144.5003725.Numerical integration computational accuracy reaches 0.1ns.The result of calculation of numerical integration is referred to as ephemeris between terrestrial time, and between terrestrial time, ephemeris forms the basis of TCG to TCB conversion, can directly utilize ephemeris between terrestrial time computed in advance to carry out interpolation when practical application.It should be noted that, containing U in integrand
eand v
ethe integral result of item can be decomposed into linear term and periodic term two parts, makes its linear term coefficient be Δ L
c, according to the analysis result of DE405, have Δ L
c=1.48082685594 × 10
-8.Order
here, L
caccurately determine it is necessary, because L
b=L
c+ L
g(vide infra) is as the linear coefficient of linear transformation between TDB and TCB.At L
cafter determining, some authors define ephemeris between terrestrial time and only include periodic term part, in this case, when time coordinate conversion calculates, utilize the interpolation results of time ephemeris must add all linear components and correct, be i.e. L
c(t
tCB-t
0) value.At present, the Tempo2 software systems of radio timing observation, adopt ephemeris between the terrestrial time that calculates based on DE405 ephemeris.Upgrade the high-precision JPL ephemeris of version if adopted, likely obtain ephemeris between more high-precision terrestrial time.Obviously, terrestrial time TT is transformed into the calculating of barycentric coordinate time TCB, should be above-mentioned formula (3) and formula (4) sum.For keep high precision, TT be transformed into TCG and TCG be transformed into TCB calculate must adopt iterative algorithm.TCB is transformed into for TCG, the time variable on formula (4) the right is TCB instead of TCG yardstick, the known TCG moment should be first adopted to calculate TCB, using the TCB that calculates first as time variable, again new TCB is calculated by formula (4) again, to obtain high-precision conversion result of calculation.In general, only need iteration once just passable.
According to International Astronomical association (IAU) resolution in 2006, barycentric dynamical time (TDB) TDB and barycentric coordinate time TCB has following relation,
TDB=TCB-L
B×(JD
TCB-T
0)×86400+TDB
0 (5)
In formula (5), T
0=2443144.5003725JD,
L
B=L
C+L
G=1.550519768×10
-8,
TDB
0=-6.55×10
-5s。
The TDB of such definition can think the time variable of JPL DE405.Needing to be utilized by known TCB formula (5) to convert TDB to, the coordinate figure exported by DE405 when reading DE405, needing to be multiplied by K again
-1, to be transformed into SI units system.Here, scale factor K=1+L
b.When time coordinate conversion calculates employing redaction JPL ephemeris, ephemeris, L between the terrestrial time of its correspondence
cand TDB
0the possibility of result Deng constant and DE405 has little difference, based on redaction JPL ephemeris terrestrial time between the analysis and research work of ephemeris and dependent constant thereof be necessary.
Suppose that the relative velocity of solar system barycenter and binary pulsar barycenter is
barycentric coordinate time TCB be transformed into double star barycentric coordinate time (for pulse list star, be then TCB be transformed into pulsar former time), need to consider special relativity time dilation, have
T in formula (6)
bBthat the pulsar pulse represented with double star barycentric coordinate time arrives the double star barycenter moment (pulse arrives the solar system barycenter moment and is transformed into same pulse and arrives the double star barycenter moment, during calculating, usually neglects pulsar distance R
0corresponding light time).T
posthat pulsar ephemeris is with reference to epoch.Owing to calculating Δ
eStime, we do not know t
bB, only know that pulse arrives the TCB moment of double star barycenter, therefore formula (6) needs to adopt iterative algorithm.First calculate Δ with the TCB moment of pulse arrival double star barycenter
eS, the TCB moment then pulse being arrived double star barycenter adds the Δ that first time calculates
eSas t
bB, then recalculate Δ
eS.General iteration is once just much of that.
Make Δ t
crepresent that time coordinate conversion is corrected, Δ t
cequal above-mentioned formula (3), (4) and (6) sum.
(4) binary pulsar orbital motion postpones
In general, binary pulsar timing model parameter provides the Keplerian orbit parameter of binary pulsar: projection semi-major axis of orbit, orbital period, excentricity, periastron time and periastron argument, on different accuracy, also provide relativity parameter of double star orbital motion simultaneously, as: the advance of periastron, orbital period derivative, pulsar and companion star's quality etc.Based on general relativity, application kepler parameters and double star mass parameter can also derive other relevant parameters.Utilize these parameters, we can calculate double star orbital motion easily and postpone Δ t
ob, Δ t
obcomprise: Roemer postpones, Einstein postpones, Shapiro postpones and aberration postpones.Document " Tempo2:Timing model and precision " summarizes binary pulsar orbital motion achievement in research, describe binary pulsar timing model and the precision thereof of the employing of Tempo2 software systems, give computing formula and the algorithm of all kinds double star orbital motion delay, Δ t can be directly applied to
obcalculating.
(5) the pulsar independent navigation measurement model of universality
The clock correction supposing satellite clock is δ t, and be arranged on X-ray detector on satellite within certain integral time, observing the pulsar pulse obtained arrive the satellite moment (TOA) is t
obs, the impulse ejection moment of pulsar rotation model prediction is expressed as t
p, pass is therebetween
t
p-t
obs=-Δ
g-Δ
Sh+Δt
C-Δt
ob+δt (7)
Formula (1) and (2) are substituted into formula (7), Wo Menyou
In formula (8), the meaning of each symbol is the same, Δ t
cthat time coordinate conversion is corrected, Δ t
obthat double star track postpones.
In actual navigation application, the apparent position vector of satellite
known,
be
error, then
Formula (9) is substituted into formula (8), and by nonlinear terms linearization, obtains
In formula (10), each symbol with " ~ " represent satellite apparent position vector or by Proximity Vector
the physical quantity calculated.
formula (3), (4) and (6) sum,
utilize observation when the moment, satellite clock was former and satellite position adopted value calculates.
be that binary pulsar orbital motion postpones, calculated by double star barycentric coordinate time, double star orbit parameter and satellite position adopted value.
unknown quantity to be asked with δ t.If navigation pulsar is single star, about the physical quantity of double star is omitted in formula (1).Formula (10) is the pulsar navigation measurement model of universality.
2. build pulsar independent navigation observation equation
Pulse TOA to be compared with full sized pules profile by observation integrated pulse profile and determines.Each x-ray photon TOA in the X-ray detector observational record lower integral time, if build integrated pulse profile by photon TOA at satellite position, due to the motion of satellite relative pulse star, can make integrated pulse profile fog.For improving signal to noise ratio (S/N ratio), we utilize the removing of formula (10) the right
other all known terms sum with beyond δ t item, is transformed into pulsar x-ray photon x time by each x-ray photon TOA, builds integrated pulse profile in pulsar position, by comparing with full sized pules profile, obtains the impulse ejection moment t observed
pobs.So have
Formula (11) is the observation equation of pulsar navigation.
Pulsar independent navigation measurement model and observation equation build for earth satellite navigation, in principle, are also applicable to the independent navigation application of deep space exploration aircraft.Such as, the artificial satellite around Mars of apply pulse star navigates, and just time coordinate conversion calculating is different from earth satellite.Need to be transformed into Mars barycentric coordinate time around during Martian satellite former, and then be transformed into solar system barycentric coordinate time (TCB).Mars barycentric coordinate time is transformed into TCB and calculates, and needs to utilize JPL DE series celestial body in solar system ephemeris, adopts numerical integrating method, precalculates and obtain Mars time ephemeris.The construction method of Mars time ephemeris, be quite analogous to geocentric coordinate time (TCG) be transformed into TCB adopt terrestrial time between ephemeris (see " reference frame conversion " trifle).
Claims (1)
1. a pulsar independent navigation measurement model construction method for universality, is characterized in that comprising the steps:
(1) geometric delay of pulsar pulse signal is calculated
Wherein, at ephemeris with reference in epoch, solar system center-of-mass system, pulsar unit vector is
relative solar system centroid distance is R
0, spatial movement speed is
observation moment satellite relative solar system barycenter vector is
observation moment and ephemeris are Δ t with reference to difference epoch, and the light velocity is c;
(2) physical delay of signal is calculated
Wherein, PB
sSbe the main celestial body number of the solar system, G is Newtonian gravitational constant, M
ka kth celestial body in solar system quality,
the vector of the relative kth of a satellite celestial body in solar system, r
kbe
mould, m
Θsolar mass, r
Θbe satellite to sun distance, ψ is the subtended angle of the sun and pulsar relative satellite;
(3) computing time, coordinate conversion corrected Δ t
c
Wherein t
0corresponding JD=2443144.5003725, t
tCGthe TCG moment,
w
0terrestrial gravitation gesture and the rotation gesture sum at Geoid place,
the position vector of observation moment satellite relative to the earth's core,
be the earth's core relative to solar system velocity vector of enter of mass, t
tCBbe the TCB moment, celestial body in solar system is in the gravitational potential at the earth's core place
m
jthe quality removing extratellurian J celestial body in solar system, r
eJthe distance that the earth arrives this celestial body,
the impact of relativity higher order term,
represent all asteroid impacts of the solar system, t
bBarrive the double star barycenter moment, t with the pulsar pulse that double star barycentric coordinate time represents
posthat pulsar ephemeris is with reference to epoch;
(4) calculate double star orbital motion and postpone Δ t
ob, Δ t
obcomprise Roemer delay, Einstein postpones, Shapiro postpones and aberration postpones;
(5) the pulsar independent navigation measurement model of universality
The clock correction supposing satellite clock is δ t, and it is t that the X-ray detector be arranged on satellite observes the pulsar pulse obtained arrive the satellite moment within integral time
obs, the impulse ejection moment of pulsar rotation model prediction is expressed as t
p, closing therebetween is t
p-t
obs=-Δ
g-Δ
sh+ Δ t
c-Δ t
ob+ δ t, obtains the pulsar navigation measurement model of universality:
(6) pulsar independent navigation observation equation is built
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410542106.4A CN104316048B (en) | 2014-10-14 | 2014-10-14 | Method for building universal pulsar-based autonomous navigation measurement model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410542106.4A CN104316048B (en) | 2014-10-14 | 2014-10-14 | Method for building universal pulsar-based autonomous navigation measurement model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104316048A true CN104316048A (en) | 2015-01-28 |
CN104316048B CN104316048B (en) | 2017-02-15 |
Family
ID=52371314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410542106.4A Expired - Fee Related CN104316048B (en) | 2014-10-14 | 2014-10-14 | Method for building universal pulsar-based autonomous navigation measurement model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104316048B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105547303A (en) * | 2016-03-16 | 2016-05-04 | 中国人民解放军国防科学技术大学 | Autonomous navigation method for libration point constellation |
CN105953791A (en) * | 2016-04-27 | 2016-09-21 | 清华大学 | Single detector X-ray pulsar navigation time-sharing observation method and device |
CN107421533A (en) * | 2017-06-22 | 2017-12-01 | 北京航空航天大学 | A kind of deep space probe X-ray pulsar TOA/DTOA Combinated navigation methods |
CN107655485A (en) * | 2017-09-25 | 2018-02-02 | 北京理工大学 | A kind of cruise section independent navigation position deviation modification method |
CN108981750A (en) * | 2018-07-10 | 2018-12-11 | 西安电子科技大学 | X-ray pulse double star photon sequence emulation mode |
CN114543817A (en) * | 2021-12-31 | 2022-05-27 | 南京航空航天大学 | Spacecraft navigation method and system |
CN116609816A (en) * | 2023-07-19 | 2023-08-18 | 山东大学 | Deep space multisource elastic fusion navigation method and system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005043189A (en) * | 2003-07-28 | 2005-02-17 | National Institute Of Information & Communication Technology | Space navigation system and method therefor |
CN101038169A (en) * | 2007-02-13 | 2007-09-19 | 北京空间飞行器总体设计部 | Navigation satellite autonomous navigation system and method based on X-ray pulsar |
CN103017774A (en) * | 2012-12-24 | 2013-04-03 | 中国科学院国家授时中心 | Pulsar navigation method with single detector |
CN103033188A (en) * | 2012-12-24 | 2013-04-10 | 中国科学院国家授时中心 | Navigation satellite autonomous time synchronization method based on synthetic aperture observation |
CN103217162A (en) * | 2013-03-21 | 2013-07-24 | 西安电子科技大学 | Pulsar accumulated pulse profile time delay measurement method by sparse representation |
CN103674032A (en) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | Satellite autonomous navigation system and method integrating pulsar radiation vector and timing observation |
-
2014
- 2014-10-14 CN CN201410542106.4A patent/CN104316048B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005043189A (en) * | 2003-07-28 | 2005-02-17 | National Institute Of Information & Communication Technology | Space navigation system and method therefor |
CN101038169A (en) * | 2007-02-13 | 2007-09-19 | 北京空间飞行器总体设计部 | Navigation satellite autonomous navigation system and method based on X-ray pulsar |
CN103674032A (en) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | Satellite autonomous navigation system and method integrating pulsar radiation vector and timing observation |
CN103017774A (en) * | 2012-12-24 | 2013-04-03 | 中国科学院国家授时中心 | Pulsar navigation method with single detector |
CN103033188A (en) * | 2012-12-24 | 2013-04-10 | 中国科学院国家授时中心 | Navigation satellite autonomous time synchronization method based on synthetic aperture observation |
CN103217162A (en) * | 2013-03-21 | 2013-07-24 | 西安电子科技大学 | Pulsar accumulated pulse profile time delay measurement method by sparse representation |
Non-Patent Citations (2)
Title |
---|
杨廷高 等: "脉冲星时间尺度及其 TOA 预报初步分析", 《时间频率学报》 * |
杨廷高: "关于脉冲星脉冲到达时间转换方程", 《时间频率学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105547303A (en) * | 2016-03-16 | 2016-05-04 | 中国人民解放军国防科学技术大学 | Autonomous navigation method for libration point constellation |
CN105547303B (en) * | 2016-03-16 | 2018-06-29 | 中国人民解放军国防科学技术大学 | A kind of autonomous navigation method of libration point constellation |
CN105953791B (en) * | 2016-04-27 | 2018-12-11 | 清华大学 | Simple detector X-ray pulsar navigation timesharing observation method and device |
CN105953791A (en) * | 2016-04-27 | 2016-09-21 | 清华大学 | Single detector X-ray pulsar navigation time-sharing observation method and device |
CN107421533A (en) * | 2017-06-22 | 2017-12-01 | 北京航空航天大学 | A kind of deep space probe X-ray pulsar TOA/DTOA Combinated navigation methods |
CN107421533B (en) * | 2017-06-22 | 2019-07-30 | 北京航空航天大学 | A kind of deep space probe X-ray pulsar TOA/DTOA Combinated navigation method |
CN107655485A (en) * | 2017-09-25 | 2018-02-02 | 北京理工大学 | A kind of cruise section independent navigation position deviation modification method |
CN107655485B (en) * | 2017-09-25 | 2020-06-16 | 北京理工大学 | Cruise section autonomous navigation position deviation correction method |
CN108981750A (en) * | 2018-07-10 | 2018-12-11 | 西安电子科技大学 | X-ray pulse double star photon sequence emulation mode |
CN108981750B (en) * | 2018-07-10 | 2020-04-28 | 西安电子科技大学 | X-ray pulse double-star photon sequence simulation method |
CN114543817A (en) * | 2021-12-31 | 2022-05-27 | 南京航空航天大学 | Spacecraft navigation method and system |
CN116609816A (en) * | 2023-07-19 | 2023-08-18 | 山东大学 | Deep space multisource elastic fusion navigation method and system |
CN116609816B (en) * | 2023-07-19 | 2023-11-03 | 山东大学 | Deep space multisource elastic fusion navigation method and system |
Also Published As
Publication number | Publication date |
---|---|
CN104316048B (en) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103017774B (en) | Pulsar navigation method with single detector | |
CN104316048A (en) | Method for building universal pulsar-based autonomous navigation measurement model | |
Liu et al. | Novel algorithm for X-ray pulsar navigation against Doppler effects | |
Emadzadeh et al. | Relative navigation between two spacecraft using X-ray pulsars | |
Deng et al. | Interplanetary spacecraft navigation using pulsars | |
CN103674032B (en) | Merge the autonomous navigation of satellite system and method for pulsar radiation vector timing observation | |
CN105203101B (en) | A kind of deep space probe capture section astronomical navigation method based on target celestial body ephemeris amendment | |
CN107024211B (en) | A kind of deep space probe angle measurement/differential speed measuring/difference ranges Combinated navigation method | |
CN104457705A (en) | Initial orbit determination method for deep space target celestial body based on space-based autonomous optical observation | |
Liu et al. | Pulsar navigation for interplanetary missions using CV model and ASUKF | |
CN103217161B (en) | A kind of pulsar navigation position and velocity joint method of estimation | |
Liu et al. | Doppler/XNAV–integrated navigation system using small-area X-ray sensor | |
CN107144283A (en) | A kind of high considerable degree optical pulsar hybrid navigation method for deep space probe | |
Wang et al. | Pulse phase estimation of X-ray pulsar with the aid of vehicle orbital dynamics | |
Kai et al. | Performance enhancement of X-ray pulsar navigation using autonomous optical sensor | |
CN113375659A (en) | Pulsar navigation method based on starlight angular distance measurement information | |
Gu et al. | Optical/radio/pulsars integrated navigation for Mars orbiter | |
Feng et al. | Autonomous orbit determination and its error analysis for deep space using X-ray pulsar | |
Liu et al. | Closed-loop EKF-based pulsar navigation for Mars explorer with Doppler effects | |
Wang et al. | X-ray pulsar/starlight Doppler deeply-integrated navigation method | |
CN101915926B (en) | Method for confirming one-way distance in satellite laser ranging (SLR) | |
CN103047986A (en) | Large-scale space-time and on-orbit dynamic effect simulation method | |
CN102607563B (en) | System for performing relative navigation on spacecraft based on background astronomical information | |
Gui et al. | A novel celestial aided time-differenced pulsar navigation method against ephemeris error of Jupiter for Jupiter exploration | |
CN116698048A (en) | Combined navigation method based on pulsar/inter-satellite ranging/landmark |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 Termination date: 20171014 |
|
CF01 | Termination of patent right due to non-payment of annual fee |