发明名称 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法
摘要 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法,本发明涉及星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法。本发明的目的是为了解决现有星敏感器标定方法不能较好地完成对星敏感器的安装矩阵误差的标定,以及不能定期对星敏感器安装矩阵和导航系统偏差进行标定和校正的问题。按以下方案实现:一:信息θ<sub>t,m</sub>和X<sub>t,m</sub>;二:建立姿态信息和轨道参数信息测量模型;三:确定<img file="DDA0000718005810000011.GIF" wi="70" he="72" />和<img file="DDA0000718005810000012.GIF" wi="112" he="64" />四:计算<img file="DDA0000718005810000013.GIF" wi="61" he="75" />和<img file="DDA0000718005810000014.GIF" wi="96" he="72" />五:求取算术平均值;六:建立星敏感器实际的姿态安装矩阵和轨道参数信息校正模型;七:确定Δθ的方向;八:对六进行校正;九:进行姿态确定和轨道参数信息确定;十:每隔N个姿态重新执行。本发明应用于卫星姿态确定技术与卫星导航技术领域。
申请公布号 CN104792340A 申请公布日期 2015.07.22
申请号 CN201510249554.X 申请日期 2015.05.15
申请人 哈尔滨工业大学 发明人 李敏;耿云海;张迎春;陈雪芹;李化义;谢成清;盛靖
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法,其特征在于,一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法具体是按照以下步骤进行的:步骤一:在卫星星下点轨迹下的固定空间坐标上建立M个地面测量接收站,卫星向M个建立好的地面测量接收站发射激光信息,激光信息包括卫星姿态信息θ<sub>t,m</sub>和轨道参数信息X<sub>t,m</sub>数据,M为正整数;步骤二:地面测量接收站对步骤一中的卫星姿态信息θ<sub>t,m</sub>和轨道参数信息X<sub>t,m</sub>进行测量,得到地面测量接收站测量到的卫星姿态信息θ<sub>t,m</sub>和地面测量接收站测量到的轨道参数信息X<sub>t,m</sub>,并建立卫星姿态信息测量模型和轨道参数信息测量模型;步骤三:利用多矢量定姿方法结合步骤二中得到的卫星姿态信息测量模型,计算出由地面测量接收站确定的卫星姿态信息估计值<img file="FDA0000718005780000011.GIF" wi="95" he="81" />利用卫星轨道动力学模型,结合滤波方法计算出由地面测量接收站确定的轨道参数信息估计值<img file="FDA0000718005780000012.GIF" wi="108" he="78" />将由地面测量接收站确定的卫星姿态信息估计值<img file="FDA0000718005780000013.GIF" wi="76" he="77" />和由地面测量接收站确定轨道参数信息估计值<img file="FDA0000718005780000014.GIF" wi="87" he="77" />进行存储;步骤四:卫星姿态确定系统利用星敏感器测量到的卫星姿态信息θ<sub>t,n</sub>、卫星姿态运动学、卫星姿态动力学模型和滤波方法计算出星敏感器测量到的卫星姿态信息的估计值<img file="FDA0000718005780000015.GIF" wi="97" he="86" />并进行存储;卫星自主导航系统利用星敏感器测量到的轨道参数信息X<sub>t,n</sub>、轨道动力学模型和滤波方法计算出星敏感器测量到的轨道参数信息的估计值<img file="FDA0000718005780000016.GIF" wi="105" he="92" />并进行存储;步骤五:对步骤三中的由地面测量接收站确定的卫星姿态信息估计值<img file="FDA0000718005780000017.GIF" wi="75" he="78" />和由地面测量接收站确定的轨道参数信息估计值<img file="FDA0000718005780000018.GIF" wi="86" he="77" />求取算术平均值,即:<img file="FDA0000718005780000019.GIF" wi="628" he="124" />消除系统随机测量误差;其中,所述<img file="FDA00007180057800000110.GIF" wi="60" he="77" />为地面测量接收站确定的卫星姿态信息估计值的算术平均值;<img file="FDA00007180057800000111.GIF" wi="74" he="77" />为地面测量接收站确定的轨道参数信息估计值的算术平均值;对步骤四中星敏感器测量到的卫星姿态信息的估计值<img file="FDA00007180057800000112.GIF" wi="68" he="86" />和星敏感器测量到的轨道参数信息估计值<img file="FDA0000718005780000021.GIF" wi="80" he="85" />求取算术平均值,即:<img file="FDA0000718005780000022.GIF" wi="608" he="124" />消除系统随机测量误差;其中,所述<img file="FDA0000718005780000023.GIF" wi="54" he="86" />为星敏感器测量到的卫星姿态信息估计值的算术平均值;<img file="FDA0000718005780000024.GIF" wi="68" he="84" />为星敏感器测量到的轨道参数信息估计值的算术平均值;步骤六:建立星敏感器实际的姿态安装矩阵和轨道参数信息校正模型;步骤七:将步骤五中得到的地面测量接收站确定的卫星姿态信息估计值的算术平均值<img file="FDA0000718005780000025.GIF" wi="59" he="78" />和地面测量接收站确定的轨道参数信息估计值的算术平均值<img file="FDA0000718005780000026.GIF" wi="76" he="77" />发送给卫星,卫星将步骤五中得到的星敏感器测量到的卫星姿态信息估计值的算术平均值<img file="FDA0000718005780000027.GIF" wi="53" he="87" />与步骤五中得到的地面测量接收站确定的卫星姿态信息估计值的算术平均值<img file="FDA0000718005780000028.GIF" wi="58" he="76" />作差,即:姿态误差角信息△θ为:<img file="FDA0000718005780000029.GIF" wi="277" he="83" />卫星将步骤五中得到的星敏感器测量到的轨道参数信息估计值的算术平均值<img file="FDA00007180057800000210.GIF" wi="68" he="85" />与步骤五中得到的地面测量接收站确定的轨道参数信息估计值的算术平均值<img file="FDA00007180057800000211.GIF" wi="70" he="76" />作差,即:轨道参数误差信息△X为:<img file="FDA00007180057800000212.GIF" wi="319" he="84" />利用步骤五中地面测量接收站确定的卫星姿态信息估计值的算术平均值<img file="FDA00007180057800000213.GIF" wi="84" he="79" />步骤五中星敏感器测量到的卫星姿态信息估计值的算术平均值<img file="FDA00007180057800000214.GIF" wi="54" he="82" />和星敏感器实际的姿态安装矩阵<img file="FDA00007180057800000217.GIF" wi="84" he="74" /><img file="FDA00007180057800000215.GIF" wi="245" he="78" />确定姿态误差角信息△θ的正负方向;步骤八:将步骤七得到的具有正负方向的姿态误差角信息△θ和轨道参数误差信息△X代入星敏感器误差安装矩阵<img file="FDA00007180057800000216.GIF" wi="70" he="85" />和步骤六中的轨道参数信息校正模型,对步骤六中的星敏感器安装矩阵和轨道参数信息进行校正;步骤九:利用卫星结合姿态动力学、运动学方程和轨道动力学方程对步骤八中校正后的星敏感器安装矩阵和轨道参数信息进行姿态确定和轨道参数信息确定;步骤十:每隔N个姿态,重新执行步骤一至步骤九。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号