[发明专利]一种星历无更新条件下天基物联网终端的轨道预报方法

权利要求书

1.一种星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，该方法包括

数据初始化，并获取预设时刻的TLE数据；

调用两行根数拟合预测模型，基于所述预测模型及预设的所述TLE数据，计算两行根数预测值；

根据所述两行根数预测值，计算卫星的入境与出境时间。

2.如权利要求1所述的星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，所述数据初始化，具体包括：

初始化预报的起始时刻、预报天数、终端的经纬度信息、卫星入境仰角阈值。

3.如权利要求2所述的星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，所述根据所述两行根数预测值，计算卫星的入境与出境时间，具体包括：

将所述两行根数预测值代入到SGP4卫星轨道预报模型中，计算得到卫星在地心地固坐标系下，从预报起始时刻到终止时刻的卫星位置矢量；

根据所述终端经纬度信息，计算终端在地心地固坐标系下，从预报起始时刻到终止时刻的终端位置矢量；

根据卫星位置矢量与终端位置矢量，计算从预报起始时刻到终止时刻的终端到卫星的仰角，并在判定所述仰角大于或等于卫星入境仰角阈值时，计算卫星的入境与出境时间。

4.如权利要求1至3中任一项所述的星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，所述两行根数预测值中包括：平近点角值，所述计算两行根数预测值，具体包括：

根据建立的平近点角拟合预测模型，计算预报起始时刻的所述平近点角值，所述平近点角值的计算公式为：

m1＝(y0+((tx-t0)*n1*360/(24*3600)))％360.0

式中，m1为平近点角，y0为参考点时刻的平近点角，tx-t0为预报起始时刻距参考点时刻的秒数，n1为每天环绕地球圈数平均值，％为取模运算。

5.如权利要求4所述的星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，所述每天环绕地球圈数平均值n1为从参考点时刻到预报起始时刻以天为间隔，对预测的卫星每天的环绕地球圈数n2进行平均值计算所得的均值，其中，预测的卫星每天的环绕地球圈数n2的计算公式为：

n2＝p11*t2^3+p12*t2^2+p13*t2+p14

式中，n2为卫星每天的环绕地球圈数，t2为当前时间距最早的样本历史TLE历时的秒数，p11、p12、p13、p14为第一卫星拟合参数。

6.如权利要求1至3中任一项所述的星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，所述两行根数预测值中包括：近地点幅角预测值，所述计算两行根数预测值，具体包括：

根据建立的近地点幅角拟合预测模型，计算预报起始时刻的所述近地点幅角预测值，所述近地点幅角预测值的计算公式为：

ω1＝(p21*t+p22)％360.0

式中，ω1为计算的过程值，ω2为所述近地点幅角预测值，t为预报起始时刻距最早的样本历史TLE历时的秒数，p21、p22为第二卫星拟合参数。

7.如权利要求6所述的星历无更新条件下天基物联网终端的轨道预报方法，其特征在于，所述方法还包括：对所述近地点幅角预测值进行修正，并根据修正后的近地点幅角预测值，计算所述卫星的入境与出境时间，其中，修正所述近地点幅角预测值的计算公式为：

式中，ω3为修正后的近地点幅角预测值，％为取模运算,x1、x2、x3为修正值，ω2_1、ω2_2、ω2_3、ω2_4为修正角度区间值。