CN109815606B - 误差链对碎片碰撞清除成功率影响分析方法、系统 - Google Patents
误差链对碎片碰撞清除成功率影响分析方法、系统 Download PDFInfo
- Publication number
- CN109815606B CN109815606B CN201910091555.4A CN201910091555A CN109815606B CN 109815606 B CN109815606 B CN 109815606B CN 201910091555 A CN201910091555 A CN 201910091555A CN 109815606 B CN109815606 B CN 109815606B
- Authority
- CN
- China
- Prior art keywords
- error
- clearing
- collision
- comprehensive
- calculating
- 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
Classifications
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种误差链对碎片碰撞清除成功率影响分析方法、系统，包括如下步骤：获取误差项目步骤：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；计算综合误差步骤：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；计算清除概率步骤：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率。本发明通过分解影响碰撞清除成功的各项误差，计算误差的量级和具体数值，并通过数值模拟，计算碰撞清除成功概率，可以论证清除器的能力作用范围，并约束各部件的技术指标要求。
Description
技术领域
本发明涉及航天空间碎片清除任务技术领域，具体地，涉及误差链对碎片碰撞清除成功率影响分析方法、系统。尤其涉及清除误差链对碎片碰撞清除成功率影响分析方法。
背景技术
空间目标与碎片是人类太空航天器面临的巨大威胁，来自外太空的彗星等接近天体是地球面临的重大威胁，因此开发一种动能清除系统既可用于对大尺寸的太空碎片进行清除，又可以对来袭的彗星等不明天体进行碰撞摧毁，以瓦解目标，或者使其改变既有运行轨迹，达到防护航天器和保护地球的目的，过程如附图2所示。进行上述任务有网捕、机械臂捕捉、激光清除或者动能清除等方式，利用空间平台飞行清除物体进行太空垃圾的清除或太空威胁天体的清除代表一种先进的未来发展方向，由于不包含末制导系统，需要控制目标跟踪瞄准、平台指向、飞行方向等一系列影响碰撞清除成功率的因素。
例如，专利文献CN107719708A(申请号：201710890821.0)公开了一种空间碎片的捕获及清除装置及其方法，包括卫星组，卫星组通过导电缆绳与捕获网连接；所述卫星组包括多个小卫星，小卫星的缆绳收放装置通过导电缆绳与捕获网连接，卫星组通过导电缆绳带动捕获网和碎片一起旋转；所述捕获网包括上层网和下层网，其中下层网为底部封闭结构，上层网上设置有能够收缩/张开的捕获网入口，导电缆绳与下层网的底部连接。该装置采用多个小卫星实现了空间碎片的捕获和清除，其成本低，捕获效率高，并且有效减少燃料消耗。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种误差链对碎片碰撞清除成功率影响分析方法、系统。
根据本发明提供的一种误差链对碎片碰撞清除成功率影响分析方法，包括如下步骤：
获取误差项目步骤：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；
计算综合误差步骤：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；
计算清除概率步骤：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率。
优选地，所述误差源类别包括：目标位置预报误差、指向误差、清除器方向误差以及时间延迟误差；
目标位置预报误差包括：测量误差、轨道模型误差；测量误差包括：测角误差、测距误差、测速误差；
指向误差包括：平台指向精度误差、平台稳定度误差；
清除器方向误差包括：清除器初始飞行角度误差、飞行速度误差；
时间延迟误差包括：清除器飞行时延；
误差源的数值包括：不同清除距离下的误差源数值。
优选地，所述误差源的综合误差包括：伴随清除模式下不同清除距离的第一综合误差、交会清除模式下不同清除距离的第二综合误差；
所述计算综合误差步骤包括：
根据获取的误差源的类别、误差源的量级以及所述不同清除距离下的误差源数值，计算伴随清除模式下不同清除距离的第一综合误差及第一综合误差相应的清除偏差：
γx表示第一综合误差；
γN表示第N种误差源下的不同清除距离下的误差源数值，其中，N为正整数；
ρ表示清除距离；
根据获取的误差源的类别、误差源的量级以及误差源的数值，计算交会清除模式下不同清除距离的第二综合误差及第二综合误差相应的清除偏差：
γy表示第二综合误差；
γx1表示飞行时延误差；
γx2表示速度误差；
优选地，第一综合误差
中，N＝7，γ1为测角误差，γ2为测距误差，γ3为测速误差，γ4为指向精度，γ5为稳定度或时延，γ6为角度误差，γ7为轨道模型误差；
飞行时延误差为：
γx1＝v·dt
速度误差为：
v表示飞行速度；
dt表示时延变量；
d表示与目标距离；
vmin表示飞行速度最小可能值；
vmax表示飞行速度最大可能值。
优选地，所述计算清除概率步骤包括：
根据伴随清除模式下不同清除距离的第一综合误差的值，计算伴随清除模式下不同清除距离的碰撞清除成功概率；
根据交会清除模式下不同清除距离的第二综合误差的值，计算伴随交会模式下不同清除距离的碰撞清除成功概率。
根据本发明提供的一种误差链对碎片碰撞清除成功率影响分析系统，包括如下模块：
获取误差项目模块：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；
计算综合误差模块：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；
计算清除概率模块：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率。
优选地，所述误差源类别包括：目标位置预报误差、指向误差、清除器方向误差以及时间延迟误差；
目标位置预报误差包括：测量误差、轨道模型误差；测量误差包括：测角误差、测距误差、测速误差；
指向误差包括：平台指向精度误差、平台稳定度误差；
清除器方向误差包括：清除器初始飞行角度误差、飞行速度误差；
时间延迟误差包括：清除器飞行时延；
误差源的数值包括：不同清除距离下的误差源数值。
优选地，所述误差源的综合误差包括：伴随清除模式下不同清除距离的第一综合误差、交会清除模式下不同清除距离的第二综合误差；
所述计算综合误差模块包括：
根据获取的误差源的类别、误差源的量级以及所述不同清除距离下的误差源数值，计算伴随清除模式下不同清除距离的第一综合误差及第一综合误差相应的清除偏差：
γx表示第一综合误差；
γN表示第N种误差源下的不同清除距离下的误差源数值，其中，N为正整数；
ρ表示清除距离；
根据获取的误差源的类别、误差源的量级以及误差源的数值，计算交会清除模式下不同清除距离的第二综合误差及第二综合误差相应的清除偏差：
γy表示第二综合误差；
γx1表示飞行时延误差；
γx2表示速度误差；
优选地，所述第一综合误差
中，N＝7，γ1为测角误差，γ2为测距误差，γ3为测速误差，γ4为指向精度，γ5为稳定度或时延，γ6为角度误差，γ7为轨道模型误差；
飞行时延误差为：
γx1＝v·dt
速度误差为：
v表示飞行速度；
dt表示时延变量；
d表示与目标距离；
vmin表示飞行速度最小可能值；
vmax表示飞行速度最大可能值。
优选地，所述计算清除概率模块包括：
根据伴随清除模式下不同清除距离的第一综合误差的值，计算伴随清除模式下不同清除距离的碰撞清除成功概率；
根据交会清除模式下不同清除距离的第二综合误差的值，计算伴随交会模式下不同清除距离的碰撞清除成功概率。
与现有技术相比，本发明具有如下的有益效果：
本发明通过分解影响碰撞清除成功的各项误差，计算误差的量级和具体数值，并通过数值模拟，计算碰撞清除成功概率，可以论证清除器的能力作用范围，并约束各部件的技术指标要求。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明提供的误差链对碎片碰撞清除成功率影响分析方法的流程步骤示意图。
图2为本发明提供的目标面碰撞清除过程示意图。
图3为本发明提供的位置误差示意图。
图4为本发明提供的影响碰撞清除成功率的误差树示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步地理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示，根据本发明提供的一种误差链对碎片碰撞清除成功率影响分析方法，包括如下步骤：
获取误差项目步骤：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；进一步地，建立影响碰撞清除成功率的误差模型方法包括：识别测量系统及执行系统的误差因素，根据误差数值及动态清除模型，模拟清除过程，利用蒙特卡洛法计算成功概率；
计算综合误差步骤：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；
计算清除概率步骤：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率。进一步地，模拟清除过程，利用蒙特卡洛法计算碰撞清除成功概率，确定每一个飞行清除器的碰撞位置；根据角度误差总偏差，利用正态分布函数生成点坐标，模拟105个点，统计碰撞数量，计算碰撞概率。
具体地，所述误差源类别包括：目标位置预报误差、指向误差、清除器方向误差以及时间延迟误差；
目标位置预报误差包括：测量误差、轨道模型误差；测量误差包括：测角误差、测距误差、测速误差；
指向误差包括：平台指向精度误差、平台稳定度误差；
清除器方向误差包括：清除器初始飞行角度误差、飞行速度误差；
时间延迟误差包括：清除器飞行时延；
误差源的数值包括：不同清除距离下的误差源数值。
具体地，所述误差源的综合误差包括：伴随清除模式下不同清除距离的第一综合误差、交会清除模式下不同清除距离的第二综合误差；进一步地，伴随清除模式下测量误差主要体现在角度测量精度，影响清除器指向，对最终碰撞清除成功率的影响随着清除距离的增加线性增加；轨道模型误差表示清除解算误差，即清除器指向角度解算的误差大小，影响最终碰撞清除成功率；平台指向误差、清除器飞行角度误差对碰撞清除成功率的影响与测量角度原理相同；平台稳定度与时间迟延误差综合影响飞行状态清除器指向；时延误差与目标相对运动速度的累积作用影响最终碰撞清除成功率，在伴随清除模式下的影响较小；清除器速度误差通过影响碰撞时间，影响最终碰撞清除成功率，在伴随清除模式下的影响较小。更进一步地，相对伴随清除模式，时延迟误差与目标相对运动速度的累积作用影响最终碰撞清除成功率，交会清除模式下目标运动速度较快，是影响碰撞清除成功率的主要因素。
所述计算综合误差步骤包括：
根据获取的误差源的类别、误差源的量级以及所述不同清除距离下的误差源数值，计算伴随清除模式下不同清除距离的第一综合误差及第一综合误差相应的清除偏差：
γx表示第一综合误差；
γN表示第N种误差源下的不同清除距离下的误差源数值，其中，N为正整数；
ρ表示清除距离；
根据获取的误差源的类别、误差源的量级以及误差源的数值，计算交会清除模式下不同清除距离的第二综合误差及第二综合误差相应的清除偏差：
γy表示第二综合误差；
γx1表示飞行时延误差；
γx2表示速度误差；
具体地，第一综合误差
中，N＝7，γ1为测角误差，γ2为测距误差，γ3为测速误差，γ4为指向精度，γ5为稳定度或时延，γ6为角度误差，γ7为轨道模型误差；
飞行时延误差为：
γx1＝v·dt
速度误差为：
v表示飞行速度；
dt表示时延变量；
d表示与目标距离；
vmin表示飞行速度最小可能值；
vmax表示飞行速度最大可能值。
具体地，所述计算清除概率步骤包括：
根据伴随清除模式下不同清除距离的第一综合误差的值，计算伴随清除模式下不同清除距离的碰撞清除成功概率；
根据交会清除模式下不同清除距离的第二综合误差的值，计算伴随交会模式下不同清除距离的碰撞清除成功概率。
本发明提供的误差链对碎片碰撞清除成功率影响分析系统，可以通过本发明给的误差链对碎片碰撞清除成功率影响分析方法的步骤流程实现。本领域技术人员可以将所述误差链对碎片碰撞清除成功率影响分析方法，理解为所述误差链对碎片碰撞清除成功率影响分析系统的一个优选例。
根据本发明提供的一种误差链对碎片碰撞清除成功率影响分析系统，包括如下模块：
获取误差项目模块：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；
计算综合误差模块：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；
计算清除概率模块：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率。
具体地，所述误差源类别包括：目标位置预报误差、指向误差、清除器方向误差以及时间延迟误差；
目标位置预报误差包括：测量误差、轨道模型误差；测量误差包括：测角误差、测距误差、测速误差；
指向误差包括：平台指向精度误差、平台稳定度误差；
清除器方向误差包括：清除器初始飞行角度误差、飞行速度误差；
时间延迟误差包括：清除器飞行时延；
误差源的数值包括：不同清除距离下的误差源数值。
具体地，所述误差源的综合误差包括：伴随清除模式下不同清除距离的第一综合误差、交会清除模式下不同清除距离的第二综合误差；
所述计算综合误差模块包括：
根据获取的误差源的类别、误差源的量级以及所述不同清除距离下的误差源数值，计算伴随清除模式下不同清除距离的第一综合误差及第一综合误差相应的清除偏差：
γx表示第一综合误差；
γN表示第N种误差源下的不同清除距离下的误差源数值，其中，N为正整数；
ρ表示清除距离；
根据获取的误差源的类别、误差源的量级以及误差源的数值，计算交会清除模式下不同清除距离的第二综合误差及第二综合误差相应的清除偏差：
γy表示第二综合误差；
γx1表示飞行时延误差；
γx2表示速度误差；
具体地，所述第一综合误差
中，N＝7，γ1为测角误差，γ2为测距误差，γ3为测速误差，γ4为指向精度，γ5为稳定度或时延，γ6为角度误差，γ7为轨道模型误差；
飞行时延误差为：
γx1＝v·dt
速度误差为：
v表示飞行速度；
dt表示时延变量；
d表示与目标距离；
vmin表示飞行速度最小可能值；
vmax表示飞行速度最大可能值。
具体地，所述计算清除概率模块包括：
根据伴随清除模式下不同清除距离的第一综合误差的值，计算伴随清除模式下不同清除距离的碰撞清除成功概率；
根据交会清除模式下不同清除距离的第二综合误差的值，计算伴随交会模式下不同清除距离的碰撞清除成功概率。
下面通过优选例，对本发明进行更为具体的说明。
优选例1：
一种清除误差链对碎片碰撞清除成功率影响分析方法,包括：
1、统计误差项目：
包括目标位置预报误差：相对测量误差，测距误差、测角误差以及测速误差及轨道模型误差；
包括指向误差，平台指向误差及稳定度误差；
包括清除器飞行的角度和速度误差；
包括时间延迟误差：清除器飞行时延。
2、计算误差源数值大小。
3、进行模拟碰撞清除成功概率仿真。
优选例2：
一种清除误差链对碎片碰撞清除成功率影响分析方法，影响碰撞清除成功率的主要误差分为四部分。具体包括目标相对测量预测误差，包括测量测角、测距、测速误差、轨道模型误差；平台指向误差和稳定度误差；清除器速度误差和角度误差；时延误差。
误差源及数量级如下所示：
相对轨道模型误差随清除距离变化情况如下：
确定了误差源的项目和指标后，需设计清除器的使用模式，这里主要分为伴随清除模式和交会清除模式，误差对碰撞清除成功率的影响在两个模式下略有不同。
伴随清除模式影响因素说明：
测量误差主要体现在角度测量精度，影响清除器指向，对最终碰撞清除成功率的影响随着清除距离的增加线性增加；轨道模型误差表示清除解算误差，影响最终碰撞清除成功率；平台指向误差、清除器飞行角度误差对碰撞清除成功率的影响与测量角度原理相同；平台稳定度与时间迟延误差综合影响飞行状态清除器指向；时延误差与目标相对运动速度的累积作用影响最终碰撞清除成功率，伴随模式下影响较小；清除器速度误差通过影响碰撞时间，影响最终碰撞清除成功率，伴随模式影响较小。
交会清除模式影响因素说明：
相对伴随清除模式，时延迟误差与目标相对运动速度的累积作用影响最终碰撞清除成功率，交会清除模式下目标运动速度较快，是影响碰撞清除成功率的主要因素。
优选例3：
一种清除误差链对碎片碰撞清除成功率影响分析方法，包括：
a)伴随清除误差统计及碰撞清除成功概率计算
在伴随清除模式下，各误差源的大小及对偏差量的贡献度如下表所示：
不同清除距离下误差源对清除偏差影响分析表(单位:m)如下：
| 清除距离 | 1km | 3km | 5km | 10km | 15km | 20km | 50km |
| 测角误差 | 1.745 | 5.235 | 8.726 | 17.453 | 26.179 | 34.906 | 87.266 |
| 测距误差 | 0.002 | 0.006 | 0.008 | 0.02 | 0.03 | 0.04 | 0.1 |
| 测速误差 | 0.005 | 0.016 | 0.025 | 0.06 | 0.09 | 0.11 | 0.25 |
| 指向精度 | 1.745 | 5.235 | 8.726 | 17.453 | 26.179 | 34.906 | 87.266 |
| 稳定度/时延 | 0.004 | 0.013 | 0.021 | 0.043 | 0.065 | 0.087 | 0.218 |
| 角度误差 | 1.745 | 5.235 | 8.726 | 17.453 | 26.179 | 34.906 | 87.266 |
| 轨道模型误差 | 0.5 | 0.8 | 1 | 1.5 | 2 | 2.5 | 5 |
| 综合误差 | 3.064 | 9.104 | 15.148 | 30.267 | 45.389 | 60.512 | 151.233 |
令各项误差为角度方向对应分别γ1到γN，则这几类误差造成的综合误差为
则相应的清除偏差为一个半径为δr＝γρ的圆，其中ρ为清除距离。根据偏离量大小，利用蒙特卡洛法，计算不同清除距离及清除模式下的碰撞清除成功概率如下表所示。
不同目标尺寸下碰撞清除成功概率如下：
b)交会清除误差统计及碰撞清除成功概率计算
相比于伴随清除，除了测量误差以外因为时延造成的误差和速度造成的误差影响较大。
不同清除距离下误差源对清除偏差影响分析(单位:m)
令角度综合误差为γ1，飞行时延误差为γ2，速度误差为γ3，则这几类误差造成的综合误差为
时间延迟上对应的偏离量为γ2＝v·dt；速度误差造成误差为/>
则相应的清除偏差为一个半径为δr＝γρ的圆，其中ρ为清除距离，如附图3所示。根据偏离量大小，利用蒙特卡洛法，计算不同清除距离及清除模式下的碰撞清除成功概率如下表所示。
对不同大小物体交会碰撞清除成功概率
通过两种模式的清除仿真，对误差链对碰撞清除成功率的影响进行偏离误差量的计算，并对碰撞清除成功的概率进行蒙特卡洛方法的仿真和计算。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。
Claims (8)
1.一种误差链对碎片碰撞清除成功率影响分析方法，其特征在于，包括如下步骤：
获取误差项目步骤：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；
计算综合误差步骤：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；
计算清除概率步骤：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率；
所述误差源类别包括：目标位置预报误差、指向误差、清除器方向误差以及时间延迟误差；
目标位置预报误差包括：测量误差、轨道模型误差；测量误差包括：测角误差、测距误差、测速误差；
指向误差包括：平台指向精度误差、平台稳定度误差；
清除器方向误差包括：清除器初始飞行角度误差、飞行速度误差；
时间延迟误差包括：清除器飞行时延；
误差源的数值包括：不同清除距离下的误差源数值。
2.根据权利要求1所述的误差链对碎片碰撞清除成功率影响分析方法，其特征在于，所述误差源的综合误差包括：伴随清除模式下不同清除距离的第一综合误差、交会清除模式下不同清除距离的第二综合误差；
所述计算综合误差步骤包括：
根据获取的误差源的类别、误差源的量级以及所述不同清除距离下的误差源数值，计算伴随清除模式下不同清除距离的第一综合误差及第一综合误差相应的清除偏差：
γx表示第一综合误差；
γN表示第N种误差源下的不同清除距离下的误差源数值，其中，N为正整数；
ρ表示清除距离；
根据获取的误差源的类别、误差源的量级以及误差源的数值，计算交会清除模式下不同清除距离的第二综合误差及第二综合误差相应的清除偏差：
γy表示第二综合误差；
γx1表示飞行时延误差；
γx2表示速度误差；
3.根据权利要求2所述的误差链对碎片碰撞清除成功率影响分析方法，其特征在于，第一综合误差
中，N＝7，γ1为测角误差，γ2为测距误差，γ3为测速误差，γ4为指向精度，γ5为稳定度或时延，γ6为角度误差，γ7为轨道模型误差；
飞行时延误差为：
γx1＝v·dt
速度误差为：
v表示飞行速度；
dt表示时延变量；
d表示与目标距离；
vmin表示飞行速度最小可能值；
vmax表示飞行速度最大可能值。
4.根据权利要求3所述的误差链对碎片碰撞清除成功率影响分析方法，其特征在于，所述计算清除概率步骤包括：
根据伴随清除模式下不同清除距离的第一综合误差的值，计算伴随清除模式下不同清除距离的碰撞清除成功概率；
根据交会清除模式下不同清除距离的第二综合误差的值，计算伴随交会模式下不同清除距离的碰撞清除成功概率。
5.一种误差链对碎片碰撞清除成功率影响分析系统，其特征在于，包括如下模块：
获取误差项目模块：获取误差源的类别、误差源的量级以及误差源的数值，建立影响碰撞清除成功率的误差模型；
计算综合误差模块：根据影响碰撞清除成功率的误差模型，计算误差源的综合误差；
计算清除概率模块：根据影响碰撞清除成功率的误差模型及误差源的综合误差，模拟清除过程，计算碰撞清除成功概率；
所述误差源类别包括：目标位置预报误差、指向误差、清除器方向误差以及时间延迟误差；
目标位置预报误差包括：测量误差、轨道模型误差；测量误差包括：测角误差、测距误差、测速误差；
指向误差包括：平台指向精度误差、平台稳定度误差；
清除器方向误差包括：清除器初始飞行角度误差、飞行速度误差；
时间延迟误差包括：清除器飞行时延；
误差源的数值包括：不同清除距离下的误差源数值。
6.根据权利要求5所述的误差链对碎片碰撞清除成功率影响分析系统，其特征在于，所述误差源的综合误差包括：伴随清除模式下不同清除距离的第一综合误差、交会清除模式下不同清除距离的第二综合误差；
所述计算综合误差模块包括：
根据获取的误差源的类别、误差源的量级以及所述不同清除距离下的误差源数值，计算伴随清除模式下不同清除距离的第一综合误差及第一综合误差相应的清除偏差：
γx表示第一综合误差；
γN表示第N种误差源下的不同清除距离下的误差源数值，其中，N为正整数；
ρ表示清除距离；
根据获取的误差源的类别、误差源的量级以及误差源的数值，计算交会清除模式下不同清除距离的第二综合误差及第二综合误差相应的清除偏差：
γy表示第二综合误差；
γx1表示飞行时延误差；
γx2表示速度误差；
7.根据权利要求6所述的误差链对碎片碰撞清除成功率影响分析系统，其特征在于，所述第一综合误差
中，N＝7，γ1为测角误差，γ2为测距误差，γ3为测速误差，γ4为指向精度，γ5为稳定度或时延，γ6为角度误差，γ7为轨道模型误差；
飞行时延误差为：
γx1＝v·dt
速度误差为：
v表示飞行速度；
dt表示时延变量；
d表示与目标距离；
vmin表示飞行速度最小可能值；
vmax表示飞行速度最大可能值。
8.根据权利要求7所述的误差链对碎片碰撞清除成功率影响分析系统，其特征在于，所述计算清除概率模块包括：
根据伴随清除模式下不同清除距离的第一综合误差的值，计算伴随清除模式下不同清除距离的碰撞清除成功概率；
根据交会清除模式下不同清除距离的第二综合误差的值，计算伴随交会模式下不同清除距离的碰撞清除成功概率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910091555.4A CN109815606B (zh) | 2019-01-30 | 2019-01-30 | 误差链对碎片碰撞清除成功率影响分析方法、系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910091555.4A CN109815606B (zh) | 2019-01-30 | 2019-01-30 | 误差链对碎片碰撞清除成功率影响分析方法、系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109815606A CN109815606A (zh) | 2019-05-28 |
| CN109815606B true CN109815606B (zh) | 2023-03-24 |
Family
ID=66605831
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910091555.4A Active CN109815606B (zh) | 2019-01-30 | 2019-01-30 | 误差链对碎片碰撞清除成功率影响分析方法、系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109815606B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111428369A (zh) * | 2020-03-26 | 2020-07-17 | 中国人民解放军32035部队 | 一种空间目标碰撞预警结果置信度计算方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104537230A (zh) * | 2014-12-23 | 2015-04-22 | 中国科学院国家天文台 | 一种航天器发射预警碰撞风险分析方法和分析装置 |
| CN105243278A (zh) * | 2015-10-20 | 2016-01-13 | 西北工业大学 | 一种基于概率理论的航天器智能预警方法 |
| CN106055780A (zh) * | 2016-05-26 | 2016-10-26 | 北京航空航天大学 | 一种基于gpu加速的空间碎片预警方法 |
| CN106845700A (zh) * | 2017-01-09 | 2017-06-13 | 西北工业大学 | 一种多碎片主动清除离线最优导引算法 |
| CN107450312A (zh) * | 2017-07-06 | 2017-12-08 | 南京航空航天大学 | 考虑航天器尺寸的防碰撞方法 |
| RU2643020C1 (ru) * | 2016-09-06 | 2018-01-29 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Омский государственный технический университет" | Способ проведения лётно-конструкторских испытаний автономного стыковочного модуля для очистки орбит от космического мусора |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130082146A1 (en) * | 2011-09-30 | 2013-04-04 | Raytheon Bbn Technologies Corp. | System and method for creating an artificial atmosphere for the removal of space debris |
| KR101537301B1 (ko) * | 2013-10-28 | 2015-07-20 | 한국항공우주연구원 | Ｃｓｍ 기반 충돌위험 분석 시스템 |
-
2019
- 2019-01-30 CN CN201910091555.4A patent/CN109815606B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104537230A (zh) * | 2014-12-23 | 2015-04-22 | 中国科学院国家天文台 | 一种航天器发射预警碰撞风险分析方法和分析装置 |
| CN105243278A (zh) * | 2015-10-20 | 2016-01-13 | 西北工业大学 | 一种基于概率理论的航天器智能预警方法 |
| CN106055780A (zh) * | 2016-05-26 | 2016-10-26 | 北京航空航天大学 | 一种基于gpu加速的空间碎片预警方法 |
| RU2643020C1 (ru) * | 2016-09-06 | 2018-01-29 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Омский государственный технический университет" | Способ проведения лётно-конструкторских испытаний автономного стыковочного модуля для очистки орбит от космического мусора |
| CN106845700A (zh) * | 2017-01-09 | 2017-06-13 | 西北工业大学 | 一种多碎片主动清除离线最优导引算法 |
| CN107450312A (zh) * | 2017-07-06 | 2017-12-08 | 南京航空航天大学 | 考虑航天器尺寸的防碰撞方法 |
Non-Patent Citations (3)
| Title |
|---|
| Analysis of the effect of attachment point bias during large space debris removal using a tethered space tug;Zhongyi Chu等;《Acta Astronautica》;20170628;全文 * |
| 空间机动飞行器清除碎片任务规划与制导控制技术;王伟林;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20190115;全文 * |
| 空间目标轨道预报误差与碰撞概率问题研究;白显宗;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20141215;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109815606A (zh) | 2019-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kersevan et al. | Recent developments of Monte-Carlo codes molflow+ and synrad+ | |
| CN103901907B (zh) | 一种软着陆避障模拟试验系统 | |
| CN104236546A (zh) | 一种卫星星光折射导航误差确定与补偿方法 | |
| CN106990401A (zh) | 基于全波形机载激光雷达数据二类高程误差修正方法 | |
| Ordaz et al. | Full-carpet design of a low-boom demonstrator concept | |
| CN109815606B (zh) | 误差链对碎片碰撞清除成功率影响分析方法、系统 | |
| CN114936471A (zh) | 一种基于并行计算的航天器碰撞预警分层快速筛选方法 | |
| CN104792299B (zh) | 一种基于测角数据的小行星轨道识别方法 | |
| CN104864873A (zh) | 一种利用人体运动特征辅助地图定位的方法 | |
| CN115446834B (zh) | 一种基于占据栅格配准的车底巡检机器人单轴重定位方法 | |
| CN107944084A (zh) | 一种深空探测器行星借力飞行轨道评估方法 | |
| CN104200030B (zh) | 一种圆参考轨道下给定边界的卫星初始相对状态确定方法 | |
| CN114780915A (zh) | 碰撞预警服务平台的数据正确性的确定方法及装置、设备 | |
| CN110231619B (zh) | 基于恩克法的雷达交接时刻预报方法及装置 | |
| Olivieri et al. | Simulations of satellites mock-up fragmentation | |
| CN102819667B (zh) | 基于时间约束库的飞行器协同攻击发射时序优化方法 | |
| CN104615008A (zh) | 一种自动紧急刹车aeb系统的测试结果分析方法和系统 | |
| CN112861373A (zh) | 近地小行星撞击轨道的生成方法和装置 | |
| Olivier | A Simulation and modeling framework for space situational awareness | |
| Feng et al. | Boosting the Robustness Verification of DNN by Identifying the Achilles's Heel | |
| CN116176869B (zh) | 一种基于先验信息的低轨机动目标初轨确定方法 | |
| CN103941593A (zh) | 低轨卫星姿态仿真方法 | |
| de la Vega et al. | Detection of outliers in sets of GNSS tracks from Volunteered Geographic Information | |
| Pekala | Atmospheric scattering of light emitted by extensive air showers | |
| Nair et al. | Sequence-based Encoding of Light Curves for Exoplanet Detection |
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 |