CN107544262A - 一种运载火箭自适应精确回收控制方法 - Google Patents
一种运载火箭自适应精确回收控制方法 Download PDFInfo
- Publication number
- CN107544262A CN107544262A CN201711027681.0A CN201711027681A CN107544262A CN 107544262 A CN107544262 A CN 107544262A CN 201711027681 A CN201711027681 A CN 201711027681A CN 107544262 A CN107544262 A CN 107544262A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- mfrac
- mtd
- 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 19
- 238000011084 recovery Methods 0.000 title claims abstract description 12
- 230000003044 adaptive effect Effects 0.000 claims abstract description 18
- 230000008859 change Effects 0.000 claims description 13
- 238000013461 design Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000004064 recycling Methods 0.000 claims description 8
- 238000009415 formwork Methods 0.000 claims description 4
- 230000035935 pregnancy Effects 0.000 claims description 4
- 230000009897 systematic effect Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000003094 perturbing effect Effects 0.000 claims description 3
- 230000001141 propulsive effect Effects 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 2
- 239000000446 fuel Substances 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 5
- 239000007787 solid Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 239000003380 propellant Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开一种运载火箭自适应精确回收控制方法,采用真实体现其动力学特性的非线性模型,建立返回段动力学模型,并引入系统参数不确定性和干扰;为估计重复使用运载火箭动力学模型中的参数不确定和提高估计的准确性,建立自适应滑模观测器以及设计参数Γ和的自适应估计律;基于反演控制思想并设计滑模动态面跟踪控制器。本发明同时考虑了气动参数的未建模动态及燃料消耗影响引起的不确定性因素,提出模型更符合垂直返回段运载火箭一子级控制系统的实际运行情况;设计的自适应参数观测器不需要不确定和干扰上界这一已知条件,在实际应用中更具有优势;不确定估计的准确性因为时变增益在线相对于干扰适当的调节而得到了改善。
Description
技术领域
本发明涉及航空航天飞行控制领域,具体涉及可重复使用运载火箭一子级垂直着陆返回的姿态与位置的精确控制问题。
背景技术
自人类开展航天活动以来,低成本、高可靠、高效率的航天发射系统一直是航天界追求的目标。可重复使用运载火箭完成预定发射任务后,全部或部分返回并安全着陆,经过检修维护与燃料加注,可再次执行发射任务,通过多次使用分摊费用来降低运载火箭的生产与发射成本,因此各航天大国都将其作为未来发展重点。而垂直返回段的安全着陆控制是重复使用运载火箭的关键技术。目前,重复使用运载火箭多采用两级入轨部分重复使用的方式,即仅对火箭一子级进行回收使用。2015年中国运载火箭技术研究院公开的专利CN201510220610.7中提到了将N台固体游动发动机沿周向均匀地安装在运载火箭的尾部,并分别与运载火箭的N台舵机相连(N≥4),采集惯性测量组合实时测量的信息,经制导计算后直接生成俯仰、偏航、滚动通道姿态控制量并驱动与所述N台舵机相连的固体游动发动机侧摆工作,最终各台游动发动机侧摆产生的侧向力合成产生箭体姿态控制力的方法,由于该方法仅限于安装有固体游动发动机的固体运载火箭,所以,不适应于一般的没有安装多台游动发动机运载火箭的姿态控制;专利CN201080036159.9提出了有关可重复使用太空运载火箭在海洋中或其它水体中的平台上着陆和重新恢复助推级和/或其它其中部分的系统和方法,实现太空运载火箭助推级垂直动力降落在预先定位的航海、平台的甲板上,它利用双向气动控制面控制助推级的轨迹,利用航海平台向助推级实时广播其位置信息,来弥补由于洋流漂移和/或其他因素带来的航海平台的位置误差,从而保证着陆时的位置精度,但是该发明没有考虑安全垂直回收着陆时的另一个关键技术——姿态控制。目前,其他公开资料及文献也未见有这方面的详细报道。
由于运载火箭一子级返回飞行的空域跨度变化大、动压变化剧烈、气动参数偏差和干扰严重,各通道呈现严重的非线性耦合特性,这使得垂直回收时姿态和位置的精确控制设计难度较大,特别是当系统内部参数改变或发生严重外部扰动时,经典控制理论很难满足现代运载火箭一子级特殊机动下的高性能控制需求。因此,必须采用现代控制理论进行可重复使用运载火箭一子级垂直回收控制方法研究。近年来,国内外许多学者针对这一问题,分别都提出了自己的非线性系统设计方法,但大多存在以下问题之一:
(1)模型较为理想化,不符合运载火箭一子级系统实际的运行情况;
(2)适用条件受限较大,工程应用难以实现;
(3)系统不确定性和干扰等对系统的控制性能影响较大;
(4)不能满足全局范围内的稳定性、鲁棒性和性能等指标要求。
发明内容
鉴于上述现有技术中的不足之处,本发明提出基于自适应滑模动态面技术的重复使用运载火箭一子级精确垂直回收控制方法,通过以下步骤实现:
步骤1、建立火箭返回段动力学模型,
其中V为飞行速度,m为火箭一子级质量,P为发动机总推力,εn为返回过程中第n次(n=1、2或3)有动力飞行段的变推力因子且0≤εn≤1,α、β分别为攻角和侧滑角,X、Y、Z为速度坐标系下的气动力分量,g为重力加速度,Rx、Ry、Rz为发射坐标系下的位置分量,r为火箭一子级质心距地心距离,θ、σ为弹道倾角和弹道偏角,为地球半径,wx、wy、wz和Vx、Vy、Vz分别为发射坐标系下的地球自转角速度和速度分量;
步骤2、将步骤1中的动力学模型转化为一般的状态空间形式,并引入系统参数不确定性和干扰,将所述非线性模型变换为如下形式,
其中,x1=[Rx,Ry,Rz]T、x2=[V,θ,σ]T为状态变量,u=[X,Y,Z]T为控制变量,Δf为系统关于大气密度未建模动态以及火箭一子级本体质量及质心变化引起的摄动项;ΔB2为一子级本体质量及质心变化引起的不确定项,d为未知干扰项;其中
f(x1,x2)=[f1,f2,f3]T,
步骤3、建立自适应滑模观测器以及设计参数Γ和的自适应估计律,具体如下:
这里Σ是一个事先确定的对角矩阵,Γ(0)>0,这里未知正数ξ为Δf+ΔB2u的上界,即|Δf+ΔB2u|≤ξ,时变增益K(t)定义见步骤4。
上述正定矩阵Σ和正实数μ可以通过如下不等式进行求解:
步骤4、利用所述滑模观测器和所述自适应估计律得到的状态和不确定项的估计值,抵消不确定性与未知干扰,设计滑模动态面跟踪控制器,具体如下:
选取如下虚拟控制律:
让x2d通过时间常数为τ2的一阶滤波器,得实际控制器为
其中δ为x2d通过一阶滤波器后得到输出量,即新的替代x2d的状态量;
步骤5、验证运载火箭一子级闭环控制系统的稳定性。
本发明同时考虑了气动参数的未建模动态及燃料消耗影响引起的不确定性因素,建立的模型更符合垂直返回段运载火箭一子级控制系统的实际运行情况;设计的自适应参数观测器不需要不确定和干扰上界这一已知条件,在实际应用中更具有优势;不确定估计的准确性因为时变增益在线相对于干扰适当的调节而得到了改善;重复使用运载火箭一子级垂直回收阶段,设计的控制律实现精确垂直着陆回收的控制目标。
附图说明
图1为本发明运载火箭自适应精确回收控制方法一实施例的框图;
图2为图1实施例弹道偏角跟踪轨迹仿真对比图;
图3为图1实施例弹道倾角跟踪轨迹仿真对比图;
图4为图1实施例x轴方向位置跟踪轨迹仿真对比图;
图5为图1实施例y轴方向位置跟踪轨迹仿真对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
为了本领域普通技术人员可以更好地理解本发明的实施,本实施例提供了利用Matlab2012b软件进行自适应滑模动态面跟踪控制与标准动态面跟踪控制的性能对比的仿真验证结果。
运载火箭自适应精确回收控制方法,包括以下步骤:
步骤一、充分考虑重复使用运载火箭的动压变化剧烈、非线性耦合等工作特点,采用真实体现其动力学特性的非线性模型,建立返回段动力学模型具体如下:
其中V为飞行速度,m为火箭一子级质量,P位发动机总推力,εn为返回过程中第n次有动力飞行段的变推力因子且0≤εn≤1,α、β分别为攻角和侧滑角,X、Y、Z为速度坐标系下的气动力分量,g为重力加速度,Rx、Ry、Rz为发射坐标系下的位置分量,r为火箭一子级质心距地心距离,θ、σ为弹道倾角和弹道偏角,为地球半径,wx、wy、wz和Vx、Vy、Vz分别为发射坐标系下的地球自转角速度和速度分量。
步骤二、为正确描述重复使用运载火箭气动参数不确定和干扰严重的情况,将返回段的动力学模型转化为一般的状态空间形式,并引入系统参数不确定性和干扰,将步骤一的非线性模型进一步修正为如下形式:
其中,x1=[Rx,Ry,Rz]T、x2=[V,θ,σ]T为状态变量,u=[X,Y,Z]T为控制变量;其中Δf为系统关于大气密度未建模动态以及火箭一子级本体质量及质心变化引起的摄动项;ΔB2为一子级本体质量及质心变化引起的不确定项,d为未知干扰项。
f(x1,x2)=[f1,f2,f3]T,
步骤三、为估计重复使用运载火箭动力学模型中的参数不确定和提高估计的准确性,建立自适应滑模观测器以及设计参数Γ和的自适应估计律,具体如下:
这里Σ是一个事先确定的对角矩阵,Γ(0)>0,这里未知正数ξ为Δf+ΔB2u的上界,即有|Δf+ΔB2u|≤ξ,时变增益K(t)定义见步骤四。
上述正定矩阵Σ和正实数μ可以通过如下不等式进行求解:
步骤四、为验证步骤三所设计的滑模观测器以及设计自适应参数估计律的正确性,需利用Lyapunov稳定性理论证明对系统状态和不确定性的估计误差是渐近稳定的,包含以下步骤:
(a)定义Lyapunov函数:
时变增益K(t)定义如下:
(i)如果有
γ1>且K(0)>0
(ii)如果有
K(t)=γ2ηTsign(η)+γ3,γ3>0
这里η是符号函数经过低通滤波后的平均值,γ3>0,γ2=K(t*)。t*是最大的时间值,即是且
(b)根据状态估计误差结果不同分两种情况分析:
(i)当时为第一种情况,此时有
(ii)当即时为第二种情况,此时有
综上,当正定矩阵∑和正实数μ满足如下不等式时,动态误差系统是渐近稳定的。
步骤五、基于反演控制思想,利用前面的滑模观测器和自适应参数估计器得到的状态和不确定项的估计值,抵消不确定性与未知干扰,并设计滑模动态面跟踪控制器,具体如下:
选取如下虚拟控制律:
让x2d通过时间常数为τ2的一阶滤波器,得实际控制器为
其中δ为x2d通过一阶滤波器后得到输出量,即新的替代x2d的状态量。
步骤六、利用Lyapunov稳定性理论,验证火箭一子级闭环控制系统对垂直回收控制指令的跟踪误差是最终一致有界的,包含以下步骤:
(a)定义Lyapunov函数:V=V2+V3
这里φ为实际闭环系统的状态量δ与外环虚拟控制量x2d之间的误差,即φ=δ-x2d。
(b)对Lyapunov函数进行求导,可得不等式:
其中
(c)选取适当的实数k1、k2和k3使其满足下列条件
此时证明了重复使用运载火箭一子级闭环跟踪控制系统对垂直回收轨迹/位置、姿态速度控制指令的跟踪误差最终一致有界,且可收敛于任意小的一个区域内。
本发明利用Matlab2012b软件进行自适应滑模动态面跟踪控制与标准动态面跟踪控制的性能对比的仿真验证:
(1)火箭一子级回收段控制系统参数选取:取k1=1,k2=1.5,μ=3,τ1=1,τ2=0.1,
γ1=2,矩阵Σ=diag(2,3,4)。
(2)初始参数选取:初始位置x1(0)为[9,6,200],单位米;初始速度166米/秒,弹道倾角初始角度θ0=85°,弹道偏角初始角度为σ0=2°。
(3)给定指令参数:着陆回收速度指令Vc=0.2米/秒,倾角回收指令角度θC=90°,弹道偏角指令为σC=0°。
(4)燃料/推进剂消耗参数选取:给定O=1.3吨/秒和O=3吨/秒两种情况。
结果说明:
如图2和图3所示,随着推进剂消耗量的增大,使用标准动态面控制方法,系统状态量与指令之间的偏差随之增大,如弹道倾角θ和弹道偏角σ明显偏离θc和σc,而滑模动态面跟踪控制方法则不会受此影响,均能准确稳定跟踪姿态控制指令,验证了姿态动态跟踪误差系统是渐近稳定的,精确垂直回收姿态控制得以实现。
由图4和图5所示,标准动态面控制方法对运载火箭一子级垂直降落阶段的落点位置控制的精确度不够,与给定指令之间的偏差较大,而本文所提方法即使在存在燃料不断消耗引起不确定性的情况下,仍保证垂直回收段火箭一子级对轨迹/位置控制指令的跟踪误差最终一致有界,且可收敛于任意小的一个区域内,实现了精确定点安全回收的控制目标。通过对比,不难发现本发明所设计的非线性系统控制方法的有效性。
本发明未详细说明部分都属于领域技术人员公知常识,以上所述仅为本发明的一个具体实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种运载火箭自适应精确回收控制方法,其特征在于,包括以下步骤:
步骤1、建立火箭返回段动力学模型,
<mrow>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mi>n</mi>
</msub>
<mi>P</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&alpha;</mi>
<mo>-</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>R</mi>
<mi>x</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>R</mi>
<mi>y</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>&theta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mi>n</mi>
</msub>
<mi>P</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>Y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>R</mi>
<mi>x</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>R</mi>
<mi>y</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>w</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mover>
<mi>&sigma;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mi>n</mi>
</msub>
<mi>P</mi>
<mi>&beta;</mi>
<mo>-</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>R</mi>
<mi>x</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>g</mi>
<mfrac>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>R</mi>
<mi>y</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>w</mi>
<mi>x</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>y</mi>
</msub>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mover>
<mi>R</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>x</mi>
<mo>=</mo>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&sigma;</mi>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mover>
<mi>R</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>y</mi>
<mo>=</mo>
<mi>V</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&sigma;</mi>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<mover>
<mi>R</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>z</mi>
<mo>=</mo>
<mo>-</mo>
<mi>V</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&sigma;</mi>
</mrow>
其中V为飞行速度,m为火箭一子级质量,P为发动机总推力,εn为返回过程中第n次(n=1、2或3)有动力飞行段的变推力因子且0≤εn≤1,α、β分别为攻角和侧滑角,X、Y、Z为速度坐标系下的气动力分量,g为重力加速度,Rx、Ry、Rz为发射坐标系下的位置分量,r为火箭一子级质心距地心距离,θ、σ为弹道倾角和弹道偏角,为地球半径,wx、wy、wz和Vx、Vy、Vz分别为发射坐标系下的地球自转角速度和速度分量;
步骤2、将步骤1中的动力学模型转化为一般的状态空间形式,并引入系统参数不确定性和干扰,将所述非线性模型变换为如下形式,
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>f</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mi>u</mi>
<mo>+</mo>
<msub>
<mi>&Delta;B</mi>
<mn>2</mn>
</msub>
<mi>u</mi>
<mo>+</mo>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,x1=[Rx,Ry,Rz]T、x2=[V,θ,σ]T为状态变量,u=[X,Y,Z]T为控制变量,Δf为运载火箭控制系统关于大气密度未建模动态以及火箭一子级本体质量及质心变化引起的摄动项;ΔB2为一子级本体质量及质心变化引起的不确定项,d为未知干扰项;其中
<mrow>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&sigma;</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&sigma;</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mi>&sigma;</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mi>m</mi>
<mi>V</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mi>m</mi>
<mi>V</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
f(x1,x2)=[f1,f2,f3]T,
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>n</mi>
</msub>
<mi>m</mi>
</mfrac>
<mi>P</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&alpha;</mi>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>x</mi>
<mi>r</mi>
</mfrac>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>y</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>n</mi>
</msub>
<mi>m</mi>
</mfrac>
<mi>P</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>x</mi>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>y</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>w</mi>
<mi>z</mi>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>n</mi>
</msub>
<mi>m</mi>
</mfrac>
<mi>P</mi>
<mi>&beta;</mi>
<mo>+</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>x</mi>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>g</mi>
<mfrac>
<mrow>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>y</mi>
</mrow>
<mi>r</mi>
</mfrac>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>w</mi>
<mi>x</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>y</mi>
</msub>
<mi>cos</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3、建立自适应滑模观测器以及设计参数Γ和的自适应估计律,具体如下:
<mrow>
<msub>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>&Sigma;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mi>u</mi>
<mo>+</mo>
<mover>
<mi>&xi;</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mi>&Gamma;</mi>
<mover>
<mover>
<mi>&xi;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>&Gamma;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mi>&mu;</mi>
<mi>&Gamma;</mi>
<mo>+</mo>
<msub>
<mi>I</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<mover>
<mover>
<mi>&xi;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<msup>
<mi>&Gamma;</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
其中Σ是一个事先确定的对角矩阵,Γ(0)>0,这里未知正数ξ为Δf+ΔB2u的上界,即|Δf+ΔB2u|≤ξ,时变增益K(t)定义如下:
(i)如果有
且K(0)>0
(ii)如果有
K(t)=γ2ηTsign(η)+γ3,γ3>0
<mrow>
<mi>&pi;</mi>
<mover>
<mi>&eta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mi>&eta;</mi>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>&pi;</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
其中η是符号函数经过低通滤波后的平均值,γ3≥0,γ2=K(t*).t*是最大的时间值,即且
上述正定矩阵Σ和正实数μ可以通过如下不等式进行求解:
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>min</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>&Sigma;</mo>
<mo>)</mo>
</mrow>
<mo>></mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
</mrow>
<mrow>
<msup>
<mi>&Gamma;</mi>
<mi>T</mi>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&Gamma;</mi>
<mi>T</mi>
</msup>
<mi>&Gamma;</mi>
<mo><</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>I</mi>
<mn>3</mn>
</msub>
<mo>;</mo>
</mrow>
步骤4、利用所述滑模观测器和所述自适应估计律得到的状态和不确定项的估计值,抵消不确定性与未知干扰,设计滑模动态面跟踪控制器,具体如下:
选取如下虚拟控制律:
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mn>1</mn>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
</mrow>
让x2d通过时间常数为τ2的一阶滤波器,得实际控制器为
<mrow>
<mi>u</mi>
<mo>=</mo>
<msubsup>
<mi>B</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mover>
<mi>&delta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mi>f</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mover>
<mi>&xi;</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msubsup>
<mi>B</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中δ为x2d通过一阶滤波器后得到输出量,即新的替代x2d的状态量;
步骤5、验证运载火箭一子级闭环控制系统的稳定性。
2.根据权利要求1所述的运载火箭自适应精确回收控制方法,其特征在于,还包括步骤3'、用来验证步骤3所述的滑模观测器以及所述自适应估计律的正确性,证明对系统状态和不确定性的估计误差是渐近稳定的,具体是:
给定李雅普诺夫函数:
<mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mover>
<mi>&xi;</mi>
<mo>~</mo>
</mover>
<mi>T</mi>
</msup>
<msup>
<mi>&Gamma;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>&Gamma;</mi>
<mover>
<mi>&xi;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<mover>
<mi>d</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>,</mo>
<mo>|</mo>
<mi>d</mi>
<mo>|</mo>
<mo>&le;</mo>
<mover>
<mi>d</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
其中通过选取合适的李雅普诺夫函数可以证明时变增益K(t)超过未知干扰上界,所以保证估计误差渐近收敛到0;一旦时变增益K(t)便由步骤3中的(ii)式调节,可减少对于时变干扰的抖振;
由于成立,所以所述滑模观测器和所述自适应估计律对系统状态和不确定性的估计误差是渐近稳定的。
3.根据权利要求1所述的运载火箭自适应精确回收控制方法,其特征在于,步骤5的过程具体是:
取李雅普诺夫函数:
V=V2+V3,
其中e1=x1d-x1,e2=x2-δ,引入一个新的状态变量δ,令x2d通过时间常数为τ2的一阶滤波器获得新的变量δ,且δ(0)=x2(0),φ=δ-x2d;推导成立,即闭环跟踪控制系统对位置控制指令x1d以及姿态速度控制指令的跟踪误差最终一致有界,证明该闭环控制系统实现了对重复使用运载火箭一子级垂直回收位置和姿态的精确回收控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711027681.0A CN107544262B (zh) | 2017-10-27 | 2017-10-27 | 一种运载火箭自适应精确回收控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711027681.0A CN107544262B (zh) | 2017-10-27 | 2017-10-27 | 一种运载火箭自适应精确回收控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107544262A true CN107544262A (zh) | 2018-01-05 |
CN107544262B CN107544262B (zh) | 2019-10-11 |
Family
ID=60968347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711027681.0A Active CN107544262B (zh) | 2017-10-27 | 2017-10-27 | 一种运载火箭自适应精确回收控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107544262B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108423196A (zh) * | 2018-04-08 | 2018-08-21 | 中国人民解放军战略支援部队航天工程大学 | 第一级能重复使用的两级入轨航天器的入轨方法 |
CN108828946A (zh) * | 2018-07-06 | 2018-11-16 | 合肥工业大学 | 一种用于火箭整流罩减震降噪的自适应控制方法 |
CN109343341A (zh) * | 2018-11-21 | 2019-02-15 | 北京航天自动控制研究所 | 一种基于深度强化学习的运载火箭垂直回收智能控制方法 |
CN110243362A (zh) * | 2019-06-27 | 2019-09-17 | 西北工业大学 | 一种中高空超声速靶标导航方法 |
CN110989650A (zh) * | 2019-12-30 | 2020-04-10 | 航天科工火箭技术有限公司 | 一种基于火箭返回的着陆控制方法和装置 |
CN111473696A (zh) * | 2020-03-24 | 2020-07-31 | 北京理工大学 | 一种基于落点估计的制导火箭垂直攻击制导方法 |
CN113408158A (zh) * | 2021-08-19 | 2021-09-17 | 中国科学院力学研究所 | 一种适用于运载火箭级间冷分离的实现方法 |
CN113479347A (zh) * | 2021-07-13 | 2021-10-08 | 北京理工大学 | 一种火箭垂直回收着陆段轨迹控制方法 |
CN113504723A (zh) * | 2021-07-05 | 2021-10-15 | 北京航空航天大学 | 一种基于逆强化学习的运载火箭减载控制方法 |
CN113701571A (zh) * | 2021-08-17 | 2021-11-26 | 航天科工火箭技术有限公司 | 一种垂直回收演示验证火箭飞行时序设计方法 |
CN114326408A (zh) * | 2022-01-06 | 2022-04-12 | 广东省智能机器人研究院 | 兼顾位置和姿态调节的移动机器人反演控制方法及装置 |
CN116594309A (zh) * | 2023-07-19 | 2023-08-15 | 东方空间技术(山东)有限公司 | 一种火箭子级垂直回收控制方法、计算设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4504031A (en) * | 1979-11-01 | 1985-03-12 | The Boeing Company | Aerodynamic braking and recovery method for a space vehicle |
US20050022500A1 (en) * | 2003-07-28 | 2005-02-03 | Buehler David Benjamin | Rocket engine passivation system |
CN103941742A (zh) * | 2014-04-29 | 2014-07-23 | 中国科学院自动化研究所 | 一种无人机滑跑纠偏控制装置及方法 |
CN104460319A (zh) * | 2014-10-20 | 2015-03-25 | 南京工业大学 | 多批次间歇反应过程二维多模型预测迭代学习控制方法 |
CN104635495A (zh) * | 2015-02-09 | 2015-05-20 | 北京航空航天大学 | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 |
CN105931130A (zh) * | 2016-04-11 | 2016-09-07 | 南京工业大学 | 一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法 |
CN106021628A (zh) * | 2015-07-03 | 2016-10-12 | 中国运载火箭技术研究院 | 一种运载火箭垂直返回弹道设计方法 |
-
2017
- 2017-10-27 CN CN201711027681.0A patent/CN107544262B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4504031A (en) * | 1979-11-01 | 1985-03-12 | The Boeing Company | Aerodynamic braking and recovery method for a space vehicle |
US20050022500A1 (en) * | 2003-07-28 | 2005-02-03 | Buehler David Benjamin | Rocket engine passivation system |
CN103941742A (zh) * | 2014-04-29 | 2014-07-23 | 中国科学院自动化研究所 | 一种无人机滑跑纠偏控制装置及方法 |
CN104460319A (zh) * | 2014-10-20 | 2015-03-25 | 南京工业大学 | 多批次间歇反应过程二维多模型预测迭代学习控制方法 |
CN104635495A (zh) * | 2015-02-09 | 2015-05-20 | 北京航空航天大学 | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 |
CN106021628A (zh) * | 2015-07-03 | 2016-10-12 | 中国运载火箭技术研究院 | 一种运载火箭垂直返回弹道设计方法 |
CN105931130A (zh) * | 2016-04-11 | 2016-09-07 | 南京工业大学 | 一种考虑量测信号丢失的改进集合卡尔曼滤波估计方法 |
Non-Patent Citations (5)
Title |
---|
NEBYLOV A,NEBYLOV V: "Reusable space planes challenges and control problems", 《INTERNATIONAL FEDERATION OF AUTOMATIC CONTROL》 * |
ZHIFENG GAO,等: "Active fault tolerant control design for reusable launch vehicle using adaptive sliding mode technique", 《JOURNAL OF THE FRANKLIN INSTITUTE》 * |
ZHIFENG GAO,等: "Adaptive Dynamic Surface Tracking Scheme Design for Rossler Chaotic Systems with Control Input Saturation", 《2015 27TH CHINESE CONTROL AND DECISION CONFERENCE》 * |
汪小卫,等: "运载火箭子级回收技术研究", 《航天返回与遥感》 * |
钱默抒: "用于舰载无人机回收的惯导测高系统研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108423196A (zh) * | 2018-04-08 | 2018-08-21 | 中国人民解放军战略支援部队航天工程大学 | 第一级能重复使用的两级入轨航天器的入轨方法 |
CN108828946A (zh) * | 2018-07-06 | 2018-11-16 | 合肥工业大学 | 一种用于火箭整流罩减震降噪的自适应控制方法 |
CN108828946B (zh) * | 2018-07-06 | 2021-04-27 | 合肥工业大学 | 一种用于火箭整流罩减震降噪的自适应控制方法 |
CN109343341A (zh) * | 2018-11-21 | 2019-02-15 | 北京航天自动控制研究所 | 一种基于深度强化学习的运载火箭垂直回收智能控制方法 |
CN109343341B (zh) * | 2018-11-21 | 2021-10-01 | 北京航天自动控制研究所 | 一种基于深度强化学习的运载火箭垂直回收智能控制方法 |
CN110243362A (zh) * | 2019-06-27 | 2019-09-17 | 西北工业大学 | 一种中高空超声速靶标导航方法 |
CN110989650B (zh) * | 2019-12-30 | 2023-10-27 | 航天科工火箭技术有限公司 | 一种基于火箭返回的着陆控制方法和装置 |
CN110989650A (zh) * | 2019-12-30 | 2020-04-10 | 航天科工火箭技术有限公司 | 一种基于火箭返回的着陆控制方法和装置 |
CN111473696A (zh) * | 2020-03-24 | 2020-07-31 | 北京理工大学 | 一种基于落点估计的制导火箭垂直攻击制导方法 |
CN111473696B (zh) * | 2020-03-24 | 2021-01-26 | 北京理工大学 | 一种基于落点估计的制导火箭垂直攻击制导方法 |
CN113504723B (zh) * | 2021-07-05 | 2023-11-28 | 北京航空航天大学 | 一种基于逆强化学习的运载火箭减载控制方法 |
CN113504723A (zh) * | 2021-07-05 | 2021-10-15 | 北京航空航天大学 | 一种基于逆强化学习的运载火箭减载控制方法 |
CN113479347A (zh) * | 2021-07-13 | 2021-10-08 | 北京理工大学 | 一种火箭垂直回收着陆段轨迹控制方法 |
CN113479347B (zh) * | 2021-07-13 | 2023-12-08 | 北京理工大学 | 一种火箭垂直回收着陆段轨迹控制方法 |
CN113701571A (zh) * | 2021-08-17 | 2021-11-26 | 航天科工火箭技术有限公司 | 一种垂直回收演示验证火箭飞行时序设计方法 |
CN113408158B (zh) * | 2021-08-19 | 2021-11-09 | 中国科学院力学研究所 | 一种适用于运载火箭级间冷分离的实现方法 |
CN113408158A (zh) * | 2021-08-19 | 2021-09-17 | 中国科学院力学研究所 | 一种适用于运载火箭级间冷分离的实现方法 |
CN114326408A (zh) * | 2022-01-06 | 2022-04-12 | 广东省智能机器人研究院 | 兼顾位置和姿态调节的移动机器人反演控制方法及装置 |
CN114326408B (zh) * | 2022-01-06 | 2024-03-19 | 广东省智能机器人研究院 | 兼顾位置和姿态调节的移动机器人反演控制方法及装置 |
CN116594309A (zh) * | 2023-07-19 | 2023-08-15 | 东方空间技术(山东)有限公司 | 一种火箭子级垂直回收控制方法、计算设备及存储介质 |
CN116594309B (zh) * | 2023-07-19 | 2023-09-29 | 东方空间技术(山东)有限公司 | 一种火箭子级垂直回收控制方法、计算设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107544262B (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107544262A (zh) | 一种运载火箭自适应精确回收控制方法 | |
CN107966156B (zh) | 一种适用于运载火箭垂直回收段的制导律设计方法 | |
CN102163059B (zh) | 推力变向无人机的姿态控制系统及控制方法 | |
CN113479347B (zh) | 一种火箭垂直回收着陆段轨迹控制方法 | |
CN109657256B (zh) | 一种高精度弹道式再入标称返回轨道仿真方法 | |
CN104022742A (zh) | 基于神经网络观测器的飞行器姿态鲁棒反演容错控制方法 | |
CN102819266B (zh) | 一种拟周期j2不变相对轨道编队飞行控制方法 | |
CN106371312A (zh) | 基于模糊控制器的升力式再入预测‑校正制导方法 | |
CN102749852A (zh) | 多源干扰系统的容错抗干扰控制方法 | |
CN102880052A (zh) | 基于时标功能分解的高超声速飞行器执行器饱和控制方法 | |
CN105159308A (zh) | 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法 | |
CN109250153A (zh) | 火星大气进入段轨迹最优跟踪制导方法 | |
CN105573337A (zh) | 一种满足再入角和航程约束的离轨制动闭路制导方法 | |
CN105807789A (zh) | 基于t-s模糊观测器补偿的uuv控制方法 | |
CN107132764A (zh) | 基于自抗扰控制技术的碎片抓捕后轨道转移控制方法 | |
CN108516104A (zh) | 多级运载航天器分离可靠性仿真平台 | |
CN103112603A (zh) | 欠驱动高速自旋卫星建立正常姿态的方法 | |
CN109683624A (zh) | 用于小型无人直升机姿态控制的非线性鲁棒控制方法 | |
CN108646554A (zh) | 一种基于指定性能的飞行器快速抗干扰纵向制导方法 | |
CN113268084A (zh) | 一种无人机编队智能容错控制方法 | |
CN105759609A (zh) | 一种基于显式非线性模型预测控制的舰载机自主着舰方法 | |
Benito et al. | Nonlinear predictive controller for drag tracking in entry guidance | |
CN104635495B (zh) | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 | |
Mao et al. | REAL-TIME ONLINE TRAJECTORY PLANNING AND GUIDANCE FOR TERMINAL AREA ENERGY MANAGEMENT OF UNMANNED AERIAL VEHICLE. | |
CN107894775A (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 |