CN111427004A - 适用于地面测站天线对卫星指向的坐标转换方法 - Google Patents

适用于地面测站天线对卫星指向的坐标转换方法 Download PDF

Info

Publication number
CN111427004A
CN111427004A CN202010198203.1A CN202010198203A CN111427004A CN 111427004 A CN111427004 A CN 111427004A CN 202010198203 A CN202010198203 A CN 202010198203A CN 111427004 A CN111427004 A CN 111427004A
Authority
CN
China
Prior art keywords
matrix
calculating
ground
satellite
conversion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010198203.1A
Other languages
English (en)
Inventor
俞航
李秀伟
徐晔
吕旺
陆启省
董瑶海
黄欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN202010198203.1A priority Critical patent/CN111427004A/zh
Publication of CN111427004A publication Critical patent/CN111427004A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
    • G01S1/02—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using radio waves
    • G01S1/08—Systems for determining direction or position line

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明提供了一种适用于地面测站天线对卫星指向的坐标转换方法，通过给定目标时刻和地面测站天线的地理位置信息，计算得到惯性系到地固系的转换矩阵和地固系到站心系的转换矩阵，用于地面测站天线对卫星指向相关的坐标转换计算。本发明不依赖于仿真软件或过多假设内容，对地球实际在轨运行的因素考虑比较全面，有效解决了地面测站天线对卫星指向控制相关坐标系转换的计算问题，而且达到了比较高的转换精度。

Description

适用于地面测站天线对卫星指向的坐标转换方法
技术领域
本发明涉及卫星轨道控制及计算领域，具体地，涉及一种适用于地面测站天线对卫星指向的坐标转换方法。
背景技术
由于卫星信号微弱而且方向性强，为了捕获运动卫星上的通讯信号，必须实时调整天线姿态与卫星的位置偏差以满足通讯需求，因而对天线的跟踪和搜索能力提出了越来越高的要求。这就要求雷达天线必须要根据指令调整指向，实时跟踪运动目标。因此，雷达天线指向过程的动态精度已经成为天线系统功能的重要指标之一，设计高精度的计算方法对天线指向过程相关坐标系之间的转换关系进行计算具有普遍的实用意义。
国内现有的星地指向算法研究，大多集中于地面测站位置固定的情况下，在卫星本体坐标系下，进行卫星对地面测站指向的优化设计，而对于地面测站天线对卫星的指向导引的研究较少。本发明针对这一实际情况，提出一种适用于地面测站天线对卫星指向的坐标转换方法，完成天线指向过程相关坐标系之间的转换关系的计算。
专利“一种数传天线指向角度的仿真分析方法”(专利号：CN105184002A)介绍了一种计算星载数传天线对地面站指向的方法，此方法通过已有的卫星轨道仿真软件STK对卫星的实际位置进行仿真求解，并计算数传天线的二维指向角。该专利的不足在于对坐标系转换的描述比较简略，没有给出转换矩阵的算法。本发明的优势在于设计了一套相关坐标系转换矩阵的详细计算流程。
文献“卫星锐波束天线指向算法及仿真”(见《中国空间科学技术》，2008年，2期)介绍了一种卫星星载锐波束天线对地表目标点指向的算法，该方法的不足之处在于，在坐标系转换的过程中忽略的影响因素较多，没有考虑类似岁差、章动等影响，指向精度较低。本发明对这些可能误差因素造成的影响进行了公式推导，使坐标系转换过程更为精确。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种适用于地面测站天线对卫星指向的坐标转换方法。
根据本发明提供的一种适用于地面测站天线对卫星指向的坐标转换方法，通过给定目标时刻和地面测站天线的地理位置信息，计算得到惯性系到地固系的转换矩阵和地固系到站心系的转换矩阵，用于地面测站天线对卫星指向相关的坐标转换计算。
优选地，包括如下步骤：
目标时刻秒计数值计算步骤：根据给定的目标时刻t1，计算历元J2000.0至给定目标时刻的秒计数值tc；
惯性系到地固系转换矩阵计算步骤：根据所述计算出的历元J2000.0至给定目标时刻的秒计数值tc，计算t1时刻惯性系到地固系的转换矩阵MECI2ECF；
地固系到站心系转换矩阵计算步骤：根据给定的地面测站天线的经纬度，计算地固系到站心系的转换矩阵MECF2CT；
优选地，所述目标时刻秒计数值计算步骤中的计算历元J2000.0至给定目标时刻的秒计数值tc的方法如下：
给定目标时刻t1，计算历元J2000.0至给定目标时刻的秒计数值tc，输入t1时刻的年year、月month、日day、时hour、分min、秒sec，计算儒略日JD：
Figure BDA0002418376880000021
其中，floor()为向下取整运算；
根据儒略日JD计算历元J2000.0至给定目标时刻的秒计数值tc：
tc＝(JD-2455197.5)×86400+315547200
优选地，所述惯性系到地固系转换矩阵计算步骤中的计算t1时刻惯性系到地固系的转换矩阵MECI2ECF的方法如下：
根据所述历元J2000.0至给定目标时刻的秒计数值tc，计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR；
再根据计算的地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR，计算惯性系到地固系的转换矩阵MECI2ECF。
优选地，所述计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR的计算方法如下：
首先计算黄经章动Δψ、黄赤交角ε和交角章动Δε：
Figure BDA0002418376880000031
Figure BDA0002418376880000032
Figure BDA0002418376880000033
其中，T2k为相对历元J2000.0的儒略世纪数：
Figure BDA0002418376880000034
地球自转矩阵ER计算方法如下：
计算赤经章动Δμ：
Δμ＝Δψ*cosε
计算格林尼治平恒星时
Figure BDA0002418376880000035
Figure BDA0002418376880000036
计算格林尼治真恒星时SG：
Figure BDA0002418376880000037
计算地球自转矩阵ER：
Figure BDA0002418376880000038
章动矩阵NR计算方法如下：
NR＝RX(-ε-Δε)RZ(-Δψ)RX(ε)
其中，
Figure BDA0002418376880000039
Figure BDA0002418376880000041
Figure BDA0002418376880000042
岁差矩阵PR计算方法如下：
计算岁差常数ζA、θA、ZA：
Figure BDA0002418376880000043
计算岁差矩阵PR：
PR＝RZ(-ZA)RY(θA)RZ(-ζA)
其中，
Figure BDA0002418376880000044
Figure BDA0002418376880000045
Figure BDA0002418376880000046
优选地，所述根据地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR，计算惯性系到地固系的转换矩阵MECI2ECF的计算方法如下：
MECI2ECF＝ER*NR*PR
优选地，所述地固系到站心系转换矩阵计算步骤中的计算地固系到站心系的转换矩阵MECF2CT的方法如下：
在站心系下，计算地固系到站心系的转换矩阵MECF2CT，描述为在地固系下一次绕Z轴的旋转和一次绕X轴的旋转：
MECF2CT＝Rx(90°-lat)Rz(90°+lon)
其中，
Figure BDA0002418376880000051
Figure BDA0002418376880000052
并根据同一坐标系下卫星与地面测站天线的相对位置关系进行一次坐标原点的平移，即完成坐标系的转换过程，可以表示为：
RwCT＝MECF2CT*(RwECF-RtECF)
其中，
RwECF为卫星在地固系下的位置；RtECF为地面测站天线在地固系下的位置；RwCT为卫星在站心系下的位置。
与现有技术相比，本发明具有如下的有益效果：
本发明不依赖于仿真软件或过多假设内容，对地球实际在轨运行的因素考虑比较全面，有效解决了地面测站天线对卫星指向控制相关坐标系转换的计算问题，而且达到了比较高的转换精度。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明提供的适用于地面测站天线对卫星指向的坐标转换方法流程示意图。
图2为本发明提供的站心系OCTXCTYCTZCT示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
根据本发明提供的一种适用于地面测站天线对卫星指向的坐标转换方法，通过给定目标时刻和地面测站天线的地理位置信息，计算得到惯性系到地固系的转换矩阵和地固系到站心系的转换矩阵，用于地面测站天线对卫星指向相关的坐标转换计算。
具体地，包括如下步骤：
目标时刻秒计数值计算步骤：根据给定的目标时刻t1，计算历元J2000.0至给定目标时刻的秒计数值tc；
惯性系到地固系转换矩阵计算步骤：根据所述计算出的历元J2000.0至给定目标时刻的秒计数值tc，计算t1时刻惯性系到地固系的转换矩阵MECI2ECF；
地固系到站心系转换矩阵计算步骤：根据给定的地面测站天线的经纬度，计算地固系到站心系的转换矩阵MECF2CT；
具体地，所述目标时刻秒计数值计算步骤中的计算历元J2000.0至给定目标时刻的秒计数值tc的方法如下：
给定目标时刻t1，计算历元J2000.0至给定目标时刻的秒计数值tc，输入t1时刻的年year、月month、日day、时hour、分min、秒sec，计算儒略日JD：
Figure BDA0002418376880000061
其中，floor()为向下取整运算；
根据儒略日JD计算历元J2000.0至给定目标时刻的秒计数值tc：
tc＝(JD-2455197.5)×86400+315547200
具体地，所述惯性系到地固系转换矩阵计算步骤中的计算t1时刻惯性系到地固系的转换矩阵MECI2ECF的方法如下：
根据所述历元J2000.0至给定目标时刻的秒计数值tc，计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR；
再根据计算的地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR，计算惯性系到地固系的转换矩阵MECI2ECF。
具体地，所述计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR的计算方法如下：
首先计算黄经章动Δψ、黄赤交角ε和交角章动Δε：
Figure BDA0002418376880000071
Figure BDA0002418376880000072
Figure BDA0002418376880000073
其中，T2k为相对历元J2000.0的儒略世纪数：
Figure BDA0002418376880000074
地球自转矩阵ER计算方法如下：
计算赤经章动Δμ：
Δμ＝Δψ*cosε
计算格林尼治平恒星时
Figure BDA0002418376880000075
Figure BDA0002418376880000076
计算格林尼治真恒星时SG：
Figure BDA0002418376880000077
计算地球自转矩阵ER：
Figure BDA0002418376880000078
章动矩阵NR计算方法如下：
NR＝RX(-ε-Δε)RZ(-Δψ)RX(ε)
其中，
Figure BDA0002418376880000079
Figure BDA00024183768800000710
Figure BDA00024183768800000711
岁差矩阵PR计算方法如下：
计算岁差常数ζA、θA、ZA：
Figure BDA0002418376880000081
计算岁差矩阵PR：
PR＝RZ(-ZA)RY(θA)RZ(-ζA)
其中，
Figure BDA0002418376880000082
Figure BDA0002418376880000083
Figure BDA0002418376880000084
具体地，所述根据地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR，计算惯性系到地固系的转换矩阵MECI2ECF的计算方法如下：
MECI2ECF＝ER*NR*PR
具体地，所述地固系到站心系转换矩阵计算步骤中的计算地固系到站心系的转换矩阵MECF2CT的方法如下：
在站心系下，计算地固系到站心系的转换矩阵MECF2CT，描述为在地固系下一次绕Z轴的旋转和一次绕X轴的旋转：
MECF2CT＝Rx(90°-lat)Rz(90°+lon)
其中，
Figure BDA0002418376880000085
Figure BDA0002418376880000091
并根据同一坐标系下卫星与地面测站天线的相对位置关系进行一次坐标原点的平移，即完成坐标系的转换过程，可以表示为：
RwCT＝MECF2CT*(RwECF-RtECF)
其中，
RwECF为卫星在地固系下的位置；RtECF为地面测站天线在地固系下的位置；RwCT为卫星在站心系下的位置。
下面通过优选例，对本发明进行更为具体地说明。
优选例1：
本发明需要用到的坐标系：所述的惯性系为J2000.0惯性坐标系，所述的地固系为WGS-84坐标系。下面给出站心系的定义。
站心系OCTXCTYCTZCT
站心系的定义为，原点OCT为地面天线原点，基本平面OCTXCTYCT面为当地水平面，OCTXCT沿当地子午圈指向正北，OCTZCT垂直基本平面指向天顶，OCTYCT按右手法则确定，如图2所示。
下面详述本发明的计算过程：
利用MATLAB对此算法进行仿真验证，地球相关参数以及站心系按上文所述设置，某型号卫星在UTC时间2018年12月3日5时30分的星历数据如下：
Figure BDA0002418376880000092
Figure BDA0002418376880000101
(1)根据给定的目标时刻t1，计算历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，具体步骤如下：
输入t1时刻(UTC时间)的年(year)、月(month)、日(day)、时(hour)、分(min)、秒(sec)，计算儒略日JD：
Figure BDA0002418376880000102
其中，floor()为向下取整运算。
根据儒略日JD计算历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc：
tc＝(JD-2455197.5)×86400+315547200
计算结果为：
tc＝597087000
(2)根据所述步骤(1)计算出的历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，计算t1时刻惯性系到地固系的转换矩阵MECI2ECF，具体步骤如下：
根据所述步骤(1)计算得到的历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，计算地球自转矩阵ER，章动矩阵NR，岁差矩阵PR，由于极移对转换矩阵的计算影响很小，本发明中不考虑此项。
首先计算黄经章动Δψ、黄赤交角ε和交角章动Δε：
Figure BDA0002418376880000103
Figure BDA0002418376880000104
Figure BDA0002418376880000111
其中，T2k为相对历元J2000.0(2000年1月1日12时)的儒略世纪数：
Figure BDA0002418376880000112
地球自转矩阵ER计算方法：
计算赤经章动Δμ：
Δμ＝Δψ*cosε
计算格林尼治平恒星时
Figure BDA0002418376880000113
Figure BDA0002418376880000114
计算格林尼治真恒星时SG：
Figure BDA0002418376880000115
计算地球自转矩阵ER：
Figure BDA0002418376880000116
章动矩阵NR计算方法：
NR＝RX(-ε-Δε)RZ(-Δψ)RX(ε)
其中，
Figure BDA0002418376880000117
Figure BDA0002418376880000118
Figure BDA0002418376880000119
岁差矩阵PR计算方法：
计算岁差常数ζA、θA、ZA：
Figure BDA0002418376880000121
计算岁差矩阵PR：
PR＝RZ(-ZA)RY(θA)RZ(-ζA)
其中，
Figure BDA0002418376880000122
Figure BDA0002418376880000123
Figure BDA0002418376880000124
根据地球自转矩阵ER、章动矩阵NR、岁差矩阵PR、计算惯性系到地固系的转换矩阵MECI2ECF：
MECI2ECF＝ER*NR*PR
计算结果为：
Figure BDA0002418376880000125
(3)根据给定的地面测站天线的经纬度，计算地固系到站心系的转换矩阵MECF2CT，具体步骤如下：
在站心系下，计算地固系到站心系的转换矩阵MECF2CT，描述为在地固系下一次绕Z轴的旋转和一次绕X轴的旋转：
MECF2CT＝Rx(90°-lat)Rz(90°+lon)
其中，
Figure BDA0002418376880000126
Figure BDA0002418376880000131
计算结果为：
Figure BDA0002418376880000132
并根据同一坐标系下卫星与地面测站天线的相对位置关系进行一次坐标原点的平移，即完成坐标系的转换过程，可以表示为：
RwCT＝MECF2CT*(RwECF-RtECF)
其中，RwECF为卫星在地固系下的位置；RtECF为地面测站天线在地固系下的位置；RwCT为卫星在站心系下的位置。
优选例2：
本发明的技术解决问题是：克服现有技术的不足，提供一种精度较高的适用于地面测站天线对卫星指向的坐标转换方法，通过给定目标时刻和地面测站天线的地理位置信息，计算得到惯性系-地固系的转换矩阵和地固系-站心系的转换矩阵，可以用于地面测站天线对卫星指向相关的坐标转换计算。
本发明提出的适用于地面测站天线对卫星指向的坐标转换方法，对坐标系转换关系的影响因素考虑全面，计算精度较高，有效解决了地面测站天线实时对卫星指向过程中相关坐标系转换计算的需求。
如图1所示，本发明提供的适用于地面测站天线对卫星指向的坐标转换方法流程示意图，本发明的技术解决方案的步骤如下：
(1)根据给定的目标时刻t1，计算历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，具体过程如下：
给定目标时刻t1，计算历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，输入t1时刻(UTC时间)的年(year)、月(month)、日(day)、时(hour)、分(min)、秒(sec)，计算儒略日JD：
Figure BDA0002418376880000141
其中，floor()为向下取整运算。
根据儒略日JD计算历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc：
tc＝(JD-2455197.5)×86400+315547200
(2)根据所述步骤(1)计算出的历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，计算t1时刻惯性系到地固系的转换矩阵MECI2ECF，具体的过程如下：
根据所述步骤(1)计算得到的历元J2000.0(2000年1月1日12时)至给定目标时刻的秒计数值tc，计算地球自转矩阵ER，章动矩阵NR，岁差矩阵PR，由于极移对转换矩阵的计算影响很小，本发明中不考虑此项。
首先计算黄经章动Δψ、黄赤交角ε和交角章动Δε：
Figure BDA0002418376880000142
Figure BDA0002418376880000143
Figure BDA0002418376880000144
其中，T2k为相对历元J2000.0(2000年1月1日12时)的儒略世纪数：
Figure BDA0002418376880000145
地球自转矩阵ER计算方法：
计算赤经章动Δμ：
Δμ＝Δψ*cosε
计算格林尼治平恒星时
Figure BDA0002418376880000146
Figure BDA0002418376880000147
计算格林尼治真恒星时SG：
Figure BDA0002418376880000151
计算地球自转矩阵ER：
Figure BDA0002418376880000152
章动矩阵NR计算方法：
NR＝RX(-ε-Δε)RZ(-Δψ)RX(ε)
其中，
Figure BDA0002418376880000153
Figure BDA0002418376880000154
Figure BDA0002418376880000155
岁差矩阵PR计算方法：
计算岁差常数ζA、θA、ZA：
Figure BDA0002418376880000156
计算岁差矩阵PR：
PR＝RZ(-ZA)RY(θA)RZ(-ζA)
其中，
Figure BDA0002418376880000157
Figure BDA0002418376880000158
Figure BDA0002418376880000161
根据地球自转矩阵ER、章动矩阵NR、岁差矩阵PR、计算惯性系到地固系的转换矩阵MECI2ECF：
MECI2ECF＝ER*NR*PR
(3)根据给定的地面测站天线的经纬度，计算地固系到站心系的转换矩阵MECF2CT，具体的过程如下：
在站心系下，计算地固系到站心系的转换矩阵MECF2CT，描述为在地固系下一次绕Z轴的旋转和一次绕X轴的旋转：
MECF2CT＝Rx(90°-lat)Rz(90°+lon)
其中，
Figure BDA0002418376880000162
Figure BDA0002418376880000163
并根据同一坐标系下卫星与地面测站天线的相对位置关系进行一次坐标原点的平移，即完成坐标系的转换过程，可以表示为：
RwCT＝MECF2CT*(RwECF-RtECF)
其中，RwECF为卫星在地固系下的位置；RtECF为地面测站天线在地固系下的位置；RwCT为卫星在站心系下的位置。
在本申请的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本申请和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本申请的限制。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (7)

1.一种适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，通过给定目标时刻和地面测站天线的地理位置信息，计算得到惯性系到地固系的转换矩阵和地固系到站心系的转换矩阵，用于地面测站天线对卫星指向相关的坐标转换计算。
2.根据权利要求1所述的适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，包括如下步骤：
目标时刻秒计数值计算步骤：根据给定的目标时刻t1，计算历元J2000.0至给定目标时刻的秒计数值tc；
惯性系到地固系转换矩阵计算步骤：根据所述计算出的历元J2000.0至给定目标时刻的秒计数值tc，计算t1时刻惯性系到地固系的转换矩阵MECI2ECF；
地固系到站心系转换矩阵计算步骤：根据给定的地面测站天线的经纬度，计算地固系到站心系的转换矩阵MECF2CT。
3.如权利要求2所述的适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，所述目标时刻秒计数值计算步骤中的计算历元J2000.0至给定目标时刻的秒计数值tc的方法如下：
给定目标时刻t1，计算历元J2000.0至给定目标时刻的秒计数值tc，输入t1时刻的年year、月month、日day、时hour、分min、秒sec，计算儒略日JD：
Figure FDA0002418376870000011
其中，floor()为向下取整运算；
根据儒略日JD计算历元J2000.0至给定目标时刻的秒计数值tc：
tc＝(JD-2455197.5)×86400+315547200。
4.如权利要求2所述的适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，所述惯性系到地固系转换矩阵计算步骤中的计算t1时刻惯性系到地固系的转换矩阵MECI2ECF的方法如下：
根据所述历元J2000.0至给定目标时刻的秒计数值tc，计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR；
再根据计算的地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR，计算惯性系到地固系的转换矩阵MECI2ECF。
5.如权利要求4所述的适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，所述计算地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR的计算方法如下：
首先计算黄经章动Δψ、黄赤交角ε和交角章动Δε：
Figure FDA0002418376870000021
Figure FDA0002418376870000022
Figure FDA0002418376870000023
其中，T2k为相对历元J2000.0的儒略世纪数：
Figure FDA0002418376870000024
地球自转矩阵ER计算方法如下：
计算赤经章动Δμ：
Δμ＝Δψ*cosε
计算格林尼治平恒星时
Figure FDA0002418376870000025
Figure FDA0002418376870000026
计算格林尼治真恒星时SG：
Figure FDA0002418376870000027
计算地球自转矩阵ER：
Figure FDA0002418376870000028
章动矩阵NR计算方法如下：
NR＝RX(-ε-Δε)RZ(-Δψ)RX(ε)
其中，
Figure FDA0002418376870000031
Figure FDA0002418376870000032
Figure FDA0002418376870000033
岁差矩阵PR计算方法如下：
计算岁差常数ζA、θA、ZA：
Figure FDA0002418376870000034
计算岁差矩阵PR：
PR＝RZ(-ZA)RY(θA)RZ(-ζA)
其中，
Figure FDA0002418376870000035
Figure FDA0002418376870000036
Figure FDA0002418376870000037
6.如权利要求4所述的适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，所述根据地球自转矩阵ER、章动矩阵NR以及岁差矩阵PR，计算惯性系到地固系的转换矩阵MECI2ECF的计算方法如下：
MECI2ECF＝ER*NR*PR。
7.如权利要求2所述的适用于地面测站天线对卫星指向的坐标转换方法，其特征在于，所述地固系到站心系转换矩阵计算步骤中的计算地固系到站心系的转换矩阵MECF2CT的方法如下：
在站心系下，计算地固系到站心系的转换矩阵MECF2CT，描述为在地固系下一次绕Z轴的旋转和一次绕X轴的旋转：
MECF2CT＝Rx(90°-lat)Rz(90°+lon)
其中，
Figure FDA0002418376870000041
Figure FDA0002418376870000042
并根据同一坐标系下卫星与地面测站天线的相对位置关系进行一次坐标原点的平移，即完成坐标系的转换过程，可以表示为：
RwCT＝MECF2CT*(RwECF-RtECF)
其中，
RwECF为卫星在地固系下的位置；RtECF为地面测站天线在地固系下的位置；RwCT为卫星在站心系下的位置。
CN202010198203.1A 2020-03-19 2020-03-19 适用于地面测站天线对卫星指向的坐标转换方法 Pending CN111427004A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010198203.1A CN111427004A (zh) 2020-03-19 2020-03-19 适用于地面测站天线对卫星指向的坐标转换方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010198203.1A CN111427004A (zh) 2020-03-19 2020-03-19 适用于地面测站天线对卫星指向的坐标转换方法

Publications (1)

Publication Number Publication Date
CN111427004A true CN111427004A (zh) 2020-07-17

Family

ID=71549839

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010198203.1A Pending CN111427004A (zh) 2020-03-19 2020-03-19 适用于地面测站天线对卫星指向的坐标转换方法

Country Status (1)

Country Link
CN (1) CN111427004A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109655064A (zh) * 2018-11-27 2019-04-19 孙秀聪 一种基于相对转动的地球固连系-惯性系快速高精度转换方法
CN112329202A (zh) * 2020-09-30 2021-02-05 北京空间飞行器总体设计部 一种火星车对环绕器天线指向算法的优化实现方法
CN115859839A (zh) * 2023-02-23 2023-03-28 中国人民解放军32035部队 基于rcs序列阶跃效应的抛物面天线载荷指向估计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6085128A (en) * 1998-02-06 2000-07-04 The United States Of America As Represented By The Secretary Of The Navy Orbit/covariance estimation and analysis (OCEAN) determination for satellites
CN104729457A (zh) * 2015-04-16 2015-06-24 哈尔滨工业大学 太阳相对近地轨道微小卫星位置的确定方法
CN104777843A (zh) * 2015-04-03 2015-07-15 上海微小卫星工程中心 一种空间飞行器对地面站高精度指向控制方法
CN105158777A (zh) * 2015-07-31 2015-12-16 上海卫星工程研究所 用于无源测向定位的数据源生成方法
CN108974395A (zh) * 2018-06-21 2018-12-11 中国人民解放军战略支援部队航天工程大学 基于天基激光平台驱动的空间目标变轨计算方法及其装置
CN110838864A (zh) * 2018-08-19 2020-02-25 南京理工大学 无人值守的卫星地面站跟踪控制系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6085128A (en) * 1998-02-06 2000-07-04 The United States Of America As Represented By The Secretary Of The Navy Orbit/covariance estimation and analysis (OCEAN) determination for satellites
CN104777843A (zh) * 2015-04-03 2015-07-15 上海微小卫星工程中心 一种空间飞行器对地面站高精度指向控制方法
CN104729457A (zh) * 2015-04-16 2015-06-24 哈尔滨工业大学 太阳相对近地轨道微小卫星位置的确定方法
CN105158777A (zh) * 2015-07-31 2015-12-16 上海卫星工程研究所 用于无源测向定位的数据源生成方法
CN108974395A (zh) * 2018-06-21 2018-12-11 中国人民解放军战略支援部队航天工程大学 基于天基激光平台驱动的空间目标变轨计算方法及其装置
CN110838864A (zh) * 2018-08-19 2020-02-25 南京理工大学 无人值守的卫星地面站跟踪控制系统

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109655064A (zh) * 2018-11-27 2019-04-19 孙秀聪 一种基于相对转动的地球固连系-惯性系快速高精度转换方法
CN109655064B (zh) * 2018-11-27 2021-04-27 孙秀聪 一种基于相对转动的地球固连系-惯性系快速高精度转换方法
CN112329202A (zh) * 2020-09-30 2021-02-05 北京空间飞行器总体设计部 一种火星车对环绕器天线指向算法的优化实现方法
CN112329202B (zh) * 2020-09-30 2023-08-15 北京空间飞行器总体设计部 一种火星车对环绕器天线指向算法的优化实现方法
CN115859839A (zh) * 2023-02-23 2023-03-28 中国人民解放军32035部队 基于rcs序列阶跃效应的抛物面天线载荷指向估计方法

Similar Documents

Publication Publication Date Title
CN111427002B (zh) 地面测控天线指向卫星的方位角计算方法
CN111427004A (zh) 适用于地面测站天线对卫星指向的坐标转换方法
Petrovski et al. Digital satellite navigation and geophysics: A practical guide with GNSS signal simulator and receiver laboratory
CN102736091B (zh) 星球表面广域探测的卫星导航方法及系统
CN105184002B (zh) 一种数传天线指向角度的仿真分析方法
CN111381256B (zh) 主动遥感卫星天线相位中心偏移误差计算的方法和系统
CN111427001A (zh) 适用于地面测站天线对卫星指向的目标定位方法
CN111427003A (zh) 地面测站天线对卫星的指向导引系统
CN103675773A (zh) 一种定标器与卫星指向对准的确定方法
CN102426017A (zh) 一种基于星敏感器确定载体相对于地理坐标系姿态的方法
CN112966211B (zh) 卫星对目标观测下视角计算方法及系统
CN105160125A (zh) 一种星敏感器四元数的仿真分析方法
Lorell et al. Gravity field of Mars from Mariner 9 tracking data
Weiss et al. Orbit and clock product generation
CN111679242A (zh) 适用于指向在轨航天器的地面天线导引方法
CN112762925A (zh) 一种基于地磁计和陀螺仪的低轨卫星定姿方法
CN111366126A (zh) 地面测站天线对卫星指向的视向量计算系统
CN111427000A (zh) 适用于地面测站天线对卫星指向的目标视向量确定方法
Jin et al. Sentinel-6A attitude modeling with dual GNSS antennas and its impact on precise orbit determination
CN105891851B (zh) 一种基于导航卫星漏信号进行定位在轨试验验证方法
Gustavsson Development of a MatLab based GPS constellation simulation for navigation algorithm developments
CN112607056B (zh) 雷达卫星对目标观测自主启动触发方法和系统
Zhang et al. 3D digital track map-based GNSS NLOS signal analytical identification method
Tsai et al. Hardware-in-the-loop Validation of GPS/GNSS Based Mission Planning for LEO Satellites
JPH1172718A (ja) 天体望遠鏡

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20200717

RJ01 Rejection of invention patent application after publication