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

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

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

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

轨道拦截预警算法研究

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

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

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

基于脉冲响应相关函数的模态辨识算法

利用不同时刻系统测量噪声具有独立同分布的特点,对脉冲响应进行相关处理,可以提高信噪比。为解决传统的基于脉冲响应的模态辨识算法的辨识精度问题,用某点脉冲响应同参考点响应之间的相关函数代替该点的脉冲响应以形成基于脉冲响应相关函数的模态辨识算法,并从理论上证明脉冲响应相关函数算法的辨识精度优于传统的脉冲响应算法。通过对一个四自由度系统模态频率和阻尼比的仿真辨识,验证了脉冲响应相关函数辨识算法的有效性。     

系统辨识在导弹建模中的应用

引入系统辨识的方法 ,确保当原始气动参数表内存在“突跳点”时 ,仍可以较为精确地得到导弹弹体三通道线性化模型。论述了这种方法的可行性、以及具体实现的步骤 ,最后给出仿真框图、仿真结果以及辨识误差的统计特性。结果表明 ,建模精度可以满足工程设计要求。     

Comparative analysis of four different single-frequency PPP models on positioning performance and atmosphere delay retrieval

Single-frequency precise point positioning (SF-PPP) has attracted increasing attention due to its high precision and cost effectiveness. With various strategies to handle the dominant error, i.e., ionosphere delay, the ionosphere-float (IF), ionosphere-free-half (IFH), ionosphere-corrected (IC), and ionosphere-weighted (IW) SF-PPP models are certain to possess different characteristics and performance levels. This study is dedicated to assessing and comparing the four models from model characteristics, positioning performance, and atmosphere delay retrieval. The model comparison shows that IC and IW models are full-rank while IF and IFH models have a rank deficiency of size one that will result in biased estimations, which means the better solvability of IC and IW models. The experiments are carried out based on the 7-day Global Positioning System (GPS) observations collected at 57 global Multi-GNSS Experiment (MGEX) stations and Global Ionosphere Map (GIM) products. The results indicate that the IW model can accelerate SF-PPP convergence and achieve higher positioning accuracy compared to the other three SF-PPP models, especially in kinematic mode. With convergence criteria of 0.25 m in horizontal and 0.5 m in vertical, the east/north/up convergence times of IW model are 0.5/15.0/25.0 min and 0.5/16.0/36.5 min for static and kinematic modes, respectively. The IW model is able to achieve an instantaneous positioning accuracy of 0.28/0.35/0.75 m. In addition, a real kinematic test also demonstrates the best positioning solutions of IW model. Regarding troposphere delay retrieval, the IF, IFH, and IW models obtain a comparable daily accuracy of 3.0 cm on average, while the IC model achieves the worst accuracy of 8.0 cm. For precise ionosphere delay estimation, IW model only needs an average initialization time of 34.3 min, but a longer initialization time of 51.6 min is required for IF model. The daily precision of ionosphere delay estimation for IW model can reach up to 10.8 cm. At the present accuracy of GIM products, it is suggested that the IW model should be adopted for SF-PPP first due to its superior performance in positioning and atmosphere delay retrieval.     

一种基于自测试的卫星快检技术研究

传统卫星设计以功能实现为主，可测试性较低，卫星综合测试基本依赖遥测遥控，测试周期长、效率低、成本高，不能很好地适应未来卫星批量化研制需求。参考民用汽车快速电检技术，从卫星可测试性与自测试设计出发，设计了基于总线调度的卫星快检系统，提高了卫星快速检修与快速发射能力，满足未来卫星低成本高效率研制需求。     

Comparison of convergence time and positioning accuracy among BDS,GPS and BDS/GPS precise point positioning with ambiguity resolution

Precise point positioning (PPP) usually takes about 30?min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60?days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0?cm, 1.5?cm, and 5.9?cm respectively, while those of the Gfix are 0.8?cm, 0.5?cm, and 1.7?cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7?cm, 0.5?cm, and 1.9?cm for the E, N, U components, respectively.     

二阶n次多项式自治系统的极限环的唯一性

主要利用文［１］中的变换，将下列二阶ｎ次多项式自治系统ｄｘｄｔ＝ｇ１（ｘ）＋ｈ１（ｘ）ｙｄｙｄｔ＝ｇ２（ｘ）＋ｈ２（ｘ）ｙ（＊）（其中，ｇ１（ｘ）＝∑ｎｉ＝０ａｉｘｉ，ｈ１（ｘ）＝∑ｎ－１ｉ＝０ｂｉｘｉ，ｇ２（ｘ）＝∑ｎｉ＝０ｃｉｘｉ，ｈ２（ｘ）＝∑ｎ－１ｉ＝０ｄｉｘｉ）变换成Ｌｉｅｎａｒｄ方程，再利用构造Ｄｕｌａｃ函数的方法和文［２］中的一个定理，得到了二阶ｎ次多项式自治系统（＊）的极限环唯一性的几个充分条件。