双星定位系统/SINS深组合导航系统研究

以卫星模拟器为基础，介绍了一种双星定位系统定姿的基本原理，解决了双星定位系统存在的位置滞后即定位实时性较差的缺陷，把双星定位系统的位置与姿态信息和捷联惯导系统的位置、姿态信息进行组合后，可有效地提高系统的精度．     

Tetrahedral satellite formation: Geomagnetic measurements exchange and interpolation

This note presents a study of a four-satellite tetrahedral formation to collect, process, and exchange multipoint measurements of geomagnetic field in a near-polar orbit. The study is conducted as a series of numerical experiments based on simulated spacecraft orbits and corresponding geomagnetic field models output. The four satellites are assumed to move in near-circular orbits specifically chosen to maintain the tetrahedron quality. The satellites exchange their simulated magnetometers readings and the collected multipoint measurements are processed on board of any of them thus creating an instantaneous interpolated map of the geomagnetic field in the interior of the tetrahedron. Interpolation is carried out with the use of Kriging algorithms, known in geostatistics for capturing spatial correlation of the data and taking into account statistical properties of the interpolated variables. We propose a concept of a servicing formation, and analyze interpolation accuracy for different formation sizes. It is then discussed how the processed multipoint measurements can be provided as a service to other nearby satellites. Finally, we show that using the existing COTS magnetometers it is possible to obtain real-time interpolation data, which are more precise at a given point and time than a conventional onboard magnetic field model, thus ensuring better attitude determination routines performance in the serviced spacecraft.     

Higher order ionospheric error correction in BDS precise orbit determination

The ionospheric error affects the accuracy of the Global Navigation Satellite Systems observation and precise orbit determination. Usually, only the first order ionospheric error is considered, which can be eliminated by the ionospheric-free linear combination observation. But the remaining higher order ionospheric error will affect the accuracy of observations and their applications. In this paper, the influence of the higher order ionospheric error have been studied by using the International Geomagnetic Reference Field 13 and the Global Ionosphere Maps model produced by the Center for Orbit Determination in Europe. Focus on ionospheric error, the experiment of paper at doy 302 in 2019, which show that the second order ionospheric error impacting BeiDou Navigation Satellite System (BDS) B1I and B3I observation is 6.3569 mm and 11.8484 mm, respectively. Whereas, the third order ionospheric error impacting BDS B1I and B3I observation is 0.1734 mm and 0.3977 mm, respectively. Due to the current measurement accuracy of BDS carrier-phase observation can reach 2 mm, the influence of high order ionospheric error on observation should be considered. For BDS precise orbit determination, the orbit overlapping results are indicated that its orbit accuracy can be improved approximately 5 mm with the higher order ionospheric error correction, which is also in agreement with the results of Satellite Laser Ranging in this work.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

惯性/偏振光组合定姿技术研究

惯性导航系统是一种全参数、全天候的自主导航系统,但其姿态误差随时间发散。太阳光在大气中传播时在天空具有稳定的偏振模式,惯性/偏振光组合定姿系统就是利用偏振光中的方位信息实现惯性导航系统姿态角误差的估计。研究了惯性/偏振光组合定姿系统模型,利用跑车试验数据对所设计的组合定姿系统进行了仿真。结果表明:惯性导航系统引入偏振光导航信息后,可以有效地抑制惯性导航姿态角误差发散的趋势,提高姿态精度,满足中等精度姿态参考系统的需求。     

Nonlinear filtering methods for spacecraft navigation based on differential algebra

The paper investigates the problem of nonlinear filtering applied to spacecraft navigation. Differential algebraic (DA) techniques are proposed as a valuable tool to implement the higher-order numerical and analytic extended Kalman filters. Working in the DA framework allows us to consistently reduce the required computational effort without losing accuracy. The performance of the proposed filters is assessed on different orbit determination problems with realistic orbit uncertainties. The case of nonlinear measurements is also considered. Numerical simulations show the good performance of the filter in case of both complex dynamics and highly nonlinear measurement problems.     

快速凝固制备的Cu—Sn—Ag—In钎料薄带的韧脆性探讨

用快速凝固法将代银钎料Ｃｕ－Ｓｎ－Ａｇ－Ｉｎ合金制备成薄带，透射电镜分析结果表明薄带组织为晶态。扫描电镜分析结果显示薄带显微组织由银或富银的面心立方固溶体ａ相和脆性相组成。     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Impact of DORIS ground antennas environment on their radio signal quality

Among the factors which may disrupt the DORIS measurements quality, the ground antennas environment is of high importance. For a set of 15 selected DORIS beacon, the differences between the effective and theoretical power received on-board the satellites (SPOT-5 and Envisat) have been analyzed in terms of spatial direction around the antenna. Such antenna maps have also been established regarding the Doppler residuals of the least-square precise orbit adjustment. Thanks to 360° views from the antennas and aerial views of the sites, the impact of the signal obstructions (trees, roofs, antennas …) on power attenuation and Doppler residuals is discussed. Depending on the nature of the obstructed object, the attenuation level can reach more than 5 dB, and the residual RMS of the orbit adjustment may be doubled from the nominal value, reaching 1 mm/s locally. The nature of the ground at the foot of the antennas has been correlated to DORIS signal quality at high elevation: reflections on flat surfaces (e.g. roofs) affect the signal more significantly than reflections on natural ground (e.g. soil). In particular, a modeling of the multipath phenomenon affecting Fairbanks site has been established and fits remarkably with the observations. Finally, an evaluation of the direct impact of obstructing objects on the orbit has also been performed. The example of a scaffolding at Kauai site displays a few millimeters error in the along-track position of the satellite.     

New improved orbit solutions for the ERS-1 and ERS-2 satellites

Improved orbit solutions of the European Remote Sensing Satellites ERS-1 and ERS-2 have been computed in the ITRF2005 terrestrial reference frame using the recent models based mainly on IERS Conventions 2003. These solutions cover the periods 3 August 1991 to 8 July 1996 for ERS-1, and 3 May 1995 to 4 July 2003 for ERS-2. For each satellite, the final orbit solution is based on a combination of three separate orbit solutions independently computed at the Delft Institute of Earth Observation and Space Systems (DEOS) of the Delft University of Technology (The Netherlands), the Navigation Support Office of the European Space Operations Centre (ESOC, Germany) and the Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences (Germany) using three different software packages for precise orbit determination, but using the same models in the same terrestrial reference frame within the European Space Agency (ESA) project ‘Reprocessing of Altimeter Products for ERS (REAPER)’. Validation using radar altimeter data indicates that the new combined orbits of ERS-1 and ERS-2 computed by us are significantly more accurate, approaching the 2–3 cm level in radial direction, than previously available orbit solutions.