Navigation satellites generally use satellite-ground and inter-satellite observation data for precise orbit determination.In orbit determination,the satellite position is often referenced to the satellite’s centroid,...Navigation satellites generally use satellite-ground and inter-satellite observation data for precise orbit determination.In orbit determination,the satellite position is often referenced to the satellite’s centroid,while the observational measurements are referenced to the satellite’s antenna phase center.The deviation between the satellite’s centroid and the antenna phase center forms the satellite antenna phase center error,which affects the precision of orbit determination.This paper takes a global navigation satellite system(GNSS)MEO satellite as an example and analyzes the actual situation of the satellite antenna phase center deviation and phase center variation based on the ground calibration data of the in-orbit satellite antenna phase center and the satellite’s in-orbit working status.The analysis shows that the antenna phase center variation caused by the satellite’s in-orbit operation is only at the centimeter level,which has a limited impact on orbit determination accuracy.The main source of precise orbit determination error is the satellite antenna phase center deviation,which can be corrected using ground calibration data.展开更多
Angles-only relative orbit determination for space non-cooperative targets based on passive sensor is subject to weakly observable problem of the relative state between two spacecraft. Previously, the evidence for ang...Angles-only relative orbit determination for space non-cooperative targets based on passive sensor is subject to weakly observable problem of the relative state between two spacecraft. Previously, the evidence for angles-only observability was found by using cylindrical dynamics, however, the solution of orbit determination is still not provided. This study develops a relative orbit determination algorithm with the cylindrical dynamics based on differential evolution. Firstly, the relative motion dynamics and line-of-sight measurement model for nearcircular orbit are established in cylindrical coordinate system.Secondly, the observability is qualitatively analyzed by using the dynamics and measurement model where the unobservable geometry is found. Then, the angles-only relative orbit determination problem is modeled into an optimal searching frame and an improved differential evolution algorithm is introduced to solve the problem. Finally, the proposed algorithm is verified and tested by a set of numerical simulations in the context of highEarth and low-Earth cases. The results show that initial relative orbit determination(IROD) solution with an appropriate accuracy in a relative short span is achieved, which can be used to initialize the navigation filter.展开更多
The HY-2 satellite carrying a satellite-borne GPS receiver is the first Chinese radar altimeter satellite, whose radial orbit determination precision must reach the centimeter level. Now HY-2 is in the test phase so t...The HY-2 satellite carrying a satellite-borne GPS receiver is the first Chinese radar altimeter satellite, whose radial orbit determination precision must reach the centimeter level. Now HY-2 is in the test phase so that the observations are not openly released. In order to study the precise orbit determination precision and procedure for HY-2 based on the satellite- borne GPS technique, the satellite-borne GPS data are simulated in this paper. The HY-2 satellite-borne GPS antenna can receive at least seven GPS satellites each epoch, which can validate the GPS receiver and antenna design. What＇s more, the precise orbit determination processing flow is given and precise orbit determination experiments are conducted using the HY-2-borne GPS data with both the reduced-dynamic method and the kinematic geometry method. With the 1 and 3 mm phase data random errors, the radial orbit determination precision can achieve the centimeter level using these two methods and the kinematic orbit accuracy is slightly lower than that of the reduced-dynamic orbit. The earth gravity field model is an important factor which seriously affects the precise orbit determination of altimeter satellites. The reduced-dynamic orbit determination experiments are made with different earth gravity field models, such as EIGEN2, EGM96, TEG4, and GEMT3. Using a large number of high precision satellite-bome GPS data, the HY-2 precise orbit determination can reach the centimeter level with commonly used earth gravity field models up to above 50 degrees and orders.展开更多
Low-Earth-Orbit（LEO） formation-flying satellites have been widely applied in many kinds of space geodesy. Precise Relative Orbit Determination（PROD） is an essential prerequisite for the LEO formation-flying satell...Low-Earth-Orbit（LEO） formation-flying satellites have been widely applied in many kinds of space geodesy. Precise Relative Orbit Determination（PROD） is an essential prerequisite for the LEO formation-flying satellites to complete their mission in space. The contribution of the BeiDou Navigation Satellite System（BDS） to the accuracy and reliability of PROD of LEO formation-flying satellites based on a Global Positioning System（GPS） is studied using a simulation method. Firstly, when BDS is added to GPS, the mean number of visible satellites increases from9.71 to 21.58. Secondly, the results show that the 3-Dimensional（3 D） accuracy of PROD, based on BDS-only, GPS-only and BDS ＋ GPS, is 0.74 mm, 0.66 mm and 0.52 mm, respectively. When BDS co-works with GPS, the accuracy increases by 29.73%. Geostationary-Earth-Orbit（GEO） satellites and Inclined Geosynchronous-Orbit（IGSO） satellites are only distributed over the Asia-Pacific region; however, they could provide a global improvement to PROD. The difference in PROD results between the Asia-Pacific region and the non-Asia-Pacific region is not apparent. Furthermore, the value of the Ambiguity Dilution Of Precision（ADOP）, based on BDS ＋ GPS, decreases by 7.50% and 8.26%, respectively, compared with BDS-only and GPS-only. Finally, if the relative position between satellites is only a few kilometres, the effect of ephemeris errors on PROD could be ignored. However, for a several-hundred-kilometre separation of the LEO satellites, the SingleDifference（SD） ephemeris errors of GEO satellites would be on the order of centimetres. The experimental results show that when IGSO satellites and Medium-Earth-Orbit（MEO） satellites co-work with GEO satellites, the accuracy decreases by 17.02%.展开更多
The precision orbit determination of challenging minisatellite payload(CHAMP) satellite was done based on position and navigation data analyst(PANDA) software which is developed in Wuhan University, using the onboard ...The precision orbit determination of challenging minisatellite payload(CHAMP) satellite was done based on position and navigation data analyst(PANDA) software which is developed in Wuhan University, using the onboard GPS data of year 2002 from day 126 to 131. The orbit accuracy was assessed by analyzing the difference from GFZ post-processed science orbits (PSO), the GPS carrier and pseudo-range data residuals and the satellite laser ranging (SLR) residuals.展开更多
At present, Global Navigation Satellite Systems（GNSS） users usually eliminate the influence of ionospheric delay of the first order items by dual-frequency ionosphere-free combination. But there is still residual io...At present, Global Navigation Satellite Systems（GNSS） users usually eliminate the influence of ionospheric delay of the first order items by dual-frequency ionosphere-free combination. But there is still residual ionospheric delay error of higher order term. The influence of the higher-order ionospheric corrections on both GPS precision orbit determination and static Precise Point Positioning（PPP） are studied in this paper. The influence of higher-order corrections on GPS precision orbit determination, GPS observations and static PPP are analyzed by neglecting or considering the higher-order ionospheric corrections by using a globally distributed network which is composed of International GNSS Service（IGS） tracking stations. Numerical experimental results show that, the root mean square（RMS） in three dimensions of satellite orbit is 36.6 mme35.5 mm. The maximal second-order ionospheric correction is 9 cm, and the maximal third-order ionospheric correction is 1 cm. Higher-order corrections are influenced by latitude and station distribution. PPP is within 3 mm in the directions of east and up. Furthermore, the impact is mainly visible in the direction of north, showing a southward migration trend, especially at the lower latitudes where the influence value is likely to be bigger than 3 mm.展开更多
Autonomous orbit determination via integration of epoch-differenced gravity gradients and starlight refraction is proposed in this paper for low-Earth-orbiting satellites operating in GPSdenied environments. Starlight...Autonomous orbit determination via integration of epoch-differenced gravity gradients and starlight refraction is proposed in this paper for low-Earth-orbiting satellites operating in GPSdenied environments. Starlight refraction compensates for the significant along-track position error that occurs from only using gravity gradients and benefits from integration in terms of improved accuracy in radial and cross-track position estimates. The between-epoch differencing of gravity gradients is employed to eliminate slowly varying measurement biases and noise near the orbit revolution frequency. The refraction angle measurements are directly used and its Jacobian matrix derived from an implicit observation equation. An information fusion filter based on a sequential extended Kalman filter is developed for the orbit determination. Truth-model simulations are used to test the performance of the algorithm, and the effects of differencing intervals and orbital heights are analyzed. A semi-simulation study using actual gravity gradient data from the Gravity field and steady-state Ocean Circulation Explorer（GOCE） combined with simulated starlight refraction measurements is further conducted, and a three-dimensional position accuracy of better than 100 m is achieved.展开更多
We study on reduced dynamic orbit determination using differenced phase in adjacent epochs for spacebome dual-frequency GPS. This method not only overcomes the shortcomings that the epoch-difference kinematic method c...We study on reduced dynamic orbit determination using differenced phase in adjacent epochs for spacebome dual-frequency GPS. This method not only overcomes the shortcomings that the epoch-difference kinematic method cannot be used when observation geometry is poor or observations are insufficient, but also avoids solving the ambiguity in the zero-difference reduced dynamic method. As the epoch-difference method is not sensitive to the impact of phase cycle slips, it can lower the difficulty of slip detection in phase observation preprocessing. In the solution strategies, we solve the high-dimensional matrix computation problems by decomposing the long observation arc into a number of short arcs. By gravity recovery and climate experiment （GRACE） satellite orbit determination and compared with GeoForschungsZentrum （GFZ） post science orbit, for epoch-difference reduced dynamic method, the root mean squares （RMSs） of radial, transverse and normal components are 1.92 cm, 3.83 cm and 3.80 cm, and the RMS in three dimensions is 5.76 cm. The solution＇s accuracy is comparable to the zero-difference reduced dynamic method.展开更多
When formation flying spacecrafts are used as platform to gain earth oriented observation, precise baselines between these spacecrafts are always essential. Gravity recovery and climate experiment （GRACE） mission is...When formation flying spacecrafts are used as platform to gain earth oriented observation, precise baselines between these spacecrafts are always essential. Gravity recovery and climate experiment （GRACE） mission is aimed at mapping the global gravity field and its variation. Accurate baseline of GRACE satellites is necessary for the gravity field modeling. The determination of kinematic and reduced dynamic relative orbits of twin satellites has been studied in this paper, and an accuracy of 2 mm for dynamic relative orbits and 5 mm for kinematic ones can be obtained, whereby most of the double difference onboard GPS ambiguities are resolved.展开更多
Solar radiation pressure is the main driving force and error source for precision orbit determination of navigation satellites.It is proportional to the solar irradiance,which is the"sun constant".In regular...Solar radiation pressure is the main driving force and error source for precision orbit determination of navigation satellites.It is proportional to the solar irradiance,which is the"sun constant".In regular calculation,the"solar constant"is regard as a constant.However,due to the existence of sunspots,flares,etc.,the solar constant is not fixed,the change in the year is about 1%.To investigate the variation of solar irradiance,we use interpolation and average segment modeling of total solar irradiance data of SORCE,establishing variance solar radiation pressure(VARSRP)model and average solar radiation pressure(AVESRP)model based on the built solar pressure model(SRPM)(constant model).According to observation data of global positioning system(GPS)and Beidou system(BDS)in 2015 and comparing the solar pressure acceleration of VARSRP,AVESRP and SRPM,the magnitude of change can reach 10-10 m/s^2.In addition,according to the satellite precise orbit determination,for GPS satellites,the results of VARSRP and AVESRP are slightly smaller than those of the SRPM model,and the improvement is between 0.1 to 0.5 mm.For geosynchronous orbit(GEO)satellites of BDS,The AVESRP and VARSRP have an improvement of 3.5 mm and 4.0 mm,respectively,based on overlapping arc,and SLR check results show the AVESRP model and the VARSRP model is improved by 2.3 mm and 3.5 mm,respectively.Moreover,the change of inclined geosynchronous orbit(IGSO)satellites and medium earth orbit(MEO)satellites is relatively small,and the improvement is smaller than 0.5 mm.展开更多
Using the FengYun-3C(FY-3C)onboard BeiDou Navigation Satellite System(BDS)and Global Positioning System(GPS)data from 2013 to 2017,this study investigates the performance and contribution of BDS to precise orbit deter...Using the FengYun-3C(FY-3C)onboard BeiDou Navigation Satellite System(BDS)and Global Positioning System(GPS)data from 2013 to 2017,this study investigates the performance and contribution of BDS to precise orbit determination(POD)for a low-Earth orbit(LEO).The overlap comparison result indicates that code bias correction of BDS can improve the POD accuracy by 12.4%.The multi-year averaged one-dimensional(1D)root mean square(RMS)of the overlapping orbit differences(OODs)for the GPS-only solution is 2.0,1.7,and 1.5 cm,respectively,during the 2013,2015,and 2017 periods.The 1D RMS for the BDS-only solution is 150.9,115.0,and 47.4 cm,respectively,during the 2013,2015,and 2017 periods,which is much worse than the GPS-only solution due to the regional system of BDS and the few BDS channels of the FY-3C receiver.For the BDS and GPS combined solution(also known as the GC combined solution),the averaged 1D RMS is 2.5,2.3,and 1.6 cm,respectively,in 2013,2015,and 2017,while the GC combined POD presents a significant accuracy improvement after the exclusion of geostationary Earth orbit(GEO)satellites.The main reason for the improvement seen after this exclusion is the unfavorable satellite tracking geometry and poor orbit accuracy of GEO satellites.The accuracy of BDS-only and GC combined solutions have gradually improved from 2013 to 2017,thanks to improvements in the accuracy of International GNSS Service(IGS)orbit and clock products in recent years,especially the availability of a high-frequency satellite clock product(30 s sampling interval)since 2015.Moreover,the GC POD(without GEO)was able to achieve slightly better accuracy than the GPS-only POD in 2017,indicating that the fusion of BDS and GPS observations can improve the accuracy of LEO POD.GC combined POD can significantly improve the reliability of LEO POD,simply due to system redundancy.An increased contribution of BDS to LEO POD can be expected with the launch of more BDS satellites and with further improvements in the accuracy of BDS satellite products in the near future.展开更多
Geomagnetic orbit determination fits for nanosatellites which pursue low cost and high-density ratio,but one of its disadvantages is the poor position accuracy introduced by magnetic bias.Here,a new method,named the f...Geomagnetic orbit determination fits for nanosatellites which pursue low cost and high-density ratio,but one of its disadvantages is the poor position accuracy introduced by magnetic bias.Here,a new method,named the fuzzy regulating unscented Kalman filter(FRUKF),is proposed.The magnetic bias is regarded as a random walk model,and a fuzzy regulator is designed to estimate the magnetic bias more accurately.The input of the regulator is the derivative of magnetic bias estimated from unscented Kalman filter(UKF).According to the fuzzy rule,the process noise covariance is adaptively determined.The FRUKF is evaluated using the real-flight data of the SWARMA.The experimental results show that the root-mean-square(RMS)position error is 3.1 km and the convergence time is shorter than the traditional way.展开更多
The Tianhui-202(TH02-02)satellite formation,as a supplement to the microwave mapping satellite system Tianhui-201(TH02-01),is the first Interferometric Synthetic Aperture Radar(InSAR)satellite formation-flying system ...The Tianhui-202(TH02-02)satellite formation,as a supplement to the microwave mapping satellite system Tianhui-201(TH02-01),is the first Interferometric Synthetic Aperture Radar(InSAR)satellite formation-flying system that supports the tracking of BeiDou global navigation Satellite system(BDS3)new B1C and B2a signals.Meanwhile,the twin TH02-02 satellites also support the tracking of Global Positioning System(GPS)L1&L2 and BDS B1I&B3I signals.As the spaceborne receiver employs two independent boards to track the Global Navigation Satellite System(GNSS)satellites,we design an orbit determination strategy by estimating independent receiver clock offsets epoch by epoch for each GNSS to realize the multi-GNSS data fusion from different boards.The performance of the spaceborne receiver is evaluated and the contribution of BDS3 to the kinematic and reduced-dynamic Precise Orbit Determination(POD)of TH02-02 satellites is investigated.The tracking data onboard shows that the average number of available BDS3 and GPS satellites are 8.7 and 9.1,respectively.The carrier-to-noise ratio and carrier phase noise of BDS3 B1C and B2a signals are comparable to those of GPS.However,strong azimuth-related systematic biases are recognized in the pseudorange multipath errors of B1C and B3I.The pseudorange noise of BDS3 signals is better than that of GPS after eliminating the multipath errors from specific signals.Taking the GPS-based reduced-dynamic orbit with single-receiver ambiguity fixing technique as a reference,the results of BDS3-only and BDS3+GPS combined POD are assessed.The Root Mean Square(RMS)of orbit comparison of BDS3-based kinematic and reduced-dynamic POD with reference orbit are better than 7 cm and 3 cm in three-Dimensional direction(3D).The POD performance based on B1C&B2a data is comparable to that based on B1I&B3I.The precision of BDS3+GPS combined kinematic orbit can reach up to 3 cm(3D RMS),which has a more than 25%improvement relative to the GPS-only solution.In addition,the consistency between the BDS3+GPS combined reduced-dynamic orbit and the GPS-based ambiguity-fixed orbit is better than 1.5 cm(3D RMS).展开更多
In recent years,the large Low Earth Orbit(LEO)constellations have become a hot topic due to their great potential to improve the Global Navigation Satellite Systems(GNSS)positioning performance.One of the important fo...In recent years,the large Low Earth Orbit(LEO)constellations have become a hot topic due to their great potential to improve the Global Navigation Satellite Systems(GNSS)positioning performance.One of the important focus is how to obtain the accurate and reliable orbits for these constellations with dozens of LEO satellites.The GNSS-based Precise Orbit Determination(POD)will be exclusively performed to achieve this goal,where the Integer Ambiguity Resolution(IAR)plays a key role in acquiring high-quality orbits.In this study,we present a comprehensive analysis of the benefit of the single-receiver IAR in LEO POD and discuss its implication for the future LEO constellations.We perform ambiguity-fixed LEO POD for four typical missions,including Gravity Recovery and Climate Experiment(GRACE)Follow-On(GRACE-FO),Swarm,Jason-3 and Sentinel-3,using the Uncalibrated Phase Delay(UPD)products generated by our GREAT(GNSS+REsearch,Application and Teaching)software.The results show that the ambiguity fixing processing can significantly improve the accuracy of LEO orbits.There are negligible differences between our UPD-based ambiguity-fixed orbits and those based on the Observable Signal Bias(OSB)and Integer Recovery Clock(IRC)products,indicating the good-quality of UPD products we generated.Compared to the float solution,the fixed solution presents a better consistency with the external precise science orbits and the largest accuracy improvement of 5 mm is achieved for GRACE-FO satellites.Meanwhile,the benefit can be observed in laser ranging residuals as well,with a Standard Deviation(STD)reduction of 3–4 mm on average for the fixed solutions.Apart from the absolute orbits,the relative accuracy of the space baseline is also improved by 20–30%in the fixed solutions.The result demonstrates the superior performance of the ambiguity-fixed LEO POD,which appears as a particularly promising technique for POD of future LEO constellations.展开更多
The two line elements(TLEs),released by the North American Aerospace Defense Command(NORAD),are chosen for CubeSats' mission operators.Unfortunately,they have errors and other accompanied problems,which cause larg...The two line elements(TLEs),released by the North American Aerospace Defense Command(NORAD),are chosen for CubeSats' mission operators.Unfortunately,they have errors and other accompanied problems,which cause large deviations in the in-track component.When a TLE value is available at a certain epoch,the dominant error is the angular error.It is proposed to correct the angular error by solving-for the mean argument of latitude at the desired epoch.A batch least squares technique and range rate measurements are used for the correction process.With the assistance of satellite tool kit(STK)software and Matlab,a simulation to verify the orbit determination(OD)technique is implemented.This paper provides an angular correction low cost OD method and presents a complete analysis for various test cases.This approach maintains high accuracy in cross-track and radial and makes great improvement in in-track at the same time,but it is exclusive for circular orbits.When it is applied to an elliptical orbit,the error will be unacceptable.Therefore,the angular error is corrected using the longitude of periapsis which totally mitigates the error at the epoch under consideration.For inclinations less than 20 o,the mean longitude is preferred for the angular correction as it provides more accuracy compared with the mean argument of latitude.展开更多
The visibility for low earth orbit（LEO） satellites provided by the BeiDou-2 system is analyzed and compared with the global positioning system（GPS）. In addition, the spaceborne receivers＇ observations are simulat...The visibility for low earth orbit（LEO） satellites provided by the BeiDou-2 system is analyzed and compared with the global positioning system（GPS）. In addition, the spaceborne receivers＇ observations are simulated by the BeiDou satellites broadcast ephemeris and LEO satellites orbits. The precise orbit determination（POD） results show that the along-track component accuracy is much better over the service area than the non-service area, while the accuracy of the other two directions keeps at the same level over different areas. However, the 3-dimensional（3D） accuracy over the two areas shows almost no difference. Only taking into consideration the observation noise and navigation satellite ephemeris errors, the 3D accuracy of the POD is about30 cm. As for the precise relative orbit determination（PROD）, the 3D accuracy is much better over the eastern hemisphere than that of the western hemisphere. The baseline length accuracy is 3.4 mm over the service area, and it is still better than 1 cm over the non-service area. This paper demonstrates that the BeiDou regional constellation could provide global service to LEO satellites for the POD and the PROD. Finally, the benefit of geostationary earth orbit（GEO） satellites is illustrated for POD.展开更多
Based on the orbit integration and orbit fitting method, the influence of the characters of the gravity model, with different precisions, on the movement of low Earth orbit satellites was studied. The way and the effe...Based on the orbit integration and orbit fitting method, the influence of the characters of the gravity model, with different precisions, on the movement of low Earth orbit satellites was studied. The way and the effect of absorbing the influence of gravity model error on CHAMP and GRACE satellite orbits, using linear and periodical empirical acceleration models and the so-called "pseudo-stochastic pulses" model, were also analyzed.展开更多
Satellite-to-Satellite tricking (SST) data can be used to determine the orbits of spacecraft in two ways. One is combined orbit determination, which combines SST data with ground-based tracking data and exploits the ...Satellite-to-Satellite tricking (SST) data can be used to determine the orbits of spacecraft in two ways. One is combined orbit determination, which combines SST data with ground-based tracking data and exploits the enhanced tracking geometry. The other is the autonomous orbit determination, which uses only SST. The latter only fits some particular circumstances since it suffers the rank defect problem in other circumstances. The proof of this statement is presented. The nature of the problem is also investigated in order to find an effective solution. Several. methods of solution are discussed. The feasibility of the methods is demonstrated by their application to a simulation.展开更多
Autonomous navigation of navigation satellite is discussed. The method of auto-orbit determination using the erosslink range and orientation parameters constraining is put forward. On the basis of the analysis of its ...Autonomous navigation of navigation satellite is discussed. The method of auto-orbit determination using the erosslink range and orientation parameters constraining is put forward. On the basis of the analysis of its feasibility, some useful conclusions are given.展开更多
For the two newly launched satellites(PRN number 27 and 28) of the future global BeiDou navigation satellite system(BDS-3), there is no available broadcast ephemeris data and other initial orbit information, but the i...For the two newly launched satellites(PRN number 27 and 28) of the future global BeiDou navigation satellite system(BDS-3), there is no available broadcast ephemeris data and other initial orbit information, but the initial orbit is the fundamental of the comprehensive analysis of the satellites and their signals. Precise orbit determination(POD) also requires determination of a priori initial value with a certain precision in order to avoid problems such as filter divergence during POD. Compared with the Newton iteration method, which relies on the initial value, this study utilizes the Bancroft algorithm to directly solve the nonlinear equations with the advantage of numerical stability. The initial orbits of these two satellites are calculated based on new code signals, and their results are analyzed and discussed. The experimental results show that, with the exception of very few epochs, when the new code signal is utilized, the median and robust variance factor of the observed residuals computed using pseudo-range observations and the solved initial orbits are less than 4 and 2 m, respectively. It also shows that this solution can be used for rapid initial orbit recovery after maneuvers of the new BeiDou satellites.展开更多
文摘Navigation satellites generally use satellite-ground and inter-satellite observation data for precise orbit determination.In orbit determination,the satellite position is often referenced to the satellite’s centroid,while the observational measurements are referenced to the satellite’s antenna phase center.The deviation between the satellite’s centroid and the antenna phase center forms the satellite antenna phase center error,which affects the precision of orbit determination.This paper takes a global navigation satellite system(GNSS)MEO satellite as an example and analyzes the actual situation of the satellite antenna phase center deviation and phase center variation based on the ground calibration data of the in-orbit satellite antenna phase center and the satellite’s in-orbit working status.The analysis shows that the antenna phase center variation caused by the satellite’s in-orbit operation is only at the centimeter level,which has a limited impact on orbit determination accuracy.The main source of precise orbit determination error is the satellite antenna phase center deviation,which can be corrected using ground calibration data.
基金supported by the National Natural Science Foundation of China (12272168)the Foundation of Science and Technology on Space Intelligent Control Laboratory (HTKJ2023KL502015)。
文摘Angles-only relative orbit determination for space non-cooperative targets based on passive sensor is subject to weakly observable problem of the relative state between two spacecraft. Previously, the evidence for angles-only observability was found by using cylindrical dynamics, however, the solution of orbit determination is still not provided. This study develops a relative orbit determination algorithm with the cylindrical dynamics based on differential evolution. Firstly, the relative motion dynamics and line-of-sight measurement model for nearcircular orbit are established in cylindrical coordinate system.Secondly, the observability is qualitatively analyzed by using the dynamics and measurement model where the unobservable geometry is found. Then, the angles-only relative orbit determination problem is modeled into an optimal searching frame and an improved differential evolution algorithm is introduced to solve the problem. Finally, the proposed algorithm is verified and tested by a set of numerical simulations in the context of highEarth and low-Earth cases. The results show that initial relative orbit determination(IROD) solution with an appropriate accuracy in a relative short span is achieved, which can be used to initialize the navigation filter.
基金supported partially by the National Natural Science Foundation of China (Nos. 40974004 and 40974016)Key Laboratory of Dynamic Geodesy of CAS, China (No. L09-01) R&I Team Support Program and the Graduate Science and Technology Foundation of SDUST, China (No. YCA110403)
文摘The HY-2 satellite carrying a satellite-borne GPS receiver is the first Chinese radar altimeter satellite, whose radial orbit determination precision must reach the centimeter level. Now HY-2 is in the test phase so that the observations are not openly released. In order to study the precise orbit determination precision and procedure for HY-2 based on the satellite- borne GPS technique, the satellite-borne GPS data are simulated in this paper. The HY-2 satellite-borne GPS antenna can receive at least seven GPS satellites each epoch, which can validate the GPS receiver and antenna design. What＇s more, the precise orbit determination processing flow is given and precise orbit determination experiments are conducted using the HY-2-borne GPS data with both the reduced-dynamic method and the kinematic geometry method. With the 1 and 3 mm phase data random errors, the radial orbit determination precision can achieve the centimeter level using these two methods and the kinematic orbit accuracy is slightly lower than that of the reduced-dynamic orbit. The earth gravity field model is an important factor which seriously affects the precise orbit determination of altimeter satellites. The reduced-dynamic orbit determination experiments are made with different earth gravity field models, such as EIGEN2, EGM96, TEG4, and GEMT3. Using a large number of high precision satellite-bome GPS data, the HY-2 precise orbit determination can reach the centimeter level with commonly used earth gravity field models up to above 50 degrees and orders.
基金supported by the National Natural Science Foundation of China (Nos. 91438202, 61370013)
文摘Low-Earth-Orbit（LEO） formation-flying satellites have been widely applied in many kinds of space geodesy. Precise Relative Orbit Determination（PROD） is an essential prerequisite for the LEO formation-flying satellites to complete their mission in space. The contribution of the BeiDou Navigation Satellite System（BDS） to the accuracy and reliability of PROD of LEO formation-flying satellites based on a Global Positioning System（GPS） is studied using a simulation method. Firstly, when BDS is added to GPS, the mean number of visible satellites increases from9.71 to 21.58. Secondly, the results show that the 3-Dimensional（3 D） accuracy of PROD, based on BDS-only, GPS-only and BDS ＋ GPS, is 0.74 mm, 0.66 mm and 0.52 mm, respectively. When BDS co-works with GPS, the accuracy increases by 29.73%. Geostationary-Earth-Orbit（GEO） satellites and Inclined Geosynchronous-Orbit（IGSO） satellites are only distributed over the Asia-Pacific region; however, they could provide a global improvement to PROD. The difference in PROD results between the Asia-Pacific region and the non-Asia-Pacific region is not apparent. Furthermore, the value of the Ambiguity Dilution Of Precision（ADOP）, based on BDS ＋ GPS, decreases by 7.50% and 8.26%, respectively, compared with BDS-only and GPS-only. Finally, if the relative position between satellites is only a few kilometres, the effect of ephemeris errors on PROD could be ignored. However, for a several-hundred-kilometre separation of the LEO satellites, the SingleDifference（SD） ephemeris errors of GEO satellites would be on the order of centimetres. The experimental results show that when IGSO satellites and Medium-Earth-Orbit（MEO） satellites co-work with GEO satellites, the accuracy decreases by 17.02%.
文摘The precision orbit determination of challenging minisatellite payload(CHAMP) satellite was done based on position and navigation data analyst(PANDA) software which is developed in Wuhan University, using the onboard GPS data of year 2002 from day 126 to 131. The orbit accuracy was assessed by analyzing the difference from GFZ post-processed science orbits (PSO), the GPS carrier and pseudo-range data residuals and the satellite laser ranging (SLR) residuals.
基金funded by the China Natural Science Funds the National Natural Science Foundation of China (41374009)Postdoctoral Applied Research Project (2015186)
文摘At present, Global Navigation Satellite Systems（GNSS） users usually eliminate the influence of ionospheric delay of the first order items by dual-frequency ionosphere-free combination. But there is still residual ionospheric delay error of higher order term. The influence of the higher-order ionospheric corrections on both GPS precision orbit determination and static Precise Point Positioning（PPP） are studied in this paper. The influence of higher-order corrections on GPS precision orbit determination, GPS observations and static PPP are analyzed by neglecting or considering the higher-order ionospheric corrections by using a globally distributed network which is composed of International GNSS Service（IGS） tracking stations. Numerical experimental results show that, the root mean square（RMS） in three dimensions of satellite orbit is 36.6 mme35.5 mm. The maximal second-order ionospheric correction is 9 cm, and the maximal third-order ionospheric correction is 1 cm. Higher-order corrections are influenced by latitude and station distribution. PPP is within 3 mm in the directions of east and up. Furthermore, the impact is mainly visible in the direction of north, showing a southward migration trend, especially at the lower latitudes where the influence value is likely to be bigger than 3 mm.
基金supported by the National Natural Science Foundation of China (No.11002008)funded in part by Ministry of Science and Technology of China (No.2014CB845303)
文摘Autonomous orbit determination via integration of epoch-differenced gravity gradients and starlight refraction is proposed in this paper for low-Earth-orbiting satellites operating in GPSdenied environments. Starlight refraction compensates for the significant along-track position error that occurs from only using gravity gradients and benefits from integration in terms of improved accuracy in radial and cross-track position estimates. The between-epoch differencing of gravity gradients is employed to eliminate slowly varying measurement biases and noise near the orbit revolution frequency. The refraction angle measurements are directly used and its Jacobian matrix derived from an implicit observation equation. An information fusion filter based on a sequential extended Kalman filter is developed for the orbit determination. Truth-model simulations are used to test the performance of the algorithm, and the effects of differencing intervals and orbital heights are analyzed. A semi-simulation study using actual gravity gradient data from the Gravity field and steady-state Ocean Circulation Explorer（GOCE） combined with simulated starlight refraction measurements is further conducted, and a three-dimensional position accuracy of better than 100 m is achieved.
基金National Natural Science Foundation of China （61002033, 60902089） Open Research Fund of State Key Laboratory of Astronautic Dynamics （2011ADL-DW0103）
文摘We study on reduced dynamic orbit determination using differenced phase in adjacent epochs for spacebome dual-frequency GPS. This method not only overcomes the shortcomings that the epoch-difference kinematic method cannot be used when observation geometry is poor or observations are insufficient, but also avoids solving the ambiguity in the zero-difference reduced dynamic method. As the epoch-difference method is not sensitive to the impact of phase cycle slips, it can lower the difficulty of slip detection in phase observation preprocessing. In the solution strategies, we solve the high-dimensional matrix computation problems by decomposing the long observation arc into a number of short arcs. By gravity recovery and climate experiment （GRACE） satellite orbit determination and compared with GeoForschungsZentrum （GFZ） post science orbit, for epoch-difference reduced dynamic method, the root mean squares （RMSs） of radial, transverse and normal components are 1.92 cm, 3.83 cm and 3.80 cm, and the RMS in three dimensions is 5.76 cm. The solution＇s accuracy is comparable to the zero-difference reduced dynamic method.
基金Supported by the National Natural Science Foundation of China (No. 40874004, No.40504002)the National 973 Program of China (No.2006CB701301)+1 种基金the National 863 Program of China (No. 2006AA12Z326, No. 2007AA12Z345)the 111 Project(No. B07037)
文摘When formation flying spacecrafts are used as platform to gain earth oriented observation, precise baselines between these spacecrafts are always essential. Gravity recovery and climate experiment （GRACE） mission is aimed at mapping the global gravity field and its variation. Accurate baseline of GRACE satellites is necessary for the gravity field modeling. The determination of kinematic and reduced dynamic relative orbits of twin satellites has been studied in this paper, and an accuracy of 2 mm for dynamic relative orbits and 5 mm for kinematic ones can be obtained, whereby most of the double difference onboard GPS ambiguities are resolved.
基金supported by the National Key Research and Development Program of China (No.2016YFB0501405)the National Natural Science Foundation of China (No.11973073)+1 种基金the Basic Project of Ministry of Science and Technology of China (No.2015FY310200)the Shanghai Key Laboratory of Space Navigation and Position Techniques (No.06DZ22101)
文摘Solar radiation pressure is the main driving force and error source for precision orbit determination of navigation satellites.It is proportional to the solar irradiance,which is the"sun constant".In regular calculation,the"solar constant"is regard as a constant.However,due to the existence of sunspots,flares,etc.,the solar constant is not fixed,the change in the year is about 1%.To investigate the variation of solar irradiance,we use interpolation and average segment modeling of total solar irradiance data of SORCE,establishing variance solar radiation pressure(VARSRP)model and average solar radiation pressure(AVESRP)model based on the built solar pressure model(SRPM)(constant model).According to observation data of global positioning system(GPS)and Beidou system(BDS)in 2015 and comparing the solar pressure acceleration of VARSRP,AVESRP and SRPM,the magnitude of change can reach 10-10 m/s^2.In addition,according to the satellite precise orbit determination,for GPS satellites,the results of VARSRP and AVESRP are slightly smaller than those of the SRPM model,and the improvement is between 0.1 to 0.5 mm.For geosynchronous orbit(GEO)satellites of BDS,The AVESRP and VARSRP have an improvement of 3.5 mm and 4.0 mm,respectively,based on overlapping arc,and SLR check results show the AVESRP model and the VARSRP model is improved by 2.3 mm and 3.5 mm,respectively.Moreover,the change of inclined geosynchronous orbit(IGSO)satellites and medium earth orbit(MEO)satellites is relatively small,and the improvement is smaller than 0.5 mm.
基金We are very grateful to the IGS,GFZ,and WHU for providing the precise orbit and clock products of GPS and BDS.Thanks also go to the EPOS-RT/PANDA software from GFZ.This study is financially supported by the National Natural Science Foundation of China(41774030,41974027,41974029,and 41505030)the Hubei Province Natural Science Foundation of China(2018CFA081)The numerical calculations in this paper were done on the supercomputing system at the Supercomputing Center of Wuhan University.
文摘Using the FengYun-3C(FY-3C)onboard BeiDou Navigation Satellite System(BDS)and Global Positioning System(GPS)data from 2013 to 2017,this study investigates the performance and contribution of BDS to precise orbit determination(POD)for a low-Earth orbit(LEO).The overlap comparison result indicates that code bias correction of BDS can improve the POD accuracy by 12.4%.The multi-year averaged one-dimensional(1D)root mean square(RMS)of the overlapping orbit differences(OODs)for the GPS-only solution is 2.0,1.7,and 1.5 cm,respectively,during the 2013,2015,and 2017 periods.The 1D RMS for the BDS-only solution is 150.9,115.0,and 47.4 cm,respectively,during the 2013,2015,and 2017 periods,which is much worse than the GPS-only solution due to the regional system of BDS and the few BDS channels of the FY-3C receiver.For the BDS and GPS combined solution(also known as the GC combined solution),the averaged 1D RMS is 2.5,2.3,and 1.6 cm,respectively,in 2013,2015,and 2017,while the GC combined POD presents a significant accuracy improvement after the exclusion of geostationary Earth orbit(GEO)satellites.The main reason for the improvement seen after this exclusion is the unfavorable satellite tracking geometry and poor orbit accuracy of GEO satellites.The accuracy of BDS-only and GC combined solutions have gradually improved from 2013 to 2017,thanks to improvements in the accuracy of International GNSS Service(IGS)orbit and clock products in recent years,especially the availability of a high-frequency satellite clock product(30 s sampling interval)since 2015.Moreover,the GC POD(without GEO)was able to achieve slightly better accuracy than the GPS-only POD in 2017,indicating that the fusion of BDS and GPS observations can improve the accuracy of LEO POD.GC combined POD can significantly improve the reliability of LEO POD,simply due to system redundancy.An increased contribution of BDS to LEO POD can be expected with the launch of more BDS satellites and with further improvements in the accuracy of BDS satellite products in the near future.
基金supported by the National Natural Science Foundation of China(No.61673212).
文摘Geomagnetic orbit determination fits for nanosatellites which pursue low cost and high-density ratio,but one of its disadvantages is the poor position accuracy introduced by magnetic bias.Here,a new method,named the fuzzy regulating unscented Kalman filter(FRUKF),is proposed.The magnetic bias is regarded as a random walk model,and a fuzzy regulator is designed to estimate the magnetic bias more accurately.The input of the regulator is the derivative of magnetic bias estimated from unscented Kalman filter(UKF).According to the fuzzy rule,the process noise covariance is adaptively determined.The FRUKF is evaluated using the real-flight data of the SWARMA.The experimental results show that the root-mean-square(RMS)position error is 3.1 km and the convergence time is shorter than the traditional way.
基金funded by the National Natural Science Foundation of China(Nos.61803018 and 41874028)the Key Laboratory Found,China(No.6142210200105)the National Key R&D Program of China(No.2020YFA0713502).
文摘The Tianhui-202(TH02-02)satellite formation,as a supplement to the microwave mapping satellite system Tianhui-201(TH02-01),is the first Interferometric Synthetic Aperture Radar(InSAR)satellite formation-flying system that supports the tracking of BeiDou global navigation Satellite system(BDS3)new B1C and B2a signals.Meanwhile,the twin TH02-02 satellites also support the tracking of Global Positioning System(GPS)L1&L2 and BDS B1I&B3I signals.As the spaceborne receiver employs two independent boards to track the Global Navigation Satellite System(GNSS)satellites,we design an orbit determination strategy by estimating independent receiver clock offsets epoch by epoch for each GNSS to realize the multi-GNSS data fusion from different boards.The performance of the spaceborne receiver is evaluated and the contribution of BDS3 to the kinematic and reduced-dynamic Precise Orbit Determination(POD)of TH02-02 satellites is investigated.The tracking data onboard shows that the average number of available BDS3 and GPS satellites are 8.7 and 9.1,respectively.The carrier-to-noise ratio and carrier phase noise of BDS3 B1C and B2a signals are comparable to those of GPS.However,strong azimuth-related systematic biases are recognized in the pseudorange multipath errors of B1C and B3I.The pseudorange noise of BDS3 signals is better than that of GPS after eliminating the multipath errors from specific signals.Taking the GPS-based reduced-dynamic orbit with single-receiver ambiguity fixing technique as a reference,the results of BDS3-only and BDS3+GPS combined POD are assessed.The Root Mean Square(RMS)of orbit comparison of BDS3-based kinematic and reduced-dynamic POD with reference orbit are better than 7 cm and 3 cm in three-Dimensional direction(3D).The POD performance based on B1C&B2a data is comparable to that based on B1I&B3I.The precision of BDS3+GPS combined kinematic orbit can reach up to 3 cm(3D RMS),which has a more than 25%improvement relative to the GPS-only solution.In addition,the consistency between the BDS3+GPS combined reduced-dynamic orbit and the GPS-based ambiguity-fixed orbit is better than 1.5 cm(3D RMS).
基金National Natural Science Foundation of China[41974027]Sino-German mobility programme[M-0054].
文摘In recent years,the large Low Earth Orbit(LEO)constellations have become a hot topic due to their great potential to improve the Global Navigation Satellite Systems(GNSS)positioning performance.One of the important focus is how to obtain the accurate and reliable orbits for these constellations with dozens of LEO satellites.The GNSS-based Precise Orbit Determination(POD)will be exclusively performed to achieve this goal,where the Integer Ambiguity Resolution(IAR)plays a key role in acquiring high-quality orbits.In this study,we present a comprehensive analysis of the benefit of the single-receiver IAR in LEO POD and discuss its implication for the future LEO constellations.We perform ambiguity-fixed LEO POD for four typical missions,including Gravity Recovery and Climate Experiment(GRACE)Follow-On(GRACE-FO),Swarm,Jason-3 and Sentinel-3,using the Uncalibrated Phase Delay(UPD)products generated by our GREAT(GNSS+REsearch,Application and Teaching)software.The results show that the ambiguity fixing processing can significantly improve the accuracy of LEO orbits.There are negligible differences between our UPD-based ambiguity-fixed orbits and those based on the Observable Signal Bias(OSB)and Integer Recovery Clock(IRC)products,indicating the good-quality of UPD products we generated.Compared to the float solution,the fixed solution presents a better consistency with the external precise science orbits and the largest accuracy improvement of 5 mm is achieved for GRACE-FO satellites.Meanwhile,the benefit can be observed in laser ranging residuals as well,with a Standard Deviation(STD)reduction of 3–4 mm on average for the fixed solutions.Apart from the absolute orbits,the relative accuracy of the space baseline is also improved by 20–30%in the fixed solutions.The result demonstrates the superior performance of the ambiguity-fixed LEO POD,which appears as a particularly promising technique for POD of future LEO constellations.
基金supported by the Research Fund for the Doctoral Program of Higher Education of China (No.20113219110025)
文摘The two line elements(TLEs),released by the North American Aerospace Defense Command(NORAD),are chosen for CubeSats' mission operators.Unfortunately,they have errors and other accompanied problems,which cause large deviations in the in-track component.When a TLE value is available at a certain epoch,the dominant error is the angular error.It is proposed to correct the angular error by solving-for the mean argument of latitude at the desired epoch.A batch least squares technique and range rate measurements are used for the correction process.With the assistance of satellite tool kit(STK)software and Matlab,a simulation to verify the orbit determination(OD)technique is implemented.This paper provides an angular correction low cost OD method and presents a complete analysis for various test cases.This approach maintains high accuracy in cross-track and radial and makes great improvement in in-track at the same time,but it is exclusive for circular orbits.When it is applied to an elliptical orbit,the error will be unacceptable.Therefore,the angular error is corrected using the longitude of periapsis which totally mitigates the error at the epoch under consideration.For inclinations less than 20 o,the mean longitude is preferred for the angular correction as it provides more accuracy compared with the mean argument of latitude.
基金co-supported by the National Natural Science Foundation of China (Nos: 61002033, 61370013)the Program for New Century Excellent Talents in University and the Scientific Research Foundation for the Returned Overseas Chinese Scholars of Ministry of Education of China
文摘The visibility for low earth orbit（LEO） satellites provided by the BeiDou-2 system is analyzed and compared with the global positioning system（GPS）. In addition, the spaceborne receivers＇ observations are simulated by the BeiDou satellites broadcast ephemeris and LEO satellites orbits. The precise orbit determination（POD） results show that the along-track component accuracy is much better over the service area than the non-service area, while the accuracy of the other two directions keeps at the same level over different areas. However, the 3-dimensional（3D） accuracy over the two areas shows almost no difference. Only taking into consideration the observation noise and navigation satellite ephemeris errors, the 3D accuracy of the POD is about30 cm. As for the precise relative orbit determination（PROD）, the 3D accuracy is much better over the eastern hemisphere than that of the western hemisphere. The baseline length accuracy is 3.4 mm over the service area, and it is still better than 1 cm over the non-service area. This paper demonstrates that the BeiDou regional constellation could provide global service to LEO satellites for the POD and the PROD. Finally, the benefit of geostationary earth orbit（GEO） satellites is illustrated for POD.
基金Funded by the Natural Science Foundation of China (No. 40504002)the 973 Program of China (No. 2006CB701301).
文摘Based on the orbit integration and orbit fitting method, the influence of the characters of the gravity model, with different precisions, on the movement of low Earth orbit satellites was studied. The way and the effect of absorbing the influence of gravity model error on CHAMP and GRACE satellite orbits, using linear and periodical empirical acceleration models and the so-called "pseudo-stochastic pulses" model, were also analyzed.
文摘Satellite-to-Satellite tricking (SST) data can be used to determine the orbits of spacecraft in two ways. One is combined orbit determination, which combines SST data with ground-based tracking data and exploits the enhanced tracking geometry. The other is the autonomous orbit determination, which uses only SST. The latter only fits some particular circumstances since it suffers the rank defect problem in other circumstances. The proof of this statement is presented. The nature of the problem is also investigated in order to find an effective solution. Several. methods of solution are discussed. The feasibility of the methods is demonstrated by their application to a simulation.
文摘Autonomous navigation of navigation satellite is discussed. The method of auto-orbit determination using the erosslink range and orientation parameters constraining is put forward. On the basis of the analysis of its feasibility, some useful conclusions are given.
基金supported by the Collaborative Precision Positioning Project funded by the Ministry of Science and Technology of China (No.2016YFB0501900)China Natural Science Funds (No.41231064,41674022,41574015)
文摘For the two newly launched satellites(PRN number 27 and 28) of the future global BeiDou navigation satellite system(BDS-3), there is no available broadcast ephemeris data and other initial orbit information, but the initial orbit is the fundamental of the comprehensive analysis of the satellites and their signals. Precise orbit determination(POD) also requires determination of a priori initial value with a certain precision in order to avoid problems such as filter divergence during POD. Compared with the Newton iteration method, which relies on the initial value, this study utilizes the Bancroft algorithm to directly solve the nonlinear equations with the advantage of numerical stability. The initial orbits of these two satellites are calculated based on new code signals, and their results are analyzed and discussed. The experimental results show that, with the exception of very few epochs, when the new code signal is utilized, the median and robust variance factor of the observed residuals computed using pseudo-range observations and the solved initial orbits are less than 4 and 2 m, respectively. It also shows that this solution can be used for rapid initial orbit recovery after maneuvers of the new BeiDou satellites.
