1. 航天器自主光学导航技术概述光学导航作为深空探测任务中的核心技术其基本原理是通过分析目标天体在成像仪上的位置信息结合航天器姿态数据实现相对轨道确定。这项技术在近年来的小行星探测任务中展现出独特优势特别是在资源受限的小型航天器平台上。在DESTINY这类小型深空探测任务中传统的地面测控导航方式面临诸多挑战。首先深空通信链路带宽有限难以为高精度导航提供充足的测量数据支持其次地面测控存在时间延迟无法满足飞越关键阶段的实时控制需求。这使得自主光学导航成为小型深空探测器的必然选择。光学导航系统的核心组件包括光学成像系统如TCAP旋转望远镜星载计算机OBC姿态确定与控制系统AOCS导航算法软件其工作流程可概括为望远镜获取目标天体图像图像处理器提取天体中心位置导航算法处理位置信息和姿态数据输出相对轨道估计结果为姿态控制系统提供引导指令关键提示在小型航天器设计中光学导航系统的性能往往受限于平台资源。望远镜尺寸、计算能力、电源预算等因素都需要在系统设计阶段仔细权衡。2. 对准误差对导航精度的影响机制2.1 望远镜对准误差的来源与分类在DESTINY任务中TCAP旋转望远镜存在三类主要对准误差详见表3中的ID6参数类型A误差ϕA,θA,ψA源于航天器本体坐标系与星敏感器(STR)之间的对准偏差主要由星敏感器安装精度和地面标定残余误差导致典型值4.6×10⁻³ deg (3σ)类型B误差ϕB,δB由望远镜旋转轴与设计方向之间的偏差引起受机械加工精度和装配工艺影响典型值0.01 deg (3σ)类型C/D误差δC,δD反映望远镜光轴与旋转轴之间的不垂直度与光学元件加工和装调精度相关典型值0.01 deg (3σ)这些误差参数共同构成了望远镜的视线(LoS)方向偏差在数学上可以表示为旋转矩阵的扰动项。当望远镜进行旋转跟踪时这些误差会产生复杂的耦合效应。2.2 误差传播的数学模型视线方向误差对导航精度的影响可以通过几何关系建模。设理想视线方向为L⁰实际视线方向L受误差影响可表示为L C(Δ)·L⁰其中C(Δ)是包含所有对准误差参数的复合旋转矩阵。在STR坐标系下这个误差会导致测量的目标天体位置出现偏差δyδy H·δx v式中H是观测矩阵δx是状态误差v是观测噪声。对准误差主要通过改变H矩阵的有效性来影响导航精度。仿真结果显示图13当存在上述对准误差时仅位置滤波器的导航精度从0.4 km(1σ)劣化到4.6 km(1σ)误差放大约11倍。这种劣化在近距离飞越阶段尤为显著因为此时微小的角度误差会对应较大的位置偏差。实践经验在任务前期通过地面测试精确标定这些对准参数非常关键。但受限于小型航天器的测试条件和成本往往难以将误差降到足够低的水平这就需要算法层面进行补偿。3. 改进的UKF导航算法设计3.1 算法核心思想传统的光学导航滤波器通常只估计相对位置或将视线误差简化为固定偏置角如偏置模型滤波器。本文提出的改进算法创新性地将7个对准误差参数纳入状态向量x [rᵀ, ϕA, θA, ψA, ϕB, δB, δC, δD]ᵀ其中r是相对位置矢量。这种设计使得滤波器能够显式建模各类对准误差的物理特性捕捉望远镜旋转时的误差变化规律实现误差参数与轨道参数的协同估计状态向量的维数从传统的3维仅位置扩展到10维显著增加了算法的复杂性但也提高了对真实系统的建模精度。3.2 UKF算法实现细节无迹卡尔曼滤波(UKF)特别适合处理这类非线性估计问题。算法的主要步骤如下sigma点生成 Xₖ [x̂ₖ, x̂ₖ±√(nλ)Pₖ]其中λ是缩放参数n是状态维数P为协方差矩阵时间更新 Xₖ₊₁⁻ f(Xₖ) x̂ₖ₊₁⁻ ΣWᵢ⁽ᵐ⁾Xₖ₊₁⁻ Pₖ₊₁⁻ ΣWᵢ⁽ᶜ⁾(Xₖ₊₁⁻-x̂ₖ₊₁⁻)(Xₖ₊₁⁻-x̂ₖ₊₁⁻)ᵀ Q观测更新 Zₖ₊₁⁻ h(Xₖ₊₁⁻) ẑₖ₊₁⁻ ΣWᵢ⁽ᵐ⁾Zₖ₊₁⁻ K PₓₓPₓₓ⁻¹ x̂ₖ₊₁⁺ x̂ₖ₊₁⁻ K(yₖ₊₁-ẑₖ₊₁⁻) Pₖ₊₁⁺ Pₖ₊₁⁻ - KSKᵀ滤波器参数设置如表4所示其中过程噪声Q和观测噪声R根据系统特性精心调整。特别值得注意的是UKF中的λ参数设置为-7这种非标准值的选择是为了在高维状态下保持数值稳定性。3.3 计算复杂度优化10维状态的UKF需要生成2n121个sigma点每个时间步长需要进行21次状态和观测方程的计算。为满足星载计算机的实时性要求我们采取了以下优化措施矩阵运算优化利用对称性减少计算量并行计算将sigma点计算分配到多个线程定点数运算在保证精度前提下使用定点数查表法预先计算常用函数值这些优化使得算法在OBC上单次迭代耗时控制在0.18秒以内图24满足1Hz的实时性要求。4. 仿真验证与结果分析4.1 数值仿真设置基于DESTINY任务飞越Phaethon小行星的场景建立了高保真仿真环境。关键参数如表3所示初始相对距离500 km相对速度3.3 km/s初始位置误差90 km (T/R轴), 130 km (S轴)观测误差0.006 deg对准误差如前所述采用蒙特卡洛方法进行100次仿真评估不同滤波器的性能。对比的算法包括仅位置滤波器偏置模型滤波器估计ϕA,ψA提出的完整模型滤波器4.2 性能对比分析仿真结果图14表6显示无误差时仅位置滤波器精度0.4 km有误差时仅位置滤波器劣化至4.6 km偏置模型和完整模型均能将精度恢复至0.11 km在接近阶段距离5分钟两种改进算法表现相当。这是因为此时望远镜旋转角较小图17视线误差主要表现为固定偏置。然而在近距离飞越阶段图19-20当望远镜快速旋转时偏置模型滤波器的误差增大至0.115 km完整模型保持0.021 km的高精度这种差异源于两种算法对误差动态特性的建模能力。完整模型通过估计底层对准参数能够准确预测旋转导致的视线变化图21而偏置模型则无法适应这种变化。4.3 硬件在环验证为验证算法的工程实用性在DESTINY的OBC样机上进行了硬件在环测试HILS。系统架构如图22所示包含实时动力学仿真计算机星载计算机LEON3处理器128MB RAMSpaceWire通信网络数据记录系统测试结果图23表明算法在真实硬件上的性能与数值仿真一致计算时间分布如图24所示90%情况下0.155s最坏情况下0.18s满足1Hz周期要求内存占用方面算法需要约20MB的常驻内存主要消耗在状态向量存储协方差矩阵中间计算结果缓冲区这在128MB的OBC上完全可行同时为其他任务留出了充足资源。5. 工程应用建议基于本研究结果对于小型深空探测任务的光学导航系统设计提出以下建议对准误差标定发射前尽可能精确标定对准参数保留星上实时估计能力作为备份算法选择对于长距离巡航段可采用简化算法节省资源关键飞越阶段应使用完整模型算法计算资源分配确保导航算法有专用计算时间片为UKF预留足够内存空间系统集成导航系统与姿态控制系统紧密耦合建立标准化的数据接口协议测试验证开展全数字仿真、硬件在环、半物理测试特别关注快速旋转工况下的性能在DESTINY任务中这套算法已经成功应用于TCAP系统的工程设计。任务规划在飞越前1小时启动高精度导航模式确保在最近距离时获得最佳定位精度。未来随着立方星等超小型探测器在深空探测中的应用增多这类资源高效的自主导航算法将变得更加重要。可能的扩展方向包括与视觉SLAM技术结合引入深度学习进行特征提取多探测器协同导航通过持续优化自主光学导航技术有望成为小型深空探测器的标准配置推动太阳系探索的民主化和常态化发展。