CN111679592A - A closed-loop semi-physical simulation system and method for spacecraft pursuit and escape game - Google Patents
A closed-loop semi-physical simulation system and method for spacecraft pursuit and escape game Download PDFInfo
- Publication number
- CN111679592A CN111679592A CN202010571039.4A CN202010571039A CN111679592A CN 111679592 A CN111679592 A CN 111679592A CN 202010571039 A CN202010571039 A CN 202010571039A CN 111679592 A CN111679592 A CN 111679592A
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- simulation
- escape
- relative
- model
- 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
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 169
- 238000000034 method Methods 0.000 title claims description 41
- 238000005259 measurement Methods 0.000 claims abstract description 61
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims description 46
- 238000004422 calculation algorithm Methods 0.000 claims description 34
- 230000008569 process Effects 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000012544 monitoring process Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000003287 optical effect Effects 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 6
- 239000003550 marker Substances 0.000 claims description 6
- 238000000691 measurement method Methods 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 claims description 2
- 239000002131 composite material Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 10
- 230000003993 interaction Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000002474 experimental method Methods 0.000 description 10
- 238000013519 translation Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Processing Or Creating Images (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开一种航天器追逃博弈闭环半实物仿真系统,包括相对测量模拟子系统、星载GNC模拟子系统、闭环运行控制子系统和相对运动学模拟子系统。根据并行仿真原理将相对测量模拟子系统、星载GNC模拟子系统和相对运动学模拟子系统形成一个并行的闭环运行整体系统,使仿真中的时间同步关系更加协调,减少时间的延迟,且可通过仿真结果反映出导航、制导和计算环节间的相互影响;该仿真系统引入了真实的测量环节、制导环节以及星载计算环节,相对纯数学仿真,置信度更高;相对物理仿真,成本低,实施简单;此外,该仿真系统对航天器轨道自然运动和追逃博弈机动采用独立并行化处理,使得仿真更加逼真。
The invention discloses a closed-loop semi-physical simulation system of a spacecraft chase and escape game, which includes a relative measurement simulation subsystem, an on-board GNC simulation subsystem, a closed-loop operation control subsystem and a relative kinematics simulation subsystem. According to the principle of parallel simulation, the relative measurement simulation subsystem, the spaceborne GNC simulation subsystem and the relative kinematics simulation subsystem form a parallel closed-loop operation overall system, which makes the time synchronization relationship in the simulation more coordinated, reduces the time delay, and can The interaction among navigation, guidance and calculation links is reflected through the simulation results; the simulation system introduces real measurement links, guidance links and onboard calculation links, which has higher confidence than pure mathematical simulation; lower cost than physical simulation , the implementation is simple; in addition, the simulation system adopts independent parallel processing for the natural motion of the spacecraft orbit and the game maneuver of pursuit and escape, which makes the simulation more realistic.
Description
技术领域technical field
本发明涉及航天器轨道控制技术领域,尤其是一种航天器追逃博弈闭环半实物仿真系统及方法。The invention relates to the technical field of spacecraft orbit control, in particular to a closed-loop semi-physical simulation system and method of a spacecraft pursuit and escape game.
背景技术Background technique
随着世界主要大国的交会对接技术的发展成熟,我国在轨航天器当前面临着被接近、捕获甚至操纵的威胁。特别当两航天器都能施加自主追踪机动或规避机动时,就会形成空间追逃动态博弈的局面。在此背景下,航天器追逃博弈问题也成为学术界研究的热点问题,一系列支撑航天器开展轨道追逃博弈的导航、制导和控制算法也相应地被发展出来,包含非合作目标仅侧角相对导航算法、改进的视线制导算法以及微分对策制导算法等。但是当前仍缺乏一套检验航天器追逃博弈导航和制导算法有效性的实验系统,特别是能对航天器追逃博弈整个闭环过程进行仿真的实验系统。由于航天器轨道追逃问题是一个综合通信、导航、控制、计算和运动的复杂过程,在真实的太空环境进行追逃博弈物理仿真的成本较高、风险较大、周期较长,所以迫切需要一种地面半实物仿真系统,对航天器追逃博弈中关键的导航环节、制导环节、星载计算环节以及所形成的闭环过程进行仿真,对轨道追逃博弈相关核心算法的效果进行检验。但目前公开文献中尚未发现国内有可用于航天器追逃博弈的闭环半实物仿真系统。With the development and maturity of the rendezvous and docking technology of the world's major powers, my country's orbiting spacecraft is currently facing the threat of being approached, captured or even manipulated. Especially when both spacecraft can perform autonomous tracking maneuvers or evasive maneuvers, a dynamic game of space pursuit and escape will be formed. In this context, the problem of the spacecraft pursuit and escape game has also become a hot issue in academic research. A series of navigation, guidance and control algorithms that support the spacecraft to carry out the orbital pursuit and escape game have also been developed accordingly, including non-cooperative targets only side by side. Angle-relative navigation algorithm, improved line-of-sight guidance algorithm and differential game guidance algorithm, etc. However, there is still a lack of an experimental system to test the effectiveness of the navigation and guidance algorithms of the spacecraft pursuit and escape game, especially an experimental system that can simulate the entire closed-loop process of the spacecraft pursuit and escape game. Since the problem of spacecraft orbital pursuit and escape is a complex process of integrated communication, navigation, control, calculation and motion, the physical simulation of the pursuit and escape game in the real space environment is costly, risky, and has a long period of time. Therefore, it is urgent to A semi-physical simulation system on the ground simulates the key navigation links, guidance links, on-board calculation links and the closed-loop process formed in the spacecraft pursuit and escape game, and tests the effect of the core algorithms related to the orbital pursuit and escape game. However, no closed-loop hardware-in-the-loop simulation system that can be used for spacecraft pursuit and escape game has been found in the open literature.
发明内容SUMMARY OF THE INVENTION
本发明提供一种航天器追逃博弈闭环半实物仿真系统及方法,用于克服现有技术中没有能够实现对航天器追逃博弈中关键进行仿真、对核心算法进行检验的方法等缺陷。The present invention provides a closed-loop semi-physical simulation system and method for a spacecraft pursuit and escape game, which is used to overcome the defects of the prior art, such as the lack of a method capable of simulating the key points in the spacecraft pursuit and escape game and checking the core algorithm.
为实现上述目的,本发明提出一种航天器追逃博弈闭环半实物仿真系统,包括:In order to achieve the above-mentioned purpose, the present invention proposes a closed-loop semi-physical simulation system of a spacecraft chase and escape game, including:
相对测量模拟子系统,用于模拟太空中的测量环境和测量手段对逃逸航天器模型和追踪航天器模型进行信息采集,获得逃逸航天器模型和追踪航天器模型的相对状态,将所述相对状态传递给星载GNC模拟子系统和闭环运行控制子系统;所述相对状态包括图像、逃逸航天器模型与追踪航天器模型之间的距离以及逃逸航天器模型和追踪航天器模型的倾角;The relative measurement simulation subsystem is used to simulate the measurement environment and measurement methods in space, collect information on the escape spacecraft model and the tracking spacecraft model, obtain the relative state of the escape spacecraft model and the tracking spacecraft model, and convert the relative state Passed to the onboard GNC simulation subsystem and the closed-loop operation control subsystem; the relative state includes the image, the distance between the escape spacecraft model and the tracking spacecraft model, and the inclination of the escape spacecraft model and the tracking spacecraft model;
星载GNC模拟子系统,用于根据所述相对状态,利用导航滤波算法进行滤波迭代,获得包含相对位置和相对速度在内的六维度状态信息,根据所述六维度状态信息,利用追逃博弈机动算法,获得追踪航天器模型和逃逸航天器模型的控制量,将所述控制量转化为第一控制指令发送给闭环运行控制子系统;The satellite-borne GNC simulation subsystem is used to perform filtering iteration by using a navigation filtering algorithm according to the relative state, to obtain six-dimensional state information including relative position and relative velocity, and use the pursuit and escape game according to the six-dimensional state information. The maneuvering algorithm obtains the control quantities of the tracking spacecraft model and the escape spacecraft model, and converts the control quantities into first control commands and sends them to the closed-loop operation control subsystem;
闭环运行控制子系统,用于根据所述相对状态对航天器轨道自然运动进行积分生成第二控制指令,对所述第一控制指令和所述第二控制指令进行加和获得综合指令,将所述综合指令解算获得逃逸航天器模型和追踪航天器模型的位置速度、姿态角和姿态角速度,将所述位置速度、姿态角和姿态角速度发送给相对运动学模拟子系统;The closed-loop operation control subsystem is used to integrate the natural motion of the spacecraft orbit according to the relative state to generate a second control command, add the first control command and the second control command to obtain a comprehensive command, and combine the Obtaining the position velocity, attitude angle and attitude angular velocity of the escape spacecraft model and the tracking spacecraft model by calculating the integrated command, and sending the position velocity, attitude angle and attitude angular velocity to the relative kinematics simulation subsystem;
相对运动学模拟子系统,用于根据所述位置速度、姿态角和姿态角速度驱动逃逸航天器模型和追踪航天器模型运动;a relative kinematics simulation subsystem for driving the escape spacecraft model and tracking the motion of the spacecraft model according to the position velocity, attitude angle and attitude angular velocity;
所述闭环运行控制子系统分别与所述相对测量模拟子系统、所述星载GNC模拟子系统和所述相对运动学模拟子系统双向电连接,将所述相对测量模拟子系统、所述星载GNC模拟子系统和所述相对运动学模拟子系统形成一个并行的闭环运行整体系统。The closed-loop operation control subsystem is electrically connected in both directions with the relative measurement simulation subsystem, the on-board GNC simulation subsystem and the relative kinematics simulation subsystem, and the relative measurement simulation subsystem, the satellite The GNC simulation subsystem and the relative kinematics simulation subsystem form a parallel closed-loop running overall system.
为实现上述目的,本发明还提出一种航天器追逃博弈闭环半实物仿真方法,利用如上述所述的航天器追逃博弈闭环半实物仿真系统进行仿真。In order to achieve the above object, the present invention also proposes a closed-loop semi-physical simulation method for the spacecraft pursuit and escape game, which uses the above-mentioned closed-loop semi-physical simulation system for the spacecraft pursuit and escape game for simulation.
与现有技术相比,本发明的有益效果有:Compared with the prior art, the beneficial effects of the present invention are:
1、本发明提供的航天器追逃博弈闭环半实物仿真系统,根据并行仿真原理将相对测量模拟子系统、星载GNC模拟子系统和相对运动学模拟子系统形成一个并行的闭环运行整体系统,使仿真中的时间同步关系更加协调,减少时间的延迟,且可通过仿真结果反映出导航、制导和计算环节间的相互影响。1. The closed-loop semi-physical simulation system of the spacecraft chase and escape game provided by the present invention forms a parallel closed-loop operation overall system according to the principle of parallel simulation, the relative measurement simulation subsystem, the on-board GNC simulation subsystem and the relative kinematics simulation subsystem, The time synchronization relationship in the simulation is more coordinated, the time delay is reduced, and the interaction between the navigation, guidance and calculation links can be reflected through the simulation results.
2、本发明提供的航天器追逃博弈闭环半实物仿真系统引入了真实的测量环节、制导环节以及星载计算环节,相对纯数学仿真,置信度更高;相对物理仿真,成本低,实施简单。2. The closed-loop semi-physical simulation system of the spacecraft chase and escape game provided by the present invention introduces real measurement links, guidance links and on-board calculation links, which is relatively pure mathematical simulation and has a higher degree of confidence; relatively physical simulation, low cost and simple implementation .
3、本发明提供的航天器追逃博弈闭环半实物仿真系统,对航天器轨道自然运动和追逃博弈机动采用独立并行化处理,通过这样的并行设计可以保证航天器模型在没有机动指令或机动规划存在时间延迟的情况下,航天器模型在仿真中可保持自由漂移的状态而不是停止运动,与实际太空中一致,使得仿真更加逼真。3. The closed-loop semi-physical simulation system of the spacecraft pursuit and escape game provided by the present invention adopts independent parallel processing for the natural motion of the spacecraft orbit and the maneuvering of the pursuit and escape game. Through such a parallel design, it can be ensured that the spacecraft model has no maneuvering instructions or maneuvers. When there is a time delay in planning, the spacecraft model can maintain a free drift state instead of stopping motion in the simulation, which is consistent with the actual space, making the simulation more realistic.
4、本发明提供的航天器追逃博弈闭环半实物仿真系统可很好的用于进行航天器导航和制导算法的验证进而开展航天器追逃博弈闭环半实物仿真。4. The closed-loop semi-physical simulation system of the spacecraft pursuit and escape game provided by the present invention can be well used to verify the navigation and guidance algorithms of the spacecraft and then carry out the closed-loop semi-physical simulation of the spacecraft pursuit and escape game.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。In order to explain the embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention, and for those of ordinary skill in the art, other drawings can also be obtained according to the structures shown in these drawings without creative efforts.
图1为本发明提供的航天器追逃博弈闭环半实物仿真系统的总体架构图;Fig. 1 is the overall architecture diagram of the closed-loop semi-physical simulation system of the spacecraft chase and escape game provided by the present invention;
图2为本发明提供的航天器追逃博弈闭环半实物仿真系统的工作原理图;Fig. 2 is the working principle diagram of the closed-loop semi-physical simulation system of the spacecraft chase and escape game provided by the present invention;
图3为本发明提供的航天器追逃博弈闭环半实物仿真系统的闭环运行示意图;Fig. 3 is the closed-loop operation schematic diagram of the closed-loop semi-physical simulation system of the spacecraft chase and escape game provided by the present invention;
图4为九自由度相对运动平台的结构图;Figure 4 is a structural diagram of a nine-degree-of-freedom relative motion platform;
图5为本发明采用的逃逸航天器模型和追踪航天器模型的结构图;Fig. 5 is the structure diagram of the escape spacecraft model and the tracking spacecraft model adopted by the present invention;
图6为本发明中逃逸航天器模型和追踪航天器模型上标志点设置图;6 is a diagram showing the setting of marker points on the escaping spacecraft model and the tracking spacecraft model in the present invention;
图7为本发明提供的航天器追逃博弈闭环半实物仿真系统中测量环节、运动环节和规划计算环节时间同步关系图;Fig. 7 is the time synchronization relation diagram of the measurement link, the movement link and the planning calculation link in the closed-loop semi-physical simulation system of the spacecraft chase and escape game provided by the present invention;
图8a为实施例1中R-bar相对距离滤波实验结果图;Fig. 8a is a graph of experimental results of R-bar relative distance filtering in Example 1;
图8b为实施例1中R-bar相对速度滤波实验结果图;Figure 8b is a graph of the experimental results of R-bar relative velocity filtering in Example 1;
图8c为实施例1中V-bar相对距离滤波实验结果图;Figure 8c is a graph of the experimental results of V-bar relative distance filtering in Example 1;
图8d为实施例1中V-bar相对速度滤波实验结果图;Figure 8d is a graph of the experimental results of V-bar relative velocity filtering in Example 1;
图8e为实施例1中H-bar相对距离滤波实验结果图;Figure 8e is a graph of experimental results of H-bar relative distance filtering in Example 1;
图8f为实施例1中H-bar相对距离滤波实验结果图;Figure 8f is a graph of experimental results of H-bar relative distance filtering in Example 1;
图9a为实施例2中R-bar相对距离的博弈实验结果图;Fig. 9a is the game experiment result diagram of R-bar relative distance in
图9b为实施例2中R-bar相对速度的博弈实验结果图;Fig. 9b is the game experiment result diagram of R-bar relative velocity in
图9c为实施例2中V-bar相对距离的博弈实验结果图;Fig. 9c is the game experiment result diagram of V-bar relative distance in
图9d为实施例2中V-bar相对速度的博弈实验结果图;Fig. 9d is the game experiment result diagram of V-bar relative velocity in
图9e为实施例2中H-bar相对距离的博弈实验结果图;Fig. 9e is the game experiment result diagram of H-bar relative distance in
图9f为实施例2中H-bar相对速度的博弈实验结果图。FIG. 9f is a graph showing the result of a game experiment of the relative velocity of H-bar in Example 2. FIG.
附图符号说明:1:X方向平动电机;2:Y方向平动电机;3:逃逸模拟端绕Y轴转动电机;4:追踪模拟端绕Y轴转动电机;5:标志点;51:远场标志点;52:近场标志点。Description of the drawing symbols: 1: X-direction translation motor; 2: Y-direction translation motor; 3: The escape analog end rotates the motor around the Y axis; 4: The tracking analog end rotates the motor around the Y axis; 5: Mark point; 51: Far-field marker point; 52: Near-field marker point.
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。The realization, functional characteristics and advantages of the present invention will be further described with reference to the accompanying drawings in conjunction with the embodiments.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
另外,本发明各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。In addition, the technical solutions between the various embodiments of the present invention can be combined with each other, but must be based on the realization by those of ordinary skill in the art. When the combination of technical solutions is contradictory or cannot be realized, it should be considered that the combination of technical solutions does not exist and is not within the scope of protection claimed by the present invention.
本发明提出一种航天器追逃博弈闭环半实物仿真系统,利用实验室的九自由度相对运动平台实现对逃逸航天器模型和追踪航天器模型的驱动。九自由度相对运动平台实如图4所示,一端为追踪模拟端,其上安装有追踪航天器模型,控制着相对运动中平台X轴的平动,最大速度为0.5m/s,最大加速度为0.05m/s2;另一端为逃逸模拟端,其上安装有逃踪航天器模型,控制着相对运动中平台YOZ平面的平动,最大速度为0.4m/s,最大加速度为0.05m/s2;此外,九自由度相对运动平台还包括导轨、电机、通信电缆、底座等。该九自由度相对运动平台可提供3个自由度的位置平动和追踪航天器模型、逃逸航天器模型各3个自由度的姿态转动,允许的最大转动角±30°,最大转动角速度为5°/s。The invention proposes a closed-loop semi-physical simulation system of a spacecraft pursuit and escape game, which utilizes a nine-degree-of-freedom relative motion platform in a laboratory to drive the escape spacecraft model and the pursuit spacecraft model. The nine-degree-of-freedom relative motion platform is shown in Figure 4. One end is the tracking simulation end, and the tracking spacecraft model is installed on it, which controls the translation of the X-axis of the platform during the relative motion. The maximum speed is 0.5m/s, and the maximum acceleration is 0.05m/s2; the other end is the escape simulation end, on which is installed the escape spacecraft model, which controls the translation of the YOZ plane of the platform in relative motion, with a maximum speed of 0.4m/s and a maximum acceleration of 0.05m/s2 ; In addition, the nine-degree-of-freedom relative motion platform also includes guide rails, motors, communication cables, bases, and the like. The 9-DOF relative motion platform can provide 3-DOF position translation and three-degree-of-freedom attitude rotation for tracking the spacecraft model and escape spacecraft model. The maximum allowable rotation angle is ±30°, and the maximum rotation angular velocity is 5 °/s.
逃逸航天器模型和追踪航天器模型均如图5所示,该逃逸航天器模型和追踪航天器模型上均安装有图像采集设备和测量设备,其中图像采集设备的安装位置分别为逃逸航天器模型和追踪航天器模型的中心轴线、正面的底线中央和正面侧边线中央。Both the escape spacecraft model and the tracking spacecraft model are shown in Figure 5. Image acquisition equipment and measurement equipment are installed on the escape spacecraft model and the tracking spacecraft model. The installation positions of the image acquisition equipment are respectively the escape spacecraft model. And trace the center axis of the spacecraft model, the center of the bottom line of the front and the center of the side line of the front.
由于当追踪航天器模型与逃逸航天器模型的距离较远时,目标在图像采集设备中的成像是一个光斑或者亮点;而当距离较近时,目标在图像采集设备中呈现出大体轮廓。因此,本发明在逃逸航天器模型上安装有两套光学标志点,分别适用于近场和远场的情况,如图6所示。Because when the distance between the tracking spacecraft model and the escape spacecraft model is far, the image of the target in the image acquisition device is a light spot or bright spot; when the distance is short, the target appears in the image acquisition device as a general outline. Therefore, in the present invention, two sets of optical marker points are installed on the escape spacecraft model, which are respectively applicable to the near field and the far field, as shown in FIG. 6 .
本发明提出一种航天器追逃博弈闭环半实物仿真系统,如图1~3所示,包括:The present invention proposes a closed-loop semi-physical simulation system for spacecraft pursuit and escape game, as shown in Figures 1 to 3, including:
(1)相对测量模拟子系统,用于模拟太空中的测量环境和测量手段对逃逸航天器模型和追踪航天器模型进行信息采集,获得逃逸航天器模型和追踪航天器模型的相对状态,将所述相对状态传递给星载GNC模拟子系统和闭环运行控制子系统;所述相对状态包括图像、逃逸航天器模型与追踪航天器模型之间的距离以及逃逸航天器模型和追踪航天器模型的倾角;(1) The relative measurement simulation subsystem is used to simulate the measurement environment and measurement methods in space to collect information on the escape spacecraft model and the tracking spacecraft model, and obtain the relative state of the escape spacecraft model and the tracking spacecraft model. The relative state is passed to the onboard GNC simulation subsystem and the closed-loop operation control subsystem; the relative state includes the image, the distance between the escape spacecraft model and the tracking spacecraft model, and the inclination of the escape spacecraft model and the tracking spacecraft model. ;
图像由设置在航天器模型上的图像采集设备采集得到;The image is acquired by the image acquisition device set on the spacecraft model;
逃逸航天器模型与追踪航天器模型之间的距离由设置在航天器模型上的激光测距仪采集得到;The distance between the escape spacecraft model and the tracking spacecraft model is acquired by the laser rangefinder set on the spacecraft model;
逃逸航天器模型和追踪航天器模型的倾角由设置在航天器模型上的测角仪采集得到。The inclination angles of the escape spacecraft model and the tracking spacecraft model are acquired by the goniometer set on the spacecraft model.
(2)星载GNC(guidance navigation and control,导航、制导与控制)模拟子系统,用于根据相对状态,利用导航滤波算法进行滤波迭代,获得包含相对位置和相对速度在内的六维度状态信息,根据六维度状态信息,利用追逃博弈机动算法,获得追踪航天器模型和逃逸航天器模型的控制量,将控制量转化为第一控制指令发送给闭环运行控制子系统;(2) On-board GNC (guidance navigation and control, navigation, guidance and control) simulation subsystem, which is used to perform filtering iteration by using the navigation filtering algorithm according to the relative state to obtain six-dimensional state information including relative position and relative velocity , according to the six-dimensional state information, use the pursuit-escape game maneuvering algorithm to obtain the control quantities of the tracking spacecraft model and the escape spacecraft model, and convert the control quantities into first control commands and send them to the closed-loop operation control subsystem;
控制量为追踪航天器模型和逃逸航天器模型在位置速度、姿态角和姿态角速度上的变化量。The control quantity is the variation of the position velocity, attitude angle and attitude angular velocity of the tracking spacecraft model and the escape spacecraft model.
星载GNC模拟子系统的功能在模拟的星载机上完成,星载GNC模拟子系统内的计算在星载计算机芯片上进行。The functions of the on-board GNC simulation subsystem are completed on the simulated on-board aircraft, and the calculations in the on-board GNC simulation subsystem are performed on the on-board computer chip.
(3)闭环运行控制子系统,用于根据相对状态对航天器轨道自然运动进行积分生成第二控制指令,对第一控制指令和述第二控制指令进行加和获得综合指令,将综合指令解算获得逃逸航天器模型和追踪航天器模型的位置速度、姿态角和姿态角速度,将位置速度、姿态角和姿态角速度发送给相对运动学模拟子系统;(3) The closed-loop operation control subsystem is used to integrate the natural motion of the spacecraft orbit according to the relative state to generate a second control command, add the first control command and the second control command to obtain a comprehensive command, and solve the comprehensive command Calculate the position velocity, attitude angle and attitude angular velocity of the escape spacecraft model and the tracking spacecraft model, and send the position velocity, attitude angle and attitude angular velocity to the relative kinematics simulation subsystem;
解算包括六自由度动力学解算、姿态运动学解算和相对运动解算。The solution includes six-degree-of-freedom dynamics solution, attitude kinematics solution and relative motion solution.
(4)相对运动学模拟子系统,用于根据位置速度、姿态角和姿态角速度驱动逃逸航天器模型和追踪航天器模型运动;(4) The relative kinematics simulation subsystem is used to drive the escape spacecraft model and track the motion of the spacecraft model according to the position velocity, attitude angle and attitude angular velocity;
相对运动学模拟子系统主要通过九自由度相对运动平台实现其功能。The relative kinematics simulation subsystem realizes its functions mainly through the nine-degree-of-freedom relative motion platform.
闭环运行控制子系统分别与相对测量模拟子系统、星载GNC模拟子系统和相对运动学模拟子系统双向电连接,将相对测量模拟子系统、星载GNC模拟子系统和相对运动学模拟子系统形成一个并行的闭环运行整体系统。The closed-loop operation control subsystem is electrically connected to the relative measurement and simulation subsystem, the spaceborne GNC simulation subsystem and the relative kinematics simulation subsystem, respectively, and connects the relative measurement and simulation subsystem, the spaceborne GNC simulation subsystem and the relative kinematics simulation subsystem. Form a parallel closed-loop operation of the overall system.
为了使仿真过程中各个子系统的时间同步关系更加协调,减少时间的延迟,在设计闭环运行整体系统时将测量环节、运动环节和规划计算环节相互并行,在tk节点到tk+1节点的一帧中,各个环节的时间同步关系如图7所示。同时,为了使得九自由度的相对运动平台运动的连续性更好,在并行的仿真实验中,设仿真步长为Δt,测量周期为Δtz,运动周期为Δtm,三者应满足以下同步关系:Δt≥Δtm>>Δtz,前一个不等式可取等号,为了保证测量数据接近仿真节点,测量频率相对仿真步长要足够高。In order to make the time synchronization relationship of each subsystem in the simulation process more coordinated and reduce the time delay, the measurement link, motion link and planning calculation link are parallel to each other when designing the closed - loop operation of the overall system. In a frame of , the time synchronization relationship of each link is shown in Figure 7. At the same time, in order to make the motion continuity of the nine-degree-of-freedom relative motion platform better, in the parallel simulation experiment, set the simulation step size as Δt, the measurement period as Δt z , and the motion period as Δt m , the three should satisfy the following synchronization Relation: Δt≥Δt m >>Δt z , the previous inequality can take the equal sign, in order to ensure that the measurement data is close to the simulation node, the measurement frequency should be high enough relative to the simulation step size.
在其中一个实施例中,航天器追逃博弈闭环半实物仿真系统还包括平台综合监控子系统,用于对闭环运行整体系统进行监视。In one of the embodiments, the closed-loop hardware-in-the-loop simulation system for the spacecraft pursuit and escape game further includes a platform integrated monitoring subsystem, which is used to monitor the entire closed-loop operation system.
平台综合监控子系统的功能是对仿真数据监视,对试验现场视频监视以及将九自由度相对运动平台仿真模拟运动反演展示。仿真数据监视包含对航天器模型仿真姿态数据监视、航天器模型仿真位置数据监视、航天器模型仿真相对位置速度监视。在仿真过程的任意地方出现异常或者超限时,平台综合监控子系统将进行及时停止指令的发送,终止仿真进程。The functions of the platform integrated monitoring subsystem are to monitor the simulation data, monitor the video of the test site, and display the simulated motion of the nine-degree-of-freedom relative motion platform. The simulation data monitoring includes the monitoring of the spacecraft model simulation attitude data, the spacecraft model simulation position data monitoring, and the spacecraft model simulation relative position speed monitoring. When an exception or overrun occurs anywhere in the simulation process, the platform integrated monitoring subsystem will send a timely stop command to terminate the simulation process.
在另一个实施例中,相对测量模拟子系统包括光学环境模拟模块和测量模块;In another embodiment, the relative measurement simulation subsystem includes an optical environment simulation module and a measurement module;
光学环境模拟模块用于模拟包括星空背景和地影遮挡的实际太空环境;对于星空背景的模拟需要用星空图片的静态场景模拟,同时要用不同功率的LED光源模拟恒星亮度。对于位置动态变化较大的恒星如太阳,采用可移动光源进行近似模拟。The optical environment simulation module is used to simulate the actual space environment including the starry sky background and ground shadow occlusion; for the simulation of the starry sky background, the static scene simulation of the starry sky image is required, and the brightness of the stars is simulated by LED light sources of different powers. For stars with large dynamic changes in position, such as the sun, a movable light source is used for approximate simulation.
测量模块用于在光学环境模拟模块提供的环境下通过测量设备实时测量逃逸航天器模型与追踪航天器模型之间的距离以及逃逸航天器模型和追踪航天器模型的倾角,通过图像采集设备实时采集逃逸航天器模型与追踪航天器模型的图像。The measurement module is used to measure the distance between the escape spacecraft model and the tracking spacecraft model and the inclination angle of the escape spacecraft model and the tracking spacecraft model in real time through the measurement equipment in the environment provided by the optical environment simulation module, and collect in real time through the image acquisition equipment. Images of the escape and pursuit spacecraft models.
本实施例中,图像采集设备选用型号为A5131MG75的相机,分辨率为1280×1024,帧率为75fps,光学尺寸为1/2″,像素尺寸为4.8μm,镜头型号为:MH0620S。有效焦距为6mm,相对孔径为F2.0,水平视场角为60.8°,垂直视场角为42.7°,畸变率小于0.1%,最小物距为0.1m,接口为C-Mount,滤镜接口M25.5×P0.5。In this embodiment, the image acquisition device selects the camera model A5131MG75, the resolution is 1280×1024, the frame rate is 75fps, the optical size is 1/2″, the pixel size is 4.8μm, and the lens model is: MH0620S. The effective focal length is 6mm, the relative aperture is F2.0, the horizontal field of view is 60.8°, the vertical field of view is 42.7°, the distortion rate is less than 0.1%, the minimum object distance is 0.1m, the interface is C-Mount, and the filter interface is M25.5 ×P0.5.
测量设备包括激光测距仪和倾角仪。激光测距仪型号为SLD-C30,量程为0.15m-30m,测量精度为2mm,输出接口为RS485。测量设备中倾角仪型号为SCA126V,测量角度范围-60°到+60°,分辨率0.01,绝对精度0.08,响应时间0.02S,输出接口RS485。Measurement equipment includes laser rangefinders and inclinometers. The laser rangefinder model is SLD-C30, the range is 0.15m-30m, the measurement accuracy is 2mm, and the output interface is RS485. The model of the inclinometer in the measuring equipment is SCA126V, the measuring angle range is -60° to +60°, the resolution is 0.01, the absolute accuracy is 0.08, the response time is 0.02S, and the output interface is RS485.
在下一实施例中,星载GNC模拟子系统包括图像处理模块、导航滤波模块和机动规划模块;In the next embodiment, the on-board GNC simulation subsystem includes an image processing module, a navigation filtering module and a maneuver planning module;
图像处理模块用于对所述图像进行处理,获得图像采集设备的测角;测角即为图像采集设备采集图像时的角度。The image processing module is used for processing the image to obtain the angle measurement of the image acquisition device; the angle measurement is the angle when the image acquisition device captures the image.
导航滤波模块用于根据测角以及逃逸航天器模型与追踪航天器模型之间的距离,利用导航滤波算法进行滤波迭代,获得包含相对位置和相对速度在内的六维度状态信息;The navigation filtering module is used to perform filtering iteration by using the navigation filtering algorithm according to the angle measurement and the distance between the escape spacecraft model and the tracking spacecraft model to obtain six-dimensional state information including relative position and relative velocity;
机动规划模块用于根据六维度状态信息,利用追逃博弈机动算法,获得追踪航天器模型和逃逸航天器模型的控制量,将控制量转化为第一控制指令发送给闭环运行控制子系统。The maneuver planning module is used to obtain the control quantities of the tracking spacecraft model and the escape spacecraft model according to the six-dimensional state information and using the pursuit-escape game maneuvering algorithm, and convert the control quantities into first control instructions and send them to the closed-loop operation control subsystem.
本实施例中,星载GNC模拟子系统内的计算在星载计算机芯片上进行,本实施例采用TMS320C6000系列DSP芯片进行数据处理,主频300MHZ,系统的外围设备集成在PCB底板上,外围设备包含LCD接口(显示终端)、外部扩展(DSP扩展槽)、程序下载与在线调试接口(JATG)、电源接口、网络通讯接口(Ethnernet)和UART串口。In this embodiment, the calculation in the on-board GNC simulation subsystem is performed on the on-board computer chip. In this embodiment, TMS320C6000 series DSP chips are used for data processing, and the main frequency is 300MHZ. The peripheral devices of the system are integrated on the PCB backplane. Including LCD interface (display terminal), external expansion (DSP expansion slot), program download and online debugging interface (JATG), power interface, network communication interface (Ethnernet) and UART serial port.
在某个实施例中,图像处理模块对图像进行处理,获得图像采集设备的测角的具体过程为:In a certain embodiment, the image processing module processes the image, and the specific process of obtaining the angle measurement of the image acquisition device is as follows:
获取图像;get image;
提取图像的灰度矩阵,在灰度矩阵中寻找疑似标志点,将疑似标志点与实际标志点进行匹配,获得图像中的标志点位置;Extract the grayscale matrix of the image, find the suspected mark point in the grayscale matrix, match the suspected mark point with the actual mark point, and obtain the position of the mark point in the image;
根据相机成像原理对标志点位置进行逆向几何关系解算,得到图像采集设备的测角。该解算通过图像长、宽大小和相机视场角的对应比例关系,反向得到图像采集设备中标志点对应的的测角。According to the principle of camera imaging, the inverse geometric relationship of the position of the marker point is solved, and the angle measurement of the image acquisition device is obtained. The solution obtains the angle measurement corresponding to the mark point in the image acquisition device in reverse through the corresponding proportional relationship between the image length and width and the camera's field of view angle.
在某个实施例中,导航滤波模块中加载的导航滤波算法为卡尔曼滤波、扩展卡尔曼滤波、无迹卡尔曼滤波和相机偏心安装下的滤波算法中的一种。In a certain embodiment, the navigation filtering algorithm loaded in the navigation filtering module is one of Kalman filtering, extended Kalman filtering, unscented Kalman filtering, and filtering algorithms under camera eccentric installation.
在另一个实施例中,机动规划模块中加载的追逃博弈机动算法为比例导引律、最优制导律和微分对策制导律中的一种。In another embodiment, the maneuvering algorithm of the chase-and-flight game loaded in the maneuver planning module is one of a proportional guidance law, an optimal guidance law and a differential game guidance law.
本发明还提出一种航天器追逃博弈闭环半实物仿真方法,利用如上述所述的航天器追逃博弈闭环半实物仿真系统进行仿真。The present invention also proposes a closed-loop semi-physical simulation method for a spacecraft chase and escape game, which uses the above-mentioned closed-loop semi-physical simulation system for a spacecraft chase and escape game for simulation.
在某个实施例中,航天器追逃博弈闭环半实物仿真方法包括:In a certain embodiment, the closed-loop hardware-in-the-loop simulation method of the spacecraft chase and escape game includes:
101:根据相似性原理,设定仿真系统的缩比因子;101: According to the similarity principle, set the scaling factor of the simulation system;
102:根据缩比因子,将逃逸航天器和追踪航天器的实际初始位置和姿态换算为逃逸航天器模型和追踪航天器模型在仿真系统中的初始位置和姿态,并驱动逃逸航天器模型和追踪航天器模型运动至初始位置和姿态;102: According to the scaling factor, convert the actual initial position and attitude of the escape spacecraft and the tracking spacecraft into the initial position and attitude of the escape spacecraft model and the tracking spacecraft model in the simulation system, and drive the escape spacecraft model and tracking. The spacecraft model moves to the initial position and attitude;
103:利用相对测量模拟子系统模拟太空中的测量环境和测量手段对逃逸航天器模型和追踪航天器模型进行信息采集,获得逃逸航天器模型和追踪航天器模型的状态信息;103: Use the relative measurement simulation subsystem to simulate the measurement environment and measurement methods in space to collect information on the escape spacecraft model and the tracking spacecraft model, and obtain the status information of the escape spacecraft model and the tracking spacecraft model;
利用星载GNC模拟子系统根据状态信息,通过导航滤波算法进行滤波迭代,获得包含相对位置和相对速度在内的六维度状态信息,根据六维度状态信息,利用追逃博弈机动算法,获得追踪航天器模型和逃逸航天器模型的控制量,将控制量转化为第一控制指令;According to the state information, the spaceborne GNC simulation subsystem is used to filter and iterate through the navigation filtering algorithm to obtain six-dimensional state information including relative position and relative velocity. The control quantity of the spacecraft model and the escape spacecraft model is converted into the first control command;
利用闭环运行控制子系统根据相对状态对航天器轨道自然运动进行积分生成第二控制指令,对第一控制指令和第二控制指令进行加和获得综合指令,将综合指令解算获得逃逸航天器模型和追踪航天器模型的位置速度、姿态角和姿态角速度;The closed-loop operation control subsystem is used to integrate the natural motion of the spacecraft orbit according to the relative state to generate a second control command, add the first control command and the second control command to obtain a comprehensive command, and solve the comprehensive command to obtain an escape spacecraft model and track the position velocity, attitude angle and attitude angular velocity of the spacecraft model;
利用相对运动学模拟子系统根据位置速度、姿态角和姿态角速度驱动逃逸航天器模型和追踪航天器模型运动;Use the relative kinematics simulation subsystem to drive the escape spacecraft model and track the motion of the spacecraft model according to the position velocity, attitude angle and attitude angular velocity;
利用闭环运行控制子系统驱动相对测量模拟子系统、星载GNC模拟子系统和相对运动学模拟子系统,使逃逸航天器模型和追踪航天器模型继续运动,从而实现航天器追逃博弈闭环半实物仿真。The closed-loop operation control subsystem is used to drive the relative measurement simulation subsystem, the on-board GNC simulation subsystem and the relative kinematics simulation subsystem, so that the escape spacecraft model and the tracking spacecraft model continue to move, so as to realize the closed-loop semi-physical object of the spacecraft pursuit and escape game. simulation.
在某个实施例中,对于步骤101,根据相似性原理,设定仿真系统的缩比因子的步骤,包括:In a certain embodiment, for step 101, according to the similarity principle, the step of setting the scaling factor of the simulation system includes:
001:为了实现轨道博弈控制地面实验过程与空间真实过程完全相似,进而达到过程再现的目的,需要使地面验证实验模型与空间原型系统满足时间相似、运动相似、几何相似、力相似等。对于相似系统,相似准则在对应点和对应时刻上都保持相同的取值。即可,根据相似π定理,获得相对测量模拟子系统的测量过程和相对运动学模拟子系统的博弈过程中的相似准则为:001: In order to realize that the orbital game control ground experiment process is completely similar to the real process in space, and thus achieve the purpose of process reproduction, it is necessary to make the ground verification experimental model and the space prototype system satisfy the time similarity, motion similarity, geometric similarity, force similarity, etc. For similar systems, the similarity criterion maintains the same value at the corresponding point and at the corresponding time. That is, according to the similarity π theorem, the similarity criterion in the measurement process of the relative measurement simulation subsystem and the game process of the relative kinematic simulation subsystem is obtained as:
(πi)p=(πi)m,i=1,2,3,...,k-r (1)(π i ) p = (π i ) m , i = 1,2,3,...,kr (1)
式中,πi为变量;k为系统中所有变量个数;r为系统的基本量纲数;下标p为原型系统;下标m为模型系统;In the formula, π i is the variable; k is the number of all variables in the system; r is the basic dimension number of the system; the subscript p is the prototype system; the subscript m is the model system;
002:根据测量过程和博弈过程,将变量的缩比因子定义为模型系统变量与原型系统变量的比值:002: According to the measurement process and the game process, the scaling factor of the variable is defined as the ratio of the model system variable to the prototype system variable:
式中,λf为变量f的缩比因子;下标p为原型系统;下标m为模型系统;In the formula, λ f is the scaling factor of the variable f; the subscript p is the prototype system; the subscript m is the model system;
通过分析知,测量过程和博弈过程涉及到的基本量纲有三个:时间、长度和质量。Through analysis, there are three basic dimensions involved in the measurement process and the game process: time, length and quality.
003:根据步骤002定义时间、长度和质量三个基本量纲的缩比因子λt、λxn和λm,即和利用λt、λxn和λm,通过相似准则,获得其他变量的缩比因子应满足的关系式:003: Define the scaling factors λ t , λ xn and λ m of the three basic dimensions of time, length and quality according to step 002, namely and Using λ t , λ xn and λ m , through the similarity criterion, the relational expressions that the scaling factors of other variables should satisfy are obtained:
式中,λv,λa,λω,λα分别为速度,加速度,角速度,角加速度对应的缩比因子;λJ为转动惯量对应的缩比因子;λμ为地心引力常数对应的缩比因子。In the formula, λ v , λ a , λ ω , λ α are the scaling factors corresponding to the velocity, acceleration, angular velocity, and angular acceleration, respectively; λ J is the scaling factor corresponding to the moment of inertia; λ μ is the gravitational constant corresponding to Scaling factor.
本某个实施例中,在实际设计缩比因子时,缩比因子应同时满足以下设计准则:In this certain embodiment, when the scaling factor is actually designed, the scaling factor should satisfy the following design criteria at the same time:
缩比因子的选取应使得缩比后的距离尺度在相对运动学模拟子系统的三维尺度范围内;The selection of the scaling factor should make the scaled distance scale within the three-dimensional scale of the relative kinematics simulation subsystem;
缩比因子的选取应使得缩比后的航天器机动能力在相对运动学模拟子系统运动能力范围内;The selection of the scaling factor should make the scaled spacecraft maneuverability within the range of the motion capability of the relative kinematics simulation subsystem;
缩比因子的选取应满足逃逸航天器模型和追踪航天器模型的运动展示效果需求。The selection of the scaling factor should meet the needs of the motion display effect of the escape spacecraft model and the tracking spacecraft model.
实施例1Example 1
本实施例主要为了验证本发明提供的航天器追逃博弈闭环半实物仿真系统中导航滤波模块加载的导航滤波算法的效果。This embodiment is mainly for verifying the effect of the navigation filtering algorithm loaded by the navigation filtering module in the closed-loop hardware-in-the-loop simulation system of the spacecraft chase and escape game provided by the present invention.
本实施例中,导航滤波模块加载的导航滤波算法为卡尔曼滤波。In this embodiment, the navigation filtering algorithm loaded by the navigation filtering module is Kalman filtering.
(1)参数配置(1) Parameter configuration
设定实验的缩比因子:长度缩比因子为10000,时间缩比因子为1。追踪航天器初始轨道根数[a,e,i,ω,Ω,f]=[7020137.0m,0.0001°,33°,18°,12°,13.2°],a,e,i,ω,Ω,f分别为半长轴、偏心率、轨道倾角、升交点赤经、近地点角距,真近点角,逃逸航天器初始相对状态为:[x,y,z,vx,vy,vz]=[65000m,-1000m,4000m,-26m/s,8m/s,-11m/s],x,y,z,vx,vy,vz分别为LVLH坐标系下x,y,z三轴上的位置和速度,实验总时间为400s,实验仿真步长为0.5s,测量频率为50Hz。初始误差的协方差矩阵为:Set the scaling factor of the experiment: the length scaling factor is 10000, and the time scaling factor is 1. Tracking spacecraft initial orbital number [a,e,i,ω,Ω,f]=[7020137.0m, 0.0001°, 33°, 18°, 12°, 13.2°], a, e, i, ω, Ω , f are the semi-major axis, eccentricity, orbital inclination, ascending node right ascension, perigee angular distance, true perigee angle, respectively. The initial relative state of the escape spacecraft is: [x,y,z,v x ,v y ,v z ]=[65000m, -1000m, 4000m, -26m/s, 8m/s, -11m/s], x, y, z, v x , v y , v z are respectively x, y, The position and velocity on the z three-axis, the total experiment time is 400s, the experimental simulation step size is 0.5s, and the measurement frequency is 50Hz. The covariance matrix of the initial error is:
(2)待检验的扩展卡尔曼滤波(2) Extended Kalman filter to be tested
追逃博弈系统的状态方程和观测方程经离散后可写为:The state equation and observation equation of the chase-and-flight game system can be written as:
式中,状态量x,y,z为三坐标轴上的相对位置,为对应的相对速度;观测量zk=[r,α,β]T,r,α,β分别为逃逸航天器模型相对追踪航天器模型的距离、俯仰角、偏航角;k为仿真节点;f(·)和h(·)分别为状态方程和观测方程;up(xk-1)为追踪航天器模型的控制量;uE(xk-1)为逃逸航天器模型的控制量;w和v分别为过程噪声和观测噪声。In the formula, the state quantity x, y, z are the relative positions on the three coordinate axes, is the corresponding relative velocity; the observed value z k =[r,α,β] T , r,α,β are the distance, pitch angle and yaw angle of the escape spacecraft model relative to the tracking spacecraft model; k is the simulation node ; f(·) and h(·) are the state equation and observation equation, respectively; u p (x k-1 ) is the control quantity of the tracking spacecraft model; u E (x k-1 ) is the control of the escape spacecraft model quantity; w and v are process noise and observation noise, respectively.
观测方程满足:The observation equation satisfies:
式中,r,α,β分别为逃逸航天器模型相对追踪航天器模型的距离、俯仰角、偏航角,x,y,z为三坐标轴上的相对位置。In the formula, r, α, β are the distance, pitch angle and yaw angle of the escape spacecraft model relative to the tracking spacecraft model, respectively, x, y, z are relative positions on the three coordinate axes.
为了将扩展卡尔曼滤波用在非线性系统中,需要对非线性观测系统进行近似,扩展卡尔曼滤波对f和h进行泰勒展开,取一阶项,得到状态转移矩阵Φk和观测矩阵Hk。In order to use the extended Kalman filter in a nonlinear system, it is necessary to approximate the nonlinear observation system. The extended Kalman filter performs Taylor expansion on f and h, and takes the first-order term to obtain the state transition matrix Φ k and the observation matrix H k .
式中,符号含义同公式(4)和(5)。In the formula, the meanings of symbols are the same as those of formulas (4) and (5).
进一步,可得滤波的状态更新方程:Further, the state update equation of the filter can be obtained:
式中,为k-1节点的状态估计值;为预测的k节点的状态估计值;Φk为状态转移矩阵;Qk为过程噪声协方差矩阵;为k-1节点的估计误差协方差矩阵;为预测的k节点的估计误差协方差矩阵;为预测的k节点的观测值。In the formula, is the estimated value of the state of the k-1 node; is the predicted state estimate value of k nodes; Φ k is the state transition matrix; Q k is the process noise covariance matrix; is the estimated error covariance matrix of k-1 nodes; is the estimated error covariance matrix of the predicted k nodes; is the observed value of the predicted k nodes.
然后可以得到滤波的量测更新方程:Then the filtered measurement update equation can be obtained:
式中,“+”和“-”号分别为量测更新后和更新前;为k节点的状态估计值;zk为k节点的状态估计值;为k节点的估计误差协方差矩阵;I为单位矩阵;Rk为测量噪声协方差矩阵;Kk为卡尔曼增益;其他符号含义同公式(8)。In the formula, the "+" and "-" signs are respectively after the measurement update and before the update; is the estimated value of the state of node k; z k is the estimated value of the state of node k; is the estimated error covariance matrix of k nodes; I is the identity matrix; R k is the measurement noise covariance matrix; K k is the Kalman gain; other symbols have the same meanings as in formula (8).
从滤波算法中得到的的值为航天器模型经过滤波处理后的包含相对位置和相对速度在内的六维度状态信息。from the filtering algorithm The value of is the six-dimensional state information including relative position and relative velocity after filtering processing of the spacecraft model.
(3)实验结果(3) Experimental results
空间非合作目标导航滤波实验结果分别以V-bar、H-bar、R-bar三个方向的位置、速度、加速度以及位置误差、速度误差、加速度误差展示,如图8a~图8f所示。由图可知,经过卡尔曼滤波(经典滤波算法)对空间非合作目标进行导航滤波,收敛速度快,精度较高。同时,通过误差分析可知,由于引入了真实测量环节,滤波误差呈现非高斯特性。The experimental results of the space non-cooperative target navigation filtering are displayed in the three directions of V-bar, H-bar, R-bar, position, velocity, acceleration, and position error, velocity error, and acceleration error, respectively, as shown in Figures 8a-8f. It can be seen from the figure that the Kalman filter (classical filtering algorithm) is used to navigate and filter the space non-cooperative target, and the convergence speed is fast and the precision is high. At the same time, through the error analysis, it can be seen that due to the introduction of the real measurement link, the filtering error presents non-Gaussian characteristics.
实施例2Example 2
本实施例主要为了验证本发明提供的航天器追逃博弈闭环半实物仿真系统中机动规划模块加载的追逃博弈机动算法的效果。This embodiment is mainly for verifying the effect of the pursuit and escape game maneuver algorithm loaded by the maneuver planning module in the spacecraft pursuit and escape game closed-loop semi-physical simulation system provided by the present invention.
本实施例中,导航滤波模块加载的导航滤波算法为微分对策制导律。In this embodiment, the navigation filtering algorithm loaded by the navigation filtering module is a differential game guidance law.
(1)参数配置(1) Parameter configuration
设定实验的缩比因子:长度缩比因子为1000,时间缩比因子为20。初始时刻,逃逸航天器在400km高度的圆轨道上,追踪航天器在逃逸航天器5km附近,双方即将开始时长为的博弈。在虚拟动点LVLH坐标系下,两航天器模型的相对状态[x,y,z,vx,vy,vz]=[250m,-5000m,-200m,0m/s,0m/s,0m/s],x,y,z,vx,vy,vz分别为LVLH坐标系下x,y,z三轴上的位置和速度。仿真参数的配置见表1,实验结果见图9a~图9f。Set the scaling factor of the experiment: the length scaling factor is 1000, and the time scaling factor is 20. At the initial moment, the escape spacecraft is on a circular orbit at a height of 400km, and the tracking spacecraft is near the escape spacecraft 5km. The two sides are about to start a game of time. In the virtual moving point LVLH coordinate system, the relative states of the two spacecraft models [x, y, z, v x , v y , v z ]=[250m,-5000m,-200m,0m/s,0m/s, 0m/s], x, y, z, v x , v y , v z are the position and velocity on the x, y, z axes in the LVLH coordinate system, respectively. The configuration of the simulation parameters is shown in Table 1, and the experimental results are shown in Figure 9a to Figure 9f.
表1追逃博弈仿真实验参数配置Table 1 Parameter configuration of the simulation experiment of the chase and escape game
(2)待检验的微分对策制导律(2) Guidance law of differential game to be tested
两航天器模型的相对状态满足下述的运动方程:The relative states of the two spacecraft models satisfy the following equations of motion:
式中,x(t)为两航天器的相对状态量,uP和uE分别为追踪航天器模型和逃逸航天器模型的机动量;A和B为相对运动方程中的系数矩阵,假设参考轨道的角速度为ω,则二者的值满足:In the formula, x(t) is the relative state quantity of the two spacecraft, u P and u E are the maneuvering quantities of the tracking spacecraft model and the escape spacecraft model, respectively; A and B are the coefficient matrices in the relative motion equation, assuming the reference The angular velocity of the orbit is ω, then the values of the two satisfy:
定义微分对策中两航天器模型的支付函数为(下标P和E分别表示追踪航天器模型和逃逸航天器模型):The payoff function of the two spacecraft models in the differential game is defined as (the subscripts P and E represent the pursuit spacecraft model and the escape spacecraft model, respectively):
式中,Q,RP,RE分别为支付函数中的权重系数,JP和JE分别为追踪航天器模型和逃逸航天器模型的支付函数。In the formula, Q, R P , and RE are the weight coefficients in the payoff function, respectively, and J P and J E are the payoff functions of the tracking spacecraft model and the escape spacecraft model, respectively.
根据极大极小值原理,可得最优控制律:According to the maximum and minimum principle, the optimal control law can be obtained:
式中,KP和KE分别为追踪航天器和逃逸航天器的反馈增益;P为黎卡提矩阵,其他符号含义同公式(10)和(11)。In the formula, K P and K E are the feedback gains of the tracking spacecraft and the escape spacecraft, respectively; P is the Riccati matrix, and the meanings of other symbols are the same as formulas (10) and (11).
在仿真实验中,连续的追逃博弈系统需要被离散化为:In the simulation experiment, the continuous pursuit and flight game system needs to be discretized as:
式中,k为仿真节点;其他符号含义同公式(10)。In the formula, k is the simulation node; the meaning of other symbols is the same as formula (10).
进而,最优控制量的解算可以被离散为:Furthermore, the solution of the optimal control variable can be discretized as:
式中,值是节点k的滤波估计值;其他符号对应公式(12)中变量在节点k的值。In the formula, value is the filter estimate for node k; other symbols correspond to the value of the variable at node k in equation (12).
(3)实验结果(3) Experimental results
从图9a~图9f可以看出,半实物仿真轨迹整体上与数学仿真的轨迹吻合度较好,两航天器模型在100秒后达到相对速度和位置都接近零的交会状态,说明了基于该种微分对策交会制导方法在太空开展轨道追逃博弈的可行性较高,基本可以达到预期的博弈效果。仔细比较发现,在x轴(R-bar方向)上的轨迹误差会比y轴(V-bar方向)和z轴(H-bar方向)上略大,y轴(V-bar方向)和z轴(H-bar方向)上速度曲线的峰值会比理论曲线峰值要略大,这也体现了数学仿真和半实物仿真的结果差异。最后需要说明的是,为避免九自由度相对运动平台上两航天器模型的碰撞,仿真中预设了30cm的安全距离,博弈到达安全距离即交会结束。It can be seen from Figures 9a to 9f that the semi-physical simulation trajectories are in good agreement with the mathematical simulation trajectories as a whole, and the two spacecraft models reach the intersection state where the relative speed and position are close to zero after 100 seconds. This kind of differential game rendezvous guidance method is highly feasible to carry out orbital pursuit and escape game in space, and can basically achieve the expected game effect. Careful comparison found that the trajectory error on the x-axis (R-bar direction) will be slightly larger than that on the y-axis (V-bar direction) and z-axis (H-bar direction), and the y-axis (V-bar direction) and z The peak value of the velocity curve on the axis (H-bar direction) will be slightly larger than the peak value of the theoretical curve, which also reflects the difference between the results of mathematical simulation and hardware-in-the-loop simulation. Finally, it should be noted that in order to avoid the collision of the two spacecraft models on the nine-degree-of-freedom relative motion platform, a safety distance of 30cm is preset in the simulation, and the game will end when the game reaches the safety distance.
总体上,通过实施例2和实施例3,演示了本发明用于检验航天器追逃过程中的导航、制导算法的有效性的过程。In general, the process of the present invention for verifying the effectiveness of the navigation and guidance algorithms in the process of chasing and escaping a spacecraft is demonstrated through
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是在本发明的发明构思下,利用本发明说明书及附图内容所作的等效结构变换,或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。The above descriptions are only the preferred embodiments of the present invention, and are not intended to limit the scope of the present invention. Under the inventive concept of the present invention, the equivalent structural transformations made by the contents of the description and drawings of the present invention, or the direct/indirect application Other related technical fields are included in the scope of patent protection of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010571039.4A CN111679592B (en) | 2020-06-22 | 2020-06-22 | Spacecraft pursuit and escape game closed-loop semi-physical simulation system and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010571039.4A CN111679592B (en) | 2020-06-22 | 2020-06-22 | Spacecraft pursuit and escape game closed-loop semi-physical simulation system and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111679592A true CN111679592A (en) | 2020-09-18 |
CN111679592B CN111679592B (en) | 2023-04-07 |
Family
ID=72436195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010571039.4A Active CN111679592B (en) | 2020-06-22 | 2020-06-22 | Spacecraft pursuit and escape game closed-loop semi-physical simulation system and method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111679592B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112363410A (en) * | 2020-11-13 | 2021-02-12 | 浙江大学 | Intelligent autonomous control research and verification system for spacecraft |
CN112666984A (en) * | 2020-12-29 | 2021-04-16 | 北京电子工程总体研究所 | Aircraft pursuit and escape game method and system |
CN113311851A (en) * | 2021-04-25 | 2021-08-27 | 北京控制工程研究所 | Spacecraft pursuit-escape intelligent orbit control method and device and storage medium |
CN113325704A (en) * | 2021-04-25 | 2021-08-31 | 北京控制工程研究所 | Spacecraft backlight approaching intelligent orbit control method and device and storage medium |
CN115390449A (en) * | 2022-08-26 | 2022-11-25 | 西北工业大学 | Downstream position grabbing track game control method based on genetic algorithm |
CN116258317A (en) * | 2022-12-09 | 2023-06-13 | 哈尔滨工业大学 | A game theory-based satellite formation anti-surveillance optimal control method, device and medium |
CN117406618A (en) * | 2022-07-08 | 2024-01-16 | 南京理工大学 | Numerical simulation system and method for spacecraft formation and star group operation control |
CN118819004A (en) * | 2024-09-19 | 2024-10-22 | 成都流体动力创新中心 | A signal-level radar semi-physical simulation method and system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050060092A1 (en) * | 2003-08-05 | 2005-03-17 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
CN101122780A (en) * | 2007-09-04 | 2008-02-13 | 北京控制工程研究所 | Lunar soft landing guidance, navigation and control semi-physical simulation test system |
CN108519958A (en) * | 2018-02-05 | 2018-09-11 | 中国人民解放军国防科技大学 | A Method of Analyzing and Constructing Spacecraft Pursuit and Escape Boundary and Judging Capture and Escape Area |
CN109238287A (en) * | 2018-09-06 | 2019-01-18 | 中国人民解放军国防科技大学 | Spacecraft escape path planning method and system |
GB201910670D0 (en) * | 2018-07-25 | 2019-09-11 | Univ Northwestern Polytechnical | Neural network-based method for pursuit-evasion game of spacecrafts |
CN110673486A (en) * | 2019-10-22 | 2020-01-10 | 北京航空航天大学 | Multi-spacecraft pursuit and escape control method based on dynamic game theory |
-
2020
- 2020-06-22 CN CN202010571039.4A patent/CN111679592B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050060092A1 (en) * | 2003-08-05 | 2005-03-17 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
CN101122780A (en) * | 2007-09-04 | 2008-02-13 | 北京控制工程研究所 | Lunar soft landing guidance, navigation and control semi-physical simulation test system |
CN108519958A (en) * | 2018-02-05 | 2018-09-11 | 中国人民解放军国防科技大学 | A Method of Analyzing and Constructing Spacecraft Pursuit and Escape Boundary and Judging Capture and Escape Area |
GB201910670D0 (en) * | 2018-07-25 | 2019-09-11 | Univ Northwestern Polytechnical | Neural network-based method for pursuit-evasion game of spacecrafts |
CN109238287A (en) * | 2018-09-06 | 2019-01-18 | 中国人民解放军国防科技大学 | Spacecraft escape path planning method and system |
CN110673486A (en) * | 2019-10-22 | 2020-01-10 | 北京航空航天大学 | Multi-spacecraft pursuit and escape control method based on dynamic game theory |
Non-Patent Citations (1)
Title |
---|
ZHEN-YU LI ETC: "A dimension-reduction solution of free-time differential games for spacecraft pursuit-evasion" * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112363410A (en) * | 2020-11-13 | 2021-02-12 | 浙江大学 | Intelligent autonomous control research and verification system for spacecraft |
CN112363410B (en) * | 2020-11-13 | 2022-09-30 | 浙江大学 | Intelligent autonomous control research and verification system for spacecraft |
CN112666984A (en) * | 2020-12-29 | 2021-04-16 | 北京电子工程总体研究所 | Aircraft pursuit and escape game method and system |
CN113311851A (en) * | 2021-04-25 | 2021-08-27 | 北京控制工程研究所 | Spacecraft pursuit-escape intelligent orbit control method and device and storage medium |
CN113325704A (en) * | 2021-04-25 | 2021-08-31 | 北京控制工程研究所 | Spacecraft backlight approaching intelligent orbit control method and device and storage medium |
CN113325704B (en) * | 2021-04-25 | 2023-11-10 | 北京控制工程研究所 | Spacecraft backlight approach intelligent orbit control method, device and storage medium |
CN117406618A (en) * | 2022-07-08 | 2024-01-16 | 南京理工大学 | Numerical simulation system and method for spacecraft formation and star group operation control |
CN115390449A (en) * | 2022-08-26 | 2022-11-25 | 西北工业大学 | Downstream position grabbing track game control method based on genetic algorithm |
CN116258317A (en) * | 2022-12-09 | 2023-06-13 | 哈尔滨工业大学 | A game theory-based satellite formation anti-surveillance optimal control method, device and medium |
CN116258317B (en) * | 2022-12-09 | 2023-10-31 | 哈尔滨工业大学 | Satellite formation anti-monitoring optimal control method, device and medium based on game theory |
CN118819004A (en) * | 2024-09-19 | 2024-10-22 | 成都流体动力创新中心 | A signal-level radar semi-physical simulation method and system |
Also Published As
Publication number | Publication date |
---|---|
CN111679592B (en) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111679592B (en) | Spacecraft pursuit and escape game closed-loop semi-physical simulation system and method | |
CN100495261C (en) | Lunar soft landing guidance, navigation and control semi-physical simulation test system | |
CN101726296B (en) | Space robot vision measurement, path planning, GNC integrated simulation system | |
CN102759358B (en) | Relative posture dynamics modeling method based on dead satellite surface reference points | |
CN106681170B (en) | Semi-physical guidance simulation method and system | |
CN114488848B (en) | Unmanned aerial vehicle autonomous flight system and simulation experiment platform for indoor building space | |
Tsiotras | ASTROS: A 5DOF experimental facility for research in space proximity operations | |
CN103344256B (en) | A laboratory test method for multi-field star sensors | |
CN101794527B (en) | Manual control rendezvous and docking semi-physical simulation testing system | |
Benninghoff et al. | Development and hardware-in-the-loop test of a guidance, navigation and control system for on-orbit servicing | |
Cho et al. | A 5-dof experimental platform for spacecraft rendezvous and docking | |
CN109445310A (en) | A kind of infrared imaging guidance hardware-in-the-loop system and method | |
CN106227066A (en) | The configurable visual display method of spacecraft space manipulation | |
CN112650076A (en) | Constellation cooperative control ground simulation system | |
Fravolini et al. | Development of modelling and control tools for aerial refueling for uavs | |
CN111872938A (en) | Spatial three-dimensional large-scale kinematics simulation system and method | |
CN110851232B (en) | Space mission presentation device based on dysmorphism screen | |
CN115752514A (en) | Autonomous optical navigation verification system for deep space small celestial body | |
CN111537000A (en) | Ground-based verification system and method for optical navigation algorithm in the landing segment of small celestial bodies in deep space | |
CN117075495A (en) | A ground semi-physical simulation system based on multi-spacecraft attitude control | |
Wilde et al. | Experimental evaluation of a cots time-of-flight camera as rendezvous sensor for small satellites | |
CN103606332A (en) | Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof | |
CN107449419A (en) | The Full Parameterized vision measuring method of the continuous kinematic parameter of body target | |
CN114859949A (en) | Satellite remote repair ground simulation device and method | |
CN114114961A (en) | Space non-cooperative target optical tracking and aiming closed-loop verification system |
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 |