飞行器多回路模型簇复合根轨迹补偿控制器设计方法技术

本发明专利技术提供了一种飞行器多回路模型簇复合根轨迹补偿控制器设计方法，该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇；根据飞行包线内的幅频裕度和相位裕度军标要求，给出了对应根轨迹描述下的闭环极点分布限制指标，通过加入多级串联滞后—超前补偿控制器并在飞行器全包线内的闭环极点分布限制指标和系统辨识中的模型辨识方法确定多级串联滞后—超前补偿控制器级数和参数值；从根轨迹描述下的闭环极点分布限制的概念出发设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。

【技术实现步骤摘要】
【专利摘要】本专利技术提供了一种，该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇；根据飞行包线内的幅频裕度和相位裕度军标要求，给出了对应根轨迹描述下的闭环极点分布限制指标，通过加入多级串联滞后—超前补偿控制器并在飞行器全包线内的闭环极点分布限制指标和系统辨识中的模型辨识方法确定多级串联滞后—超前补偿控制器级数和参数值；从根轨迹描述下的闭环极点分布限制的概念出发设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。【专利说明】
本专利技术涉及一种飞行器控制器设计方法，特别涉及，属于测控技术和飞行力学等范畴。
技术介绍
飞行器起降过程的控制对飞行安全有重要作用；由于飞行器起降过程中飞行速度变化大，即使按照纵向模型也会面临强非线性问题；另一方面，飞行器的操纵舵存在饱和、死区等现象；从飞行安全考虑，超低空飞行(如飞机起飞/着陆)时，控制器必须保证系统具有一定的稳定裕度、无超调和平稳性，这样，就使得超低空飞行控制器设计非常复杂，不能直接套用现有控制理论进行飞行器控制的设计。在现代实际飞行控制器的设计中，一少部分采用状态空间法进行设计，而大多数仍然采用以PID为代表的经典频域法和逆Nyquist阵列法为代表的现代频率法进行控制器设计。现代控制理论以状态空间法为特征、以解析计算为主要手段、以实现性能指标为最优的现代控制理论，而后有发展了最优控制方法、模型参考控制方法、自适应控制方法、动态逆控制方法，反馈线性化方法、直接非线性优化控制、变增益控制法、神经网络控制方法，模糊控制方法，鲁棒控制法以及多种方法组合控制等一系列控制器设计方法，发表的学术论文数以万计，例如2011年Ghasemi A设计了自适应模糊滑模控制的再入飞行器(GhasemiA，Moradi Mj Menhaj M B.Adaptive Fuzzy Sliding Mode Control Design for a Low-LiftReentry Vehicle.Journal of Aerospace Engineering, 2011，25 (2):210-216)，2013年Babaei A R为非最小相位和非线性飞行器设计了模糊滑模控制自动驾驶仪(Babaei A R, Mortazavi M, Moradi M H.Fuzzy sliding mode autopilot designfor nonminimum phase and nonlinear UAV.Journal of Intelligent and FuzzySystems, 2013，24(3):499_509)，很多研究仅仅停留在理想化的仿真研究阶段；而且这种设计存在三个问题:(I)由于无法进行飞行器超低空操纵稳定性试验，难以得到精确的被控对象的数学模型；(2)对于军标规定的稳定裕度等评价飞行控制系统的重要性能指标，状态空间法远不像经典频率法那样能以明显的形式表达出来；(3)控制器结构过于复杂、没有考虑实际控制器和飞行状态的约束，设计的控制器物理上不可实现。英国的学者Rosenbrock系统地、开创性地研究了如何将频域法推广到多变量系统的设计中去，利用矩阵对角优势概念，把多变量问题转化为能用人们熟知的古典方法的单变量系统的设计问题，以后相继出现了 Mayne序列回差法，MacFarlane特征轨迹法、Owens并矢展开法等方法，共同特点是把多输入一多输出、回路间严重关联的多变量系统的设计，化为一系列单变量系统的设计问题，进而可选用某一种古典方法(Nyquist和Bode的频率响应法，Evans的根轨迹法等)完成系统的设计，上述这些方法保留和继承了古典图形法的优点，不要求特别精确的数学模型，容易满足工程上的限制。特别是当采用有图形显示终端的人一机对话式的计算机辅助设计程序实现时，可以充分发挥设计者的经验和智慧，设计出既满足品质要求，又是物理上可实现的、结构简单的控制器；国内外对多变量频率法进行了改进研究(高大远，罗成，沈辉，胡德文，挠性卫星姿态解藕控制器多变量频率域设计方法，宇航学报，2007，Vol.28(2)，pp442-447 ;熊柯，夏智勋，郭振云，倾斜转弯高超声速巡航飞行器多变量频域法解耦设计，弹箭与制导学报，2011，Vol.31(3)，pp25-28)但是，这种设计方法可考虑系统不确定问题时保守性过大，在飞行器操纵舵限制情况下不能得到合理的设计结果。综上所述，目前的控制方法还不能在飞行器模型变化、按照全飞行包线内的稳定裕度指标设计出超调量小、平稳的低空飞行控制器。
技术实现思路
为了克服现有方法不能在飞行器在全飞行包线内模型变化大的情况下设计出符合全飞行包线内的稳定裕度指标的超调量小、平稳低空飞行控制器的技术缺陷，本专利技术提供了一种，该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇；根据飞行包线内的幅频裕度和相位裕度军标要求，给出了对应根轨迹描述下的闭环极点分布限制指标，通过加入多级串联滞后一超前补偿控制器并在飞行器全包线内的闭环极点分布限制指标和系统辨识中的模型辨识方法确定多级串联滞后一超前补偿控制器级数和参数值；从根轨迹描述下的闭环极点分布限制的概念出发设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。 本专利技术解决其技术问题所采用的技术方案:一种，其特点是包括以下步骤:步骤1、给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的操纵舵面与飞行高度之模型簇，并且能够跨越飞行包线获得飞行器的颤振频率，得到对应的飞行器操纵舵面与飞行高度之间开环传递函数模型簇矩阵为:【权利要求】1.一种，其特点是包括以下步骤: 步骤1、给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的操纵舵面与飞行高度之模型簇，并且能够跨越飞行包线获得飞行器的颤振频率，得到对应的飞行器操纵舵面与飞行高度之间开环传递函数模型簇矩阵为: 【文档编号】G05B13/04GK103823374SQ201410069798【公开日】2014年5月28日 申请日期:2014年2月28日 优先权日:2014年2月28日 【专利技术者】史忠科 申请人:西安费斯达自动化工程有限公司本文档来自技高网...

【技术保护点】
一种飞行器多回路模型簇复合根轨迹补偿控制器设计方法，其特点是包括以下步骤：步骤1、给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的操纵舵面与飞行高度之模型簇，并且能够跨越飞行包线获得飞行器的颤振频率，得到对应的飞行器操纵舵面与飞行高度之间开环传递函数模型簇矩阵为：其中，G为m×m方阵，m>1为正整数，s为拉普拉氏变换的自变量，h为飞行器飞行高度，M为马赫数，P为m×m单模方阵，D为m×m多项式对角矩阵，Q为m×m单模方阵，为多项式，n>1为正整数；选取满足条件：以及其中，GE为m×m方阵，PE为m×m单模方阵，DE为m×m多项式对角矩阵，di,E为DE的第IE行、第IE列元素，为D的第IE行、第IE列元素，IE=1,2,…,m,QE为m×m单模方阵，为多项式，arg为相角数学符号；飞行器多回路系统的控制器设为：GCA(s)=QE-1(s)Ga0(s)PE-1(s)]]>其中，GCA(s)为m×m方阵，Ga0(s)=diag[Gc,1(s),Gc,2(s),…,Gc,m(s)]为m×m对角矩阵；为Ga0(s)的第IE行、第IE列元素，IE=1,2,…,m；步骤2、控制器IE=1,2,…,m的设计过程如下：(1)令具体表达形式为：GIE(s)=e-σ(h,M)sK(h,M)A(h,M,s)B(h,M,s)]]>其中A(h,M,s)=sm+am?1(h,M)sm?1+am?2(h,M)sm?2+…+a1(h,M)s+a0(h,M)、B(h,M,s)=sn+bn?1(h,M)sn?1+bn?2(h,M)sn?2+…+b1(h,M)s+b0(h,M)为多项式，s为传递函数中常用的拉普拉斯变化后的变量，h,M分别为飞行高度和马赫数，σ(h,M)是俯仰回路的延迟时间，K(h,M)为随h,M变化的增益，al(h,M),l=0,1,2,…,m?1为多项式A(h,M,s)中随h,M变化的系数簇，bi(h,M),i=0,1,2,…,n?1为多项式B(h,M,s)中随h,M变化的系数簇；(2)候选多级串联滞后—超前补偿环节的传递函数为：Gc,IE(s)=kcΠi=1NTD1(i)s+1a(i)TL1(i)s+1·TD2(i)s/a(i)+1TL2(i)s+1]]>式中，kc为待确定的常数增益，N为整数，表示待确定的滞后—超前补偿环节的级数，TD1(i)、TL1(i)、TD2(i)、TL2(i),i=1,2,…,N为待确定的时间常数，a(i)>1,i=1,2,…,N为待确定的参数；加入多级串联滞后—超前补偿环节后，整个系统的开环传递函数为：GIE(s)Gc,IE(s)=kce-σ(h,M)sK(h,M)A(h,M,s)B(h,M,s)Πi=1NTD1(i)s+1a(i)TL1(i)s+1·TD2(i)s/a(i)+1TL2(i)s+1]]>对应的根轨迹方程为：e-σ(h,M)sK(h,M)A(h,M,s)Πi=1N{[TD1(i)s+1]·[TD2(i)s/a(i)+1}}+Πi=1N{[a(i)TL1(i)s+1]·[TL2(i)s+1]}·B(h,M,s)=0]]>(3)设s=σ+jω，其中：σ为s的实部，ω为s的虚部，j为虚部符号；系统的稳定裕度指标设定为：σ≤?ζ2，其中，ζ为非零实数，ξ给定数；这样，系统的稳定裕度指标可以转化为：根据{e-σ(h,M)sK(h,M)A(h,M,s)Πi=1N{[TD1(i)s+1]·[TD2(i)s/a(i)+1}}+Πi=1N{[a(i)TL1(i)s+1]·[TL2(i)s+1]}·B(h,M,s)}s=σ+jω=0]]>或Re{{e-σ(h,M)sK(h,M)A(h,M,s)Πi=1N{[TD1(i)s+1]·[TD2(i)s/a(i)+1}}+Πi=1N{[a(i)TL1(i)s+1]·[TL2(i)s+1]}·B(h,M,s)}s=σ+jω=0Im{{e-σ(h,M)sK(h,M)A(h,M,s)Πi=1N{[TD1(i)s+1]&CenterDot...

【技术特征摘要】

【专利技术属性】
技术研发人员：史忠科，
申请(专利权)人：西安费斯达自动化工程有限公司，
类型：发明
国别省市：