[发明专利]一种下行链路最差干扰快速估算方法

权利要求书

1.一种下行链路最差干扰快速估算方法，其特征在于，包括以下步骤：

S101：初始化NGSO卫星的纬度幅角；

S102：计算NGSO卫星在J2000坐标系的三维位置矢量、赤经和J2000坐标系到轨道坐标系的坐标转换矩阵；

S103：初始化GSO卫星的赤经；

S104：计算GSO卫星在J2000坐标系的位置矢量；

S105：计算GSO卫星到NGSO卫星的延长线在轨道坐标系的单位矢量；

S106：判断GSO卫星到NGSO卫星的延长线是否与地球表面相交，如果相交，转到S107，不相交转到S111；

S107：计算GSO卫星到NGSO卫星的延长线与地球表面的交点在J2000坐标系的矢量表达式；

S108：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线方向的发射增益；

S109：计算NGSO卫星空间电台的发射在GSO卫星到NGSO卫星的延长线与地球表面的交点处GSO地球站的EPFD；

S110：存储EPFD，以及对应NGSO卫星纬度幅角、GSO卫星赤经以及GSO卫星到NGSO卫星的延长线与地球表面的交点的坐标；

S111：更新GSO卫星的赤经，经判断后转到S104，直至完成所有赤经的计算；

S112：更新NGSO卫星的纬度幅角，经判断后转到S102，直至完成所有纬度幅角的计算；

S113：完成所述计算，获得所有NGSO卫星纬度幅角、所有GSO卫星赤经对应的所有EPFD值，根据EPFD限值格式，每种场景都存在一个“100％”不能超过的限值，即完全不能超过的值，进一步计算所有EPFD值的最大值，即为最差干扰值。

2.根据权利要求1所述的方法，其特征在于，所述初始化NGSO卫星的纬度幅角包括：

S1011：将NGSO卫星的纬度幅角取为u＝u₀，其中，u₀为NGSO卫星的初始纬度幅角u₀＝0°。

3.根据权利要求2所述的方法，其特征在于，所述计算NGSO卫星在J2000坐标系的位置矢量、赤经以及J2000坐标系到轨道坐标系的坐标转换矩阵包括：

S1021：根据纬度幅角u和近地点幅角ω，计算真近点角：

θ＝u-ω

S1022：根据所述真近点角计算NGSO卫星轨道半径为：

其中，a是NGSO卫星半长轴，e是NGSO卫星偏心率；

S1023：NGSO卫星在J2000坐标系的位置矢量为：

其中，为单位向量，计算公式为：

其中，Ω是NGSO卫星升交点赤经，i是NGSO卫星轨道倾角；

S1024：将NGSO卫星在J2000坐标系的位置矢量的三维分量记为x_N、y_N和z_N，NGSO卫星在J2000坐标系的赤经为：

S1025：J2000坐标系到轨道坐标系的坐标转换矩阵为：

其中，M_x()、M_y()和M_z()分别为绕X、Y和Z轴的坐标转换矩阵计算公式为：

4.根据权利要求3所述的方法，其特征在于，所述初始化GSO卫星的赤经包括：

S1031：初始化GSO卫星的赤经：

L_G＝L_N-ΔL_G

其中，ΔL_G为GSO卫星赤经区间宽度的一半。

5.根据权利要求4所述的方法，其特征在于，所述计算GSO卫星在J2000坐标系的位置矢量包括：

S1041：GSO卫星在J2000坐标系的位置矢量r_G为：

其中，R_e是地球平均半径，B_G为GSO卫星的赤纬，H_G为GSO卫星的高度。