发明名称 飞行器多回路模型簇复合根轨迹多级PID鲁棒控制器设计方法
摘要 本发明提供了一种飞行器多回路模型簇复合根轨迹多级PID鲁棒控制器设计方法,该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇;根据飞行包线内的幅频裕度和相位裕度军标要求,给出了对应根轨迹描述下的闭环极点分布限制指标,通过加入多级PID控制器并在飞行器全包线内的闭环极点分布限制指标和系统辨识中的模型辨识方法确定多级PID鲁棒控制器级数和参数值;从根轨迹描述下的闭环极点分布限制的概念出发设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。
申请公布号 CN103809433A 申请公布日期 2014.05.21
申请号 CN201410070271.4 申请日期 2014.02.28
申请人 西安费斯达自动化工程有限公司 发明人 史忠科
分类号 G05B11/42(2006.01)I;G05D1/10(2006.01)I 主分类号 G05B11/42(2006.01)I
代理机构 代理人
主权项 1.一种飞行器多回路模型簇复合根轨迹多级PID鲁棒控制器设计方法,其特点是包括以下步骤: 步骤1、给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的操纵舵面与飞行高度之模型簇,并且能够跨越飞行包线获得飞行器的颤振频率,得到对应的飞行器操纵舵面与飞行高度之间开环传递函数模型簇矩阵为: <img file="FDA0000470855960000011.GIF" wi="1372" he="69" />其中,G为m×m方阵,m&gt;1为正整数,s为拉普拉氏变换的自变量,h为飞行器飞行高度,M为马赫数,Δ为不确定向量,P为m×m单模方阵,D为m×m多项式对角矩阵,Q为m×m单模方阵,<img file="FDA0000470855960000012.GIF" wi="45" he="57" />为多项式,n&gt;1为正整数;选取<img file="FDA0000470855960000013.GIF" wi="1209" he="71" />满足条件:<img file="FDA0000470855960000014.GIF" wi="1351" he="251" />以及<img file="FDA0000470855960000015.GIF" wi="1696" he="94" />其中,G<sub>E</sub>为m×m方阵,P<sub>E</sub>为m×m单模方阵,D<sub>E</sub>为m×m多项式对角矩阵,d<sub>i,E</sub>为D<sub>E</sub>的第I<sub>E</sub>行、第I<sub>E</sub>列元素,<img file="FDA0000470855960000019.GIF" wi="63" he="70" />为D的第I<sub>E</sub>行、第I<sub>E</sub>列元素,I<sub>E</sub>=1,2,…,m,Q<sub>E</sub>为m×m单模方阵,<img file="FDA0000470855960000018.GIF" wi="66" he="54" />为多项式,arg为相角数学符号;飞行器多回路系统的控制器设为: <img file="FDA0000470855960000016.GIF" wi="629" he="83" />其中,G<sub>CA</sub>(s)为m×m方阵,G<sub>a0</sub>(s)=diag[G<sub>c,1</sub>(s),G<sub>c,2</sub>(s),…,G<sub>c,m</sub>(s)]为m×m对角矩阵;<img file="FDA00004708559600000110.GIF" wi="198" he="87" />为G<sub>a0</sub>(s)的第I<sub>E</sub>行、第I<sub>E</sub>列元素,I<sub>E</sub>=1,2,…,m;步骤2、控制器<img file="FDA00004708559600000111.GIF" wi="227" he="92" />I<sub>E</sub>=1,2,…,m的设计过程如下:(1)令<img file="FDA0000470855960000017.GIF" wi="1047" he="94" />具体表达形式为:<img file="FDA0000470855960000021.GIF" wi="969" he="142" />其中 A(h,M,s)=s<sup>m</sup>+a<sub>m-1</sub>(h,M)s<sup>m-1</sup>+a<sub>m-2</sub>(h,M)s<sup>m-2</sup>+…+a<sub>1</sub>(h,M)s+a<sub>0</sub>(h,M)、 B(h,M,s)=s<sup>n</sup>+b<sub>n-1</sub>(h,M)s<sup>n-1</sup>+b<sub>n-2</sub>(h,M)s<sup>n-2</sup>+…+b<sub>1</sub>(h,M)s+b<sub>0</sub>(h,M)为多项式,s为传递函数中常用的拉普拉斯变化后的变量,h,M分别为飞行高度和马赫数,σ(h,M)是俯仰回路的延迟时间,K(h,M)为随h,M变化的增益,a<sub>l</sub>(h,M),l=0,1,2,…,m-1为多项式A(h,M,s)中随h,M变化的系数簇,b<sub>i</sub>(h,M),i=0,1,2,…,n-1为多项式B(h,M,s)中随h,M变化的系数簇,△k(s)为模型的不确定项; (2)PID控制器的传递函数为: <img file="FDA0000470855960000022.GIF" wi="798" he="152" />式中,k<sub>c</sub>为待确定的常数增益,N为整数,表示待确定的多级PID控制器的级数,k<sub>P</sub>(i)、k<sub>I</sub>(i)、k<sub>D</sub>(i)i=1,2,…,N为待确定的常数; 加入多级PID控制器后,整个系统的开环传递函数为: <img file="FDA0000470855960000023.GIF" wi="1698" he="156" />对应的根轨迹方程为: <img file="FDA0000470855960000024.GIF" wi="1762" he="144" />(3)其中:σ为s的实部,ω为s的虚部,j为虚部符号;系统的稳定裕度指标设定为:σ≤-ζ<sup>2</sup>,<img file="FDA0000470855960000025.GIF" wi="198" he="149" />其中,ζ为非零实数,ξ给定数;根据飞行试验或风洞试验建立模型不确定项造成的滞后相角<img file="FDA0000470855960000026.GIF" wi="410" he="82" />弧度,幅值<img file="FDA0000470855960000027.GIF" wi="352" he="89" />系统的稳定裕度指标调整为:<img file="FDA0000470855960000028.GIF" wi="284" he="82" />和<img file="FDA0000470855960000029.GIF" wi="378" he="136" />其中, △<sub>M</sub>和△<sub>a</sub>均为整实数;这样,系统的稳定裕度指标可以转化为:根据 <img file="FDA0000470855960000031.GIF" wi="1700" he="153" />或<img file="FDA0000470855960000032.GIF" wi="1824" he="324" />所得到的根轨迹必须满足<img file="FDA0000470855960000033.GIF" wi="284" he="82" />和<img file="FDA0000470855960000034.GIF" wi="380" he="148" />根据该指标和极大似然准则或其它准则共同约束下,可以根据系统模型结构辨识中的极大似然方法或辨识方法确定多级PID控制器的级数N、常数k<sub>P</sub>(i)、k<sub>I</sub>(i)、k<sub>D</sub>(i)i=1,2,…,N。
地址 710075 陕西省西安市高新区科技路金桥国际广场12101号