ÈçÌ⣬LMIÖÐÒªÇó¶Ô½ÇÏߵľØÕó±ØÐëÊÇ·½Õó²ÅÄÜÇó½â£¬¿ÉÊÇÎÒÒªÇóµÄ¾ØÕó²»ÊÇ·½ÕóÔõô°ì£¿
ÎÊÌâÈçͼ1.
Z=[98.235;50.862;49.144;48.756;25.397;25.514;75.375];
A=[1 -1 -1 0 0 0 0;0 0 1 -1 0 0 0;0 1 0 0 -1 -1 0;0 0 0 1 1 0 -1];
e1=[0.5;0.5;0.5;0.5;0.5;0.5;0.5];
e2=[0.4;0.4;0.4;0.4;0.4;0.4;0.4];
r1=[0.08;0.08;0.08;0.08;0.08;0.08;0.08];
r2=[0.06;0.06;0.06;0.06;0.06;0.06;0.06];
AΪ4X7¾ØÕó£»X£¬Z£¬e1£¬e2£¬r1£¬r2¾ùΪ7x1¾ØÕó
´úÂëÈçÏ£º
Z1=z-e1; Z2=-(z+e2);
setlmis([]);
X=lmivar(2,[7 1]);
lmiterm([1 1 1 X],-1,1);
lmiterm([1 1 1 0], Z1);
lmiterm([2 1 1 X],1,1);
lmiterm([2 1 1 0],Z2);
lmiterm([3 1 1 X], A,1);
lmiterm([3 1 1 0],-r1);
lmiterm([4 1 1 X],-A,1);
lmiterm([4 1 1 0],r2);
lmis=getlmis;
ÔËÐкóÈçÏ´íÎó£º
lmiterm([1 1 1 X],-1,1);
lhs of LMI #1, block (1,1): terms in diagonal blocks must be square
MatlabÒªÇóÿһ¸ö²»µÈʽ±ØÐëÊÇ·½Õ󣬼´X£¬AXÒ²±ØÐëÊÇ·½Õó¡£Ôõô°ì£¿²»¿ÉÄÜÿ´Î¶¼ÕýºÃÊÇ·½Õóѽ¡£
Èç¹ûת»¯³ÉͼƬ2ÔÙÇó½â£¬
³ÌÐòÈçÏÂ
Z1=z-e1; Z2=-(z+e2);
setlmis([]);
X=lmivar(2,[7 1]);
lmiterm([1 1 1 X],-1,1);
lmiterm([1 1 1 0],Z1);
lmiterm([1 2 2 X],1,1);
lmiterm([1 2 2 0],Z2);
lmiterm([1 3 3 X], A,1);
lmiterm([1 3 3 0],r1);
lmiterm([1 4 4 X],A,-1);
lmiterm([1 4 4 0],r2);
lmis=getlmis;
»¹ÊdzöÏÖͬÑùµÄÎÊÌâ¡£
LMIÒªÇó²»µÈʽб¶Ô½ÇÏßÉϵľØÕóΪ·½Õ󣬼´nxn¾ØÕó¡£¶øÎÒÕâÀïÒªÇó½âµÄ²»µÈʽ¾ØÕó£¨1£©-£¨4£©¶¼²»ÊÇ·½ÕóÔõô°ì£¿ÒªÔõô´¦Àí£¿
ÓÃfminconº¯Êý¿ÉÒԽ⣬µ«ÊǾÍÊÇÏëѧ»áʹÓÃLMIº¯Êý£¬Èç¹û¶¼ÊÇ·½Õó²Å¿ÉÒԽ⣬ÄÇÆñ²»ÊÇLMIÕâ¸ö¹¤¾ßÏäºÜ¾ÖÏÞ£¿
![100RMBÇóMATLAB¼òµ¥µÄLMIÎÊÌ⣡]()
111111.png
![100RMBÇóMATLAB¼òµ¥µÄLMIÎÊÌ⣡-1]()
22222.png |