基于误差发散规律的低轨卫星大气阻力系数计算方法

低轨卫星轨道预报精度受到大气模型和大气阻力系数精度的制约,给一些高精度的空间和航天任务带来不利影响.提出了一种基于沿迹方向误差发散规律的大气阻力系数计算新方法.首先通过理论推导给出低轨卫星轨道预报中沿迹误差发散的分析表达式,定量描述初值误差和模型误差对沿迹误差的综合影响;提出利用定轨段的基本信息,优选预报段所采用的阻力系数,抑制沿迹误差的发散速率,从而降低沿迹方向预报误差的最大值,提高短期预报精度.以400 km附近的GRACE-A卫星的全弧段星载GPS高精度资料为基础,检验了方法的精度和成功率.结果表明:相对于传统的定轨预报方法,新方法能提高24 h短期预报精度约45%,成功率约71%,总体有效率约86%;方法对低、中、高等3种太阳辐射水平均有效,对于中低等级的地磁扰动也有效,具备较好的应用价值.     

天文大气折射的较差测量方法及试观测结果

受到大气折射的影响,天文观测上通常回避仰角15°以下的目标的观测,但作为大气折射的完整理论研究,低仰角下的大气折射仍然是值得分析探究的.特别是对某些工程应用方面,低仰角的目标有时必须要观测.提出了一套新的利用较差方法测定大气折射的思路.利用一台较大视场的望远镜从天顶开始,在不同高度上对星空作一系列观测,计算不同天顶距处大气折射函数的各阶导数,最后经数值积分可给出大气折射实测值.该方法不依赖于严格的地方参数和复杂精密的观测仪器,并且观测原理相对简单. 2007年底,利用一台简易的大视场望远镜在兴隆观测站进行了试验观测,根据较差方法实测得到真天顶距44.8°至87.5°的大气折射值,初步证明了大气折射较差测量方法的可行性.受到观测条件的限制,本次实测结果精度有限,偶然误差最大约为6",并且存在一定的系统差.在天顶距84°时,与普尔科沃大气折射表的差值约为15".如何消除因积分模型误差引入的累积误差是今后需要解决的关键问题.     

基于单站短弧段光学观测的低轨卫星轨道预报算法

基于单站短弧段的目标跟踪预报技术是保证高精度光电经纬仪在非常规环境下正常跟踪捕获目标的重要途径.构建了基于非线性滤波器的跟踪预报算法,能够在正常情况下为闭环跟踪提供引导数据,同时构建了基于非线性变换的目标预报算法,可以在无有效观测数据情况下为经纬仪提供轨迹引导,保证目标不丢失.证明了非线性滤波在单站短弧段跟踪预报算法中比常规扩展卡尔曼滤波(Extended Kalman Filter,EKF)更有效.计算表明,本文设计的非线性滤波器可以作为光学跟踪设备的引导算法,引导精度同经纬仪的随机测量精度等量级,在设备系统误差达到50"时,预报60 s的精度可达到20",仍然满足跟踪设备的视场要求.     

由极迹半径的外推预报地极坐标

本文根据极迹半径的变化规律提出了一种通过外推极迹半径预报地极坐标的方法。本方法的特点是预报时地极坐标的x和y两个分量结合成整体运算。一般情况下半年的预报度在80％以上,但在极迹半径最小的时期预测误差太大。     

卫星与空间碎片碰撞预警的快速算法

提出了一种空间目标与空间碎片碰撞预警的快速计算方法，该算法不仅计算速度快，而且减缓了沿迹误差的影响，延长了预警的有效期，利用该方法，检测一个卫星与5800个空间碎片在一天内是否碰撞，计算时间小于3秒；用2天前的资料，预报计算法国CERISE卫星和空间碎片的碰撞时间误差约为0．02秒，空间位置误差约为500—600米．     

区域网GPS卫星独立定轨的研究

根据我国广域差分GPS系统对GPS卫星定轨的要求,通过对国内GPS跟踪网实测数据的处理,分析和讨论了区域网定轨的数据处理方法和可能达到的定轨精度。为了提高所定轨道的稳定性和先进实时预报的精度,通过对所定轨道的误差分析提出在轨道沿迹方向引入经验加速度计算方案。计算结果表明,采用此方法后GPS卫星的定轨精度有了显提高。既论证了利用我国区域GPS网和广播星历进行独立定轨的可行性,也阐述了提高轨道预报精度的方法。     

德令哈13.7m望远镜重力变形研究

借助于数字摄影测量结果调整天线面板,使德令哈13.7 m望远镜在仰角52°时获得最佳反射面面形,从而使天线效率在观测仰角范围内得到整体优化.与之前基于经纬仪测量的面板调整结果相比,天线口径效率提高约1倍.依据不同俯仰姿态下的测量结果,得到了天线的重力变形模型,包括副面偏移和倾斜、主面焦距和面形偏差随仰角变化的规律.根据不同仰角的面形偏差测量数据反演反射面重力变形模型时,采用了数据拟合方法,这样可以减小测量误差对模型精度的影响.     

辛积分器中沿迹误差的一种补偿方法

辛积分器严格描述了一摄动Ｈａｍｉｌｔｏｎ系统的流，因而导致天体轨道的沿迹误差随时间呈线性增长趋势。本文利用这一特点，提出了一种对其沿迹误差进行估算的数值方法，从而达到了对数值结果进行沿迹误差补偿的目的，数值结果证实了此方法在较大积分步长和较长积分时间的数值计算中是有效的。     

流星雷达测距精度的改进

针对目前流星雷达测距误差大的问题，本文提出了提高流星雷达测距精度的新方法，即提高采样速率，用相关分析确定回波脉冲参考点的方法。该方法使流星雷达的测距精度提高一个数量级，测距误差降到±１４ｍ，使流星雷达不仅可以用来观测研究流星，还可用于监测飞机、火箭的飞行等，扩大流星雷达的应用。     

卫星定位精度对在轨甚低频干涉测量影响分析

由于地球电离层的阻挡以及其他干扰,在地面难以进行有效的甚低频天文观测,而使用搭载于绕地或绕月轨道卫星上的甚低频天线进行干涉观测会大大提高观测灵敏度和角分辨率.卫星定位精度会影响观测数据的处理结果,进而影响成图质量,并且会大幅影响飞行项目的复杂度和总成本.首先分析卫星姿态控制精度和星间基线测量精度对绕地轨道甚低频干涉观测的影响,之后对干涉观测中信号延时误差进行仿真,研究延时误差与数据相关处理中条纹搜索范围之间的关系,并对数据计算速度需求进行估算.分析和仿真结果可以辅助相关空间甚低频观测项目总体方案的制定和工程指标的优化.     

A Constant Elevation Search Method for Narrow Beam Radar to Capture Space Targets

Position prediction is the basis of guiding radar to seek, capture, and track space targets, although there’s a risk of insufficient accuracy for narrow beam radar due to the influence of orbit determination and propagation error. Since the along-track error is the main factor that affects the accuracy of position prediction, this paper proposed a Constant Elevation Search (CES) method for narrow beam radar, based on the idea of along-track error compensation, to seek and capture transit targets. With along-track error estimation, all possible positions on a specified elevation can be observed, thus the success rate of capturing transit targets will be improved. Simulation shows that when the position prediction method fails, the CES method still works.     

An Orbital Prediction Algorithm for LEO Satellites Based on Optical Observations of Short Arcs at Single Station

The forecasting technique of the target tracking based on the short arcs at single station is an important way to guarantee that high-precision photoelectric theodolites can normally track and capture the targets in unconventional environments. We construct the tracking prediction algorithm based on nonlinear filter, which can provide the guiding data for the closed loop tracking under normal circumstances. At the same time we also construct the target prediction algorithm based on the nonlinear transformation, without valid observational data, which can provide a track guidance for the theodolite and ensure that the targets will not be lost. It is demonstrated that the nonlinear filtering is more effective than the EKF (extended Kalman filter) in the tracking prediction algorithm of the short arcs at single station. The results indicate that the nonlinear filter designed in this paper can be used as the guiding algorithm for the optical tracking equipments. And its guiding accuracy is in the same order of magnitude of the theodolite's random measurement accuracy. When the systematic error of the equipments reaches 50″, the accuracy can achieve 20″ for predictions in 60 s. This still satisfies the requirement of the field of view of the tracking equipments.     

基于历史TLE的空间目标轨道预报误差演化规律研究

北美防空司令部(North American Aerospace Defense Command, NORAD)发布的双行根数(Two Line Element, TLE)是广大航天工作者最常用的轨道根数,与其对应的轨道模型是SGP4/SDP4 (Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4)解析模型.由于TLE中并没有包含相应的轨道精度信息,编目轨道的应用范围受到很大的限制.基于Space-Track网站发布的历史TLE数据和配套的SGP4/SDP4动力学模型,采用定轨标预报的方法统计并生成了大量目标轨道的预报误差,通过对预报轨道的时间区间划分给出了每个目标的预报误差随预报时间变化的拟合系数,并进一步对不同类型轨道预报误差的演化规律和特征进行了分类讨论,给出了4种轨道类型目标的轨道预报误差随时间演化的平均解析模型,为拓展双行根数的应用提供有价值的参考.     

基于光学辅助指向测量的太赫兹望远镜指向误差模型研究

针对太赫兹波段天文点源目标较少, 指向测量相对困难的特点, 研究了利用与太赫兹天线共轴的小型光学望远镜来辅助太赫兹望远镜指向测量以及建立指向误差修正模型的方法. 依托紫金山天文台1.2 m斜轴式太赫兹天线开展了光学辅助指向测量的实验研究, 利用一台安装在天线背架上的100mm口径折射式光学望远镜获得了优于2$''$的指向测量精度. 此外, 通过对斜轴天线的结构分析以及大气折射和本地恒星时(Local Sidereal Time, LST)偏差等误差来源的分析, 建立了包含23个误差项的斜轴式光学指向修正模型, 实现了约3$''$的拟合精度. 最后, 借助高精度数字摄影测量对光电轴一致性进行了标定, 并针对其对指向模型精度的影响进行了讨论. 研究成果将为南极5 m太赫兹望远镜(The 5m Dome A Terahertz Explorer, DATE5)及其他太赫兹望远镜提供指向测量和指向修正模型方面的技术参考.     

日月轨道计算对TLE深空编目目标预报精度的影响

双行根数(Two Line Element, TLE)作为一类广泛使用的空间物体编目数据, 其预报精度和误差特性是TLE编目 在空间碎片研究中所要关注的问题之一. TLE编目需要配合SGP4/SDP4 (Simplified General Perturbations 4/Simplified Deep Space 4)模型进行轨道预报, 对深空物体来说, 主要考虑带谐项$J_2$、$J_3$、$J_4$摄动、 第三体日月摄动和特殊轨道共振问题修正等. 其中, SGP4/SDP4模型第三体摄动计算时, 对日月轨道近似采用了长期进动根数和 简单平运动的方式, 在外推10d时存在约2$^\circ$--3${^\circ     

The meteor radar response function: Theory and application to narrow beam MST radar

The meteor radar response function is an important tool for analyzing meteor backscatter observed by radar systems. We extend previous work on the development of the response function to include a non-uniform meteor ionization profile, provided by meteor ablation theory, in contrast to what has been assumed in the past. This has the advantage that the height distribution of meteors expected to be observed by a radar meteor system may be accurately modeled. Such modeling leads to meteor height distributions that have implications for the composition of those meteoroids ablating at high altitudes which may be observed by “non-traditional” meteor radars operating at MF/HF. The response function is then employed to investigate meteor backscatter observed by narrow beam MST radars which in recent years have been used increasingly to observe meteors.     

Combination Scattering by Anisotropic Langmuir Turbulence with Application to Solar Radar Experiments

We develop a theory for radar signal scattering by anisotropic Langmuir turbulence in the solar corona due to a t+l ⇄ t process. Langmuir turbulence is assumed to be generated within a cone by a narrow type III burst electron beam. Using wave-kinetic theory we obtain expressions for the frequency shift, scattering cross-section of the turbulence, coefficient of absorption (due to scattering) and optical depth. On the basis of those expressions we give some estimates for an echo spectrum. We show that the minimum radar echo frequency shift is determined by the minimal phase velocity of the Langmuir waves, the maximum shift is determined by the electron beam velocity, but in any case it can not exceed −wt/2 (decay) and wt (coalescence), where wt is the frequency of a radar signal. The angular characteristics of the scattered signal differ dramatically for the cases of coalescence and decay. The signal is scattered into a narrow cone high above the specular reflection point (wp ≪ wt), but in the vicinity of wp ∼ wt/2 the red-shifted echo is scattered isotropically, while the blue-shifted echo is scattered into a even narrower cone. We show that absorption (due to scattering) increases with increasing radar frequency. The dependence of the absorption on the local plasma frequency is strongly determined by the Langmuir turbulence spectrum. Our theory shows that the role of the nonlinear scattering process t+l ⇄ t is essential and that such process can be used for radar studies of the spectral energy density of anisotropic Langmuir turbulence.     

An Improved Grey Model and Its Application Research on the Prediction of Real-time GPS Satellite Clock Errors

In the work on the real-time GPS precise point positioning, the realtime and reliable prediction of the satellite clock error is one of the keys to the realization of the GPS real-time high accuracy point positioning. The satelliteborne GPS atomic clock has high frequency, is very sensitive and extremely easy to be influenced by the outside world and its own factors. Therefore, it is very difficult for one to know well its complicated and detailed law of change, with these attributes being in accordance with the characteristics of the theory of grey system. Thus, it is considered that the process of variation of the clock error is regarded as a grey system. On the basis of the exploration of the limitations of the quadratic polynomial and grey model satellite clock error predictions, the research on the real-time prediction of the GPS satellite clock error by taking advantage of the improved grey model is proposed. Finally, the materials of the GPS satellite clock error of 3 different time intervals are used to make the accuracy analysis of the clock error prediction of different sampling intervals, to study the relation between the grey model exponential coefficient and the prediction accuracy and to make the analysis of the comparison of the prediction accuracy with that of the quadratic polynomial method. The general relation between the different types of satellite clock errors and the model exponential coefficients is summarized and compared with the IGS final clock error ephemeris product to test and verify the feasibility and availability of the improved prediction model proposed in the present article so as to provide the higher-accuracy satellite clock error products for the real-time GPS dynamic precise point positioning.     

Modeling and calibration of pointing errors with alt-az telescope

This paper presents a new model for improving the pointing accuracy of a telescope. The Denavit–Hartenberg (D–H) convention was used to perform an error analysis of the telescope's kinematics. A kinematic model was used to relate pointing errors to mechanical errors and the parameters of the kinematic model were estimated with a statistical model fit using data from two large astronomical telescopes. The model illustrates the geometric errors caused by imprecision in manufacturing and assembly processes and their effects on the pointing accuracy of the telescope. A kinematic model relates pointing error to axis position when certain geometric errors are assumed to be present in a telescope. In the parameter estimation portion, the semi-parametric regression model was introduced to compensate for remaining nonlinear errors. The experimental results indicate that the proposed semi-parametric regression model eliminates both geometric and nonlinear errors, and that the telescope's pointing accuracy significantly improves after this calibration.     

Determination of stellar positions with a cross-correlation procedure

A technique of cross-correlation for the determination of stellar positions is proposed. It is a technique suitable for the recording of stellar positions on transit instruments, photoelectric astrolabes, and also on the meridian circles with reversed observing method. An advantage of this process is that there is no hypothesis required for the modelling image profile, and for the noise. It is closer to unbiased estimation. Results by this new method on a transit instrument are presented. For the brighter stars, there is no great difference in accuracy between the median position and the cross-correlation value according to the results.