线性矩阵不等式(LMI)在控制理论中的应用

news2024/12/23 20:41:47

目录

(一)Matlab中的LMI处理工具包

 (二)为什么LMI成为控制理论领域重要工具?

(三)LMI在与Lyapunov不等式的关系

(1)线性矩阵不等式

 (2)线性矩阵不等式系统

(3)舒尔(Schur)补

(四)LMI中常见引理

引理2(广义KYP引理[4])

推论1(广义KYP引理推论[4])

引理3(射影定理[1])

引理4 (Jensen不等式[5,6]

引理5 (Finsler's引理[7]):

参考文献


(一)Matlab中的LMI处理工具包

        matlab中有专门求解线性矩阵不等式的工具包YALMIP,可以在官网上下载安装,可参考yalmip安装教程。yalmip只提供了一些基本的LMI求解方法,有更复杂的不等式求解需求可以再安装cplex工具包。关于怎么使用yalmip工具包可参考yalmip代码编写教程。

 (二)为什么LMI成为控制理论领域重要工具?

        线性矩阵不等式(LMI)技术是分析和综合控制系统的重要工具,尤其是鲁棒控制领域,主要因素有以下三个[1]:

  • 对于经典的控制方法,LMI技术的优点就是操作简单。在LMI技术出现以前,人们通过求解Ricaati方程来设计最优控制器,但是Ricaati方程求解很难。而LMI技术仅需少量的概念和基本原理就能开发出实用的工具(现在我们用YALMIP工具包就能很方便的求解LMI问题)。
  • LMI技术提供了控制问题的广泛前景,包括鲁棒性分析,标称H∞,H2和鲁棒控制综合,多目标综合,线性参数变化综合,其中一些无法在古典控制领域中解决处理。
  • LMI技术是利用凸优化的强有效的数值工具,并在理论体系上附加上有效的软件工具。

(三)LMI在与Lyapunov不等式的关系

(1)线性矩阵不等式

考虑线性矩阵不等式表达如下[2]:

上式中 ,x_1,....x_m为决策变量,特别地,在上述一般形式中的函数 为实对称矩阵。上式中的F(x)<0,表示矩阵F(x)是负定的,即对所有的非零向量 v\epsilon R^m, v^TF(x)v<0或者F(x)的最大特征值小于0。
        在许多系统与控制问题中,问题的变量以矩阵的形式出现,如Lyapunov矩阵不等式:

其中,矩阵 为已知的具有合适的维度的常数矩阵,且Q 为已知的对称矩阵,
为对称矩阵变量。假定 E_1,E_2,....,E_MR^n中的一组基,则对任意对称 存在x_1,x_2,....,x_M使得

         因此有:

        通过这种转化,就得到了一个更为一般的线性矩阵不等式的表达方式。

 (2)线性矩阵不等式系统

假设存在多个矩阵不等式:

 总体构成称线性矩阵不等式系统。引进F(x)=diag\left \{F_1(x),....F_k(x) \right \},则F_1(x)<0,...F_k(x)<0同时成立当且仅当F(x)<0。因此一个线性矩阵不等式系统也可以用一个单一的线性矩阵不等式表示。

(3)舒尔(Schur)补

        在许多一些非线性矩阵不等式转化成线性矩阵不等式的问题中,我们常常用到矩阵的Schur补性质。考虑一个矩阵,并将S分块:

其中S_{11}r \times r维的。假定S_{11}是非奇异的,则称为“S_{11}在S中的Sehar补”。以下引理给出了矩阵的Schur补性质。

 引理1 Schur补性质

对给定的对称矩阵:

其中S_{11}r \times r维的。以下3个条件是等价的:

(1)S<0

(2)S_{11}<0,S_{22}-S{12}^TS_{11}^{-1}S_{12}<0

(3)S_{22}<0,S_{11}-S{12}S_{22}^{-1}S_{12}^T<0

(证明方法可参考 俞立的《鲁棒控制——线性矩阵不等式处理方法》 p8页)

        在一些控制问题中,经常遇到二次型矩阵不等式[3]:

 其中A,B,Q=Q^T>0,R=R^T>0是给定的适当维数的常数矩阵,P是对称矩阵变量,则应用引理1,可以将上面的矩阵不等式的可行性问题转化成一个等价的矩阵不等式:

的可行性问题,而后者是一个关于矩阵变量P的线性矩阵不等式。

        因此在控制问题中,我们经常要设计Lyapunov函数V(t),并为了保证系统稳定性要求\dot{V}(t)<0,那我们就可以把\dot{V}(t)<0这个不等式转化为线性矩阵不等式的形式,并用matlab中的YALMIP直接求解出即可。

(四)LMI中常见引理

引理2(广义KYP引理[4])

给定矩阵\ThetaF以及\Phi\Psi,且N_w代表T_wF的零空间,其中T_w=\left [ \begin{matrix} I & -jwI \end{matrix} \right ],则不等式

N^{*}_w\Theta N_w<0,w\epsilon \left [ \begin{matrix} \varpi _1 & \varpi _2 \end{matrix} \right ]

成立,当且仅当存在对称矩阵P以及Q>0,使得成立,其中

其中*表示矩阵的共轭转置,j是虚数单位,\Phi \otimes P表示右克罗内克积,即

\Phi \otimes P= \left [ \begin{matrix} 0 & P \\ P & 0 \end{matrix} \right ]

推论1(广义KYP引理推论[4])

针对线性系统G(jw)为系统由扰动到控制输出的传递函数,则对于给定的对称矩阵\Pi,下面两个陈述是等价的:

1)有限频域不等式

成立。

2)存在对称矩阵Р和Q满足Q>0,使得

成立,其中

表示矩阵\Pi的右上块和右下块,矩阵中*表示其对应块的转置。

引理3(射影定理[1])

对于给定的标量\Gamma ,\Lambda ,\Theta,矩阵F满足\Gamma F\Lambda +(\Gamma F\Lambda )^T+\Theta <0,当且仅当下面两个条件成立:

引理3(反射影定理[1]):P是给定的正定对称矩阵,不等式\Psi +S+S^T<0等价于下面的线性矩阵不等式(LMI)求解问题:

式中符号[W]_s用以表示矩阵W与其转置的加和,即[W]_s=W+W^T

引理4 (Jensen不等式[5,6]

对于任意正定对称常数矩阵M\epsilon R^{n \times n},标量r满足r>0,并且存在一个向量w:[0,r]\rightarrow R^n,那么以下不等式成立:

引理5 (Finsler's引理[7]):

x\epsilon R^n,p\epsilon S^n,以及H\epsilon R^{m \times n}满足H的秩小于n,rank(H) = r <n.则下列两式等价:

参考文献

【1】Apkarian P,Tuan H D,Bernussou J.Continuous-Time analysis,eigenstructure as-signment,and H2synthesis with enhanced Linear Matrix Inequalities(LMI)char-acterizations[J].IEEE Transactions on Automatic Control,2001,42(12):1941–1946.

【2】《鲁棒控制理论及应用》王娟 张涛 徐国凯

【3】《鲁棒控制——线性矩阵不等式处理方法》俞立

【4】Iwasaki T,Hara S.Generalized KYP Lemma:unified frequency domain inequal-ities with design applications[J].IEEE Transactions on Automatic Control,2005,50(1):41–59.

【5】Wu J,Chen X,Gao H.H∞filtering with stochastic sampling[J].Signal Proces-siong,2010,90(4):1131–1145.

【6】Gao H,Wu J,Shi P.Robust sampled-data H∞control with stochastic sampling[J].Automatica,2009,45(7):1729–1736.

【7】Qiu J,Feng G,Yang J.New results on robust energy-to-peakfiltering for discrete-time switched polytopic linear systems with time-varying delay[J].IET ControlTheory and Applications,2008,2(9):795–806.

·

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1006640.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

如何进行有效测试?——基于MEAP诞生的软件测试实用指南告诉你

近年来出现了一 些新的出版方式&#xff0c;MEAP(Manning Early Access Program)就是其中的一种&#xff0c;把开源运动扩展到出版行业。在MEAP中&#xff0c;读者可在图书出版前逐章阅读早期版本。在作者写作过程中&#xff0c;读者可以及时提供反馈&#xff0c;帮助作者写出更…

苹果发布会,华为成热搜第一?

文 | 琥珀消研社 作者 | April叶 琥珀消研社快评&#xff1a;雷蒙多给华为手机“代言”后&#xff0c;苹果库克又给华为“代言”了&#xff0c;苹果发布会当晚&#xff0c;华为竟然登顶热搜第一#苹果发布会 #华为 #Mate 60 Pro#手机#科技 苹果开发布会&#xff0c;华为却冲上…

Matplotlib | 高阶绘图案例【3】- 五大战区高校排名

文章目录 &#x1f3f3;️‍&#x1f308; 1. 导入模块&#x1f3f3;️‍&#x1f308; 2. 数据处理2.1 高效数据2.2 学校排名 &#x1f3f3;️‍&#x1f308; 3. 绘图3.1 绘制图布&#xff0c;设置极坐标系3.2 绘制学校排名柱状图3.3 绘制五大战区扇形区域3.4 添加战区、学校…

javaSE___语法基础(关键字、标识符、数据类型、字符、字符串、运算符)

语法基础 一、关键字的概述和使用二、标识符的概述和组成规则三、标识符中常见的命名规则四、数据类型的概述和分类五、数据类型的概述和分类六、数据类型转换之强制转换七、数据类型转换之强制转换超出范围八、字符和字符串参与运算九、算术运算符的基本用法十、算术运算符和-…

闭包的理解

1.什么是闭包&#xff1f; 变量的私有化。一个函数内的变量,随着函数的执行完毕,对于的变量也会随着销毁,闭包可以让变量在函数执行完毕之后不必销毁,通常将这个变量通过匿名函数的形式return出去,这个变量只能被访问,不能被修改。 2.证明变量执行玩被销毁 (1)函数体没有被包…

((type *)0)->member的用法

问题缘由&#xff0c;在学习 rt-thread 内核的时候遇到了这么一行代码&#xff1a; to_thread rt_list_entry(rt_thread_priority_table[0].next,struct rt_thread, tlist); 而 rt_list_entry 的宏定义如下&#xff1a; /* 已知一个结构体里面的成员的地址&#xff0c;反推…

利用 Python 中的地理空间数据与 GeoPandas

推荐&#xff1a;使用 NSDT编辑器快速搭建3D应用场景 空间数据的真正潜力在于它能够连接数据点及其各自的位置&#xff0c;为高级分析创造无限的可能性。地理空间数据科学是数据科学中的一个新兴领域&#xff0c;旨在利用地理空间信息并通过空间算法和机器学习或深度学习等先进…

Android:使用命令行发现keytool不是内部命令解决办法

一、前言&#xff1a;最近在搞引入高德地图的SDK&#xff0c;发现需要给app签名打包。记录一下。 二、当我在命令行中输入keytool的时候说keytool不是内部命令 解决方案&#xff1a; 找到系统属性--------高级----------点击环境变量 双击点开 找到java\jre1.8.0_202\bin新建…

GLTF在线编辑器

GLTF在线编辑器提供了一个内置的模型查看器&#xff0c;可以加载和预览 glTF/glb 文件。用户可以在不用安装任何插件的情况下直接在浏览中快速查看和编辑器3D模型。 它的功能特点如下&#xff1a; 1、打开GLTF模型 用户可以在GLTF编辑器中拖入GLB/GLTF模型或者选择打开本地GL…

Bootstarp4 设计网页轮播组件

很多网站都有广告轮播功能&#xff0c;可使用bootstrap4中的carousel组件非常简单的实现。 目录 下载bootstrap4 轮播功能实现 简单实现轮播组件 增加标识图标 增加标题和说明 切换淡入淡出 设置数据间隔 总结 下载bootstrap4 下载 Bootstrap Bootstrap v4 中文文档 …

Android13 下拉菜单栏中添加快捷截图按钮

Android 13 原生系统下拉状态栏中是没有快捷截图按钮,现在需要添加快捷截图功能。 添加快捷截图功能后的效果图: 涉及修改的文件如下: modified: vendor/mediatek/proprietary/packages/apps/SystemUI/res/values/config.xml modified: vendor/mediatek/proprietary/…

软件测试的流程有哪些?

一、需求调研 1. 委托方提供资料 A. 填写测试委托申请表 B. 操作手册 C. 开发需求规格说明书 D. 开发合同及招标文件等 2. 双方技术沟通确定测试具体内容&#xff0c;如功能性测试、性能效率测试、信息安全性测试、兼容性测试、可靠性测试等。 3. 我方给出测试方案及报价…

SpringMVC学习|Servlet回顾、理解SpringMVC小demo、SpringMVC原理

Servlet回顾 创一个空的maven父工程&#xff0c;导入相关依赖&#xff0c;测试的、spring的、servlet的、jsp以及jstl标签的。 创建一个空的maven子项目&#xff0c;并添加web支持 为了保险起见&#xff0c;在子项目中添加servlet和jsp的依赖 编写一个Servlet类&#xff0c;实…

【Python】从入门到上头—常用内置模块基础应用(13)

文章目录 datetimebase64hashlibhmacurllibXMLHTMLParserrandom小结 datetime datetime是Python处理日期和时间的标准库。 获取当前日期时间 from datetime import datetime now datetime.now() # 获取当前datetime print(now) #2023-09-13 10:28:48.621343 print(type(no…

JSP ssm 网上求职管理系统myeclipse开发mysql数据库springMVC模式java编程计算机网页设计

一、源码特点 JSP ssm 网上求职管理系统是一套完善的web设计系统&#xff08;系统采用SSM框架进行设计开发&#xff0c;springspringMVCmybatis&#xff09;&#xff0c;对理解JSP java编程开发语言有帮助&#xff0c;系统具有完整的源代码和数据库&#xff0c;系统主要采…

3900页手册415集视频426G资料迅为RK3568开发板

资料优势 专为3568编写|迅为原创|拒绝网络拼凑 20个手册2800页手册进行结构分层适用于学习与开发 为了方便大家清晰快速的学习&#xff0c;迅为iTOP-3568开发板手册资料全面升级&#xff0c;对手册内容进行了结构分层&#xff0c;共计20个文档&#xff0c;超2800页的资料专为…

sql注入之盲注总结

死在山野的风里&#xff0c;活在自由的梦里 sql注入之盲注总结 布尔 - 逻辑判断时间 - 延时判断报错 - 报错回显&#xff08;强制性报错 &#xff09;以DWVA靶场为例 盲注就是在注入的过程中&#xff0c;获取的数据不能显示到前端页面&#xff0c;此时&#xff0c;我们需要利用…

2023年武汉初级、中级、副高工程师职称评审条件是什么?甘建二告诉你

武汉中级职称评审通知已经出来了&#xff0c;武汉市中级、副高申报时间是8月31--9月18号&#xff0c;那么评这些职称需要满足什么条件呢&#xff1f;今天甘建二给大家分享一下评职称学历年限、专业等是怎么计算的&#xff0c;点赞收藏慢慢看。 一、武汉初级职称评审条件&#x…

【MySQL进阶】SQL性能分析

一、SQL性能分析 1.SQL执行频率 MySQL 客户端连接成功后&#xff0c;通过 show [session|global] status 命令可以提供服务器状态信 息。通过如下指令&#xff0c;可以查看当前数据库的 INSERT 、 UPDATE 、 DELETE 、 SELECT 的访问频次&#xff1a; -- session 是查看当…

jvm深入研究文档--整体概念

阿丹&#xff1a; 精通JVM对于一个java工程师非常重要&#xff0c;要是深入了解了jvm就可以有效的面对下面的问题 程序调优&#xff1a;JVM的配置和调优对于程序的运行有着至关重要的影响。不同的业务场景需要不同的JVM配置&#xff0c;比如设置不同的垃圾收集器、调整新生代和…