How to find the MADB h based on LMI?
A=[-2,0;0,-0.9];
B=[-1,0;-1,-1];
D=[-0.2;-0.3];
M=[0.3,0.1];
N=[0.1,0.2];
G=[0.1,0;0,0.1];
Ed=[0;0];
Ea=[1,0;0,1];
Eb=[1,0;0,1];
K1=0.2;
K2=0.5;
hd=9;
setlmis([]);
Q1=lmivar(1,[2,1]);
Z1=lmivar(1,[2,1]);
Q2=lmivar(1,[2,1]);
Z2=lmivar(1,[2,1]);
Q1=lmivar(1,[2,1]);
P=lmivar(1,[2,1]);
I=lmivar(1,[2,1]);
epi=lmivar(1, [1,0]);
miu=lmivar(1, [1,0]);
h=lmivar(1, [1,0]);
lmiterm([1,1,1,P],A’,1);
lmiterm([1,1,1,P],1,A);
lmiterm([1,1,1,P],1,-D*K1*M);
lmiterm([1,1,1,P],-M’*K1’*D’,1);
lmiterm([1,1,1,-Z1],-1,1);
lmiterm([1,1,1,Q2],1,1);
lmiterm([1,1,2,Z1],1,1);
lmiterm([1,1,3,P],1,B);
lmiterm([1,1,3,P],-1,D*K1*N);
lmiterm([1,1,4,P],1,D);
lmiterm([1,1,4,epi],1,-M’*K2′);
lmiterm([1,1,4,epi],1,M’*K1′);
lmiterm([1,1,5,Z1],(h/2)*A’,1);
lmiterm([1,1,5,Z2],(h/2)*A’,1);
lmiterm([1,1,5,Z1],-M’*K1’*D’,1);
lmiterm([1,1,5,Z2],-M’*K1’*D’,1);
lmiterm([1,1,6,P],G’,1);
lmiterm([1,1,7,miu],1,Ea);
lmiterm([1,1,7,miu],1,-Ed*K1*M);
lmiterm([1,2,2,Z1],-1,1);
lmiterm([1,2,2,Z2],-1,1);
lmiterm([1,2,2,Q1],1,1);
lmiterm([1,2,2,Q2],-1,1);
lmiterm([1,2,3,Z2],1,1);
lmiterm([1,3,3,Z2],-1,1);
lmiterm([1,3,3,Q1],-1,1);
lmiterm([1,3,4,epi],1,-N’*K2′);
lmiterm([1,3,4,epi],1,N’*K1′);
lmiterm([1,3,5,Z1],(h/2)*B’,1);
lmiterm([1,3,5,Z2],(h/2)*B’,1);
lmiterm([1,3,5,Z1],-N’*K1’*D’,1);
lmiterm([1,3,5,Z2],-N’*K1’*D’,1);
lmiterm([1,3,7,miu],1,Eb);
lmiterm([1,3,7,miu],1,-Ed*K1*N);
lmiterm([1,4,4,0],-1);
lmiterm([1,4,5,Z1],(h/2)*D’,1);
lmiterm([1,4,5,Z2],(h/2)*D’,1);
lmiterm([1,4,7,miu],1,Ed);
lmiterm([1,5,5,Z1],-1,1);
lmiterm([1,5,5,Z2],-1,1);
lmiterm([1,5,6,Z1],(h/2)*G’,1);
lmiterm([1,5,6,Z2],(h/2)*G’,1);
lmiterm([1,6,6,0],-1);
lmiterm([1,7,7,0],-1);
lmiterm([-2,1,1,Q1],1,1);
lmiterm([-3,1,1,Z1],1,1);
lmiterm([-4,1,1,Q2],1,1);
lmiterm([-5,1,1,Z2],1,1);
lmiterm([-6,1,1,P],1,1);
lmiterm([-7,1,1,epi],1,1);
lmiterm([-8,1,1,miu],1,1);
lmiterm([-9,1,1,h],1,1);
lmisys=getlmis;
[tmin, xfeas] =feasp(lmisys) ;
if(tmin<0)
string=[‘tmin= ‘,num2str(tmin)];
disp(string)
options=[1e-5, 0, 0, 0, 0];
c=mat2dec(lmisys, eye(1,1) ) ;
[copt,xopt]=mincx(lmisys , c, options);
-h==dec2mat(lmisys, xopt, -h);
endA=[-2,0;0,-0.9];
B=[-1,0;-1,-1];
D=[-0.2;-0.3];
M=[0.3,0.1];
N=[0.1,0.2];
G=[0.1,0;0,0.1];
Ed=[0;0];
Ea=[1,0;0,1];
Eb=[1,0;0,1];
K1=0.2;
K2=0.5;
hd=9;
setlmis([]);
Q1=lmivar(1,[2,1]);
Z1=lmivar(1,[2,1]);
Q2=lmivar(1,[2,1]);
Z2=lmivar(1,[2,1]);
Q1=lmivar(1,[2,1]);
P=lmivar(1,[2,1]);
I=lmivar(1,[2,1]);
epi=lmivar(1, [1,0]);
miu=lmivar(1, [1,0]);
h=lmivar(1, [1,0]);
lmiterm([1,1,1,P],A’,1);
lmiterm([1,1,1,P],1,A);
lmiterm([1,1,1,P],1,-D*K1*M);
lmiterm([1,1,1,P],-M’*K1’*D’,1);
lmiterm([1,1,1,-Z1],-1,1);
lmiterm([1,1,1,Q2],1,1);
lmiterm([1,1,2,Z1],1,1);
lmiterm([1,1,3,P],1,B);
lmiterm([1,1,3,P],-1,D*K1*N);
lmiterm([1,1,4,P],1,D);
lmiterm([1,1,4,epi],1,-M’*K2′);
lmiterm([1,1,4,epi],1,M’*K1′);
lmiterm([1,1,5,Z1],(h/2)*A’,1);
lmiterm([1,1,5,Z2],(h/2)*A’,1);
lmiterm([1,1,5,Z1],-M’*K1’*D’,1);
lmiterm([1,1,5,Z2],-M’*K1’*D’,1);
lmiterm([1,1,6,P],G’,1);
lmiterm([1,1,7,miu],1,Ea);
lmiterm([1,1,7,miu],1,-Ed*K1*M);
lmiterm([1,2,2,Z1],-1,1);
lmiterm([1,2,2,Z2],-1,1);
lmiterm([1,2,2,Q1],1,1);
lmiterm([1,2,2,Q2],-1,1);
lmiterm([1,2,3,Z2],1,1);
lmiterm([1,3,3,Z2],-1,1);
lmiterm([1,3,3,Q1],-1,1);
lmiterm([1,3,4,epi],1,-N’*K2′);
lmiterm([1,3,4,epi],1,N’*K1′);
lmiterm([1,3,5,Z1],(h/2)*B’,1);
lmiterm([1,3,5,Z2],(h/2)*B’,1);
lmiterm([1,3,5,Z1],-N’*K1’*D’,1);
lmiterm([1,3,5,Z2],-N’*K1’*D’,1);
lmiterm([1,3,7,miu],1,Eb);
lmiterm([1,3,7,miu],1,-Ed*K1*N);
lmiterm([1,4,4,0],-1);
lmiterm([1,4,5,Z1],(h/2)*D’,1);
lmiterm([1,4,5,Z2],(h/2)*D’,1);
lmiterm([1,4,7,miu],1,Ed);
lmiterm([1,5,5,Z1],-1,1);
lmiterm([1,5,5,Z2],-1,1);
lmiterm([1,5,6,Z1],(h/2)*G’,1);
lmiterm([1,5,6,Z2],(h/2)*G’,1);
lmiterm([1,6,6,0],-1);
lmiterm([1,7,7,0],-1);
lmiterm([-2,1,1,Q1],1,1);
lmiterm([-3,1,1,Z1],1,1);
lmiterm([-4,1,1,Q2],1,1);
lmiterm([-5,1,1,Z2],1,1);
lmiterm([-6,1,1,P],1,1);
lmiterm([-7,1,1,epi],1,1);
lmiterm([-8,1,1,miu],1,1);
lmiterm([-9,1,1,h],1,1);
lmisys=getlmis;
[tmin, xfeas] =feasp(lmisys) ;
if(tmin<0)
string=[‘tmin= ‘,num2str(tmin)];
disp(string)
options=[1e-5, 0, 0, 0, 0];
c=mat2dec(lmisys, eye(1,1) ) ;
[copt,xopt]=mincx(lmisys , c, options);
-h==dec2mat(lmisys, xopt, -h);
end A=[-2,0;0,-0.9];
B=[-1,0;-1,-1];
D=[-0.2;-0.3];
M=[0.3,0.1];
N=[0.1,0.2];
G=[0.1,0;0,0.1];
Ed=[0;0];
Ea=[1,0;0,1];
Eb=[1,0;0,1];
K1=0.2;
K2=0.5;
hd=9;
setlmis([]);
Q1=lmivar(1,[2,1]);
Z1=lmivar(1,[2,1]);
Q2=lmivar(1,[2,1]);
Z2=lmivar(1,[2,1]);
Q1=lmivar(1,[2,1]);
P=lmivar(1,[2,1]);
I=lmivar(1,[2,1]);
epi=lmivar(1, [1,0]);
miu=lmivar(1, [1,0]);
h=lmivar(1, [1,0]);
lmiterm([1,1,1,P],A’,1);
lmiterm([1,1,1,P],1,A);
lmiterm([1,1,1,P],1,-D*K1*M);
lmiterm([1,1,1,P],-M’*K1’*D’,1);
lmiterm([1,1,1,-Z1],-1,1);
lmiterm([1,1,1,Q2],1,1);
lmiterm([1,1,2,Z1],1,1);
lmiterm([1,1,3,P],1,B);
lmiterm([1,1,3,P],-1,D*K1*N);
lmiterm([1,1,4,P],1,D);
lmiterm([1,1,4,epi],1,-M’*K2′);
lmiterm([1,1,4,epi],1,M’*K1′);
lmiterm([1,1,5,Z1],(h/2)*A’,1);
lmiterm([1,1,5,Z2],(h/2)*A’,1);
lmiterm([1,1,5,Z1],-M’*K1’*D’,1);
lmiterm([1,1,5,Z2],-M’*K1’*D’,1);
lmiterm([1,1,6,P],G’,1);
lmiterm([1,1,7,miu],1,Ea);
lmiterm([1,1,7,miu],1,-Ed*K1*M);
lmiterm([1,2,2,Z1],-1,1);
lmiterm([1,2,2,Z2],-1,1);
lmiterm([1,2,2,Q1],1,1);
lmiterm([1,2,2,Q2],-1,1);
lmiterm([1,2,3,Z2],1,1);
lmiterm([1,3,3,Z2],-1,1);
lmiterm([1,3,3,Q1],-1,1);
lmiterm([1,3,4,epi],1,-N’*K2′);
lmiterm([1,3,4,epi],1,N’*K1′);
lmiterm([1,3,5,Z1],(h/2)*B’,1);
lmiterm([1,3,5,Z2],(h/2)*B’,1);
lmiterm([1,3,5,Z1],-N’*K1’*D’,1);
lmiterm([1,3,5,Z2],-N’*K1’*D’,1);
lmiterm([1,3,7,miu],1,Eb);
lmiterm([1,3,7,miu],1,-Ed*K1*N);
lmiterm([1,4,4,0],-1);
lmiterm([1,4,5,Z1],(h/2)*D’,1);
lmiterm([1,4,5,Z2],(h/2)*D’,1);
lmiterm([1,4,7,miu],1,Ed);
lmiterm([1,5,5,Z1],-1,1);
lmiterm([1,5,5,Z2],-1,1);
lmiterm([1,5,6,Z1],(h/2)*G’,1);
lmiterm([1,5,6,Z2],(h/2)*G’,1);
lmiterm([1,6,6,0],-1);
lmiterm([1,7,7,0],-1);
lmiterm([-2,1,1,Q1],1,1);
lmiterm([-3,1,1,Z1],1,1);
lmiterm([-4,1,1,Q2],1,1);
lmiterm([-5,1,1,Z2],1,1);
lmiterm([-6,1,1,P],1,1);
lmiterm([-7,1,1,epi],1,1);
lmiterm([-8,1,1,miu],1,1);
lmiterm([-9,1,1,h],1,1);
lmisys=getlmis;
[tmin, xfeas] =feasp(lmisys) ;
if(tmin<0)
string=[‘tmin= ‘,num2str(tmin)];
disp(string)
options=[1e-5, 0, 0, 0, 0];
c=mat2dec(lmisys, eye(1,1) ) ;
[copt,xopt]=mincx(lmisys , c, options);
-h==dec2mat(lmisys, xopt, -h);
end lmi, optimization MATLAB Answers — New Questions