一种利于交班的舰空导弹最优中制导律设计

针对舰空导弹超视距拦截低空飞行目标问题,在复合制导的中制导段,建立了导弹拦截目标的相对运动模型,依据制导平台无线电指令修正信息,考虑舰空导弹拦截低空目标时的飞行速度特性,给出了预测遭遇点的实时解算模型,并确定了利于中末制导交班的中制导终端约束条件。基于拦截预测遭遇点,引入伪控制量的概念,应用线性二次型最优控制理论,设计了一种基于拦截预测遭遇点的利于制导交班的最优中制导律,仿真结果验证了预测遭遇点实时解算模型和最优中制导律的正确性与有效性。     

Libration point orbits for lunar global positioning systems

With the development of lunar exploration, a lunar global positioning system (LGPS) is demanded for both on-ground and in-flight lunar exploration missions. The traditional configuration of constellation requires at least eighteen satellites to cover the whole lunar surface continuously. In this paper, the configurations of the libration point orbits (LPOs) constellations are investigated. By using the constellations on the Earth–Moon _L1$L_1$ and _L2$L_2$ LPOs, the basic functions of the LGPS can be realized by using eight to fourteen satellites. First, the LPO and the combinations of LPOs, which can be used in the constellations of the LGPS, are investigated. The criteria and procedures of the configuration design are introduced. Second, the configurations of LPOs constellations are investigated in the Earth–Moon circular-restricted three-body problem (CR3BP). The size of the LPOs and the distribution of the satellites on these LPOs are determined by using an exhaustive algorithm and a global optimization method, respectively. The key performance parameters of these constellations are computed. Third, the constellations with good performance in the CR3BP are redesigned in the more accurate Earth–Moon based Sun-perturbed bicircular four-body problem (B4BP). Moreover, in order to avoid the ground coverage problem caused by the perturbation of the Sun, some modifications are implemented, and the configuration of the no blind area LGPS in the B4BP is obtained.     

Symmetry and asymmetry of ionospheric weather at magnetic conjugate points for two midlatitude observatories

Variations of the ionospheric weather W-index for two midlatitude observatories, namely, Grahamstown and Hermanus, and their conjugate counterpart locations in Africa are studied for a period from October 2010 to December 2011. The observatories are located in the longitude sector, which has consistent magnetic equator and geographic equator so that geomagnetic latitudes of the line of force are very close to the corresponding geographic latitudes providing opportunity to ignore the impact of the difference of the gravitational field and the geomagnetic field at the conjugate points on the ionosphere structure and dynamics. The ionosondes of Grahamstown and Hermanus provide data of the critical frequency (foF2), and Global Ionospheric Maps (GIM) provide the total electron content (TECgps) along the magnetic field line up to the conjugate point in the opposite hemisphere. The global model of the ionosphere, International Reference Ionosphere, extended to the plasmasphere altitude of 20,200 km (IRI-Plas) is used to deliver the F2 layer peak parameters from TECgps at the magnetic conjugate area. The evidence is obtained that the electron gas heated by day and cooled by night at the summer hemisphere as compared with the opposite features in the conjugate winter hemisphere testifies on a reversal of plasma fluxes along the magnetic field line by the solar terminator. The ionospheric weather W-index is derived from NmF2 (related with foF2) and TECgps data. It is found that symmetry of W-index behavior in the magnetic conjugate hemispheres is dominant for the equinoxes when plasma movement along the magnetic line of force is imposed on symmetrical background electron density and electron content. Asymmetry of the ionospheric storm effects is observed for solstices when the plasma diffuse down more slowly into the colder winter hemisphere than into the warmer summer hemisphere inducing either plasma increase (positive phase) or decrease (negative phase of W-index) in the ionospheric and plasmaspheric plasma density.     

一种弹目遭遇点预测方法

为充分利用舰空导弹中制导段所获信息并减小导弹中制导段的弹道曲率,给出一种弹目预测遭遇点的解算方法,并据此设计比例导引制导律。在初始发射坐标系中,假设目标由当前位置以当前速度大小沿当前速度方向匀速运动至遭遇点,导弹由当前位置以当前速度大小按照一定的导引规律匀速运动至遭遇点。在该假设条件下,导弹到达遭遇点的总航路大小一方面...     

小推力航天器的地月低能转移轨道

 在限制性四体模型下研究基于小推力方式的地月低能转移问题,通过借助于平动点轨道的相空间结构来揭示小推力转移的机理。重点研究了小推力转移自由飞行段的构造：经由LL₁点穿越获得最小能量的低能转移;而经由LL₁点Halo轨道穿越,得到(M,N)圈穿越轨道;由于Halo轨道相对于平动点增加了一维度的选择,根据(2,2)圈穿越轨道构造该转移的自由飞行段。在地球势阱逃逸和月球势阱捕获段,分别设计了合适的小推力的控制律及发动机开/关机时间,成功实施近地球段的小推力加速和近月球段的减速。尽管未对所得到的结果进行优化,所得转移轨道的燃料消耗也与类似边界条件的SMART-1轨道基本一致。     

基于标准弹道点的最优制导方法

由于基于需要速度的显式制导方法要确定虚拟目标点,虚拟目标点的确定要对地球扁率和再入阻力的影响进行修正,存在着百米以上的修正误差,为此提出了基于标准弹道点的最优制导方法,即将标准弹道的某点作为目标点进行需要椭圆轨道计算及最优入轨控制。提供了采用弹上迭代制导方法实现最优控制,同时给出了预测关机参数的需要椭圆轨道终端约束校正方法。该制导方法无需进行地球扁率和再入阻力的影响修正,可以减少方法误差,仿真计算证明了这一点。     

结构可靠性分析的拟蒙特卡罗方法

由于蒙特卡罗(MC)方法具有程序结构简单,收敛速度与问题维数无关等优点,故其在结构可靠性分析中得到了广泛应用。但是对于小失效概率等问题,计算效率低这一主要缺点限制了该方法的应用范围。为了解决MC方法存在的问题,通过引进单位超立方体中不同的低偏差点集代替伪随机数序列,并结合了重要抽样技术建立了结构可靠性分析的拟蒙特卡罗(QMC)方法。该方法不但可以大幅度减少抽样点数目,还能够得到确定性的估计值, 避免传统MC方法只能得到概率意义下误差的缺陷。通过数值算例可以看出本文方法具有较高的计算精度和效率,因此适用于结构可靠性分析问题。     

基于CDMA蜂窝网的飞行器三维定位算法

针对中低空飞行器提出了一种将码分多址(CDMA)蜂窝网导频信号被动定位与惯性导航系统进行组合的导航定位方法。该定位方法是根据CDMA蜂窝网导频信号的到达时差（TDOA）进行两次加权最小二乘解算,并利用飞行器惯性组合的原始信息进行解模糊,从而实现中低空飞行器的三维实时定位。蒙特卡罗仿真证明：在高斯白噪声条件下,该算法运算速度快,精度高,其均方误差逼近克拉美罗下界（CRLB）,为组合导航实时性要求提供了理论参考。     

轨道拦截预警算法研究

针对空间轨道拦截作战模式,提出了一种拦截预警的计算方法.推导了求解轨道交点位置的表达式,然后结合拦截卫星的末端威胁距离给出了判断拦截是否发生的地心距判据和时间差判据,最后通过仿真算例验证了交点表达式以及判据的正确性,并讨论了拦截卫星的威胁距离与时间差门限的关系.     

基于ICP算法的手术导航三维配准技术

针对计算机辅助手术三维导航技术中术前CT图像与术中实际空间的配准问题,提出一种基于最近点迭代(ICP,Iterative Closest Point)算法的特征点云配准技术.利用医学图像空间和实际空间特征区域的两片点云坐标进行三维配准.对CT图像进行重建、分割及交互式操作得到医学图像特征点云;利用光学定位仪实时采集实际空间中对应区域的点云;通过主元分析(PCA,Principal Component Analysis)获取两组点云数据的特征向量进行初配准;进行最近点迭代使配准矩阵收敛到一个最优解,其中采用k-d tree寻找邻近点加速迭代过程.以塑料脊柱模型骨为对象进行了脊柱手术导航配准精度实验,进一步对实验中的点云数据加入高斯噪声以进行误差分析.结果表明这种配准方法简单可靠,在模型骨情况下配准精度在1mm以内.