问题遇到的现象和发生背景 1、要求:A=[-11.6028 7.1632 0 0;6.4909 -27.837 0 0;10 0 -10 0;0 10 0 -10];B=[3.086;5.4458;0;0];C=[0 0 1 0;0 0 0 1];
其中的W∈Rn×n,V∈Rp×p为设计的对称正定矩阵,X=inv(P),求X的值。
下面是我的代码:(运行报错)
clc;clear;close all
A=[-11.6028 7.1632 0 0;6.4909 -27.837 0 0;10 0 -10 0;0 10 0 -10];
B=[3.086;5.4458;0;0];
C=[0 0 1 0;0 0 0 1];
c=[1 0 0 0 1 0 0 1 0 1];
W=eye();
V=eye();
setlmis([]);
X=lmivar(1,[2 1;2 1]);
P=inv(X);
lmiterm([1 1 1 P],1,A);
lmiterm([1 1 1 P],A',1);
lmiterm([1 1 1 inv(V)],-C',C);
lmiterm([1 1 2 P],1,1);
lmiterm([1 2 1 P],1,1);
lmiterm([1 2 2 inv(W)],-1,1);
lmisys=getlmis;
lmiterm([1 1 1 P],-1,1);
lmiterm([1 1 2 1],1,1);
lmiterm([1 2 1 1],1,1);
lmiterm([1 2 2 X],-1,1);
lmis=getlmis;
options=[1e-5,0,0,0,0];
[copt,xopt]=mincx(lmisys,c,option);
X=dec2mat(lmisys,xopt,X)
copt
xopt