CN102929135B - 飞行器时间滞后时变模型的控制器设计方法 - Google Patents
飞行器时间滞后时变模型的控制器设计方法 Download PDFInfo
- Publication number
- CN102929135B CN102929135B CN201210380994.5A CN201210380994A CN102929135B CN 102929135 B CN102929135 B CN 102929135B CN 201210380994 A CN201210380994 A CN 201210380994A CN 102929135 B CN102929135 B CN 102929135B
- Authority
- CN
- China
- Prior art keywords
- matrix
- time
- controller
- aircraft
- formula
- 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.)
- Expired - Fee Related
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种飞行器时间滞后时变模型的控制器设计方法,用于解决现有的鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术问题。技术方案是给出时变滞后系统分段鲁棒稳定可解条件,直接利用线性系统状态反馈的闭环期望极点选择,并根据所有闭环期望极点的实部全部为负数的特点,给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员对风洞或飞行试验得到的含有时间滞后时变不确定性飞行器模型直接设计飞行控制器,解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
Description
技术领域
本发明涉及一种控制器设计方法,特别是涉及一种飞行器时间滞后时变模型的控制器设计方法。
背景技术
飞机鲁棒控制是目前国际航空界研究的重点课题之一,在高性能飞机控制器设计时,必须考虑鲁棒稳定性和鲁棒控制问题;实际飞行器模型是很复杂的未知模型结构的非线性微分方程式,为了描述这种复杂的非线性,人们通常采用风洞和飞行试验得到按离散数据描述的试验模型;为了减少风险并降低试验成本,通常按照不同高度、马赫数进行飞行机动试验,这样,描述飞行器试验模型的离散数据并不是很多,这种模型对静稳定性较好的飞行器很实用。然而,现代和未来的战斗机为了提高“机敏性”都放宽了对静态稳定性的限制,战斗机通常要求在开环临界稳定点附近工作;这样就要求飞行控制系统能良好地处理模型不确定性问题;在实际飞行控制系统设计中要考虑以下几个主要问题:(1)将试验得到离散数据用某一逼近模型来描述,模型中存在未建模动态;(2)风洞试验不能进行全尺寸模型自由飞、存在约束,飞行试验离散点选择、初始飞行状态、机动飞行的输入动作选择等不可能将所有的非线性充分激励,采用系统辨识所得模型存在各种各样的误差;(3)飞行环境与试验环境有区别,流场变化和干扰等使得实际气动力、力矩模型与试验模型有区别;(4)执行部件与控制元件存在制造容差,系统运行过程中也存在老化、磨损等现象,与飞行试验的结果不相同;(5)在实际工程问题中,需要控制器比较简单、可靠,通常需要对数学模型人为地进行简化,去掉一些复杂的因素;因此,在研究现代飞机的控制问题时,就必须考虑鲁棒性问题;特别是飞行器迎角、侧滑角测量和许多物理、化工过程中存在着不同程度的时间滞后不确定性,如果在系统的分析或设计过程忽略这些时间滞后,就可能出现错误的结果或引起系统的不稳定。
1980年后,国际上开展了多种不确定系统的控制理论研究,特别是由加拿大学者Zames提出的H-infinit理论,Zames认为,基于状态空间模型的LQG方法之所以鲁棒性不好,主要是因为用白噪声模型表示不确定的干扰是不现实的;因此,在假定干扰属于某一已知信号集的情况下,Zames提出用其相应灵敏度函数的范数作为指标,设计目标是在可能发生的最坏干扰下使系统的误差在这种范数意义下达到极小,从而将干扰抑制问题转化为求解使闭环系统稳定;从此,国内外很多学者展开了H-infinit控制方法研究;在航空界,该方法一直处于探索阶段,美国NASA,德国宇航研究院、荷兰等国都对鲁棒控制方法进行了研究,取得了很多仿真和实验结果;国内的航空院校也对飞机鲁棒控制方法进行了一系列的研究,如文献(史忠科、吴方向等,《鲁棒控制理论》,国防工业出版社,2003年1月;苏宏业.《鲁棒控制基础理论》,科学出版社,2010年10月)介绍,但这些结果与实际应用的距离还相差甚大,难以直接对实际飞行控制器进行设计并应用;特别是很多研究仅仅根据李雅普诺夫定理给出了不确定时间滞后系统鲁棒稳定性条件,但对于这些不等式解的存在条件等问题涉及较少,不能得到具体实现时间滞后鲁棒控制器设计步骤,没有解决直接设计鲁棒飞行控制器的技术问题。
发明内容
为了克服现有鲁棒控制理论缺乏设计步骤难以直接设计飞行控制器的技术不足,本发明提供一种飞行器时间滞后时变模型的控制器设计方法;该方法提供了实际时间滞后时变系统鲁棒稳定控制器的设计条件,直接利用线性系统状态反馈的闭环期望极点选择,并根据所有闭环期望极点的实部全部为负数的特点,给出了限定条件不等式直接设计反馈矩阵,可以对风洞或飞行试验得到的含有时间滞后时变不确定性飞行器模型直接设计飞行控制器,解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
本发明解决其技术问题所采用的技术方案是:一种飞行器时间滞后时变模型的控制器设计方法,其特点是包括以下步骤:
步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为:
式中,x(t)∈Rn,u(t)∈Rm分别为状态和输入向量,A0(t),Aτ(t),B(t)为已知常系数矩阵,τ为未知延迟时间,ΔA0(t)为系数矩阵未知部分;
按照已知A(t),B(t)的变化范围分类,即在不同时间段将A0(t),B0(t)表达成:
式中,A0i、B0i为已知的常数矩阵,ΔA0iΔB0i为未知矩阵,tij、tij为时间常数,r、p为正整数,i、j为下标,不同时间段的A(t),B(t)表达式形式相同;在时间段tij≤t<tij+Tij内,飞行控制器为:u(t)=-Kix(t)
式中,Ki为常数反馈矩阵;
带入(1)式中,有:
步骤二、选取(A0i-B0iK0i)的特征值各不相同且实部为负,设计反馈矩阵Ki使得满足条件:
该控制器使得
式中,Mi为线性变换矩阵,
σik,ωik(k=1,2,…,n)为实数,jωik(k=1,2,…,n)表示虚数,diag为对角矩阵符号,
ΔA0i-ΔB0iKi通常假设为ΔA0i-ΔB0iKi=HiFiWi,Hi,Wi均假设为已经矩阵,0<Fi≤I,I=diag[1,1,…,1]为单位阵。
本发明的有益效果是:通过本发明提供的时变滞后系统分段鲁棒稳定可解条件,直接利用线性系统状态反馈的闭环期望极点选择,并根据所有闭环期望极点的实部全部为负数的特点,给出了限定条件不等式直接设计反馈矩阵。使得本研究领域的工程技术人员对风洞或飞行试验得到的含有时间滞后时变不确定性飞行器模型直接设计飞行控制器,解决了当前研究只给出鲁棒稳定性不等式而无法直接设计飞行控制器的技术问题。
下面结合具体实施方式对本发明作详细说明。
具体实施方式
本发明飞行器时间滞后时变模型的控制器设计方法具体步骤如下:
1、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为:
式中,x(t)∈Rn,u(t)∈Rm分别为状态和输入向量,A0(t),Aτ(t),B(t)为已知常系数矩阵,τ为未知延迟时间,ΔA0(t)为系数矩阵未知部分;全文符号相同;
按照已知A(t),B(t)的变化范围分类,即在不同时间段将A0(t),B0(t)表达成:
式中,A0i、B0i为已知的常数矩阵,ΔA0iΔB0i为未知矩阵,tij、Tij为时间常数,r、p为正整数,i、j为下标,不同时间段的A(t),B(t)表达式形式可以相同;
在时间段tij≤t<tij+Tij内,飞行控制器为:u(t)=-Kix(t)
式中,Ki为常数反馈矩阵;
带入(1)式中,有:
2、选取(A0i-B0iK0i)的特征值各不相同且实部为负,设计反馈矩阵Ki使得满足条件:
该控制器使得
式中,Mi为线性变换矩阵,
σik,ωik(k=1,2,…,n)为实数,jωik(k=1,2,…,n)表示虚数,diag为对角矩阵符号,
ΔA0i-ΔB0iKi通常假设为ΔA0i-ΔB0iKi=HiFiWi,Hi,Wi均假设为已经矩阵,0<Fi≤I,I=diag[1,1,…,1]为单位阵;
取飞行器纵向运动状态变量为x=[qαθ]T,输入变量为u=δe,其中q为俯仰角速度,α为气流迎角,θ为俯仰角,δe为升降舵偏角;在时间段20≤t<100内,状态方程系数矩阵为:
不确定部分为:
0<F1≤I,0≤λ1<1,
选择闭环期望极点即A01-B01K1的特征值σ(A01-B01K1)=diag[-0.5,-1,-2],可得:
控制器为:K1=[-0.37941.5382-0.6231]。
Claims (1)
1.一种飞行器时间滞后时变模型的控制器设计方法,其特征在于包括以下步骤:
步骤一、在给定高度、马赫数条件下通过风洞或飞行试验得到含有时间滞后不确定性的飞行器时变模型为:
式中,x(t)∈Rn,u(t)∈Rm分别为状态和输入向量,A0(t),Aτ(t),B(t)为已知常系数矩阵,τ为未知延迟时间,ΔA0(t)为系数矩阵未知部分;
按照已知A(t),B(t)的变化范围分类,即在不同时间段将A0(t),B0(t)表达成:
tij≤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)式中,有:
步骤二、选取(A0i-B0iKi)的特征值各不相同且实部为负,设计反馈矩阵Ki使得满足条件:
该控制器使得
鲁棒稳定;
式中,Mi为线性变换矩阵,
σik,ωik为实数,jωik表示虚数,其中k=1,2,…,n;diag为对角矩阵符号,
ΔA0i-ΔB0iKi通常假设为ΔA0i-ΔB0iKi=HiFiWi,Hi,Wi均假设为已知矩阵,0<Fi≤I,I=diag[1,1,…,1]为单位阵;
取飞行器纵向运动状态变量为输入变量为u=δe,其中q为俯仰角速度,α为气流迎角,为俯仰角,δe为升降舵偏角;在时间段20≤t<100内,状态方程系数矩阵为:
不确定部分为:
0<F1≤I,0≤λ1<1,
选择闭环期望极点即A01-B01K1的特征值σ(A01-B01K1)=diag[-0.5,-1,-2],得:
控制器为:K1=[-0.3794 1.5382 -0.6231]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210380994.5A CN102929135B (zh) | 2012-10-10 | 2012-10-10 | 飞行器时间滞后时变模型的控制器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210380994.5A CN102929135B (zh) | 2012-10-10 | 2012-10-10 | 飞行器时间滞后时变模型的控制器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102929135A CN102929135A (zh) | 2013-02-13 |
CN102929135B true CN102929135B (zh) | 2015-05-27 |
Family
ID=47643964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210380994.5A Expired - Fee Related CN102929135B (zh) | 2012-10-10 | 2012-10-10 | 飞行器时间滞后时变模型的控制器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102929135B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108427430A (zh) * | 2018-03-30 | 2018-08-21 | 南京航空航天大学 | 基于网络控制的四旋翼飞行器控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102436437A (zh) * | 2011-11-17 | 2012-05-02 | 西北工业大学 | 基于角速度的飞行器极限飞行时四元数傅里埃近似输出方法 |
CN102495829A (zh) * | 2011-11-17 | 2012-06-13 | 西北工业大学 | 基于角速度的飞行器极限飞行时四元数沃尔什近似输出方法 |
CN102591212A (zh) * | 2012-03-01 | 2012-07-18 | 北京航空航天大学 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
CN102707624A (zh) * | 2012-05-31 | 2012-10-03 | 西北工业大学 | 基于飞行器常规模型的纵向控制器区域设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8359320B2 (en) * | 2007-10-31 | 2013-01-22 | At&T Intellectual Property I, Lp | Metadata repository and methods thereof |
-
2012
- 2012-10-10 CN CN201210380994.5A patent/CN102929135B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102436437A (zh) * | 2011-11-17 | 2012-05-02 | 西北工业大学 | 基于角速度的飞行器极限飞行时四元数傅里埃近似输出方法 |
CN102495829A (zh) * | 2011-11-17 | 2012-06-13 | 西北工业大学 | 基于角速度的飞行器极限飞行时四元数沃尔什近似输出方法 |
CN102591212A (zh) * | 2012-03-01 | 2012-07-18 | 北京航空航天大学 | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 |
CN102707624A (zh) * | 2012-05-31 | 2012-10-03 | 西北工业大学 | 基于飞行器常规模型的纵向控制器区域设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102929135A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102929142B (zh) | 飞行器不确定时变模型的控制器设计方法 | |
CN102929141B (zh) | 飞行器时间滞后时变模型逼近及控制器设计方法 | |
CN102929138B (zh) | 含有非线性的飞行器控制器设计方法 | |
CN102929140B (zh) | 飞行器时间滞后模型逼近及控制器设计方法 | |
CN102929139B (zh) | 飞行器多区间模型的控制器设计方法 | |
CN102929134A (zh) | 飞行器多个时间滞后模型的控制器设计方法 | |
CN102929130B (zh) | 一种鲁棒飞行控制器设计方法 | |
CN102929144B (zh) | 飞行器区间模型的控制器设计方法 | |
CN102929128B (zh) | 飞行器不确定模型的控制器设计方法 | |
CN102081349B (zh) | 多输入-多输出等价气动伺服弹性鲁棒稳定性的飞行试验确定方法 | |
CN101364085A (zh) | 非线性微分黄金分割自适应控制方法 | |
CN105676852A (zh) | 小型无人直升机无动力学模型结构自适应姿态控制方法 | |
Blonigan | Least squares shadowing for sensitivity analysis of large chaotic systems and fluid flows | |
CN102929143B (zh) | 飞行器时间滞后模型的控制器设计方法 | |
CN102929129B (zh) | 飞行器多个时间滞后时变模型的控制器设计方法 | |
CN102929135B (zh) | 飞行器时间滞后时变模型的控制器设计方法 | |
Li et al. | Nonlinear robust control of hypersonic aircrafts with interactions between flight dynamics and propulsion systems | |
CN102929145B (zh) | 飞行器多个时间滞后时变模型逼近及控制器设计方法 | |
CN102929136B (zh) | 飞行器多个时间滞后模型逼近及控制器设计方法 | |
CN102929137B (zh) | 含有分布参数项的飞行器控制器设计方法 | |
CN103777523B (zh) | 飞行器多回路模型簇复合pid鲁棒控制器设计方法 | |
Chang et al. | Lateral control for ultra-low altitude airdrop based on the L1 adaptive control augmentation | |
CN103809449B (zh) | 飞行器多回路模型簇颤振抑制复合pid鲁棒控制器设计方法 | |
CN103792848A (zh) | 纵向飞行模型簇人机闭环复合根轨迹多级pid鲁棒控制器设计方法 | |
CN103809433A (zh) | 飞行器多回路模型簇复合根轨迹多级pid鲁棒控制器设计方法 |
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: 20150527 Termination date: 20211010 |
|
CF01 | Termination of patent right due to non-payment of annual fee |