当前位置: 首页 > news >正文

线性矩阵不等式(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.

·

相关文章:

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

目录 &#xff08;一&#xff09;Matlab中的LMI处理工具包 &#xff08;二&#xff09;为什么LMI成为控制理论领域重要工具&#xff1f; &#xff08;三&#xff09;LMI在与Lyapunov不等式的关系 &#xff08;1&#xff09;线性矩阵不等式 &#xff08;2&#xff09;线性矩阵…...

如何在Python爬虫程序中使用HTTP代理?

在进行网络爬虫时&#xff0c;我们经常需要使用代理服务器来隐藏自己的真实IP地址&#xff0c;以避免被目标网站封禁或限制访问。本文将介绍如何将HTTP代理配置到Python爬虫程序中使用。 什么是HTTP代理&#xff1f; HTTP代理是一种网络代理&#xff0c;它充当客户端和服务器之…...

ARM架构指令集--专用指令

四、状态寄存器专用指令 CPSR寄存器-N Z C V T为0时 为ARM状态 F为0时 为开启FIQ状态 I为0时 为开启IRQ状态 图1 图2 一开始都是SVC指令&#xff0c;因为在操作系统启动的时候&#xff0c;在做一些初始化的操作&#xff0c;不允许被打断 图3 复位后CPSR寄存器为0xD3--…...

免费IP类api接口:含ip查询、ip应用场景查询、ip代理识别、IP行业查询...

免费IP类api接口&#xff1a;含ip查询、ip应用场景查询、ip代理识别… IP归属地-IPv6区县级&#xff1a;根据IP地址&#xff08;IPv6版本&#xff09;查询归属地信息&#xff0c;包含国家、省、市、区县和运营商等信息。IP归属地-IPv6城市级&#xff1a;根据IP地址&#xff08…...

Android Studio 中AGP ,Gradle ,JDK,SDK都是什么?

当进行 Android 开发时&#xff0c;以下是关键概念和工具的解释&#xff1a; Android Gradle Plugin&#xff1a; Android Gradle Plugin 是一个由 Google 提供的构建工具&#xff0c;它与 Gradle 配合使用来构建和打包 Android 应用。它提供了一组任务和功能&#xff0c;使开发…...

算法通关18关 | 回溯模板如何解决复原IP问题

18关的前几篇文章看过之后&#xff0c;对回溯的模板问题基本解题思路就知道了&#xff0c;就是固定的for循环问题&#xff0c;外层for循环控制横向&#xff0c;递归控制纵向&#xff0c;还要考虑撤销操作和元素是否能被重复利用问题&#xff0c;重复利用的情景较少&#xff0c;…...

Layui快速入门之第五节 导航

目录 一&#xff1a;基本概念 导航依赖element模块 API 渲染 属性 事件 二&#xff1a;水平导航 常规用法&#xff1a; 三&#xff1a;垂直导航 四&#xff1a;侧边垂直导航 五&#xff1a;导航主题 六&#xff1a;加入徽章等元素 七&#xff1a;面包屑导航 ps&a…...

使用分支——Git Checkout

这篇文章写的挺好&#xff1b; https://zhuanlan.zhihu.com/p/465954849 这里要注意&#xff0c;git 新的命令&#xff0c;通过 git switch 切换分支&#xff0c;虽然git checkout 分支 还可以用&#xff1b; 游离状态的HEADS 在我们已经见识到git checkout命令对于分支的三…...

【2023】数据挖掘课程设计:基于TF-IDF的文本分类

目录 一、课程设计题目 基于TF-IDF的文本分类 二、课程设计设置 1. 操作系统 2. IDE 3. python 4. 相关的库 三、课程设计目标 1. 掌握数据预处理的方法&#xff0c;对训练集数据进行预处理&#xff1b; 2. 掌握文本分类建模的方法&#xff0c;对语料库的文档进行建模…...

java.lang.NoSuchMethodError: java.lang.reflect.Field.trySetAccessible()Z

java.lang.NoSuchMethodError: java.lang.reflect.Field.trySetAccessible()Z 将JDK升级为11即可。 File --Project Structure – SDK Location --Gradle Setting --Gradle JDK 选择11...

如何使用SQL系列 之 如何在MySQL中使用存储过程

引言 通常&#xff0c;当使用关系型数据库时&#xff0c;你直接在应用程序代码中发出单独的结构化查询语言(SQL)查询来检索或操作数据&#xff0c;如SELECT、INSERT、UPDATE或DELETE。这些语句直接作用于并操作底层数据库表。如果相同的语句或一组语句中使用多个应用程序访问同…...

用 Github Codespaces 免费搭建本地开发测试环境

如何丝滑地白嫖一个本地开发环境&#xff1f;怎么新建一个代码空间&#xff1f; 1&#xff1a;通过Github网页新建2&#xff1a;通过VSCode插件新建 为代码创建相应的开发测试环境 如何丝滑地白嫖一个本地开发环境&#xff1f; 使用Codespaces为开发者解决这样的痛点&#xf…...

PyTorch实战-实现神经网络图像分类基础Tensor最全操作详解(一)

目录 前言 一、PyTorch数据结构-Tensor 1.什么是Tensor 2.数据Tensor使用场景 3.张量形态 标量&#xff08;0D 张量&#xff09; 向量&#xff08;1D 张量&#xff09; 矩阵(2D张量) 3D 张量与高维张量 二、Tensor的创建 1. 从列表或NumPy数组创建 2. 使用特定的初始…...

第29章_瑞萨MCU零基础入门系列教程之改进型环形缓冲区

本教程基于韦东山百问网出的 DShanMCU-RA6M5开发板 进行编写&#xff0c;需要的同学可以在这里获取&#xff1a; https://item.taobao.com/item.htm?id728461040949 配套资料获取&#xff1a;https://renesas-docs.100ask.net 瑞萨MCU零基础入门系列教程汇总&#xff1a; ht…...

如何搭建一个react项目(详细介绍)

要搭建一个基本的 React 项目&#xff0c;你需要执行以下步骤。在开始之前&#xff0c;请确保你已经安装了 Node.js 和 npm&#xff08;Node 包管理器&#xff09;。 搭建一个React项目 1&#xff0c;创建项目目录2&#xff0c;初始化项目3&#xff0c;安装 React 和 ReactDOM4…...

ActiveMQ用法

ActiveMQ 和 JMS的关系&#xff1f; ActiveMQ是流行的开源消息中间件&#xff0c;JMS是Java平台定义的一种消息传递的标准。ActiveMQ实现了JMS规范&#xff0c;因此可以使用JMS API来与ActiveMQ进行交互。 JMS定义了一种标准的API。API包括了一些接口和类&#xff0c;用于创建…...

TouchGFX之缓存位图

位图缓存是专用RAM缓冲区&#xff0c;应用可将位图保存&#xff08;或缓存&#xff09;在其中。 如果缓存了位图&#xff0c;在绘制位图时&#xff0c;TouchGFX将自动使用RAM缓存作为像素来源。位图缓存在许多情况下十分有用。 从RAM读取数据通常比从闪存读取要快&#xff08;特…...

线性代数的本质(十)——矩阵分解

文章目录 矩阵分解LU分解QR分解特征值分解奇异值分解奇异值分解矩阵的基本子空间奇异值分解的性质矩阵的外积展开式 矩阵分解 矩阵的因式分解是把矩阵表示为多个矩阵的乘积&#xff0c;这种结构更便于理解和计算。 LU分解 设 A A A 是 m n m\times n mn 矩阵&#xff0c;…...

vue实现鼠标拖拽div左右移动的功能

直接代码&#xff1a; <template><div class"demo"><div class"third-part" id"发展历程"><div class"title">发展历程</div><div class"content" id"nav" v-if"dataList…...

基于Python和mysql开发的商城购物管理系统分为前后端(源码+数据库+程序配置说明书+程序使用说明书)

一、项目简介 本项目是一套基于Python和mysql开发的商城购物管理系统&#xff0c;主要针对计算机相关专业的正在做毕设的学生与需要项目实战练习的Python学习者。 包含&#xff1a;项目源码、项目文档、数据库脚本等&#xff0c;该项目附带全部源码可作为毕设使用。 项目都经过…...

MySQL内外连接、索引特性

目录 内连接 外连接 索引特性 理解索引 删除索引 MySQL内外连接是一种用于联接两个或多个表的操作。内连接只返回满足连接条件的行&#xff0c;外连接返回满足条件和不满足条件的行。 内连接 SQL如下&#xff1a; SELECT ... FROM t1 INNER JOIN t2 ON 连接条件 [INNER …...

滚动条设置

不同浏览器滚动条样式及滚动定位 是否可以滚动 overflow&#xff1a;scroll overflow&#xff1a;autooverflow:scroll – 只有超出了盒子才会有滚动条 overflow:auto – 一直有滚动的盒子&#xff0c;只是超出了盒子才会出现滚动条滑块&#xff0c;可以滚动 谷歌浏览器滚动…...

【AI】机器学习——感知机

文章目录 4.1 感知机基本概念4.2 策略4.2.1 数据集的线性可分性4.2.2 学习策略目标损失函数的构造关于距离的解释 4.3 算法4.3.1 原始形式损失函数的梯度下降法 4.3.2 PLA例题4.3.3 算法收敛性 4.4 PLA对偶形式4.4.1 原始PLA分析4.4.2 PLA对偶形式4.4.3 优点 4.1 感知机基本概念…...

蓝牙遥控器在T2-U上的应用

文章目录 简介优势使用流程示例代码遥控器命令表遥控器代码实现开启遥控器配对功能运行 简介 Tuya beacon 协议是基于 BLE 广播通信技术&#xff0c;完善配对解绑、组包拆包、群组管理、加密解密、安全策略&#xff0c;形成的一种轻量、安全的可接入涂鸦云的蓝牙协议。 蓝牙 …...

数据驱动的数字营销与消费者运营

引言&#xff1a;基于海洋馆文旅企业在推广宣传中&#xff0c;如何通过指标体系量化分析广告收益对业务带来的收益价值的思考&#xff1f; 第一部分:前链路引流投放的策略与实战 1.1 动态广告的实现: 偶然与必然 动态广告是一种基于实时数据和用户行为的广告形式&#xff0c;它…...

Qt点亮I.MX6U开发板的一个LED

本篇开始将会介绍与开发版相关的Qt项目&#xff0c;首先从点亮一个LED开始。I.MX6U和STM32MP157的相关信息都会用到&#xff0c;但是后期还是将I.MX6U的学习作为重点。当然其他开发版的开发也可以参考本博文。 文章目录 1. Qt是如何操控开发板上的一个LED2. 出厂内核设备树中注…...

网络摄像头-流媒体服务器-视频流客户端

取电脑的视频流 当涉及交通事件检测算法和摄像头视频数据处理时&#xff0c;涉及的代码案例可能会非常复杂&#xff0c;因为这涉及到多个组件和技术。以下是一个简单的Python代码示例&#xff0c;演示如何使用OpenCV库捕获摄像头视频流并进行实时车辆检测&#xff0c;这是一个…...

Django05_反向解析

Django05_反向解析 5.1 反向解析概述 随着功能的不断扩展&#xff0c;路由层的 url 发生变化&#xff0c;就需要去更改对应的视图层和模板层的 url&#xff0c;非常麻烦&#xff0c;不便维护。这个时候我们可以通过反向解析&#xff0c;将 url解析成对应的 试图函数 通过 path…...

基于HTML、CSS和JavaScript制作一个中秋节倒计时网页

&#x1f482; 个人网站:【工具大全】【游戏大全】【神级源码资源网】&#x1f91f; 前端学习课程&#xff1a;&#x1f449;【28个案例趣学前端】【400个JS面试题】&#x1f485; 寻找学习交流、摸鱼划水的小伙伴&#xff0c;请点击【摸鱼学习交流群】 这个项目使用HTML、CSS和…...

富斯I6刷10通道固件

使用USB转串口模块刷写10通道固件 一、下载固件 1. 十通道英文固件 下载地址: https://github.com/benb0jangles/FlySky-i6-Mod-/tree/master 选择 FlySky-i6-Mod–master \ 10ch Mod i6 Updater \ 10ch_MOD_i6_Programmer_V1 路径下的文件,亲测可用。 2. 原版六通道中…...

wordpress 表单/网络营销工资一般多少

那么&#xff0c;ELK 到底是什么呢&#xff1f; “ELK”是三个开源项目的首字母缩写&#xff0c;这三个项目分别是&#xff1a;Elasticsearch、Logstash 和 Kibana。Elasticsearch 是一个搜索和分析引擎。Logstash 是服务器端数据处理管道&#xff0c;能够同时从多个来源采集数…...

企业网站建设预算方案/网站访问量

iOS8之后&#xff0c;苹果推出了WebKit这个框架&#xff0c;用来替换原有的UIWebView&#xff0c;新的控件优点多多&#xff0c;不一一叙述。由于一直在适配iOS7&#xff0c;就没有去替换&#xff0c;现在仍掉了iOS7&#xff0c;以为很简单的就替换过来了&#xff0c;然而在替换…...

专门做正品的网站/搜索引擎优化的方法包括

前言 Redis 是我们目前大规模使用的缓存中间件&#xff0c;由于他强调高效而又便捷的功能&#xff0c;得到了广泛的使用。单节点的Redis已经达到了很高的性能&#xff0c;为了提高可用性我们可以使用Redis 集群。本文参考了Rdis的官方文档和使用Redis官方提供的Redis Cluster工…...

网站设计公司哪里好/超级推荐的关键词怎么优化

参考:https://docs.opencv.org/master/da/d13/tutorial_aruco_calibration.html 说明:从视频中挑选标定样本,键盘输入“c”:选择该帧作为样本;键盘输入“d”:不选择该帧作为样本,进行下一帧挑选;键盘输入“esc”:样本挑选完毕; #include <opencv2/highgui.hpp&…...

wordpress 日期/杭州seo代理公司

Sed简介sed 是一种在线编辑器&#xff0c;它一次处理一行内容。处理时&#xff0c;把当前处理的行存储在临时缓冲区中&#xff0c;称为“模式空间”&#xff08;pattern space&#xff09;&#xff0c;接着用sed命令处理缓冲区中的内容&#xff0c;处理完成后&#xff0c;把缓冲…...

初中信息科技怎么自己做网站/网站建设报价单模板

-- 学习C可分为4个层次&#xff1a; 第一层次&#xff0c;C基础&#xff1a;挑选一本入门书籍&#xff0c;如《C Primer》、《C大学教程》、或Stroustrup撰写的经典《C程序设计语言》或他一年半前的新作《C程序设计原理与实践》&#xff0c;而一般C课程也止于此&#xff0c;另…...