看看https://wk.baidu.com/view/2eda0c054b2fb4daa58da0116c175f0e7cd119f6
clear
clc
A = [0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 1
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0];
B = [0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 1]';
beta = 3;
setlmis([])
P = lmivar(2, [6 6]);
lmiterm([1 1 1 P], 1, A, 's')
lmiterm([1 1 1 P], beta, 1)
lmiterm([1 1 2 P], 1, -B)
lmiterm([1 2 1 P], B', 1)
lmiterm([1 2 2 0], -1)
lmis1 = getlmis;
% solve
[tmin, xfeas] = feasp(lmis1);
P = dec2mat(lmis1, xfeas, P)