基于位置观测的北斗广播星历速度精度分析

利用广播星历计算导航卫星的速度向量是GNSS高精度实时测速的必要条件.本文分析了仅以卫星位置向量为观测量的北斗广播星历的速度计算精度.从广播星历拟合过程出发,推导了北斗18参数模型的速度向量计算公式.基于北斗13颗在轨卫星一年的实际轨道数据,分析了全年广播星历计算卫星速度向量的精度.结果表明,利用18参数模型计算的速度误差最大在10-4m·-1量级;在相同拟合时段条件下,地球静止轨道（GEO）和倾斜地球同步轨道（IGSO）卫星的速度精度相当,高于中圆地球轨道（MEO）卫星.通过对位置残差序列分析,得出位置残差误差较小且变化趋势平稳是广播星历计算速度精度较高的原因.分析和计算结果验证了仅用位置观测量拟合北斗广播星历算法的有效性.      

卫星位置偏差对双星/GIS组合定位的影响及组合定位试验

建立了双星定位和双星／GIS组合定位的计算模型.在此基础上研究了星历误差对双星定位结果的 影响,卫星位置偏差对组合定位结果的影响,在组合定位中,在经、纬度方向上,分别给卫星位置加入均方差 为1°的位置偏差.结果表明,给卫星位置加入均方差为1°的经、纬度方向位置偏差,并不会降低组合定位的 精度.对这一重要结论,从几何上给出了解释.最后进行了双星／GIS组合定位的物理试验.试验结果表明, 双星／GIS组合定位确实可以大大提高双星定位的精度.     

倾斜同步轨道卫星交叉点位置演化及保持

倾斜同步轨道(IGSO)卫星交叉点位置在各种轨道摄动因素和入轨偏差的共同影响下,其交叉点位置会在赤道上发生漂移,这种漂移将导致服务区域的漂移,从而影响系统的性能,因此IGSO卫星需要进行定期的交叉点位置保持来保证服务的稳定。对IGSO卫星交叉点位置的演化和保持问题进行了深入的研究,首先,讨论了IGSO卫星交叉点位置的演化规律,给出了在地球非球形引力、日月引力摄动以及卫星入轨偏差作用下的交叉点位置的演化公式,并分析了公式的精度;接着,在对交叉点位置演化规律的分析基础上,进一步讨论了IGSO卫星交叉点位置保持策略。     

静止轨道卫星东西位置保持控制参数的优化方法

为提高静止轨道卫星东西位置保持的轨道控制精度,提出一种优化轨道控制参数和推进剂消耗量估算的方法.这种方法考虑卫星轨道机动期间其姿态控制喷气对卫星轨道的摄动力,据此对位置保持参数中的轨控推力器点火时间进行补偿,并根据所有推力器的喷气时间估算推进剂消耗量,以修正轨控策略,实现更加精准的轨道控制效果.以在轨运行的静止轨道卫星的多次向西位置保持的变轨参数优化过程为实例,验证本方法的有效性和正确性,可应用于后续静止轨道卫星东西位置的实际操作中.     

圆轨道星座位置保持控制的仿真与研究

星座构形可由卫星相对地球的位置和卫星相对卫星的位置联合描述。由于星座中卫星所受摄动,星座构形会在运行中发生改变,需要加以控制。针对同一圆轨道内的两颗卫星,推导了描述星座位置运动的数学模型,通过仿真验证了模型的精度;进而利用二次型最优控制理论给出了基于相对位置的控制规律,提出了星座相对位置保持、绝对位置可移动的星座控制方法;论证了相对位置保持控制所需燃料要比绝对位置保持控制所需燃料少,仿真结果验证了方法的有效性。     

静止轨道卫星连续式等间隔脉冲推力东西位置保持策略

参考地球静止轨道卫星传统东西位置保持控制原理,提出连续式等间隔脉冲推力东西位置保持策略.在摄动分析基础上,给出该策略的控制流程以及控制参数计算公式,进行了轨道参数、姿态控制和角动量管理三方面的影响分析,仿真验证了该策略作用下卫星半长轴、偏心率和平经度的变化趋势,结果表明新策略可以长期有效地控制卫星平经度在漂移环范围内,并具有控制流程简单,位保精度高等优势.     

基于能量色散特征计算磁层顶重联位置的方法研究

磁层顶磁场重联是太阳风向磁层输入能量的主要方式.重联如何触发一直是空间物理研究的难点,其机制仍然有待深入研究.由于卫星穿越磁层顶时,很难恰好穿越重联发生的区域,因此难以观测到重联的触发条件.本文利用THEMIS卫星观测,确立了反演磁层顶重联点的方法.当重联刚开始发生时,卫星能够观测到离子的能量色散特征,可利用其计算卫星到重联发生位置的距离.沿着磁力线模型追踪该距离即可反演出磁层顶发生重联的位置.与其他方法进行了对比分析,结果显示本文方法比其他方法具有更高的精度.      

基于电推力器进行南北位置保持的一种地球同步轨道注入参数方法

摘要： 电推力器在静止轨道卫星上应用越来越广泛,特别是基于电推力器进行南北位置保持,可以有效节省推进剂.提出改进的GPS星历参数解析算法,在此基础上考虑包含电推力模型在内多摄动项模型进行地面精密轨道计算,采用微分修正法,提出一种地球同步轨道注入参数方法,该方法可应用于星上自主完成基于电推力器的南北位置保持.仿真算例表明使用该方法得到的轨道注入参数,卫星能够在保证姿态确定精度的同时,完成南北位置保持任务.     

钟差和动力学参数联合约束下的北斗卫星轨道快速确定

受地球非球形引力、第三体摄动和太阳光压等摄动因素的影响,导航卫星位置存在长周期变化趋势,需要定期对导航卫星进行轨道机动,以保持卫星轨位和导航服务区.导航卫星机动后的定轨,特别是GEO卫星,其频繁轨控后的轨道快速确定问题,是制约卫星可用度和导航系统服务性能的重要因素.在基于伪距相位数据的轨道测定中,轨道与钟差的统计相关是制约卫星轨道快速确定的关键因素,特别是在观测弧段短的情况下,待估参数之间的相关性更强,动力学参数估计结果严重失真会导致轨道预报精度衰减明显.当卫星钟差与测站钟差通过外部手段高精度测定后,可以减少待估参数的估计,同时利用长弧定轨的动力学与运动学参数先验信息,对短弧定轨模式进行参数约束,卫星定轨精度将有很大的提升空间.通过钟差与力学参数的联合约束,实现了北斗卫星短弧快速定轨,解决了卫星机动后的轨道快速确定问题,SLR评估的卫星机动后4 h定轨外符视向精度优于0.71 m,比常规方法提高了3倍,预报1 h轨道视向精度为1.89 m,用户等效距离误差(UERE)精度达到1.85 m.     

同步卫星的轨道摄动及定点位置保持

未来大功率通信、广播卫星将采用大面积太阳翼,较强的太阳辐射压会引起轨道偏心率摄动量很大,影响东西向位置精度.本文着重讨论此种情况下的东西向位置保持修正策略;给出一种平经度漂移环方式;导出漂移环参数、最长修正周期及最小修正速度的计算公式并讨论了偏心率修正速度与修正周期的关系.比较了切向速度脉冲两脉冲控制方式及径向切向速度脉冲两脉冲控制方式所需的修正速度及允许的最长修正周期.     

Application of CCD drift-scan photoelectric technique on monitoring GEO satellites

Geosynchronous Earth Orbit (GEO) satellites are widely used because of their unique characteristics of high-orbit and remaining permanently in the same area of the sky. Precise monitoring of GEO satellites can provide a key reference for the judgment of satellite operation status, the capture and identification of targets, and the analysis of collision warning. The observation using ground-based optical telescopes plays an important role in the field of monitoring GEO targets. Different from distant celestial bodies, there is a relative movement between the GEO target and the background reference stars, which makes the conventional observation method limited for long focal length telescopes. CCD drift-scan photoelectric technique is applied on monitoring GEO targets. In the case of parking the telescope, the good round images of the background reference stars and the GEO target at the same sky region can be obtained through the alternating observation of CCD drift-scan mode and CCD stare mode, so as to improve the precision of celestial positioning for the GEO target. Observation experiments of GEO targets were carried out with 1.56-meter telescope of Shanghai Astronomical Observatory. The results show that the application of CCD drift-scan photoelectric technique makes the precision of observing the GEO target reach the level of 0.2″, which gives full play to the advantage of the long focal length of the telescope. The effect of orbit improvement based on multi-pass of observations is obvious and the prediction precision of extrapolating to 72-h is in the order of several arc seconds in azimuth and elevation.     

木星伽利略卫星互掩互食天象的测光技术

对行星天然卫星互掩互食天象进行测光可以得到非常准确的天体测量信息, 这些信息可以用来进一步改进行卫星的运动理论和历表. 对2003年在云南天文台1米望远镜上实际观测到的互掩互食图像进行测光处理, 分析表明, 在对木星去晕和测光孔径为1.5~2.0倍半高全宽条件下, 所得光变曲线有较好的稳定性.      

基于单站CCD漂移扫描光电设备的同步轨道目标测定轨

地基光电观测在同步轨道目标监测领域具有重要作用.为评估单站光电设备对同步轨道目标的实际测定轨能力,利用上海天文台佘山站1.56m望远镜,采用CCD漂移扫描光电技术,对3颗北斗同步卫星开展试验观测,基于卫星精密星历评估目标的测定轨外符精度.结果表明:同步轨道目标的天文定位在方位和俯仰方向上的外符精度均好于0.3";在单圈次观测情况下,尽管轨道预报精度较低,约为数千米量级,但是观测弧段内定轨精度可优于百米;在多圈次观测情况下,轨道改进效果显著,定轨精度优于50m,外推至4d的轨道预报精度为百米量级.此外,定量评估了每晚不同观测时间跨度下同步轨道目标的测定轨精度,为单站光电设备实际应用提供了参考.      

Differential code bias estimation based on uncombined PPP with LEO onboard GPS observations

Differential Code Bias (DCB) is an essential correction that must be provided to the Global Navigation Satellite System (GNSS) users for precise position determination. With the continuous deployment of Low Earth Orbit (LEO) satellites, DCB estimation using observations from GNSS receivers onboard the LEO satellites is drawing increasing interests in order to meet the growing demands on high-quality DCB products from LEO-based applications, such as LEO-based GNSS signal augmentation and space weather research. Previous studies on LEO-based DCB estimation are usually using the geometry-free combination of GNSS observations, and it may suffer from significant leveling errors due to non-zero mean of multipath errors and short-term variations of receiver code and phase biases. In this study, we utilize the uncombined Precise Point Positioning (PPP) model for LEO DCB estimation. The models for uncombined PPP-based LEO DCB estimation are presented and GPS observations acquired from receivers onboard three identical Swarm satellites from February 1 to 28, 2019 are used for the validation. The results show that the average Root Mean Square errors (RMS) of the GPS satellite DCBs estimated with onboard data from each of the three Swarm satellites using the uncombined PPP model are less than 0.18 ns when compared to the GPS satellite DCBs obtained from IGS final daily Global Ionospheric Map (GIM) products. Meanwhile, the corresponding average RMS of GPS satellite DCBs estimated with the conventional geometry-free model are 0.290, 0.210, 0.281 ns, respectively, which are significantly larger than those obtained with the uncombined PPP model. It is also noted that the estimated GPS satellite DCBs by Swarm A and C satellites are highly correlated, likely attributed to their similar orbit type and space environment. On the other hand, the Swarm receiver DCBs estimated with uncombined PPP model, with Standard Deviation (STD) of 0.065, 0.037 and 0.071 ns, are more stable than those obtained from the official Swarm Level 2 products with corresponding STD values of 0.115, 0.101, and 0.109 ns, respectively. The above indicates that high-quality DCB products can be estimated based on uncombined PPP with LEO onboard observations.     

Analysis of total visual and CCD V-broadband observations of Comet C/1995 O1 (Hale-Bopp): 1995–2001

We developed a general method for determination of water production rates from groundbased visual observations and applied it to Comet Hale–Bopp. Our main objective is to extend the method to include total visual magnitude observations obtained with CCD detector and V filter in the analysis of total visual magnitudes. We compare the CCD V-broadband careful observations of Liller [Liller, W. Pre-perihelion CCD photometry of Comet 1995 O1 (Hale-Bopp). Planet. Space Sci. 45, 1505–1513, 1997; Liller, W. CCD photometry of Comet C/1995 O1 (Hale-Bopp): 1995–2000. Int. Comet Quart. 23(3), 93–97, 2001] with the total visual magnitude observations from experienced international observers found in the International Comet Quarterly (ICQ) archive. A data set of ∼400 CCD observations covering about the same 6 years time span of the ∼12,000 ICQ selected total visual magnitude observations were used in the analysis. A least-square method applied to the water production rates, yields power laws as a function of the heliocentric distances for the pre- and post-perihelion phases. The average dimension of the nucleus as well as its effective active area is determined and compared with values published in the literature.     

Long-term photometric signature study of two GEO satellites

Geostationary earth orbit satellites have been extensively used for unique high-orbit stationary characteristics. Long-term precise investigation is an important issue in the observation of GEO satellites, since it can provide valuable information on the satellites’ operation state, discrimination and early warning analysis. Ground-based optical-electronic devices play a significant role in the observation. 4-month photometric signature variation of two satellites is presented based on the successive observations using the 1.56-meter telescope of Shanghai Astronomical Observatory (SHAO). It can be concluded that the long-term brightness change mainly results from the sun declination angle and regular orbit maneuver. Moreover, the solar panel offsets of the two satellites are analyzed and found to be approximately 4 degrees. Estimation of photometric accuracy reaches 0.15 mag for the application of CCD drift-scan optical-electronic technique.     

Operation and Application of A, B Satellites for Environment and Disaster Monitoring and Forecasting

Environment and disaster monitoring and forecasting small satellite constellation A and B satellites (HJ-1-A, B) are called "environment and disaster reduction satellites A and B' for short. The constellation adopts a 10:30 LT sun-synchronous circular orbit, with orbit altitude of 649 km. HJ-1-A and HJ-1-B are distributed with a phase difference of 180o in the same orbital plane, so as to enhance the time resolution of earth observation. The satellites have orbit maintenance capability, the lifetime is 3 years. Both satellites adopt CAST968 platforms. Two wide-coverage multispectral CCD cameras with resolution 30 m and width 700 km, a super-spectral imager with resolution 100 m and width 50 km as well as a data transmission subsystem of 120 Mbit/s are deployed on HJ-1-A, which also carries Ka communication testing equipment of Thailand. HJ-1-B has two wide-coverage multispectral CCD cameras (the same as satellite A), one infrared camera with resolution 150 m and width 720 km and a data transmission subsystem of 60 Mbit/s. The coverage period of the wide-coverage multispectral CCD camera is 48 hours. The revisit period of super-spectral imager is 96 hours and the coverage period of infrared camera is 96 hours.      

CCD 遥感图像几何校正的实时处理算法初探

为向星上处理迈进,实时完成ＣＣＤ遥感图像１Ｂ级几何校正,迫切需要建立一个准确而又能快速进行运算的反映校正前后几何位置对应关系的函数表达式和重采样算法。通过论证把校正前后一维位置对应关系高度精确地表示为一元三次多项式关系,从而改用递推极大地简化计算;在此基础上提出了一维灰度重采样新算法——加权平均法及其快速算法。当整幅图每行均按此处理之后,再根据获得的ＧＰＳ定轨定姿数据加注出图像的地理经纬网,从而得到完整的１Ｂ级几何校正产品。初步微机仿真实验表明,只需适当引入多处理器,本算法达到实时完全可行。     

Stereoscopic observations from meteorological satellites

The capability of making stereoscopic observations of clouds from meteorological satellites is a new basic analysis tool with a broad spectrum of applications. Stereoscopic observations from satellites were first made using the early vidicon tube weather satellites (e.g., Ondrejka and Conover [1]). However, the only high quality meteorological stereoscopy from low orbit has been done from Apollo and Skylab, (e.g., Shenk $\text{et al.}$ [2] and Black [3], [4]). Stereoscopy from geosynchronous satellites was proposed by Shenk [5] and Bristor and Pichel [6] in 1974 which allowed Minzner $\text{et al.}$ [7] to demonstrate the first quantitative cloud height analysis. In 1978 Bryson [8] and desJardins [9] independently developed digital processing techniques to remap stereo images which made possible precision height measurement and spectacular display of stereograms (Hasler $\text{et al.}$ [10], and Hasler [11]). In 1980 the Japanese Geosynchronous Satellite (GMS) and the U.S. GOES-West satellite were synchronized to obtain stereo over the central Pacific as described by Fujita and Dodge [12] and in this paper. Recently the authors have remapped images from a Low Earth Orbiter (LEO) to the coordinate system of a Geosynchronous Earth Orbiter (GEO) and obtained stereoscopic cloud height measurements which promise to have quality comparable to previous all GEO stereo. It has also been determined that the north-south imaging scan rate of some GEOs can be slowed or reversed. Therefore the feasibility of obtaining stereoscopic observations world wide from combinations of operational GEO and LEO satellites has been demonstrated.Stereoscopy from satellites has many advantages over infrared techniques for the observation of cloud structure because it depends only on basic geometric relationships. Digital remapping of GEO and LEO satellite images is imperative for precision stereo height measurement and high quality displays because of the curvature of the earth and the large angular separation of the two satellites. A general solution for accurate height computation depends on precise navigation of the two satellites. Validation of the geosynchronous satellite stereo using high altitude mountain lakes and vertically pointing aircraft lidar leads to a height accuracy estimate of ± 500 m for typical clouds which have been studied. Applications of the satellite stereo include: 1) cloud top and base height measurements, 2) cloud-wind height assignment, 3) vertical motion estimates for convective clouds (Mack $\text{et al.}$ [13], [14]), 4) temperature vs. height measurements when stereo is used together with infrared observations and 5) cloud emissivity measurements when stereo, infrared and temperature sounding are used together (see Szejwach $\text{et al.}$ [15]).When true satellite stereo image pairs are not available, synthetic stereo may be generated. The combination of multispectral satellite data using computer produced stereo image pairs is a dramatic example of synthetic stereoscopic display. The classic case uses the combination of infrared and visible data as first demonstrated by Pichel $\text{et al.}$ [16]. Hasler $\text{et at.}$ [17], Mosher and Young [18] and Lorenz [19], have expanded this concept to display many channels of data from various radiometers as well as real and simulated data fields.A future system of stereoscopic satellites would be comprised of both low orbiters (as suggested by Lorenz and Schmidt [20], [19]) and a global system of geosynchronous satellites. The low earth orbiters would provide stereo coverage day and night and include the poles. An optimum global system of stereoscopic geosynchronous satellites would require international standarization of scan rate and direction, and scan times (synchronization) and resolution of at least 1 km in all imaging channels. A stereoscopic satellite system as suggested here would make an extremely important contribution to the understanding and prediction of the atmosphere.     

Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK

The overlapping carrier frequencies L1/E1, L5/E5a and B2/E5b from GPS/Galileo/BDS allow inter-system double-differencing of observations, which shows a clear advantage over differencing of the observations of each constellation independently. However, the inter-system biases destroy the integer nature of the inter-system double-differencing ambiguities. Two methods of direct rounding and parameter estimation are used to determine the ISB value. By analyzing data collected from Curtin University from 2015 to 2018, the phase and code inter-system bias (ISB) are related to the receiver type, firmware version and the selected overlapping frequency. Upgrade of receiver firmware version results in changes of ISB values. For example, the upgrade of Javad firmware in Dec, 15, 2017 causes the difference of 0.5 cycles ISB between BDS GEO and non-GEO satellites. By comparing the three dynamic models which include white noise process, random walk process, and random constant in the parameter estimation method, the ISB determined by the random constant model is consistent with the value obtained by the direct rounding method. After the calibration of ISBs, the performances of tightly combined positioning are assessed. The success rate of ambiguity resolution and accuracy of positioning for the tight combination (TC) are significantly improved in comparison with that for the loose combination (LC) over short baselines. For L5/E5a, on which only few satellites can be observed, the maximum increase in success rates of ambiguity resolution can reach 31.7%, i.e., from 54.9% of LC to larger than 86.6% of TC, and the positioning accuracies can even be increased by 0.13 m, i.e., from 0.208 of LC to 0.074 m of TC in East direction for the mix-receiver TRIMBLE NETR9-SEPT POLARX4 in 2018.