LMI实例2

%1．定义线性矩阵不等式约束
setlmis([]);
A=[-1 -2 1;3 2 1;1 -2 -1];
B=[1;0;1];
Q=[1 -1 0;-1 -3 -12;0 -12 -36];
X=lmivar(1,[3 1]); % 变量X，满对称的
lmiterm([1 1 1 X],1,A,'s');
lmiterm([1 1 1 0],Q);
lmiterm([1 2 2 0],-1);
lmiterm([1 2 1 X],B',1);
LMIs=getlmis;
c=mat2dec(LMIs,eye(3));%将目标函数写成，其中)(TraceXxcTx是矩阵变量X中的独立元所构成的向量。
                        %由于引进向量的目的是要选择cX的对角元，因此它可以作为相应于IX=的决策向量得到，即
options=[1e-5,0,0,0,0];
>> [copt,xopt]=mincx(LMIs,c,options);%调用mincx计算最小值xopt，目标函数的全局最小值copt=c’*xopt,
                        %其中1e-5给定了所要求的关于copt的计算精度。作为求解器mincx运行的结果，以下的信息将出现在屏幕上
                         Solver for linear objective minimization under LMI constraints 

 Iterations   :    Best objective value so far 
 
     1
     2                  -8.511476
     3                 -13.063640
***                 new lower bound:   -34.023978
     4                 -15.768450
***                 new lowe