CN103217161B - A kind of pulsar navigation position and velocity joint method of estimation - Google Patents
A kind of pulsar navigation position and velocity joint method of estimation Download PDFInfo
- Publication number
- CN103217161B CN103217161B CN201310062534.2A CN201310062534A CN103217161B CN 103217161 B CN103217161 B CN 103217161B CN 201310062534 A CN201310062534 A CN 201310062534A CN 103217161 B CN103217161 B CN 103217161B
- Authority
- CN
- China
- Prior art keywords
- interval
- sub
- pulse
- observation
- centerdot
- 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.)
- Expired - Fee Related
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
A kind of pulsar navigation position and velocity joint method of estimation, belong to Spacecraft Autonomous Navigation field.First observation interval is divided into multiple sub-interval by the method.Then, in the sub-interval of each observation, by pulsar radiation period superimposed pulse signal to obtain the sub-profile of pulse accumulation, adopt the sub-profile of accurate maximum Likelihood process pulse accumulation to obtain the pulse arrival time at every sub-interval.Finally, according to these pulse arrival times, the least square estimation method is utilized to estimate position and the speed of spacecraft.Location of the present invention and constant speed precision very high, close to the U.S. labor lower bound of carat, and calculated amount is very little.Therefore, the present invention has important practical significance to the Spacecraft Autonomous Navigation based on pulsar.
Description
Technical field
The invention belongs to spacecraft navigation field, particularly a kind of pulsar navigation position and velocity joint method of estimation.
Background technology
In the last few years, in order to fight for space resources, spacefaring nation carries out survey of deep space activity one after another, and accurate navigator fix is prerequisite and the basis of all survey of deep space activities.At present, spacecraft navigation system mainly comprises following several: land station's navigational system, GPS (Global Position System, GPS), gravity navigation system, earth-magnetic navigation system and celestial navigation system (CelestialNavigation System, CNS) etc.But these technology are all not suitable for survey of deep space, and subject matter is as follows: land station's navigational system, GPS, gravity navigation system and earth-magnetic navigation system can only provide navigation information for the spacecraft of the aircraft of earth surface and terrestrial space; CNS installs optical sensitive device on spacecraft, realizes location by the relative position measuring fixed star and nearly celestial body, and its positioning precision affects by the distance between spacecraft and nearly celestial body, cannot meet the requirement of interplanetary mission hi-Fix.
X-ray pulsar navigation is a kind of emerging Spacecraft Autonomous Navigation mode, can provide the navigation informations such as long-time, high-precision location and constant speed in whole space for spacecraft.X-ray pulsar is a kind of magnetic neutron star of high speed rotation, can continuous external stable radiation, foreseeable, unique X ray signal.Spaceborne X-ray detector receives the pulse signal of a period of time (5 ~ 10 minutes), to accumulate and process the time t that can obtain and arrive spacecraft to it by the recurrence interval
sC.And this pulse arrives solar system barycenter (solar system barycenter, SSB) time t
bthen obtain by pulse timing model prediction.Time delay t
b-t
sCit is then the basic observation of pulsar navigation system.Can obtain multiple observed quantity from many pulsar direction of visual lines, recycling Navigation Filter can estimate Space Vehicle position.
The acquisition of pulse arrival time need by long-time (5 ~ 10 minutes) pile-up pulse star radiation signal.But because spacecraft is motion, and cannot obtain long-time, high-precision velocity information, the pulse signal that X-ray detector receives will inevitably be subject to the impact of Doppler effect, and this effect cannot effectively be eliminated.In pulse signal accumulation, if accumulate by the pulsar radiation signal natural period, pulse accumulation profile there will be " passivation ", and positioning precision also can significantly decline, and pulsar navigation system even can be caused time serious normally to work.In upper " On pulse phase estimation and tracking of variable celestial X-ray sources (impulse phase about variable celestial X-ray source is estimated and the follows the tracks of) " literary composition announced of Institute of Navigation 63rd Annual Meeting (the 63rd annual meeting of navigation association), Golshan, according to the distortion of pulse profile, utilizes maximum Likelihood to determine position and speed.The precision of the method is close to the U.S. labor lower bound of carat.But because this technology obtains velocity estimation value by the maximal value of gridding method search likelihood function, the calculating of each net point all relates to whole impulse radiation signal raw data.Under prior art conditions, detector resolution is at musec order, and accumulated time is 5 ~ 10 minutes, then original data volume is 10
8~ 10
9magnitude, very huge.This causes calculated amount, and (only the calculating of a net point just comprises 10 very greatly
8~ 10
9magnitude additive operation and 10
8~ 10
9magnitude multiplying), so this algorithm cannot requirement of real time.At periodical " Science China:Physics, Mechanics & Astronomy (Chinese science: physics mechanics uranology) " the 6th phase in 2011 " Modeling and Doppler measurement of X-ray pulsar (X-ray pulsar modeling and Doppler measurement) " literary composition in, pulsar signal is decomposed into multiple Gaussian function, propose the concept of pulse accumulation profile entropy, with the error between full sized pules profile and accumulation profile for objective function, by adjustment time delay and Doppler speed, make objective function minimum, thus obtain the velocity estimation value of degree of precision.These two kinds of methods also need the raw data repeatedly calculating impulse radiation signal, and calculated amount is very large.
Summary of the invention
The present invention is directed to existing pulsar navigation location and constant speed method calculated amount this defect large, provide one pulsar navigation position and velocity joint method of estimation fast, to realize calculating in real time in-orbit.
Technical scheme of the present invention is a kind of pulsar navigation position and velocity joint method of estimation, comprises the following steps:
Step 1, total observation interval (t
0, t
f) be divided into the individual equal sub-interval of observation of m, (t
0+ (i-1) T
obs/ m, t
0+ iT
obs/ m), i=1,2 ... m, wherein, observation time is T
obs=t
f-t
0;
Step 2, in the sub-interval of each observation, by pulsar cycle superimposed pulse signal, obtains the sub-profile of pulse accumulation, then utilizes accurate maximum Likelihood to estimate impulse phase estimated value
according to gained impulse phase estimated value
conversion obtains pulse arrival time t
i;
Step 3, with Least Square in Processing step 2 gained pulse arrival time t
i, obtain Space Vehicle position and speed.
And, in step 2, describedly utilize accurate maximum Likelihood to estimate impulse phase estimated value
as shown in the formula,
Wherein, function h () is normalization full sized pules profile,
be i-th observation sub-interval pile-up pulse profile, θ and
represent phase place.
And, in step 2, by gained impulse phase estimated value
conversion obtains pulse arrival time t
iimplementation be, by impulse phase estimated value
as the impulse phase at i-th sub-interval of observation
following formula is adopted to change,
Wherein, n
ithe pulse number of cycles between spacecraft and solar system barycenter,
it is solar system barycenter
the phase place that moment is corresponding, P
0for the recurrence interval.
And, in step 3, obtain the estimated value of Space Vehicle position and speed by following formula
Wherein, c is the light velocity.
The present invention's advantage is compared with prior art:
(1) the present invention carrys out estimating Doppler speed according to the distortion of pulse accumulation profile, but utilizes the change of pulse arrival time to carry out estimating speed.Therefore, calculate without the need to relating to a large amount of impulse radiation photon, calculated amount greatly reduces.
(2) estimated accuracy provided by the invention is close to the U.S. labor lower bound of carat, location and constant speed precision high.
Accompanying drawing explanation
Fig. 1 is embodiments of the invention process flow diagrams.
Embodiment
Technical solution of the present invention is described in detail below in conjunction with drawings and Examples.The technical program can adopt the automatic operational scheme of computer software technology.
Embodiment is using Crab pulsar as navigation pulsar, and see Fig. 1, the flow process of embodiment comprises the following steps:
Step 1: embodiment utilizes X-ray detector to collect pulsar x-ray photon, by the observation interval segmentation in such a way that the pulsed photonic sequence pair that detector obtains is answered.Total observation interval (t
0, t
f) be divided into the individual equal sub-interval (t of observation of m
0+ (i-1) T
obs/ m, t
0+ iT
obs/ m), i=1,2 ... m, wherein, t
0for observation initial time, t
ffor the observation end time, observation time is T
obs=t
ft
0.T is set in embodiment
0=0s, t
f=300s, m=30.
The value of m need meet
wherein, c is the light velocity, and v is spacecraft speed, and δ is pulse signal resolution, and its value equals X-ray detector temporal resolution.δ=1ms, v=10km/s in embodiment.
Step 2: in the sub-interval of each observation, by pulsar radiation period superimposed pulse signal, obtains the sub-profile of pulse accumulation, then utilizes quick accurate maximum Likelihood to estimate impulse phase, by the impulse phase estimated value obtained
be converted into pulse arrival time t
i.As shown in fig. 1, there is the sub-profile 1 of accumulation, accumulate sub-profile 2 ... accumulate sub-profile m, process respectively.
Because accurate maximum Likelihood of the prior art only estimates impulse phase, its calculated amount is only relevant with profile length, instead of pulse signal photon numbers, so calculated amount is less.Can see document: Rinauro S, Colonnese S, Scarano G.Fast near-maximum likelihood phase estimation of X-ray pulsars (quick accurate maximum likelihood pulsar phase estimation) .Signal Processing (signal transacting), 2013,93 (1): 326-331. the present invention utilize accurate maximum Likelihood to improve arithmetic speed.
The present invention utilizes accurate maximum Likelihood to estimate phase place as shown in Equation 1:
Wherein, function h () is normalization full sized pules profile,
it is the pulse accumulation profile at i-th sub-interval of observation.As in Fig. 1, the 1st, 2 ... the pulse accumulation profile at m the sub-interval of observation is designated as the sub-profile 1,2 of accumulation respectively ... m.θ and
represent phase place, be variable, wherein θ is integration variable,
for independent variable.During concrete enforcement, full sized pules profile observes acquisition for a long time by land station, and then normalization makes profile integration be 1, can obtain normalization full sized pules profile.EPN database (The EuropeanPulsar Network Data Archive, Europe observations of pulsar grid database) disclose full sized pules profile and the correlation parameter of pulsar, China technician obtains full sized pules profile by this database.
Due to the motion of spacecraft, i-th sub-interval (t of observation of acquisition
0+ (i-1) T
obs/ m, t
0+ iT
obs/ m) the corresponding initial time of phase value be not t
0.Below, the present invention investigates i-th initial time observing the phase value at sub-interval corresponding of acquisition.The present invention supposes that the recurrence interval is P
0, the speed of spacecraft on pulsar direction of visual lines is v, and first impulse phase at i-th sub-interval of observation is
wherein,
therefore, a kth impulse phase at i-th sub-interval of observation is
Wherein, c is the light velocity.
The present invention supposes there be K pulse in a sub-interval of observation, and the span of k is 0,1 ... K-1.When adopting nearly Likelihood estimation, the impulse phase at i-th sub-interval of observation of acquisition
for
average, can calculate by formula 3:
This impulse phase equals the phase place of (K-1)/2 pulse, the phase value at i-th the sub-interval of observation namely obtained
corresponding to the intermediate time at i-th sub-interval of observation,
By the impulse phase at i-th sub-interval of observation
be converted into pulse arrival time t
i, as shown in Equation 4:
Wherein, n
ithe pulse number of cycles between spacecraft and solar system barycenter,
at solar system barycenter
the impulse phase that moment is corresponding.P
0for the recurrence interval.
In embodiment, estimate that such as formula utilizing nearly maximum Likelihood shown in (1) phase place obtains impulse phase estimated value
then such as formula the impulse phase estimated value will obtained (4) Suo Shi
as the impulse phase at i-th sub-interval of observation
be converted into pulse arrival time t
i.In the present embodiment, the recurrence interval P of Crab pulsar
0for 0.0334s.
Step 3: with Least Square in Processing step 2 gained pulse arrival time t
i, Space Vehicle position and speed can be obtained.
According to step 2, the present invention obtains a series of data pair,
i=1,2 ... m.
with pulse arrival time t
ibetween relation as follows:
Wherein,
with
for spacecraft is at t
0the position in moment and speed, c is the light velocity.
Therefore, position and velocity estimation value as follows:
Specific embodiment described herein is only to the explanation for example of the present invention's spirit.Those skilled in the art can make various amendment or supplement or adopt similar mode to substitute to described specific embodiment, but can't depart from spirit of the present invention or surmount the scope that appended claims defines.
Claims (3)
1. pulsar navigation position and a velocity joint method of estimation, is characterized in that: comprise the following steps,
Step 1, total observation interval (t
0, t
f) be divided into the individual equal sub-interval of observation of m, (t
0+ (i-1) T
obs/ m, t
0+ iT
obs/ m), i=1,2...m, wherein, observation time is T
obs=t
f-t
0;
The value of m need meet
wherein, c is the light velocity, and v is spacecraft speed, and δ is pulse signal resolution, and its value equals X-ray detector temporal resolution;
Step 2, in the sub-interval of each observation, by pulsar cycle superimposed pulse signal, obtains the sub-profile of pulse accumulation, then utilizes accurate maximum Likelihood to estimate impulse phase estimated value
according to gained impulse phase estimated value
conversion obtains pulse arrival time t
i;
Step 3, with Least Square in Processing step 2 gained pulse arrival time t
i, obtain Space Vehicle position and speed; Comprise the estimated value simultaneously obtaining Space Vehicle position and speed by following formula
Wherein, c is the light velocity.
2. pulsar navigation position and velocity joint method of estimation according to claim 1, is characterized in that: in step 2, describedly utilizes accurate maximum Likelihood to estimate impulse phase estimated value
as shown in the formula,
Wherein, function h () is normalization full sized pules profile,
be i-th observation sub-interval pile-up pulse profile, θ and
represent phase place.
3. pulsar navigation position and velocity joint method of estimation according to claim 2, is characterized in that: in step 2, by gained impulse phase estimated value
conversion obtains pulse arrival time t
iimplementation be, by impulse phase estimated value
as the impulse phase at i-th sub-interval of observation
following formula is adopted to change,
Wherein, n
ithe pulse number of cycles between spacecraft and solar system barycenter,
it is solar system barycenter
the phase place that moment is corresponding, P
0for the recurrence interval.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310062534.2A CN103217161B (en) | 2013-02-27 | 2013-02-27 | A kind of pulsar navigation position and velocity joint method of estimation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310062534.2A CN103217161B (en) | 2013-02-27 | 2013-02-27 | A kind of pulsar navigation position and velocity joint method of estimation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103217161A CN103217161A (en) | 2013-07-24 |
CN103217161B true CN103217161B (en) | 2015-10-21 |
Family
ID=48815165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310062534.2A Expired - Fee Related CN103217161B (en) | 2013-02-27 | 2013-02-27 | A kind of pulsar navigation position and velocity joint method of estimation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103217161B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103776454B (en) * | 2014-01-21 | 2016-08-17 | 西安电子科技大学 | Maximum likelihood phase estimation method based on X-ray pulsar |
CN103900577B (en) * | 2014-04-14 | 2016-08-17 | 武汉科技大学 | A kind of Relative Navigation towards formation flight tests the speed and Combinated navigation method |
CN103954279B (en) * | 2014-05-19 | 2016-08-17 | 武汉科技大学 | Doppler's differential speed measuring model and the method with X-ray pulsar integrated navigation |
CN104535067B (en) * | 2015-01-14 | 2017-05-24 | 中国人民解放军国防科学技术大学 | Method for quickly calculating arrival time of pulse signal based on sector search |
CN104567937B (en) * | 2015-01-27 | 2017-08-29 | 中国空间技术研究院 | A kind of high-precision X-ray pulsar signal imitation method |
CN109186612B (en) * | 2018-09-06 | 2021-09-24 | 武汉科技大学 | Rapid pulsar period estimation method based on compressed sensing and Clamei bound |
CN111649735B (en) * | 2020-06-12 | 2021-11-16 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
CN113375697B (en) * | 2021-06-23 | 2022-09-23 | 西安电子科技大学 | X-ray pulsar signal time delay estimation method based on maximum posterior estimation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101038169A (en) * | 2007-02-13 | 2007-09-19 | 北京空间飞行器总体设计部 | Navigation satellite autonomous navigation system and method based on X-ray pulsar |
CN102175246A (en) * | 2010-12-31 | 2011-09-07 | 北京航空航天大学 | Spacecraft navigation system for equivalent device of X pulsar detector |
-
2013
- 2013-02-27 CN CN201310062534.2A patent/CN103217161B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101038169A (en) * | 2007-02-13 | 2007-09-19 | 北京空间飞行器总体设计部 | Navigation satellite autonomous navigation system and method based on X-ray pulsar |
CN102175246A (en) * | 2010-12-31 | 2011-09-07 | 北京航空航天大学 | Spacecraft navigation system for equivalent device of X pulsar detector |
Non-Patent Citations (4)
Title |
---|
X射线脉冲星自主导航的光子到达时间转换;费保俊等;《空间科学学报》;20101231;第30卷(第1期);第85-90页 * |
基于泊松模型的X射线脉冲星信号的最大似然TOA估计;李建勋等;《天文学报》;20100731;第51卷(第3期);第263-270页 * |
基于非齐次Poisson过程的脉冲星到达信号的最大似然相位估计;李建勋等;《信号处理》;20100831;第26卷(第8期);第1252-1256页 * |
李建勋等.基于脉冲星定时模型的自主导航定位方法.《中国科学 G辑:物理学 力学 天文学》.2009,第39卷(第2期),第311-316页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103217161A (en) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103217161B (en) | A kind of pulsar navigation position and velocity joint method of estimation | |
CN103017774B (en) | Pulsar navigation method with single detector | |
Emadzadeh et al. | Relative navigation between two spacecraft using X-ray pulsars | |
Brunt et al. | Assessment of NASA airborne laser altimetry data using ground-based GPS data near Summit Station, Greenland | |
WO2010047875A2 (en) | Method and apparatus for autonomous, in-receiver prediction of gnss ephemerides | |
CN103776454A (en) | Maximum likelihood phase estimation method based on X-ray pulsar | |
CN105022045A (en) | Multi-source data fusion-based three-dimensional ionosphere chromatographic method | |
CN110208831A (en) | A method of realizing No. three Satellite Orbit Determinations of Beidou and time synchronization | |
CN104848862B (en) | The punctual method and system in a kind of ring fire detector precision synchronous location | |
Den Ouden et al. | Stand-alone single-frequency GPS ice velocity observations on Nordenskiöldbreen, Svalbard | |
CN104316048B (en) | Method for building universal pulsar-based autonomous navigation measurement model | |
Konrad et al. | Determining the age distribution of Colle Gnifetti, Monte Rosa, Swiss Alps, by combining ice cores, ground-penetrating radar and a simple flow model | |
CN103017772A (en) | Optical and pulsar fusion type self-navigating method based on observability analysis | |
CN102645229A (en) | Spaceborne laser altimeter in-orbit moon elevation detection uncertainty evaluation method | |
Emadzadeh et al. | A new relative navigation system based on X-ray pulsar measurements | |
Weng et al. | Advanced data assimilation for cloud-resolving hurricane initialization and prediction | |
CN104990548B (en) | Processing method based on the dynamic pulse star signal that epoch converts into | |
CN116774264A (en) | Moving target positioning method based on low orbit satellite opportunistic signal Doppler | |
CN106986049B (en) | Precise parallel optimization design method for deep space gravity-assist track | |
Wang et al. | X-ray pulsar/starlight Doppler deeply-integrated navigation method | |
CN106597421A (en) | Prediction model-based delay and delay rate fast search method applied to very long baseline antenna array | |
CN102607563B (en) | System for performing relative navigation on spacecraft based on background astronomical information | |
CN104316056B (en) | Pulsar-based navigation constellation time synchronization and orientation parameter measurement method | |
CN103940424B (en) | A kind of X-ray pulsar navigation signal integer ambiguity detection and method of estimation | |
CN106153052B (en) | A kind of construction method of the template pulse of pulsar |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151021 Termination date: 20160227 |