[未答复] feasp求解的可行解代入验证发现不可行

[复制链接]
zyl2333 发表于 2022-5-5 12:37:58
A1_1=[0.8 0.35;-0.4 3.1];
A2_1=[-0.1 0.4;0.6 0.25];
B2_1=[0.1 0;0 0];
B1_1=[-0.42 0;0.15 0.4];
C1_1=[0 0.5;-0.5 0.2];
C2_1=[-0.8 0;0.5 0];
D1_1=[0 0;0.5 0];
D2_1=[0 0;0 0];
U=[0 0;0 0];
V=[-1 -1;-1 -1];
p12=0.43;

setlmis([]);     
K=lmivar(1,[1 0]);
Q=lmivar(1,[2 1]);
P1=lmivar(1,[2 1]);
P2=lmivar(1,[2 1]);
R1=lmivar(1,[2 1]);
G1=lmivar(2,[2 2]);
Y1=lmivar(2,[2 2]);
M=lmivar(1,[2 1]);
%%
lmiterm([1,1,1,R1],-1,1);
lmiterm([1,1,2,G1'],1,C1_1'*C2_1);
lmiterm([1,1,2,Y1'],1,D1_1'*C2_1);
lmiterm([1,1,3,G1'],1,C1_1'*D2_1);
lmiterm([1,1,3,Y1'],1,D1_1'*D2_1);
lmiterm([1,1,4,G1'],1,A1_1');
lmiterm([1,1,4,Y1'],1,B1_1');
lmiterm([1,1,5,G1'],1,1);
lmiterm([1,1,6,G1'],1,C1_1');
lmiterm([1,1,6,Y1'],1,D1_1');
lmiterm([1,2,2,Q],-1,1);
lmiterm([1,2,2,0],C2_1'*C2_1);
lmiterm([1,2,3,0],C2_1'*D2_1);  
lmiterm([1,2,4,0],A2_1');
lmiterm([1,3,3,K],-1,1);  
lmiterm([1,3,3,0],D2_1'*D2_1);
lmiterm([1,3,4,0],B2_1');
lmiterm([1,4,4,G1'],-1,1);
lmiterm([1,4,4,G1],-1,1);
lmiterm([1,4,4,P2],p12,1);
lmiterm([1,4,4,P1],(1-p12),1);
lmiterm([1,5,5,M],-1,1);
lmiterm([1,6,6,0],-1);
%% 矩阵变量约束      
lmisys=getlmis;
[tmin,xfeasp]=feasp(lmisys);
PP1=dec2mat(lmisys,xfeasp,P1);
PP2=dec2mat(lmisys,xfeasp,P2);
RR1=dec2mat(lmisys,xfeasp,R1);
KK=dec2mat(lmisys,xfeasp,K);
QQ=dec2mat(lmisys,xfeasp,Q);
MM=dec2mat(lmisys,xfeasp,M);
GG1=dec2mat(lmisys,xfeasp,G1);
YY1=dec2mat(lmisys,xfeasp,Y1);
Z1=[-RR1  GG1'*C1_1'*C2_1+YY1'*D1_1'*C2_1  GG1'*C1_1'*D2_1+YY1'*D1_1'*D2_1 GG1'*A1_1'+YY1'*B1_1'  GG1'  GG1'*C1_1'+YY1'*D1_1';
    C2_1'*C1_1*GG1+C2_1'*D1_1*YY1  -QQ+C2_1'*C2_1  C2_1'*D2_1  A2_1'  U  U;   
    D2_1'*C1_1*GG1+D2_1'*D1_1*YY1  D2_1'*C2_1  -KK+D2_1'*D2_1  B2_1'  U  U;
    A1_1*GG1+B1_1*YY1  A2_1  B2_1  -GG1'-GG1+p12*PP2+(1-p12)*PP1  U  U;
    GG1  U  U  U  -MM  U;
    C1_1*GG1+D1_1*YY1  U  U  U  U  V;
];
E1=eig(Z1); %求Z的特征值
if all(E1)<0
    disp('ok');
end

求出来的特征值并不都是小于零的。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表