公开/公告号CN102176160A
专利类型发明专利
公开/公告日2011-09-07
原文格式PDF
申请/专利权人 北京控制工程研究所;
申请/专利号CN201010623823.1
发明设计人 孟斌;
申请日2010-12-31
分类号G05D1/00;
代理机构中国航天科技专利中心;
代理人杨虹
地址 100080 北京市2729信箱
入库时间 2023-12-18 03:08:57
法律状态公告日
法律状态信息
法律状态
2012-09-26
授权
授权
2011-11-16
实质审查的生效 IPC(主分类):G05D1/00 申请日:20101231
实质审查的生效
2011-09-07
公开
公开
技术领域
本发明涉及一种高超声速飞行器控制方法,特别是涉及一种基于偏量特征模型的高超声速飞行器自适应学习控制方法,属于飞行器控制技术领域。
背景技术
高超声速飞行器技术是航空史上继发明飞机、突破声障之后的第三个划时代的里程碑,同时也开辟了进入太空的新模式。由于高超声速飞行器的独特优势,作为武器将在未来战争的快速打击、远程突防中占据重要的地位,因此高超声速飞行器的研究受到各军事大国政府的极大重视。高超声速飞行器必须采用特殊的飞行模式,使得其飞行包线更加复杂,成为一个参数和状态相互耦合、参数大范围变化、多通道耦合、具有严重不确定性的复杂非线性系统。这对控制理论和方法的研究提出了新的挑战。
基于特征模型的黄金分割控制方法是吴宏鑫院士提出的,经过20多年的研究,在理论和应用上均取得了重要进展,形成了一套实用性很强的自适应控制理论和方法。该方法需要辨识的参数少,可以保证闭环系统的暂态性能和稳态性能,具有强鲁棒性和自适应性。基于特征模型的全系数自适应控制方法的基本思想是,首先建立系统的特征模型,然后按照特征模型参数设计全系数自适应控制律。所建立的系统的特征模型为下述二阶时变差分方程:
y(k+2)=f1(k)y(k+1)+f2(k)y(k)+f3(k)u(k)
其中u(k),y(k)分别表示系统的输入和输出;系数fi(k),i=1,2,3,是慢时变的且其范围可事先确定。
对于上述二阶时变差分方程,使用加权最小二乘方法等辨识特征模型参数fi,i=1,2,3,然后将其投影到给定范围内。辨识结果记为i=1,2,3。
进而,设计全系数自适应控制律u=u0+ug+ui+ud,其中,
ui(k)=ui(k-1)-kie(k)
ud(k)=-kde(k)
l1=0.382,l2=0.618,e(k)=y(k)-yr(k),ki1>>ki2>0,或者,ki1,ki2,cd,ld为所需调整参数,yr(k)为跟踪目标函数。u0,ug,ui,和ud分别称为维持/跟踪控制律,黄金分割控制,逻辑积分控制和逻辑微分控制。黄金分割自适应控制律是该方法中特有的一种反馈控制律。
由于高超声速飞行器速度快,参数变化范围大,当设计上述基于特征模型的全系数自适应控制律,出现了执行机构舵抖动过快的问题。对于本发明所给出基于偏量特征模型的高超声速飞行器自适应学习控制方法,没有公开的具有完整实用意义的方法。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种能抑制高超声速飞行器舵抖动的偏量特征模型的高超声速飞行器自适应学习控制方法。
本发明的技术解决方案是:基于偏量特征模型的高超声速飞行器自适应学习控制方法,通过以下步骤实现:
第一步,根据公式组(1)的高超声速飞行器动力学方程,建立公式(2)的偏量特征模型,
其中
Δyi(k+2)=ai1(k)Δyi(k+1)+ai2(k)Δyi(k)+bi(k)Δui(k),i=1,2,3 (2),
其中y1(k)为第k个采样周期的攻角α,y2(k)为第k个采样周期的侧滑角β,y3(k)为第k个采样周期的高超声速飞行器的滚动角φ,k为自然数,Δyi(k)=yi(k+1)-yi(k),Δui(k)=ui(k+1)-ui(k),u1(k)为第k个采样周期的高超声速飞行器的升降舵偏转角δe,u2(k)为第k个采样周期的高超声速飞行器的方向舵偏转角δr,u3(k)为第k个采样周期的高超声速飞行器的左右升降舵差动角δa,ai1、ai2和bi为偏量特征模型的系数,且ai1∈[2-h2b,2+h2b],ai2∈[-1-h2b,-1+h2b],bi∈[-h2b,h2b],h为采样周期,b为大于0的常数;
第二步,利用公式(3)的对偶辨识法对第一步中得到的偏量特征模型的系数ai1、ai2和bi进行辨识,
i=1,2,3 (3)
其中μi为系数修正参数且μi>0,q-1为滞后一步算子,l1=0.382,l2=0.618,和分别为ai1、ai2和bi的辨识值,为的转置,
第三步,根据第二步得到的ai1、ai2和bi辨识值和利用公式(4)得到黄金分割自适应学习控制律Δui(k),
其中,Δui(k)=ui(k+1)-ui(k),u1(k)为第k个采样周期的高超声速飞行器的升降舵偏转角δe,u2(k)为第k个采样周期的高超声速飞行器的方向舵偏转角δr,u3(k)为第k个采样周期的高超声速飞行器的左右升降舵差动角δa,u1,u2,u3分别表示δe,δr,δa,λi为控制律修正参数且λi>0,ei(k)=Δyi(k)-Δyir(k),yir(k)为跟踪目标函数,Δyir(k)=yir(k+1)-yir(k);
第四步,利用第三步得到黄金分割自适应学习控制律Δui(k)返回到公式(1)的高超声速飞行器动力学方程,控制高超声速飞行器的方向舵偏转角δr、升降舵偏转角δe和左右升降舵差动角δa。
本发明与现有技术相比有益效果为:
(1)本发明引入了两个可变的修正参数,可以通过适当调整修正参数,对每一步控制增量进行限制,解决了执行机构舵抖动过快的问题;
(2)本发明建立了偏量特征模型,为黄金分割自适应学习控制律的设计奠定了基础;
(3)本发明给出了对偶参数辨识方法,可以减小系统控制的暂态误差;
(4)本发明给出了基于偏量特征模型的黄金分割自适应学习控制设计方法,从而可以通过调整控制律修正参数,限制控制增量的大小,降低了特征模型系数和状态的相关度,解决了特征模型的系数辨识问题;
(5)本发明针对高超声速飞行器,设计基于偏量特征模型的黄金分割自适应学习控制律和对偶辨识算法,减小了暂态误差和舵的抖动。
附图说明
图1为本发明流程图。
具体实施方式
本发明如图1所示,通过以下步骤实现:
1、根据公式组的高超声速飞行器动力学方程,建立偏量特征模型。
忽略地球自转、曲率以及风的影响,且忽略气动舵面对飞行器气动力的影响,利用线性化的气动系数,高超声速飞行器姿态动力学为:
其中
下面对被控输出α,β和φ求2阶导数,由公式组(1)可得
其中,
lay=-(I1cosα+I5sinα+cosα)tanβωz-(I2cosα+I6sinα-sinα)tanβωx+kyy
lβx=(I2sinα-I5cosα+cosα)ωy-2sinαcosαtanβωz-cos2αtanβωx+(kxxsinα-kzxcosα)
lβz=(I1sinα-I6cosα+sinα)ωy-sin2αtanβωz+(kxzsinα-kzzcosα)
lφx=(I2+I5cosφtanθ+cosφtanθ)ωy+(I3sinφtanθ-sinφtanθ)ωz+
(I4sinφtanθ)ωx+(kxx+kzxcosφtanθ)
利用小角度假设:
sinα≈α,cosα≈1,tanβ≈β
以及
与公式组(1)中的第三个公式联立,得到
将(1-4)~(1-7)代入(1-1)~(1-3),得到
其中,Δα,Δβ,Δφ是建模误差,
将(1-8)~(1-10)记为:
其中,y1,y2,y3分别表示α,β,φ,u1,u2,u3分别表示δe,δr,δa,Δ1,Δ2,Δ3分别表示Δα,Δβ,Δφ,fi(·)表示非线性函数。
将(1-11)离散化,可得
yi(k+2)-2yi(k+1)+yi(k)=h2fi(k)+h2Δi(k)
其中,h表示采样周期,
fi(k)=fi(y1(k),y1(k+1),y2(k),y2(k+1),y3(k),y3(k+1),ui(k),ui(k+1))
进而可得:
Δyi(k+2)-2Δyi(k+1)+Δyi(k)=h2fi(k+1)-h2fi(k)+h2Δi(k+1)-h2Δi(k)
其中,Δyi(k)=yi(k+1)-yi(k),由上式易得:
其中,符号表示在y(k)和y(k+1)之间某一点的值,ξi表示误差。由于Δui≠0,则存在ηi,使得
ξi=ηiΔui
进而有
将上式记为:
Δyi(k+2)=ai1(k)Δyi(k+1)+ai2(k)Δyi(k)+bi(k)Δui(k),i=1,2,3 (2)
其中,
由的上界为b,(b为大于0的常数)易得:
ai1∈[2-h2b,2+h2b],a2∈[-1-h2b,-1+h2b],bi∈[-h2b,h2b]
公式(2)即为所建立的偏量特征模型。
其中,ai1、ai2和bi为偏量特征模型的系数,y1(k)为第k个采样周期的攻角α,y2(k)为第k个采样周期的侧滑角β,y3(k)为第k个采样周期的高超声速飞行器的滚动角φ,k为自然数,Δyi(k)=yi(k+1)-yi(k),Δui(k)=ui(k+1)-ui(k),u1(k)为第k个采样周期的高超声速飞行器的升降舵偏转角δe,u2(k)为第k个采样周期的高超声速飞行器的方向舵偏转角δr,u3(k)为第k个采样周期的高超声速飞行器的左右升降舵差动角δa。
2、利用对偶辨识法对偏量特征模型的系数ai1、ai2和bi进行辨识。
选取指标
J(Δui(k))=[Pi(q-1)Δyi(k+2)-Ri(q-1)Δyir(k+1)]2+λi(Δui(k))2,i=1,2,3 (3-1)
其中,yir(k)为跟踪目标函数,
Pi(q-1)=1-l2ai1q-1-l1ai2q-2,Ri(q-1)=l1ai1+l2ai2q-1
由
易得:
其中,ei(k)=Δyi(k)-Δyir(k),上式即为所建立的黄金分割自适应学习控制。
与(3-1)对偶选取参数辨识指标:
其中,μi为系数修正参数且μi>0,减小系数修正参数μi可以增强观测量对参数估计值的修正作用,增强算法对时变参数的跟踪能力,θi(k)=[ai1(k) ai2(k) bi(k)],表示θi(k)的辨识值,利用矩阵求逆引理,
可得
由于
设
则
则
i=1,2,3 (3-3)
公式(3-3)即为所建立的对偶参数辨识方法。
根据自适应控制确定性等价原理,使用时将公式(3-3)改为公式(3),
i=1,2,3 (3)
其中μi为系数修正参数且μi>0,q-1为滞后一步算子,l1=0.382,l2=0.618,和分别为ai1、ai2和bi的辨识值,为的转置,
3、将得到的ai1、ai2和bi辨识值和利用公式(4)得到黄金分割自适应学习控制律Δui(k),
根据自适应控制理论中的确定性等价原理,利用辨识参数代替系统参数,公式(3-2)可化为公式(4)。
其中,Δui(k)=ui(k+1)-ui(k),u1(k)为第k个采样周期的高超声速飞行器的升降舵偏转角δe,u2(k)为第k个采样周期的高超声速飞行器的方向舵偏转角δr,u3(k)为第k个采样周期的高超声速飞行器的左右升降舵差动角δa,`u1,u2,u3分别表示δe,δr,δa,λi为控制律修正参数且λi>0,减小算法中的参数λi可以增强对控制的修正作用,ei(k)=Δyi(k)-Δyir(k),yir(k)为跟踪目标函数,Δyir(k)=yir(k+1)-yir(k)。
4、利用黄金分割自适应学习控制律Δui(k)返回到公式组(1)的高超声速飞行器动力学方程,控制高超声速飞行器的方向舵偏转角δr、升降舵偏转角δe和左右升降舵差动角δa。
本发明未详细说明部分属本领域技术人员公知常识。
机译: 高超声速滑行飞行器的坡峰轨迹
机译: 高超声速飞行器的集成的向内转向进气口和喷嘴
机译: 高超声速飞行器的集成的向内转向进气口和喷嘴