CN102862686A - 再入飞行器的最优积分滑模姿态控制方法及控制器 - Google Patents
再入飞行器的最优积分滑模姿态控制方法及控制器 Download PDFInfo
- Publication number
- CN102862686A CN102862686A CN2012103672887A CN201210367288A CN102862686A CN 102862686 A CN102862686 A CN 102862686A CN 2012103672887 A CN2012103672887 A CN 2012103672887A CN 201210367288 A CN201210367288 A CN 201210367288A CN 102862686 A CN102862686 A CN 102862686A
- Authority
- CN
- China
- Prior art keywords
- omega
- control
- nominal
- module
- attitude
- 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
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种再入飞行器的最优积分滑模姿态控制方法及控制器,属于飞行器控制技术领域。首先针对飞行器的标称模型设计了SDRE标称姿态控制律,使标称***的性能满足提出的最优指标。然后,考虑***的不确定性,在SDRE标称控制律的基础上设计积分滑模控制律,使***在满足性能指标要求的同时,对不确定性具有鲁棒性。为了减弱抖振,引入二阶滑模设计思想,使控制器输出较光滑。本发明设计的姿态控制器不仅能保证期望的指标,而且具有较好的鲁棒性。
Description
技术领域
本发明涉及一种再入飞行器的最优积分滑模姿态控制方法及控制器,属于飞行器控制技术领域。
背景技术
飞行器无动力再入飞行过程中,空气密度和飞行器速度变化较大,动力学参数变化剧烈,通道间的耦合作用非常严重,表现出强烈的多变量耦合和非线性,且伴随着其他未知干扰和不确定,因此,针对快时变、强不确定性的***设计鲁棒姿态控制器非常关键。
目前,针对再入飞行的特点,许多学者提出了不同的姿态控制器,例如鲁棒控制器、自适应控制器、最优控制器、滑模控制器等。基于状态依赖黎卡提方程(SDRE)的方法给出了一种非线性***的最优控制律设计方法,它通过构造仿射非线性***的状态依赖参数(SDC)形式,将***的调节器设计问题转换为“LQR(线性二次型调节器)”问题,且充分保留了***中的非线性特性。这种方法实现简单,可通过调整权值矩阵有效折衷控制量和***的动态性能,并且,经过参数化后的***矩阵以及设计的权值矩阵与***状态相关,***具有很好的设计灵活性。
滑模控制对匹配的参数不确定性和外部扰动具有良好的鲁棒性,且具有快速的动态响应能力,因此考虑将SDRE方法与滑模方法结合,设计最优滑模控制器,在发挥SDRE性能优势的同时,保证***的鲁棒性。将最优控制理论与滑模控制理论相结合的一种思想在于将线性二次型最优控制理论应用于滑模面的设计,该方法被应用于不确定线性***以提高LQR的性能,然而对于非线性***,在进行滑模面优化时,会导致难以求解的两点边值问题;另一种思想是用滑变结构控制理论对最优控制器进行鲁棒化设计,采用积分滑模的概念,控制律的一部分是针对标称线性***确定的LQR控制律,另一部分是积分滑模控制,保证了***鲁棒性,然而,由于使用了符号函数,此类控制律在提高***鲁棒性的同时引入了滑模的抖振。采用边界层可以减弱抖振,但是此时滑模面被限制在一个较小的区域内,而无法证明滑模的可达性,使得***性能偏离最优指标。因此,需要设计控制器,既保证***在存在不确定性时的***性能,又能减弱***抖振。
发明内容
本发明针对飞行器再入段的姿态控制问题,结合积分滑模与状态依赖黎卡提方程(SDRE)方法,设计了一种最优积分滑模(OISMC)姿态控制方法。首先针对飞行器的标称模型设计了SDRE标称姿态控制律,使标称***的性能满足提出的最优指标。然后,考虑***的不确定性,在SDRE标称控制律的基础上设计积分滑模控制律,使***在满足性能指标要求的同时,对不确定性具有鲁棒性。为了减弱抖振,引入二阶滑模设计思想,使控制器输出较光滑。
本发明的最优积分滑模姿态控制方法具体包括以下步骤:
步骤1,生成飞行器的状态向量。
结合飞行器的实际姿态角Ω=[α,β,μ]T,姿态角速度ω=[ωx,ωy,ωz]T,以及速度V,组成状态向量x:x=[V,α,β,μ,ωx,ωy,ωz]T。
步骤2,建立再入飞行器的数学模型:
y=h(x)
其中,状态向量x=[V α β μ ωx ωy ωz]T,控制力矩u=[Mx My Mz]T,输出变量y=[α β μ]T,f(x)=[f1(x) f2(x)...f7(x)]T。
f1(x)=(-X-mgsinγ)/m
f2(x)=ωz+tanβ(ωysinα-ωxcosα)-(Y-mgcosγcosμ)/(mVcosβ)
f3(x)=ωxsinα+ωycosα+(Z+mgcosγsinμ)/(mV)
f4(x)=secβ(ωxcosα-ωysinα)+[(tanβ+tanγsinμ)(Y-mgcosγcosμ)
+(Z+mgcosγsinμ)tanγcosμ]/(mV)
h(x)=[α β μ]T
其中α,β,μ分别为攻角、侧滑角和速度倾侧角,X,Y,Z为速度坐标系下阻力、升力和侧力,V为飞行器的速度。m为飞行器质量,Ixx,Iyy,Izz,Ixy为飞行器对机体坐标系各轴的转动惯量以及惯量积,ωx,ωy,ωz分别是滚转角速率、偏航角速率和俯仰角速率,γ为弹道倾角。Mx,My,Mz为俯仰、偏航、滚转三个方向的力矩。d(x)表示包括参数摄动、外部扰动以及未建模动态等聚合不确定性,由于再入过程中速度快,大气环境变化剧烈,d(x)无法忽略。
步骤3,针对飞行器的标称模型(d(x)=0),将步骤2建立的再入飞行器模型转化为状态依赖参数(SDC)形式:
A(x),B的表达式如下。
式中,
步骤4,针对再入飞行器的标称模型,根据SDRE方法理论计算标称控制律u*。
给定最优指标 Q(x)和R(x)是权值系数,Q7×7(x)≥0,R3×3(x)>0,根据控制量和***的动态性能调整Q(x)和R(x)权阵。
解如下代数Ricatti方程得到P(x):
AT(x)P(x)+P(x)A(x)+Q(x)-P(x)B(x)R-1(x)BTP(x)=0
计算标称控制律u*:
u*=-R(x)-1B(x)TP(x)[x-xc] (2)
式中xc=[0,αc,βc,μc,0,0,0]T,αc,βc,μc为制导***给出的姿态角指令。
步骤5,设计积分滑模面s如下:
s=Cx+z (3)
sig(s)和sign(s)定义如下:
sig(s)=[sig(s1),sig(s2),sig(s3)]T
sign(s)=[sign(s1),sign(s2),sign(s3)]T
步骤6,计算最优积分滑模控制力矩u。
u由标称控制律u*以及积分滑模切换项usw组成:
u=u*+usw (5)
其中u*为步骤4得到的标称控制律;
usw=-(CB)-1·[k1sig(s)+k2∫sign(s)]
s是步骤5中设计的积分滑模面。k1,k2为常值参数,满足k2≥4L。
步骤7,控制分配,得到舵偏角指令:
δ=G-1u
其中δ=[δe,δa,δr]T,δe,δa,δr分别为升降舵、副翼、方向舵的偏角,u=[Mx,My,Mz]T为步骤6得到的控制力矩,G是转换矩阵,由气动参数决定。
步骤8,将步骤7得到的舵偏角指令输入飞行器,对其进行姿态控制;同时,飞行器输出当前飞行器的各个状态ω,V,X,Y,Z,作为姿态控制的输入,重复步骤1-步骤8。
从而使得飞行器实现利用实际的姿态角Ω=[α β μ]T跟踪制导***给出的姿态角指令Ωc=[αc,βc,μc]T的目的。
根据上述最优积分滑模姿态控制方法,提出一种最优积分滑模姿态控制器,包括状态向量生成模块、状态依赖参数化模块、控制参数选择模块、SDRE求解模块、滑模面计算模块、标称控制律计算模块、最优积分滑模控制律计算模块和控制分配器。其中,状态向量生成模块与标称控制律计算模块连接,状态依赖参数化模块的输出分别连至SDRE求解模块、标称控制律计算模块、滑模面计算模块;控制参数选择模块的输出分别连至SDRE求解模块、滑模面计算模块;SDRE求解模块连接标称控制律计算模块;滑模面计算模块、标称控制律计算模块和滑模面计算模块输出至最优积分滑模控制律计算模块;最优积分滑模控制律计算模块连接控制分配器。
状态向量生成模块接收飞行器的速度V、姿态角速度ω=[ωx,ωy,ωz]T、制导***的姿态角指令Ωc=[αc,βc,μc]T、飞行器实际的姿态角Ω=[α β μ]T,生成状态向量x=[V,α,β,μ,ωx,ωv,ωz]T。
状态依赖参数化模块接收飞行器实际的姿态角Ω=[α β μ]T、速度V、姿态角速度ω=[ωx,ωy,ωz]T、速度坐标系下阻力、升力和侧力X,Y,Z,将再入飞行器模型转化为状态依赖参数(SDC)形式。
SDRE求解模块接收状态依赖参数化模块的输出A(x),B,结合选择的控制器参数,通过求解SDRE方程得到P(x)。
标称控制律计算模块接收状态依赖参数矩阵A(x),B,SDRE方程的解P(x),以及状态x,得到标称控制量u*。
滑模面计算模块接收状态依赖参数矩阵A(x),B,参数C,以及标称控制律u*,得到滑模面s。
最优积分滑模控制律计算模块接收标称控制律u*,滑模面s,以及控制参数k1,k2,得到控制力矩Mx,My,Mz。
控制分配模块将得到的控制力矩输出u=[Mx,My,Mz]T分配至舵面执行机构,得到舵偏角指令δ=[δe,δa,δr]T。
有益效果
1.标称控制律u*决定了***的响应,当存在不确定性时,***依然能够实现与标称***同样的相应。
2.***能够实现期望的最优指标,通过调整Q,R权阵,可有效折衷控制量和***的动态性能,且SDRE标称控制律完全保留了飞行控制***的非线性。
3.积分滑模面的引入,抑制了再入飞行中不确定性,保证了***的鲁棒性。
4.采用二阶滑模的思想设计控制律,有效减弱了抖振。
本发明设计的姿态控制器不仅能保证期望的指标,而且具有较好的鲁棒性。
附图说明
图1为本发明的最优积分滑模姿态控制方法的流程图;
图2为本发明的再入飞行器姿态控制***结构图;
图3为本发明的最优积分滑模姿态控制器结构图;
图4为具体实施方式中攻角跟踪误差;
图5为具体实施方式中侧滑角跟踪误差;
图6为具体实施方式中倾侧角跟踪误差;
图7为具体实施方式中舵面响应曲线。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实施例加以进一步说明。
本发明设计的控制器目的是在获得非线性姿态控制***最优指标的前提下,保证***的鲁棒性。
根据再入飞行器的模型,不考虑不确定性,再入飞行器的标称模型可以表示为:
y=h(x)
首先将再入飞行器的标称模型转化为了状态依赖参数形式:
然后,选择指标函数
根据SDRE理论,可得到使上述指标满足的标称控制器:
u*=-R(x)-1B(x)TP(x)[x-xc]
该控制器即为“SDRE标称控制律计算”中得到的标称控制律。针对再入飞行器的标称***,该控制器能够保证***的局部最优。考虑到再入过程中环境变化剧烈,***存在不确定性,若只采用上述标称控制律,***性能会偏离最优指标,因此通过引入积分滑模的方法保证控制器的鲁棒性。
在标称控制律的基础上增加滑模切换项,形成最优积分滑模控制律如式(5)所示。
对(3)所示的滑模面求导并将控制量代入,
上式可写作
ueq=[CB]-1[CBu*-Cd(x')]
代入不确定***(1)得:
因此,不确定***在发明的控制器控制下的性能与SDRE控制律u*控制下的性能完全一致,也即能实现最优性能指标。
本实施例通过在再入飞行器姿态控制仿真平台上进行实验,以验证本发明提出的最优积分滑模具有良好的性能。根据实施方案中的设计步骤,在飞行器姿态控制仿真平台上根据设计的最优积分滑模姿态控制方法搭建姿态控制器如图2所示,控制器的输入为制导***给出的姿态角信号Ωc、飞行器的实际姿态角Ω、速度坐标系下阻力、升力和侧力X,Y,Z、飞行器的速度V、以及飞行器的姿态角速度ω=[ωx,ωy,ωz]T,输出为舵面偏角。根据图2搭建再入飞行器姿态控制***。
在本实施例中,飞行器姿态角初始状态如下:攻角初值α=0°,侧滑角初值β=2°,速度倾侧角μ=0°,给定常值姿态角跟踪目标αc=4°,βc=0°,μc=10°。控制器参数C=[03×4 I3],k1=0.1,k2=0.05。SDRE标称控制律设计时,选择权值系数Q=diag{0 500 1000 500 0 0 0},R=diag{1 0.1 0.01},可以看出,由于只需要对状态α,β,μ进行控制,所以Q阵中相关的项系数较大,对于无需控制的状态,对应项取0。通过调节Q和R矩阵,可以对控制量和控制性能进行权衡,既保证良好的控制性能,又避免控制量饱和。
由于再入飞行器飞行条件大范围变化,且常常具有气动参数摄动等不确定性,因此对于再入飞行器的姿态控制问题,不仅要检验标称情况下的控制性能,还需要检验控制器在环境参数剧烈变化和***具有较强不确定性的情况下,能否进行鲁棒、精确地控制。为检验控制***的鲁棒性,本实施例中对气动系数、大气环境进行拉偏实验以模拟恶劣的环境,将大气密度负向拉偏20%,气动力系数负向拉偏10%,力矩系数负向拉偏30%,并且对xyz方向的力矩系数分别施加-0.001,-0.001,-0.0001的常值偏差。
采用发明的最优积分滑模对飞行器进行姿态控制,为了进行对比,与采用SDRE控制器对再入飞行器进行姿态控制的结果进行对比,姿态角跟踪误差如图4-图6所示。图中分别给出了采用SDRE控制方法对标称模型及不确定模型的控制结果,以及发明的最优积分滑模控制器对不确定模型的控制结果。可以看出,没有参数不确定性及外部干扰时,SDRE方法跟踪的效果较好,姿态角跟踪误差约3s即可达到0,同时超调较小,整个动态响应过程良好。然而当存在不确定性时,***的跟踪性能明显变差,飞行器的姿态偏离了给定值,出现了明显的静差。而采用最优积分滑模时,与只采用SDRE方法相比,明显具有更好的动态品质和跟踪性能,且最优积分滑模控制律对于不确定性***的控制效果与SDRE控制律对于标称***的控制律几乎重合,提出的姿态控制算法具有较强的鲁棒性。
发明的控制器输出升降舵副翼、方向舵的舵面偏角δe,δa,δr,由于舵偏范围受实际硬件条件限制:-30°≤δe,δa,δr≤30°。图7给出了采用SDRE控制器控制标称***、采用SDRE控制器控制不确定***、以及采用最优积分滑模姿态控制器控制不确定***时的舵面偏角。由仿真结果可以看出,由于选取了合适的Q、R参数,舵面没有达到饱和,而且发明的控制器输出控制量光滑,避免了抖振带来的不良影响。
Claims (2)
1.再入飞行器的最优积分滑模姿态控制方法,其特征在于:包括以下步骤:
步骤1,生成飞行器的状态向量;
结合飞行器的实际姿态角Ω=[α,β,μ]T,姿态角速度ω=[ωx,ωy,ωz]T,以及速度V,组成状态向量x:x=[V,α,β,μ,ωx,ωy,ωz]T;
步骤2,建立再入飞行器的数学模型:
y=h(x)
其中,状态向量x=[V α β μ ωx ωy ωz]T,控制力矩u=[Mx My Mz]T,输出变量y[α β μ]T,f(x)=[f1(x) f2(x)...f7(x)]T;
f1(x)=(-X-mgsinγ)/m
f2(x)=ωz+tanβ(ωysinα-ωxcosα)-(Y-mgcosγcosμ)/(mVcosβ)
f3(x)=ωxsinα+ωycosα+(Z+mgcosγsinμ)/(mV)
f4(x)=secβ(ωxcosα-ωysinα)+[(tanβ+tanγsinμ)(Y-mgcosγcosμ)
+(Z+mgcosγsinμ)tanγcosμ]/(mV)
h(x)=[α β μ]T
其中α,β,μ分别为攻角、侧滑角和速度倾侧角,X,Y,Z为速度坐标系下阻力、升力和侧力,V为飞行器的速度;m为飞行器质量,Ixx,Iyy,Izz,Ixy为飞行器对机体坐标系各轴的转动惯量以及惯量积,ωx,ωy,ωz分别是滚转角速率、偏航角速率和俯仰角速率,γ为弹道倾角;Mx,My,Mz为俯仰、偏航、滚转三个方向的力矩;d(x)表示包括参数摄动、外部扰动以及未建模动态等聚合不确定性;
步骤3,针对d(x)=0的飞行器标称模型,将步骤2建立的再入飞行器模型转化为状态依赖参数形式:
A(x),B的表达式如下;
式中,
步骤4,针对再入飞行器的标称模型,根据SDRE方法理论计算标称控制律u*;
给定最优指标 Q(x)和R(x)是权值系数,Q7×7(x)≥0,R3×3(x)>0,根据控制量和***的动态性能调整Q(x)和R(x)权阵;
解如下代数Ricatti方程得到P(x):
AT(x)P(x)+P(x)A(x)+Q(x)-P(x)B(x)R-1(x)BTP(x)=0
计算标称控制律u*:
u*=-R(x)-1B(x)TP(x)[x-xc] (2)
式中xc=[0,αc,βc,μc,0,0,0]T,αc,βc,μc为制导***给出的姿态角指令;
步骤5,设计积分滑模面s:
s=Cx+z (3)
sig(s)和sign(s)定义如下:
sig(s)=[sig(s1),sig(s2),sig(s3)]T
sign(s)=[sign(s1),sign(s2),sign(s3)]T
步骤6,计算最优积分滑模控制力矩u;
u由标称控制律u*以及积分滑模切换项usw组成:
u=u*+usw (5)
其中u*为步骤4得到的标称控制律;
usw=-(CB)-1·[k1 sig(s)+k2∫sign(s)]
步骤7,控制分配,得到舵偏角指令:
δ=G-1u
其中δ=[δe,δa,δr]T,δe,δa,δr分别为升降舵、副翼、方向舵的偏角,u=[Mx,My,Mz]T为步骤6得到的姿态控制力矩,G是转换矩阵,由气动参数决定。
步骤8,将步骤7得到的舵偏角指令输入飞行器,对其进行姿态控制;同时,飞行器输出当前飞行器的各个状态ω,V,X,Y,Z,作为姿态控制的输入,重复步骤1-步骤8;从而使得飞行器实现利用实际的姿态角Ω=[α β μ]T跟踪制导***给出的姿态角指令Ωc=[αc,βc,μc]T的目的。
2.根据权利要求1所述的方法设计的一种再入飞行器的最优积分滑模姿态控制器,其特征在于:包括状态向量生成模块、状态依赖参数化模块、控制参数选择模块、SDRE求解模块、滑模面计算模块、标称控制律计算模块、最优积分滑模控制律计算模块和控制分配器;其中,状态向量生成模块与标称控制律计算模块连接,状态依赖参数化模块的输出分别连至SDRE求解模块、标称控制律计算模块、滑模面计算模块;控制参数选择模块的输出分别连至SDRE求解模块、滑模面计算模块;SDRE求解模块连接标称控制律计算模块;滑模面计算模块、标称控制律计算模块和滑模面计算模块输出至最优积分滑模控制律计算模块;最优积分滑模控制律计算模块连接控制分配器;
状态向量生成模块接收飞行器的速度V、姿态角速度ω=[ωx,ωy,ωz]T、制导***的姿态角指令Ωc=[αc,βc,μc]T、飞行器实际的姿态角Ω=[α β μ]T,生成状态向量x=[V,α,β,μ,ωx,ωy,ωz]T;
状态依赖参数化模块接收飞行器实际的姿态角Ω=[α β μ]T、速度V、姿态角速度ω=[ωx,ωy,ωz]T、速度坐标系下阻力、升力和侧力X,Y,Z,将再入飞行器模型转化为状态依赖参数形式;
SDRE求解模块接收状态依赖参数化模块的输出A(x),B,结合选择的控制器参数,通过求解SDRE方程得到P(x);
标称控制律计算模块接收状态依赖参数矩阵A(x),B,SDRE方程的解P(x),以及状态x,得到标称控制量u*;
滑模面计算模块接收状态依赖参数矩阵A(x),B,参数C,以及标称控制律u*,得到滑模面s;
最优积分滑模控制律计算模块接收标称控制律u*,滑模面s,以及控制参数k1,k2,得到控制力矩Mx,My,Mz;
控制分配模块将得到的控制力矩输出u=[Mc,My,Mz]T分配至舵面执行机构,得到舵偏角指令δ=[δe,δa,δr]T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210367288.7A CN102862686B (zh) | 2012-09-28 | 2012-09-28 | 再入飞行器的最优积分滑模姿态控制方法及控制器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210367288.7A CN102862686B (zh) | 2012-09-28 | 2012-09-28 | 再入飞行器的最优积分滑模姿态控制方法及控制器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102862686A true CN102862686A (zh) | 2013-01-09 |
CN102862686B CN102862686B (zh) | 2014-08-06 |
Family
ID=47441881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210367288.7A Expired - Fee Related CN102862686B (zh) | 2012-09-28 | 2012-09-28 | 再入飞行器的最优积分滑模姿态控制方法及控制器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102862686B (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699131A (zh) * | 2013-12-26 | 2014-04-02 | 北京控制工程研究所 | 一种卫星控制***离散积分滑模容错控制方法 |
CN103853157A (zh) * | 2014-03-19 | 2014-06-11 | 湖北蔚蓝国际航空学校有限公司 | 一种基于自适应滑模的飞行器姿态控制方法 |
CN104134008A (zh) * | 2014-08-08 | 2014-11-05 | 北京航天自动控制研究所 | 飞行器姿态运动通道间运动耦合特性的交联影响评估方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
CN104462810A (zh) * | 2014-12-05 | 2015-03-25 | 哈尔滨工业大学 | 一种适用于轮控卫星姿态机动与跟踪控制的sdre参数调节方法 |
CN104503233A (zh) * | 2014-11-27 | 2015-04-08 | 哈尔滨工业大学 | 适用于卫星姿态控制的干扰力矩辨识方法 |
CN104634182A (zh) * | 2014-12-16 | 2015-05-20 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN104699918A (zh) * | 2015-03-27 | 2015-06-10 | 北京润科通用技术有限公司 | 一种计算舵偏角的方法及*** |
CN104898431A (zh) * | 2015-06-10 | 2015-09-09 | 北京理工大学 | 一种基于扰动观测器的再入飞行器有限时间控制方法 |
CN104898685A (zh) * | 2015-05-27 | 2015-09-09 | 中国人民解放军海军航空工程学院 | 一种基于概率最优的飞行器俯仰通道攻角跟踪控制方法 |
CN104932531A (zh) * | 2015-06-24 | 2015-09-23 | 南京航空航天大学 | 一种基于滑模控制的四旋翼飞行器的最优抗输入饱和控制方法 |
CN104950671A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 基于自适应模糊的再入飞行器pid型滑模姿态控制方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN104950672A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种最优积分滑模控制方法 |
CN105138001A (zh) * | 2015-09-10 | 2015-12-09 | 中国人民解放军国防科学技术大学 | 一种四旋翼飞行器姿态控制方法 |
CN105404304A (zh) * | 2015-08-21 | 2016-03-16 | 北京理工大学 | 基于归一化神经网络的航天器容错姿态协同跟踪控制方法 |
CN106248082A (zh) * | 2016-09-13 | 2016-12-21 | 北京理工大学 | 一种飞行器自主导航***及导航方法 |
CN106950982A (zh) * | 2017-02-16 | 2017-07-14 | 北京临近空间飞行器***工程研究所 | 再入飞行器姿控动力***高空力矩特性辨识方法 |
CN108287476A (zh) * | 2018-02-06 | 2018-07-17 | 哈尔滨工业大学 | 基于高阶滑模控制和扰动观测器的空间翻滚非合作目标自主交会制导方法 |
CN109085848A (zh) * | 2018-08-02 | 2018-12-25 | 西北工业大学 | 空空导弹直接力/气动力有限时间抗饱和控制方法 |
CN110109473A (zh) * | 2019-04-29 | 2019-08-09 | 北京理工大学 | 可应用于旋转飞行器的非线性滚转稳定控制***及方法 |
CN110618608A (zh) * | 2018-06-20 | 2019-12-27 | 河南科技大学 | 一种复合制导跟踪控制方法与装置 |
CN110687801A (zh) * | 2019-11-22 | 2020-01-14 | 深圳信息职业技术学院 | 一种基于前向制导的自适应积分滑模三维制导律设计方法 |
CN111796568A (zh) * | 2020-07-02 | 2020-10-20 | 北京航空航天大学 | 面向控制***的积分滑模控制方法、装置及设备 |
CN113900440A (zh) * | 2021-07-21 | 2022-01-07 | 中国电子科技集团公司电子科学研究院 | 一种无人机控制律设计方法、装置及可读存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364085A (zh) * | 2008-09-12 | 2009-02-11 | 北京控制工程研究所 | 非线性微分黄金分割自适应控制方法 |
CN102139767A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行过程中的横向调整方法 |
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
-
2012
- 2012-09-28 CN CN201210367288.7A patent/CN102862686B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364085A (zh) * | 2008-09-12 | 2009-02-11 | 北京控制工程研究所 | 非线性微分黄金分割自适应控制方法 |
CN102139767A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行过程中的横向调整方法 |
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
Non-Patent Citations (2)
Title |
---|
张军等: "一种高超声速飞行器的非线性再入姿态控制方法", 《空间控制技术与应用》 * |
张军等: "基于SDRE方法的挠性航天器姿态控制", 《宇航学报》 * |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699131A (zh) * | 2013-12-26 | 2014-04-02 | 北京控制工程研究所 | 一种卫星控制***离散积分滑模容错控制方法 |
CN103699131B (zh) * | 2013-12-26 | 2016-05-04 | 北京控制工程研究所 | 一种卫星控制***离散积分滑模容错控制方法 |
CN103853157A (zh) * | 2014-03-19 | 2014-06-11 | 湖北蔚蓝国际航空学校有限公司 | 一种基于自适应滑模的飞行器姿态控制方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
CN104134008A (zh) * | 2014-08-08 | 2014-11-05 | 北京航天自动控制研究所 | 飞行器姿态运动通道间运动耦合特性的交联影响评估方法 |
CN104503233A (zh) * | 2014-11-27 | 2015-04-08 | 哈尔滨工业大学 | 适用于卫星姿态控制的干扰力矩辨识方法 |
CN104462810B (zh) * | 2014-12-05 | 2017-06-16 | 哈尔滨工业大学 | 一种适用于轮控卫星姿态机动与跟踪控制的sdre参数调节方法 |
CN104462810A (zh) * | 2014-12-05 | 2015-03-25 | 哈尔滨工业大学 | 一种适用于轮控卫星姿态机动与跟踪控制的sdre参数调节方法 |
CN104634182A (zh) * | 2014-12-16 | 2015-05-20 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN104634182B (zh) * | 2014-12-16 | 2016-02-10 | 北京控制工程研究所 | 一种跳跃式再入标准弹道在线修正的跟踪制导方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN104656450B (zh) * | 2015-01-20 | 2017-04-19 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN104699918A (zh) * | 2015-03-27 | 2015-06-10 | 北京润科通用技术有限公司 | 一种计算舵偏角的方法及*** |
CN104699918B (zh) * | 2015-03-27 | 2017-08-11 | 北京润科通用技术有限公司 | 一种计算舵偏角的方法及*** |
CN104898685B (zh) * | 2015-05-27 | 2017-09-19 | 中国人民解放军海军航空工程学院 | 一种基于概率最优的飞行器俯仰通道攻角跟踪控制方法 |
CN104898685A (zh) * | 2015-05-27 | 2015-09-09 | 中国人民解放军海军航空工程学院 | 一种基于概率最优的飞行器俯仰通道攻角跟踪控制方法 |
CN104898431A (zh) * | 2015-06-10 | 2015-09-09 | 北京理工大学 | 一种基于扰动观测器的再入飞行器有限时间控制方法 |
CN104950672A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种最优积分滑模控制方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN104950671A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 基于自适应模糊的再入飞行器pid型滑模姿态控制方法 |
CN104950672B (zh) * | 2015-06-10 | 2017-09-08 | 北京理工大学 | 一种最优积分滑模控制方法 |
CN104932531A (zh) * | 2015-06-24 | 2015-09-23 | 南京航空航天大学 | 一种基于滑模控制的四旋翼飞行器的最优抗输入饱和控制方法 |
CN105404304A (zh) * | 2015-08-21 | 2016-03-16 | 北京理工大学 | 基于归一化神经网络的航天器容错姿态协同跟踪控制方法 |
CN105404304B (zh) * | 2015-08-21 | 2018-05-25 | 北京理工大学 | 基于归一化神经网络的航天器容错姿态协同跟踪控制方法 |
CN105138001A (zh) * | 2015-09-10 | 2015-12-09 | 中国人民解放军国防科学技术大学 | 一种四旋翼飞行器姿态控制方法 |
CN105138001B (zh) * | 2015-09-10 | 2017-08-25 | 中国人民解放军国防科学技术大学 | 一种四旋翼飞行器姿态控制方法 |
CN106248082B (zh) * | 2016-09-13 | 2019-06-04 | 北京理工大学 | 一种飞行器自主导航***及导航方法 |
CN106248082A (zh) * | 2016-09-13 | 2016-12-21 | 北京理工大学 | 一种飞行器自主导航***及导航方法 |
CN106950982A (zh) * | 2017-02-16 | 2017-07-14 | 北京临近空间飞行器***工程研究所 | 再入飞行器姿控动力***高空力矩特性辨识方法 |
CN108287476A (zh) * | 2018-02-06 | 2018-07-17 | 哈尔滨工业大学 | 基于高阶滑模控制和扰动观测器的空间翻滚非合作目标自主交会制导方法 |
CN110618608A (zh) * | 2018-06-20 | 2019-12-27 | 河南科技大学 | 一种复合制导跟踪控制方法与装置 |
CN110618608B (zh) * | 2018-06-20 | 2022-04-15 | 河南科技大学 | 一种复合制导跟踪控制方法与装置 |
CN109085848A (zh) * | 2018-08-02 | 2018-12-25 | 西北工业大学 | 空空导弹直接力/气动力有限时间抗饱和控制方法 |
CN109085848B (zh) * | 2018-08-02 | 2021-05-07 | 西北工业大学 | 空空导弹直接力/气动力有限时间抗饱和控制方法 |
CN110109473A (zh) * | 2019-04-29 | 2019-08-09 | 北京理工大学 | 可应用于旋转飞行器的非线性滚转稳定控制***及方法 |
CN110687801A (zh) * | 2019-11-22 | 2020-01-14 | 深圳信息职业技术学院 | 一种基于前向制导的自适应积分滑模三维制导律设计方法 |
CN111796568A (zh) * | 2020-07-02 | 2020-10-20 | 北京航空航天大学 | 面向控制***的积分滑模控制方法、装置及设备 |
CN113900440A (zh) * | 2021-07-21 | 2022-01-07 | 中国电子科技集团公司电子科学研究院 | 一种无人机控制律设计方法、装置及可读存储介质 |
CN113900440B (zh) * | 2021-07-21 | 2023-03-14 | 中国电子科技集团公司电子科学研究院 | 一种无人机控制律设计方法、装置及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102862686B (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102862686B (zh) | 再入飞行器的最优积分滑模姿态控制方法及控制器 | |
CN106444799B (zh) | 基于模糊扩张状态观测器和自适应滑模的四旋翼无人机控制方法 | |
CN104950898B (zh) | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 | |
CN105416276B (zh) | 基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 | |
Huang et al. | Fixed-time autonomous shipboard landing control of a helicopter with external disturbances | |
CN102929151A (zh) | 一种基于指数时变二阶滑模的再入飞行姿态控制方法 | |
CN103105850B (zh) | 一种近空间飞行器故障诊断与容错控制方法 | |
CN103645636B (zh) | 一种pid控制器参数优化方法 | |
CN108594639B (zh) | 一种基于强化学习的全垫升气垫船航迹跟踪控制方法 | |
CN110347170B (zh) | 可重复使用运载器再入段鲁棒容错制导控制***及工作方法 | |
CN104950671A (zh) | 基于自适应模糊的再入飞行器pid型滑模姿态控制方法 | |
CN105785762A (zh) | 一种基于自适应反步滑模的两轴惯性稳定平台高精度控制方法 | |
CN105629734B (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN105607473B (zh) | 小型无人直升机的姿态误差快速收敛自适应控制方法 | |
CN101937233A (zh) | 近空间高超声速飞行器非线性自适应控制方法 | |
CN103558857A (zh) | 一种btt飞行器的分布式复合抗干扰姿态控制方法 | |
CN105786024A (zh) | 一种基于模型误差补偿的机载光电平台高精度跟踪控制器及其跟踪控制方法 | |
CN103197543B (zh) | 基于运动状态综合识别的高速飞行器自适应控制方法 | |
CN102929283A (zh) | 基于sdre的再入飞行器自适应最优滑模姿态控制方法 | |
CN105759609B (zh) | 一种基于显式非线性模型预测控制的舰载机自主着舰方法 | |
CN107450323B (zh) | 高超声速飞行器再入段神经网络复合学习控制方法 | |
CN111045440B (zh) | 一种高超声速飞行器俯冲段快速滚转控制方法 | |
CN102654772A (zh) | 一种基于控制力受限情况下飞行器航迹倾角反演控制方法 | |
CN108919824A (zh) | 舰载无人机纵向着舰控制方法 | |
CN111142550B (zh) | 民用飞机辅助驾驶控制方法、***及飞行品质评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140806 Termination date: 20190928 |
|
CF01 | Termination of patent right due to non-payment of annual fee |