本发明专利技术公开了一种飞行器时间滞后时变模型逼近及控制器设计方法，用于解决现有的鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术问题。技术方案是给出时变逼近系统分段鲁棒稳定可解条件，直接利用线性系统状态反馈的闭环期望极点选择，并根据所有闭环期望极点的实部全部为负数的特点，给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员对风洞或飞行试验得到的含有时间滞后时变不确定性飞行器模型直接设计飞行控制器，解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
【技术实现步骤摘要】
本专利技术涉及一种控制器设计方法，特别是涉及一种。
技术介绍
飞机鲁棒控制是目前国际航空界研究的重点课题之一，在高性能飞机控制器设计时，必须考虑鲁棒稳定性和鲁棒控制问题；实际飞行器模型是很复杂的未知模型结构的非线性微分方程式，为了描述这种复杂的非线性，人们通常采用风洞和飞行试验得到按离散数据描述的试验模型；为了减少风险并降低试验成本，通常按照不同高度、马赫数进行飞行机动试验，这样，描述飞行器试验模型的离散数据并不是很多，这种模型对静稳定性较好的飞行器很实用。然而，现代和未来的战斗机为了提高“机敏性”都放宽了对静态稳定性的限制，战斗机通常要求在开环临界稳定点附近工作；这样就要求飞行控制系统能良好地处理·模型不确定性问题；在实际飞行控制系统设计中要考虑以下几个主要问题(I)将试验得到离散数据用某一逼近模型来描述，模型中存在未建模动态；(2)风洞试验不能进行全尺寸模型自由飞、存在约束，飞行试验离散点选择、初始飞行状态、机动飞行的输入动作选择等不可能将所有的非线性充分激励，采用系统辨识所得模型存在各种各样的误差；(3)飞行环境与试验环境有区别，流场变化和干扰等使得实际气动力、力矩模型与试验模型有区别；(4)执行部件与控制元件存在制造容差，系统运行过程中也存在老化、磨损等现象，与飞行试验的结果不相同；(5)在实际工程问题中，需要控制器比较简单、可靠，通常需要对数学模型人为地进行简化，去掉一些复杂的因素；因此，在研究现代飞机的控制问题时，就必须考虑鲁棒性问题；特别是飞行器迎角、侧滑角测量和许多物理、化工过程中存在着不同程度的时间滞后不确定性，如果在系统的分析或设计过程忽略这些时间滞后，就可能出现错误的结果或引起系统的不稳定。1980年后，国际上开展了多种不确定系统的控制理论研究，特别是由加拿大学者Zames提出的H-infinit理论，Zames认为,基于状态空间模型的LQG方法之所以鲁棒性不好，主要是因为用白噪声模型表示不确定的干扰是不现实的；因此，在假定干扰属于某一已知信号集的情况下，Zames提出用其相应灵敏度函数的范数作为指标，设计目标是在可能发生的最坏干扰下使系统的误差在这种范数意义下达到极小，从而将干扰抑制问题转化为求解使闭环系统稳定；从此，国内外很多学者展开了 H-infinit控制方法研究；在航空界，该方法一直处于探索阶段，美国NASA，德国宇航研究院、荷兰等国都对鲁棒控制方法进行了研究，取得了很多仿真和实验结果；国内的航空院校也对飞机鲁棒控制方法进行了一系列的研究，如文献(史忠科、吴方向等，《鲁棒控制理论》，国防工业出版社，2003年I月；苏宏业.《鲁棒控制基础理论》，科学出版社，2010年10月)介绍，但这些结果与实际应用的距离还相差甚大，难以直接对实际飞行控制器进行设计并应用；特别是很多研究仅仅根据李雅普诺夫定理给出了不确定时间滞后系统鲁棒稳定性条件，但对于这些不等式解的存在条件等问题涉及较少，不能得到具体实现时间滞后鲁棒控制器设计步骤，没有解决直接设计鲁棒飞行控制器的技术问题。
技术实现思路
为了克服现有鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术不足，本专利技术提供一种；该方法提供了实际时间滞后时变系统鲁棒稳定控制器的分段逼近设计条件，直接利用线性系统状态反馈的闭环期望极点选择，并根据所有闭环期望极点的实部全部为负数的特点，给出了限定条件不等式直接设计反馈矩阵，可以对风洞或飞行试验得到的含有时间滞后时变不确定性飞行器模型直接设计飞行控制器，解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。 本专利技术解决其技术问题所采用的技术方案是一种，其特点是包括以下步骤步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为x(t) = [A0(t) + AA0(0]x + Ar(t)x(t-τ) + B(0u(i)( I)式中，x(t) e Rn，u(t) e Rni分别为状态和输入向量，AQ(t)，AT (t)，B(t)为已知常系数矩阵，τ为未知延迟时间，AA0(t)为系数矩阵未知部分；全文符号相同；ΓΝ( ,τ) = Φ(/-τ,/,τ)联立j.pt、,迭代求解 N(t, τ ), [mm {xreal (t) - Φ( , 0, TMt0) — j。Φ0， , τ)Β{ )η{ ) }'式中Φ (t，t0, τ )为状态转移矩阵φ{ , η, τ) = [A0(O + Ar(/)N(r, Γ)]φ( , t0, T)Xreal⑴为实际系统状态响应；等价后系统的系数矩阵为A(t) =A0(t)+Aτ (t)N(t, τ ) + AA(t)式中:AA(t)为系数矩阵未知部分；按照已知A(t)，B(t)的变化范围分类，即在不同时间段将A(t)，B(t)表达成 Α(0 = A0i +AA0iR 丄=j = \,2,- P)[Β(0 = Β0 +ΔΒ0,式中，AQi、Btli为已知的常数矩阵，AAtli ΔΒΜ为未知矩阵，tu、Tij为时间常数，r、P为正整数，i、j为下标，不同时间段的A(t)，B(t)表达式形式可以相同；在时间段tij ( 内，飞行控制器为式中，Ki为常数反馈矩阵；带入(O式中，有i(/) = [(A , — B 0,.Κ ,) + (ΔΑ0 — AB0K f)]x(t) + ΔΑΓ(/)χ(/-τ)式中ΛΑτ (t)为系统等价后未知残余部分；步骤二、选取(Acii-BtliKtli)的特征值各不相同且实部为负，设计反馈矩阵Ki使得满足条件Ai >Μ[[(ΔΑ0,-ΔΒ() Κ )τΜτΜ1(ΔΑ0 -ΔΒ0 ΚΙ.)+ΔΑ (ΟΜ-τΜ1ΔΑΓ(Ο]Μ ；该控制器使得χ( O = [(A0f — B0fKi)+ (Δ A 0 — Δ B 0iK t )]χ (O + Δ A r (t)ii(t — τ)鲁棒稳定；式中，Mi为线性变换矩阵，Μ Γ1 (Αο,- - Β ο,κ ,-)Μ , = diag [σ + Jmn, σ + jo)i2,…，σin + jmin ],σ ik，C0ik(k=l,2,…，η)为实数，j ω ik(k=l, ζ, ···, η)表示虚数，diag 为对角矩阵符号，·Λ. = diag [σ,21； σ ···， σ1 ；Λ Atli-ΛBtliKi 通常假设为 Λ Atli-ΛBtliKi=HiFiWi,Hi,Wi 均假设为已经矩阵，(KFi ( I，I=diag[l, I,…，I]为单位阵。本专利技术的有益效果是通过本专利技术提供的时变逼近系统分段鲁棒稳定可解条件，直接利用线性系统状态反馈的闭环期望极点选择，并根据所有闭环期望极点的实部全部为负数的特点，给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员对风洞或飞行试验得到的含有时间滞后时变不确定性飞行器模型直接设计飞行控制器，解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。下面结合具体实施方式对本专利技术作详细说明。具体实施例方式本专利技术具体步骤如下I、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为x(t) = [A0 (/) + AA0(/本文档来自技高网...
【技术保护点】
一种飞行器时间滞后时变模型逼近及控制器设计方法，其特征在于包括以下步骤：步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为：x·(t)=[A0(t)+ΔA0(t)x+Aτ(t)x(t-τ)+B(t)u(t)---(1)式中，x(t)∈Rn，u(t)∈Rm分别为状态和输入向量，A0(t)，Aτ(t)，B(t)为已知常系数矩阵，τ为未知延迟时间，ΔA0(t)为系数矩阵未知部分；全文符号相同；联立迭代求解N(t，τ)，式中：Φ(t，t0，τ)为状态转移矩阵Φ·(t,t0,τ)=[A0(t)+Aτ(t)N(t,τ)]Φ(t,t0,τ)xreal(t)为实际系统状态响应；等价后系统的系数矩阵为A(t)＝A0(t)+Aτ(t)N(t，τ)+ΔA(t)式中：ΔA(t)为系数矩阵未知部分；按照已知A(t)，B(t)的变化范围分类，即在不同时间段将A(t)，B(t)表达成：A(t)=A0i+ΔA0iB(t)=B0i+ΔB0itij≤t<tij+Tij(i=1,2,...,r,j=1,2,...,p)式中，A0i、B0i为已知的常数矩阵，ΔA0i？ΔB0i为未知矩阵，tij、Tij为时间常数，r、p为正整数，i、j为下标，不同时间段的A(t)，B(t)表达式形式相同；在时间段tij≤t＜tij+Tij内，飞行控制器为：u(t)＝？Kix(t)式中，Ki为常数反馈矩阵；带入(1)式中，有：x·(t)=[(A0i-B0iKi)+(ΔA0i-ΔB0iKi)]x(t)+ΔAτ(t)x(t-τ)式中：ΔAτ(t)为系统等价后未知残余部分；步骤二、选取(A0i？B0iK0i)的特征值各不相同且实部为负，设计反馈矩阵Ki使得满足条件：Λi>MiT[(ΔA0i-ΔB0iKi)TMi-TMi-1(ΔA0i-ΔB0iKi)+ΔAτT(t)Mi-TMi-1ΔAτ(t)]Mi;该控制器使得x·(t)=[(A0i-B0iKi)+(ΔA0i-ΔB0iKi)]x(t)+ΔAτ(t)x(t-τ)鲁棒稳定；式中，Mi为线性变换矩阵，Mi-1(A0i-B0iKi)Mi=diag[σi1+jωi1,σi2+jωi2,...,σin+jωin],σik，ωik(k＝1，2，…，n)为实数，jωik(k＝1，2，…，n)表示虚数，diag为对角矩阵符号，Λi=diag[σi12,σi22,...,σin2];ΔA0i？ΔB0iKi通常假设为ΔA0i？ΔB0iKi＝HiFiWi，Hi，Wi均假设为已经矩阵，0＜Fi≤I，I＝diag[1，1，…，1]为单位阵。FDA00002234812200012.jpg...
【技术特征摘要】
1. 一种飞行器时间滞后时变模型逼近及控制器设计方法，其特征在于包括以下步骤步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为 x(0 = [A0(t) + AA0(t)]x+Ar-τ) + B(/)u(/)( I) 式中，X(t) e Rn,U(t) e Rni分别为状态和输入向量，AtlUhA,⑴，B(t)为已知常系数矩阵，τ为未知延迟时间，AAtlU)为系数矩阵未知部分；全文符号相同； .Ν(/,τ) = Φ( -τ,/,Γ) 联立j . f、，、(>, 、D,、( , , 17迭代求解 N(t, τ ), [mm {xreal (/) - Φ( , t0, TMt0) — j。Φ(Χ ζ, T)B(g)u(g)dg} ~ 式中Φ (t, t0, τ )为状态转移矩阵 Φ( ,/0, τ) = [A0(O + AT(t)N{t, γ)]Φ(/, /0, τ) XreaI (t)为实际系统状态响应； 等价后系统的系数矩阵为 A (t) = Α0( )+Ατ (t)N(t, τ ) + Δ A (t) 式中AA(t)为系数矩阵未知部分； 按照已知A(t)，B(t)的变化范围分类，即在不同时间段将A(t)，B(t)表达成 Γ A(0 = Ani + AAni Β(/) = Β0 .+ΔΒ0, 13 式中，A0i、B0i为已知的常数矩阵，ΔΑμ AB0i为未知矩阵，tu、Tij...
【专利技术属性】
技术研发人员：史忠科，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。