【Matlab】LMI 求解器快速入门，常用指令

Zhao-Jichao

已于 2022-10-06 17:35:09 修改

阅读量4.2k

点赞数 7

分类专栏： Simulink/Matlab 文章标签： Nash NN Optimal

于 2022-03-25 19:37:12 首次发布

本文链接：https://blog.csdn.net/weixin_36815313/article/details/123743540

版权

Simulink/Matlab 专栏收录该内容

49 篇文章

订阅专栏

本文通过一个具体的例子介绍了如何使用Matlab中的LMI工具箱来解决矩阵不等式问题，包括设置矩阵变量、定义不等式约束以及求解过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >


1	【Matlab】LMI 求解器快速入门，常用指令
2	【Matlab】利用 LMI 解矩阵不等式方程

例题

首先以段广仁老师的一个例题作为问题切入点。

问题

给定矩阵
$A_1 = \left[\begin{matrix} -1 & 2 \\ 1 & -3 \\ \end{matrix}\right], A_2 = \left[\begin{matrix} -0.8 & 1.5 \\ 1.3 & -2.7 \\ \end{matrix}\right], A_3 = \left[\begin{matrix} -1.4 & 0.9 \\ 0.7 & -2.0 \\ \end{matrix}\right],$

求满足如下关系的矩阵 $P$
$\left\{\begin{aligned} &A_1^\text{T} P + P A_1 < 0 \\ &A_2^\text{T} P + P A_2 < 0 \\ &A_3^\text{T} P + P A_3 < 0 \\ &I < P \\ \end{aligned}\right.$

解决代码

% 段广仁老师 LMI P420
% Example B.1
clear
clc

A1 =[-1  2 
      1 -3];
A2 =[-0.8  1.5 
      1.3 -2.7];
A3 =[-1.4  0.9 
      0.7 -2.0];

setlmis([])
P = lmivar(1, [2 1]);

% #1  A_1' P + P A_1
lmiterm([1 1 1 P], 1, A1, 's');
% #2  A_2' P + P A_2
lmiterm([2 1 1 P], 1, A2, 's');
% #3  A_3' P + P A_3
lmiterm([3 1 1 P], 1, A3, 's');
% #4  P
lmiterm([-4 1 1 P], 1, 1);
% #4  I
lmiterm([ 4 1 1 0], 1);

lmis1 = getlmis;

% solve
[tmin, xfeas] = feasp(lmis1);
P = dec2mat(lmis1, xfeas, P)