测量噪声下的高超声速飞行器组合观测鲁棒控制.pdf
《测量噪声下的高超声速飞行器组合观测鲁棒控制.pdf》由会员分享,可在线阅读,更多相关《测量噪声下的高超声速飞行器组合观测鲁棒控制.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 45 卷第 2 期 2024 年 2 月宇 航 学 报Journal of AstronauticsNo.22024FebruaryVol.45测量噪声下的高超声速飞行器组合观测鲁棒控制张睿1,2,李世华1,魏振岩3,许斌2(1.东南大学自动化学院,南京 210096;2.西北工业大学自动化学院,西安 710129;3.北京机电工程研究所,北京 100074)摘要:针对高超声速飞行器大包络飞行时,风场等外部干扰和测量噪声导致控制系统性能变差甚至失稳的问题,提出了一种基于状态观测的高超声速飞行器鲁棒控制策略。考虑噪声和外部干扰影响下系统状态受扰的问题,设计了扩张状态观测器(ESO)和自适应K
2、alman滤波器(AKF)的特殊组合,在统计特性未知的噪声和外部干扰影响下准确估计系统状态。其中,自适应Kalman滤波器用于滤除噪声,为扩张状态观测器提供不受噪声影响的状态估计,扩张状态观测器连续为自适应Kalman滤波器提供外部干扰估计。基于系统状态估计设计线性高斯二次型控制器,实现高超声速飞行器精确轨迹跟踪。最后,通过仿真验证了该方法在测量噪声及风干扰影响下的鲁棒跟踪能力。关键词:高超声速飞行器;扩张状态观测器;自适应Kalman滤波器;鲁棒控制;测量噪声中图分类号:V249.1 文献标识码:A 文章编号:1000-1328(2024)02-0203-09 DOI:10.3873/j.i
3、ssn.1000-1328.2024.02.005Robust Control of Hypersonic Flight Vehicle Based on Combined State Observation Under Measurement NoisesZHANG Rui1,2,LI Shihua1,WEI Zhenyan3,XU Bin2(1.School of Automation,Southeast University,Nanjing 210096,China;2.School of Automation,Northwestern Polytechnical University,
4、Xian 710129,China;3.Beijing Research Institute of Mechanical&Electrical Technology,Beijing 100074,China)Abstract:In response to the poor performance and even instability of the control system caused by existed external disturbances and measurement noises,a robust control strategy based on state obse
5、rvation is proposed.Considering the unstable controller caused by using inaccurate system states under disturbances and noises,the special combination of extended state observer(ESO)and adaptive Kalman filter(AKF)is proposed to accurately estimate the system states which are contaminated by disturba
6、nces and noises that with unknown statistical properties.The AKF is used to offset noises and to provide the ESO with noiseless state estimation,while the ESO continuously provides the estimation of wind disturbances for the AKF.To achieve good tracking performance,the linear quadratic Gaussian cont
7、rol using the states estimation of observer is designed for HFV.Simulations are implemented to verify that robust tracking is obtained under the proposed method.Key words:Hypersonic flight vehicle;Extended state observer;Adaptive Kalman filter;Robust control;Measurement noises0引言高超声速飞行器(Hypersonic f
8、light vehicle,HFV)能够以马赫数超过5的高速度巡航飞行,在民用和军事领域发挥着至关重要的作用,越来越受到人们的重视。当前已经成功进行了 X-43A、X-51A和猎鹰HTV-2等型号的试飞验证。HFV的机身与推力系统设计呈现高度一体化,机身外形纤细,材料轻收稿日期:2023-08-10;修回日期:2023-11-19基金项目:江苏省卓博计划(2022ZB117);陕西省自然科学基础研究计划(2023JC-XJ-08);国家自然科学基金重点项目(61933010)宇航学报第 45 卷薄,气动、推力、弹性三者之间的耦合突出。文献1-2 分别给出了HFV锥体加速器模型和曲线拟合得到的面
9、向控制模型。考虑到飞行环境、气动参数和外界扰动的大范围变化,HFV控制过程中存在扰动和噪声干扰。为了使HFV在各种干扰下获得良好的跟踪性能,已有大量的先进控制方案研究。文献 3-4 将面向控制的HFV纵向通道模型分解为速度子系统和高度子系统,分别设计控制器实现参考轨迹跟踪。其中,针对速度子系统,设计了鲁棒自适应控制器。针对高度子系统,提出了动态逆控制,避免了传统反步控制的“计算爆炸”问题。文献 5 在反馈线性化控制的框架下,采用滑模控制来提高HFV控制系统的鲁棒性。文献 6 分别设计了有限时间终端滑模控制和高阶滑模控制以获得更快的跟踪误差收敛速度和更高的控制精度。考虑HFV大包络飞行时复杂飞行
10、环境导致的系统未知动力学,文献 7-9 分别采用神经网络对未知非线性动力学进行逼近,并将逼近结果进一步前馈到控制器中,以获得更高的跟踪控制精度。考虑风等外部干扰引起的扰动,当扰动存在上界时,文献 10 设计了扰动观测器。然而,上述研究忽略了实际HFV系统中攻角、航迹角等系统状态难以测量的问题,所设计的控制器难以工程实现。针对这一问题,文献 11 同时考虑无法测量的状态和干扰,将不确定性和外部干扰视为一个新的系统状态,与原有系统状态方程联立,设计扩张状态观测器(Extended state observer,ESO)同时对不可测系统状态、不确定和干扰构成的集总不确定进行估计。文献 12-13 在
11、此基础上进一步面向HFV控制系统提出了高增益观测器和变增益观测器,以提高ESO的收敛速度和估计精度。文献 14 研究了有限时间观测器,以实现更快的状态估计。值得注意的是,在上述工作中,HFV控制器都是在假设系统无噪声的情况下设计的。这种假设是不现实的。由于模型构建误差、极端飞行环境以及飞行器超高速飞行产生的电磁干扰等因素将影响传感器的测量精度,系统状态中必然存在测量噪声。如果将被噪声污染的系统状态直接应用到控制器设计中,会导致跟踪性能下降,甚至失控。因此,必须在控制器设计中降低噪声的影响。文献 15 利用二阶命令滤波器来处理系统噪声。文献 16 提出了线性高斯二次型(Linear quadra
12、tic Gaussian,LQG)/回路传递恢复(Loop transfer recovery,LTR)控制器,实现噪声影响下的HFV高精度全通道姿态控制。然而,上述HFV控制研究虽然在控制器设计中增加了噪声处理环节,却忽略了外部干扰的影响,导致控制系统抗噪能力增强,抗扰能力有限。同时考虑噪声和外部干扰的影响,文献 17 面向线性系统提出了ESO和Kalman滤波的组合观测策略,但是该方法不仅要求被控对象为线性系统,更要求传感器测量噪声的统计特性提前已知,这一点实际工程中难以保证,在呈现非线性且面临复杂多变飞行环境的HFV控制系统中更是完全不能实现。如何获取不受干扰和测量噪声影响的系统状态,并
13、基于该系统状态设计HFV控制器实现对高度参考信号和速度参考信号的高精度鲁棒跟踪仍待解决。针对HFV大包络飞行时,飞行器模型和传感器测量受风干扰和噪声影响,导致控制系统性能变差甚至失稳的问题,本文提出了一种基于状态观测的鲁棒控制策略。首先,采用小扰动线性化理论,推导HFV的线性纵向通道模型;其次,基于该模型设计“AKF+ESO”组合状态观测器,在未知测量噪声及外部干扰的影响下准确重构系统状态;再次,基于“AKF+ESO”估计的系统状态,设计LQG控制器实现HFV对高度参考信号和速度参考信号的高精度跟踪;最后,仿真验证本文HFV组合观测控制策略的有效性。1高超声速飞行器动力学模型1.1高超声速飞行
14、器的纵向通道模型考虑风场干扰的影响,HFV的纵向通道动力学模型为3-4:V=Tcos-Dm-gsinh=Vsin =L+TsinmV-gcosV+d(t)=q-q=MyyIyy+dq(t)(1)式中:V,h,q分别为速度、高度、航迹角、攻角、俯仰角速率;m为飞行器质量;g为重力加速度;Myy,Iyy分别为俯仰转动力矩和俯仰转动惯量;d(t),dq(t)为风场造成的外部干扰;T,D,L分别为飞行器204第 2 期张睿等:测量噪声下的高超声速飞行器组合观测鲁棒控制发动机推力、阻力、升力,并且满足:T=T()+T0()q S()C3T3+C2T2+CT+C0T+q S()C3T3+C2T2+CT+C
15、0TD q S()C2D2+CD+C0DL=L0+L q SC0L+q SCLMyy=MT+M0()+MeezTT+q Sc()C2M2+CM+C0M+q Sc CeMe(2)式中:,e分别为燃油当量、舵面偏角;q 为动压;S为参考面积;c 为平均气动弦;C()T,C()T,C()D,C()L,C()M为气动参数;zT为推力力矩耦合系数。1.2动力学模型线性化采用小扰动线性化理论线性化动力学模型(1)。在平衡点(0,u0)处进行泰勒展开,保留泰勒展开一次项,忽略高阶项,动力学模型(1)可改写为如下线性时不变状态方程,即:()t=A()t+Bu()t+Bd()ty()t=C()ty()t=C0(
16、)t(3)式中:(t)=V,h,qT为系统状态向量;u(t)=,eT为控制输入;(t)=d(t),dq(t)T为外部干扰;y(t)R5为测量输出;y(t)R2为受控输出。定义f1=V,f2=h,f3=,f4=,f5=q,fi为fi(i=1,2,3,4,5)对(t)的偏导数,ufi为fi对u(t)的偏导数,则A和B为偏导数矩阵在平衡点(0,u0)处的取值,有:A=f1,f2,f3,f4,f5T=0,u=u0B=uf1,uf2,uf3,uf4,uf5T=0,u=u0(4)具体形式为:A=00-gcoscos-Tsin-m0sin0Vcos00-LmV20gsinVL+sin+TcosmV0LmV2
17、0-gsinV-L+sin+TcosmV1000zT+qIyy0(5)B=T()cosm0T()sinmV T()sinmV00000MeT(6)式(5)和式(6)中,系统状态向量选取为平衡点处的系统状态0,系统输入选取为平衡点处的输入量u0,且:L=L+Tsin-mgcosq=q Sc()2C2M+CM=q S()3C3T2+2C2T+CT+q S()3C3T2+2C2T+CT=q S()2C2D+CD(7)模型(3)的其他参数为:Bd=001-1000001TC=I5 5C0=1000001000(8)假设采样时间为Ts,采用zero-order-hold(ZOH)方法对HFV线性动力学模
18、型(3)进行离散化处理。记k=(kTs),uk=u(kTs),k=(kTs),yk=y(kTs),y k=y(kTs),k为当前时刻,动力学模型(3)离散化后得:k+1=Ak+Buk+Bkyk=Cky k=C0k(9)式中:A=I5 5+ATs,B=BTs,B=BdTs。考虑测量噪声的影响,离散动力学模型(9)可重新写为:k+1=Ak+Buk+Bkyk=Ck+ky k=C0k(10)式中:k R5为测量噪声。记dk=d(kTs),dqk=dq(kTs),k=dk,dqkT。将k引 入 系 统 状 态,定 义 新 的 系 统 状 态 为k=Vk,hk,k,k,qk,dk,dqkT。HFV 的动力
19、学模型(10)可重新描述为如下扩张状态方程:k+1=Ak+Buk+Bkyk=Ck+k(11)式中:A=AB02 502 2,B=B02 2,B=05 2I2 2,C=205宇航学报第 45 卷C05 2,k=d(k)dq(k)。假设1.测量噪声、外部干扰及其导数有界。因此,存在未知正常数m,dm,dqm,dm,dqm使得k m,|dk|dm,|dqk|dqm,|dk|dm,|dqk|dqm。假设2.定义k为高斯白噪声,定义扰动估计误差为过程噪声vk=Bk-Bk,且vk为高斯白噪声,满足:Ek=0,EkTj=RkkjE vk=0,EvkvTj=QkkjETkvk=0(12)其中,kj=0,k j
20、1,k=j(13)式中:k为k的估计,由ESO(18)得到;k和j为两个计算时刻;Rk R5和Qk R5为正定矩阵,在实际系统中Rk和Qk通常未知。注1.线性动力学模型(3)中的矩阵A和B为非线性动力学函数fi在平衡点(0,u0)处的雅可比矩阵,均为时不变的常值矩阵。这种在平衡点处展开的HFV小扰动线性化方法在文献 16 中得到了详细论证。2基于AKF和ESO的状态组合观测考虑被噪声污染的测量信号如果直接用于HFV控制可能激发执行器额外的高频控制量,降低控制性能甚至导致系统失控,必须在控制器设计中使用不受噪声和干扰影响的系统状态。本节提出了“AKF+ESO”的状态组合观测方法,在统计特性未知的
21、测量噪声和外部干扰影响下获取准确的系统状态估计。面向HFV离散动力学模型中统计特性未知的测量噪声,设计AKF实时估计噪声的统计特性并获取不受噪声影响的系统状态,随后将降噪后的系统状态输入ESO实现对外部干扰的估计。值得注意的是,ESO估计的外部干扰又被引入了AKF,以提升AKF的抗干扰能力。这种自增强拓扑使组合观测结构能够在外部干扰、噪声并存的环境下精确地重构系统状态。考虑系统(10),设计AKF为:k+1/k=Ak+Buk+BkPk+1/k=APkAT+QkKk+1=Pk+1/kCT()CPk+1/kCT+Rk-1Pk+1=()I5 5-Kk+1C Pk+1/kZk+1=yk+1-Ck+1/
22、kk+1=k+1/k+Kk+1Zk+1(14)式中:k为k的估计值;Kk为 Kalman 滤波增益;k+1/k为状态一步预测;Pk+1/k为状态一步预测均方误差阵;Pk为状态预测均方误差阵;Zk为量测新息;Qk为Qk的估计;Rk为Rk的估计。Qk和Rk的更新律为:Qk=()1-kQk-1+k()KkZkZTkKTk-APk-1ATRk=()1-kRk-1+k()ZkZTk-CPk-1CT(15)式中:k=k-1/()k-1+b,k的初值设置为 0=1,b为遗忘因子且通常取值为b(0.95,0.999)。其中,Qk更新律中包含负定项-APk-1AT,Rk更新律中包含负定项-CPk-1CT。负定项
23、导致系统可能出现发散,为了保证系统稳定性,可舍弃负定项,损失部分状态估计精度,具体形式如下:Qk=()1-kQk-1+kKkZkZTkKTkRk=()1-kRk-1+kZkZTk(16)基于扩张状态方程(11),设计ESO为:k+1=Ak+Buk+Ks(yk-Ck)(17)式中:k为k的估计值;Ks R7 5为ESO增益。为了提升ESO的抗噪能力,将式(17)中的yk替换为 AKF估计预测的系统输出Ck,因此,式(17)可重新写为:k+1=Ak+Buk+Ks(Ck-Ck)(18)定义AKF估计误差为:k+1=k+1-k+1(19)定义ESO估计误差为:k+1=k+1-k+1(20)联立k+1和
24、k+1,可得“AKF+ESO”组合观测的系统状态误差为:k+1=k+1k+1=k+1k+1-k+1k+1(21)进一步计算可得:206第 2 期张睿等:测量噪声下的高超声速飞行器组合观测鲁棒控制k+1=Ak+Buk+BkAk+Buk+Bk-Ak+Buk+Bk+Kk+1()yk+1-Ck+1/kAk+Buk+Ks()Ck-Ck=Ak+BkAk+Bk-Kk+1()C()Ak+Buk+Bk+k+1-Ck+1Ks()Ck-Ck=Ak+BkAk+Bk-Kk+1C()Ak+Bk+Kk+1k+1KsC()k-k-KsC()k-k=Ak+BkAk+Bk-Kk+1C()Ak+Bk+Kk+1k+1-KsCk+Ks
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 测量 噪声 高超 声速 飞行器 组合 观测 鲁棒控制
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。