CN112966333A - 一种绳系飞网捕获非合作目标的碰撞中心估计方法 - Google Patents
一种绳系飞网捕获非合作目标的碰撞中心估计方法 Download PDFInfo
- Publication number
- CN112966333A CN112966333A CN202110235471.0A CN202110235471A CN112966333A CN 112966333 A CN112966333 A CN 112966333A CN 202110235471 A CN202110235471 A CN 202110235471A CN 112966333 A CN112966333 A CN 112966333A
- Authority
- CN
- China
- Prior art keywords
- net
- representing
- collision
- rope
- unit
- 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
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000008859 change Effects 0.000 claims abstract description 18
- 238000013528 artificial neural network Methods 0.000 claims abstract description 17
- 239000002245 particle Substances 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 108050006002 RNA polymerase sigma factor FliA Proteins 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000009975 flexible effect Effects 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G4/00—Tools specially adapted for use in space
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Artificial Intelligence (AREA)
- Aviation & Aerospace Engineering (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Automation & Control Theory (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Remote Sensing (AREA)
- Financial Or Insurance-Related Operations Such As Payment And Settlement (AREA)
Abstract
本发明公开了一种绳系飞网捕获非合作目标的碰撞中心估计方法,首先建立空间绳系飞网机器人系统动力学模型,接下来设计可机动单元系绳拉力观测器,然后计算碰撞时刻可机动单元拉力突变峰值,最后采用基于BP神经网络的算法估计碰撞中心点位置。本发明不依靠非合作目标相对速度、质量等信息,可以快速直接通过输入可机动单元系绳拉力精确估计碰撞中心点。
Description
技术领域
本发明属于航天技术领域,具体涉及一种绳系飞网碰撞中心估计方法。
背景技术
空间绳系飞网具有灵活性和大包络性等优点,近年来,被广泛应用于空间碎片清除、废弃卫星捕获等在轨服务任务。其中,碰撞抓捕环节是决定任务成败的关键。由于绳系飞网柔性特质和特殊的弹射发射方式,以及目标运动状态通常无法测量或预知,整个抓捕过程中,飞网与目标的碰撞是随机的。碰撞中心的位置直接影响捕获收口的期望位置,对于后续的轨迹规划及控制具有重要意义。
然而,非合作目标的运动状态和姿态都是未知的,目标与空间绳系飞网之间没有信息交流,故碰撞抓捕过程具有任意性,从而使飞网的收拢问题成为研究的难点。针对这一问题,有很多专利和文章给出了一些解决方法。例如,申请号为CN109063307A的中国发明专利,提出了一种用于计算可机动飞网抓捕目标的碰撞位移和碰撞力计算方法,通过可机动单元的测量力、位置和速度信息得到飞网中心的碰撞位移和碰撞力。但这种方法并未直接给出碰撞中心点的位置信息。本发明在前有发明的基础上,在非合作目标的相对速度、质量和体积等未知条件下,仅通过突变时刻可机动单元系绳拉力峰值来估计碰撞中心点。
发明内容
为了克服现有技术的不足,本发明提供了一种绳系飞网捕获非合作目标的碰撞中心估计方法,首先建立空间绳系飞网机器人系统动力学模型,接下来设计可机动单元系绳拉力观测器,然后计算碰撞时刻可机动单元拉力突变峰值,最后采用基于BP神经网络的算法估计碰撞中心点位置。本发明不依靠非合作目标相对速度、质量等信息,可以快速直接通过输入可机动单元系绳拉力精确估计碰撞中心点。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:建立空间绳系飞网机器人系统动力学模型;
空间绳系飞网机器人系统包括四个可机动单元;
空间绳系飞网机器人系统采用弹簧-质杆动力学模型,空间绳系飞网的每个系绳节点的质量表示为:
其中,i表示飞网中的节点,mCM表示可机动单元质量,ρ是飞网绳段材料的密度,A是飞网绳段的横截面积,l是飞网中未伸长的网格边长;
飞网相邻质点间的系绳拉力矢量表示为:
其中,是从第i个节点到第j个节点的单位向量,ri是飞网相邻质点间实际长度;Fi=(k(l-l0)+βkv)表示飞网相邻质点间的系绳拉力大小,l0是每个飞网网格绳段的实际长度,k是飞网绳段的弹性系数,k=EA/l0,E是飞网绳段材料的杨氏模量,β是飞网绳段的阻尼系数,ξ是飞网绳段材料的阻尼比,ν是网格末节点所在轴向相对速度的映射;
根据牛顿第二定律得到飞网中每个质点的动力学方程为:
miai=∑Fi+Fext (3)
其中,ai是质点运动加速度,∑Fi是每个质点系绳拉力之和,Fext表示重力和碰撞力之和;
四个可机动单元的动力学模型为:
mCMaCM=Fj1+Fj2+Fext (4)
其中,Fj1和Fj2分别表示与可机动单元相连接的两段绳子的拉力;
步骤2:设计可机动单元系绳拉力观测器;
可机动单元非线性系统定义为:
设z为系统扰动估计误差,则
进一步得到四个可机动单元构成的非线性系统:
x2=e2+αq (9)
设计扰动估计自适应律w及额外补偿项β(x1,x2),使流形保持不变,对系统扰动估计误差z求导得:
定义扰动估计自适应律w为:
将式(12)带入式(11)中得:
步骤3:计算碰撞时刻可机动单元拉力突变峰值;
根据步骤2中观测得出的可机动单元系绳突变拉力,得到干扰力与拉力峰值分离表达式为:
Ff=f(xtd)=Fd+F (14)
Fd=f(x0) (15)
F=Ff-Fd (16)
其中,Ff表示可机动单元突变时刻峰值,f(xtp)表示拉力观测器,Fd表示耦合干扰力,F表示突变时刻拉力,f(x0)表示非碰撞时刻系绳拉力观测器;
步骤4:基于BP神经网络飞网碰撞训练过程如下:
假定空间绳系飞网为正方形,大小为Q*Q,飞网网格边长为S;在飞网中取P个样本点;
训练样本集Fp为:
Fp=[F1k,F2k,F3k,F4k] (17)
其中,F1k,F2k,F3k,F4k分别表示样本集中碰撞时刻四个可机动单元拉力突变峰值,k代表样本序号,k=1,2,...P;
神经网络参数为:
A=[F11,F21,F31,F41] (18)
Hinj=Fp*W1+b1 (19)
Houtj=Hinjg(x) (20)
X=Houtj*W2+b2 (21)
Y=Houtj*W2+b2
其中,A为抽取的样本矩阵,Hinj为隐藏层输入矩阵,Houtj为隐藏层输出矩阵,W1表示输入层到隐藏层权值,b1表示输入层到隐藏层参数,W2表示隐藏层到输出层权值,b2表示隐藏层到输出层参数,为激活函数,X,Y表示样本输出层,即碰撞中心点的位置坐标;
误差反向传播调整各层权值:
其中,E为误差,xt为样本期望坐标,Xt为实际输出坐标;
采用梯度下降法使误差趋近极小值:
训练样本集权值更新:
其中,Wji表示第i层到第j层权值输出权值,△Wji表示输出层权值调整量;
步骤5:空间绳系飞网与刚性目标碰撞时,向训练完成的BP神经网络输入步骤3计算得到的可机动单元系绳拉力突变峰值,BP神经网络通过迭代计算,输出碰撞中心点位置坐标(X,Y)。
优选地,所述Q=10,S=20,P=121。
本发明的有益效果如下:
1、本发明设计了一种针对柔性飞网与刚性目标碰撞系绳拉力突变的干扰观测器,通过此观测器计算出突变时刻拉力,以确定BP神经网络的输入层。
2、本发明提出了一种基于BP神经网络的飞网碰撞模中心估计算法,不依靠非合作目标相对速度、质量等信息,可以快速直接通过输入可机动单元系绳拉力精确估计碰撞中心点。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明的目的是提出一种绳系飞网捕获非合作目标碰撞中心点的估计方法,设计可机动单元系绳拉力观测器,计算碰撞时刻可机动单元拉力突变峰值,训练一种BP神经网络,使用该网络可直接输出碰撞中心点坐标。
如图1所示,一种绳系飞网捕获非合作目标的碰撞中心估计方法,包括以下步骤:
步骤1:建立空间绳系飞网机器人系统动力学模型;
空间绳系飞网机器人系统包括四个可机动单元;
空间绳系飞网机器人系统采用弹簧-质杆动力学模型,空间绳系飞网的每个系绳节点的质量表示为:
其中,i表示飞网中的节点,mCM表示可机动单元质量,ρ是飞网绳段材料的密度,A是飞网绳段的横截面积,l是飞网中未伸长的网格边长;
飞网相邻质点间的系绳拉力矢量表示为:
其中,是从第i个节点到第j个节点的单位向量,ri是飞网相邻质点间实际长度;Fi=(k(l-l0)+βkv)表示飞网相邻质点间的系绳拉力大小,l0是每个飞网网格绳段的实际长度,k是飞网绳段的弹性系数,k=EA/l0,E是飞网绳段材料的杨氏模量,β是飞网绳段的阻尼系数,ξ是飞网绳段材料的阻尼比,ν是网格末节点所在轴向相对速度的映射;
根据牛顿第二定律得到飞网中每个质点的动力学方程为:
miai=∑Fi+Fext (3)
其中,ai是质点运动加速度,∑Fi是每个质点系绳拉力之和,Fext表示重力和碰撞力之和;
四个可机动单元的动力学模型为:
mCMaCM=Fj1+Fj2+Fext (4)
其中,Fj1和Fj2分别表示与可机动单元相连接的两段绳子的拉力;
步骤2:设计可机动单元系绳拉力观测器;
可机动单元非线性系统定义为:
设z为系统扰动估计误差,则
进一步得到四个可机动单元构成的非线性系统:
x2=e2+αq (9)
设计扰动估计自适应律w及额外补偿项β(x1,x2),使流形保持不变,对系统扰动估计误差z求导得:
定义扰动估计自适应律w为:
将式(12)带入式(11)中得:
步骤3:计算碰撞时刻可机动单元拉力突变峰值;
根据步骤2中观测得出的可机动单元系绳突变拉力,得到干扰力与拉力峰值分离表达式为:
Ff=f(xtd)=Fd+F (14)
Fd=f(x0) (15)
F=Ff-Fd (16)
其中,Ff表示可机动单元突变时刻峰值,f(xtp)表示拉力观测器,Fd表示耦合干扰力,F表示突变时刻拉力,f(x0)表示非碰撞时刻系绳拉力观测器;
步骤4:基于BP神经网络飞网碰撞训练过程如下:
以10*10m正方形飞网为例,网格边长为20cm,将每个节点进行碰撞动力学仿真实验,共取121个样本点。
训练样本集Fp为:
Fp=[F1k,F2k,F3k,F4k] (17)其中,F1k,F2k,F3k,F4k分别表示样本集中碰撞时刻四个可机动单元拉力突变峰值,k代表样本序号,k=1,2,...121;
神经网络参数为:
A=[F11,F21,F31,F41] (18)
Hinj=Fp*W1+b1 (19)
Houtj=Hinjg(x) (20)
X=Houtj*W2+b2 (21)
Y=Houtj*W2+b2
其中,A为抽取的样本矩阵,Hinj为隐藏层输入矩阵,Houtj为隐藏层输出矩阵,W1表示输入层到隐藏层权值,b1表示输入层到隐藏层参数,W2表示隐藏层到输出层权值,b2表示隐藏层到输出层参数,为激活函数,X,Y表示样本输出层,即碰撞中心点的位置坐标;
误差反向传播调整各层权值:
其中,E为误差,xt为样本期望坐标,Xt为实际输出坐标;
采用梯度下降法使误差趋近极小值:
训练样本集权值更新:
其中,Wji表示第i层到第j层权值输出权值,△Wji表示输出层权值调整量;
步骤5:空间绳系飞网与刚性目标碰撞时,向训练完成的BP神经网络输入步骤3计算得到的可机动单元系绳拉力突变峰值,BP神经网络通过迭代计算,输出碰撞中心点位置坐标(X,Y)。
Claims (2)
1.一种绳系飞网捕获非合作目标的碰撞中心估计方法,其特征在于,包括以下步骤:
步骤1:建立空间绳系飞网机器人系统动力学模型;
空间绳系飞网机器人系统包括四个可机动单元;
空间绳系飞网机器人系统采用弹簧-质杆动力学模型,空间绳系飞网的每个系绳节点的质量表示为:
其中,i表示飞网中的节点,mCM表示可机动单元质量,ρ是飞网绳段材料的密度,A是飞网绳段的横截面积,l是飞网中未伸长的网格边长;
飞网相邻质点间的系绳拉力矢量表示为:
其中,是从第i个节点到第j个节点的单位向量,ri是飞网相邻质点间实际长度;Fi=(k(l-l0)+βkv)表示飞网相邻质点间的系绳拉力大小,l0是每个飞网网格绳段的实际长度,k是飞网绳段的弹性系数,k=EA/l0,E是飞网绳段材料的杨氏模量,β是飞网绳段的阻尼系数,ξ是飞网绳段材料的阻尼比,ν是网格末节点所在轴向相对速度的映射;
根据牛顿第二定律得到飞网中每个质点的动力学方程为:
miai=∑Fi+Fext (3)
其中,ai是质点运动加速度,∑Fi是每个质点系绳拉力之和,Fext表示重力和碰撞力之和;
四个可机动单元的动力学模型为:
mCMaCM=Fj1+Fj2+Fext (4)
其中,Fj1和Fj2分别表示与可机动单元相连接的两段绳子的拉力;
步骤2:设计可机动单元系绳拉力观测器;
可机动单元非线性系统定义为:
设z为系统扰动估计误差,则
进一步得到四个可机动单元构成的非线性系统:
x2=e2+αq (9)
设计扰动估计自适应律w及额外补偿项β(x1,x2),使流形保持不变,对系统扰动估计误差z求导得:
定义扰动估计自适应律w为:
将式(12)带入式(11)中得:
步骤3:计算碰撞时刻可机动单元拉力突变峰值;
根据步骤2中观测得出的可机动单元系绳突变拉力,得到干扰力与拉力峰值分离表达式为:
Ff=f(xtd)=Fd+F (14)
Fd=f(x0) (15)
F=Ff-Fd (16)
其中,Ff表示可机动单元突变时刻峰值,f(xtp)表示拉力观测器,Fd表示耦合干扰力,F表示突变时刻拉力,f(x0)表示非碰撞时刻系绳拉力观测器;
步骤4:基于BP神经网络飞网碰撞训练过程如下:
假定空间绳系飞网为正方形,大小为Q*Q,飞网网格边长为S;在飞网中取P个样本点;
训练样本集Fp为:
Fp=[F1k,F2k,F3k,F4k] (17)其中,F1k,F2k,F3k,F4k分别表示样本集中碰撞时刻四个可机动单元拉力突变峰值,k代表样本序号,k=1,2,...P;
神经网络参数为:
A=[F11,F21,F31,F41] (18)
Hinj=Fp*W1+b1 (19)
Houtj=Hinjg(x) (20)
X=Houtj*W2+b2 (21)
Y=Houtj*W2+b2
其中,A为抽取的样本矩阵,Hinj为隐藏层输入矩阵,Houtj为隐藏层输出矩阵,W1表示输入层到隐藏层权值,b1表示输入层到隐藏层参数,W2表示隐藏层到输出层权值,b2表示隐藏层到输出层参数,为激活函数,X,Y表示样本输出层,即碰撞中心点的位置坐标;
误差反向传播调整各层权值:
其中,E为误差,xt为样本期望坐标,Xt为实际输出坐标;
采用梯度下降法使误差趋近极小值:
训练样本集权值更新:
其中,Wji表示第i层到第j层权值输出权值,△Wji表示输出层权值调整量;
步骤5:空间绳系飞网与刚性目标碰撞时,向训练完成的BP神经网络输入步骤3计算得到的可机动单元系绳拉力突变峰值,BP神经网络通过迭代计算,输出碰撞中心点位置坐标(X,Y)。
2.根据权利要求1所述的一种绳系飞网捕获非合作目标的碰撞中心估计方法,其特征在于,所述Q=10,S=20,P=121。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110235471.0A CN112966333B (zh) | 2021-03-03 | 2021-03-03 | 一种绳系飞网捕获非合作目标的碰撞中心估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110235471.0A CN112966333B (zh) | 2021-03-03 | 2021-03-03 | 一种绳系飞网捕获非合作目标的碰撞中心估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112966333A true CN112966333A (zh) | 2021-06-15 |
CN112966333B CN112966333B (zh) | 2022-09-13 |
Family
ID=76276324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110235471.0A Active CN112966333B (zh) | 2021-03-03 | 2021-03-03 | 一种绳系飞网捕获非合作目标的碰撞中心估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112966333B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113703468A (zh) * | 2021-08-06 | 2021-11-26 | 浙江大学 | 空间绳系机器人位姿一体化控制执行机构 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108363837A (zh) * | 2018-01-18 | 2018-08-03 | 西北工业大学 | 一种绳系飞网系统碰撞力的计算方法 |
CN109063307A (zh) * | 2018-07-26 | 2018-12-21 | 西北工业大学 | 一种可机动飞网抓捕目标的碰撞位移和碰撞力计算方法 |
-
2021
- 2021-03-03 CN CN202110235471.0A patent/CN112966333B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108363837A (zh) * | 2018-01-18 | 2018-08-03 | 西北工业大学 | 一种绳系飞网系统碰撞力的计算方法 |
CN109063307A (zh) * | 2018-07-26 | 2018-12-21 | 西北工业大学 | 一种可机动飞网抓捕目标的碰撞位移和碰撞力计算方法 |
Non-Patent Citations (1)
Title |
---|
刘彤等: "空间绳系飞网运动建模与仿真", 《计算机仿真》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113703468A (zh) * | 2021-08-06 | 2021-11-26 | 浙江大学 | 空间绳系机器人位姿一体化控制执行机构 |
CN113703468B (zh) * | 2021-08-06 | 2023-11-14 | 浙江大学 | 空间绳系机器人位姿一体化控制执行机构 |
Also Published As
Publication number | Publication date |
---|---|
CN112966333B (zh) | 2022-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109063307B (zh) | 一种可机动飞网抓捕目标的碰撞位移和碰撞力计算方法 | |
CN106446442B (zh) | 一种火星伞舱组合体的降落伞展开过程稳定性评估方法 | |
CN103955611B (zh) | 一种空间绳系机器人的通用动力学模型的建立方法 | |
CN104777844B (zh) | 一种高超声速临近空间飞行器航迹跟踪方法 | |
CN114077258B (zh) | 一种基于强化学习ppo2算法的无人艇位姿控制方法 | |
CN109828602B (zh) | 一种基于观测补偿技术的航迹回路非线性模型变换方法 | |
CN112966333B (zh) | 一种绳系飞网捕获非合作目标的碰撞中心估计方法 | |
CN108427322A (zh) | 一种大柔性飞行器基于在线辨识的建模方法 | |
CN110069821B (zh) | 一种空间飞网抛射展开参数确定方法 | |
CN112683261B (zh) | 一种基于速度预测的无人机鲁棒性导航方法 | |
CN111737897A (zh) | 一种深海网箱高密度养殖鱼群的数值模拟方法 | |
CN112668252A (zh) | 一种系留气球动态响应特性计算方法 | |
CN108363837B (zh) | 一种绳系飞网系统碰撞力的计算方法 | |
CN113607378B (zh) | 一种绳系支撑飞行器模型强迫自由角运动模拟与抑制方法 | |
CN103984812A (zh) | 一种适用于平流层卫星系统的空间绳系的仿真方法 | |
CN104715133B (zh) | 一种待辨识对象的运动学参数在轨辨识方法和装置 | |
Liu et al. | Continuous dynamic simulation for morphing wing aeroelasticity | |
CN116933381A (zh) | 一种柔性绳索飞行器系统动力学仿真分析方法 | |
Azouz et al. | Dynamic analysis of airships with small deformations | |
Chin et al. | A numerical model of a towed cable-body system | |
CN112380692A (zh) | 一种运载火箭的大气层内在线轨迹规划方法 | |
Hembree et al. | Efficient tether dynamic model formulation using recursive rigid-body dynamics | |
CN115016291B (zh) | 一种基于风场信息的飞行器抗干扰姿态控制系统和方法 | |
CN118484652B (zh) | 基于elm与ukf的飞行器气动参数在线辨识方法 | |
CN113110571A (zh) | 一种基于降维状态观测器估计飞行攻角的方法 |
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 |