申请/专利权人：南京信息工程大学

申请日：2023-03-28

公开（公告）日：2024-01-30

公开（公告）号：CN116774255B

主分类号：G01S19/27

分类号：G01S19/27;G01S19/37

优先权：

专利状态码：有效-授权

法律状态：2024.01.30#授权;2023.10.10#实质审查的生效;2023.09.19#公开

摘要：本发明公开了一种IGSO导航卫星每日重现周期预报方法，通过使用历史的广播星历数据，计算IGSO卫星以往一段时间的每日重现周期值，并构建相应IGSO卫星每日重现周期值时间序列，根据观察时间序列的特性，确定IGSO卫星每日重现周期值随时间变化的数学模型，模型包括趋势项部分和周期项部分，根据时间序列特征确定周期项中的部分参数，模型其余参数根据时间序列值进行最优估计得到，通过所建立的模型，准确预测未来数个星期的IGSO卫星每日重现周期值，可服务于预估卫星覆盖区域的偏移情况、卫星的稳定性特征分析、实时多路径消除应用以及为卫星轨道调整提供科学依据等方面。

主权项：1.一种IGSO导航卫星每日重现周期预报方法，其特征在于，包括以下步骤：步骤1：根据历史广播星历数据计算IGSO导航卫星以往的每日重现周期，通过在广播星历中得到的轨道摄动参数和轨道长半径得到某天的IGSO卫星重现周期值，即日平均值；步骤2：构建IGSO卫星每日重现周期值的时间序列，并根据连续数天的IGSO卫星周期值时间序列进行绘图，显示其随时间变化的基本特征；步骤3：确定IGSO卫星每日重现周期值随时间变化的数学模型及确定部分模型参数：根据上述步骤2中得到的IGSO卫星每日重现周期值的时间序列，得到关于时间序列的趋势项、周期项，通过该趋势项和周期项确定重现周期随时间变化的模型，并且在周期项中通过观察实现序列得出其角频率；其中时间序列呈现明显的变大或减小趋势，则趋势项如下： 3；公式（3）中，为时间序列的开始时间，为模型计算值所对应的时间，为待确定的斜率参数，为待确定的截距参数；若时间序列不存在趋势项，则可令为0；若时间序列呈现明显的周期特征，则周期项如下： 4；公式（4）中，为振幅，为角频率，为初始相位；若时间序列中存在多个频率不同的周期项，则可由多个正弦函数组成，若时间序列中不存在周期项，则可令为0；最终的IGSO卫星每日重现周期随时间变化的数学模型可表示成由趋势项和周期项合成的函数，即 5；步骤4：根据IGSO卫星每日重现周期值时间序列估计数学模型参数：根据所述趋势项与周期项，得出初始相位和振幅的最优参数；其中，该趋势项和一个周期项，可以表示为： 7；所述步骤3确定了角频率，公式（7），，，为四个待估参数，将公式（7）变形为： 8；公式（8）中、为两个中间变量，为，为，通过该式得到模型参数和；将所述步骤2中获取的IGSO卫星每日重现周期时间序列值代入公式（8），左边可写成向量形式，为： 9；右边的方程系数矩阵为：（10）；则参数，，，估计如下：（11）；公式（11）中，，，，为，，，的最优估计值，初始相位的最优估计值为： （12）；振幅的最优估值为：（13）；步骤5：根据所建立的数学模型预测IGSO卫星未来的每日重现周期值；根据步骤4的最优参数确定IGSO卫星每日重现周期值随时间变化的数学模型，并将某一时刻的数值带入得到当时的重现周期预报值。

全文数据：

权利要求：

百度查询： 南京信息工程大学 一种IGSO导航卫星每日重现周期预报方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD