[发明专利]基于低轨道地球卫星的地影时刻预报的星上确定方法

摘要

本发明公开了一种基于低轨道地球卫星的地影时刻预报的星上确定方法，该地影时刻预报运行在星载计算机中。本发明首先输入一个预报时刻，然后判断预报时刻是否位于时间节点，若位于则采用牛顿下山法计算预报时刻的进出地影的纬度幅角，若位于轨道参数冻结范围内，则采用轨道参数解析算法计算预报时刻的进出地影的纬度幅角；最后利用纬度幅角－时刻关系反解得到预报时刻对应的进出地影时刻。本发明是以轨道要素表征的变换矩阵作为信息输入，通过拟定的判定角与进出地影关系，并利用星载计算机中设置的时间节点、上注星历，获得上注星历精确预报卫星在每个时间节点上的地影时刻。

主权项

一种基于低轨道地球卫星的地影时刻预报的星上确定方法，所述低轨道地球卫星中的星载计算机用于负责星上数据与程序的存储、处理以及各分系统的协调管理；利用星载计算机中的星上数据并结合卫星是否位于地影区域的角度来得到卫星进出地影的时刻，从而使得星载计算机中的地影时刻自主预报更加精确、迅速；其特征在于：所述的地影时刻自主预报包括有下列步骤；步骤一，获取地心赤道惯性坐标系下的位置矢量；(A)获取卫星V在地心赤道惯性坐标系O‑XiYiZi中的位置矢量，记为卫星‑位置矢量即卫星V在坐标系Oi中各个轴上的分量分别为：u为纬度幅角，单位为度；Ω为轨道升交点赤经，单位为度；i为轨道倾角，单位为度；(B)获取太阳A在地心赤道惯性坐标系O‑XiYiZi中的位置矢量，记为太阳‑位置矢量即太阳A在坐标系Oi中各个轴上的分量分别为：Λ为太阳黄经，单位为度；ε为黄赤交角，单位为度；(C)获取太阳与卫星的连线AV在地心赤道惯性坐标系O‑XiYiZi中的矢量，记为卫星‑太阳矢量且利用地球在太阳光照射下的几何阴影关系，设定太阳光为平行光照射地球的情况下，则有步骤二，获取地心赤道惯性坐标系O‑XiYiZi到星‑地‑日坐标系O‑XcYcZc的变换矩阵依据由两矢量的分量列阵求坐标变换矩阵方法对卫星‑位置矢量和太阳‑位置矢量进行转换处理，得到变换矩阵其中：K11表示卫星V在地心赤道惯性坐标系O‑XiYiZi的Xi轴上的位置矢量，即K12表示卫星V在地心赤道惯性坐标系O‑XiYiZi的Yi轴上的位置矢量，即K13表示卫星V在地心赤道惯性坐标系O‑XiYiZi的Zi轴上的位置矢量，即K31表示卫星V在地心赤道惯性坐标系O‑XiYiZi的Xi轴上的位置矢量与太阳A在地心赤道惯性坐标系O‑XiYiZi的Xi轴上的位置矢量的向量积，即u为纬度幅角，单位为度；Ω为轨道升交点赤经，单位为度；i为轨道倾角，单位为度；Λ为太阳黄经，单位为度；ε为黄赤交角，单位为度；K32表示卫星V在地心赤道惯性坐标系O‑XiYiZi的Yi轴上的位置矢量与太阳A在地心赤道惯性坐标系O‑XiYiZi的Yi轴上的位置矢量的向量积，即K32＝sinusinicosΛ‑cosucosΩsinΛsinε+sinucosisinΩsinΛsinε；K33表示卫星V在地心赤道惯性坐标系O‑XiYiZi的Zi轴上的位置矢量与太阳A在地心赤道惯性坐标系O‑XiYiZi的Zi轴上的位置矢量的向量积，即K21表示K31与K11的向量积，即K22表示K32与K12的向量积，即K23表示K33与K13的向量积，即步骤三，获取卫星‑太阳矢量在星‑地‑日坐标系O‑XcYcZc中的投影，即卫星‑太阳‑投影点卫星‑太阳‑投影点在O‑XcYcZc的Xc轴上的位置分量记为Yc轴上的位置分量记为Zc轴上的位置分量记为则E1为与轨道要素关联的第一系数，E2为与轨道要素关联的第二系数，F1为与轨道要素关联的第三系数，F2为与轨道要素关联的第四系数，F3为与轨道要素关联的第五系数；E1＝‑cosicosΛsinΩ+cosisinΛcosεcosΩ+sinisinΛsinε；E2＝sinΛcosεsinΩ+cosΛcosΩ；F1=14(cos2isin2Λcosϵ-sin2Λcosϵ)sin2Ω+14(cos2icos2Λ-cos2isin2Λcos2ϵ-cos2Λ+sin2Λcos2ϵ)cos2Ω+14sin2isin2ΛsinϵsinΩ-14sin2isin2Λsin2ϵcosΩ+14(2sin2icos2Λ+cos2icos2icos2Λ+2sin2Λsin2icos2ϵ+cos2Λ+2sin2Λcos2isin2ϵ+cos2isin2Λcos2ϵ+2sin2Λsin2ϵ+sin2Λcos2ϵ);]]>F2=12(cosicos2Λ-cosisin2Λcos2ϵ)sin2Ω-12cosisin2Λcosϵcos2Ω-12sinisin2Λsin2ϵsinΩ-12sinisin2ΛsinϵcosΩ;]]>F3=-14(cos2isin2Λcosϵ+sin2Λcosϵ)sin2Ω+14(cos2icos2Λ+cos2isin2Λcos2ϵ-cos2Λ+sin2Λcos2ϵ)cos2Ω-14sin2isin2ΛsinϵsinΩ+14sin2isin2Λsin2ϵcosΩ+14(-2sin2icos2Λ-cos2icos2Λ-2sin2Λsin2icos2ϵ+cos2Λ-2sin2Λcos2isin2ϵ-cos2isin2Λcos2ϵ+2sin2Λsin2ϵ+sin2Λcos2ϵ);]]>步骤四，依据卫星‑太阳‑投影点在Yc轴上的位置分量与Xc轴上的位置分量的比值，来表征判定角β的正切角关系在地球视为均匀球体且卫星轨道的偏心率e最小的情况下，卫星恰好处于进地影位置或者出地影位置，此时的判定角β设为初始常值β0，且其中π取值为3.14；Re为地球平均半径，单位为米；a为卫星的轨道半长轴，单位为米；在卫星轨道确定、以及任意一时间节点tw的初始常值β0确定的情况下，卫星恰好进出地影的关系记为f(utw)=F2sin2u+F3cos2u-tanβ0×E1sinu-tanβ0×E2cosu+F1=0;]]>步骤五，获取卫星进出地影的时刻；步骤(5‑1)，通过星载计算机的界面输入一个预报时刻tq，所述预报时刻tq的形式为年月日时分秒；步骤(5‑2)，判断所述预报时刻tq是否位于时间节点tw；(A)若预报时刻tq位于时间节点tw上，即tq＝tw，则采用牛顿下山法对进出地影关系f(u)进行迭代运算，并以前一个时间节点tw‑1的进出地影的纬度幅角作为迭代初值，得到预报时刻tq的进出地影的纬度幅角对于时间节点tw对应的卫星进出地影的纬度幅角的牛顿下山迭代关系为由于预报时刻tq位于时间节点tw上，即tq＝tw，能够得到为时间节点tw上的进地影纬度幅角迭代值；为时间节点上的出地影纬度幅角迭代值；为前一个时间节点tw‑1上的进地影纬度幅角迭代值；为前一个时间节点tw‑1上的出地影纬度幅角迭代值；δ为下山因子；为的进出地影的函数值；为的导数值；为的进出地影的函数值；为的导数值；(B)若预报时刻tq不位于时间节点tw上，且位于轨道参数冻结范围内，则采用轨道参数解析算法对内的时间节点tw的进出地影的纬度幅角进行计算，得到预报时刻tq的进出地影的纬度幅角轨道参数解析算法是指：先选取出属于轨道参数冻结范围内的时间节点tw；然后计算预报时刻tq对应的进地影纬度幅角即计算预报时刻tq对应的出地影纬度幅角即为tq处进地影纬度幅角的变化量；为tq处出地影纬度幅角的变化量；ΔE1＝‑(cosicosΛcosΩ+cosisinΛcosεsinΩ)WΩΔt；ΔE2＝(sinΛcosεcosΩ‑cosΛsinΩ)WΩΔt；ΔF1=12(cos2isin2Λcosϵ-sin2Λcosϵ)cos2ΩWΩΔt-12(cos2icos2Λ-cos2isin2Λcos2ϵ-cos2Λ+sin2Λcos2ϵ)cos2ΩWΩΔt+14sin2isin2ΛsinϵsinΩWΩΔt+14sin2isin2Λsin2ϵcosΩWΩΔt;]]>ΔF2=(cosicos2Λ-cosisin2Λcos2ϵ)cos2ΩWΩΔt+cosisin2Λcosϵsin2ΩWΩΔt-12sinisin2Λsin2ϵcosΩWΩΔt+12sinisin2ΛsinϵsinΩWΩΔt;]]>ΔF3=-12(cos2isin2Λcosϵ+sin2Λcosϵ)cos2ΩWΩΔt-12(cos2icos2Λ-cos2isin2Λcos2ϵ-cos2Λ+sin2Λcos2ϵ)sin2ΩWΩΔt-14sin2isin2ΛsinϵsinΩWΩΔt-14sin2isin2Λsin2ϵsinΩWΩΔt;]]>E1为与轨道要素关联的第一系数，E2为与轨道要素关联的第二系数，F1为与轨道要素关联的第三系数，F2为与轨道要素关联的第四系数，F3为与轨道要素关联的第五系数；ΔE1为与轨道要素关联的第一系数的变化量，ΔE2为与轨道要素关联的第二系数的变化量，ΔF1为与轨道要素关联的第三系数的变化量，ΔF2为与轨道要素关联的第四系数的变化量，ΔF3为与轨道要素关联的第五系数的变化量；Δt为tq相对于时间节点tw的时间间隔，则Δt＝tq‑tw；WΩ为升交点赤经Ω的平均变化率，且J2为地球引力势的二阶谐系数，J2＝1.0826300×10‑3，Re为地球平均半径，μ为地球引力常数，i为轨道倾角，a为轨道半长轴，e为轨道偏心率；步骤(5‑3)，利用纬度幅角－时刻关系对进行反解，得到对应的进地影时刻对应的进地影时刻TΩ为卫星轨道的交点周期；tp为卫星运行的当前时刻；tp+1为卫星运行的下一时刻；为tp+1时刻所对应的纬度幅角；为tp时刻所对应的纬度幅角。