我没有得到预期的P和Q正定矩阵。请帮我找到错误后拖lmi)

1视图(30天)
“* P + P * - 2 * P * * B 'P < 0, P > 0 - > (1)
A_new‘* Q + Q * A_new - 2 * C_new”* C_new < 0, Q > 0 - > (2)
(0 = 1;1 0]
B = (1, 1)
A_new = [0 1 1;1 0 1;0 0 0)
C_new = [1 0 0]
我写了第一LMI的代码
setlmis ([])
P = lmivar([大小(1,1)]);
lmiterm ([1 1 1 P], 1, 1);% P > 0
lmiterm ([2 1 1 P], ', 1,“年代”);
lmiterm ([2 1 1 P] 2 B * B * P);
LMIs1 = getlmis;
[TMIN1, XFEAS1] = feasp (LMIs1);
P = dec2mat (LMIs1 XFEAS1 P)
我写了第二LMI的代码
setlmis ([])
Q = lmivar(1,[大小(A_new, 1) 1]);
lmiterm ([3 1 1 Q], 1, 1);% Q > 0
lmiterm ([4 1 1 Q], A_new ', 1,“年代”);
lmiterm ([4 1 1 1] 2 C_new”* C_new);
LMIs2 = getlmis;
[TMIN2, XFEAS2] = feasp (LMIs2);
Q = dec2mat (LMIs2, XFEAS2 Q)
1评论
Zhengqi王
Zhengqi王 2021年7月4日
我已经运行您的代码。我认为没有问题。结果P和Q正定因为所有矩阵的特征值是积极的。你可以测试它eig (P)和eig (Q)。

登录置评。

答案(0)

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!