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

权利要求书

1.一种基于低轨道地球卫星的地影时刻预报的星上确定方法，所述低轨道地球卫星中的星载计算机用于负责星上数据与程序的存储、处理以及各分系统的协调管理；利用星载计算机中的星上数据并结合卫星是否位于地影区域的角度来得到卫星进出地影的时刻，从而使得星载计算机中的地影时刻自主预报更加精确、迅速；

其特征在于：所述的地影时刻自主预报包括有下列步骤；

步骤一，获取地心赤道惯性坐标系下的位置矢量；

(A)获取卫星V在地心赤道惯性坐标系O-X_iY_iZ_i中的位置矢量，记为卫星-位置矢量即卫星V在坐标系O_i中各个轴上的分量分别为：xOiV=cosucosΩ-sinucosisinΩ,yOiV=cosusinΩ+sinucosicosΩ,zOiV=sinusini;]]>

u为纬度幅角，单位为度；

Ω为轨道升交点赤经，单位为度；

i为轨道倾角，单位为度；

(B)获取太阳A在地心赤道惯性坐标系O-X_iY_iZ_i中的位置矢量，记为太阳-位置矢量即太阳A在坐标系Oi中各个轴上的分量分别为：yOiA=sinΛcosϵ,zOiA=sinΛsinϵ;]]>

Λ为太阳黄经，单位为度；

ε为黄赤交角，单位为度；

(C)获取太阳与卫星的连线AV在地心赤道惯性坐标系O-X_iY_iZ_i中的矢量，记为卫星-太阳矢量且DV-AOi=DAOi-DVOi,xOiV-A=xOiA-xOiV,yOiV-A=yOiA-yOiV,zOiV-A=zOiA-zOiV,]]>DV-AOi=(xOiV-A,yOiV-A,zOiV-A);]]>

利用地球在太阳光照射下的几何阴影关系，设定太阳光为平行光照射地球的情况下，DVOi=(xOiV,yOiV,zOiV)=(0,0,0),]]>则有DV-AOi=DAOi,DV-AOi=(xOiV-A,yOiV-A,zOiV-A)=DAOi=(xOiA,yOiA,zOiA);]]>

步骤二，获取地心赤道惯性坐标系O-X_iY_iZ_i到星-地-日坐标系O-X_cY_cZ_c的变换矩阵

依据由两矢量的分量列阵求坐标变换矩阵方法对卫星-位置矢量和太阳-位置矢量进行转换处理，得到变换矩阵LOcOi=K11K12K13K21K22K23K31K32K33,]]>其中：

K₁₁表示卫星V在地心赤道惯性坐标系O-X_iY_iZ_i的X_i轴上的位置矢量，即K11=xOiV=cosucosΩ-sinucosisinΩ;]]>

K₁₂表示卫星V在地心赤道惯性坐标系O-X_iY_iZ_i的Y_i轴上的位置矢量，即K12=yOiV=cosusinΩ+sinucosicosΩ;]]>

K₁₃表示卫星V在地心赤道惯性坐标系O-X_iY_iZ_i的Z_i轴上的位置矢量，即K13=zOiV=sinusini;]]>

K₃₁表示卫星V在地心赤道惯性坐标系O-X_iY_iZ_i的X_i轴上的位置矢量与太阳A在地心赤道惯性坐标系O-X_iY_iZ_i的X_i轴上的位置矢量的向量积，即K₃₁＝-sinusinisinΛcosε

+cosusinΩsinΛsinε；

+sinucosicosΩsinΛsinε

u为纬度幅角，单位为度；

Ω为轨道升交点赤经，单位为度；

i为轨道倾角，单位为度；

Λ为太阳黄经，单位为度；

ε为黄赤交角，单位为度；

K₃₂表示卫星V在地心赤道惯性坐标系O-X_iY_iZ_i的Y_i轴上的位置矢量与太阳A在地心赤道惯性坐标系O-X_iY_iZ_i的Y_i轴上的位置矢量的向量积，即K₃₂＝sinusinicosΛ-cosucosΩsinΛsinε+sinucosisinΩsinΛsinε；

K₃₃表示卫星V在地心赤道惯性坐标系O-X_iY_iZ_i的Z_i轴上的位置矢量与太阳A在地心赤道惯性坐标系O-X_iY_iZ_i的Z_i轴上的位置矢量的向量积，即K₃₃＝-cosusinΩcosΛ-sinucosicosΩcosΛ

                                                ；

+cosucosΩsinΛcosε-sinucosisinΩsinΛcosε

K₂₁表示K₃₁与K₁₁的向量积，即K₂₁＝cos²usin²ΩcosΛ+sinucosucosisin2ΩcosΛ

-sinucosucosisinΛcosεcos2Ω-cos²usinΩcosΩsinΛcosε

+sin²ucos²icos²Ω+sin²ucos²isinΩcosΩsinΛcosε；

+sin²usin²icosΛ-sinucosusinicosΩsinΛsinε

+sin²usinicosisinΩsinΛsinε

K₂₂表示K₃₂与K₁₂的向量积，即K₂₂＝-cos²usinΩcosΩcosΛ-sinucosucosicos2ΩcosΛ

+cos²ucos²ΩsinΛcosε+sin²usin²isinΛcosε

-sinucosusinisinΩsinΛsinε+sin²ucos²isinΩcosΩcosΛ；

-sinucosucosisin2Ω-sin²usinicosicosΩsinΛsinε

+sin²ucos²isin²ΩsinΛcosε

K₂₃表示K₃₃与K₁₃的向量积，即K₂₃＝-sinucosusinicosΩcosΛ+cos²usinΛsinε

+sin²usinicosisinΩcosΛ-sinucosusinisinΩsinΛcosε；

-sin²usinicosicosΩsinΛcosε+sin²ucos²isinΛsinε

步骤三，获取卫星-太阳矢量在星-地-日坐标系O-X_cY_cZ_c中的投影，即卫星-太阳-投影点DV-AOc=LOcOi×DV-AOi;]]>

卫星-太阳-投影点在O-X_cY_cZ_c的X_c轴上的位置分量记为Y_c轴上的位置分量记为Z_c轴上的位置分量记为则DV-AOc=(xOc,yOc,zOc);xOc=E1sinu+E2cosu,]]>yOc=F1+F2sin2u+F3cos2u;]]>

E₁为与轨道要素关联的第一系数，E₂为与轨道要素关联的第二系数，F₁为与轨道要素关联的第三系数，F₂为与轨道要素关联的第四系数，F₃为与轨道要素关联的第五系数；

E₁＝-cosicosΛsinΩ+cosisinΛcosεcosΩ+sinisinΛsinε；

E₂＝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Λ+cos2icos2Λ+2sin2Λsin2icos2ϵ+cos2Λ+;]]>

2sin²Λcos²isin²ε+cos²isin²Λcos²ε+2sin²Λsin²ε+sin²Λcos²ε)

F2=12(cosicos2Λ-cosisin2Λcos2ϵ)sin2Ω-12cosisin2Λcosϵcos2Ω-12sinisin2Λsin2ϵsinΩ-12sinisin2ΛsinϵcosΩ;]]>

F3=-(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ϵ);]]>

步骤四，依据卫星-太阳-投影点在Y_c轴上的位置分量与X_c轴上的位置分量的比值，来表征判定角β的正切角关系

在地球视为均匀球体且卫星轨道的偏心率e最小的情况下，卫星恰好处于进地影位置或者出地影位置，此时的判定角β设为初始常值β₀，且其中π取值为3.14；R_e为地球平均半径，单位为米；a为卫星的轨道半长轴，单位为米；

在卫星轨道确定、以及任意一时间节点t_w的初始常值β₀确定的情况下，卫星恰好进出地影的关系记为

f(utw)=F2sin2u+F3cos2u-tanβ0×E1sinu-tanβ0×E2cosu+F1=0;]]>

步骤五，获取卫星进出地影的时刻；

步骤(5-1)，通过星载计算机的界面输入一个预报时刻t_q，所述预报时刻t_q的形式为年月日时分秒；

步骤(5-2)，判断所述预报时刻t_q是否位于时间节点t_w；

(A)若预报时刻t_q位于时间节点t_w上，即t_q＝t_w，则采用牛顿下山法对进出地影关系f(u)进行迭代运算，并以前一个时间节点t_w-1的进出地影的纬度幅角作为迭代初值，得到预报时刻t_q的进出地影的纬度幅角

对于时间节点tw对应的卫星进出地影的纬度幅角的牛顿下山迭代关系为由于预报时刻t_q位于时间节点t_w上，即t_q＝t_w，能够得到

为时间节点t_w上的进地影纬度幅角迭代值；

为时间节点上的出地影纬度幅角迭代值；

为前一个时间节点t_w-1上的进地影纬度幅角迭代值；

为前一个时间节点t_w-1上的出地影纬度幅角迭代值；

δ为下山因子；

为的进出地影的函数值；

为的导数值；

为的进出地影的函数值；

为的导数值；

(B)若预报时刻t_q不位于时间节点t_w上，且位于轨道参数冻结范围内，则采用轨道参数解析算法对内的时间节点t_w的进出地影的纬度幅角进行计算，得到预报时刻t_q的进出地影的纬度幅角

轨道参数解析算法是指：先选取出属于轨道参数冻结范围内的时间节点t_w；然后计算预报时刻t_q对应的进地影纬度幅角即计算预报时刻t_q对应的出地影纬度幅角即

为t_q处进地影纬度幅角的变化量；

为t_q处出地影纬度幅角的变化量；

ΔE₁＝-(cosicosΛcosΩ+cosisinΛcosεsinΩ)W_ΩΔt；

ΔE₂＝(sinΛcosεcosΩ-cosΛsinΩ)W_ΩΔt；

ΔF1=12(cos2isin2Λcosϵ-sin2Λcosϵ)cos2ΩWΩΔt-12(cos2icos2Λ-cos2isin2Λcos2ϵ-cos2Λ+sin2Λcos2ϵ)sin2ΩWΩΔt+14sin2isin2ΛsinϵcosΩWΩΔt+14sin2isin2Λsin2ϵsinΩWΩΔ;]]>

Δ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ϵcosΩWΩΔt-14sin2isin2Λsin2ϵsinΩWΩΔt;]]>

E₁为与轨道要素关联的第一系数，E₂为与轨道要素关联的第二系数，F₁为与轨道要素关联的第三系数，F₂为与轨道要素关联的第四系数，F₃为与轨道要素关联的第五系数；

ΔE₁为与轨道要素关联的第一系数的变化量，ΔE₂为与轨道要素关联的第二系数的变化量，ΔF₁为与轨道要素关联的第三系数的变化量，ΔF₂为与轨道要素关联的第四系数的变化量，ΔF₃为与轨道要素关联的第五系数的变化量；

Δt为t_q相对于时间节点t_w的时间间隔，则Δt＝t_q-t_w；

W_Ω为升交点赤经Ω的平均变化率，且J₂为地球引力势的二阶谐系数，J₂＝1.0826300×10^-3，R_e为地球平均半径，μ为地球引力常数，i为轨道倾角，a为轨道半长轴，e为轨道偏心率；

步骤(5-3)，利用纬度幅角－时刻关系对进行反解，得到对应的进地影时刻对应的进地影时刻

T_Ω为卫星轨道的交点周期；

t_p为卫星运行的当前时刻；

t_p+1为卫星运行的下一时刻；

为t_p+1时刻所对应的纬度幅角；

为t_p时刻所对应的纬度幅角。

2.根据权利要求1所述的基于低轨道地球卫星的地影时刻预报的星上确定方法，其特征在于：进出地影关系中的纬度幅角u在一个轨道周期t_周期内必定是成对出现的，将进地影纬度幅角记为u_进、出地影纬度幅角记为u_出，则有：