CN104536302B - 一种卫星仿真***和方法 - Google Patents
一种卫星仿真***和方法 Download PDFInfo
- Publication number
- CN104536302B CN104536302B CN201410497678.5A CN201410497678A CN104536302B CN 104536302 B CN104536302 B CN 104536302B CN 201410497678 A CN201410497678 A CN 201410497678A CN 104536302 B CN104536302 B CN 104536302B
- Authority
- CN
- China
- Prior art keywords
- simulation
- satellite
- coordinate system
- subsystem
- orbit
- 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.)
- Expired - Fee Related
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 420
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims description 39
- 230000005358 geomagnetic field Effects 0.000 claims description 28
- 230000014509 gene expression Effects 0.000 claims description 24
- 230000007613 environmental effect Effects 0.000 abstract description 7
- 238000005457 optimization Methods 0.000 abstract 1
- 108091092878 Microsatellite Proteins 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 230000005484 gravity Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013079 data visualisation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供一种卫星仿真***和方法,涉及卫星仿真领域,该***操作灵活、便于改进、易于优化,该***包括:仿真控制子***、环境模拟子***、视景轨道仿真子***和卫星模拟子***,其中,该仿真控制子***分别与该环境模拟子***、该视景轨道仿真子***和该卫星模拟子***连接,该环境模拟子***与该视景轨道仿真子***连接,本发明实施例用于卫星仿真。
Description
技术领域
本发明涉及卫星仿真领域,尤其涉及一种卫星仿真***和方法。
背景技术
近现代卫星的发展,以微小卫星领域最为激烈,尤其是近地轨道和种地轨道的微小卫星的研制和开发,是各个航天大国的空间技术发展的焦点和热点。微小卫星质量轻,体积小,功耗低,寿命长,能满足人们各种功能需求,如地球资源探测、卫星通讯、军事侦察和遥感测绘等等。微小卫星成本低,使用灵活,既能单独使用,也能数十颗甚至更多连成一个网络,在军事方面能实现覆盖式无死角的监视,在通讯方面能实现全球高质量全方位的通讯,在农林业方面,能更精准的预报灾害,减小损失。微小卫星具有如此广阔的应用前景,因此,对于微小卫星的研究是当今乃至未来卫星空间技术领域发展的重要方向之一。
卫星姿态与轨道控制***是卫星上各大保障***中最复杂、技术性最突出的***之一,对卫星上的有效载荷任务的完成起着非常重要的作用,如果该***出现故障或异常,将直接导致卫星不能完成既定任务,甚至使其变成一个空间垃圾。随着新一代姿态与轨道控制***变得日益先进与复杂,面对不断新增的新型号和新任务,对地面卫星仿真***的灵活性、通用性与可扩展性以及实现测试数据可视化的要求不断提高。
现有技术中,主要通过地面卫星仿真***完成卫星仿真的过程,该地面卫星仿真***通过获取磁场信息和轨道信息并通过计算完成卫星仿真,但是,地面卫星仿真***在卫星仿真过程中的计算量较大,因此造成该***在卫星仿真过程中消耗的时间较长。
发明内容
本发明的实施例提供一种卫星仿真***和方法,该***能够减少卫星仿真的时间,从而提高卫星仿真的效率。
为达到上述目的,本发明的实施例采用如下技术方案:
一方面,提供一种卫星仿真***,包括:仿真控制子***、环境模拟子***、视景轨道仿真子***和卫星模拟子***,其中,所述仿真控制子***分别与所述环境模拟子***、所述视景轨道仿真子***和所述卫星模拟子***连接,所述环境模拟子***与所述视景轨道仿真子***连接,
所述仿真控制子***,用于生成模拟卫星,并获取当前时刻所述模拟卫星对应的磁场信息和轨道信息,并将所述磁场信息和所述轨道信息发送至与所述磁场信息和轨道信息对应的所述卫星模拟子***;
所述卫星模拟子***,用于接收所述仿真控制子***发送的所述磁场信息和所述轨道信息,并根据所述磁场信息和所述轨道信息以及配置的力矩得到控制信息,并将所述控制信息发送至所述仿真控制子***;
所述仿真控制子***,还用于根据所述卫星模拟子***发送的控制信息得到所述模拟卫星在所述下一时刻的姿态信息,并根据所述下一时刻的姿态信息调整所述模拟卫星的姿态;
所述环境模拟子***,用于获取所述模拟卫星的磁场信息,并将所述磁场信息发送至所述仿真控制子***;
所述视景轨道仿真子***,用于获取所述模拟卫星的轨道信息,并将所述轨道信息发送至所述仿真控制子***。
进一步地,所述仿真控制子***具体用于,根据所述控制信息得到所述模拟卫星当前时刻的横滚角、俯仰角和偏航角,并根据所述模拟卫星当前时刻的横滚角、俯仰角和偏航角得到所述模拟卫星下一时刻的横滚角、俯仰角和偏航角,并根据所述下一时刻的横滚角、俯仰角和偏航角调整所述模拟卫星的姿态。
更进一步地,所述视景轨道仿真子***还用于,展示所述模拟卫星的姿态和对应的轨道信息,并将所述轨道信息发送至所述环境模拟子***;
所述环境模拟子***具体用于,根据所述轨道信息得到磁场信息;
所述环境模拟子***还用于,将所述磁场信息发送至所述视景轨道仿真子***;
所述视景轨道仿真子***还用于,展示所述环境模拟子***发送的磁场信息。
更进一步地,所述仿真控制子***具体用于,根据所述当前时刻的横滚角、所述俯仰角和所述偏航角通过公式:
得到所述模拟卫星当前时刻的姿态四元数;
其中,Q为所述模拟卫星当前时刻的姿态四元数,
为所述模拟卫星当前时刻的横滚角,θ为所述模拟卫星当前时刻的俯仰角,ψ为所述模拟卫星当前时刻的偏航角;
所述仿真控制子***具体用于,根据所述姿态四元数通过公式:
得到所述模拟卫星从轨道坐标系到本体坐标系的姿态旋转矩阵;其中,为所述姿态旋转矩阵;
所述仿真控制子***具体用于,根据所述控制信息得到配置的力矩和所述模拟卫星的转动惯量以及所述模拟卫星在当前时刻所述本体坐标系相对于惯性坐标系的角速度,并根据所述力矩和所述转动惯量以及所述旋转矩阵通过公式:
得到下一时刻所述本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式,
其中,为下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式,I为所述模拟卫星的转动惯量,为当前时刻所述本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影,Tb为所述力矩在本体坐标系中的投影;
所述仿真控制子***具体用于,对所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式积分得到所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度,并根据所述控制信息得到轨道坐标系相对于惯性坐标系的旋转矩阵,并根据所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度和所述轨道坐标系相对于惯性坐标系的旋转矩阵通过公式:
得到下一时刻所述模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影;
其中,为下一时刻所述模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影,为所述姿态旋转矩阵,为下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中投影,为所述模拟卫星在轨道坐标系相对于惯性坐标系的角速度在所述轨道坐标系中的投影,ωo为所述模拟卫星的轨道角速度;
所述仿真控制子***具体用于,根据所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影和所述模拟卫星的转动惯量通过方程:
得到q0的微分表达式和的微分表达式,并对所述q0的微分表达式和的微分表达式分别积分得到所述模拟卫星下一时刻的姿态四元数,并根据所述模拟卫星下一时刻的姿态四元数得到所述模拟卫星下一时刻的横滚角、俯仰角和偏航角;
其中,为q0的微分表达式,为的微分表达式,为下一时刻所述模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影, I为所述模拟卫星的转动惯量。
更进一步地,所述环境模拟子***具体用于,根据所述轨道信息得到经度和余纬,并根据所述经度和余纬通过公式:
得到地球磁位势;
其中,V为所述地球磁位势,RE为地球半径,r为所述模拟卫星到地心的距离,λ为所述经度,为所述余纬,和分别为n次m阶高阶系数,为n次m阶缔合勒让德函数;
所述环境模拟子***具体用于,根据所述地球磁位势通过公式:
分别得到在地心赤道坐标系中X轴、Y轴和Z轴的地磁场强度;
其中,为X轴的地磁场强度,为Y轴的地磁场强度,为Z轴的地磁场强度;
所述环境模拟子***具体用于,通过公式
将在地心赤道坐标系中的地磁场强度转换为在地心轨道坐标系中的地磁场强度;
其中,为在地心轨道坐标系中的地磁场强度,为在地心赤道坐标系中的地磁场强度,Cz(ω)为在Z轴方向上旋转α角度的旋转矩阵,α为近地点幅角,Cx(i)为在X轴方向上旋转i角度的旋转矩阵,i为轨道倾角,Cz(Ω-β)为在Z轴方向上旋转角度Ω-β的旋转矩阵,Ω为升交点赤经,β为所述下一时刻的俯仰角;
所述环境模拟子***具体用于,通过公式
将在地心轨道坐标系中的地磁场强度转换为在轨道系坐标中的地磁场强度;
其中,为在轨道系坐标中的地磁场强度,为在地心轨道坐标系中的地磁场强度,Cx(π/2)为在X轴方向上旋转π/2角度的旋转矩阵,CZ(f+π/2)为在Z轴方向上旋转(f+π/2)角度的旋转矩阵,f为真近点角。
更进一步地,所述仿真控制子***与所述卫星模拟子***通过RS232总线连接,所述仿真控制子***与所述环境模拟子***和所述视景轨道仿真子***分别通过以太网连接,所述环境模拟子***和所述视景轨道仿真子***通过以太网连接。
另一方面,提供一种卫星仿真的方法,包括:
仿真控制子***生成模拟卫星,并获取所述环境模拟子***发送的所述模拟卫星当前时刻对应的磁场信息和所述视景轨道仿真子***发送所述模拟卫星当前时刻对应的轨道信息;
所述仿真控制子***将所述磁场信息和所述轨道信息发送至与所述磁场信息和轨道信息对应的所述卫星模拟子***;
所述卫星模拟子***根据所述磁场信息和所述轨道信息以及配置的力矩得到控制信息,并将所述控制信息发送至所述仿真控制子***;
所述仿真控制子***根据所述卫星模拟子***发送的控制信息得到所述模拟卫星在所述下一时刻的姿态信息,并根据所述下一时刻的姿态信息调整所述模拟卫星的姿态。
进一步地,所述仿真控制子***根据所述卫星模拟子***发送的控制信息得到所述模拟卫星在所述下一时刻的姿态信息,并根据所述下一时刻的姿态信息调整所述模拟卫星的姿态包括:
所述仿真控制子***根据所述控制信息得到所述模拟卫星当前时刻的横滚角、俯仰角和偏航角,并根据所述模拟卫星当前时刻的横滚角、俯仰角和偏航角得到所述模拟卫星下一时刻的横滚角、俯仰角和偏航角,并根据所述下一时刻的横滚角、俯仰角和偏航角调整所述模拟卫星的姿态。
更进一步地,在所述根据所述下一时刻的姿态信息调整所述模拟卫星的姿态后,所述方法还包括:
所述视景轨道仿真子***展示所述模拟卫星调整后的姿态和对应的轨道信息,并将所述模拟卫星调整后的姿态对应的轨道信息发送至所述环境模拟子***;
所述环境模拟子***根据所述模拟卫星调整后的姿态对应的轨道信息得到所述模拟卫星调整后的姿态对应的磁场信息,并将所述模拟卫星调整后的姿态对应的磁场信息发送至所述视景轨道仿真子***;
所述视景轨道仿真子***展示所述模拟卫星调整后的姿态对应的磁场信息。
本发明提供一种卫星仿真***和方法,该***包括仿真控制子***、环境模拟子***、视景轨道仿真子***和卫星模拟子***,其中,该仿真控制子***分别与该环境模拟子***、该视景轨道仿真子***和该卫星模拟子***连接,该环境模拟子***与该视景轨道仿真子***连接,该***主要模拟真实空间环境,该卫星仿真***采用分布式仿真,将整个仿真***的功能分配到各个子模块,能够很好的进行***负载的分配,通过数据链对数据的分发与维护实现仿真***有条不紊的运行,该仿真***涉及到的多个模块,根据自己的任务与载荷的要求,合理分配***资源,使得***能够以最快的速度完成该***的任务,弥补单个仿真***实时性差、资源紧张、同步性差等固有的不足之处,因此,本***具有资源利用率高,实时性好,仿真***模块化程度高,同步性强,可维护等优点,同时,***实时性和同步性的增强,使得仿真***能够极大的模拟出真实的空间环境,大大提高了***的仿真效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种卫星仿真***的结构示意图;
图2为本发明实施例提供的一种卫星仿真方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种卫星仿真***,如图1所示,该***包括:
仿真控制子***10、环境模拟子***11、视景轨道仿真子***12和卫星模拟子***13,其中,该仿真控制子***10分别与该环境模拟子***11、该视景轨道仿真子***12和该卫星模拟子***13连接,该环境模拟子***13与该视景轨道仿真子***12连接,
具体地,该仿真控制子***与该卫星模拟子***通过RS232总线连接,该仿真控制子***与该环境模拟子***和该视景轨道仿真子***分别通过以太网连接,该环境模拟子***和该视景轨道仿真子***通过以太网连接。
其中,该仿真控制子***和该环境模拟子***以及该视景轨道仿真子***之间通过TCP/IP协议进行数据传输。
需要说明的是,该仿真控制子***通过控制指令控制该环境模拟子***和该视景轨道仿真子***以及卫星模拟子***的开启。
示例地,该***中各子***之间的IP分配可以如下表:
为了保证数据传输的通用性,各子***间传输的数据格式可以为:
上述IP分配以及数据格式只是举例说明,不作限定。
该仿真控制子***10,用于生成模拟卫星,并获取当前时刻该模拟卫星对应的磁场信息和轨道信息,并将该磁场信息和该轨道信息发送至与该磁场信息和轨道信息对应的该卫星模拟子***。
需要说明的是,该仿真控制子***需要将获取的磁场信息和轨道信息转换为该卫星模拟子***对应的信息,并发送至该卫星模拟子***。
在本发明一种可能的实现方式中,该仿真控制子***接收该环境模拟子***发送的该当前时刻的磁场信息和该视景轨道仿真子***发送的该当前时刻的轨道信息。
其中,该磁场信息可以是地球磁位势,该轨道信息可以是经度和余纬。
该卫星模拟子***13,用于接收该仿真控制子***10发送的该磁场信息和该轨道信息,并根据该磁场信息和该轨道信息以及配置的力矩得到控制信息,并将该控制信息发送至该仿真控制子***。
示例地,该卫星模拟子***可以采用MSP430-5438A芯片作为主控芯片。
该仿真控制子***10,还用于根据该卫星模拟子***13发送的控制信息得到该模拟卫星在下一时刻的姿态信息,并根据该下一时刻的姿态信息调整该模拟卫星的姿态。
该环境模拟子***12,用于获取该模拟卫星的磁场信息,并将该磁场信息发送至该仿真控制子***10。
该视景轨道仿真子***12,用于获取该模拟卫星的轨道信息,并将该轨道信息发送至该仿真控制子***。
其中,该视景轨道仿真子***包括STK卫星仿真软件,STK特有的轨道模拟模块可以生成多种轨道信息。
在本发明一种可能的实现方式中,生成的轨道信息如下表所示:
用户还可以通过该视景轨道仿真子***配置轨道场景参数,该轨道场景参数包括场景名称、起止时间、时间***、场景历元、动画历元、轨道外推和轨道历元以及卫星名称等。
进一步地,该仿真控制子***10具体用于,根据该控制信息得到该模拟卫星当前时刻的横滚角、俯仰角和偏航角,并根据该模拟卫星当前时刻的横滚角、俯仰角和偏航角得到该模拟卫星下一时刻的横滚角、俯仰角和偏航角,并根据该下一时刻的横滚角、俯仰角和偏航角调整该模拟卫星的姿态。
具体地,该仿真控制子***10具体用于,根据该当前时刻的横滚角、该俯仰角和该偏航角通过公式:
得到该模拟卫星当前时刻的姿态四元数;
其中,Q为该模拟卫星当前时刻的姿态四元数,
为该模拟卫星当前时刻的横滚角,θ为该模拟卫星当前时刻的俯仰角,ψ为该模拟卫星当前时刻的偏航角。
该仿真控制子***10具体用于,根据该姿态四元数通过公式:
得到该模拟卫星从轨道坐标系到本体坐标系的姿态旋转矩阵;其中,为该姿态旋转矩阵;
该仿真控制子***10具体用于,根据该控制信息得到配置的力矩和该模拟卫星的转动惯量以及该模拟卫星在当前时刻该本体坐标系相对于惯性坐标系的角速度,并根据该力矩和该转动惯量以及该旋转矩阵通过公式:
得到下一时刻该本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式。
其中,为下一时刻该模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式,I为该模拟卫星的转动惯量,为当前时刻该本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影,Tb为该力矩在本体坐标系中的投影。
需要说明的是,该力矩可以是控制力矩,重力梯度力矩,大气干扰力矩等,当配置的力矩包括多个力矩时,则上述公式中的Tb实际为所有力矩在本体坐标系中的投影之和,例如,若该配置的力矩包括控制力矩,重力梯度力矩,大气干扰力矩,则上述
具体为:
其中,为控制力矩在本体坐标系中的投影,为重力梯度力矩在本体坐标系中的投影,为大气干扰力矩在本体坐标系中的投影;这里配置的力矩包括控制力矩,重力梯度力矩,大气干扰力矩只是举例说明,还可以包括其他力矩,本发明对此不作限定。
该仿真控制子***10具体用于,对该下一时刻该模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式积分得到该下一时刻该模拟卫星在本体坐标系相对于惯性坐标系的角速度,并根据该控制信息得到轨道坐标系相对于惯性坐标系的旋转矩阵,并根据该下一时刻该模拟卫星在本体坐标系相对于惯性坐标系的角速度和该轨道坐标系相对于惯性坐标系的旋转矩阵通过公式:
得到下一时刻该模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影;
其中,为下一时刻该模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影,为该姿态旋转矩阵,为下一时刻该模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中投影,为该模拟卫星在轨道坐标系相对于惯性坐标系的角速度在该轨道坐标系中的投影, ωo为该模拟卫星的轨道角速度。
该仿真控制子***10具体用于,根据该下一时刻该模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影和该模拟卫星的转动惯量通过方程:
得到q0的微分表达式和的微分表达式,并对该q0的微分表达式和的微分表达式分别积分得到该模拟卫星下一时刻的姿态四元数,并根据该模拟卫星下一时刻的姿态四元数得到该模拟卫星下一时刻的横滚角、俯仰角和偏航角。
其中,为q0的微分表达式,为的微分表达式,为下一时刻该模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影, I为该模拟卫星的转动惯量。
进一步地,该视景轨道仿真子***12还用于,展示该模拟卫星的姿态和对应的轨道信息,并将该轨道信息发送至该环境模拟子***。
该环境模拟子***11具体用于,根据该轨道信息得到磁场信息。
该环境模拟子***11还用于,将该磁场信息发送至该视景轨道仿真子***。
该视景轨道仿真子***12还用于,展示该环境模拟子***发送的磁场信息。
具体地,该环境模拟子***11具体用于,根据该轨道信息得到经度和余纬,并根据该经度和余纬通过公式:
得到地球磁位势;
其中,V为该地球磁位势,RE为地球半径,r为该模拟卫星到地心的距离,λ为该经度,为该余纬,和分别为n次m阶高阶系数,为n次m阶缔合勒让德函数;
该环境模拟子***11具体用于,根据该地球磁位势通过公式:
分别得到在地心赤道坐标系中X轴、Y轴和Z轴的地磁场强度;
其中,为X轴的地磁场强度,为Y轴的地磁场强度,为Z轴的地磁场强度;
该环境模拟子***11具体用于,通过公式
将在地心赤道坐标系中的地磁场强度转换为在地心轨道坐标系中的地磁场强度;
其中,为在地心轨道坐标系中的地磁场强度,为在地心赤道坐标系中的地磁场强度,Cz(ω)为在Z轴方向上旋转α角度的旋转矩阵,α为近地点幅角,Cx(i)为在X轴方向上旋转i角度的旋转矩阵,i为轨道倾角,Cz(Ω-β)为在Z轴方向上旋转角度Ω-β的旋转矩阵,Ω为升交点赤经,β为该下一时刻的俯仰角;
该环境模拟子***11具体用于,通过公式
将在地心轨道坐标系中的地磁场强度转换为在轨道系坐标中的地磁场强度;
其中,为在轨道系坐标中的地磁场强度,为在地心轨道坐标系中的地磁场强度,Cx(π/2)为在X轴方向上旋转π/2角度的旋转矩阵,CZ(f+π/2)为在Z轴方向上旋转(f+π/2)角度的旋转矩阵,f为真近点角。
需要说明的是,环境模拟子***在根据下一时刻的姿态信息调整该模拟卫星的姿态后,该视景轨道仿真子***展示该模拟卫星调整后的姿态和对应的轨道信息,并将该模拟卫星调整后的姿态对应的轨道信息发送至该环境模拟子***,该环境模拟子***根据该模拟卫星调整后的姿态对应的轨道信息得到该模拟卫星调整后的姿态对应的磁场信息,并将该模拟卫星调整后的姿态对应的磁场信息发送至该视景轨道仿真子***,该视景轨道仿真子***展示该模拟卫星调整后的姿态对应的磁场信息,同时,该环境模拟子***还将该模拟卫星调整后的姿态对应的磁场信息发送至仿真控制子***,该视景轨道仿真子***将该模拟卫星调整后的姿态对应的轨道信息发送至仿真控制子***,以便该仿真控制子***根据该模拟卫星调整后的姿态对应的磁场信息和该模拟卫星调整后的姿态对应的轨道信息继续通过卫星模拟子***获取下一时刻的姿态信息。
另外,在该卫星仿真***启动后,通过该视景轨道仿真子***配置初始轨道信息和初始磁场信息,该视景轨道仿真子***将该初始轨道信息和初始磁场信息发送至仿真控制子***,该仿真控制子***将该初始轨道信息和初始磁场信息发送至对应的卫星模拟子***,卫星模拟子***根据该初始轨道信息和初始磁场信息以及配置的力矩得到控制信息,并将该控制信息发送至仿真控制子***,以便该仿真控制子***根据该控制信息得到下一时刻的姿态信息。
该卫星仿真***采用分布式仿真,将整个仿真***的功能分配到各个子模块,能够很好的进行***负载的分配,通过数据链对数据的分发与维护实现仿真***有条不紊的运行,该仿真***涉及到的多个模块,根据自己的任务与载荷的要求,合理分配***资源,使得***能够以最快的速度完成该***的任务,弥补单个仿真***实时性差、资源紧张、同步性差等固有的不足之处,因此,本***具有资源利用率高,实时性好,仿真***模块化程度高,同步性强,可维护等优点,同时,***实时性和同步性的增强,使得仿真***能够极大的模拟出真实的空间环境,大大提高了***的仿真效率,另外,这种逼真的仿真场景不但可以实现卫星GNC***的半物理验证,同时也可以实现卫星的演示与教学。
本发明实施例提供一种卫星仿真的方法,如图2所示,该方法包括:
S201、仿真控制子***生成模拟卫星,并获取该环境模拟子***发送的该模拟卫星当前时刻对应的磁场信息和该视景轨道仿真子***发送该模拟卫星当前时刻对应的轨道信息。
其中,该仿真控制子***需要将获取的磁场信息和轨道信息转换为该卫星模拟子***对应的信息,并发送至该卫星模拟子***。
S202、该仿真控制子***将该磁场信息和该轨道信息发送至与该磁场信息和轨道信息对应的该卫星模拟子***。
S203、该卫星模拟子***根据该磁场信息和该轨道信息以及配置的力矩得到控制信息,并将该控制信息发送至该仿真控制子***。
示例地,该卫星模拟子***可以采用MSP430-5438A芯片作为主控芯片。
S204、该仿真控制子***根据该卫星模拟子***发送的控制信息得到该模拟卫星在该下一时刻的姿态信息,并根据该下一时刻的姿态信息调整该模拟卫星的姿态。
具体地,该仿真控制子***根据该控制信息得到该模拟卫星当前时刻的横滚角、俯仰角和偏航角,并根据该模拟卫星当前时刻的横滚角、俯仰角和偏航角得到该模拟卫星下一时刻的横滚角、俯仰角和偏航角,并根据该下一时刻的横滚角、俯仰角和偏航角调整该模拟卫星的姿态。
需要说明的是,步骤S204的具体实现方式可以参考上述实施例中关于仿真控制子***的描述,此处不再赘述了。
进一步地,在该仿真控制子***根据该下一时刻的姿态信息调整该模拟卫星的姿态后,该视景轨道仿真子***展示该模拟卫星调整后的姿态和对应的轨道信息,并将该模拟卫星调整后的姿态对应的轨道信息发送至该环境模拟子***,该环境模拟子***根据该模拟卫星调整后的姿态对应的轨道信息得到该模拟卫星调整后的姿态对应的磁场信息,并将该模拟卫星调整后的姿态对应的磁场信息发送至该视景轨道仿真子***,该视景轨道仿真子***展示该模拟卫星调整后的姿态对应的磁场信息。
具体地,该环境模拟子***根据该轨道信息得到经度和余纬,并根据该经度和余纬通过公式:
得到地球磁位势;
其中,V为该地球磁位势,RE为地球半径,r为该模拟卫星到地心的距离,λ为该经度,为该余纬,和分别为n次m阶高阶系数,为n次m阶缔合勒让德函数;
该环境模拟子***根据该地球磁位势通过公式:
分别得到在地心赤道坐标系中X轴、Y轴和Z轴的地磁场强度;
其中,为X轴的地磁场强度,为Y轴的地磁场强度,为Z轴的地磁场强度;
该环境模拟子***通过公式
将在地心赤道坐标系中的地磁场强度转换为在地心轨道坐标系中的地磁场强度;
其中,为在地心轨道坐标系中的地磁场强度,为在地心赤道坐标系中的地磁场强度,Cz(ω)为在Z轴方向上旋转α角度的旋转矩阵,α为近地点幅角,Cx(i)为在X轴方向上旋转i角度的旋转矩阵,i为轨道倾角,Cz(Ω-β)为在Z轴方向上旋转角度Ω-β的旋转矩阵,Ω为升交点赤经,β为该下一时刻的俯仰角;
该环境模拟子***通过公式
将在地心轨道坐标系中的地磁场强度转换为在轨道系坐标中的地磁场强度;
其中,为在轨道系坐标中的地磁场强度,为在地心轨道坐标系中的地磁场强度,Cx(π/2)为在X轴方向上旋转π/2角度的旋转矩阵,CZ(f+π/2)为在Z轴方向上旋转(f+π/2)角度的旋转矩阵,f为真近点角。
需要说明的是,环境模拟子***在根据下一时刻的姿态信息调整该模拟卫星的姿态后,该视景轨道仿真子***展示该模拟卫星调整后的姿态和对应的轨道信息,并将该模拟卫星调整后的姿态对应的轨道信息发送至该环境模拟子***,该环境模拟子***根据该模拟卫星调整后的姿态对应的轨道信息得到该模拟卫星调整后的姿态对应的磁场信息,并将该模拟卫星调整后的姿态对应的磁场信息发送至该视景轨道仿真子***,该视景轨道仿真子***展示该模拟卫星调整后的姿态对应的磁场信息,同时,该环境模拟子***还将该模拟卫星调整后的姿态对应的磁场信息发送至仿真控制子***,该视景轨道仿真子***将该模拟卫星调整后的姿态对应的轨道信息发送至仿真控制子***,以便该仿真控制子***根据该模拟卫星调整后的姿态对应的磁场信息和该模拟卫星调整后的姿态对应的轨道信息继续通过卫星模拟子***获取下一时刻的姿态信息。
另外,在该卫星仿真***启动后,通过该视景轨道仿真子***配置初始轨道信息和初始磁场信息,该视景轨道仿真子***将该初始轨道信息和初始磁场信息发送至仿真控制子***,该仿真控制子***将该初始轨道信息和初始磁场信息发送至对应的卫星模拟子***,卫星模拟子***根据该初始轨道信息和初始磁场信息以及配置的力矩得到控制信息,并将该控制信息发送至仿真控制子***,以便该仿真控制子***根据该控制信息得到下一时刻的姿态信息。
该方法应用于卫星仿真***,通过采用分布式仿真,将整个仿真***的功能分配到各个子模块,能够很好的进行***负载的分配,通过数据链对数据的分发与维护实现仿真***有条不紊的运行,该仿真***涉及到的多个模块,根据自己的任务与载荷的要求,合理分配***资源,使得***能够以最快的速度完成该***的任务,弥补单个仿真***实时性差、资源紧张、同步性差等固有的不足之处,因此,本***具有资源利用率高,实时性好,仿真***模块化程度高,同步性强,可维护等优点,同时,***实时性和同步性的增强,使得仿真***能够极大的模拟出真实的空间环境,大大提高了***的仿真效率,另外,这种逼真的仿真场景不但可以实现卫星GNC***的半物理验证,同时也可以实现卫星的演示与教学。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (4)
1.一种卫星仿真***,其特征在于,包括:仿真控制子***、环境模拟子***、视景轨道仿真子***和卫星模拟子***,其中,所述仿真控制子***分别与所述环境模拟子***、所述视景轨道仿真子***和所述卫星模拟子***连接,所述环境模拟子***与所述视景轨道仿真子***连接,
所述仿真控制子***,用于生成模拟卫星,并获取当前时刻所述模拟卫星对应的磁场信息和轨道信息,并将所述磁场信息和所述轨道信息发送至与所述磁场信息和轨道信息对应的所述卫星模拟子***;
所述卫星模拟子***,用于接收所述仿真控制子***发送的所述磁场信息和所述轨道信息,并根据所述磁场信息和所述轨道信息以及配置的力矩得到控制信息,并将所述控制信息发送至所述仿真控制子***;
所述仿真控制子***,还用于根据所述卫星模拟子***发送的控制信息得到所述模拟卫星在下一时刻的姿态信息,并根据所述下一时刻的姿态信息调整所述模拟卫星的姿态;
所述环境模拟子***,用于获取所述模拟卫星的磁场信息,并将所述磁场信息发送至所述仿真控制子***;
所述视景轨道仿真子***,用于获取所述模拟卫星的轨道信息,并将所述轨道信息发送至所述仿真控制子***;
所述视景轨道仿真子***还用于,展示所述模拟卫星的姿态和对应的轨道信息,并将所述轨道信息发送至所述环境模拟子***;
所述环境模拟子***具体用于,根据所述轨道信息得到磁场信息;
所述环境模拟子***还用于,将所述磁场信息发送至所述视景轨道仿真子***;
所述视景轨道仿真子***还用于,展示所述环境模拟子***发送的磁场信息;
所述仿真控制子***具体用于,根据所述控制信息得到所述模拟卫星当前时刻的横滚角、俯仰角和偏航角,并根据所述模拟卫星当前时刻的横滚角、俯仰角和偏航角得到所述模拟卫星下一时刻的横滚角、俯仰角和偏航角,并根据所述下一时刻的横滚角、俯仰角和偏航角调整所述模拟卫星的姿态;
所述仿真控制子***具体用于,根据所述当前时刻的横滚角、所述俯仰角和所述偏航角通过公式:
得到所述模拟卫星当前时刻的姿态四元数;
其中,Q为所述模拟卫星当前时刻的姿态四元数,
为所述模拟卫星当前时刻的横滚角,θ为所述模拟卫星当前时刻的俯仰角,ψ为所述模拟卫星当前时刻的偏航角;
所述仿真控制子***具体用于,根据所述姿态四元数通过公式:
得到所述模拟卫星从轨道坐标系到本体坐标系的姿态旋转矩阵;其中,为所述姿态旋转矩阵;
所述仿真控制子***具体用于,根据所述控制信息得到配置的力矩和所述模拟卫星的转动惯量以及所述模拟卫星在当前时刻所述本体坐标系相对于惯性坐标系的角速度,并根据所述力矩和所述转动惯量以及所述旋转矩阵通过公式:
得到下一时刻所述本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式,
其中,为下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式,I为所述模拟卫星的转动惯量,为当前时刻所述本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影,Tb为所述力矩在本体坐标系中的投影;
所述仿真控制子***具体用于,对所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影的微分表达式积分得到所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度,并根据所述控制信息得到轨道坐标系相对于惯性坐标系的旋转矩阵,并根据所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度和所述轨道坐标系相对于惯性坐标系的旋转矩阵通过公式:
得到下一时刻所述模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影;
其中,为下一时刻所述模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影,为所述姿态旋转矩阵,为下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中投影,为所述模拟卫星在轨道坐标系相对于惯性坐标系的角速度在所述轨道坐标系中的投影,ωo为所述模拟卫星的轨道角速度;
所述仿真控制子***具体用于,根据所述下一时刻所述模拟卫星在本体坐标系相对于惯性坐标系的角速度在本体坐标系中的投影和所述模拟卫星的转动惯量通过方程:
得到q0的微分表达式和的微分表达式,并对所述q0的微分表达式和的微分表达式分别积分得到所述模拟卫星下一时刻的姿态四元数,并根据所述模拟卫星下一时刻的姿态四元数得到所述模拟卫星下一时刻的横滚角、俯仰角和偏航角;
其中,为q0的微分表达式,为的微分表达式,为下一时刻所述模拟卫星在本体坐标系相对于轨道坐标系的角速度在本体坐标系中的投影,I为所述模拟卫星的转动惯量。
2.根据权利要求1所述的***,其特征在于,所述环境模拟子***具体用于,根据所述轨道信息得到经度和余纬,并根据所述经度和余纬通过公式:
得到地球磁位势;
其中,V为所述地球磁位势,RE为地球半径,r为所述模拟卫星到地心的距离,λ为所述经度,为所述余纬,和分别为n次m阶高阶系数,为n次m阶缔合勒让德函数;
所述环境模拟子***具体用于,根据所述地球磁位势通过公式:
分别得到在地心赤道坐标系中X轴、Y轴和Z轴的地磁场强度;
其中,为X轴的地磁场强度,为Y轴的地磁场强度,为Z轴的地磁场强度;
所述环境模拟子***具体用于,通过公式
将在地心赤道坐标系中的地磁场强度转换为在地心轨道坐标系中的地磁场强度;
其中,为在地心轨道坐标系中的地磁场强度,为在地心赤道坐标系中的地磁场强度,Cz(ω)为在Z轴方向上旋转α角度的旋转矩阵,α为近地点幅角,Cx(i)为在X轴方向上旋转i角度的旋转矩阵,i为轨道倾角,Cz(Ω-β)为在Z轴方向上旋转角度Ω-β的旋转矩阵,Ω为升交点赤经,β为所述下一时刻的俯仰角;
所述环境模拟子***具体用于,通过公式
将在地心轨道坐标系中的地磁场强度转换为在轨道系坐标中的地磁场强度;
其中,为在轨道系坐标中的地磁场强度,为在地心轨道坐标系中的地磁场强度,Cx(π/2)为在X轴方向上旋转π/2角度的旋转矩阵,CZ(f+π/2)为在Z轴方向上旋转(f+π/2)角度的旋转矩阵,f为真近点角。
3.根据权利要求1或2所述的***,其特征在于,所述仿真控制子***与所述卫星模拟子***通过RS232总线连接,所述仿真控制子***与所述环境模拟子***和所述视景轨道仿真子***分别通过以太网连接,所述环境模拟子***和所述视景轨道仿真子***通过以太网连接。
4.一种卫星仿真的方法,其特征在于,包括:
仿真控制子***生成模拟卫星,并获取环境模拟子***发送的所述模拟卫星当前时刻对应的磁场信息和视景轨道仿真子***发送所述模拟卫星当前时刻对应的轨道信息;
所述仿真控制子***将所述磁场信息和所述轨道信息发送至与所述磁场信息和轨道信息对应的卫星模拟子***;
所述卫星模拟子***根据所述磁场信息和所述轨道信息以及配置的力矩得到控制信息,并将所述控制信息发送至所述仿真控制子***;
所述仿真控制子***根据所述卫星模拟子***发送的控制信息得到所述模拟卫星在下一时刻的姿态信息,并根据所述下一时刻的姿态信息调整所述模拟卫星的姿态;
所述仿真控制子***根据所述卫星模拟子***发送的控制信息得到所述模拟卫星在所述下一时刻的姿态信息,并根据所述下一时刻的姿态信息调整所述模拟卫星的姿态包括:
所述仿真控制子***根据所述控制信息得到所述模拟卫星当前时刻的横滚角、俯仰角和偏航角,并根据所述模拟卫星当前时刻的横滚角、俯仰角和偏航角得到所述模拟卫星下一时刻的横滚角、俯仰角和偏航角,并根据所述下一时刻的横滚角、俯仰角和偏航角调整所述模拟卫星的姿态;
在所述根据所述下一时刻的姿态信息调整所述模拟卫星的姿态后,所述方法还包括:
所述视景轨道仿真子***展示所述模拟卫星调整后的姿态和对应的轨道信息,并将所述模拟卫星调整后的姿态对应的轨道信息发送至所述环境模拟子***;
所述环境模拟子***根据所述模拟卫星调整后的姿态对应的轨道信息得到所述模拟卫星调整后的姿态对应的磁场信息,并将所述模拟卫星调整后的姿态对应的磁场信息发送至所述视景轨道仿真子***;
所述视景轨道仿真子***展示所述模拟卫星调整后的姿态对应的磁场信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410497678.5A CN104536302B (zh) | 2014-09-25 | 2014-09-25 | 一种卫星仿真***和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410497678.5A CN104536302B (zh) | 2014-09-25 | 2014-09-25 | 一种卫星仿真***和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104536302A CN104536302A (zh) | 2015-04-22 |
CN104536302B true CN104536302B (zh) | 2017-06-06 |
Family
ID=52851846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410497678.5A Expired - Fee Related CN104536302B (zh) | 2014-09-25 | 2014-09-25 | 一种卫星仿真***和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104536302B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105487405B (zh) * | 2015-12-17 | 2019-01-29 | 西安测绘研究所 | 低低跟踪重力测量卫星半物理仿真*** |
CN105867163A (zh) * | 2016-03-29 | 2016-08-17 | 上海卫星工程研究所 | 基于Labwindows和STK的航天任务分析方法 |
CN109343081A (zh) * | 2018-10-10 | 2019-02-15 | 中国人民解放军国防科技大学 | 一种gps信号动态接收环境仿真方法及*** |
CN110321644B (zh) * | 2019-07-09 | 2021-03-26 | 中国科学院电子学研究所 | 一种面向卫星***的分布式仿真***及其仿真方法 |
CN111460184B (zh) * | 2020-03-05 | 2023-04-11 | 哈尔滨工业大学 | 遥感卫星载荷数据动态模拟装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320524A (zh) * | 2008-04-22 | 2008-12-10 | 北京航空航天大学 | 多处理器实时仿真平台 |
CN102354123A (zh) * | 2011-07-18 | 2012-02-15 | 北京航空航天大学 | 一种跨平台可扩展的卫星动态仿真测试*** |
CN102411312A (zh) * | 2011-11-16 | 2012-04-11 | 清华大学 | 敏感器模拟器及具有该敏感器模拟器的卫星闭环仿真*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2674870B2 (ja) * | 1990-10-04 | 1997-11-12 | 日本電気株式会社 | 人工衛星追跡管制設備シミュレータ |
-
2014
- 2014-09-25 CN CN201410497678.5A patent/CN104536302B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320524A (zh) * | 2008-04-22 | 2008-12-10 | 北京航空航天大学 | 多处理器实时仿真平台 |
CN102354123A (zh) * | 2011-07-18 | 2012-02-15 | 北京航空航天大学 | 一种跨平台可扩展的卫星动态仿真测试*** |
CN102411312A (zh) * | 2011-11-16 | 2012-04-11 | 清华大学 | 敏感器模拟器及具有该敏感器模拟器的卫星闭环仿真*** |
Non-Patent Citations (2)
Title |
---|
卫星姿轨控实时仿真***及串口通信建模;朱承元等;《哈尔滨工业大学学报》;20050131;第37卷(第1期);全文 * |
基于RTX的卫星姿轨控***地面实时仿真***;夏红伟等;《计算机仿真》;20060930;第23卷(第9期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104536302A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104536302B (zh) | 一种卫星仿真***和方法 | |
CN102354123B (zh) | 一种跨平台可扩展的卫星动态仿真测试*** | |
CN106530894B (zh) | 一种飞行训练器虚拟平显方法及*** | |
CN107544467B (zh) | 整星条件下的双星编队控制闭环测试***及方法 | |
You | Space microsystems and micro/nano satellites | |
CN103149846A (zh) | 飞行机器人控制***半物理仿真平台 | |
GB2548109A (en) | Apparatus and method for satellite payload development | |
CN112650076B (zh) | 一种星群协同控制地面仿真*** | |
CN103646127A (zh) | 卫星轨道姿态可视化三维显示方法 | |
CN106227066B (zh) | 航天器空间操控可配置视景显示方法 | |
CN102841966A (zh) | Vpp-STK卫星仿真开发与运行平台*** | |
CN105466477A (zh) | 一种面向卫星目标和恒星目标的天基观测模拟***及方法 | |
Reid et al. | The Artemis Jr. rover: Mobility platform for lunar ISRU mission simulation | |
RU102823U1 (ru) | Комплексное обучающее устройство | |
CN105487405B (zh) | 低低跟踪重力测量卫星半物理仿真*** | |
CN114529602B (zh) | 一种空间多目标态势监视方法及装置 | |
CN116227112A (zh) | 面向遥感卫星星座体系协同的半物理仿真***与仿真方法 | |
CN113359510A (zh) | 北斗卫星导航***信号模拟器数据实时仿真***及方法 | |
CN103941592A (zh) | 飞行机器人动力学模型的在线建模方法 | |
Bi et al. | The state of the art of information integration in space applications | |
Diston | Computational modelling and simulation of aircraft and the environment, volume 1: Platform kinematics and synthetic environment | |
CN202453676U (zh) | 飞行机器人控制***半物理仿真平台 | |
Basmadji et al. | Development of ground station for a terrain observer-hardware in the loop simulations | |
CN115113639B (zh) | 一种无人机飞行控制与模拟训练方法及装置 | |
CN105487402B (zh) | 一种星敏与陀螺组合定姿全物理仿真试验方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170606 Termination date: 20210925 |
|
CF01 | Termination of patent right due to non-payment of annual fee |