一种适用于长弧段的初轨计算方法

根据单位矢量法测轨原理，在人造卫星初轨计算的单位矢量法基础上，给出了一种适用于长弧段的初轨计算方法．该方法既适用于长弧段，也适用于短弧段；适用于各种不同类型的观测资料和任意偏心率、任意轨道倾角的人造地球卫星；有利于提高初轨测定精度；并改善整个计算收敛性．特别需要指出的是，该方法与有摄初轨计算的单位矢量法相结合，为单位矢量法从初轨计算推广到轨道改进打下了坚实的基础．     

单站单圈测距资料初轨计算的单位矢量法

根据单位矢量法测轨原理，在已知轨道面参数i，Ω的前提下，给出了一种适用于单站单圈测距资料的初轨计算方法．计算结果表明，方法是有效的．在没有测角资料或测角资料精度显著低于测距精度的情况下，有应用价值．     

初轨计算方法中的常数系统差修正

在单位矢量法的基础上给出了一种常数系统误差修正方法．大量的数据验算结果表明,当测量数据中含有显著常数系统误差时,在初轨计算中修正常数系统差是有效的,可以提高测轨精度．     

论初轨计算的最佳精度及二重解

从最优估计的角度,研究短弧测角资料初轨计算的最佳精度及准最佳精度算法．按照统计理论,对一般的参数估计问题,其估计误差的方差存在Cram’e4-Rao下界,这就是初轨理论的最佳精度;通过大量的模拟计算表明,文［1］中提出参考矢量法,其精度已充分接近C·R下界,且模拟计算轨道统计精度与观测误差成正比,可以认为是一种初轨计算的准最佳精度算法;对于小偏心率情况,发现一个重要的事实：初轨计算有二重解,更准确地说,在M0的左右两个半平面内有真假两个解,偏心率愈小,平近点角M0愈接近±π／2,给出的假解的概率愈大．针对这一情况,提出一种算法,可以给出2个解,其中必定有一个是真解,它有较高的精度.     

移动站的运动对轨道确定影响的分析

介绍了测站移动速度对四元素短弧初轨确定的影响．从测站速度矢量的影响因素着手，简化地对移动站的移动速度与影响因素之间的关系进行了理论推导．以单位矢量方法中的3种方法为例，修改了轨道的条件方程，并做了大量的数据验算，最后给出了两圈外测数据的两种不同情况的单圈定轨结果，就此初步分析了移动站的移动速度对定轨结果的影响及规律．这种在定轨方法中补充考虑移动站的移动速度影响的方法，计算模型简洁，容易理解，便于运用．     

数值积分和数据融合在单位矢量法中的应用

单位矢量法是紫金山天文台提出并得到广泛应用的一种系列轨道确定算法.该算法通过投影变换,对不同类型的观测数据构成不同的条件方程,从而便于不同类型数据的加权处理,充分发挥了高精度数据在定轨中的作用,显著提高了定轨精度.     

一种有根数的约束的单位矢量法

本文在单位矢量法测轨原理的基础上,给出了一种根据已知的部分粗略轨道根数进行约束求解的实用测轨方法。该方法已经在实际工作中得到应用。大量的实测计算表明,该方法在观测资料较差时,能够较好地提高初轨计算收敛率,测轨精度也优于一般比典的初轨计算方法。     

关于改进的Encke方法中参考轨道的选择和校正

对于改进的Ｅｎｃｋｅ方法，选择适当的参考轨道是一个关键．然而，对于人造地球卫星长弧轨道计算，目前所给出的几种参考轨道均需要逐段校正，这将给定轨问题带来附加的复杂性．本文将仔细探讨如何选择参考轨道和减少校正次数．     

初轨计算中的病态分析

本文对现有初轨计算方法进行病态性分析与误差分析；研究结果表明：病态对现有初轨算法的影响，主要来源于法方程系数中包含观测误差．系数行列式愈大，定轨精度的损失愈多，当■被随机误差项△μ淹盖时，现有初轨算法将失效．此外，仿真结果还显示：■与△μ的大小还极大地依赖观测弧段的空间位置，当观测弧段包含近站点作为中点时，■最大，而■小，此时定轨精度较高；当观测弧段位于近站点的某一侧时，■小，而■大，此时定轨精度较低，观测弧段愈偏离近站点，病态影响愈大；因而在观测时，应尽量使观测弧段与近站点对称（此时μ值较大），这是提高短弧定轨的一种有效途径．     

人卫长弧定轨中的摄动计算问题

人卫精度定轨中摄动星历表的计算，常采用数值方法，对于长弧定轨，将会遇到数值稳定问题，有人提出采用Ｅｎｃｋｅ特别摄动法，但要保持数值稳定，关键在于坐标摄动法数值计算中参考轨道的选取，而不是用不用Ｅｎｃｋｅ变换，更有效的方法是一般数值方法中采用稳定化措施。     

A Hamiltonian approach to the planar optimization of mid-course corrections

Lawden’s primer vector theory gives a set of necessary conditions that characterize the optimality of a transfer orbit, defined accordingly to the possibility of adding mid-course corrections. In this paper a novel approach is proposed where, through a polar coordinates transformation, the primer vector components decouple. Furthermore, the case when transfer, departure and arrival orbits are coplanar is analyzed using a Hamiltonian approach. This procedure leads to approximate analytic solutions for the in-plane components of the primer vector. Moreover, the solution for the circular transfer case is proven to be the Hill’s solution. The novel procedure reduces the mathematical and computational complexity of the original case study. It is shown that the primer vector is independent of the semi-major axis of the transfer orbit. The case with a fixed transfer trajectory and variable initial and final thrust impulses is studied. The acquired related optimality maps are presented and analyzed and they express the likelihood of a set of trajectories to be optimal. Furthermore, it is presented which kind of requirements have to be fulfilled by a set of departure and arrival orbits to have the same profile of primer vector.     

初轨的稳健估计算法

讨论单测角资料(数据中包含异常值)初轨的稳健估计算法．该算法既可用于提供可靠的初轨,又可用作资料预处理,剔除异常值．稳健估计算法的崩溃点为16％至25％,可剔除6σ以上的异常值．     

Impulsive time-free transfers between halo orbits

A methodology is developed to design optimal time-free impulsive transfers between three-dimensional halo orbits in the vicinity of the interior L ₁ libration point of the Sun-Earth/Moon barycenter system. The transfer trajectories are optimal in the sense that the total characteristic velocity required to implement the transfer exhibits a local minimum. Criteria are established whereby the implementation of a coast in the initial orbit, a coast in the final orbit, or dual coasts accomplishes a reduction in fuel expenditure. The optimality of a reference two-impulse transfer can be determined by examining the slope at the endpoints of a plot of the magnitude of the primer vector on the reference trajectory. If the initial and final slopes of the primer magnitude are zero, the transfer trajectory is optimal; otherwise, the execution of coasts is warranted. The optimal time of flight on the time-free transfer, and consequently, the departure and arrival locations on the halo orbits are determined by the unconstrained minimization of a function of two variables using a multivariable search technique. Results indicate that the cost can be substantially diminished by the allowance for coasts in the initial and final libration-point orbits.An earlier version was presented as Paper AIAA 92-4580 at the AIAA/AAS Astrodynamics Conference, Hilton Head Island, SC, U.S.A., August 10–12, 1992.     

Optimal transfers between unstable periodic orbits using invariant manifolds

This paper presents a method to construct optimal transfers between unstable periodic orbits of differing energies using invariant manifolds. The transfers constructed in this method asymptotically depart the initial orbit on a trajectory contained within the unstable manifold of the initial orbit and later, asymptotically arrive at the final orbit on a trajectory contained within the stable manifold of the final orbit. Primer vector theory is applied to a transfer to determine the optimal maneuvers required to create the bridging trajectory that connects the unstable and stable manifold trajectories. Transfers are constructed between unstable periodic orbits in the Sun–Earth, Earth–Moon, and Jupiter-Europa three-body systems. Multiple solutions are found between the same initial and final orbits, where certain solutions retrace interior portions of the trajectory. All transfers created satisfy the conditions for optimality. The costs of transfers constructed using manifolds are compared to the costs of transfers constructed without the use of manifolds. In all cases, the total cost of the transfer is significantly lower when invariant manifolds are used in the transfer construction. In many cases, the transfers that employ invariant manifolds are three times more efficient, in terms of fuel expenditure, than the transfer that do not. The decrease in transfer cost is accompanied by an increase in transfer time of flight.     

A new approach to impulsive rendezvous near circular orbit

A new approach is presented for the problem of planar optimal impulsive rendezvous of a spacecraft in an inertial frame near a circular orbit in a Newtonian gravitational field. The total characteristic velocity to be minimized is replaced by a related characteristic-value function and this related optimization problem can be solved in closed form. The solution of this problem is shown to approach the solution of the original problem in the limit as the boundary conditions approach those of a circular orbit. Using a form of primer-vector theory the problem is formulated in a way that leads to relatively easy calculation of the optimal velocity increments. A certain vector that can easily be calculated from the boundary conditions determines the number of impulses required for solution of the optimization problem and also is useful in the computation of these velocity increments. Necessary and sufficient conditions for boundary conditions to require exactly three nonsingular non-degenerate impulses for solution of the related optimal rendezvous problem, and a means of calculating these velocity increments are presented. A simple example of a three-impulse rendezvous problem is solved and the resulting trajectory is depicted. Optimal non-degenerate nonsingular two-impulse rendezvous for the related problem is found to consist of four categories of solutions depending on the four ways the primer vector locus intersects the unit circle. Necessary and sufficient conditions for each category of solutions are presented. The region of the boundary values that admit each category of solutions of the related problem are found, and in each case a closed-form solution of the optimal velocity increments is presented. Similar results are presented for the simpler optimal rendezvous that require only one-impulse. For brevity degenerate and singular solutions are not discussed in detail, but should be presented in a following study. Although this approach is thought to provide simpler computations than existing methods, its main contribution may be in establishing a new approach to the more general problem.     

已知e的先验值的初轨寻优算法

本文讨论由短弧测角资料,已知偏心率e的先验值时,初轨的寻优算法。把偏心率的先验值看作是对e的约束条件,再固定观测弧段两端点的观测量与作为约束条件,本文导出了一维(对M_0)非线性寻优算法;针对观测资料中可能存在异常值的情形,本文提出了一种初轨的稳健寻优算法。仿真结果表明,当e的先验值有足够高的精度时,能明显提高周期P的估计精度;稳健算法还可用作剔除观测资料的异常值。     

Navigating to the Moon along low-energy transfers

This paper presents a navigation strategy to fly to the Moon along a Weak Stability Boundary transfer trajectory. A particular strategy is devised to ensure capture into an uncontrolled relatively stable orbit at the Moon. Both uncertainty in the orbit determination process and in the control of the thrust vector are included in the navigation analysis. The orbit determination process is based on the definition of an optimal filtering technique that is able to meet accuracy requirements at an acceptable computational cost. Three sequential filtering techniques are analysed: an extended Kalman filter, an unscented Kalman filter and a Kalman filter based on high order expansions. The analysis shows that only the unscented Kalman filter meets the accuracy requirements at an acceptable computational cost. This paper demonstrates lunar weak capture for all trajectories within a capture corridor defined by all the trajectories in the neighbourhood of the nominal one, in state space. A minimum Δv strategy is presented to extend the lifetime of the spacecraft around the Moon. The orbit determination and navigation strategies are applied to the case of the European Student Moon Orbiter.     

A novel algorithm for generating libration point orbits about the collinear points

This paper presents a numerical algorithm that can generate long-term libration points orbits (LPOs) and the transfer orbits from the parking orbits to the LPOs in the circular-restricted three-body problem (CR3BP) and the full solar system model without initial guesses. The families of the quasi-periodic LPOs in the CR3BP can also be constructed with this algorithm. By using the dynamical behavior of LPO, the transfer orbit from the parking orbit to the LPO is generated using a bisection method. At the same time, a short segment of the target LPO connected with the transfer orbit is obtained, then the short segment of LPO is extended by correcting the state towards its adjacent point on the stable manifold of the target LPO with differential evolution algorithm. By implementing the correction strategy repeatedly, the LPOs can be extended to any length as needed. Moreover, combining with the continuation procedure, this algorithm can be used to generate the families of the quasi-periodic LPOs in the CR3BP.     

基于CEI测量的GEO目标快速轨道恢复

连线干涉测量(Connected Element Interferometry, CEI)是一种全天时全天候的被动测角技术, 已用于空间目标的跟踪监视. 地球静止轨道(Geostationary Earth Orbit, GEO)卫星需要频繁机动以保持轨位或完成其他任务, 其机动后的快速轨道恢复能力对于监视预警极为重要. 针对基于CEI的GEO短弧定轨和预报, 分析了定轨算法的形亏和数亏, 在附加先验轨道约束的短弧定轨基础上, 提出了轨道半长轴初值的自适应优化方法. 利用亚太七号卫星的CEI仿真和实测数据进行了短弧定轨和预报, 实验结果表明, 采用优化后的半长轴初值, 30min短弧定轨和10min预报的卫星位置分量精度均优于4km, 能够满足非合作GEO目标机动后快速轨道恢复的需求.