一种基于模型参数估计的星光矢量蒙气差修正方法技术

本发明专利技术公开了一种基于模型参数估计的星光矢量蒙气差修正方法，该方法可用于大气内星敏感器参考星矢量蒙气差修正。该方法将关于外部大气特征的蒙气差模型参数设为待估计值，根据星光矢量夹角不变原理，构造了蒙气差参数估计量测方程，通过扩展卡尔曼滤波算法循环迭代实时星图数据，估计出稳定且跟随环境变化的蒙气差模型参数，能够实时对参考星矢量进行蒙气差修正，进而提高大气内星敏感器姿态测量准确度。本发明专利技术实现了蒙气差模型参数的自主估计而不依赖其他外部传感器，采用扩展卡尔曼滤波算法实时估计出的蒙气差模型参数，不受传感器及经验模型误差影响。经验模型误差影响。经验模型误差影响。

【技术实现步骤摘要】
一种基于模型参数估计的星光矢量蒙气差修正方法


[0001]本专利技术属于天文导航、大气光学领域，具体涉及一种基于模型参数估计的星光矢量蒙气差修正方法。

技术介绍

[0002]星敏感器是一种以恒星作为参考，经星图捕获、星点质心提取、星图识别/跟踪、姿态解算等处理步骤，获得空间姿态的嵌入式传感设备。作为目前绝对姿态测量精度最高的传感器，在大气层外及临近空间的航行体上已得到广泛的应用。然而对于航行于大气层内，尤其在近地面空间的载体平台来说，使用星敏感器作为姿态测量器件，首先需要克服星光大气折射(蒙气差)导致的测量误差问题。
[0003]对于蒙气差修正技术，常用方法是获得观测星矢量/参考星矢量的天顶角以及观测环境的大气参数，根据蒙气差经验模型计算每个星光矢量的蒙气差经验值，再对星光矢量进行修正。所以这种方法就依赖于大气参数传感器(如温度、湿度、气压传感器等)，并受传感器误差及经验公式误差的影响(在修正观测星矢量时，还需依赖外部姿态传感器，并受该传感器准确度的影响。如具有漂移特点的陀螺仪等)。

技术实现思路

[0004]本专利技术的目的在于克服以往星光矢量蒙气差修正方法的不足，提出一种不依赖其他外部传感器并自适应获得蒙气差经验模型参数的方法。该方法利用星光矢量间角距不变的特点，建立蒙气差参数估计量测方程，再采用扩展卡尔曼滤波(Extended Kalman Filter,EKF)算法，建立蒙气差参数滤波方程，实现不依赖外部传感器而获得随大气变化的蒙气差经验模型参数，可以简单有效地进行参考星矢量蒙气差修正，进而提高大气内星敏感器姿态解算准确度。
[0005]本专利技术采用的技术方案为：一种基于模型参数估计的星光矢量蒙气差修正方法，包括以下步骤：
[0006]步骤一：星敏感器拍摄观测星图，获得成像面上星点坐标，根据星敏感器内参、畸变模型及小孔成像模型获得星敏感器坐标系下的若干观测星矢量，再进行星图识别，得到观测星与导航星表上参考星的匹配关系；
[0007]步骤二：根据星敏感器拍摄时间和地点以及地球自转等信息，基于国际天文联合会(International Astronomical Union,IAU)基础天文学标准(Standards Of Fundamental Astronomy,SOFA)提供的高精度坐标系转换链，可将步骤一获得的参考星从导航星表采用的天球坐标系转换到星敏感器观测位置的地理坐标系，即获得地理坐标系下无蒙气差的参考星矢量，进而可以获得参考矢量的真天顶距(无蒙气差的参考星矢量和天顶矢量的夹角)；
[0008]步骤三：选择满足工程需要的蒙气差经验模型，该模型将星光矢量的蒙气差大小表示为关于天顶距的三阶函数，各阶系数都是关于大气环境特征的函数，在本方法中设为
待估计的参数。在步骤二基础上，计算所有参考星矢量的蒙气差大小并对其沿仰角增加的方向做蒙气差修正，所得修正的参考星矢量都以待估计参数表示；
[0009]步骤四：根据星光矢量夹角不变原理，星敏感器坐标系下任意两观测矢量的夹角等于观测点地理坐标系下两对应的修正后参考矢量的夹角。所以可构造蒙气差参数估计量测方程，即以所有观测矢量对的夹角与对应参考矢量对的夹角之差表示观测变量列，该观测变量列可以构造为其关于待估计参数的雅可比矩阵和参数代差的乘积形式。
[0010]步骤五：在步骤四的基础上采用扩展卡尔曼滤波(Extended Kalman Filter,EKF)算法解算蒙气差模型参数：1、初始化蒙气差模型参数及误差协方差矩阵；2、计算步骤四所述观测变量列及观测变量列雅可比矩阵；3、计算卡尔曼滤波增益；4、计算参数代差；5、更新蒙气差模型参数及误差协方差矩阵；
[0011]步骤六：重复上述步骤一到步骤五，直至获得稳定的参数估计结果即可认为蒙气差模型参数估计值可用于修正参考星矢量。代入步骤三得到修正的参考星矢量。
[0012]其中，所述步骤一中，星敏感器拍摄观测星图，所得第i颗观测星在星敏感器坐标系下的观测星矢量记为v
i
。所述星图识别方法可采用子图同构类算法(如三角形算法、金字塔算法等)，也可采用模式识别类算法(如栅格算法、采用径向和环向特征的识别算法等)，还可以采用其他类型的识别算法(如基于神经网络的星图识别算法等)。所述星图识别算法不唯一，需结合实际星敏条件选择合适的算法，但需保证识别结果正确。
[0013]其中，所述步骤二中，参考星坐标系转换时，星敏感器拍摄时间为曝光中心时刻，表示方法可选UTC时间，精确到毫秒；地点为曝光中心时刻星敏感器的经度、纬度、高度，表示方法可选WGS84椭球参考系，精确到秒和米。地球自转信息包括地球当日极移和DUT1时间，可由国际地球自转服务(International Earth Rotation Service,IERS)提供。基于国际天文联合会(International Astronomical Union,IAU)基础天文学标准(Standards Of Fundamental Astronomy，SOFA)提供的高精度坐标系转换链进行参考星坐标系转换，还需依赖导航星表上参考星的赤经、赤纬、自行、视差、径向速度等参数。观测点地理坐标系即为东北天坐标系，所得地理坐标系下第i颗参考星无蒙气差的矢量记为w
i
，天顶矢量U＝[0，0，1]T
，真天顶距为：
[0014]Z
i
＝arccos(U
·
w
i
)
ꢀꢀ
(1)
[0015]其中，所述步骤三中，蒙气差经验模型如下：
[0016]ρ＝a1tanZ+a2tan3Z+a3tan5Z
ꢀꢀ
(2)
[0017]其中，a1，a2，a3即为待估计参数，ρ为星光矢量的蒙气差大小，Z为星光矢量的真天顶距；计算所有参考星矢量的蒙气差大小ρ
i
并对其沿仰角增加的方向做蒙气差修正，所得修正的参考星矢量都是关于待估计参数a1，a2，a3的函数，修正后的参考星矢量表示为
[0018][0019]其中，I表示单位矩阵，n
i
表示参考星矢量w
i
增加仰角ρ
i
时的转轴矢量：
[0020]n
i
＝w
i
×
U
ꢀꢀ
(4)
[0021]n
i
∧表示由矢量元素构成的反对称矩阵，n
i
(1)、n
i
(2)、n
i
(3)分别表示矢量n
i
的第1、2、3个元素：
[0022][0023]其中，所述步骤四中，根据星光矢量夹角不变原理，星敏感器坐标系下任意两观测矢量的夹角等于观测点地理坐标系下两对应的修正后参考矢量的夹角，可得：
[0024][0025]若有n对匹配的星光矢量，则可构造个观测变量r
ij
：
[0026][0027本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于模型参数估计的星光矢量蒙气差修正方法，其特征在于：所述方法包括以下步骤：步骤一，由星敏感器拍摄星图获得星敏感器坐标系下的观测星矢量，再经过星图识别获得观测星与参考星匹配关系；步骤二，将参考星从天球坐标系转换到星敏感器观测时的地理坐标系，并计算参考星矢量的真天顶距；步骤三，根据蒙气差经验模型和参考星矢量真天顶距计算星光矢量蒙气差值，再对参考星矢量沿仰角增加的方向做蒙气差修正，保留蒙气差模型参数为待估计值，得到修正的参考星矢量；步骤四，根据星光矢量夹角不变原理，构造蒙气差参数估计量测方程，并计算方程中观测变量列及其关于参数的雅可比矩阵；步骤五，采用扩展卡尔曼滤波算法解算并更新蒙气差模型参数和误差协方差矩阵；步骤六，重复上述步骤一到步骤五，直至获得稳定的参数估计结果并用于修正参考星矢量。2.如权利要求1所述的一种基于模型参数估计的星光矢量蒙气差修正方法，其特征在于：所述步骤一中，星敏感器拍摄观测星图，所得第i颗观测星在星敏感器坐标系下的观测星矢量记为v
i
。3.如权利要求1所述的一种基于模型参数估计的星光矢量蒙气差修正方法，其特征在于：所述步骤二中，参考星坐标系转换时，星敏感器拍摄时间为曝光中心时刻，表示方法可选UTC时间，精确到毫秒；地点为曝光中心时刻星敏感器的经度、纬度、高度，表示方法可选WGS84椭球参考系，精确到秒和米；地球自转信息包括地球当日极移和DUT1时间，可由国际地球自转服务(International Earth Rotation Service，IERS)提供；基于国际天文联合会(International Astronomical Union，IAU)基础天文学标准(Standards Of Fundamental Astronomy，SOFA)提供的高精度坐标系转换链进行参考星坐标系转换，还需依赖导航星表上参考星的赤经、赤纬、自行、视差、径向速度等参数。观测点地理坐标系即为东北天坐标系，所得地理坐标系下第i颗参考星无蒙气差的矢量记为w
i
，天顶矢量U＝[0，0，1]
T
，第i颗参考星的真天顶距为：Z
i
＝arccos(U
·
w
i
)
ꢀꢀꢀꢀ
(1)。4.如权利要求1所述的一种基于模型参数估计的星光矢量蒙气差修正方法，其特征在于，所述步骤三中，蒙气差经验模型如下：ρ＝a1tanZ+a2tan3Z+a3tan5Z
ꢀꢀꢀꢀ
(2)其中，a1，a2，a3即为待估计参数，ρ为星光矢量的蒙气差大小，Z为星光矢量的真天顶距；计算所有参考星矢量的蒙气差大小ρ
i
并对其沿仰角增加的方向做蒙气差修正，所得修正的参考星矢量都是关于待估计参数a1，a2，a3的函数，修正后的参考星矢量表示为的函数，修正后的参考星矢量表示为其中I表示单位矩阵，n
i
表示参考星矢量w
i
增加仰角ρ
i
时的转轴矢量：n
i
＝w
i
×
U
ꢀꢀꢀꢀ
(4)n
i∧
表示由矢...

【专利技术属性】
技术研发人员：江洁，张晗，杨季三，张广军，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：