你好,
lmiterm我认为有一些错误,我更喜欢使用gui编辑器”lmiedit“它有助于避免错误,用它我可以产生相同的结果在书中。
%的定义矩阵
%的工作区
clc
清除所有
%指定参数矩阵
一个= [1 0 1;
2 0 1;
2 0 3)
广告= [1 0 1;
2 1 1;
0 0 1]
B = [1 1;
1 - 2;
0 1]
d = 0.1;
% LMI问题定义
% %决策变量的定义
setlmis ([]);
X = lmivar (1、1 [3]);
W = lmivar (2 [2、3]);
β= lmivar ([1]);
% LMI的定义
% LMI # 1
% LMI # 1 LMI (1,1)
lmiterm ([1 1 1 X], 1, ', ' s ');% LMI # 1: X * + * X
lmiterm ([1 1 1 X] 1广告”、“s”);% LMI # 1: X *广告' +广告* X
lmiterm ([1 1 1 W), B, 1, ' s ');% LMI # 1: B * W + W * B的
lmiterm ([1 1 1 0] d * * Ad”);% LMI # 1: d * *广告广告”
lmiterm ([1 2 1 X]、d *, 1);% LMI # 1: d * * X
lmiterm ([1 2 1 W] d B * 1);% LMI # 1: d * B * W
lmiterm(1 2 2β,5 * d,黑眼圈(3),' s ');% LMI # 1:β- d * *眼(3)(非对称吗?)
lmiterm (3 1 X [1], d *广告,1);% LMI # 1: d *广告* X
lmiterm(1 3 3β,5 * d,眼睛(3),' s ');% LMI # 1: d *β*眼(3)(非对称吗?)
lmiterm ([1 3 3 0] - d *眼(3));% LMI # 1: - d *眼(3)
lmiterm ([2 1 1 X], 1, 1);% LMI # 2: X
lmiterm(β3 1 1,1,1);% LMI # 3:β
lmiterm(β4 1 1,1,1);% LMI # 4:β
lmiterm ([4 1 1 0), 1);% LMI # 4: 1
stabilz2 = getlmis;
% LMI的解决方案
[tmin, xfeas] = feasp (stabilz2);
Xvalue = dec2mat (stabilz2 xfeas X)
Wvalue = dec2mat (stabilz2 xfeas W)
Betavalue = dec2mat (stabilz2 xfeasβ)
k = Wvalue *发票(Xvalue)
