CN113158341A - 航天器中继终端天线跟踪预报方法、装置和存储介质 - Google Patents
航天器中继终端天线跟踪预报方法、装置和存储介质 Download PDFInfo
- Publication number
- CN113158341A CN113158341A CN202110417513.2A CN202110417513A CN113158341A CN 113158341 A CN113158341 A CN 113158341A CN 202110417513 A CN202110417513 A CN 202110417513A CN 113158341 A CN113158341 A CN 113158341A
- Authority
- CN
- China
- Prior art keywords
- target
- angle
- relay terminal
- terminal antenna
- relay
- 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
Links
- 238000013277 forecasting method Methods 0.000 title claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 39
- 239000013598 vector Substances 0.000 claims description 34
- 230000000295 complement effect Effects 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 description 35
- 238000010586 diagram Methods 0.000 description 18
- 239000011159 matrix material Substances 0.000 description 14
- 238000004891 communication Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 230000009466 transformation Effects 0.000 description 7
- 206010019345 Heat stroke Diseases 0.000 description 5
- 208000007180 Sunstroke Diseases 0.000 description 5
- 230000036544 posture Effects 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000005034 decoration Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 210000001503 joint Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241001183967 Isodon Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computing Systems (AREA)
- Radio Relay Systems (AREA)
Abstract
本发明公开了一种航天器中继终端天线跟踪预报方法、装置和存储介质。其中,该方法包括:采集航天器的形状信息;基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上;获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件;基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。本发明解决了由于无法确定中继终端天线与中继卫星之间是否遮挡,导致中继终端天线跟踪预报准确率低的技术问题。
Description
技术领域
本发明涉及航天领域,具体而言,涉及一种航天器中继终端天线跟踪预报方法、装置和存储介质。
背景技术
目前,在我国载人航天工程任务中,航天器中继跟踪测量是最重要的天基测控手段之一,在航天器飞控任务规划与设计时,需要给出航天器所安装的中继终端天线受遮挡区域的定量计算结果作为飞控任务测控弧段设计和在轨操作的输入。而由于大型航天器组合体姿态种类更多、构型更加复杂,尤其是对于载人飞船径向对接的三目标或四目标组合体惯性飞行情况以及在径向交会对接过程中,会存在中继终端天线被大型舱体以及太阳帆板遮挡的问题,而现有的航天器中继跟踪预报方法未考虑到这一遮挡问题,导致中继终端天线对中继卫星的跟踪预报准确率低的问题。
针对上述遮挡影响下的中继终端天线跟踪预报准确率低的问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种航天器中继终端天线跟踪预报方法、装置和存储介质,以至少解决由于无法确定中继终端天线与中继卫星之间是否遮挡,导致中继终端天线跟踪预报准确率低的技术问题。
根据本发明实施例的一个方面,提供了一种航天器中继终端天线跟踪预报方法,包括:采集航天器的形状信息;基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上;获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件;基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。
可选地,获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值确定中继终端天线与地心二者之间的第一目标距离;基于第二目标坐标值确定中继卫星与地心二者之间的第二目标距离;基于第一目标距离与第二目标距离确定第一目标角度的角度值和第二目标角度的角度值,其中,第一目标角度为第一目标距离和第二目标距离分别与同一地球半径的夹角和,第二目标角度为第一距离与第二距离之间的夹角;基于第一目标角度的角度值和第二目标角度的角度值确定跟踪条件。
可选地,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在第一目标角度的角度值大于第二目标角度的角度值的情况下,确定中继终端天线与中继卫星未发生遮挡;在第一目标角度的角度值小于等于第二目标角度的角度值的情况下,确定中继终端天线与中继卫星发生遮挡。
可选地,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值和第二目标坐标值确定第三目标角度,其中,第三目标角度用于表征中继终端天线的框架角;基于第三目标角度确跟踪条件。
可选地,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在第三目标角度的角度值大于第一目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡;在第三目标角度的角度值小于等于第一目标阈值的情况下,确定中继终端天线与中继卫星未发生遮挡。
可选地,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值和第二目标坐标值确定中继终端天线与中继卫星二者之间的第一矢量;基于第一目标坐标值确定中继终端天线与太阳二者之间的第二矢量;基于第一矢量与第二矢量确定中继终端天线与中继卫星二者之间的第四目标角度;基于第四目标角度确定跟踪条件。
可选地,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡;在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线与中继卫星未发生遮挡。
可选地,在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡,包括:在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,确定中继终端天线发生日凌现象;在中继终端天线发生日凌现象的情况下,确定中继终端天线与中继卫星之间发生遮挡。
可选地,在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线与中继卫星间未相互遮挡,包括:在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线未发生日凌现象;在中继终端天线未发生日凌现象的情况下,确定中继终端天线与中继卫星之间未发生遮挡。
可选地,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值和第二目标坐标值确定中继终端天线与中继卫星之间的目标射线;确定目标射线与航天器之间的相对位置关系;基于相对位置关系确定跟踪条件。
可选地,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在相对位置关系为目标射线与航天器之间相交的情况下,确定中继终端天线与中继卫星发生遮挡;在相对位置关系为目标射线与航天器之间未相交的情况下,确定中继终端天线与中继卫星未发生遮挡。
根据本发明实施例的另一方面,还提供了一种航天器中继终端天线跟踪预报装置,包括:采集单元,用于采集航天器的形状信息;转换单元,用于基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上;处理单元,用于获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件;确定单元,用于基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。
根据本发明实施例的另一方面,还提供了一种计算机可读存储介质,该计算机可读存储介质包括存储的程序,其中,在该程序运行时控制计算机可读存储介质所在设备执行本发明实施例的航天器中继终端天线跟踪预报方法。
根据本发明实施例的另一方面,还提供了一种处理器,该处理器用于运行程序,其中,该程序运行时执行本发明实施例的航天器中继终端天线跟踪预报方法。
在本发明实施例中,采用采集航天器的形状信息;基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上;获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件;基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。也就是说,本申请通过航天器的星历数据和姿态数据将航天器的形状信息进行转换,从而得到安装于航天器上的中继终端天线的第一目标坐标值,然后获取中继卫星的第二目标坐标值,并将中继终端天线与中继卫星的坐标值输入至预设判断条件组中,从而通过得到的跟踪条件确定中继终端天线与中继卫星之间是否发生了遮挡,进而解决了由于无法确定中继终端天线与中继卫星之间是否遮挡,导致中继终端天线跟踪预报准确率低的技术问题,达到了通过确定中继终端天线与中继卫星之间的遮挡关系,提高中继终端天线跟踪预报准确率的技术效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的一种航天器中继终端天线跟踪预报方法的流程图;
图2是根据本发明实施例的一种航天器的简化几何模型示意图;
图3是根据本发明实施例的一种组合体的简化几何模型示意图;
图4是根据本发明实施例的一种航天器本体坐标系与中继终端天线坐标系的关系示意图;
图5是根据本发明实施例的一种中继终端天线与中继卫星之间存在地球遮挡的关系示意图:
图6是根据本发明实施例的一种中继终端天线坐标系示意图;
图7是根据本发明实施例的一种中继终端天线与中继卫星的太阳相对几何关系示意图:
图8是根据本发明实施例的另一种中继终端天线与中继卫星的太阳相对几何关系示意图;
图9是根据本发明实施例的一种射线与矩形平板相交计算的关系示意图;
图10是根据本发明实施例的一种射线与圆柱体相交计算的关系示意图;
图11是根据本发明实施例的一种航天器组合体中继跟踪预报计算方法的流程图;
图12是根据本发明实施例的一种航天器中继终端天线跟踪预报装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
实施例1
根据本发明实施例,提供了一种航天器中继终端天线跟踪预报方法的实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种航天器中继终端天线跟踪预报方法的流程图。如图1所示,该方法可以包括如下步骤:
步骤S102,采集航天器的形状信息。
在本发明上述步骤S102提供的技术方案中,可以设定航天器模型,然后通过根据工业部门提供的航天器三维模型文件以及相应尺寸参数文件,使用三维建模软件对上述航天器模型进行模型编辑和处理,以采集航天器的形状信息,也即,将大型航天器分解为多个圆柱体、平板等子模型,并获取各子模型对应的等效形状的选定位置坐标,并将组合体简化为多个圆柱体与多个平板的组合体,对于不同的等效形状选定位置不同,以代表等效形状的轮廓信息或者关键点参数信息等。
举例而言,图2是根据本发明实施例的一种航天器的简化几何模型示意图,图3是根据本发明实施例的一种组合体的简化几何模型示意图,如图2和图3所示,在该实施例中,可以将航天器的舱体等效成多个圆柱体,其中圆柱体为航天器某舱段的最大包络圆柱、圆柱体的直径为航天器该舱段(包括外部载荷部件)的最大直径尺寸,可以将航天器上的太阳电池板简化为矩形平板,选定特征参数信息为四个顶点位置坐标。
步骤S104,基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上。
在本发明上述步骤S104提供的技术方案中,在采集到航天器的形状信息之后,由于中继终端天线是安装在航天器上的,因此可以基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,从而可以得到中继终端天线的第一目标坐标值。
在该实施例中,可以以航天器及组合体的质心和真实本体坐标系为基础,建立简化航天器模型的本体坐标系,简化航天器模型的本体坐标系的质心为真实航天器的质心,坐标系与真实本体坐标系一致,然后设计和明确航天器姿态数据所使用的本体坐标系的定义,根据姿态定义在外部提供的航天器模型上绘制坐标轴。
在上述实施例中,构成航天器三维模型的简单几何体的关键参数信息可以为:每个圆柱体可以采用圆柱底面圆心坐标、顶面圆心坐标、圆柱体半径和圆柱高度描述,矩形平板即太阳帆板采用四个顶点的坐标参数描述;中继终端天线坐标原点在本体坐标系中的坐标与安装角度。
在该实施例中,可以通过航天器本体坐标系与中继终端天线坐标系之间的关系来确定中继终端天线在航天器上的安装位置,图4是根据本发明实施例的一种航天器本体坐标系与中继终端天线坐标系的关系示意图,如图4所示,中继终端天线可以安装在航天器本体坐标系的第III象限偏第IV象限α角的位置,其中,航天器本体坐标系到中继终端天线坐标系的旋转矩阵可以通过如下公式表示:
Lrb=RY(α-180)RZ(90)
其中,Lrb可以用于表示旋转矩阵,RY可以用于表示绕Y轴旋转某一角度的旋转矩阵,RZ可以用于表示绕Z轴旋转某一角度的旋转矩阵。
可选地,在实际应用中,可以通过与工业设计部门确认框架角的定义中继终端天线坐标,以选定和设置中继终端天线参数的设定方法。中继终端天线坐标定义的原则是进行坐标选择,将该实施例中的方位轴旋转到中继终端天线的方位轴,中继终端天线坐标绕方位轴的旋转角度可以不再旋转,通过调整可见区间使其与中继终端天线实际波束范围一致即可。
步骤S106,获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件。
在本发明上述步骤S106提供的技术方案中,可以建立中继终端天线坐标系,然后通过中继终端天线坐标系获取中继卫星的第二目标坐标值,并将中继终端天线的第一目标坐标值与中继卫星的第二目标坐标值输入至预设判断条件组,从而得到跟踪条件,该跟踪条件也即为可跟踪的判断条件。
在上述实施例中,预设判断条件组可以包括至少四个条件,其中,第一条件为用于判断中继终端天线和中继卫星之间是否有地球遮挡,第二条件为用于通过中继终端天线的框架角判断中继终端天线与中继卫星之间是否发生遮挡,第三条件为用于通过确定中继终端天线是否发生日凌现象判断中继终端天线与中继卫星之间是否发生遮挡,第四条件为用于判断中继终端天线与中继卫星之间是否有多个帆板、舱体遮挡,若通过其中任意一个条件即可判断中继终端天线与中继卫星之间发生遮挡,即可确定中继终端天线与中继卫星之间不可见,而不必继续其他条件的判断。
步骤S108,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。
在本发明上述步骤S108提供的技术方案中,跟踪条件可以为上述预设条件组中的任意一个条件,在得到跟踪条件之后,可以通过跟踪条件来确定中继终端天线与中继卫星之间是否发生遮挡。
通过本申请上述步骤S102至步骤S108,采集航天器的形状信息;基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上;获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件;基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。也就是说,该实施例通过航天器的星历数据和姿态数据将航天器的形状信息进行转换,从而得到安装于航天器上的中继终端天线的第一目标坐标值,然后获取中继卫星的第二目标坐标值,并将中继终端天线与中继卫星的坐标值输入至预设判断条件组中,从而通过得到的跟踪条件确定中继终端天线与中继卫星之间是否发生了遮挡,进而解决了由于无法确定中继终端天线与中继卫星之间是否遮挡,导致中继终端天线跟踪预报准确率低的技术问题,达到了通过确定中继终端天线与中继卫星之间的遮挡关系,提高中继终端天线对中继卫星的跟踪预报准确率的技术效果。
下面对该实施例的上述方法进行进一步介绍。
作为一种可选的实施方式,获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值确定中继终端天线与地心二者之间的第一目标距离;基于第二目标坐标值确定中继卫星与地心二者之间的第二目标距离;基于第一目标距离与第二目标距离确定第一目标角度的角度值和第二目标角度的角度值,其中,第一目标角度为第一目标距离和第二目标距离分别与同一地球半径的夹角和,第二目标角度为第一距离与第二距离之间的夹角;基于第一目标角度的角度值和第二目标角度的角度值确定跟踪条件。
在该实施例中,可以通过中继终端天线的第一目标坐标值确定中继终端天线到地球地心之间的第一目标距离,然后通过中继卫星的第二目标坐标值确定中继卫星到地球地心之间的第二目标距离,从而可以通过第一目标距离、第二目标距离以及地球半径计算出第一目标角度的角度值和第二目标角度的角度值。
在上述实施例中,图5是根据本发明实施例的一种中继终端天线与中继卫星之间存在地球遮挡的关系示意图,如图5所示,设地球半径为RE,中继终端天线到地心的距离为r1,中继卫星到地心的距离为r2,计算第一目标角度的角度值和第二目标角度的角度值的计算公式可以如下所示:
其中,φ可以用于表示第一目标角度,ε可以用于表示第二目标角度。
在上述实施例中,第一目标角度φ可以为第一目标距离和第二目标距离分别与同一地球半径的夹角和,第二目标角度为第一距离与第二距离之间的夹角。在确定出第一目标角度和第二目标角度之后,可以通过第一目标角度的角度值和第二目标角度的角度值确定跟踪条件,该跟踪条件可以为确定中继终端天线与中继卫星之间是否有地球进行遮挡。
作为一种可选的实施方式,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在第一目标角度的角度值大于第二目标角度的角度值的情况下,确定中继终端天线与中继卫星未发生遮挡;在第一目标角度的角度值小于等于第二目标角度的角度值的情况下,确定中继终端天线与中继卫星发生遮挡。
在该实施例中,在确定出第一目标角度和第二目标角度之后,可以判断第一目标角度和第二目标角度的角度值的大小,在第一目标角度的角度值大于第二目标角度的角度值的情况下,确定中继终端天线与中继卫星未发生遮挡,也即,当φ>ε时,中继终端天线与中继卫星之间无地球遮挡,此时中继终端天线与中继卫星未发生遮挡,相互之间可见;在第一目标角度的角度值小于等于第二目标角度的角度值的情况下,确定中继终端天线与中继卫星发生遮挡,也即,当φ≤ε时,中继终端天线与中继卫星之间存在地球遮挡情况,此时,中继终端天线与中继卫星发生遮挡,相互之间不可见。
作为一种可选的实施方式,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值和第二目标坐标值确定第三目标角度,其中,第三目标角度用于表征中继终端天线的框架角;基于第三目标角度确跟踪条件。
在该实施例中,图6是根据本发明实施例的一种中继终端天线坐标系示意图,如图6所示,可以通过中继终端天线的第一目标坐标值和中继卫星的第二目标坐标值计算出第三目标角度,也即中继终端天线的框架角,并通过该框架角确定跟踪条件,该跟踪条件可以为通过框架角的大小确定中继终端天线与中继卫星之间发生遮挡。
在上述实施例中,可以在中继终端天线坐标系中计算出框架角,则计算中继终端天线的框架角的计算公式可以如下所示:
其中,中继卫星的坐标可以为Pr=(xr,yr,zr),以中继终端天线的第一目标坐标为中继终端天线坐标系的原点,则框架角可以用a和b来表示,且-180°<a<180°,-90°<b<90°,r可以用于表示中继终端天线与中继卫星之间的距离。
在上述实施例中,中继终端天线安装在航天器上,且中继终端天线的基座与航天器本体之间可以固连,其安装方式可以有多种,通常用绕2个轴的旋转角度表示,即框架AB角,也称框架角,可以通过坐标旋转采用统一的模型表示。如图4所示,本申请中可以按先绕x轴旋转a角,然后绕y轴旋转b角的方式定义中继终端天线的框架角,中继终端天线的轴心是中继终端天线坐标系的Z轴,即零位为正Z,X轴是方位轴,按右手定则可以确定正方向,b角是中继终端天线-中继卫星连线与YZ平面的夹角,朝向+X为正。
作为一种可选的实施方式,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在第三目标角度的角度值大于第一目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡;在第三目标角度的角度值小于等于第一目标阈值的情况下,确定中继终端天线与中继卫星未发生遮挡。
在该实施例中,在确定好第三目标角度,也即中继终端天线的框架角之后,在第三目标角度的角度值大于第一目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡,也即在|a|>90°或|b|>90°的情况下,中继终端天线与中继星之间不可见,在第三目标角度的角度值小于等于第一目标阈值的情况下,确定中继终端天线与中继卫星未发生遮挡,也即在其它情况下,中继终端天线与中继星之间可见。
可选地,在该实施例中,在实际情况下,航天器上安装的中继终端天线的指向会受物理条件约束。
作为一种可选的实施方式,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值和第二目标坐标值确定中继终端天线与中继卫星二者之间的第一矢量;基于第一目标坐标值确定中继终端天线与太阳二者之间的第二矢量;基于第一矢量与第二矢量确定中继终端天线与中继卫星二者之间的第四目标角度;基于第四目标角度确定跟踪条件。
在该实施例中,在获取到中继终端天线的第一目标坐标值和中继卫星的第二目标坐标值之后,可以确定中继终端天线与中继卫星二者之间的第一矢量,以及中继终端天线与太阳二者之间的第二矢量,通过第一矢量和第二矢量计算中继终端天线与中继卫星二者之间的第四目标角度,也即,第一矢量与第二矢量之间的夹角,并通过该夹角确定跟踪条件,该跟踪条件可以为中继终端天线是否发生日凌现象。
在上述实施例中,图7是根据本发明实施例的一种中继终端天线与中继卫星的太阳相对几何关系示意图,图8是根据本发明实施例的另一种中继终端天线与中继卫星的太阳相对几何关系示意图,如图7和图8所示,可以通过航天器姿态轨道运动仿真来获得在特定姿态轨道条件下的相关矢量,例如,在地心J2000惯性坐标系下,设中继卫星的坐标矢量为中继终端天线的坐标矢量为太阳的坐标矢量为则中继终端天线到中继卫星的第一矢量、中继终端天线到太阳的第二矢量的计算公式可以如下所示:
则第一矢量与第二矢量之间的夹角β的计算公式可以如下所示:
作为一种可选的实施方式,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡;在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线与中继卫星未发生遮挡。
在该实施例中,在确定好第一矢量与第二矢量之间的夹角之后,可以在该夹角的角度值小于第二目标阈值,或该夹角的补角的角度值小于第二目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡,也即,β<δ或者180°-β<δ时,确定中继终端天线与中继卫星之间不可见;在该夹角的角度值大于等于第二目标阈值,或该夹角的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线与中继卫星未发生遮挡,也即,中继终端天线与中继卫星之间可见。可选地,上述δ是日凌现象发生的门限值,与天线具体参数有关,由天线制造部门提供。
作为一种可选的实施方式,在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,确定中继终端天线与中继卫星发生遮挡,包括:在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,确定中继终端天线发生日凌现象;在中继终端天线发生日凌现象的情况下,确定中继终端天线与中继卫星之间发生遮挡。
在该实施例中,由于太阳电磁辐射会对中继终端天线与中继卫星之间的通信和测量造成干扰,当太阳矢量与中继终端天线到中继卫星的连线的夹角小于一定值时干扰噪声将超过信号强度,会对通信和测量造成严重干扰甚至通信中断,这种现象就是日凌现象。对于近地大型航天器轨道,日凌现象发生时间不规则,可以采用数值计方法进行计算和预报。
在上述实施例中,由于航天测控通常是双向的,对于中继卫星的天线同样存在日凌问题,因此,日凌现象发生的判据可以为,在第四目标角度的角度值小于第二目标阈值,或第四目标角度的补角的角度值小于第二目标阈值的情况下,也即β<δ或者180°-β<δ的情况下,可以确定中继终端天线与中继卫星的天线发生了日凌现象,此时,中继终端天线与中继卫星之间由于日凌现象无法进行通信,可以确定中继终端天线与中继卫星之间发生遮挡。
作为一种可选的实施方式,在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线与中继卫星间未相互遮挡,包括:在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,确定中继终端天线未发生日凌现象;在中继终端天线未发生日凌现象的情况下,确定中继终端天线与中继卫星之间未发生遮挡。
在该实施例中,在第四目标角度的角度值大于等于第二目标阈值,或第四目标角度的补角的角度值大于等于第二目标阈值的情况下,可以确定中继终端天线未发生日凌现象。此时,由于中继终端天线与中继卫星之间的通信未被日凌现象干扰,因此,中继终端天线可以跟踪中继卫星,也即中继终端天线与中继卫星之间未发生遮挡,相互可见。
作为一种可选的实施方式,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:基于第一目标坐标值和第二目标坐标值确定中继终端天线与中继卫星之间的目标射线;确定目标射线与航天器之间的相对位置关系;基于相对位置关系确定跟踪条件。
在该实施例中,可以将航天器的帆板和舱体简化为矩形平板和圆柱体,在获取到中继终端天线的第一目标坐标值和中继卫星的第二目标坐标值之后,可以获取中继终端天线与中继卫星之间的目标射线,通过判断目标射线与矩形平板、圆柱体的相对位置关系,确定跟踪条件,该跟踪条件可以为目标射线与矩形平板、圆柱体之间是否相交。
在该实施例中,可以将中继天线坐标系内中继终端天线-中继星之间的连线与多个帆板、舱体之间的遮挡计算问题,转换为计算机图形学中射线与多个忽略厚度的矩形平板的相交检测、与多个圆柱体的相交检测问题。
在上述实施例中,图9是根据本发明实施例的一种射线与矩形平板相交计算的关系示意图,如图9所示,设在J2000地心惯性坐标系下中继终端天线位置(中继终端天线坐标系原点)为PA、中继卫星位置为Pz,首先根据坐标系转换矩阵,计算出矩形平板的四个顶点P1、P2、P3、P4在J2000地心惯性坐标系中的坐标值,然后分别计算目标射线PAPZ与三角形P1P2P3的相对位置关系、目标射线PAPZ与三角形P1P4P3的相对位置关系。
在上述实施例中,图10是根据本发明实施例的一种射线与圆柱体相交计算的关系示意图,如图10所示,首先根据圆柱体特征点参数,建立以圆柱体方向为轴的坐标系,其中圆柱体底面圆心为坐标系原点,计算航天器本体坐标系到圆柱体坐标系的旋转矩阵,计算中继终端位置PA和中继卫星位置Pz在该圆柱体坐标系下的坐标值,之后进行相交判断计算:首先判断线段方向是否与y轴方向平行,如果是,判断是否是圆柱体的上下圆面相交;其次判断线段是否与圆柱面相交,先假设圆柱是无限高,判断的方法是求解二次方程;假设线段和圆柱相交,则必定存在线段上存在某一点t,满足该点到y轴的距离为圆柱的半径;可以先通过判断根的个数,通过求解方程可以得到射线与圆柱体之间的相对位置关系。
作为一种可选的实施方式,基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡,包括:在相对位置关系为目标射线与航天器之间相交的情况下,确定中继终端天线与中继卫星发生遮挡;在相对位置关系为目标射线与航天器之间未相交的情况下,确定中继终端天线与中继卫星未发生遮挡。
在该实施例中,在相对位置关系为目标射线与矩形平板或圆柱体之间相交,则确定中继终端天线与中继卫星发生遮挡,也即中继终端天线与中继卫星之间不可见,在相对位置关系为目标射线与矩形平板或圆柱体之间不相交,则确定中继终端天线与中继卫星未发生遮挡,也即中继终端天线与中继卫星之间可见。
作为一种可选的实施方式,在得到中继终端天线跟踪预报计算结果之后,可以对该计算结果进行验证。
在该实施例中,可以基于系统仿真工具包(商业软件STK)进行中继跟踪预报结果的比对与验证工作,该商业软件STK是美国Analytical Graphics公司开发的一款在航天领域处于领先地位的商业分析软件。通过该STK进行验证的方法可以如下:
首先建立两个航天器或组合体的太阳电池阵三维模型和航天器三维本体模型,然后使用格式转换工具,转换为可以导入到STK中*.mdl格式模型文件;然后基于STK建立航天器轨道分析环境,分别导入两个航天器的太阳电池阵模型和航天器本体模型;之后设置STK中传感器对象的安装位置和指向角度限制;最后分别设置两个航天器的轨道参数和姿态参数,在某时间段选择若干特征时刻点,比对STK中中继跟踪可见性结果与使用该实施例中的方法得到的预报结果是否一致,在结果一致的情况下,确定使用该实施例中的方法得到的预报结果正确,在结果不一致的情况下,对使用该实施例中的方法得到的预报结果进行再次计算,直到结果一致为止。
该实施例的中继终端天线跟踪预报方法,通过航天器的星历数据和姿态数据将航天器的形状信息进行转换,从而得到安装于航天器上的中继终端天线的第一目标坐标值,然后获取中继卫星的第二目标坐标值,并将中继终端天线与中继卫星的坐标值输入至预设判断条件组中,从而通过得到的跟踪条件确定中继终端天线与中继卫星之间是否发生了遮挡,进而解决了由于无法确定中继终端天线与中继卫星之间是否遮挡,导致中继终端天线跟踪预报准确率低的技术问题,达到了通过确定中继终端天线与中继卫星之间的遮挡关系,提高中继终端天线的跟踪预报准确率的技术效果。
实施例2
下面结合优选的实施方式对本发明实施例的技术方案进行举例说明。
图11是根据本发明实施例的一种航天器组合体中继跟踪预报计算方法的流程图,如图11所示,该方法可以包括如下步骤:
步骤S1101,航天器及其组合体三维模型分解与简化建模。
步骤S1102,根据航天器星历数据和姿态数据计算各转换矩阵和坐标值。
步骤S1103,地球相对几何关系计算与地球遮挡判断。
步骤S1104,中继终端天线框架角计算与约束条件判断。
步骤S1105,太阳相对几何关系与日凌现象判断。
步骤S1106,多舱体与帆板相对几何关系计算与遮挡判断。
步骤S1107,输出中继终端天线的框架角和跟踪弧段预报结果。
在该实施例中,针对空间站任务期间航天器中继终端天线将会面临较为复杂的自身及其他航天器舱体与太阳帆板遮挡的问题,提出了一种航天器组合体中继跟踪预报计算方法,包括中继终端与中继卫星之间的地球遮挡判断、框架角约束条件计算与判断、太阳运动导致日凌现象判断以及航天器舱体及帆板的遮挡计算判断方法,在保证满足最大空间包络且空间位置关系准确符合航天器实际构型的前提下,可以将复杂大型航天器简化为多个圆柱体与矩形平板的组合体,通过处理和计算多个不同几何体的特征参数,实现中继终端天线与中继卫星之间的遮挡判断计算与中继跟踪预报功能。
在上述实施例中,鉴于中继跟踪预报的前提条件可以为中继天线中心与中继星几何可见(无地球遮挡)、无航天器舱体或帆板遮挡、中继框架角满足约束条件、无日凌现象,因此,该实施例可以先进行中继天线与中继星之间的地球遮挡判断,之后进行中继跟踪框架角计算和判断,然后进行日凌判断计算、最后进行航天器舱体或太阳帆板对中继星以及对太阳的遮挡判断计算,其中如任意一项计算中满足遮挡判据条件,即中断计算流程、返回遮挡状态判别结果,认为中继终端(中继终端天线或中继天线)与中继卫星之间不可见。
在该实施例中,可以设计航天器中继跟踪计算模型,在该模型中,主要涉及J2000地心惯性坐标系、航天器轨道坐标系、航天器本体坐标系、中继终端天线坐标系以及航天器姿态角等,通过各种坐标系之间的转换计算可以得到中继终端天线的框架角参数。
在上述实施例中,J2000地心惯性坐标系的原点为地心,x轴指向公元2000年的平春分点,z轴垂直于平赤道,y轴由右手定则确定;航天器轨道坐标系是航天器的姿态参考系,坐标系原点在飞行器质心,z轴指向地心,x轴与z垂直指向飞行方向,y轴垂直轨道面;航天器本体坐标系的原点位于航天器的质心或几何中心,x轴由坐标原点指向航天器的某一特征轴方向,z轴由坐标原点指向航天器的另一特征轴方向,y轴与x轴、z轴构成右手系。
在上述实施例中,设J2000地心惯性坐标系到航天器本体坐标系的变换矩阵为Lbi,轨道坐标系到本体坐标系的变换矩阵为Lbo,惯性坐标系与轨道坐标系的变换矩阵为Loi,则三者之间的关系为:
Lbi=Lbo·Loi
在上述实施例中,航天器姿态可以采用欧拉角描述,表示的是飞行器本体系相对当地轨道系的转动关系,飞行器姿态角的转系为3-1-2转系,即按偏航、滚动、俯仰角的顺序定义。航天器轨道坐标系到航天器本体坐标系的坐标变换矩阵为:
其中,Ry(θ)可以用于表示绕y轴转动角度θ的旋转矩阵,Rx(φ)可以用于表示绕x轴转动角度φ的旋转矩阵,Rz(ψ)可以用于表示绕z轴转动角度ψ的旋转矩阵。
在上述实施例中,在获取航天器的简化几何体模型参数之后,可以根据航天器在某时刻的位置与姿态,对中继终端天线进行旋转变换,比如,根据当前时刻的航天器姿态(包括舱体转动角度、太阳翼转动角度等)进行旋转变换。
在相关技术中,在空间实验室任务期间,两大型航天器组合体为三轴对地或连续偏航两种姿态,中继天线遮挡计算模型比较简单。随着我国空间站的建设与发展、飞控测控任务设计要求越来越精细化。大型航天器组合体姿态种类更多、构型更加复杂,,中继天线存在较为复杂的舱体以及太阳帆板遮挡的问题。
在上述相关技术中,由于航天器表面通常难以用解析表达式来描述,因此要进行中继跟踪遮挡判断计算,首先要将航天器三维模型表面细分为许多个微小平面面元并获取每个小面元的顶点坐标和法向矢量,但是面元之间相互遮挡关系的判断过程比较复杂,需要多次循环计算和判断比较,这一过程也可被称为目标的消隐处理。单纯地通过编程执行循环计算和判断来实现目标消隐计算,比较费时且精度较低;尤其对于复杂结构大型航天器,小面元划分的越细则面元数量越多,导致遮挡判断计算时间呈指数级增加,从而导致中继终端天线跟踪预报的准确率较低的问题。
而该实施例中的航天器组合体中继跟踪预报计算方法,相比于相关技术中的仅利用几何可见或框架角约束条件进行预报计算,无法针对多航天器动态遮挡计算方法,通过将航天器简化建模、将相关矢量与模型投影到中继天线坐标系进行遮挡判断,实现了动态、准确、定量、快速地遮挡判断与计算,可为航天器中继终端在轨利用提供可靠依据。
实施例3
本发明实施例还提供了一种航天器中继终端天线跟踪预报装置。需要说明的是,该实施例的航天器中继终端天线跟踪预报装置可以用于执行本发明实施例的航天器中继终端天线跟踪预报方法。
图12是根据本发明实施例的一种航天器中继终端天线跟踪预报装置的示意图。如图12所示,该航天器中继终端天线跟踪预报装置1200可以包括:采集单元1201、转换单元1202、处理单元1203和确定单元1204。
采集单元1201,用于采集航天器的形状信息。
转换单元1202,用于基于航天器的星历数据和姿态数据对航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,中继终端天线位于航天器上。
处理单元1203,用于获取中继卫星的第二目标坐标值,将第一目标坐标值和第二目标坐标值输入至预设判断条件组,得到跟踪条件。
确定单元1204,用于基于跟踪条件确定中继终端天线与中继卫星是否发生遮挡。
该实施例中的航天器中继终端天线跟踪预报装置,通过航天器的星历数据和姿态数据将航天器的形状信息进行转换,从而得到安装于航天器上的中继终端天线的第一目标坐标值,然后获取中继卫星的第二目标坐标值,并将中继终端天线与中继卫星的坐标值输入至预设判断条件组中,从而通过得到的跟踪条件确定中继终端天线与中继卫星之间是否发生了遮挡,进而解决了由于无法确定中继终端天线与中继卫星之间是否遮挡,导致中继终端天线跟踪预报准确率低的技术问题,达到了通过确定中继终端天线与中继卫星之间的遮挡关系,提高中继终端天线对中继卫星的跟踪预报准确率的技术效果。
实施例4
根据本发明实施例,还提供了一种计算机可读存储介质,该计算机可读存储介质包括存储的程序,其中,该程序执行实施例1中所述的航天器中继终端天线跟踪预报方法。
实施例5
根据本发明实施例,还提供了一种处理器,该处理器用于运行程序,其中,所述程序运行时执行实施例1中所述的航天器中继终端天线跟踪预报方法。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如所述单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模块的间接耦合或通信连接,可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (14)
1.一种航天器中继终端天线跟踪预报方法,其特征在于,包括:
采集航天器的形状信息;
基于所述航天器的星历数据和姿态数据对所述航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,所述中继终端天线位于所述航天器上;
获取中继卫星的第二目标坐标值,将所述第一目标坐标值和所述第二目标坐标值输入至预设判断条件组,得到跟踪条件;
基于所述跟踪条件确定所述中继终端天线与所述中继卫星是否发生遮挡。
2.根据权利要求1所述的方法,其特征在于,获取中继卫星的第二目标坐标值,将所述第一目标坐标值和所述第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:
基于所述第一目标坐标值确定所述中继终端天线与地心二者之间的第一目标距离;
基于所述第二目标坐标值确定所述中继卫星与所述地心二者之间的第二目标距离;
基于所述第一目标距离与所述第二目标距离确定第一目标角度的角度值和第二目标角度的角度值,其中,所述第一目标角度为所述第一目标距离和所述第二目标距离分别与同一地球半径的夹角和,所述第二目标角度为所述第一距离与所述第二距离之间的夹角;
基于所述第一目标角度的角度值和所述第二目标角度的角度值确定所述跟踪条件。
3.根据权利要求2所述的方法,其特征在于,基于所述跟踪条件确定所述中继终端天线与所述中继卫星是否发生遮挡,包括:
在所述第一目标角度的角度值大于所述第二目标角度的角度值的情况下,确定所述中继终端天线与所述中继卫星未发生遮挡;
在所述第一目标角度的角度值小于等于所述第二目标角度的角度值的情况下,确定所述中继终端天线与所述中继卫星发生遮挡。
4.根据权利要求1所述的方法,其特征在于,将所述第一目标坐标值和所述第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:
基于所述第一目标坐标值和所述第二目标坐标值确定第三目标角度,其中,所述第三目标角度用于表征所述中继终端天线的框架角;
基于所述第三目标角度确所述跟踪条件。
5.根据权利要求4所述的方法,其特征在于,基于所述跟踪条件确定所述中继终端天线与所述中继卫星是否发生遮挡,包括:
在所述第三目标角度的角度值大于第一目标阈值的情况下,确定所述中继终端天线与所述中继卫星发生遮挡;
在所述第三目标角度的角度值小于等于第一目标阈值的情况下,确定所述中继终端天线与所述中继卫星未发生遮挡。
6.根据权利要求1所述的方法,其特征在于,将所述第一目标坐标值和所述第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:
基于所述第一目标坐标值和所述第二目标坐标值确定所述中继终端天线与所述中继卫星二者之间的第一矢量;
基于所述第一目标坐标值确定所述中继终端天线与太阳二者之间的第二矢量;
基于所述第一矢量与所述第二矢量确定所述中继终端天线与所述中继卫星二者之间的第四目标角度;
基于所述第四目标角度确定所述跟踪条件。
7.根据权利要求6所述的方法,基于所述跟踪条件确定所述中继终端天线与所述中继卫星是否发生遮挡,包括:
在所述第四目标角度的角度值小于第二目标阈值,或所述第四目标角度的补角的角度值小于所述第二目标阈值的情况下,确定所述中继终端天线与所述中继卫星发生遮挡;
在所述第四目标角度的角度值大于等于第二目标阈值,或所述第四目标角度的补角的角度值大于等于所述第二目标阈值的情况下,确定所述中继终端天线与所述中继卫星未发生遮挡。
8.根据权利要求7所述的方法,其特征在于,在所述第四目标角度的角度值小于第二目标阈值,或所述第四目标角度的补角的角度值小于所述第二目标阈值的情况下,确定所述中继终端天线与所述中继卫星发生遮挡,包括:
在所述第四目标角度的角度值小于第二目标阈值,或所述第四目标角度的补角的角度值小于所述第二目标阈值的情况下,确定所述中继终端天线发生日凌现象;
在所述中继终端天线发生日凌现象的情况下,确定所述中继终端天线与所述中继卫星之间发生遮挡。
9.根据权利要求7所述的方法,其特征在于,在所述第四目标角度的角度值大于等于第二目标阈值,或所述第四目标角度的补角的角度值大于等于所述第二目标阈值的情况下,确定所述中继终端天线与所述中继卫星间未相互遮挡,包括:
在所述第四目标角度的角度值大于等于第二目标阈值,或所述第四目标角度的补角的角度值大于等于所述第二目标阈值的情况下,确定所述中继终端天线未发生日凌现象;
在所述中继终端天线未发生日凌现象的情况下,确定所述中继终端天线与所述中继卫星之间未发生遮挡。
10.根据权利要求1所述的方法,其特征在于,将所述第一目标坐标值和所述第二目标坐标值输入至预设判断条件组,得到跟踪条件,包括:
基于所述第一目标坐标值和所述第二目标坐标值确定所述中继终端天线与所述中继卫星之间的目标射线;
确定所述目标射线与所述航天器之间的相对位置关系;
基于所述相对位置关系确定所述跟踪条件。
11.根据权利要求10所述的方法,其特征在于,基于所述跟踪条件确定所述中继终端天线与所述中继卫星是否发生遮挡,包括:
在所述相对位置关系为所述目标射线与所述航天器之间相交的情况下,确定所述中继终端天线与所述中继卫星发生遮挡;
在所述相对位置关系为所述目标射线与所述航天器之间未相交的情况下,确定所述中继终端天线与所述中继卫星未发生遮挡。
12.一种航天器中继终端天线跟踪预报装置,其特征在于,包括:
采集单元,用于采集航天器的形状信息;
转换单元,用于基于所述航天器的星历数据和姿态数据对所述航天器的形状信息进行转换,得到中继终端天线的第一目标坐标值,其中,所述中继终端天线位于所述航天器上;
处理单元,用于获取中继卫星的第二目标坐标值,将所述第一目标坐标值和所述第二目标坐标值输入至预设判断条件组,得到跟踪条件;
确定单元,用于基于所述跟踪条件确定所述中继终端天线与所述中继卫星是否发生遮挡。
13.一种计算机可读存储介质,其特征在于,所述存储介质包括存储的程序,其中,在所述程序运行时控制所述存储介质所在设备执行权利要求1至11中任意一项所述的航天器中继终端天线跟踪预报方法。
14.一种处理器,其特征在于,所述处理器用于运行程序,其中,所述程序运行时执行权利要求1至11中任意一项所述的航天器中继终端天线跟踪预报方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110417513.2A CN113158341A (zh) | 2021-04-19 | 2021-04-19 | 航天器中继终端天线跟踪预报方法、装置和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110417513.2A CN113158341A (zh) | 2021-04-19 | 2021-04-19 | 航天器中继终端天线跟踪预报方法、装置和存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113158341A true CN113158341A (zh) | 2021-07-23 |
Family
ID=76868488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110417513.2A Pending CN113158341A (zh) | 2021-04-19 | 2021-04-19 | 航天器中继终端天线跟踪预报方法、装置和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113158341A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113595618A (zh) * | 2021-07-29 | 2021-11-02 | 东方红卫星移动通信有限公司 | 一种用于卫星通信的日凌干扰角度和日凌时间的预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160170029A1 (en) * | 2014-12-10 | 2016-06-16 | Worldvu Satellites Limited | User Terminal Having a Linear Array Antenna with Electronic and Mechanical Actuation System |
CN107944094A (zh) * | 2017-11-06 | 2018-04-20 | 中国航天空气动力技术研究院 | 一种复杂外形航天器投影面积的确定方法及系统 |
CN110161493A (zh) * | 2019-04-23 | 2019-08-23 | 中国西安卫星测控中心 | 多约束条件下航天器跟踪预报方法 |
RU2706636C1 (ru) * | 2018-10-17 | 2019-11-19 | Публичное акционерное общество "Ракетно-космическая корпорация "Энергия" имени С.П. Королёва" | Способ определения координат космического аппарата по сигналам навигационных спутников и устройство определения координат космического аппарата по сигналам навигационных спутников |
CN111008481A (zh) * | 2019-12-18 | 2020-04-14 | 北京空天技术研究所 | 航天器的气动分析方法及装置 |
CN112540390A (zh) * | 2020-11-26 | 2021-03-23 | 陕西星邑空间技术有限公司 | 一种航天器的跟踪预报计算方法及装置 |
-
2021
- 2021-04-19 CN CN202110417513.2A patent/CN113158341A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160170029A1 (en) * | 2014-12-10 | 2016-06-16 | Worldvu Satellites Limited | User Terminal Having a Linear Array Antenna with Electronic and Mechanical Actuation System |
CN107944094A (zh) * | 2017-11-06 | 2018-04-20 | 中国航天空气动力技术研究院 | 一种复杂外形航天器投影面积的确定方法及系统 |
RU2706636C1 (ru) * | 2018-10-17 | 2019-11-19 | Публичное акционерное общество "Ракетно-космическая корпорация "Энергия" имени С.П. Королёва" | Способ определения координат космического аппарата по сигналам навигационных спутников и устройство определения координат космического аппарата по сигналам навигационных спутников |
CN110161493A (zh) * | 2019-04-23 | 2019-08-23 | 中国西安卫星测控中心 | 多约束条件下航天器跟踪预报方法 |
CN111008481A (zh) * | 2019-12-18 | 2020-04-14 | 北京空天技术研究所 | 航天器的气动分析方法及装置 |
CN112540390A (zh) * | 2020-11-26 | 2021-03-23 | 陕西星邑空间技术有限公司 | 一种航天器的跟踪预报计算方法及装置 |
Non-Patent Citations (7)
Title |
---|
YAO ZHANG 等: "An improved reaching law of SMC for spacecraft tracking system", 2015 IEEE INTERNATIONAL CONFERENCE ON INFORMATION AND AUTOMATION, 1 October 2015 (2015-10-01) * |
刘勇 等: "引入姿态误差的中继跟踪预报", 飞行器测控学报, vol. 31, no. 04, 14 August 2012 (2012-08-14), pages 11 - 14 * |
吴刚 等: "航空器目标的星间天线指向实时监视系统设计", 测试技术学报, vol. 24, no. 05, 15 September 2010 (2010-09-15), pages 385 - 389 * |
周倜 等: "航天器天线的建模与仿真", 系统仿真学报, vol. 28, no. 09, 8 September 2016 (2016-09-08), pages 2042 - 2048 * |
李于衡 等: "中继卫星与用户星双向跟踪关系及策略", 西安电子科技大学学报, no. 01, 20 February 2007 (2007-02-20), pages 6 - 10 * |
李于衡 等: "在轨跟踪与数据中继卫星测控关键技术(上)", 上海航天, no. 04, 25 August 2006 (2006-08-25), pages 1 - 7 * |
李于衡;易克初;田红心;: "跟踪与数据中继卫星(TDRS)跟踪用户星的条件分析", 空间科学学报, no. 05, 30 September 2006 (2006-09-30), pages 377 - 381 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113595618A (zh) * | 2021-07-29 | 2021-11-02 | 东方红卫星移动通信有限公司 | 一种用于卫星通信的日凌干扰角度和日凌时间的预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109781060B (zh) | 星载点波束天线对地指向精度评估方法 | |
CN110125936A (zh) | 一种空间机器人的柔顺控制方法和地面实验验证系统 | |
CN110231029B (zh) | 一种水下机器人多传感器融合数据处理方法 | |
CN107014398B (zh) | 卫星模拟太阳敏感器故障检测方法及装置 | |
CN109781059A (zh) | 星载点波束天线对地指向精度评估系统 | |
CN103644918A (zh) | 卫星对月探测数据定位处理方法 | |
CN113158341A (zh) | 航天器中继终端天线跟踪预报方法、装置和存储介质 | |
CN114510673A (zh) | 一种基于欧拉角转换实时计算卫星测控角的方法 | |
Qiao et al. | Asteroid approaching orbit optimization considering optical navigation observability | |
CN117113645A (zh) | 一种大尺度空间试训任务场景快速生成系统 | |
Daly et al. | Problems in radiation environment models at low altitudes | |
CN111008481A (zh) | 航天器的气动分析方法及装置 | |
CN112046794B (zh) | 基于混合高斯模型的固定时间约束的航天器集群控制方法 | |
CN111337031B (zh) | 一种基于姿态信息的航天器地标匹配自主位置确定方法 | |
CN113722958A (zh) | 一种不规则形状小天体引力场高效建模方法 | |
Darling et al. | Uncertainty propagation of correlated quaternion and Euclidean states using the Gauss-Bingham density | |
CN107065547B (zh) | 一种基于零空间方法的非合作目标自主交会策略 | |
CN108287352B (zh) | 一种卫星圆锥传感器作用范围的确认方法及系统 | |
Ticona et al. | Attitude determination and control system for nadir pointing and detumbling using magnetorquer for 1u bolivian cubesat | |
Huang et al. | Nonlinear Robust $ H_ {\infty} $ Control for Spacecraft Body-Fixed Hovering Around Noncooperative Target Via Modified $\theta-D $ Method | |
Zhang et al. | A diverse space target dataset with multidebris and realistic on-orbit environment | |
Zhou et al. | The general method of the tanker car mouth pose measurement | |
CN111880404A (zh) | 基于非奇异终端滑膜控制的航天器集群控制方法 | |
Shi et al. | Research on Starlink constellation simulation and target area visibility algorithm | |
CN115859479B (zh) | 一种用于航天器轨道设计的gvct构建仿真方法及装置 |
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 |