關于定理1的求解，以下是我的代碼，求出來沒有可行解，請問哪里出問題了？
clc
clear
A=[0 1 0;0 0 1;-1 -5 -3];
Ad=0.6*eye(3);
G=0.1*eye(3);
B=[0.5 0.3 0;0 0.9 0.1;-0.5 -0.2 0.1];
E=eye(3);
A1=[zeros(6,3),[E;A]];
Ad1=[zeros(6,3),[zeros(3);Ad]];
% Ad1=0.6*eye(6);
G1=[zeros(3); G];
B1=[zeros(3,3);B];
H1=[-eye(3);zeros(3)];
KN=[-1.4443 0.4449 -0.3640 -2.2578 -0.6375 -0.9807
-0.2795 -0.4155 1.0127 -0.1189 -0.7128 0.9251
-0.5218 -1.4436 -1.1944 -0.8794 -0.0251 -0.6915];
KN=-0.5*ones(3,6);
mu=0.4;
Z1=0.1*eye(3);
Z1=[0.1 0 0;0 0.2 0;0 0 0.3];
Z2=0.1*eye(3);
Z2=[0.3 0 0;0 0.4 0;0 0 0.5];
Z=[Z1,zeros(3);zeros(3),Z2];
R=0.1*eye(3);
% x1=[1 1 1 1 1 1]';
gamma=1;
%%
setlmis([])
kappa=lmivar(1,[1,1]);
P=lmivar(2,[6,6]);
Q=lmivar(2,[6,6]);
lmiterm([-2 1 1 P],1,1);
lmiterm([-3 1 1 Q],1,1);
lmiterm([-4 1 1 kappa],1,1);
lmiterm([1 1 1 P],1,A1+B1*KN);
lmiterm([1 1 1 -P],(A1+B1*KN)',1);
lmiterm([1 1 1 Q],1,1);
lmiterm([1 1 1 0],Z+KN'*R*KN);
lmiterm([1 1 1 kappa],1,1);
lmiterm([1 1 2 P],1,Ad1);
lmiterm([1 1 3 P],1,H1);
lmiterm([1 2 2 Q],-(1-mu),1);
lmiterm([1 2 3 0],0);
lmiterm([1 3 3 0],-gamma);
%
lmis=getlmis;
[tmin,xfea]=feasp(lmis);
P=dec2mat(lmis,xfea,P)
Q=dec2mat(lmis,xfea,Q)
kappa=dec2mat(lmis,xfea,kappa)
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/66708.html
標籤：其他
上一篇：LeetCode Week 6：第 51 ~ 55 題
下一篇：系統清理