[发明专利]动态星图拖尾星像的数值模拟方法

说明书

(一)技术领域

本发明涉及一种动态拖尾星像的数值模拟方法，属于图像分析领域，是一项动态星图模拟的关键技术。

(二)背景技术

作为飞行器主要的姿态测量装置之一，星敏感器定姿精度最高，误差不随时间积累，应用越来越广泛。星敏感器工程应用最大的瓶颈性制约因素是载体机动导致的动态成像拖尾，星敏感器成像需要一定的曝光时间，在曝光时间内，如果捷联载体存在姿态机动，星光不再单纯成像为能量聚集的单个星像点。而是一条能量沿拖尾轨迹延展、峰值灰度显著下降、信噪比衰减的星像，称为拖尾星像。星模拟器是星敏感器的试验测试装置，星敏感器动态图像的处理能力的检验测试，需要星模拟器生成具有拖尾特征的动态星图，星像动态拖尾已严重影响到星敏感器的功能，不可以再利用静态星图序列检验捷联星敏感器算法的动态适应能力。

当星敏感器相对惯性系绝对静止时，静态星像呈近似圆形，能量也近似为二维高斯分布，如图1所示。如果捷联载体相对惯性系存在姿态机动，则在星敏感器本体产生相应的三轴角速率，导致动态成像拖尾，如图2所示，图中星k相对星敏感器本体具有相对角运动。理想拖尾星像覆盖范围内任一像素的能量，是幅值衰减的高斯点扩散函数沿拖尾轨迹的线积分。但程序中做线积分时间长，难以保障实时性，需要探索高实时性的拖尾星像数值模拟方法。

(三)发明内容

本发明的目的在于提供一种动态拖尾星像的数值模拟方法。

本发明的目的是通过以下技术方案实现的：

鉴于飞行器与星敏感器的捷联安装方式，利用安装矩阵可以实现飞行器本体和星敏感器本体姿态的相互转换。利用安装矩阵，将飞行器本体的三轴角速率转换为星敏感器本体的三轴角速率，进而可以得到星敏感器的姿态四元数运动学方程。

(1)确定龙格库塔积分步长及姿态四元数序列解算。星像点二维高斯灰度分布的均方差作为星敏感器的一个内参数，定义为星像点的高斯半径。根据高斯半径确定出时间步长。基于星敏感器的姿态四元数运动学方程并利用星敏感器四元数初值，进行龙格库塔数值积分解算，依次解出在时段内的星敏感器姿态四元数序列。

(2)基于四元数序列求解视场内所有星在星敏感器本体系下的星光矢量序列。

(3)依据光学系统模型，将星光矢量序列映射到星敏感器成像阵列面上，依次得到视场中每颗星的映射位置点序列，作为动态拖尾星像中心轨迹的离散点，参看图2动态拖尾成像机理及数值模拟方法示意图。

(4)以各映射位置点序列为中心，建立二维高斯灰度点扩散函数，对不同位置点处的高斯点扩散函数作叠加，即可模拟出每颗星的动态拖尾星像。图3为一条动态拖尾星像的仿真图。

(四)附图说明

图1为静态星像的灰度二维高斯分布仿真图；

图2为动态拖尾成像机理及数值模拟方法示意图；

图3为一条动态拖尾星像的仿真图。

图4为应用本发明方法模拟动态拖尾星像的实施方式流程图。

(五)具体实施方式

下面以视场内某颗星k为例，其赤经赤纬为(α_k,δ_k)，结合附图对本发明方法作进一步的描述和说明：

(1)星敏感器姿态四元数定义为：Q＝q₀+q₁i+q₂j+q₃k，(i，j，k分别为单位矢量，q₀为标量)，亦可表示为向量Q＝[q₀,q₁,q₂,q₃]^T。基于四元数Q的星敏感器姿态四元数运动学方程如下：