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

目录

(一)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.mzph.cn/news/78073.shtml

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

相关文章

Redis 缓存过期淘汰策略

面试题&#xff1a; 生产上你们的redis内存设置多少&#xff1f;如何配置、修改redis的内存大小如果内存满了你怎么办redis清理内存的方式&#xff1f;定期出测除和惰性曲除了解过吗redis缓存淘太策略有哪些&#xff1f;分别是什么&#xff1f;你用那个&#xff1f;redis的LRU…

做机器视觉工程师,其实挺没意思的

3.康耐视VisionPro高级脚本系列教程-3.脚本编辑错误和运行错误调试方法&#xff0c;break和Contitinuee的差别_哔哩哔哩_bilibili 其实人生就是“有时有意思&#xff0c;有时没意思”。 心里有太多的不甘心&#xff0c;太多的苦水&#xff0c;是没法再吃学习的苦&#xff0c…

分库分表---理论

目录 一、垂直切分 1、垂直分库 2、垂直分表 3、垂直切分优缺点 二、水平切分 1、水平分库 2、水平分表 3、水平切分优缺点 三、数据分片规则 1、Hash取模分表 2、数值Range分表 3、一致性Hash算法 四、分库分表带来的问题 1、分布式事务问题 2、跨节点关联查询…

指针,动态内存分配

什么是指针 指针就是地址&#xff0c;地址就是内存单元的编号 指针变量&#xff0c;存放其他变量地址的变量 &#xff0c;存放地址的变量 第2句话&#xff0c;修改i的值也不影响p的值&#xff0c;可以看下面的例子&#xff0c;p放的是地址&#xff0c;p的地址修改了&#xff…

阿里云通义千问向全社会开放,近期将开源更大参数规模大模型

9月13日&#xff0c;阿里云宣布通义千问大模型已首批通过备案&#xff0c;并正式向公众开放&#xff0c;广大用户可登录通义千问官网体验&#xff0c;企业用户可以通过阿里云调用通义千问API。 通义千问在技术创新和行业应用上均位居大模型行业前列。IDC最新的AI大模型评估报告…

TypeScript高级类型

&#x1f3ac; 岸边的风&#xff1a;个人主页 &#x1f525; 个人专栏 :《 VUE 》 《 javaScript 》 ⛺️ 生活的理想&#xff0c;就是为了理想的生活 ! 目录 映射类型&#xff08;Mapped Types&#xff09; 1. Readonly 2. Partial 3. Pick 4. Record 条件类型&#xf…

Mysql中group by 使用中发现的问题

当使用 MySQL 的 GROUP BY 语句时&#xff0c;根据指定的列对结果进行分组。在 GROUP BY 分组时&#xff0c;如果某个字段在分组中有多个不同的值&#xff0c;那么就会出现你提到的该字段一直在变化的情况。 这种情况通常是由于在 GROUP BY 中选择的字段与其他非聚合字段不兼容…

互联网医院App开发:构建医疗服务的技术指南

互联网医院App的开发是一个复杂而具有挑战性的任务&#xff0c;但它也是一个充满潜力的领域&#xff0c;可以为患者和医疗专业人员提供更便捷的医疗服务。本文将引导您通过一些常见的技术步骤来构建一个简单的互联网医院App原型&#xff0c;以了解该过程的基本概念。 技术栈选…

打印日志遇到的问题,logback与zookeeper冲突

在做项目时需要打印日志引入了logback打印日志&#xff0c;但是一直无法打印&#xff0c;于是一路查找原因。发现zookeeper中默认带的有个logback和我自己引入的logback版本冲突了&#xff0c;这样直接使用exclusions标签将zookeeper中自带的日志框架全部排除即可 按理说到这一…

Windows下防火墙端口配置

在电脑或者服务器上部署某个应用后&#xff0c;如果需要对外提供服务可能就需要在主机防火墙上设置开启需要的端口&#xff0c;那么具体怎样操作呢 1.打开windows防火墙 2.设置防火墙入站规则 如下图“高级安全Windows Defender 防火墙”页面&#xff0c;点击左侧“入站规则”…

离散数学 学习 之 一阶逻辑基本概念 (一 )

个体词可以 理解为主语 &#xff0c; 3 不是偶数 &#xff0c;3 就是 个体常项 凡整数都能被 2 整除 &#xff0c; 整数就是 个体变项 上面的谓词是 &#xff08; 1 &#xff09; 是无理数 &#xff08; 2 &#xff09; 是有理数 &#xff08; 3 &#xff09; 与 同岁 &#xf…

机器人任务挖掘与智能超级自动化技术解析

本文为上海财经大学教授、安徽财经大学学术副校长何贤杰出席“会计科技Acctech应对不确定性挑战”高峰论坛时的演讲内容整理。何贤杰详细介绍了机器人任务挖掘与智能超级自动化技术的发展背景、关键技术和应用场景。 从本质来说&#xff0c;会计是非常适合智能化、自动化的。会…

JWT安全

JWT是什么&#xff1f; JSON Web Token (JWT)是一个开放标准(RFC 7519)&#xff0c;它定义了一种紧凑的、自包含的方式&#xff0c;用于作为JSON对象在各方之间安全地传输信息。 JWT全称JSON Web Token&#xff0c;是一种传输信息的标准&#xff0c;用于JSON对象在各方之间安全…

SQL7 查找年龄大于24岁的用户信息

描述 题目&#xff1a;现在运营想要针对24岁以上的用户开展分析&#xff0c;请你取出满足条件的设备ID、性别、年龄、学校。 用户信息表&#xff1a;user_profile iddevice_idgenderageuniversityprovince12138male21北京大学Beijing23214male复旦大学Shanghai36543female20…

Android性能优化之应用瘦身(APK瘦身)

关于作者&#xff1a;CSDN内容合伙人、技术专家&#xff0c; 从零开始做日活千万级APP。 专注于分享各领域原创系列文章 &#xff0c;擅长java后端、移动开发、人工智能等&#xff0c;希望大家多多支持。 目录 一、导读二、概览2.1 apk组成 三、优化方向3.1 源代码3.1.1 代码混…

GDB的C++调试方法

本文记录基础的GDB调试过程&#xff0c;包含指令如下&#xff1a; 文章目录 准备编译文件GDB启动GDB开启代码行设置断点运行程序查看pc的指令查看监视的变量以及断点设置快照checkpoint实验1实验2 nextnextistepx/i $pcfinishinfo break 和 delete Numrefbreak col. if condit…

【操作系统】进程控制

进程控制&#xff1a;创建新进程&#xff0c;撤销已有进程&#xff0c;实现进程状态转换等。 原语&#xff1a;进程控制用的程序段。执行期间不允许中断&#xff0c;用&#xff02;关中断&#xff02;和&#xff02;开中断&#xff02;指令&#xff08;特权指令&#xff09;实…

创建java文件 自动添加作者、时间等信息 – IDEA 技巧

2023 09 亲测 文章目录 效果修改位置配置信息 效果 每次创建文件的时候&#xff0c;自动加上作者、时间等信息 修改位置 打开&#xff1a;File —> Settings —> Editor —> File and Code Templates —> includes —> FileHeader 配置信息 /*** author : Java…

zookeeper/HA集群配置

1.zookeep配置 1.1 安装4台虚拟机 &#xff08;1&#xff09;按照如下设置准备四台虚拟机&#xff0c;其中三台作为zookeeper&#xff0c;配置每台机器相应的IP&#xff0c;hostname&#xff0c;下载vim&#xff0c;ntpdate配置定时器定时更新时间&#xff0c;psmisc&#xff…

C++与C编译后符号表对比(一百九十二)

简介&#xff1a; CSDN博客专家&#xff0c;专注Android/Linux系统&#xff0c;分享多mic语音方案、音视频、编解码等技术&#xff0c;与大家一起成长&#xff01; 优质专栏&#xff1a;Audio工程师进阶系列【原创干货持续更新中……】&#x1f680; 人生格言&#xff1a; 人生…