适用于地面测站天线对卫星指向的目标定位方法技术

技术编号:24936084 阅读:30 留言:0更新日期:2020-07-17 20:34
本发明专利技术提供了一种适用于地面测站天线对卫星指向的目标定位方法,通过给定目标时刻的卫星星历数据,经过卫星轨道的相关计算以及相关坐标系的转换计算,得到卫星在地固系下的位置,从而完成卫星的定位过程。本发明专利技术考虑卫星实际运行情况计算地面测站天线与卫星的位置关系,不依赖于仿真软件或过多假设内容,有效解决了用于地面站接收天线对卫星指向的的卫星定位问题,而且达到了比较高的定位精度。

【技术实现步骤摘要】
适用于地面测站天线对卫星指向的目标定位方法
本专利技术涉及卫星轨道计算与控制领域,具体地,涉及一种适用于地面测站天线对卫星指向的目标定位方法。
技术介绍
雷达已应用到制导和超视距探测等多个重要领域,在科技建设中起着非常重要的作用,随着外空间目标探测、控制的需求,地面测站雷达天线的跟踪和搜索能力提出了越来越高的要求,因为星地指向的偏差会导致链路传输信息的信噪比下降,若是超过了最大站位容限,甚至会出现信号丢失的现象。这就要求对卫星的定位要更加快速和精确。因此,设计精度高的卫星定位计算方法具有普遍的实用意义。卫星的定位计算,主要是依据卫星的轨道信息及时间信息,计算出卫星在惯性系下的位置,进而计算卫星在地固系下的位置。国内现有的相关算法研究,大多集中于地面测站位置固定的情况下,在卫星本体坐标系下,进行卫星对地面测站指向的优化设计,而对于地面测站天线对卫星的指向定位的研究较少。本专利技术针对这一实际情况,提出一种用于地面测站天线的精度较高的卫星实时定位方法,实现地面测站天线对卫星指向的目标定位。专利“一种绕月卫星双轴天线对地指向控制方法”(专利号:CN101204994A)介绍了一种计算绕月卫星对地心指向的方法,根据地面上注星历数据推算卫星位置,计算卫星对地球的可见区域,并计算出双轴天线的指向角度。该专利为对地心指向,主要结合月球相关坐标系进行计算。本专利技术与其不同之处在于主要结合地球及地表位置相关坐标系进行计算,完成了地面测站天线对卫星的指向定位计算。专利“一种深空探测器天线指向的设计方法”(专利号:CN105184002A)介绍了一种深空探测器天线对地心指向的方法,此方法是用于实现深空探测器天线对地心的定向。该专利针对天线对地心的定向而并非是地表的给定位置,而且直接给出了探测器天线对地心的指向矢量,没有通过轨道参数计算卫星位置的算法。本专利技术与其不同之处在于,设计了通过给定时刻的卫星轨道参数计算地固系卫星位置的计算过程。
技术实现思路
针对现有技术中的缺陷,本专利技术的目的是提供一种适用于地面测站天线对卫星指向的目标定位方法。根据本专利技术提供的一种适用于地面测站天线对卫星指向的目标定位方法,通过给定目标时刻的卫星星历数据,经过卫星轨道的相关计算以及相关坐标系的转换计算,得到卫星在地固系下的位置,从而完成卫星的定位过程。优选地,包括如下步骤:惯性系卫星位置计算步骤:根据给定的目标时刻t1的星历信息,计算t1时刻卫星在惯性系下的位置RwECI;目标时刻秒计数值计算步骤:根据给定的目标时刻t1,计算历元J2000.0至给定目标时刻的秒计数值tc;地固系卫星位置计算步骤:根据所述计算得到的历元J2000.0至给定目标时刻的秒计数值tc,计算t1时刻卫星在地固系下的位置RwECF。优选地,所述惯性系卫星位置计算步骤中的计算t1时刻卫星在惯性系下的位置RwECI的方法如下:输入目标时刻t1的星历信息,包括:轨道半长轴a、轨道偏心率e、轨道倾角i、升交点赤经Ω、近地点幅角ω以及平近点角M;计算t1时刻卫星在惯性系下的位置RwECI:RwECI=Q*rp其中,旋转矩阵Q按照3-1-3旋转顺序进行描述:矢量rp:其中,M1为真近点角:优选地,所述目标时刻秒计数值计算步骤中的计算历元J2000.0至给定目标时刻的秒计数值tc的方法如下:输入t1时刻的年year、月month、日day、时hour、分min、秒sec,计算儒略日JD:其中,floor()为向下取整运算;根据儒略日JD计算历元J2000.0至给定目标时刻的秒计数值tc:tc=(JD-2455197.5)×86400+315547200优选地,所述地固系卫星位置计算步骤中的计算t1时刻卫星在地固系下的位置RwECF的方法如下:计算地球自转矩阵ER,章动矩阵NR,岁差矩阵PR;根据计算得到的地球自转矩阵ER,章动矩阵NR,岁差矩阵PR,计算惯性系到地固系的转换矩阵MECI2ECF;根据所述计算得到的t1时刻卫星在惯性系下的位置RwECI,计算t1时刻卫星在地固系下的位置RwECF。优选地,所述计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR的计算方法如下:根据所述计算得到的历元J2000.0(至给定目标时刻的秒计数值tc,计算黄经章动Δψ、黄赤交角ε和交角章动Δε:其中,T2k为相对历元J2000.0的儒略世纪数:地球自转矩阵ER计算方法:计算赤经章动Δμ:Δμ=Δψ*cosε计算格林尼治平恒星时计算格林尼治真恒星时SG:计算地球自转矩阵ER:章动矩阵NR计算方法:NR=RX(-ε-Δε)RZ(-Δψ)RX(ε)其中,岁差矩阵PR计算方法:计算岁差常数ζA、θA、ZA:计算岁差矩阵PR:PR=RZ(-ZA)RY(θA)RZ(-ζA)其中,优选地,所述根据计算得到的地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR,计算惯性系到地固系的转换矩阵MECI2ECF的计算方法如下:MECI2ECF=ER*NR*PR优选地,所述根据计算得到的t1时刻卫星在惯性系下的位置RwECI,计算t1时刻卫星在地固系下的位置RwECF的计算方法如下:RwECF=MECI2ECF*RwECI与现有技术相比,本专利技术具有如下的有益效果:本专利技术考虑卫星实际运行情况计算地面测站天线与卫星的位置关系,不依赖于仿真软件或过多假设内容,有效解决了用于地面站接收天线对卫星指向的的卫星定位问题,而且达到了比较高的定位精度。附图说明通过阅读参照以下附图对非限制性实施例所作的详细描述,本专利技术的其它特征、目的和优点将会变得更明显:图1为本专利技术提供的地面测站天线对卫星指向的目标定位方法流程示意图。图2为本专利技术提供的地面测站天线与卫星位置示意图。图3为本专利技术提供的惯性系下卫星位置变化曲线示意图。图4为本专利技术提供的地固系下卫星位置变化曲线示意图。具体实施方式下面结合具体实施例对本专利技术进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本专利技术,但不以任何形式限制本专利技术。应当指出的是,对本领域的普通技术人员来说,在不脱离本专利技术构思的前提下,还可以做出若干变化和改进。这些都属于本专利技术的保护范围。根据本专利技术提供的一种适用于地面测站天线对卫星指向的目标定位方法,通过给定目标时刻的卫星星历数据,经过卫星轨道的相关计算以及相关坐标系的转换计算,得到卫星在地固系下的位置,从而完成卫星的定位过程。具体地,包括如下步骤:惯性系卫星位置计算步骤:根本文档来自技高网
...

【技术保护点】
1.一种适用于地面测站天线对卫星指向的目标定位方法,其特征在于,通过给定目标时刻的卫星星历数据,经过卫星轨道的相关计算以及相关坐标系的转换计算,得到卫星在地固系下的位置,从而完成卫星的定位过程。/n

【技术特征摘要】
1.一种适用于地面测站天线对卫星指向的目标定位方法,其特征在于,通过给定目标时刻的卫星星历数据,经过卫星轨道的相关计算以及相关坐标系的转换计算,得到卫星在地固系下的位置,从而完成卫星的定位过程。


2.根据权利要求1所述的适用于地面测站天线对卫星指向的目标定位方法,其特征在于,包括如下步骤:
惯性系卫星位置计算步骤:根据给定的目标时刻t1的星历信息,计算t1时刻卫星在惯性系下的位置RwECI;
目标时刻秒计数值计算步骤:根据给定的目标时刻t1,计算历元J2000.0至给定目标时刻的秒计数值tc;
地固系卫星位置计算步骤:根据所述计算得到的历元J2000.0至给定目标时刻的秒计数值tc,计算t1时刻卫星在地固系下的位置RwECF。


3.如权利要求2所述的适用于地面测站天线对卫星指向的目标定位方法,其特征在于,所述惯性系卫星位置计算步骤中的计算t1时刻卫星在惯性系下的位置RwECI的方法如下:
输入目标时刻t1的星历信息,包括:轨道半长轴a、轨道偏心率e、轨道倾角i、升交点赤经Ω、近地点幅角ω以及平近点角M;
计算t1时刻卫星在惯性系下的位置RwECI:
RwECI=Q*rp
其中,旋转矩阵Q按照3-1-3旋转顺序进行描述:



矢量rp:



其中,M1为真近点角:




4.如权利要求2所述的适用于地面测站天线对卫星指向的目标定位方法,其特征在于,所述目标时刻秒计数值计算步骤中的计算历元J2000.0至给定目标时刻的秒计数值tc的方法如下:
输入t1时刻的年year、月month、日day、时hour、分min、秒sec,计算儒略日JD:



其中,floor()为向下取整运算;
根据儒略日JD计算历元J2000.0至给定目标时刻的秒计数值tc:
tc=(JD-2455197.5)×86400+315547200。


5.如权利要求2所述的适用于地面测站天线对卫星指向的目标定位方法,其特征在于,所述地固系卫星位置计算步骤中的计算t1时刻卫星在地固系下的位置RwECF的方法如下:<...

【专利技术属性】
技术研发人员:吕旺郑峰司力琼顾强俞航王田野徐晔杨珺
申请(专利权)人:上海卫星工程研究所
类型:发明
国别省市:上海;31

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1