【转载】Matlab中LMI(线性矩阵不等式)工具箱使用教程

@TOC
转载
原文地址：https://www.cnblogs.com/Hand-Head/articles/5412511.html

这一段被老板逼着论文开题，自己找方向比较着急，最后选择了供应链控制理论的一个方向。我要写的论文，用到了Matlab的LMI工具，以及某篇论文中的H-inf稳定定理。自己好好研究了好长时间，怎么也无法实现该论文当中的算例。研究了一个多月，自己简直就快崩溃了，也搞不定问题。我很是怀疑自己的选题是不是正确，并且怀疑自己是不是选的太难了。如果连论文中的算例都无法实现，如何把该模型应用到自己论文当中呢？功夫不负有心人，昨日我加入了Mathworks的Matlab的Newsgroup，结果碰见一牛人Johan，立即就把论文中的算例给写成程序。但是他做出的结果和论文仍然有差别，我仍有点不甘心，人家的论文发表在Automatica上，如果连这种期刊都水的要命，那么就没有什么学术水平了。

今天中午，仍然不甘心，老爸给我打了电话让我看红场阅兵，于是我边看PPMate边漫无边际的核对着自己的程序。终于做出了和算例一致的结果。

我搜出来的都是一些简单的算例，并且机会没有中文教程，我在这里就斗胆把自己的体会写出来，试着给大家提供一点参考。

LMI：Linear Matrix Inequality，就是线性矩阵不等式。
在Matlab当中，我们可以采用图形界面的lmiedit命令，来调用GUI接口，但是我认为采用程序的方式更方便（也因为我不懂这个lmiedit的GUI）。
对于LMI Lab， 其中有三种求解器（solver）： feasp，mincx和gevp。
每个求解器针对不同的问题：
feasp：解决可行性问题（feasibility problem），例如：A(x)

mincx：在线性矩阵不等式的限制下解决最小化问题（Minimization of a linear objective under LMI constraints），例如最小化c’x，在限制条件A(x) < B(x)下。

gevp：解决广义特征值最小化问题。例如：最小化lambda，在0

要解决一个LMI问题，首要的就是要把线性矩阵不等式表示出来。

对于以下类型的任意的LMI问题

N’ * L(X1, . . . , XK) * N < M’ * R(X1, . . . , XK) * M

其中X1, . . . , XK是结构已经事先确定的矩阵变量。左侧和右侧的外部因子（outer factors）N和M是给定的具有相同维数的矩阵。
左侧和右侧的内部因子（inner factors）L(.)和R(.)是具有相同结构的对称块矩阵。每一个块由X1, . . . , XK以及它们的转置组合而成形成的。

解决LMI问题的步骤有两个：

1、定义维数以及每一个矩阵的结构，也就是定义X1, . . . , XK。

2、描述每一个LMI的每一项内容（Describe the term content of each LMI）

此处介绍两个术语：

矩阵变量（Matrix Variables）：例如你要求解X满足A(x)

变量。

项（Terms）：项是常量或者变量（Terms are either constant or variable）。

常项（Constant Terms）是确定的矩阵。可变项（Variable Terms)是哪些含有矩阵变

量的项，例如：XA, XC’。如果是XA + XC’，那么记得要把它当成两项来处理。

好了废话不说了，让我们来看个例子吧（下面是一线性时滞系统）。

针对这个式子，如果存在满足如下LMI的正矩阵(positive-define)的Q，S1，S2和矩阵M，那么我们就称作

该系统为H-inf渐进稳定的，并且gammar是上限。

该论文的地址为：论文原文地址
H∞ state feedback control for generalized continuous/discrete time-delay system

该论文的算例为：

我们要实现的就利用LMI进行求解，验证论文结果。

首先我们要用setlmis([])命令初始化一个LMI系统。

接下来，我们就要设定矩阵变量了。采用函数为lmivar

语法：X = lmivar(type,struct)

type=1: 定义块对角的对称矩阵。 每一个对角块或者是全矩阵<任意对称矩阵>，标量<单位矩阵的乘积>，或者是零阵。

如果X有R个对角块，那么后面这个struct就应该是一个Rx2阶的的矩阵，在此矩阵中，struct(r,1)表示第r个块的大小，struct(r,2) 表示第r个块的类型<1–全矩阵，0–标量，-1–零阵）。

比如一个矩阵有两个对角块，其中一个是2x2的全对称矩阵，第二个是1x1的一个标量，那么该矩阵变量应该表示为X = lmivar(1, [2 1; 1 0]) 。

type=2: mxn阶的矩阵，只需要写作struct = [m,n]即可。

type=3: 其它类型。针对类型3，X的每一个条目（each entry of X）被定义为0或者是+(-)xn，此处xn代表了第n个决策变量。

那么针对我们的例子，我们如此定义变量：

% Q is a symmetric matrix, has a block size of 2 and this block is symmetric
Q = lmivar(1, [2 1]);

% S1 a symmeric matrix, size 2
S1 = lmivar(1, [2 1]);

% S2 is 1 by 1 matrix
S2 = lmivar(1, [1 0]);

% Type of 2, size 1 by 2
M = lmivar(2, [1 2]);

定义完成变量之后，我们就该用lmiterm来描述LMI中的每一个项了。Matlab的官方文档提示我们，如果要描述一个LMI只需要描述上三角或者下三角元素就可以了，否则会描述成另一个LMI。

When describing an LMI with several blocks, remember to specify only the terms in the blocks on or below the diagonal (or equivalently, only the terms in blocks on or above the diagonal).

语法为：lmiterm(termID,A,B,flag)

termID是一个四维整数向量，来表示该项的位置和包含了哪些矩阵变量。

termID(1)可以为+p或者-p，+p代表了这个项位于第p个线性矩阵不等式的左边，-p代表了这个项位于第p个线性矩阵不等式的右边。注意：按照惯例来讲，左边通常指较小的那边。

termID(2:3):

1、对于外部变量来说，取值为[0,0];

2、对于左边或者右边的内部变量来说，如果该项在(i,j)位置，取值[i,j]

termID(4):

1、对于外部变量，取值为0

2、对于AXB，取值X

3、对于A*X’*B，取值-X

flag(可选，值为s）:

因为：(AXB) + (AXB)T = AXB + B’*X’*A’，所以采用s来进行简写。

比如：针对A*X + X’*A’

我们采用笨方法：

lmiterm([1 1 1 X],A,1)
lmiterm([1 1 1 -X],1,A’)

那么简写就是lmiterm([1 1 1 X],A,1,‘s’)

接下来我们就看该论文中的算例吧：(1,1)位置是

-Q+BdS2Bd’+AdS1Ad’;

我们应该表示为：

% pos in (1, 1)
lmiterm([1 1 1 Q], -1, 1);
lmiterm([1 1 1 S2], Bd, Bd’);
lmiterm([1 1 1 S1], Ad, Ad’);

其它位置仿照写就行了，不懂了多看帮助文档。

把每一个项都定义以后，要记得

lmis = getlmis;

[tmin, feas] = feasp(lmis)

getlmis:是在完成定义变量和项之后，LMI系统的内部表示就可以通过此命令获得（After completing the description of a given LMI system with lmivar and lmiterm, its internal representation lmisys is obtained with the command）。

feasp是调用feasp求解器，看有没有可行解。feas就是可行解。

下面我把代码贴上去，那些常数矩阵都在此源程序中定义了。

A = [2 1; 0 1];
Ad = [0.2 0.1; 0 0.1];
B1 = [0.1 0.1]’;
B2 = [1 1]’;
Bd = [0.1 0.1]’;

C = [1, 1];
Cd = [0.1, 0.1];

D11 = 0.1;
D12 = 1;
Dd = 0.1;

gammar = 1;

% Initial a LMI system
setlmis([]);

% Define Variables

% Q is a symmetric matrix, has a block size of 2 and this block is symmetric
Q = lmivar(1, [2 1]);

% S1 a symmeric matrix, size 2
S1 = lmivar(1, [2 1]);

% S2 is 1 by 1 matrix
S2 = lmivar(1, [1 0]);

% Type of 2, size 1 by 2
M = lmivar(2, [1 2]);

% Q, S1, S2 > 0
lmiterm([-2 1 1 Q], 1, 1);

lmiterm([-3 1 1 S1], 1, 1);

lmiterm([-4 1 1 S2], 1, 1);

% pos in (1, 1)
lmiterm([1 1 1 Q], -1, 1);
lmiterm([1 1 1 S2], Bd, Bd’);
lmiterm([1 1 1 S1], Ad, Ad’);

% pos (1, 2)
lmiterm([1 1 2 Q], A, 1);
lmiterm([1 1 2 M], B2, 1);

% pos(1, 3)
lmiterm([1 1 3 0], B1);

% pos(1, 4)
lmiterm([1 1 4 S2], Bd, Dd’);
lmiterm([1 1 4 S1], Ad, Cd’);

% pos(2, 2)
lmiterm([1 2 2 Q], -1, 1);

% pos(2, 4)
lmiterm([1 2 4 Q], 1, C’);

lmiterm([1 2 4 -M], 1, D12’);

% pos(2, 5)
lmiterm([1 2 5 -M], 1, 1);

% pos(2, 6)
lmiterm([1 2 6 Q], 1, 1);

% pos(3, 3)
lmiterm([1 3 3 0], -(gammar^2));

% pos(3, 4)
lmiterm([1 3 4 0], D11’);

% pos(4, 4)
lmiterm([1 4 4 0], -1);
lmiterm([1 4 4 S1], Cd, Cd’);
lmiterm([1 4 4 S2], Dd, Dd’);

lmiterm([1 5 5 S2], -1, 1);

lmiterm([1 6 6 S1], -1, 1);

lmis = getlmis;

[tmin, feas] = feasp(lmis)

运行后，就调用dec2mat把决策变量转化为矩阵形式。

Q = dec2mat(lmis, feas, Q)

Q =

1.9253   -2.2338

-2.2338 9.1054

可以看到，和论文中的一样。

欢迎使用Markdown编辑器

你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Markdown编辑器, 可以仔细阅读这篇文章，了解一下Markdown的基本语法知识。

新的改变

我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客：

1. 全新的界面设计 ，将会带来全新的写作体验；
2. 在创作中心设置你喜爱的代码高亮样式，Markdown 将代码片显示选择的高亮样式 进行展示；
3. 增加了 图片拖拽 功能，你可以将本地的图片直接拖拽到编辑区域直接展示；
4. 全新的 KaTeX数学公式 语法；
5. 增加了支持甘特图的mermaid语法¹ 功能；
6. 增加了 多屏幕编辑 Markdown文章功能；
7. 增加了 焦点写作模式、预览模式、简洁写作模式、左右区域同步滚轮设置 等功能，功能按钮位于编辑区域与预览区域中间；
8. 增加了 检查列表 功能。

功能快捷键

撤销：Ctrl/Command + Z
重做：Ctrl/Command + Y
加粗：Ctrl/Command + B
斜体：Ctrl/Command + I
标题：Ctrl/Command + Shift + H
无序列表：Ctrl/Command + Shift + U
有序列表：Ctrl/Command + Shift + O
检查列表：Ctrl/Command + Shift + C
插入代码：Ctrl/Command + Shift + K
插入链接：Ctrl/Command + Shift + L
插入图片：Ctrl/Command + Shift + G
查找：Ctrl/Command + F
替换：Ctrl/Command + G

合理的创建标题，有助于目录的生成

直接输入1次#，并按下space后，将生成1级标题。
输入2次#，并按下space后，将生成2级标题。
以此类推，我们支持6级标题。有助于使用TOC语法后生成一个完美的目录。

如何改变文本的样式

强调文本 强调文本

加粗文本 加粗文本

标记文本

删除文本

引用文本

H₂O is是液体。

2¹⁰ 运算结果是 1024.

插入链接与图片

链接: link.

图片:

带尺寸的图片:

居中的图片:

居中并且带尺寸的图片:

当然，我们为了让用户更加便捷，我们增加了图片拖拽功能。

如何插入一段漂亮的代码片

去博客设置页面，选择一款你喜欢的代码片高亮样式，下面展示同样高亮的 代码片.

// An highlighted block
var foo = 'bar';

生成一个适合你的列表

• 项目
  • 项目
    • 项目

1. 项目1
2. 项目2
3. 项目3

• 计划任务
• 完成任务

创建一个表格

一个简单的表格是这么创建的：

项目	Value
电脑	$1600
手机	$12
导管	$1

设定内容居中、居左、居右

使用:---------:居中
使用:----------居左
使用----------:居右

第一列	第二列	第三列
第一列文本居中	第二列文本居右	第三列文本居左

SmartyPants

SmartyPants将ASCII标点字符转换为“智能”印刷标点HTML实体。例如：

TYPE	ASCII	HTML
Single backticks	'Isn't this fun?'	‘Isn’t this fun?’
Quotes	"Isn't this fun?"	“Isn’t this fun?”
Dashes	-- is en-dash, --- is em-dash	– is en-dash, — is em-dash

创建一个自定义列表

Markdown

Text-to- HTML conversion tool

Authors

John

Luke

如何创建一个注脚

一个具有注脚的文本。²

注释也是必不可少的

Markdown将文本转换为 HTML。

KaTeX数学公式

您可以使用渲染LaTeX数学表达式 KaTeX:

Gamma公式展示 $\Gamma (n)=(n-1)!\forall n\in \mathbb{N}$ 是通过欧拉积分

$$\Gamma (z)={\int }_0^{\infty }{t}^{z-1}{e}^{-t}dt.$$

你可以找到更多关于的信息 LaTeX 数学表达式here.

新的甘特图功能，丰富你的文章

Mon 06 Mon 13 Mon 20 已完成 进行中 计划一 计划二 现有任务 Adding GANTT diagram functionality to mermaid

• 关于 甘特图 语法，参考 这儿,

UML 图表

可以使用UML图表进行渲染。 Mermaid. 例如下面产生的一个序列图：

张三 李四 王五 你好！李四, 最近怎么样? 你最近怎么样，王五？ 我很好，谢谢! 我很好，谢谢! 李四想了很长时间, 文字太长了 不适合放在一行. 打量着王五... 很好... 王五, 你怎么样? 张三 李四 王五

这将产生一个流程图。:

链接

长方形

圆

圆角长方形

菱形

• 关于 Mermaid 语法，参考 这儿,

FLowchart流程图

我们依旧会支持flowchart的流程图：

Created with Raphaël 2.2.0 开始 我的操作 确认？ 结束 yes no

• 关于 Flowchart流程图 语法，参考 这儿.

导出与导入

导出

如果你想尝试使用此编辑器, 你可以在此篇文章任意编辑。当你完成了一篇文章的写作, 在上方工具栏找到 文章导出 ，生成一个.md文件或者.html文件进行本地保存。

导入

如果你想加载一篇你写过的.md文件，在上方工具栏可以选择导入功能进行对应扩展名的文件导入，
继续你的创作。

---

1. mermaid语法说明 ↩︎

2. 注脚的解释 ↩︎