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

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

例题

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

问题

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

求满足如下关系的矩阵 $P$
$$\{\begin{array}{rl}& 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{array}$$

解决代码

% 段广仁老师 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)

代码解析

setlmis([])

定义一个矩阵不等式系统

lmivar()

定义一个矩阵变量

lmiterm()

描述不等式的某一项

getlmis

矩阵不等式系统定义完毕

feasp / mincx / gevp

求解器

dec2mat()

参考

1. MATLAB——LMI工具箱学习笔记
2. MATLAB——LMI工具箱学习笔记（二）
3. matlab–LMI工具箱使用学习1（lmivar()、lmiterm()函数）
4. matlab–LMI工具箱使用学习2（LMI求解）
5. Minimize Linear Objectives under LMI Constraints-MathWorks
6. LMI Applications-MathWorks