CN115758705B - 卫星南北保持策略模型的建模方法、系统、获取方法 - Google Patents
卫星南北保持策略模型的建模方法、系统、获取方法 Download PDFInfo
- Publication number
- CN115758705B CN115758705B CN202211408069.9A CN202211408069A CN115758705B CN 115758705 B CN115758705 B CN 115758705B CN 202211408069 A CN202211408069 A CN 202211408069A CN 115758705 B CN115758705 B CN 115758705B
- Authority
- CN
- China
- Prior art keywords
- satellite
- inclination angle
- moment
- strategy
- south
- 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.)
- Active
Links
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及航空航天领域,提供一种基于DoubleDQN的卫星南北保持策略模型的建模方法、系统、获取方法、设备、介质,建模方法包括:S1:初始化模型;S2:得到所有倾角控制行为及Q值;S3:获得当前时刻卫星倾角和卫星执行的倾角控制行为;S4:执行倾角控制行为,得到下一时刻卫星倾角和下次实际控制时刻;根据下次实际控制时刻、南北保持策略奖励函数得到奖励;S5:将卫星组合状态数据组存入经验池中;S6:计算目标值;S7:计算误差,并更新当前神经网络权值参数;S8:更新Q值;S9:重复执行步骤S3‑S8;每重复执行指定迭代次数的步骤S3‑S8后,更新目标神经网络权值参数;S10:重复执行步骤S2‑S9。本方案能够得到最优的决策策略,减少卫星燃料的消耗。
Description
技术领域
本发明涉及航空航天技术领域,特别是涉及一种基于Double DQN的卫星南北保持策略模型的建模方法、系统、获取方法、设备、介质。
背景技术
随着人类航天活动的不断开展,越来越多的遥感卫星为人们的日常生活提供帮助。
GEO卫星在运行的过程中会受到日、月引力和地球非球形摄动的影响,从而导致在南北纬度方向上发生漂移,所以对GEO三轴稳定卫星南北位置保持(倾角保持)控制对于航天领域有至关重要的作用。
现有的技术首先通过对卫星在轨道运行过程中受到的地球形状、日月引力、太阳光压等各种摄动力的分析建立轨道机动时的动力学模型,然后制定南北保持的短期和长期策略。这种方法通过卫星在轨道运行过程中受到的各种摄动力进行复杂的建模,然而由于空间受力的复杂性以及卫星本身参数的不确定性,导致无法对卫星进行精确建模,参数多、计算复杂,进而影响卫星倾角控制的精度,可能会消耗更多的燃料。并且,现有强化学习方法无法解决状态和动作空间维度很高的问题。
因此,亟需开发一种基于Double DQN的卫星南北保持策略模型的建模方法、系统、获取方法、设备、介质,降低建模难度,精确计算南北保持策略。
发明内容
本发明的目的是提供一种基于Double DQN的卫星南北保持策略模型的建模方法、系统、获取方法、设备、介质,对GEO三轴稳定卫星进行南北位置保持时不需要进行复杂的建模,不需要考虑空间受力的复杂性以及卫星本身参数的不确定性,强化学习具有很强的行为决策能力,能够得到最优的决策策略,减少卫星燃料的消耗。
为解决上述技术问题,作为本发明的一个方面,提供了一种基于Double DQN的卫星南北保持策略模型的建模方法,包括以下步骤:
S1:初始化模型,获取多组卫星训练状态数据组,每组卫星训练状态数据组包括卫星的初始状态、多个预期轨控时刻和预期轨控次数;每个卫星的初始状态包括初始时刻卫星倾角;
S2:向模型中输入一组卫星训练状态数据组的初始时刻卫星倾角,得到初始时刻之后所有倾角控制行为及所对应输出的Q值;
S3:获取当前时刻卫星倾角,根据贪心策略获得卫星执行的倾角控制行为;
S4:执行倾角控制行为,得到下一时刻卫星倾角和下次实际控制时刻;根据下次实际控制时刻、南北保持策略奖励函数得到奖励;所述南北保持策略奖励函数采用公式1:
其中,rt为当前时刻卫星进行倾角控制行为所获得的奖励,Δt为距离下次实际控制时刻的时间间隔,根据当前时刻及下次实际控制时刻获得;t0为与当前时刻距离最近的预期轨控时刻;t为当前时刻;
S5:将当前时刻卫星倾角、卫星执行的倾角控制行为、奖励和下一时刻卫星倾角作为一组卫星组合状态数据组存入经验池中;
S6:从经验池中取出多组卫星组合状态数据组,根据目标神经网络权值参数计算每个卫星组合状态数据组的目标值;
S7:根据损失函数计算误差,并更新当前神经网络权值参数;
S8:根据值函数更新Q值;将下一时刻卫星倾角作为当前时刻卫星倾角;
S9:重复执行步骤S3-S8,执行步骤S3-S8的次数等于该组卫星训练状态数据组的轨控次数;每重复执行指定迭代次数的步骤S3-S8后,根据当前神经网络权值参数更新目标神经网络权值参数;
S10:重复执行步骤S2-S9,直到输入完所有卫星训练状态数据组的数据。
根据本发明一示例实施方式,步骤S1中,初始化模型包括定义损失函数。
根据本发明一示例实施方式,模型的输入为卫星倾角,模型的输出为卫星执行倾角控制行为的回报值(Q值)。
根据本发明一示例实施方式,步骤S1中,卫星倾角为在轨道上的二维倾角,该二维倾角根据卫星轨道倾角和升交点赤经获得:
s=(ix,iy);
其中,s表示卫星在轨道上的二维倾角,i表示卫星轨道倾角,Ω表示升交点赤经。
根据本发明一示例实施方式,所述二维倾角为矢量数据。
根据本发明一示例实施方式,步骤S3中,初次循环时,当前时刻卫星倾角为初始时刻卫星倾角。
根据本发明一示例实施方式,步骤S3中,所述根据贪心策略获得卫星执行的倾角控制行为的方法包括:卫星以第一指定概率随机选择下一时刻的倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为;第一指定概率和第二指定概率之和等于1。
根据本发明一示例实施方式,步骤S6中,根据目标神经网络权值参数计算每个卫星组合状态数据组的目标值的方法采用公式2:
其中,yj表示目标值,γ为折扣值,θ’为目标神经网络权值参数,θ为当前神经网络权值参数,表示一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后得到的Q值,表示选择能使一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后得到的最大Q值所对应的倾角控制行为,sj+1表示一组卫星组合状态数据组中的下一时刻卫星倾角,a表示卫星执行的倾角控制行为,rj表示一组卫星组合状态数据组中的奖励。
根据本发明一示例实施方式,步骤S7中,所述损失函数采用公式3:
其中,yj表示目标值,θ为当前神经网络权值参数,Q(sj,aj;θ)表示一组卫星组合状态数据组中的当前时刻卫星执行倾角控制行为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表示当前时刻的下一时刻。
作为本发明的第二个方面,提供了一种基于Double DQN的卫星南北保持策略系统,采用所述基于Double DQN的卫星南北保持策略模型的建模方法建立系统。
作为本发明的第三个方面,提供一种卫星南北保持最优策略的获取方法,采用所述基于Double DQN的卫星南北保持策略模型的建模方法建立基于Double DQN的卫星南北保持策略模型;
根据该模型得到最优策略;
所述根据该模型得到最优策略的方法采用公式5:
其中,π表示卫星进行倾角控制的策略,π*表示通过模型学习到的最优的倾角控制策略,即卫星在初始时刻卫星倾角为s的情况下经过策略π*的倾角控制行为a下产生最大的回报。
作为本发明的第四个方面,提供一种电子设备,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现所述基于Double DQN的卫星南北保持策略模型的建模方法。
作为本发明的第五个方面,提供一种计算机可读介质,其上存储有计算机程序,所述程序被处理器执行时实现所述基于Double DQN的卫星南北保持策略模型的建模方法。
本发明的有益效果是:
本方案通过神经网络进行建模,利用当前卫星倾角数据进行深度强化学习和决策,不需要利用卫星在轨道运行过程中受到的各种摄动力进行复杂的建模,不仅能够得到最优的南北控制策略,而且能够减少卫星燃料的消耗,这对实际的航天应用具有重要的意义和价值。
附图说明
图1示意性示出了基于Double DQN的卫星南北保持策略模型的建模方法的步骤图。
图2示意性示出了一种电子设备的框图。
图3示意性示出了一种计算机可读介质的框图。
具体实施方式
现在将参考附图更全面地描述示例实施例。然而,示例实施例能够以多种形式实施,且不应被理解为限于在此阐述的实施例;相反,提供这些实施例使得本申请将全面和完整,并将示例实施例的构思全面地传达给本领域的技术人员。在图中相同的附图标记表示相同或类似的部分,因而将省略对它们的重复描述。
此外,所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施例中。在下面的描述中,提供许多具体细节从而给出对本申请的实施例的充分理解。然而,本领域技术人员将意识到,可以实践本申请的技术方案而没有特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知方法、装置、实现或者操作以避免模糊本申请的各方面。
附图中所示的方框图仅仅是功能实体,不一定必须与物理上独立的实体相对应。即,可以采用软件形式来实现这些功能实体,或在一个或多个硬件模块或集成电路中实现这些功能实体,或在不同网络和/或处理器装置和/或微控制器装置中实现这些功能实体。
附图中所示的流程图仅是示例性说明,不是必须包括所有的内容和操作/步骤,也不是必须按所描述的顺序执行。例如,有的操作/步骤还可以分解,而有的操作/步骤可以合并或部分合并,因此实际执行的顺序有可能根据实际情况改变。
应理解,虽然本文中可能使用术语第一、第二、第三等来描述各种组件,但这些组件不应受这些术语限制。这些术语乃用以区分一组件与另一组件。因此,下文论述的第一组件可称为第二组件而不偏离本申请概念的教示。如本文中所使用,术语“及/或”包括相关联的列出项目中的任一个及一或多者的所有组合。
本领域技术人员可以理解,附图只是示例实施例的示意图,附图中的模块或流程并不一定是实施本申请所必须的,因此不能用于限制本申请的保护范围。
本方案基于深度学习较强的感知能力来从环境中获取观测信息,基于强化学习较强的决策能力获得预期回报值来评估顶座价值。整个学习过程可以描述为:在某一时刻,卫星与飞行的环境交互来获取观测信息,通过神经网络将当前状态信息映射为相应的动作(倾角控制行为),环境对此动作做出反应,得到对应的奖励值和下一个观测信息,并把完整的交互信息存入经验池中。通过不断循环以上过程,最终可以得到实现目标的最优策略。
本方案所述的卫星为GEO三轴稳定卫星。地球同步轨道(GEO),是指人造卫星在地球赤道上空约36000公里处围绕地球运行的圆形轨道。因为卫星绕地球运行周期与地球自转同步,卫星与地球之间处于相对静止的状态,因此在这种轨道上运行的卫星被简称为“同步卫星”,又称为“静止卫星”或“固定卫星”。三轴稳定就是卫星不旋转,本体在X、Y、Z三个方向上均稳定,换言之就是与地球保持一定的姿态关系。
深度Q网络(Deep Q Networks,简称DQN)算法是深度强化学习中的一种网络,是深度学习与Q学习的结合。由于其综合了强化学习和深度学习的优势,所以目前在各个领域都已经得到广泛的应用。
深度强化学习作为人工智能领域的一个新的研究热点,将深度学习和强化学习相结合,通过一种端对端的学习方式实现从原始输入到输出的直接控制和决策。由于深度学习是基于神经网络结构的,所以它对环境具有较强的感知能力,但是缺乏一定的决策控制能力;而强化学习恰好具有很强的行为决策能力。因此,深度强化学习将深度学习的感知能力和强化学习的决策能力相结合,优势互补,可以直接从高维原始数据中学习控制策略。深度强化学习方法自提出以来,在许多需要感知高维度原始输入数据和决策控制的任务中,已经取得了实质性的突破,由于深度学习端到端的学习优势,故深度强化学习可以解决建模难、规划难的问题。
DQN算法在计算目标值和计算当前值用的是同一个网络,即目标值的计算使用了当前要训练的Q网络中的参数,而又用目标值来更新网络的参数,这样二者循环依赖,不利于算法的收敛。相比较于DQN,Nature DQN增加了目标网络,通过双网络结构来减少目标Q值计算与待更新Q网络参数之间的依赖关系,综合了强化学习和深度学习的优势,从而极大地提高了DQN算法的稳定性。
Nature DQN通过使用两个独立的但完全相同的Q网络(一个网络作为当前Q网络,另一个作为目标Q网络),减少计算目标网络的目标值和当前网络参数之间的相关性。通过复制当前网络的权重参数给目标网络,使得当前网络每隔一定步长C更新目标网络,双网络结构让目标Q值在一段时间里使保持不变,降低了计算目标Q值和当前网络参数的相关性,进而提高了算法收敛性和稳定性。
DQN和Nature DQN算法在进行值函数优化目标时,动作(倾角控制行为)选择和策略评估都是基于相同的值函数来实现的。无论是DQN还是Nature DQN,都无法克服Q-Learning本身固有的缺陷——过估计,从而使得估计的值函数比真实值大。
过估计是指估计的值函数比真实值函数要大,其根源主要在于Q-Learning中的最大化操作,即目标值其中的max操作使得估计的值函数比值函数的真实值大(注:对于真实的策略来说并在给定的状态下并不是每次都选择使得Q值最大的倾角控制行为,因为一般真实的策略都是随机性策略,所以在这里目标值直接选择倾角控制行为最大的Q值往往会导致目标值要高于真实值)。
不同于这两种算法,DoubleDQN(简称DDQN)则是基于NatureDQN双网络架构,将动作选择和策略评估分开进行,使用当前神经网络权值参数θ来选择最优动作,使用目标神经网络权值参数θ'来评估最优倾角控制行为,从而解决DQN和NatureDQN算法过度估计的问题。DDQN算法和NatureDQN的算法步骤的区别就在目标Q值的计算方式不同。DDQN能够估计出更加准确的Q值,可获得更稳定有效的策略。
作为本发明的第一个实施方式,提供一种基于DoubleDQN的卫星南北保持策略模型的建模方法,如图1所示,包括以下步骤:
S1:初始化模型,获取多组卫星训练状态数据组,每组卫星训练状态数据组包括卫星的初始状态、多个预期轨控时刻和预期轨控次数;每个卫星的初始状态包括初始时刻卫星倾角。
初始化模型的方法包括:定义损失函数;初始化经验池的容量为N,经验池用于存储训练的样本;初始化网络模型的当前神经网络权值参数θ和目标神经网络权值参数θ′,θ′=θ;初始化任务训练的指定迭代次数T1;初始化网络的输入为卫星倾角s,计算出的网络输出为卫星执行倾角控制行为后的回报值Q。
卫星在某一时刻的运动状态可由开普勒轨道六根数表示:半长轴、偏心率、升交点赤经、近地点幅角、轨道倾角、平近点角。
在南北保持策略问题中,由于卫星在运行的过程中会受到日、月引力和地球非球形摄动的影响,会使得轨道倾角i发生漂移。卫星轨道平面在空间的方位通常用倾角i和升交点赤经Ω两个轨道要素来描述。但是在倾角很小的情况下,为了避免奇点,使用以下轨道要素来代替i和Ω:
卫星倾角的数据为卫星在轨道上的二维倾角,该二维倾角为矢量数据,卫星在轨道上的二维倾角矢量可表示为:
s=(ix,iy)。
因此,卫星在轨道上的二维倾角根据卫星轨道倾角和升交点赤经获得:
s=(ix,iy);
其中,s表示卫星在轨道上的二维倾角,i表示卫星轨道倾角,Ω表示升交点赤经。
多组卫星训练状态数据组组成数据集,该数据集中的卫星状态的数据大于或等于100组,卫星状态的数据越多,模型训练出来的结果越准确。
多组卫星训练状态数据组的数据为训练集的数据,可以采用仿真数据,也可以采用仿真数据和真实数据相结合。一个时间段内的时间线上包括多个时间点,每个时间点上卫星的状态都不同,在不同的时间点执行轨控策略会得到不同的效果。本方案的通过多组卫星训练状态数据组,每组卫星的初始时刻卫星倾角对应一个时间点的卫星倾角,每组卫星训练状态数据组的初始时刻对应的时间点均不同,即每组卫星训练状态数据组的初始时刻均不同。
S2:向模型中输入一组卫星训练状态数据组的初始时刻卫星倾角,得到初始时刻之后所有倾角控制行为及所对应输出的Q值。
初始时刻卫星执行倾角控制行为后,得到下一时刻卫星倾角。下一时刻卫星执行倾角控制行为后,得到下一时刻的下一时刻卫星倾角。以此类推,得到下多个时刻的倾角控制行为。
S3:获取当前时刻卫星倾角,根据贪心策略获得卫星执行的倾角控制行为。
初次循环时,当前时刻卫星倾角为初始时刻卫星倾角。
根据贪心策略获得卫星执行的控制行为的方法包括:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为;第一指定概率和第二指定概率之和等于1。
如果第一指定概率大于第二指定概率,则根据贪心策略获得卫星执行的倾角控制行为的方法采用:卫星以第一指定概率随机选择倾角控制行为;
如果第二指定概率大于第一指定概率,则根据贪心策略获得卫星执行的倾角控制行为的方法采用:卫星以第二指定概率执行最大Q值所对应的倾角控制行为;
如果第一指定概率等于第二指定概率,则根据贪心策略获得卫星执行的倾角控制行为的方法择其一:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为。
贪心策略为ε-greedy策略。
第一指定概率为ε,ε随迭代次数的增加而减小。
当前时刻卫星执行的倾角控制行为为at。
S4:执行倾角控制行为,得到下一时刻卫星倾角及下次实际控制时刻;根据下次实际控制时刻、南北保持策略奖励函数得到奖励。下次实际控制时刻即为下次实际需要进行倾角控制行为的时刻。
由于日、月引力和地球非球形摄动的影响,GEO三轴稳定卫星轨道倾角和升交点赤经会不断变化,每年倾角矢量漂移率的变化为0.75°~0.95°。为使卫星的南北位置误差保持在规定的范围内,必须定期实施南北位置保持(轨道倾角控制)。卫星的南北位置保持是在轨运行时最复杂的测控工作,在燃料量和倾角保持精度要求一定的前提下,选择目标轨道倾角,使卫星两次控制时间间隔最长,且消耗燃料最省是南北保持控制的关键。
卫星的倾角保持是将倾角控制在偏差为Δi的范围内,让倾角持续漂移,直至接近偏差为Δi的倾角边界时执行倾角机动。考虑到每次对卫星进行控制都需要耗费大量的人工作业,所以控制的次数越少,风险就越低,这就要求对卫星进行倾角保持的周期越长。
在每次控制消耗的燃料等量的前提下,t时刻(当前时刻)倾角的变化量决定了下一次进行轨道倾角控制的时间。卫星标称轨道的倾角为i,在某一时刻对卫星进行控制后,根据外推的星历得到卫星倾角在超出i±Δi范围的时刻。基于卫星倾角保持周期长的目的,设计了t时刻的奖励策略。t时刻的奖励策略即为南北保持策略奖励函数采用公式1:
其中,rt为当前时刻卫星进行倾角控制行为所获得的奖励,Δt为距离下次实际控制时刻的时间间隔,Δt根据当前时刻及下次实际控制时刻获得,即Δt=t1-t;t为当前时刻,t1为下次实际控制时刻,t0为与当前时刻距离最近的预期轨控时刻。
下次实际控制时刻为外推出来的要进行的倾角控制的时刻,即倾角超出预定范围的时刻。
t+1时刻为t时刻(当前时刻)的下一时刻。t时刻执行倾角控制行为后,t+1时刻的倾角发生相应变化,t时刻的下多个时刻倾角也发生相应变化。Δt时间越长,保持周期越长,奖励越高。
S5:将当前时刻卫星倾角、卫星执行的倾角控制行为、奖励和下一时刻卫星倾角作为一组卫星组合状态数据组存入经验池中。
S6:从经验池中取出多组卫星组合状态数据组,根据目标神经网络权值参数计算每个卫星组合状态数据组的目标值。
卫星组合状态数据组的数量为m,m为大于0的自然数,且m小于卫星训练状态数据组的数量。m组卫星组合状态数据组为小批量的卫星组合状态数据组。卫星组合状态数据组的数量根据卫星训练状态数据组的数量确定。
所述根据目标神经网络权值参数计算每个卫星组合状态数据组的目标值的方法采用公式2:
其中,yj表示目标值,γ为折扣值(衰减因子),θ为目标神经网络权值参数,θ为当前神经网络权值参数,表示一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后得到的Q值,表示选择能使一组卫星组合状态数据组中的下一时刻卫星执行倾角控制行为a后得到的最大Q值所对应的倾角控制行为,sj+1表示一组卫星组合状态数据组中的下一时刻卫星倾角,a表示卫星执行的倾角控制行为,rj表示一组卫星组合状态数据组中的奖励。
指的是一个倾角控制行为,该倾角控制行为是该组卫星组合状态数据组的所有下一时刻卫星执行倾角控制行为中的一个,并且执行该倾角控制行为后,得到的Q值比该组卫星组合状态数据组的所有下一时刻卫星执行倾角控制行为得到的Q值大。
公式2将动作(倾角控制行为)选择和策略评估分开进行,使用当前神经网络权值参数θ来选择最优动作,使用目标神经网络权值参数θ′来评估最优倾角控制行为,从而解决DQN和Nature DQN算法过度估计的问题。
模型收敛的条件为:损失函数计算出来的误差在指定范围内。
迭代完成的条件为:执行完所有步骤。
S7:根据损失函数计算误差,并更新当前神经网络权值参数。
还根据目标值计算误差。
损失函数采用公式3:
其中,yj表示目标值,θ为当前神经网络权值参数,Q(sj,aj;θ)表示一组卫星组合状态数据组中的当前时刻卫星执行倾角控制行为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表示当前时刻卫星倾角(t时刻卫星倾角),at表示当前时刻卫星执行的倾角控制行为,st+1表示当前时刻的下一时刻卫星倾角,rt表示奖励。
其中α和γ的范围都在0~1之间。
S9:重复执行步骤S3-S8,执行步骤S3-S8的次数等于该组卫星训练状态数据组的轨控次数;每重复执行指定迭代次数的步骤S3-S8后,根据当前神经网络权值参数更新目标神经网络权值参数。
完成指定迭代次数的迭代后,将目标神经网络权值参数更新为当前神经网络权值参数。通过这种方式,减少计算目标网络的目标值和当前网络参数之间的相关性。通过复制当前网络的权重参数给目标网络,使得当前网络每隔一定步长C更新目标网络,双网络结构让目标Q值在一段时间里使保持不变,降低了计算目标Q值和当前网络参数的相关性,进而提高了算法收敛性和稳定性。
S10:重复执行步骤S2-S9,直到输入完所有卫星训练状态数据组的数据。本方案的建模方法以卫星倾角数据作为神经网络模型的输入,产生的回报值作为输出,采用DoubleDQN神经网络,不需要利用卫星在轨道运行过程中受到的各种摄动力进行复杂的建模,直接采用深度强化学习来进行学习和决策,基于DQN算法进行改进适用于训练大规模的神经网络且极大提高了DQN算法的稳定性,不仅能够得到最优的南北控制策略,而且能够减少卫星燃料的消耗这对实际的航天应用具有重要的意义和价值。
根据本发明的第二个具体实施方式,本发明提供一种基于Double DQN的卫星南北保持策略系统,采用第一个实施方式的基于Double DQN的卫星南北保持策略模型的建模方法建立系统。
根据本发明的第三个具体实施方式,本发明提供一种卫星南北保持最优策略的获取方法,采用第一个实施方式的基于Double DQN的卫星南北保持策略模型的建模方法建立基于Double DQN的卫星南北保持策略模型;
根据该模型得到最优策略。
根据该模型得到最优策略的方法采用公式5:
其中,π表示卫星进行倾角控制的策略,π*表示通过模型学习到的最优的倾角控制策略,即卫星在初始时刻卫星倾角为s的情况下经过策略π*的倾角控制行为a下产生最大的回报。
根据本发明的第四个具体实施方式,本发明提供一种电子设备,如图2所示,图2是根据一示例性实施例示出的一种电子设备的框图。
下面参照图2来描述根据本申请的这种实施方式的电子设备200。图2显示的电子设备200仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。
如图2所示,电子设备200以通用计算设备的形式表现。电子设备200的组件可以包括但不限于:至少一个处理单元210、至少一个存储单元220、连接不同系统组件(包括存储单元220和处理单元210)的总线230、显示单元240等。
其中,所述存储单元存储有程序代码,所述程序代码可以被所述处理单元210执行,使得所述处理单元210执行本说明书中描述的根据本申请各种示例性实施方式的步骤。例如,所述处理单元210可以执行如图1中所示的步骤。
所述存储单元220可以包括易失性存储单元形式的可读介质,例如随机存取存储单元(RAM)2201和/或高速缓存存储单元2202,还可以进一步包括只读存储单元(ROM)2203。
所述存储单元220还可以包括具有一组(至少一个)程序模块2205的程序/实用工具2204,这样的程序模块2205包括但不限于:操作系统、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。
总线230可以为表示几类总线结构中的一种或多种,包括存储单元总线或者存储单元控制器、外围总线、图形加速端口、处理单元或者使用多种总线结构中的任意总线结构的局域总线。
电子设备200也可以与一个或多个外部设备200’(例如键盘、指向设备、蓝牙设备等)通信,使得用户能与该电子设备200交互的设备通信,和/或该电子设备200能与一个或多个其它计算设备进行通信的任何设备(例如路由器、调制解调器等等)通信。这种通信可以通过输入/输出(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.一种基于Double DQN的卫星南北保持策略模型的建模方法,其特征在于,包括以下步骤:
S1:初始化模型,获取多组卫星训练状态数据组,每组卫星训练状态数据组包括卫星的初始状态、多个预期轨控时刻和预期轨控次数;每个卫星的初始状态包括初始时刻卫星倾角;
S2:向模型中输入一组卫星训练状态数据组的初始时刻卫星倾角,得到初始时刻之后所有倾角控制行为及所对应输出的Q值;
S3:获取当前时刻卫星倾角,根据贪心策略获得卫星执行的倾角控制行为;
S4:执行倾角控制行为,得到下一时刻卫星倾角和下次实际控制时刻;根据下次实际控制时刻、南北保持策略奖励函数得到奖励;所述南北保持策略奖励函数采用公式1:
其中,rt为当前时刻卫星进行倾角控制行为所获得的奖励,Δt为距离下次实际控制时刻的时间间隔,根据当前时刻及下一次实际控制时刻获得;t0为与当前时刻距离最近的预期轨控时刻;t为当前时刻;
S5:将当前时刻卫星倾角、卫星执行的倾角控制行为、奖励和下一时刻卫星倾角作为一组卫星组合状态数据组存入经验池中;
S6:从经验池中取出多组卫星组合状态数据组,根据目标神经网络权值参数计算每个卫星组合状态数据组的目标值;
S7:根据损失函数计算误差,并更新当前神经网络权值参数;
S8:根据值函数更新Q值;将下一时刻卫星倾角作为当前时刻卫星倾角;
S9:重复执行步骤S3-S8,执行步骤S3-S8的次数等于该组卫星训练状态数据组的预期轨控次数;每重复执行指定迭代次数的步骤S3-S8后,根据当前神经网络权值参数更新目标神经网络权值参数;
S10:重复执行步骤S2-S9,直到输入完所有卫星训练状态数据组的数据。
3.根据权利要求1所述的基于Double DQN的卫星南北保持策略模型的建模方法,其特征在于,步骤S3中,所述根据贪心策略获得卫星执行的倾角控制行为的方法包括:卫星以第一指定概率随机选择倾角控制行为、或以第二指定概率执行最大Q值所对应的倾角控制行为;第一指定概率和第二指定概率之和等于1。
6.根据权利要求1所述的基于Double DQN的卫星南北保持策略模型的建模方法,其特征在于,步骤S8中,所述根据值函数更新Q值的方法采用公式4:
Q(st,at)←Q(st,at)+α[rt+γmaxQ(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.一种基于Double DQN的卫星南北保持策略系统,其特征在于,采用权利要求1-6中任一项所述的建模方法建立系统。
9.一种电子设备,其特征在于,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-6中任一所述的方法。
10.一种计算机可读介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现如权利要求1-6中任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211408069.9A CN115758705B (zh) | 2022-11-10 | 2022-11-10 | 卫星南北保持策略模型的建模方法、系统、获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211408069.9A CN115758705B (zh) | 2022-11-10 | 2022-11-10 | 卫星南北保持策略模型的建模方法、系统、获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115758705A CN115758705A (zh) | 2023-03-07 |
CN115758705B true CN115758705B (zh) | 2023-05-05 |
Family
ID=85369193
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211408069.9A Active CN115758705B (zh) | 2022-11-10 | 2022-11-10 | 卫星南北保持策略模型的建模方法、系统、获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115758705B (zh) |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112668235B (zh) * | 2020-12-07 | 2022-12-09 | 中原工学院 | 基于离线模型预训练学习的ddpg算法的机器人控制方法 |
CN114613169B (zh) * | 2022-04-20 | 2023-02-28 | 南京信息工程大学 | 一种基于双经验池dqn的交通信号灯控制方法 |
-
2022
- 2022-11-10 CN CN202211408069.9A patent/CN115758705B/zh active Active
Non-Patent Citations (1)
Title |
---|
采用偏心率和倾角矢量联合隔离法实现双星共位;马煦;孙海燕;黄显安;李永财;;电讯技术(第07期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115758705A (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3688675B1 (en) | Distributional reinforcement learning for continuous control tasks | |
Boain | AB-Cs of sun-synchronous orbit mission design | |
US11801948B2 (en) | Machine learning system and method for orbital trajectory planning | |
Sullivan et al. | Using reinforcement learning to design a low-thrust approach into a periodic orbit in a multi-body system | |
Chan et al. | Autonomous imaging and mapping of small bodies using deep reinforcement learning | |
Zuo et al. | A case learning-based differential evolution algorithm for global optimization of interplanetary trajectory design | |
CN115758707B (zh) | 卫星东西保持策略模型的建模方法、系统、获取方法 | |
Guzzetti | Reinforcement learning and topology of orbit manifolds for stationkeeping of unstable symmetric periodic orbits | |
CN115795817B (zh) | 卫星东西保持策略模型的建模方法、系统、获取方法 | |
CN115795816B (zh) | 卫星东西保持策略模型的建模方法、模型、获取方法 | |
Bonasera et al. | Designing impulsive station-keeping maneuvers near a sun-earth l2 halo orbit via reinforcement learning | |
Herrmann et al. | Reinforcement learning for the agile earth-observing satellite scheduling problem | |
CN115758705B (zh) | 卫星南北保持策略模型的建模方法、系统、获取方法 | |
CN115865167B (zh) | 卫星南北保持策略模型的建模方法、系统、获取方法 | |
CN115865166B (zh) | 卫星南北保持策略模型的建模方法、系统、获取方法 | |
CN115758706B (zh) | 卫星东西保持策略模型的建模方法、模型、获取方法 | |
CN115758704B (zh) | 卫星南北保持策略模型的建模方法、系统、获取方法 | |
Yang et al. | A station-keeping control method for GEO spacecraft based on autonomous control architecture | |
CN115806062B (zh) | 卫星相对相位保持策略模型的建模方法、系统、获取方法 | |
CN115892516B (zh) | 卫星相对相位保持策略模型的建模方法、模型、获取方法 | |
CN115806060B (zh) | 卫星相对相位保持策略模型的建模方法、模型、获取方法 | |
CN115806061B (zh) | 卫星相对相位保持策略模型的建模方法、模型、获取方法 | |
Herrmann et al. | Autonomous spacecraft tasking using monte carlo tree search methods | |
Adams | Theory and Applications of Gram-Scale Spacecraft | |
Geng et al. | Impulsive guidance of optimal pursuit with conical imaging zone for the evader |
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 |