CN115758704A - 卫星南北保持策略模型的建模方法、模型、获取方法 - Google Patents

卫星南北保持策略模型的建模方法、模型、获取方法 Download PDF

Info

Publication number
CN115758704A
CN115758704A CN202211408067.XA CN202211408067A CN115758704A CN 115758704 A CN115758704 A CN 115758704A CN 202211408067 A CN202211408067 A CN 202211408067A CN 115758704 A CN115758704 A CN 115758704A
Authority
CN
China
Prior art keywords
satellite
inclination
inclination angle
state data
strategy
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.)
Granted
Application number
CN202211408067.XA
Other languages
English (en)
Other versions
CN115758704B (zh
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.)
Emposat Co Ltd
Original Assignee
Emposat Co Ltd
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 Emposat Co Ltd filed Critical Emposat Co Ltd
Priority to CN202211408067.XA priority Critical patent/CN115758704B/zh
Publication of CN115758704A publication Critical patent/CN115758704A/zh
Application granted granted Critical
Publication of CN115758704B publication Critical patent/CN115758704B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及航空航天领域,提供一种卫星南北保持策略模型的建模方法、模型、获取方法、设备、介质,建模方法包括:S1:初始化模型;S2:向模型中输入初始时刻卫星倾角,得到倾角控制行为及Q值;S3:获取当前时刻卫星倾角,获得卫星执行的倾角控制行为;S4:执行倾角控制行为,得到下一时刻卫星倾角;根据下一时刻卫星倾角、南北保持策略奖励函数得到奖励;S5:将卫星组合状态数据组存入经验池;S6:计算目标值;S7:计算误差,更新神经网络权值参数;S8:更新Q值;将下一时刻卫星倾角作为当前时刻卫星倾角;S9:重复执行步骤S3‑S8;S10:重复执行步骤S2‑S9,直到输入完所有卫星训练状态数据组的数据。本方案能得到最优的决策策略,减少卫星燃料消耗。

Description

卫星南北保持策略模型的建模方法、模型、获取方法
技术领域
本发明涉及航空航天技术领域,特别是涉及一种卫星南北保持策略模型的建模方法、模型、获取方法、设备、介质。
背景技术
随着人类航天活动的不断开展,越来越多的遥感卫星为人们的日常生活提供帮助。
GEO卫星在运行的过程中会受到日、月引力和地球非球形摄动的影响,从而导致在南北纬度方向上发生漂移,所以对GEO三轴稳定卫星南北位置保持(倾角保持)控制对于航天领域有至关重要的作用。
现有的技术首先通过对卫星在轨道运行过程中受到的地球形状、日月引力、太阳光压等各种摄动力的分析建立轨道机动时的动力学模型,然后制定南北保持的短期和长期策略。这种方法通过卫星在轨道运行过程中受到的各种摄动力进行复杂的建模,然而由于空间受力的复杂性以及卫星本身参数的不确定性,导致无法对卫星进行精确建模,参数多、计算复杂,进而影响卫星倾角控制的精度,可能会消耗更多的燃料。并且,现有强化学习方法无法解决状态和动作空间维度很高的问题。
因此,亟需开发一种卫星南北保持策略模型的建模方法、模型、获取方法、设备、介质,降低建模难度,精确计算南北保持策略。
发明内容
本发明的目的是提供一种卫星南北保持策略模型的建模方法、模型、获取方法、设备、介质,对GEO三轴稳定卫星进行南北位置保持时不需要进行复杂的建模,不需要考虑空间受力的复杂性以及卫星本身参数的不确定性,强化学习具有很强的行为决策能力,能够得到最优的决策策略,减少卫星燃料的消耗。
为解决上述技术问题,作为本发明的一个方面,提供了一种基于深度Q网络的卫星南北保持策略模型的建模方法,包括以下步骤:
S1:初始化模型,获取多组卫星训练状态数据组,每组卫星训练状态数据组包括卫星的初始状态、多个预期轨控时刻和预期轨控次数;每个卫星的初始状态包括初始时刻卫星倾角;
S2:向模型中输入一组卫星训练状态数据组的初始时刻卫星倾角,得到初始时刻之后所有倾角控制行为及所对应输出的Q值;
S3:获取当前时刻卫星倾角,根据贪心策略获得卫星执行的倾角控制行为;
S4:执行倾角控制行为,得到下一时刻卫星倾角;根据下一时刻卫星倾角、南北保持策略奖励函数得到奖励;所述南北保持策略奖励函数采用公式1:
Figure BDA0003937348540000021
其中,rt为当前时刻卫星进行倾角控制行为所获得的奖励,s0为标称轨道的倾角,t为当前时刻,t0为与当前时刻距离最近的预期轨控时刻,Δrs为倾角保持圆半径,st+1为当前时刻的下一时刻卫星倾角,当前时刻的下一时刻卫星倾角差为|st+1-s0|;当当前时刻的下一时刻卫星倾角不在倾角保持圆的范围内就给予惩罚;
S5:将当前时刻卫星倾角、卫星执行的倾角控制行为、奖励和下一时刻卫星倾角作为一组卫星组合状态数据组存入经验池中;
S6:从经验池中取出多组卫星组合状态数据组,根据神经网络权值参数计算每个卫星组合状态数据组的目标值;
S7:根据损失函数计算误差,并更新神经网络权值参数;
S8:根据值函数更新Q值;将下一时刻卫星倾角作为当前时刻卫星倾角;
S9:重复执行步骤S3-S8,执行步骤S3-S8的次数等于该组卫星训练状态数据组的预期轨控次数;
S10:重复执行步骤S2-S9,直到输入完所有卫星训练状态数据组的数据。
根据本发明一示例实施方式,步骤S1中,初始化模型包括定义损失函数。
根据本发明一示例实施方式,模型的输入为卫星倾角,模型的输出为卫星执行倾角控制行为后产生的回报值(Q值)。
根据本发明一示例实施方式,步骤S1中,卫星倾角为在轨道上的二维倾角,该二维倾角根据卫星轨道倾角和升交点赤经获得:
s=(ix,iy);
Figure BDA0003937348540000031
其中,s表示卫星在轨道上的二维倾角,i表示卫星轨道倾角,Ω表示升交点赤经。
根据本发明一示例实施方式,所述二维倾角为矢量数据。
根据本发明一示例实施方式,步骤S3中,初次循环时,当前时刻卫星倾角为初始时刻卫星倾角。
根据本发明一示例实施方式,步骤S3中,所述根据贪心策略获得卫星执行的倾角控制行为的方法包括:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为;第一指定概率和第二指定概率之和等于1。
根据本发明一示例实施方式,步骤S6中,所述根据神经网络权值参数计算每个卫星组合状态数据组的目标值的方法采用公式2:
Figure BDA0003937348540000032
其中,yj表示目标值,γ为折扣值,w为神经网络权值参数,
Figure BDA0003937348540000033
表示一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后的最大Q值,sj+1表示一组卫星组合状态数据组中的下一时刻卫星倾角,a表示一组卫星组合状态数据组中的当前时刻卫星执行的倾角控制行为,rj表示一组卫星组合状态数据组中的奖励。
根据本发明一示例实施方式,步骤S7中,所述损失函数采用公式3:
Figure BDA0003937348540000034
其中,yj表示目标值,w为神经网络权值参数,Q(sj,aj;w)表示一组卫星组合状态数据组中的当前时刻卫星执行倾角控制行为aj后的Q值,sj表示一组卫星组合状态数据组中的当前时刻卫星倾角,aj表示卫星执行的倾角控制行为,m为卫星组合状态数据组的数量。
根据本发明一示例实施方式,步骤S8中,所述根据值函数更新Q值的方法采用公式4:
Q(st,at)←Q(st,at)+α[rt+γmax Q(st+1,at)-Q(st,at)] (4);
其中,箭头左侧的Q(st,at)表示更新后的当前时刻的卫星执行倾角控制行为at后的Q值,箭头右侧的Q(st,at)表示更新前的当前时刻的卫星执行倾角控制行为at后的Q值,Q(st+1,at)表示更新前的当前时刻的下一时刻卫星执行倾角控制行为at后的Q值,α为权重,γ为折扣值,st表示当前时刻卫星倾角,at表示当前时刻卫星执行的倾角控制行为,st+1表示当前时刻的下一时刻卫星倾角,rt表示奖励。
t表示当前时刻,t+1表示当前时刻的下一时刻。
作为本发明的第二个方面,提供了一种基于深度Q网络的卫星南北保持策略模型,采用所述基于深度Q网络的卫星南北保持策略模型的建模方法建立模型。
作为本发明的第三个方面,提供一种卫星南北保持最优策略的获取方法,采用所述基于深度Q网络的卫星南北保持策略模型的建模方法建立基于深度Q网络的卫星南北保持策略模型;
根据该模型得到最优策略;
所述根据该模型得到最优策略的方法采用公式5:
Figure BDA0003937348540000041
其中,π表示卫星进行倾角控制的策略,π*表示通过模型学习到的最优的倾角控制策略,即卫星在初始卫星倾角数据为s的情况下经过策略π*的倾角控制行为a下产生最大的回报。
作为本发明的第四个方面,提供一种电子设备,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现所述基于深度Q网络的卫星南北保持策略模型的建模方法。
作为本发明的第五个方面,提供一种计算机可读介质,其上存储有计算机程序,所述程序被处理器执行时实现所述基于深度Q网络的卫星南北保持策略模型的建模方法。
本发明的有益效果是:
本方案通过神经网络进行建模,利用当前卫星倾角数据进行深度强化学习和决策,不需要利用卫星在轨道运行过程中受到的各种摄动力进行复杂的建模,不仅能够得到最优的南北控制策略,而且能够减少卫星燃料的消耗,这对实际的航天应用具有重要的意义和价值。
附图说明
图1示意性示出了基于深度Q网络的卫星南北保持策略模型的建模方法的步骤图。
图2示意性示出了一种电子设备的框图。
图3示意性示出了一种计算机可读介质的框图。
具体实施方式
现在将参考附图更全面地描述示例实施例。然而,示例实施例能够以多种形式实施,且不应被理解为限于在此阐述的实施例;相反,提供这些实施例使得本申请将全面和完整,并将示例实施例的构思全面地传达给本领域的技术人员。在图中相同的附图标记表示相同或类似的部分,因而将省略对它们的重复描述。
此外,所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施例中。在下面的描述中,提供许多具体细节从而给出对本申请的实施例的充分理解。然而,本领域技术人员将意识到,可以实践本申请的技术方案而没有特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知方法、装置、实现或者操作以避免模糊本申请的各方面。
附图中所示的方框图仅仅是功能实体,不一定必须与物理上独立的实体相对应。即,可以采用软件形式来实现这些功能实体,或在一个或多个硬件模块或集成电路中实现这些功能实体,或在不同网络和/或处理器装置和/或微控制器装置中实现这些功能实体。
附图中所示的流程图仅是示例性说明,不是必须包括所有的内容和操作/步骤,也不是必须按所描述的顺序执行。例如,有的操作/步骤还可以分解,而有的操作/步骤可以合并或部分合并,因此实际执行的顺序有可能根据实际情况改变。
应理解,虽然本文中可能使用术语第一、第二、第三等来描述各种组件,但这些组件不应受这些术语限制。这些术语乃用以区分一组件与另一组件。因此,下文论述的第一组件可称为第二组件而不偏离本申请概念的教示。如本文中所使用,术语“及/或”包括相关联的列出项目中的任一个及一或多者的所有组合。
本领域技术人员可以理解,附图只是示例实施例的示意图,附图中的模块或流程并不一定是实施本申请所必须的,因此不能用于限制本申请的保护范围。
本方案基于深度学习较强的感知能力来从环境中获取观测信息,基于强化学习较强的决策能力获得预期回报值来评估顶座价值。整个学习过程可以描述为:在某一时刻,卫星与飞行的环境交互来获取观测信息,通过神经网络将当前状态信息映射为相应的动作(控制行为),环境对此动作做出反应,得到对应的奖励值和下一个观测信息,并把完整的交互信息存入经验池中。通过不断循环以上过程,最终可以得到实现目标的最优策略。
本方案所述的卫星为GEO三轴稳定卫星。地球同步轨道(GEO),是指人造卫星在地球赤道上空约36000公里处围绕地球运行的圆形轨道。因为卫星绕地球运行周期与地球自转同步,卫星与地球之间处于相对静止的状态,因此在这种轨道上运行的卫星被简称为“同步卫星”,又称为“静止卫星”或“固定卫星”。三轴稳定就是卫星不旋转,本体在X、Y、Z三个方向上均稳定,换言之就是与地球保持一定的姿态关系。
深度Q网络(DeepQNetworks,简称DQN)算法是深度强化学习中的一种网络,是深度学习与Q学习的结合。由于其综合了强化学习和深度学习的优势,所以目前在各个领域都已经得到广泛的应用。
深度强化学习作为人工智能领域的一个新的研究热点,将深度学习和强化学习相结合,通过一种端对端的学习方式实现从原始输入到输出的直接控制和决策。由于深度学习是基于神经网络结构的,所以它对环境具有较强的感知能力,但是缺乏一定的决策控制能力;而强化学习恰好具有很强的行为决策能力。因此,深度强化学习将深度学习的感知能力和强化学习的决策能力相结合,优势互补,可以直接从高维原始数据中学习控制策略。深度强化学习方法自提出以来,在许多需要感知高维度原始输入数据和决策控制的任务中,已经取得了实质性的突破,由于深度学习端到端的学习优势,故深度强化学习可以解决建模难、规划难的问题。
作为本发明的第一个实施方式,提供一种基于深度Q网络的卫星南北保持策略模型的建模方法,如图1所示,包括以下步骤:
S1:初始化模型,获取多组卫星训练状态数据组,每组卫星训练状态数据组包括卫星的初始状态、多个预期轨控时刻和预期轨控次数;每个卫星的初始状态包括初始时刻卫星倾角。
初始化模型的方法包括:定义损失函数;初始化经验池的容量为N,经验池用于存储训练的样本;初始化网络模型的神经网络权值参数w;初始化模型网络的输入为卫星倾角s,计算出的网络输出为卫星执行倾角控制行为后的回报值Q。
卫星在某一时刻的运动状态可由开普勒轨道六根数表示:半长轴、偏心率、升交点赤经、近地点幅角、轨道倾角、平近点角。
在南北保持策略问题中,由于卫星在运行的过程中会受到日、月引力和地球非球形摄动的影响,会使得轨道倾角i发生漂移。卫星轨道平面在空间的方位通常用倾角i和升交点赤经Ω两个轨道要素来描述。但是在倾角很小的情况下,为了避免奇点,使用以下轨道要素来代替i和Ω:
Figure BDA0003937348540000071
卫星倾角的数据为卫星在轨道上的二维倾角,该二维倾角为矢量数据,卫星在轨道上的二维倾角矢量可表示为:
s=(ix,iy)。
因此,卫星在轨道上的二维倾角根据卫星轨道倾角和升交点赤经获得:
s=(ix,iy);
Figure BDA0003937348540000081
其中,s表示卫星在轨道上的二维倾角,i表示卫星轨道倾角,Ω表示升交点赤经。
多组卫星训练状态数据组组成数据集,该数据集中的卫星训练用的状态的数据大于或等于100组,卫星状态的数据越多,模型训练出来的结果越准确。
多组卫星训练状态数据组的数据为训练集的数据,可以采用仿真数据,也可以采用仿真数据和真实数据相结合。一个时间段内的时间线上包括多个时间点,每个时间点上卫星的状态都不同,在不同的时间点执行轨控策略会得到不同的效果。本方案的通过多组卫星训练状态数据组,每组卫星的初始时刻卫星倾角对应一个时间点的卫星倾角,每组卫星训练状态数据组的初始时刻对应的时间点均不同,即每组卫星训练状态数据组的初始时刻均不同。
S2:向模型中输入一组卫星训练状态数据组的初始时刻卫星倾角,得到初始时刻之后所有倾角控制行为及所对应输出的Q值。
初始时刻卫星执行倾角控制行为后,得到下一时刻卫星倾角。下一时刻卫星执行倾角控制行为后,得到下一时刻的下一时刻卫星倾角。以此类推,得到下多个时刻的倾角控制行为。
S3:获取当前时刻卫星倾角,根据贪心策略获得卫星执行的倾角控制行为。
初次循环时,当前时刻卫星倾角为初始时刻卫星倾角。
根据贪心策略获得卫星执行的倾角控制行为的方法包括:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为;第一指定概率和第二指定概率之和等于1。
如果第一指定概率大于第二指定概率,则根据贪心策略获得卫星执行的倾角控制行为的方法采用:卫星以第一指定概率随机选择倾角控制行为;
如果第二指定概率大于第一指定概率,则根据贪心策略获得卫星执行的倾角控制行为的方法采用:卫星以第二指定概率执行最大Q值所对应的控制行为;
如果第一指定概率等于第二指定概率,则根据贪心策略获得卫星执行的倾角控制行为的方法择其一:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为。
贪心策略为ε-greedy策略。
第一指定概率为ε,ε随迭代次数的增加而减小。
当前时刻卫星执行的倾角控制行为为at
S4:执行倾角控制行为,得到下一时刻卫星倾角;根据下一时刻卫星倾角、南北保持策略奖励函数得到奖励。
卫星的倾角保持的方法是设置以Δrs为半径的倾角保持圆,让倾角持续漂移,直至接近半径为Δrs的倾角保持圆的上边界时执行倾角机动,使倾角矢量跳变至倾角保持圆内的下边界,倾角机动方向基本沿日、月摄动的反方向。
卫星南北保持策略问题的目标是在保持轨道倾角的同时还要尽可能最小化燃料的消耗量,所以奖励策略的定义要同时考虑到倾角的变动和燃料的消耗。卫星初始质量一定,总燃耗取决于每次控制速度增量绝对值之和,而后者又取决于每次控制产生的倾角矢量改变量绝对值之和。
假设对于卫星控制频率固定(即固定经过一段时间后进行轨控),期望本次控制量可以保证在下一次控制时候倾角在保持圆内且控制量尽可能小,即在下一时刻轨道的倾角不仅满足在倾角保持圆内,而且还要尽可能接近标称轨道(理论轨道)。所以某一时刻(t时刻)倾角的变化量决定外推的某一时刻的下一时刻(t时刻的下一时刻)的轨道倾角的状态。基于这个目的,设计了t时刻的奖励策略。标称轨道的倾角为s0,倾角保持圆半径为Δrs,t时刻的奖励策略即为南北保持策略奖励函数采用公式1:
Figure BDA0003937348540000091
其中,rt为当前时刻卫星进行倾角控制行为所获得的奖励,s0为标称轨道的倾角,t为当前时刻,t0为与当前时刻距离最近的预期的轨控时刻,Δrs为倾角保持圆半径,st+1为当前时刻的下一时刻的倾角,当前时刻的下一时刻的倾角差为|st+1-s0|;当当前时刻的下一时刻的倾角不在倾角保持圆的范围内就给予惩罚。
t+1时刻为t时刻(当前时刻)的下一时刻。外推的t+1时刻的倾角为t时刻的下一时刻的倾角。
S5:将当前时刻卫星倾角、卫星执行的倾角控制行为、奖励和下一时刻卫星倾角作为一组卫星组合状态数据组存入经验池中。
S6:从经验池中取出多组卫星组合状态数据组,根据神经网络权值参数计算每个卫星组合状态数据组的目标值。
卫星组合状态数据组的数量为m,m为大于0的自然数,且m小于卫星训练状态数据组的数量。m组卫星组合状态数据组为小批量的卫星组合状态数据组。卫星组合状态数据组的数量根据卫星训练状态数据组的数量确定。
所述根据神经网络权值参数计算每个卫星组合状态数据组的目标值的方法采用公式2:
Figure BDA0003937348540000101
其中,yj表示目标值,γ为折扣值(衰减因子),w为神经网络权值参数,
Figure BDA0003937348540000102
表示一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后的最大Q值,sj+1表示一组卫星组合状态数据组中的下一时刻卫星倾角,a表示卫星执行的倾角控制行为,rj表示一组卫星组合状态数据组中的奖励。
使任务终止即为模型收敛或迭代完成。当sj+1为模型收敛或迭代完成时,yi等于rj;当sj+1不为模型收敛或迭代完成时,yi等于
Figure BDA0003937348540000103
模型收敛的条件为:损失函数计算出来的误差在指定范围内。
迭代完成的条件为:执行完所有步骤。
S7:根据损失函数计算误差,并更新神经网络权值参数。
损失函数采用公式3:
Figure BDA0003937348540000104
其中,yj表示目标值,w为神经网络权值参数,Q(sj,aj;w)表示一组卫星组合状态数据组中的当前时刻卫星执行倾角控制行为aj后的Q值,sj表示一组卫星组合状态数据组中的当前时刻卫星倾角,aj表示卫星执行的倾角控制行为,rj表示一组卫星组合状态数据组中的奖励;m为卫星组合状态数据组的数量。
误差即为损失函数采用公式3的计算结果。
神经网络权值参数通过随机梯度下降法(SGD)进行更新。
rt、at、st、st+1表示卫星训练状态数据组的数据集中的样本,rj、aj、sj、sj+1表示经验池中的样本。
步骤S5-S7对模型的参数进行调整,能够让模型的计算准确度更高。
S8:根据值函数更新Q值,将下一时刻卫星倾角作为当前时刻卫星倾角。
根据值函数更新Q值的方法采用公式4:
Q(st,at)←Q(st,at)+α[rt+γmax Q(st+1,at)-Q(st,at)] (4);
其中,箭头左侧的Q(st,at)表示更新后的当前时刻的卫星执行倾角控制行为at后的Q值,箭头右侧的Q(st,at)表示更新前的当前时刻的卫星执行倾角控制行为at后的Q值,Q(st+1,at)表示更新前的当前时刻的下一时刻的卫星执行倾角控制行为at后的Q值,α为权重,γ为折扣值(衰减因子),st表示当前时刻卫星倾角,at表示当前时刻卫星执行的倾角控制行为,st+1表示当前时刻的下一时刻卫星倾角,rt表示奖励。
其中α和γ的范围都在0~1之间。
S9:重复执行步骤S3-S8,执行步骤S3-S8的次数等于该组卫星训练状态数据组的预期轨控次数。
S10:重复执行步骤S2-29,直到输入完所有卫星训练状态数据组的数据。
本方案的建模方法以卫星倾角数据作为神经网络模型的输入,产生的回报值作为输出,采用深度神经网络,不需要利用卫星在轨道运行过程中受到的各种摄动力进行复杂的建模,直接采用深度强化学习来进行学习和决策,不仅能够得到最优的南北控制策略,而且能够减少卫星燃料的消耗这对实际的航天应用具有重要的意义和价值。
根据本发明的第二个具体实施方式,本发明提供一种基于深度Q网络的卫星南北保持策略模型,采用第一个实施方式的基于深度Q网络的卫星南北保持策略模型的建模方法建立模型。
根据本发明的第三个具体实施方式,本发明提供一种卫星南北保持最优策略的获取方法,采用第一个实施方式的基于深度Q网络的卫星南北保持策略模型的建模方法建立基于深度Q网络的卫星南北保持策略模型;
根据该模型得到最优策略。
根据该模型得到最优策略的方法采用公式5:
Figure BDA0003937348540000121
其中,π表示卫星进行倾角控制的策略,π*表示通过模型学习到的最优的倾角控制策略,即卫星在初始卫星倾角数据为s的情况下经过策略π*的倾角控制行为a下产生最大的回报。
根据本发明的第四个具体实施方式,本发明提供一种电子设备,如图2所示,图2是根据一示例性实施例示出的一种电子设备的框图。
下面参照图2来描述根据本申请的这种实施方式的电子设备200。图2显示的电子设备200仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。
如图2所示,电子设备200以通用计算设备的形式表现。电子设备200的组件可以包括但不限于:至少一个处理单元210、至少一个存储单元220、连接不同系统组件(包括存储单元220和处理单元210)的总线230、显示单元240等。
其中,所述存储单元存储有程序代码,所述程序代码可以被所述处理单元210执行,使得所述处理单元210执行本说明书中描述的根据本申请各种示例性实施方式的步骤。例如,所述处理单元210可以执行如图1中所示的步骤。
所述存储单元220可以包括易失性存储单元形式的可读介质,例如随机存取存储单元(RAM)3201和/或高速缓存存储单元2202,还可以进一步包括只读存储单元(ROM)2203。
所述存储单元220还可以包括具有一组(至少一个)程序模块2205的程序/实用工具2204,这样的程序模块2205包括但不限于:操作系统、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。
总线230可以为表示几类总线结构中的一种或多种,包括存储单元总线或者存储单元控制器、外围总线、图形加速端口、处理单元或者使用多种总线结构中的任意总线结构的局域总线。
电子设备200也可以与一个或多个外部设备200’(例如键盘、指向设备、蓝牙设备等)通信,使得用户能与该电子设备200交互的设备通信,和/或该电子设备300能与一个或多个其它计算设备进行通信的任何设备(例如路由器、调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口250进行。并且,电子设备200还可以通过网络适配器260与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。网络适配器260可以通过总线230与电子设备200的其它模块通信。应当明白,尽管图中未示出,可以结合电子设备200使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
通过以上的实施方式的描述,本领域的技术人员易于理解,这里描述的示例实施方式可以通过软件实现,也可以通过软件结合必要的硬件的方式来实现。
因此,根据本发明的第五个具体实施方式,本发明提供一种计算机可读介质。如图3所示,根据本发明实施方式的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中或网络上,包括若干指令以使得一台计算设备(可以是个人计算机、服务器、或者网络设备等)执行根据本发明实施方式的上述方法。
所述软件产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以为但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
所述计算机可读存储介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了可读程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。可读存储介质还可以是可读存储介质以外的任何可读介质,该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。可读存储介质上包含的程序代码可以用任何适当的介质传输,包括但不限于无线、有线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言的任意组合来编写用于执行本发明操作的程序代码,所述程序设计语言包括面向对象的程序设计语言—诸如Java、C++等,还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中,远程计算设备可以通过任意种类的网络,包括局域网(LAN)或广域网(WAN),连接到用户计算设备,或者,可以连接到外部计算设备(例如利用因特网服务提供商来通过因特网连接)。
上述计算机可读介质承载有一个或者多个程序,当上述一个或者多个程序被一个该设备执行时,使得该计算机可读介质实现第一个实施方式的功能。
本领域技术人员可以理解上述各模块可以按照实施例的描述分布于装置中,也可以进行相应变化唯一不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。
通过以上的实施例的描述,本领域的技术人员易于理解,这里描述的示例实施例可以通过软件实现,也可以通过软件结合必要的硬件的方式来实现。因此,根据本发明实施例的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中或网络上,包括若干指令以使得一台计算设备(可以是个人计算机、服务器、移动终端、或者网络设备等)执行根据本发明实施例的方法。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于深度Q网络的卫星南北保持策略模型的建模方法,其特征在于,包括以下步骤:
S1:初始化模型,获取多组卫星训练状态数据组,每组卫星训练状态数据组包括卫星的初始状态、多个预期轨控时刻和预期轨控次数;每个卫星的初始状态包括初始时刻卫星倾角;
S2:向模型中输入一组卫星训练状态数据组的初始时刻卫星倾角,得到初始时刻之后所有倾角控制行为及所对应输出的Q值;
S3:获取当前时刻卫星倾角,根据贪心策略获得卫星执行的倾角控制行为;
S4:执行倾角控制行为,得到下一时刻卫星倾角;根据下一时刻卫星倾角、南北保持策略奖励函数得到奖励;所述南北保持策略奖励函数采用公式1:
Figure FDA0003937348530000011
其中,rt为当前时刻卫星进行倾角控制行为所获得的奖励,s0为标称轨道的倾角,t为当前时刻,t0为与当前时刻距离最近的预期轨控时刻,Δrs为倾角保持圆半径,st+1为当前时刻的下一时刻卫星倾角,当前时刻的下一时刻卫星倾角差为|st+1-s0|;当当前时刻的下一时刻卫星倾角不在倾角保持圆的范围内就给予惩罚;
S5:将当前时刻卫星倾角、卫星执行的倾角控制行为、奖励和下一时刻卫星倾角作为一组卫星组合状态数据组存入经验池中;
S6:从经验池中取出多组卫星组合状态数据组,根据神经网络权值参数计算每个卫星组合状态数据组的目标值;
S7:根据损失函数计算误差,并更新神经网络权值参数;
S8:根据值函数更新Q值;将下一时刻卫星倾角作为当前时刻卫星倾角;
S9:重复执行步骤S3-S8,执行步骤S3-S8的次数等于该组卫星训练状态数据组的预期轨控次数;
S10:重复执行步骤S2-S9,直到输入完所有卫星训练状态数据组的数据。
2.根据权利要求1所述的基于深度Q网络的卫星南北保持策略模型的建模方法,其特征在于,步骤S1中,卫星倾角的数据为卫星在轨道上的二维倾角,该二维倾角根据卫星轨道倾角和升交点赤经获得;
s=(ix,iy);
Figure FDA0003937348530000021
其中,s表示卫星在轨道上的二维倾角,i表示卫星轨道倾角,Ω表示升交点赤经。
3.根据权利要求1所述的基于深度Q网络的卫星南北保持策略模型的建模方法,其特征在于,步骤S3中,所述根据贪心策略获得卫星执行的倾角控制行为的方法包括:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为;第一指定概率和第二指定概率之和等于1。
4.根据权利要求1所述的基于深度Q网络的卫星南北保持策略模型的建模方法,其特征在于,步骤S6中,所述根据神经网络权值参数计算每个卫星组合状态数据组的目标值的方法采用公式2:
Figure FDA0003937348530000022
其中,yj表示目标值,γ为折扣值,w为神经网络权值参数,
Figure FDA0003937348530000023
表示一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后的最大Q值,sj+1表示一组卫星组合状态数据组中的下一时刻卫星倾角,a表示一组卫星组合状态数据组中的当前时刻卫星执行的倾角控制行为,rj表示一组卫星组合状态数据组中的奖励。
5.根据权利要求1所述的基于深度Q网络的卫星南北保持策略模型的建模方法,其特征在于,步骤S7中,所述损失函数采用公式3:
Figure FDA0003937348530000031
其中,yj表示目标值,w为神经网络权值参数,Q(sj,aj;w)表示一组卫星组合状态数据组中的当前时刻卫星执行倾角控制行为aj后的Q值,sj表示一组卫星组合状态数据组中的当前时刻卫星倾角,aj表示卫星执行的倾角控制行为,m为卫星组合状态数据组的数量。
6.根据权利要求1所述的基于深度Q网络的卫星南北保持策略模型的建模方法,其特征在于,步骤S8中,所述根据值函数更新Q值的方法采用公式4:
Q(st,at)←Q(st,at)+α[rt+γmax Q(st+1,at)-Q(st,at)] (4);
其中,箭头左侧的Q(st,at)表示更新后的当前时刻的卫星执行倾角控制行为at后的Q值,箭头右侧的Q(st,at)表示更新前的当前时刻的卫星执行倾角控制行为at后的Q值,Q(st+1,at)表示更新前的当前时刻的下一时刻卫星执行倾角控制行为at后的Q值,α为权重,γ为折扣值,st表示当前时刻卫星倾角,at表示当前时刻卫星执行的倾角制行为,st+1表示当前时刻的下一时刻卫星倾角,rt表示奖励。
7.一种基于深度Q网络的卫星南北保持策略模型,其特征在于,采用权利要求1-6中任一项所述的建模方法建立模型。
8.一种卫星南北保持最优策略的获取方法,其特征在于,根据权利要求1-6中任一项所述的建模方法建立基于深度Q网络的卫星南北保持策略模型;
根据该模型得到最优策略;
所述根据该模型得到最优策略的方法采用公式5:
Figure FDA0003937348530000032
其中,π表示卫星进行倾角控制的策略,π*表示通过模型学习到的最优的倾角控制策略,即卫星在初始卫星倾角为s的情况下经过策略π*的倾角控制行为a下产生最大的回报。
9.一种电子设备,其特征在于,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-6中任一所述的方法。
10.一种计算机可读介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现如权利要求1-6中任一所述的方法。
CN202211408067.XA 2022-11-10 2022-11-10 卫星南北保持策略模型的建模方法、系统、获取方法 Active CN115758704B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211408067.XA CN115758704B (zh) 2022-11-10 2022-11-10 卫星南北保持策略模型的建模方法、系统、获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211408067.XA CN115758704B (zh) 2022-11-10 2022-11-10 卫星南北保持策略模型的建模方法、系统、获取方法

Publications (2)

Publication Number Publication Date
CN115758704A true CN115758704A (zh) 2023-03-07
CN115758704B CN115758704B (zh) 2023-05-05

Family

ID=85369190

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211408067.XA Active CN115758704B (zh) 2022-11-10 2022-11-10 卫星南北保持策略模型的建模方法、系统、获取方法

Country Status (1)

Country Link
CN (1) CN115758704B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111563326A (zh) * 2020-04-29 2020-08-21 深圳星地孪生科技有限公司 卫星轨道运行参数学习系统、设备及系统运行方法
CN113815902A (zh) * 2021-11-24 2021-12-21 北京航天驭星科技有限公司 一种用于卫星轨道控制效果标定的方法
CN114362810A (zh) * 2022-01-11 2022-04-15 重庆邮电大学 一种基于迁移深度强化学习的低轨卫星跳波束优化方法
CN114889848A (zh) * 2022-03-30 2022-08-12 中国科学院软件研究所 一种卫星姿态的控制方法、控制装置、计算机设备和介质
CN114928401A (zh) * 2022-05-17 2022-08-19 重庆邮电大学 一种基于多智能体强化学习的leo星间链路的动态规划方法
CN114933028A (zh) * 2022-07-21 2022-08-23 北京航天驭星科技有限公司 双星轨控策略控制方法、装置、电子设备及存储介质
WO2022198748A1 (zh) * 2021-03-23 2022-09-29 北京航空航天大学 一种面向araim应用的低轨卫星增强系统星座配置优化方法
CN115173923A (zh) * 2022-07-04 2022-10-11 重庆邮电大学 一种低轨卫星网络能效感知路由优化方法和系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111563326A (zh) * 2020-04-29 2020-08-21 深圳星地孪生科技有限公司 卫星轨道运行参数学习系统、设备及系统运行方法
WO2022198748A1 (zh) * 2021-03-23 2022-09-29 北京航空航天大学 一种面向araim应用的低轨卫星增强系统星座配置优化方法
CN113815902A (zh) * 2021-11-24 2021-12-21 北京航天驭星科技有限公司 一种用于卫星轨道控制效果标定的方法
CN114362810A (zh) * 2022-01-11 2022-04-15 重庆邮电大学 一种基于迁移深度强化学习的低轨卫星跳波束优化方法
CN114889848A (zh) * 2022-03-30 2022-08-12 中国科学院软件研究所 一种卫星姿态的控制方法、控制装置、计算机设备和介质
CN114928401A (zh) * 2022-05-17 2022-08-19 重庆邮电大学 一种基于多智能体强化学习的leo星间链路的动态规划方法
CN115173923A (zh) * 2022-07-04 2022-10-11 重庆邮电大学 一种低轨卫星网络能效感知路由优化方法和系统
CN114933028A (zh) * 2022-07-21 2022-08-23 北京航天驭星科技有限公司 双星轨控策略控制方法、装置、电子设备及存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵义平;李峰;经姚翔;李大伟;: "GEO电推进卫星轨道漂移策略研究", 航天器工程 *

Also Published As

Publication number Publication date
CN115758704B (zh) 2023-05-05

Similar Documents

Publication Publication Date Title
Boain AB-Cs of sun-synchronous orbit mission design
Zuo et al. A case learning-based differential evolution algorithm for global optimization of interplanetary trajectory design
Carrara An open source satellite attitude and orbit simulator toolbox for Matlab
CN108021138B (zh) 一种地磁场模型简化设计方法
JP7250221B2 (ja) 準衛星軌道
CN110816896B (zh) 一种卫星星上简易轨道外推方法
Ocampo et al. Theoretical foundation of Copernicus: a unified system for trajectory design and optimization
CN115795816B (zh) 卫星东西保持策略模型的建模方法、模型、获取方法
CN115758707B (zh) 卫星东西保持策略模型的建模方法、系统、获取方法
CN115795817B (zh) 卫星东西保持策略模型的建模方法、系统、获取方法
CN113310496A (zh) 一种确定月地转移轨道的方法及装置
CN115758704B (zh) 卫星南北保持策略模型的建模方法、系统、获取方法
CN115865166B (zh) 卫星南北保持策略模型的建模方法、系统、获取方法
CN115758705B (zh) 卫星南北保持策略模型的建模方法、系统、获取方法
CN115865167B (zh) 卫星南北保持策略模型的建模方法、系统、获取方法
CN115892516B (zh) 卫星相对相位保持策略模型的建模方法、模型、获取方法
CN115758706B (zh) 卫星东西保持策略模型的建模方法、模型、获取方法
Golikov THEONA—a numerical-analytical theory of motion of artificial satellites of celestial bodies
CN115806062B (zh) 卫星相对相位保持策略模型的建模方法、系统、获取方法
CN115806060B (zh) 卫星相对相位保持策略模型的建模方法、模型、获取方法
CN115806061B (zh) 卫星相对相位保持策略模型的建模方法、模型、获取方法
KR102056719B1 (ko) 정지궤도 관측 자기장 자료를 이용한 우주기상 요소와 Dst 지수 예측 시스템 및 방법
Miñan et al. Manoeuvre planning algorithm for satellite formations using mean relative orbital elements
Bowen On-board orbit determination and 3-axis attitude determination for picosatellite applications
Wei et al. Redesign of high-precision reference orbit for interferometric SAR satellite with injection error

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
GR01 Patent grant
GR01 Patent grant