买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网获悉上海卫星工程研究所申请的专利基于GNSS绝对定位数据的卫星姿态导引角计算方法及系统获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115892517B 。

龙图腾网通过国家知识产权局官网在2025-08-29发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211486970.8，技术领域涉及：B64G1/24；该发明授权基于GNSS绝对定位数据的卫星姿态导引角计算方法及系统是由栗双岭;温俊健;王赟;李楠;孙英梅;韩东升;仲兆宇设计研发完成，并于2022-11-24向国家知识产权局提交的专利申请。

本基于GNSS绝对定位数据的卫星姿态导引角计算方法及系统在说明书摘要公布了：本发明提供了一种基于GNSS绝对定位数据的卫星姿态导引角计算方法及系统，包括：步骤S1：根据GNSS实时绝对定位数据中的时间、位置、速度转换到惯性坐标系下，得到惯性坐标系下的位置、速度；步骤S2：将GNSS实时绝对定位的时间下位置速度数值积分递推至当前星时，得到当前星时时刻的位置和速度；步骤S3：将当前时刻的位置速度信息转换得到对应时刻的轨道参数；步骤S4：根据当前星时时刻轨道参数，计算当前时刻偏航导引角和俯仰导引角。本方法可以实时高精度计算卫星姿态导引角，有效支持高精度的姿态控制，有效补偿地球自转引起的图像形变，提高成像质量。

本发明授权基于GNSS绝对定位数据的卫星姿态导引角计算方法及系统在权利要求书中公布了：1.一种基于GNSS绝对定位数据的卫星姿态导引角计算方法，其特征在于，包括： 步骤S1：根据GNSS实时绝对定位数据中的时间、位置、速度转换到惯性坐标系下，得到惯性坐标系下的位置、速度； 步骤S2：将GNSS实时绝对定位的时间下位置速度数值积分递推至当前星时，得到当前星时时刻的位置和速度； 步骤S3：将当前星时时刻的位置速度信息转换得到对应时刻的轨道参数； 步骤S4：根据当前星时时刻轨道参数，计算当前星时时刻偏航导引角和俯仰导引角； 在所述步骤S1中： 步骤S1.1： 根据时间T0计算J2000地心惯性系到地固坐标系的坐标变换矩阵： MECI2ECF＝EP·ER·NR·PR 其中，EP为极移矩阵，简化为3×3单位阵；ER为地球自转矩阵，NR为章动矩阵，PR为岁差矩阵； 步骤S1.2： 将GNSS绝对定位数据位置速度信息由WGS84坐标系转换至J2000惯性坐标系，得到位置rECI0、速度vECI0为： rECI0＝MECI2ECF-1·rECF0 vECI0＝MECI2ECF-1vECF0+ωe ×rECF0 其中，rECFO为T0时刻WGS84坐标系下的位置，vECFO为T0时刻WGS84坐标系下的速度；上标T表示转置，上标·表示对时间求导； 在所述步骤S2中： 根据T0时刻位置速度，采用四阶龙格库塔公式，数值积分递推至当前星时T，得到T时刻的位置rECI、速度vECI； 在所述步骤S3中： 步骤S3.1：计算过程参数： 位置标量r＝|rECI|，速度标量v＝|vECI| 轨道坐标系到J2000惯性坐标系转换矩阵： 其中，Y＝vECI×rECI，X＝rECI×Y，均为3×1列向量；|X|表示X的模；轨道能量常量： 其中，地球引力常数μ＝398600km3s2； 动量矩： H＝rECI×vECI 动量矩的模： H＝|H| 半通径： p＝H2μ 拉普拉斯矢量： 其中，Lx、Ly、Lz分别为拉普拉斯矢量L在惯性坐标系中的三个坐标轴的分量；步骤S3.2：计算T时刻轨道根数： 半长轴： 偏心率： 轨道倾角： i＝arccosHziH 其中，Hzi为动量矩H在Z向分量； 纬度幅角： u＝arctan2-Morb2ECI33,Morb2ECI31 其中，Morb2ECI33为矩阵Morb2ECI第3行第3列的数值，Morb2ECI31为矩阵Morb2ECI第3行第1列的数值； 升交点赤经： Ω＝arctan2-Morb2ECI12,Morb2ECI22 其中，Morb2ECI12为矩阵Morb2ECI第1行第2列的数值，Morb2ECI22为矩阵Morb2ECI第2行第2列的数值； 近地点幅角： 真近点角： f＝u-ω。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人上海卫星工程研究所，其通讯地址为：200240 上海市闵行区元江路3666号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。