线性矩阵不等式LMI与李雅普诺夫Lyapunov稳定性

文章目录

  • 线性矩阵不等式(Linear Matrix Inequality,LMI)
      • 例子
  • Lyapunov稳定性
  • Schur Complement
      • 定义
      • Schur Complement作用/性质
      • 利用Schur Complement将LMI和Lyapunov联系起来

线性矩阵不等式(Linear Matrix Inequality,LMI)

形式为
LMI ( y ) = A 0 + A 1 y 1 + A 2 y 2 + ⋯ ≥ 0 \text{LMI}(y)=A_0+A_1y_1+A_2y_2+\cdots \geq 0 LMI(y)=A0+A1y1+A2y2+0
其中 A 0 , A 1 , A 2 , . . . A_0,A_1,A_2,... A0,A1,A2,...为对称方阵。

例子


LMI ( y ) = [ y 1 + y 2 y 2 + 1 y 1 + 1 y 3 ] , \text{LMI}(y)=\left[ \begin{matrix} y_1+y_2 & y_2+1 \\ y_1+1&y_3\end{matrix} \right], LMI(y)=[y1+y2y1+1y2+1y3],
则对应
A 0 = [ 0 1 1 0 ] , A 1 = [ 1 0 1 0 ] , A 2 = [ 1 1 0 0 ] , A 3 = [ 0 0 0 1 . ] A_0=\left[ \begin{matrix} 0 & 1 \\1&0\end{matrix} \right], A_1=\left[ \begin{matrix} 1& 0\\ 1&0\end{matrix} \right],A_2=\left[ \begin{matrix} 1 &1 \\ 0&0\end{matrix} \right],A_3=\left[ \begin{matrix} 0 & 0 \\ 0&1\end{matrix} .\right] A0=[0110],A1=[1100],A2=[1010],A3=[0001.]
随着解决线性矩阵不等式的内点法的提出、以及 MATLAB 软件中 LMI 工具箱的推出,线性矩阵不等式这一工具越来越受到人们的注意和重视。

Lyapunov稳定性

假设可以找到一个正定的Lyapunov函数 V V V(即 V > 0 V>0 V>0)且 V ˙ < 0 \dot{V}<0 V˙<0,则可以证明系统是稳定的。以线性系统为例:
x ˙ = A x + B u . \dot{x}=Ax+Bu. x˙=Ax+Bu.
假设反馈控制
u = − K x . u=-Kx. u=Kx.
取Lyapunov函数为
V ( x ) = x T P x , V(x)=x^{T}Px, V(x)=xTPx,
其中 P P P正定且对称,即 P ≻ 0 , P = P T P\succ0,P=P^{T} P0,P=PT。Lyapunov的导数为
V ˙ ( x ) = x T P x ˙ + x ˙ T P x = x T P ( A − B K ) x + x T ( A − B K ) T P x = − x T Q x , \begin{aligned} \dot{V}(x)= & x^TP\dot{x}+\dot{x}^TPx \\ =&x^TP(A-BK)x+x^T(A-BK)^TPx\\ =&-x^TQx, \end{aligned} V˙(x)===xTPx˙+x˙TPxxTP(ABK)x+xT(ABK)TPxxTQx,
其中
Q = − ( A T P + P A − P B K − K T B T P ) . Q=-(A^TP+PA-PBK-K^TB^TP). Q=(ATP+PAPBKKTBTP).
若能证明 Q ≻ 0 Q \succ 0 Q0,则该系统渐近稳定。

最优控制中常取
K = − 1 2 R − 1 B T P T , K=-\frac{1}{2}R^{-1}B^TP^T, K=21R1BTPT,
其中,前提矩阵 R R R满足 R = R T ≻ 0 R=R^T \succ 0 R=RT0 R − 1 R^{-1} R1存在且有界,于是,
Q = − ( A T P + P A − P B R − 1 B T P T ) . (1) Q=-(A^TP+PA-PBR^{-1}B^TP^T). \tag{1} Q=(ATP+PAPBR1BTPT).(1)

Schur Complement

Schur Complement可用于对一个块矩阵进行等价转换。

定义

假设一个 n × n n \times n n×n的矩阵 M M M可以写成一个块矩阵形式:
M = [ A B C D ] . M=\left[ \begin{matrix} A & B \\ C & D \end{matrix} \right]. M=[ACBD].

  1. D D D是可逆的,则 D D D M M M中的舒尔补存在且为
    A − B D − 1 C ; A-BD^{-1}C; ABD1C;

  2. A A A是可逆的,则 A A A M M M中的舒尔补存在且为
    D − C A − 1 B . D-CA^{-1}B. DCA1B.
    “来历”:对方程
    [ A B C D ] [ x y ] = [ p q ] , \left[ \begin{matrix} A & B \\ C & D \end{matrix} \right] \left[ \begin{matrix} x \\ y \end{matrix} \right]=\left[ \begin{matrix} p \\ q \end{matrix} \right], [ACBD][xy]=[pq],
    使用高斯消元法,由 D D D可逆有
    ( A − B D − 1 C ) x = p − B D − 1 q . (A-BD^{-1}C)x=p-BD^{-1}q. (ABD1C)x=pBD1q.
    A A A可逆有
    ( D − C A − 1 B ) y = q − C A − 1 p . (D-CA^{-1}B)y=q-CA^{-1}p. (DCA1B)y=qCA1p.
    未知数前面的系数即为舒尔补。

Schur Complement作用/性质

  1. M M M分别变为上三角或者下三角矩阵:若 D D D可逆,则
    M = [ A B C D ] = [ I B D − 1 0 I ] [ A − B D − 1 C 0 0 D ] [ I 0 D − 1 C I ] ; M=\left[ \begin{matrix} A & B \\ C & D \end{matrix} \right]=\left[ \begin{matrix} I & BD^{-1} \\ 0 & I \end{matrix} \right]\left[ \begin{matrix} A-BD^{-1}C & 0 \\ 0 & D \end{matrix} \right]\left[ \begin{matrix} I & 0 \\ D^{-1}C & I \end{matrix} \right]; M=[ACBD]=[I0BD1I][ABD1C00D][ID1C0I];
    A A A可逆,则
    M = [ A B C D ] = [ I 0 C A − 1 I ] [ A 0 0 D − C A − 1 B ] [ I A − 1 B 0 I ] . M=\left[ \begin{matrix} A & B \\ C & D \end{matrix} \right]=\left[ \begin{matrix} I & 0\\ CA^{-1} & I \end{matrix} \right]\left[ \begin{matrix} A & 0 \\ 0 & D-CA^{-1}B \end{matrix} \right]\left[ \begin{matrix} I & A^{-1}B \\ 0 & I \end{matrix} \right]. M=[ACBD]=[ICA10I][A00DCA1B][I0A1BI].
    利用该性质可以快速求解矩阵 M M M的逆。

  2. 特殊性质:若 M M M是对称的,即
    M = [ A B B T C ] , M=\left[ \begin{matrix} A & B \\ B^T & C \end{matrix} \right], M=[ABTBC],
    C C C可逆,则有下列性质:

  3. M ≻ 0 M \succ 0 M0,则有且仅有 C ≻ 0 C \succ 0 C0 A − B C − 1 B T ≻ 0 A-BC^{-1}B^T \succ 0 ABC1BT0

  4. C ≻ 0 C \succ 0 C0,则 M ≻ 0 M \succ 0 M0有且仅有 A − B C − 1 B T ≻ 0 A-BC^{-1}B ^T\succ 0 ABC1BT0

利用Schur Complement将LMI和Lyapunov联系起来

利用舒尔补的特殊性质,式 ( 1 ) (1) (1)大于0等效为
[ − A T P − P A P B B T P T R ] ≻ 0. \left[ \begin{matrix} -A^TP-PA & PB \\ B^TP^T&R \end{matrix} \right] \succ 0. [ATPPABTPTPBR]0.
Lyapunov稳定性的判定条件转化为线性形式,从而方便用软件包数值求解。

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

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

相关文章

BUG:docker启动之后直接退出问题

示例如下&#xff1a; 问题排查&#xff1a; 启动命令 sudo docker run --privilegedtrue --runtimenvidia --shm-size80g -v /mmm_data_center:/mmm_data_center -v /imagecenter_new/:/imagecenter_new -v /data1:/data1 -v /mnt/offline_data/:/mnt/offline_data/ --neth…

cesium-测量高度垂直距离

cesium做垂直测量 完整代码 <template><div id"cesiumContainer" style"height: 100vh;"></div><div id"toolbar" style"position: fixed;top:20px;left:220px;"><el-breadcrumb><el-breadcrumb-i…

【51单片机】直流电机实验和步进电机实验

目录 直流电机实验直流电机介绍ULN2003 芯片介绍硬件设计软件设计实验现象 步进电机实验步进电机简介步进电机的工作原理步进电机极性区分双极性步进电机驱动原理单极性步进电机驱动原理细分驱动原理 28BYJ-48 步进电机简介软件设计 橙色 直流电机实验 在未学习 PWM 之前&…

理解Jetpack Compose中的`remember`和`mutableStateOf`

理解Jetpack Compose中的remember和mutableStateOf 在现代Android开发中&#xff0c;Jetpack Compose已经成为构建原生UI的首选工具。它引入了一种声明式的编程模式&#xff0c;极大地简化了UI开发。在Compose的世界里&#xff0c;remember和mutableStateOf是两个非常关键的函…

智慧城市:打造低碳未来,引领城市数字化转型新篇章

在“万物皆可数字化”的新时代浪潮下&#xff0c;智慧城市作为未来城市发展的先锋方向&#xff0c;正在以前所未有的速度和规模重塑我们的城市面貌。 智慧城市不仅是一个技术革新的标志&#xff0c;更是城市治理、民生服务等领域全面升级的重要引擎。 一、智慧城市的多元应用领…

玩家笔记:幻兽帕鲁搭建服务器开服教程

玩转幻兽帕鲁服务器&#xff0c;阿里云推出新手0基础一键部署幻兽帕鲁服务器教程&#xff0c;傻瓜式一键部署&#xff0c;3分钟即可成功创建一台Palworld专属服务器&#xff0c;成本仅需26元&#xff0c;阿里云服务器网aliyunfuwuqi.com分享2024年新版基于阿里云搭建幻兽帕鲁服…

高频一体式读写器现场应用

一体式读写器将所有组件集成在一个设备内&#xff0c;具有设计紧凑、安装简单、一体读写的特点&#xff0c;可在生产、仓库等领域中应用。高频一体式读写器现场应用 一体式读写器可以在在工业生产线上应用&#xff0c;读取产线设备的各种信息并记录&#xff0c;实现自动化控制和…

从零开始 TensorRT(4)命令行工具篇:trtexec 基本功能

前言 学习资料&#xff1a; TensorRT 源码示例 B站视频&#xff1a;TensorRT 教程 | 基于 8.6.1 版本 视频配套代码 cookbook 参考源码&#xff1a;cookbook → 07-Tool → trtexec 官方文档&#xff1a;trtexec 在 TensorRT 的安装目录 xxx/TensorRT-8.6.1.6/bin 下有命令行…

kubectl命令

kubenetes部署服务的流程 以部署一个nginx服务来说明kubernetes系统各个组件调用关系&#xff1a; 1. 首先要明确&#xff0c;一旦kubernetes环境启动之后&#xff0c;master和node都会将自身的信息存储到etcd数据库中 2. 一个nginx服务的安装请求会首先被发送到master节点的ap…

C++ dfs 与图有关的知识(四十七)【第七篇】

今天我们接着来学习树上搜索&#xff08;dfs深度优先搜索&#xff09; 1.树的深度与子树大小 树的深度&#xff1a;规定根结点是树的第一层&#xff0c;树根的孩子结点是树的第二层&#xff0c;以此类推&#xff0c;树的深度就是结点的最大层数。 根据定义&#xff0c;如果我们…

c语言实现greedy snake(贪吃蛇)

##第一个小项目 大一学生寒假项目 最终实现效果如图 一.以C语言实现个人小项目 在我们快速学完了一个高级编程语言&#xff0c;就应该写一个小项目来加以巩固自己的学习成果。 所以今天&#xff0c;我们来尝试写一写greedy snake&#xff0c;对于大学生来说也是可以加强能…

Gateway API 实践之(七)FSM Gateway 的负载均衡算法

FSM Gateway 流量管理策略系列&#xff1a; 故障注入黑白名单访问控制限速重试会话保持健康检查负载均衡算法TLS 上游双向 TLS 在微服务和 API 网关架构中&#xff0c;负载均衡是至关重要的&#xff0c;它确保每个服务实例都能平均地处理请求&#xff0c;同时也为高可用性和故…

python 基础知识点(蓝桥杯python科目个人复习计划34)

今日复习内容&#xff1a;以做题为主 例题1&#xff1a;Alice 和 Bob的爱恨情仇 题目描述&#xff1a; Alice和Bob最近正在学习博弈论&#xff0c;为了学以致用&#xff0c;他们找来了一大堆的小饼干 &#xff0c;并通过博弈的方式来吃掉这些小饼干。他们轮流对这些饼干进行…

随机森林超参数的网格优化(机器学习的精华--调参)

随机森林超参数的网格优化&#xff08;机器学习的精华–调参&#xff09; 随机森林各个参数对算法的影响 影响力参数⭐⭐⭐⭐⭐几乎总是具有巨大影响力n_estimators&#xff08;整体学习能力&#xff09;max_depth&#xff08;粗剪枝&#xff09;max_features&#xff08;随机…

mysql学习打卡day22

今日成果&#xff1a; select * from employees where salary > (select avg(salary) from employees); -- 查询超过平均工资的员工select * from clients where client_id not in (select distinct client_id from invoices); -- 查询没有发票的用户 感谢各位读者查阅&…

ASP.NET Core 预防开放式重定向攻击

写在前面 为预防钓鱼网站的常用套路&#xff0c;在进行 Web 应用程序的开发时&#xff0c;原则上应该将所有由用户提交的数据视为不可信。如果应用程序中包含了基于 URL 内容重定向的功能&#xff0c;需要确保这种类型的重定向操作只能在应用本地完成&#xff0c;或者明确判断…

[技术杂谈]如何下载vscode历史版本

网站模板&#xff1a; https://code.visualstudio.com/updates/v1_85 如果你想下载1.84系列可以访问https://code.visualstudio.com/updates/v1_84​​​​​​ 然后看到&#xff1a; 选择对应版本下载即可&#xff0c;我是windows x64系统选择x64即可开始下载

Oracle11g安装配置详细教程

Oracle 11g的安装和配置是一个详尽的过程&#xff0c;涉及多个步骤。以下是基于Oracle 11g R2&#xff08;第二个发布版本&#xff09;在Windows环境下的安装和基本配置的概述。请注意&#xff0c;Oracle 11g在不同的操作系统上的安装步骤可能会有所不同&#xff0c;此处以Wind…

MQTT在linux下服务端和客户端的应用

MQTT&#xff08;Message Queuing Telemetry Transport&#xff09;是一种轻量级、开放标准的消息传输协议&#xff0c;设计用于受限设备和低带宽、不稳定网络的通信。 MQTT的一些关键特点和概念&#xff1a; 发布/订阅模型&#xff1a; MQTT采用发布/订阅&#xff08;Publ…

QCustomplot实现灰度曲线图

从 QCustomplot官网 https://www.qcustomplot.com/index.php/download 下载支持文件。首页有些demo可以进行参考学习。 新建一个Qt工程&#xff0c;将下载得到的qcustomplot.h和qcustomplot.cpp文件加入到当前工程。pro文件中加上 printsupport 在ui界面中&#xff0c;添加一…