CN103823375A - 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 - Google Patents
纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 Download PDFInfo
- Publication number
- CN103823375A CN103823375A CN201410069966.0A CN201410069966A CN103823375A CN 103823375 A CN103823375 A CN 103823375A CN 201410069966 A CN201410069966 A CN 201410069966A CN 103823375 A CN103823375 A CN 103823375A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- sigma
- flight
- controller
- frequency
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000013461 design Methods 0.000 title claims abstract description 30
- 239000002131 composite material Substances 0.000 title abstract 2
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 230000003028 elevating effect Effects 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 5
- 238000012850 discrimination method Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009510 drug design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000013643 reference control Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种纵向飞行模型簇复合根轨迹多级PID鲁棒控制器设计方法,该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇;根据飞行包线内的幅频裕度和相位裕度军标要求,给出了对应根轨迹描述下的闭环极点分布限制指标,通过加入多级PID控制器并在飞行器全包线内的闭环极点分布限制指标和系统辨识中的模型辨识方法确定多级PID鲁棒控制器级数和参数值;从根轨迹描述下的闭环极点分布限制的概念出发设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。
Description
技术领域
本发明涉及一种飞行器控制器设计方法,特别涉及纵向飞行模型簇复合根轨迹多级PID鲁棒控制器设计方法,属于测控技术和飞行力学等范畴。
背景技术
飞行器起降过程的控制对飞行安全有重要作用;由于飞行器起降过程中飞行速度变化大,即使按照纵向模型也会面临强非线性问题;另一方面,飞行器的操纵舵存在饱和、死区等现象;从飞行安全考虑,超低空飞行(如飞机起飞/着陆)时,控制器必须保证系统具有一定的稳定裕度、无超调和平稳性,这样,就使得超低空飞行控制器设计非常复杂,不能直接套用现有控制理论进行飞行器控制的设计。
在现代实际飞行控制器的设计中,一少部分采用状态空间法进行设计,而大多数仍然采用以PID为代表的经典频域法和逆Nyquist阵列法为代表的现代频率法进行控制器设计。现代控制理论以状态空间法为特征、以解析计算为主要手段、以实现性能指标为最优的现代控制理论,而后有发展了最优控制方法、模型参考控制方法、自适应控制方法、动态逆控制方法,反馈线性化方法、直接非线性优化控制、变增益控制法、神经网络控制方法,模糊控制方法,鲁棒控制法以及多种方法组合控制等一系列控制器设计方法,发表的学术论文数以万计,例如2011年Ghasemi A设计了自适应模糊滑模控制的再入飞行器(GhasemiA,Moradi M,Menhaj M B.Adaptive Fuzzy Sliding Mode Control Design fora Low-Lift Reentry Vehicle[J].Journal of Aerospace Engineering,2011,25(2):210-216),2013年Babaei A R为非最小相位和非线性飞行器设计了模糊滑模控制自动驾驶仪(Babaei A R,Mortazavi M,Moradi M H.Fuzzy slidingmode autopilot design for nonminimum phase and nonlinear UAV[J].Journalof Intelligent and Fuzzy Systems,2013,24(3):499-509),很多研究仅仅停留在理想化的仿真研究阶段;而且这种设计存在三个问题:(1)由于无法进行飞行器超低空操纵稳定性试验,难以得到精确的被控对象的数学模型;(2)对于军标规定的稳定裕度等评价飞行控制系统的重要性能指标,状态空间法远不像经典频率法那样能以明显的形式表达出来;(3)控制器结构过于复杂、没有考虑实际控制器和飞行状态的约束,设计的控制器物理上不可实现。
英国的学者Rosenbrock系统地、开创性地研究了如何将频域法推广到多变量系统的设计中去,利用矩阵对角优势概念,把多变量问题转化为能用人们熟知的古典方法的单变量系统的设计问题,以后相继出现了Mayne序列回差法,MacFarlane特征轨迹法、Owens并矢展开法等方法,共同特点是把多输入一多输出、回路间严重关联的多变量系统的设计,化为一系列单变量系统的设计问题,进而可选用某一种古典方法(Nyquist和Bode的频率响应法,Evans的根轨迹法等)完成系统的设计,上述这些方法保留和继承了古典图形法的优点,不要求特别精确的数学模型,容易满足工程上的限制。特别是当采用有图形显示终端的人一机对话式的计算机辅助设计程序实现时,可以充分发挥设计者的经验和智慧,设计出既满足品质要求,又是物理上可实现的、结构简单的控制器;国内外对多变量频率法进行了改进研究(高大远,罗成,沈辉,胡德文,挠性卫星姿态解藕控制器多变量频率域设计方法,宇航学报,2007,Vol.28(2),pp442-447;熊柯,夏智勋,郭振云,倾斜转弯高超声速巡航飞行器多变量频域法解耦设计,弹箭与制导学报,2011,Vol.31(3),pp25-28)但是,这种设计方法可考虑系统不确定问题时保守性过大,在飞行器操纵舵限制情况下不能得到合理的设计结果。
综上所述,目前的控制方法还不能在飞行器模型变化、按照全飞行包线内的稳定裕度指标设计出超调量小、平稳的低空飞行控制器。
发明内容
为了克服现有方法不能在飞行器在全飞行包线内模型变化大的情况下设计出符合全飞行包线内的稳定裕度指标的超调量小、平稳低空飞行控制器的技术缺陷,本发明提供了一种纵向飞行模型簇复合根轨迹多级PID鲁棒控制器设计方法,该方法在给定不同高度、马赫数条件下通过扫频飞行试验直接确定获得全包线内的幅频和相频特性构成的模型簇;根据飞行包线内的幅频裕度和相位裕度军标要求,给出了对应根轨迹描述下的闭环极点分布限制指标,通过加入多级PID控制器并在飞行器全包线内的闭环极点分布限制指标和系统辨识中的模型辨识方法确定多级PID鲁棒控制器级数和参数值;从根轨迹描述下的闭环极点分布限制的概念出发设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。
本发明解决其技术问题所采用的技术方案:一种纵向飞行模型簇复合根轨迹多级PID鲁棒控制器设计方法,其特点是包括以下步骤:
1、给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的升降舵与飞行高度之模型簇,对应的飞行器升降舵与飞行高度之间开环传递函数簇描述为:
其中
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变化的系数簇,△k(s)为模型的不确定项;
2、候选多级PID控制器的传递函数为:
式中,kc为待确定的常数增益,N为整数,表示待确定的多级PID控制器的级数,kP(i)、kI(i)、kD(i)i=1,2,…,N为待确定的常数;
加入多级PID控制器后,整个系统的开环传递函数为:
对应的根轨迹方程为:
3、设s=σ+jω,其中:σ为s的实部,ω为s的虚部,j为虚部符号;系统的稳定裕度指标设定为:σ≤-ζ2,其中,ζ为非零实数,ξ给定数;根据飞行试验或风洞试验建立模型不确定项造成的滞后相角弧度,幅值系统的稳定裕度指标调整为:和其中,△M和△a均为整实数;
这样,系统的稳定裕度指标可以转化为:根据
或 所得到的根轨迹必须满足和根据该指标和极大似然准则或其它准则共同约束下,可以根据系统模型结构辨识中的极大似然方法或辨识方法确定多级PID控制器的级数N、常数kP(i)、kI(i)、kD(i)i=1,2,…,N。
本发明的有益效果是:从根轨迹描述下的闭环极点分布限制的概念出发,通过加入多级PID控制器,在全飞行包线内按照符合给定闭环极点分布限制要求和模型辨识方法确定多级PID鲁棒控制器的参数,设计出符合全飞行包线的超调量小、平稳的低空飞行控制器。
下面结合实施例对本发明作详细说明。
具体实施方式
1、给定不同高度、马赫数下使用线性扫频信号(f0为起始频率,f1为截止频率,r=(f1-f0)/T,T为扫频时间)或对数扫频信号f(t)=A(t)sin{2πf0/r·[exp(rt)-1]}(f0为起始频率,f1为截止频率,r=ln(f1/f0)/T,T为扫频时间)对飞机激励,可直接得到允许飞行的全包线内的幅频和相频特性,构成飞行器全包线内的升降舵与飞行高度之模型簇,对应的飞行器升降舵与飞行高度之间开环传递函数簇描述为:
其中
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变化的系数簇,△k(s)为模型的不确定项;
2、候选多级PID控制器的传递函数为:
式中,kc为待确定的常数增益,N为整数,表示待确定的多级PID控制器的级数,kP(i)、kI(i)、kD(i)i=1,2,…,N为待确定的常数;
加入多级PID控制器后,整个系统的开环传递函数为:
对应的根轨迹方程为:
3、设s=σ+jω,其中:σ为s的实部,ω为s的虚部,j为虚部符号;系统的稳定裕度指标设定为:σ≤-ζ2,其中,ζ为非零实数,ξ给定数;根据飞行试验或风洞试验建立模型不确定项造成的滞后相角弧度,幅值系统的稳定裕度指标调整为:和其中,△M和△a均为整实数;
这样,系统的稳定裕度指标可以转化为:根据
或
所得到的根轨迹必须满足和根据该指标和极大似然准则或其它准则共同约束下,可以根据系统模型结构辨识中的极大似然方法或辨识方法确定多级PID控制器的级数N、常数kP(i)、kI(i)、kD(i)i=1,2,…,N。
Claims (1)
1.一种纵向飞行模型簇复合根轨迹多级PID鲁棒控制器设计方法,其特点是包括以下步骤:
1)给定不同高度、马赫数下通过扫频飞行试验直接由允许飞行的全包线内的幅频和相频特性构成飞行器全包线内的升降舵与飞行高度之模型簇,对应的飞行器升降舵与飞行高度之间开环传递函数簇描述为:
其中
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变化的系数簇,△k(s)为模型的不确定项;
2)候选多级PID控制器的传递函数为:
式中,kc为待确定的常数增益,N为整数,表示待确定的多级PID控制器的级数,kP(i)、kI(i)、kD(i)i=1,2,…,N为待确定的常数;
加入多级PID控制器后,整个系统的开环传递函数为:
对应的根轨迹方程为:
3)设s=σ+jω,其中:σ为s的实部,ω为s的虚部,j为虚部符号;系统的稳定裕度指标设定为:σ≤-ζ2,其中,ζ为非零实数,ξ给定数;根据飞行试验或风洞试验建立模型不确定项造成的滞后相角弧度,幅值系统的稳定裕度指标调整为:M和其中,△M和△a均为整实数;
这样,系统的稳定裕度指标可以转化为:根据
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410069966.0A CN103823375B (zh) | 2014-02-28 | 2014-02-28 | 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410069966.0A CN103823375B (zh) | 2014-02-28 | 2014-02-28 | 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103823375A true CN103823375A (zh) | 2014-05-28 |
CN103823375B CN103823375B (zh) | 2016-04-06 |
Family
ID=50758493
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410069966.0A Expired - Fee Related CN103823375B (zh) | 2014-02-28 | 2014-02-28 | 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103823375B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109240327A (zh) * | 2018-09-11 | 2019-01-18 | 陕西千山航空电子有限责任公司 | 一种固定翼飞机飞行阶段识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5396415A (en) * | 1992-01-31 | 1995-03-07 | Honeywell Inc. | Neruo-pid controller |
JP3493340B2 (ja) * | 2000-12-20 | 2004-02-03 | 川崎重工業株式会社 | 自動制御システムの制御則決定方法 |
CN102176118A (zh) * | 2011-01-27 | 2011-09-07 | 西北工业大学 | 多输入-多输出飞行器等价稳定裕度鲁棒确定方法 |
CN102279564A (zh) * | 2011-04-29 | 2011-12-14 | 南京航空航天大学 | 应用智能pid控制器的飞行仿真转台控制系统及方法 |
US20120035747A1 (en) * | 2002-04-18 | 2012-02-09 | Cleveland State University | Scaling and parameterizing a controller |
CN102722176A (zh) * | 2012-06-18 | 2012-10-10 | 中国航天空气动力技术研究院 | 一种可变形无人机飞行控制方法 |
-
2014
- 2014-02-28 CN CN201410069966.0A patent/CN103823375B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5396415A (en) * | 1992-01-31 | 1995-03-07 | Honeywell Inc. | Neruo-pid controller |
JP3493340B2 (ja) * | 2000-12-20 | 2004-02-03 | 川崎重工業株式会社 | 自動制御システムの制御則決定方法 |
US20120035747A1 (en) * | 2002-04-18 | 2012-02-09 | Cleveland State University | Scaling and parameterizing a controller |
CN102176118A (zh) * | 2011-01-27 | 2011-09-07 | 西北工业大学 | 多输入-多输出飞行器等价稳定裕度鲁棒确定方法 |
CN102279564A (zh) * | 2011-04-29 | 2011-12-14 | 南京航空航天大学 | 应用智能pid控制器的飞行仿真转台控制系统及方法 |
CN102722176A (zh) * | 2012-06-18 | 2012-10-10 | 中国航天空气动力技术研究院 | 一种可变形无人机飞行控制方法 |
Non-Patent Citations (1)
Title |
---|
史忠科: "飞行器模型簇描述及鲁棒控制器设计", 《控制与决策》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109240327A (zh) * | 2018-09-11 | 2019-01-18 | 陕西千山航空电子有限责任公司 | 一种固定翼飞机飞行阶段识别方法 |
CN109240327B (zh) * | 2018-09-11 | 2021-10-12 | 陕西千山航空电子有限责任公司 | 一种固定翼飞机飞行阶段识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103823375B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103809446A (zh) | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN103809433A (zh) | 飞行器多回路模型簇复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809448A (zh) | 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 | |
CN103792848B (zh) | 纵向飞行模型簇人机闭环复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809434B (zh) | 纵向飞行模型簇复合根轨迹多级pid控制器设计方法 | |
CN103809442A (zh) | 飞行器多回路模型簇复合频率鲁棒控制器设计方法 | |
CN103777523A (zh) | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 | |
CN103823375A (zh) | 纵向飞行模型簇复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809452A (zh) | 纵向飞行模型簇颤振抑制复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103823376B (zh) | 纵向飞行模型簇复合pid控制器设计方法 | |
CN103823364A (zh) | 飞行器多回路模型簇复合根轨迹补偿鲁棒控制器设计方法 | |
CN103823374A (zh) | 飞行器多回路模型簇复合根轨迹补偿控制器设计方法 | |
CN103823378A (zh) | 纵向飞行模型簇颤振抑制复合pid鲁棒控制器设计方法 | |
CN103809449A (zh) | 飞行器多回路模型簇颤振抑制复合pid鲁棒控制器设计方法 | |
CN103853049A (zh) | 纵向飞行模型簇复合频率鲁棒控制器设计方法 | |
CN103823367A (zh) | 纵向飞行模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN103809453A (zh) | 纵向飞行模型簇人机闭环复合根轨迹补偿鲁棒控制器设计方法 | |
CN103809445A (zh) | 飞行器多回路模型簇复合pid控制器设计方法 | |
CN103823377A (zh) | 纵向飞行模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 | |
CN103823366A (zh) | 纵向飞行模型簇复合根轨迹补偿控制器设计方法 | |
CN103809455A (zh) | 纵向飞行模型簇复合根轨迹补偿鲁棒控制器设计方法 | |
CN103823365A (zh) | 纵向飞行模型簇复合pid鲁棒控制器设计方法 | |
CN103809444A (zh) | 飞行器多回路模型簇人机闭环pid鲁棒控制器设计方法 | |
CN103809443A (zh) | 飞行器多回路模型簇人机闭环复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809457B (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: 20160406 |
|
CF01 | Termination of patent right due to non-payment of annual fee |