1、第 卷第期 年月系统工程与电子技术 文章编号：（）网址：收稿日期：；修回日期：；网络优先出版日期：。网络优先出版地址：基金项目：国家自然科学基金（，）；空间智能控制技术重点实验室开放基金（，）资助课题通讯作者引用格式：龚柏春，金鑫，张仁勇，等基于相对角动量的空间非合作目标机动仅测角检测方法系统工程与电子技术，（）：，（）：基于相对角动量的空间非合作目标机动仅测角检测方法龚柏春，金鑫，张仁勇，陈修桥（南京航空航天大学航天学院，江苏 南京 ；中国科学院空间应用工程与技术中心，北京 ；中国人民解放军 部队，北京 ）摘要：针对传统方法在空间非合作目标小脉冲轨道机动检测过程中性能不佳的问题，在构建了一种
2、基于相对角动量的新型待检特征量的基础上，提出了利用无源光学相机仅测角信息的轨道机动检测新方法。首先，在轨道坐标系下建立了基于 方程的相对动力学模型和相机偏离质心安装时的相对视线角测量模型；然后，给出了相对角动量的定义并量化分析了以相对角动量作为机动检测特征标志的优越性；接着，设计了基于二次滑窗方差比值的机动检测算法框架；最后，通过数值仿真对所提算法进行了有效性验证和性能分析。仿真结果表明，该方法可以实现对近圆轨道目标的快速轨道机动检测，机动时间点检测误差在 以内。关键词：空间非合作目标；轨道机动检测；相对角动量；仅测角中图分类号：文献标志码：，（.，；.，；.，）：，：；第期龚柏春等：基于相对
3、角动量的空间非合作目标机动仅测角检测方法 引言近年来，航天技术快速发展，各国的航天活动日益频繁，近地轨道上卫星数量急剧增长，拥挤的太空环境给卫星的安全带来了威胁 。根据空间安全威胁的来源和性质，可将其分为两类：第一类源于自然环境因素，如小行星碎片撞击；第二类则来自人为威慑和威胁，如卫星武器 。因此，实现对空间非合作目标的态势感知与跟踪，对保障我国现役卫星的安全而言具有重要意义。基于星载测量传感器实现对非合作目标的轨道机动检测已经成为了天基态势感知与跟踪领域需要突破的关键技术，已经得到了国内外的广泛关注和研究 。根据技术实现途径，现有的空间非合作目标轨道机动检测方法主要可分为两类，现综述如下。第
4、一类研究侧重于机动特征量的选取 ，不对测量方法有过高或特定的要求，如 等 根据卫星轨道半长轴的长期变化，提出了地球同步轨道、倾斜地球同步轨道和中地球轨道卫星面内轨道机动的检测方法，以及地球同步轨道卫星面外轨道机动的检测方法。等 使用半长轴和偏心率作为检测量来判断目标是否机动，但其方法仅适用于检测在近地点或远地点沿速度方向的脉冲机动。等 提出了使用公开的地球同步经度数据库进行神经网络训练来预测卫星状态，以星下点的经纬度和卫星高度作为特征量进行检测，但该方法仅适用于地球同步轨道，具有较高的局限性。第二类研究侧重于对卫星状态参数测量及数据处理方法进行探讨 ，以此来检测采用传统测量方法因精度不够而无法
5、识别的轨道机动，或者直接通过特定测量计算量的代价忽略机动存在以实现稳定跟踪。王庆瑞等 将距离变化率设定为检验参数，提出了基于概率判决的轨道机动检验方法，基于 准则，能够生成机动的判决门限，从而实现自适应的检验轨道机动，最终可以实现一般情况下正确率大于 的机动检测。等 提出使用颗相互协作的卫星来解决卡尔曼滤波法在目标机动后检测性能降低的问题，但该方法涉及到测量卫星相互协作的问题，过程复杂，且需要考虑到某一个或多个测量卫星无法观测时对系统的调整问题。等 提出了联合滤波的方法，其利用主、辅助滤波器的配合消除了机动对目标跟踪的影响。可见，第二类方法计算量较大，给卫星增加了负担，因此如何用普通的跟踪滤波
6、器对空间机动目标进行跟踪，是一个难点 。为了解决上述问题，本文针对空间非合作目标近程态势感知任务，在构建一种新型机动检测特征量的基础上，提出利用无源光学相机仅测量角度信息的机动检测和跟踪方法。下文将在第节和第节中建立相对动力学模型和传感器测量模型，在第节中设计基于新型特征量的机动检测算法，在第节中进行数值仿真实验与性能分析。相对运动动力学模型如图所示，建立当地垂直当地水平（，）参考坐标系，原点在感知卫星的质心，轴由地心指向卫星质心，轴垂直于轨道，轴满足右手定则。图地心坐标系（左）和 系（右）（）（）在近圆轨道、两个航天器之间的距离足够小且不考虑额外摄动的情况下，可以采用线性化的 方程来建模航天
7、器之间的相对运动 。时刻的相对轨道状态可以用时刻的状态向量表示：（）式中：为状态转移矩阵，具体形式参见文献 。传感器测量模型本文假设感知卫星自身的轨道和姿态已知，星上仅有无源光学相机可以用来对非合作目标进行相对测量。根据工程实际约束情况，相机偏离卫星质心安装，偏心测量如图所示。因此，可以直接在 参考坐标系下建立俯仰角和偏航角的测量模型：（）（）（）（）式中：，是目标时刻的相对位置；，为相机的偏置量；，为测量误差，设其为零均值高斯白噪声。图相机偏置测量示意图 系统工程与电子技术第 卷基于相对角动量的机动检测算法本文所研究的非合作目标轨道机动仅测角检测方法的基本逻辑关系如图所示。首先以式（）所示的
8、视线角测量量结合式（）所示的动力学模型进行滤波估计，获得非合作目标的相对轨道参数，然后基于相对轨道参数进行轨道机动待检量的构造，最后通过对待检量进行检验判别，来确定机动时间，并下达滤波系统重新初始化指令。图算法逻辑框图 由于本文的重点在于机动检测算法设计，所以对图中的仅测角定轨算法不过多论述，而是直接采用 等在文献 中提出的方法。视线角与相对轨道参量的映射关系如下：，（，）（）下文中的检测算法设计与分析过程将直接采用和进行。下面将重点先阐述基于相对角动量的新型机动待检特征量的构建及其敏感性分析，然后给出采用该特征量的轨道机动自动检测算法。机动待检特征量建模及其敏感性分析传统的轨道机动检测特征量
9、存在敏感性差、适用性窄等问题，难以满足空间对抗对态势感知灵敏性、检测快速性的要求。因此，本文提出一种基于相对角动量概念的新型机动检测特征量。下面首先阐述该待检特征量的定义，然后对其敏感性进行分析、证明。定义相对角动量为 （）式中：，（）卫星的机动会直接影响卫星的位置和速度矢量，而角动量的计算与物体的位置和速度矢量都具有相关性，其通过将速度和位置矢量进行叉乘的处理方式，不仅可以反映位置和速度矢量大小的变化，也可以直接通过角动量数值的变化反映位置和速度矢量方向上的变化，因此相对角动量对卫星机动更加敏感。此外，相对角动量参数对卫星所处轨道以及轨道上的位置没有要求，所以也将更具普适性。同时，与地心惯性
10、系下的轨道角动量相比，由式（）定义的轨道系下的相对角动量对目标机动更具敏感性，现将原因分析如下。卫星携带能量有限，所以卫星机动的冲量较小 ，这使其机动后的短时间内的速度变化量也很小，而卫星在低地球轨道运行时速度要达到 量级。因此，卫星在机动后短时间内的惯性位置、速度改变量相较于其机动前的位置大小、速度大小是极小的。而以地心为参考点的角动量由相对于地心的位置和速度贡献，这使惯性系下的角动量变化也很难在噪声中被发现。然而，若取测量卫星为参考点，轨道系下的相对角动量由相对位置与相对速度贡献，大大减小了变化量与原始基础量之间的差距，因此可以有效减小相对角动量的变化量被噪声淹没的风险。下面对上述定性分析
11、进行数学证明。根据式（）可以得到机动前时刻和机动完成后时刻的相对角动量分别为 （）当机动为脉冲机动时，由于机动时间很短，即，其中为机动时间。因此，（）则 （）（）记 （）与式（）式（）的过程同理，可以得到地心惯性系下的角动量为（）（）式中：，；，。为 参考坐标系到地心惯性坐标系的变化矩阵。由于常 见 的 单 个 脉 冲 所 能 提 供 的 加 速 度 大 小 约 为 ，而卫星相对于地心的运行速度至少达到 量级，所以有（）从而有（）式中：（）。因此，以地心惯性系下的角动量为检验标准，将很难在短时间内从噪声中辨别出是否机动。相反地，（）从而有 （）由此可以看出，相对角动量作为检验标准，将具有更优良
12、的特性。当考虑相对轨道估计误差时，式（）变为 （）（）（）（）第期龚柏春等：基于相对角动量的空间非合作目标机动仅测角检测方法 （）式中：和分别为相对位置和速度误差。通常，估计误差远远小于状态量，即 （）结合式（）和式（）可得 （）式（）表明了由轨道机动造成的相对角动量的变化将显著地区别于由相对轨道估计误差引起的变化。因此，只要相对轨道估计误差满足式（）所示的条件，相对角动量对目标机动是非常敏感的。轨道机动点自动判别算法若目标卫星所受到的合力并非指向感知卫星，则其相对角动量会发生变化，变化规律为 （）式中：为目标卫星受到的相对于感知卫星质心的单位质量力矩。在目标卫星没有施加机动时，仅为地心引力以
13、及相关摄动力提供，因而的变化很平缓，即 的变化是连续的。换言之，如果 出现了突变的情况，那么一定是目标卫星进行了机动（或被碰撞）。因此，检测 的突变就能判别机动。本文设计了基于二次滑窗方差比值的自主寻找机动点方法，算法流程如图所示。图机动点检测算法框图 在每完成一次滑窗内的计算后，滑窗向后移动一个格，重新计算并记录结果。将通过对相对角动量进行两次方差计算后的数据记作 。通过计算方差的方式将相对角动量突变转化为数值大小上的突变。为了使数值突变更加明显，本文采用二次滑窗的方式来处理相对角动量。设置一个宽度为、步长为的数据滑窗（包括待检测点时刻，前后再各取个数据），的选取与定轨精度相关。滑窗内的方差
14、为（）（）（）（）式中：，分别表示时刻、时刻、时刻的各方向上的相对角动量；表示各方向上的平均相对角动量。虽然可以从序列的突变来检测机动，但由于目标星的相对位置和机动方式不同，的范围差异很大，无法设置固定的阀值来判断是否机动。为了解决这个问题，首先在滤波平稳后，计算该点到检测点 的平均值：（）（）（）式中：表示检测点时刻的 。然后，为了消除各种情况间的差异，计算各个时刻 与 （）的比值：（）（）同时，为了缩小未机动时由滤波造成的浮动，并凸显机动时刻值的变化，可以将的首位替换为相对于其的无穷大，如 。为了进一步扩大机动后的值与未机动时的值的差异，可以将所有的值乘以同一倍数，如 ，则有，（）式中：下
15、标表示由 参考坐标系下、或方向计算得到，为由实验得到的阀值。若满足式（），则说明滑窗内包含机动时刻数据，认为即为机动时刻。数值仿真验证假设目标星和感知卫星都在近圆轨道上飞行，二者之间的初始相对轨道如表所示。测角相机的安装矢量为，安装误差为，测角误差为零均值高斯白噪声，其均方差将在仿真算例中具体说明，感知卫星的姿态测量包含 的常值偏差和均方差为 的白噪声。参考轨道的积分步长为。表目标卫星相对于测量卫星的初始相对轨道参数 相对轨道初始相对位置 初始相对速度（）静止 共椭圆飞越 椭圆绕飞 面外振荡 为了验证该方法的有效性，首先设置目标卫星在 时刻分别沿 参考坐标系下轴、轴和轴添加 的脉冲。之后，分别
16、在 时刻和 时刻设置不同的脉冲机动。系统工程与电子技术第 卷下面针对相对角动量对轨道机动的敏感性、敏感器精度对检测的影响、机动检测相对误差和双脉冲机动下的检测结果与持续跟踪观测效果进行仿真实验验证和性能分析。相对角动量对机动反应的仿真结果为了验证本文所提出的相对角动量对目标卫星机动的敏感性，设置相机的测角均方差为 ，并以两种不同相对轨迹为例进行验证。图、图分别为 静止轨道添加方向脉冲、跳跃轨道添加方向脉冲后相对角动量的变化。脉冲添加点均在 时刻，脉冲大小为。图 静止轨道添加方向脉冲后的相对角动量变化 图跳跃轨道添加方向脉冲后的相对角动量的变化 在图、图中，前期为滤波收敛过程造成的数值波动，在滤波收敛之后到脉冲机动之前，各个轴上的相对角动量基本稳定。但是非常明显地，在机动脉冲施加点 之后，图中轴、轴和图中轴的相对角动量发生了突变。这就表明了相对角动量对轨道机动具有优良的敏感能力。测角精度对脉冲机动检测的影响结果该仿真用于验证仅测角导航情况下的机动检测有效性。首先，在测角均方差为 的情况下得到合适的阀值；接下来，在取阀值为的情况下，进行敏感器测角均方差分别为 、，脉冲机动大小分别为 和 ，
