CN105094114B - 一种确定lpv控制器的稳定性的方法 - Google Patents
一种确定lpv控制器的稳定性的方法 Download PDFInfo
- Publication number
- CN105094114B CN105094114B CN201510350455.0A CN201510350455A CN105094114B CN 105094114 B CN105094114 B CN 105094114B CN 201510350455 A CN201510350455 A CN 201510350455A CN 105094114 B CN105094114 B CN 105094114B
- Authority
- CN
- China
- Prior art keywords
- theta
- lpv
- centerdot
- sigma
- controller
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
- G05B23/0254—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model based on a quantitative model, e.g. mathematical relationships between inputs and outputs; functions: observer, Kalman filter, residual calculation, Neural Networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种确定LPV控制器的稳定性的方法,该方法包括:根据系统建模误差和LPV参数在线测量偏差,通过模型转换得到参数不确定的LPV系统的线性系统的控制器求解问题的标准形式;将线性系统的控制器求解问题转化为求解一个线性正矩阵不等式的凸优化问题;当所述线性正矩阵不等式有解时,LPV控制器存在且稳定。通过使用本发明所提供的方法,可以确定LPV控制器的稳定性,得到控制器存在的充分条件。
Description
技术领域
本发明涉及航空航天技术,特别涉及一种确定线性参数变化(LPV,LinearParameterVarying)控制器的稳定性的方法。
背景技术
高超声速飞行器的动力学特性具有高度的非线性,并且系统参数在大范围内剧烈变化,因此一般的线性控制系统显然不能满足系统的性能要求;另外,由于计算的复杂性等不可避免的缺点,现有的非线性控制系统设计方法也不适于工程实现,因此,对于高超声速飞行器这类系统一个最常用的工程设计方法就是变增益控制。而LPV理论的出现,则弥补了传统变增益控制的不足,而且LPV系统的稳定性可以从理论上得到证明。
在现有技术中,传统的三回路自动驾驶仪通常采用的方法是增益调度技术,在工作点附近用线性时不变系统近似非线性系统,根据高超声速飞行器的飞行参数,利用插值的形式自适应地调整控制器参数。上述的现有技术虽然能够保证系统在每个工作点上满足要求的性能,但却无法从理论上保证系统在整个参数区域内的控制性能和稳定性。而高超声速飞行器在飞行过程中,其高度、速度、攻角、动压等参数将随着飞行轨迹的变化而发生剧烈的变化,系统状态距离线性化的平衡点较远,系统参数变化较快,外界扰动强,因此传统的增益预置方法的控制器存在明显的缺陷。
发明内容
有鉴于此,本发明提供一种确定LPV控制器的稳定性的方法,从而可以确定LPV控制器的稳定性。
本发明的技术方案具体是这样实现的:
一种确定LPV控制器的稳定性的方法,该方法包括:
A、根据系统建模误差和LPV参数在线测量偏差,通过模型转换得到参数不确定的LPV系统的线性系统的控制器求解问题的标准形式;
B、将线性系统的控制器求解问题转化为求解一个线性正矩阵不等式的凸优化问题;
C、当所述线性正矩阵不等式有解时,LPV控制器存在且稳定。
较佳的,所述步骤A包括:
考虑参数不确定性,对系统建模误差和LPV参数在线测量误差中的可确定信息进行进一步的提取,得到如下形式的LPV系统:
其中,参变矩阵均为LPV参数向量q(t)的测量值的仿射函数,B1、D1、E1、E2、F1、F2均是确定的时不变矩阵;
用z表示表示所述LPV系统的可测量输出,则有:
因此,将所述LPV系统进一步写成:
与所述LPV系统相对应的控制器K表示为:
其中,是的函数;使得所述由LPV系统与控制器K组成的闭环系统满足如下条件:
1)闭环系统参数依赖稳定;
2)被控输出y和扰动w满足指标是γ的H∞性能指标:
令:
则对于相对应的第二LPV系统:
则设该第二LPV系统的控制器为K1:
根据 得控制器K:
Ak=Ak1-Bk1(I+D11Dk1)-1D11Gk1
Bk=Bk1(I+D11Dk1)-1。
Ck=(I+Dk1D11)-1Ck1
Dk=Dk1(I+D11Dk1)-1
较佳的,所述步骤B包括:
将由所述第二LPV系统和控制器K1组成的闭环系统表示为:
其中:
进一步,所述闭环系统等价于:
令与所述闭环系统对应的Lyapunov函数为V(x)=xTPx,满足系统稳定要求:
P>0且
得到第一不等式:
因为p=Δq,ΔTΔ≤I,即pTp≤qTq,得到第二不等式:
qTq-pTp≥0;
当所述第一不等式和第二不等式同时成立时,等价于存在常数δ>0使得以下第三不等式成立:
令P=δP,则所述第三不等式成立等价于:
即下述LMI不等式成立:
较佳的,所述步骤C包括:
在确定LPV控制器是否存在且稳定时,使用如下所述的定理1:
定理1、存在控制器K1使得第二LPV系统稳定的充要条件是:存在P>0使得由所述第二LPV系统和控制器K1组成的闭环系统满足所述LMI不等式。
然后,将上述的LMI不等式写成:
H0+VKTUT+UKVT<0;
其中:
因此,如果记U⊥,V⊥,使得成立,则:
其中:
应用上述的定理1可得:
定理2、由所述第二LPV系统和控制器K1组成的闭环系统稳定等价于下述不等式组成立:
其中,所述不等式组中的第一个LMI不等式由计算得到,第二个LMI不等式由计算得到,第三个LMI不等式保证P>0成立;
由于 因此,将T取为使得 的可逆矩阵,则有:
因此,定理2中控制器存在的条件等价为如下的第二不等式组:
当所述第二不等式组中的线性正矩阵不等式有解时,则LPV控制器存在且稳定。
如上可见,在本发明中的确定LPV控制器的稳定性的方法中,针对高超声速飞行器的参数不确定性LPV模型,将线性系统的控制器求解转化为求解一个线性正矩阵不等式的凸优化问题,而当该线性正矩阵不等式有解时,则确定控制器存在且稳定,从而可以确定LPV控制器的稳定性,得到了控制器存在的充分条件。
附图说明
图1为本发明实施例中的确定LPV控制器的稳定性的方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明进一步详细说明。
本实施例提供了一种确定LPV控制器的稳定性的方法,该方法适用于高超声速飞行器的LPV系统。
图1为本发明实施例中的确定LPV控制器的稳定性的方法的流程示意图。如图1所示,本发明实施例中的确定LPV控制器的稳定性的方法包括:。
步骤101,根据系统建模误差和LPV参数在线测量偏差,通过模型转换得到参数不确定的LPV系统的线性系统的控制器求解问题的标准形式。
在本发明的技术方案中,可以先考虑系统建模误差和LPV参数在线测量偏差,建立一个LPV模型(即LPV系统)。
例如,较佳的,在本发明的具体实施例中,上述的步骤101包括:
考虑参数不确定性,对系统建模误差和LPV参数在线测量误差中的可确定信息进行进一步的提取,可以得到如下形式的适用于高超声速飞行器的LPV系统:
其中,参变矩阵均为LPV参数向量q(t)的测量值的仿射函数,B1、D1、E1、E2、F1、F2均是确定的时不变矩阵。
用z表示表示上述LPV系统的可测量输出,则有:
因此,上述待研究的高超声速飞行器的LPV系统可进一步写成:
与上述LPV系统相对应的控制器K可以表示为:
其中,是的函数,但不一定是仿射函数。使得上述由LPV系统与控制器K组成的闭环系统满足如下条件:
1、闭环系统参数依赖稳定;
2、被控输出y和扰动w满足指标是γ的H∞性能指标:
令:
则对于相对应的第二LPV系统:
则可设该第二LPV系统的控制器为K1:
那么,根据可得控制器K:
Ak=Ak1-Bk1(I+D11Dk1)-1D11Ck1
Bk=Bk1(I+D11Dk1)-1(6)
Ck=(I+Dk1D11)-1Ck1
Dk=Dk1(I+D11Dk1)-1
也就是说,如果求出控制器K1那么就可以确定控制器K。
因此,通过上述的模型转换,即可得到参数不确定的LPV系统鲁棒控制问题(即线性系统的控制器求解问题)的标准形式。
步骤102,将线性系统的控制器求解问题转化为求解一个线性正矩阵不等式的凸优化问题。
在本发明的技术方案中,为了确定LPV控制器是否存在且稳定,可以先将上述线性系统的控制器求解问题转化为求解一个线性正矩阵不等式的凸优化问题。
例如,较佳的,在本发明的具体实施例中,可以先将由所述第二LPV系统和控制器K1组成的闭环系统表示为:
其中:
进一步,上述闭环系统可以等价于:
令与上述闭环系统对应的Lyapunov函数为V(x)=xTPx,满足系统稳定要求:
P>0且
即可得第一不等式:
又因为p=Δq,ΔTΔ≤I,即pTp≤qTq,也即可得第二不等式:
qTq-pTp≥0(10)
当上述第一不等式(9)和第二不等式(10)同时成立时,等价于存在常数δ>0使得以下第三不等式成立:
进一步令P=δP,则上述第三不等式成立等价于:
也即下述LMI不等式成立:
因此可知,可以将线性系统的控制器求解问题转化为求解一个线性正矩阵不等式的凸优化问题。
步骤103,当所述线性正矩阵不等式有解时,LPV控制器存在且稳定。
在本发明的技术方案中,可以使用多种具体的实施方式来实现上述的步骤103。以下将以其中的一种实现方式为例,对本发明的技术方案进行详细的介绍。
例如,较佳的,在本发明的具体实施例中,所述步骤103包括:
在确定LPV控制器是否存在且稳定时,可以使用如下所述的定理1:
定理1、存在控制器K1使得第二LPV系统稳定的充要条件是:存在P>0使得由所述第二LPV系统和控制器K1组成的闭环系统满足上述的LMI不等式(12)。
然后,可以将上述的LMI不等式(12)写成:
H0+VKTUT+UKVT<0(13)
其中:
因此,如果记U⊥,V⊥,使得成立,则:
其中:
应用上述的定理1可得:
定理2、由所述第二LPV系统和控制器K1组成的闭环系统稳定等价于下述不等式组成立:
其中,上述不等式组中的第一个LMI不等式由计算得到,第二个LMI不等式由计算得到,第三个LMI不等式保证P>0成立。
进一步分析可知: 因此可以把T取为使得 的可逆矩阵,则有:
因此,定理2中控制器存在的条件可以等价为如下所述的第二不等式组:
当上述公式(14)中的线性正矩阵不等式有解时,则LPV控制器存在且稳定,可以由MATLAB中的LMI工具箱判断控制器存在条件是否成立。
由此可知,当上述线性正矩阵不等式有解时,LPV控制器存在且稳定。
综上可知,在本发明中的确定LPV控制器的稳定性的方法中,针对高超声速飞行器的参数不确定性LPV模型,将线性系统的控制器求解转化为求解一个线性正矩阵不等式的凸优化问题,而当该线性正矩阵不等式有解时,则确定控制器存在且稳定,从而可以确定LPV控制器的稳定性,得到了控制器存在的充分条件。通过使用该方法,可以对参数测量存在误差和非线性系统转化成LPV系统过程中存在建模误差的高超声速飞行器进行系统建模,提高LPV模型的建模精确程度。另外,由于本发明中的上述方法具有一定的通用性,因此也可以应用于不同外形的高超声速飞行器。此外,在实际工程应用中,考虑参数不确定性的LPV系统建模结果对于此类飞行器姿态控制的设计与研究也具有显著的参考价值。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (4)
1.一种确定LPV控制器的稳定性的方法,其特征在于,该方法包括:
A、根据系统建模误差和LPV参数在线测量偏差,通过模型转换得到参数不确定的LPV系统的线性系统的控制器求解问题的标准形式;
B、将线性系统的控制器求解问题转化为求解一个线性正矩阵不等式的凸优化问题;
C、当所述线性正矩阵不等式有解时,LPV控制器存在且稳定。
2.根据权利要求1所述的方法,其特征在于,步骤A包括:
考虑参数不确定性,对系统建模误差和LPV参数在线测量偏差中的可确定信息进行进一步的提取,得到如下形式的LPV系统:
其中,参变矩阵均为LPV参数向量q(t)的测量值的仿射函数,B1、D1、E1、E2、F1、F2均是确定的时不变矩阵;
用z表示表示所述LPV系统的可测量输出,则有:
因此,将所述LPV系统进一步写成:
与所述LPV系统相对应的控制器K表示为:
其中,是的函数;使得由LPV系统与控制器K组成的闭环系统满足如下条件:
1)闭环系统参数依赖稳定;
2)被控输出y和扰动w满足指标是γ的H∞性能指标:
令:
则对于相对应的第二LPV系统:
则设该第二LPV系统的控制器为K1:
根据得控制器K:
Ak=Ak1-Bk1(I+D11Dk1)-1D11Ck1
Bk=Bk1(I+D11Dk1)-1。
Ck=(I+Dk1D11)-1Ck1
Dk=Dk1(I+D11Dk1)-1
3.根据权利要求2所述的方法,其特征在于,步骤B包括:
将由所述第二LPV系统和控制器K1组成的闭环系统表示为:
其中:
进一步,所述闭环系统等价于:
令与所述闭环系统对应的Lyapunov函数为V(x)=xTPx,满足系统稳定要求:
P>0且
得到第一不等式:
因为p=Δq,ΔTΔ≤I,即pTp≤qTq,得到第二不等式:
qTq-pTp≥0;
当所述第一不等式和第二不等式同时成立时,等价于存在常数δ>0使得以下第三不等式成立:
令P=δP,则所述第三不等式成立等价于:
即下述LMI不等式成立:
4.根据权利要求3所述的方法,其特征在于,步骤C包括:
在确定LPV控制器是否存在且稳定时,使用如下定理1:
定理1、存在控制器K1使得第二LPV系统稳定的充要条件是:存在P>0使得由所述第二LPV系统和控制器K1组成的闭环系统满足所述LMI不等式;
然后,将上述的LMI不等式写成:
H0+VKTUT+UKVT<0;
其中:
因此,如果记U⊥,V⊥,使得成立,则:
其中:
应用上述的定理1可得:
定理2、由所述第二LPV系统和控制器K1组成的闭环系统稳定等价于下述不等式组成立:
其中,所述不等式组中的第一个LMI不等式由计算得到,第二个LMI不等式由计算得到,第三个LMI不等式保证P>0成立;
由于因此,将T取为使得的可逆矩阵,则有:
因此,定理2中控制器存在的条件等价为如下的第二不等式组:
当所述第二不等式组中的线性正矩阵不等式有解时,则LPV控制器存在且稳定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510350455.0A CN105094114B (zh) | 2015-06-23 | 2015-06-23 | 一种确定lpv控制器的稳定性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510350455.0A CN105094114B (zh) | 2015-06-23 | 2015-06-23 | 一种确定lpv控制器的稳定性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105094114A CN105094114A (zh) | 2015-11-25 |
CN105094114B true CN105094114B (zh) | 2016-06-01 |
Family
ID=54574791
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510350455.0A Expired - Fee Related CN105094114B (zh) | 2015-06-23 | 2015-06-23 | 一种确定lpv控制器的稳定性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105094114B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109992004B (zh) * | 2019-05-08 | 2022-04-22 | 哈尔滨理工大学 | 一种lpv系统异步切换状态反馈控制器设计方法 |
CN110850844B (zh) * | 2019-11-11 | 2020-09-22 | 清华大学深圳国际研究生院 | 一种多维执行器广义最小可检测故障的计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1126163A1 (en) * | 2000-02-16 | 2001-08-22 | Turbowinds N.V./S.A. | Blade pitch angle control device for wind turbine |
CN102678452B (zh) * | 2012-05-22 | 2013-10-30 | 江南大学 | 基于lpv变增益的风力机被动容错控制方法 |
-
2015
- 2015-06-23 CN CN201510350455.0A patent/CN105094114B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105094114A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104635490B (zh) | 一种单出杆液压缸位置伺服系统的输出反馈控制方法 | |
CN103617816B (zh) | 反应堆堆芯功率分布的测量方法 | |
CN105975645A (zh) | 一种基于多步的含激波区域飞行器流场快速计算方法 | |
CN104360635A (zh) | 一种电机位置伺服系统的抗干扰控制方法 | |
CN106383443B (zh) | 抗干扰控制方法及系统 | |
CN104898428A (zh) | 一种基于干扰估计的电液伺服系统自适应鲁棒控制方法 | |
CN104898550A (zh) | 动态伺服系统基于滑模扩张状态观测器的复合控制方法 | |
CN105094114B (zh) | 一种确定lpv控制器的稳定性的方法 | |
CN105093932B (zh) | 一种确定lpv变增益控制器的鲁棒性的方法 | |
CN104794332B (zh) | 一种高层建筑风致响应分析模型的不确定性分析方法 | |
CN105893699A (zh) | 一种含混合不确定性参数的连杆机构鲁棒综合方法 | |
CN110457852A (zh) | 基于迭代法的综合回弹补偿方法 | |
CN104572575B (zh) | 一种特大变形动态网格生成方法 | |
Sirisha et al. | Numerical integration of singularly perturbed delay differential equations using exponential integrating factor | |
CN107153755A (zh) | 一种页岩气井数值模拟的求解方法 | |
Mulleners et al. | Impact of surface roughness on the turbulent wake flow of a turbine blade | |
CN105093933B (zh) | 一种确定lpv变增益控制器的方法 | |
Herty et al. | Feedback controls for continuous priority models in supply chain management | |
CN104991445A (zh) | 一种全局稳定的电机伺服系统自适应输出反馈鲁棒控制方法 | |
CN108052003A (zh) | 基于光电平台精确模型的自抗扰控制器设计系统 | |
CN110909511B (zh) | 一种无曲面体积分的无粘低速绕流数值模拟方法 | |
KR101562863B1 (ko) | 래티스 볼츠만 이론을 이용한 유체 유동 시뮬레이션 방법 및 이를 실현하기 위한 기록 매체 | |
CN111008492B (zh) | 一种基于无雅克比矩阵的高阶单元欧拉方程数值模拟方法 | |
Fleischmann et al. | A systematic LPV/LFR modelling approach optimized for linearised gain scheduling control synthesis | |
Yao et al. | Nonlinear modeling by assembling piecewise linear models |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160601 Termination date: 20200623 |