【正文】 % LMI 1: inv(ew)*Ilmiterm([2 1 1 x],1,1)。 % LMI 2: xlmiterm([3 1 1 ew],1,1)。 % LMI 3: ewss=getlmis。[lp kp]=feasp(ss)。x=dec2mat(ss,kp,x)y=dec2mat(ss,kp,y)ew=dec2mat(ss,kp,ew)p=inv(x)k=y*inv(x)附錄 2MATLAB程序如下所示：a=[ 。 ]。a1=[ 。 ]。q=[1 0。0 1]。b=[。]。e=[ 。0 ]。f=[ 。 ]。I=[1 0。0 1]。setlmis([])。x=lmivar(1,[2 1])。y=lmivar(2,[1 2])。ew=lmivar(2,[1 1])。 lmiterm([1 1 1 x],1,a39。,39。s39。)。 % LMI 1: x*a39。+a*xlmiterm([1 1 1 y],b,1,39。s39。)。 % LMI 1: b*y+y39。*b39。lmiterm([1 2 1 0],a139。)。 % LMI 1: a139。lmiterm([1 2 2 0],q)。 % LMI 1: qlmiterm([1 3 1 0],e39。)。 % LMI 1: e39。lmiterm([1 3 3 ew],.5*1,I,39。s39。)。 % LMI 1: ew*I (NON SYMMETRIC?)lmiterm([1 4 1 x],f,1)。 % LMI 1: f*xlmiterm([1 4 4 0],inv(ew)*I)。 % LMI 1: inv(ew)*Ilmiterm([1 5 1 x],1,1)。 % LMI 1: xlmiterm([1 5 5 0],inv(q))。 % LMI 1: inv(q)lmiterm([2 1 1 x],1,1)。 % LMI 2: xlmiterm([3 1 1 ew],1,1)。 % LMI 3: ewss1=getlmis。[lp kp]=feasp(ss1)。x=dec2mat(ss1,kp,x)y=dec2mat(ss1,kp,y)ew=dec2mat(ss1,kp,ew)p=inv(x)k=y*inv(x)附錄 3MATLAB程序如下所示：a=[ 。 ]。a1=[ 。 ]。q=[1 0。0 1]。b=[。 ]。e=[ 。0 ]。f1=[ 。 ]。f2=[ 。 ]。 I=[1 0。 0 1]。setlmis([])。x=lmivar(1,[2 1])。y=lmivar(2,[1 2])。ew=lmivar(2,[1 1])。lmiterm([1 1 1 x],1,a39。,39。s39。)。 % LMI 1: x*a39。+a*xlmiterm([1 1 1 y],b,1,39。s39。)。 % LMI 1: b*y+y39。*b39。lmiterm([1 2 1 0],a139。)。 % LMI 1: a139。lmiterm([1 2 2 0],q)。 % LMI 1: qlmiterm([1 3 1 x],f1,1)。 % LMI 1: f1*xlmiterm([1 3 3 0],inv(ew)*I)。 % LMI 1: inv(ew)*Ilmiterm([1 4 1 x],1,1)。 % LMI 1: xlmiterm([1 4 4 0],inv(q))。 % LMI 1: inv(q)lmiterm([1 5 1 0],e39。)。 % LMI 1: e39。lmiterm([1 5 5 ew],.5*,I,39。s39。)。 % LMI 1: *ew*I (NON SYMMETRIC?) lmiterm([1 6 2 0],f2)。 % LMI 1: f2lmiterm([1 6 6 0],inv(ew)*I)。 % LMI 1: inv(ew)*Ilmiterm([2 1 1 x],1,1)。 % LMI 2: xlmiterm([3 1 1 ew],1,1)。 % LMI 3: ewss2=getlmis。 [lp kp]=feasp(ss2)。x=dec2mat(ss2,kp,x)y=dec2mat(ss2,kp,y)ew=dec2mat(ss2,kp,ew)p=inv(x)k=y*inv(x)附錄 4MATLAB程序如下所示：a=[ 。 ]。a1=[ 。 ]。q=[1 0。0 1]。b=[。 ]。e=[ 。0 ]。f1=[ 。 ]。f2=[ 。 ]。 I=[1 0。 0 1]。setlmis([])。x=lmivar(1,[2 1])。y=lmivar(2,[1 2])。ew1=lmivar(2,[1 1])。ew2=lmivar(2,[1 1])。ew3=lmivar(2,[1 1])。ew4=lmivar(2,[1 1])。ew5=lmivar(2,[1 1])。ew6=lmivar(2,[1 1])。p1=lmivar(1,[2 1])。d1=lmivar(2,[1 1])。 lmiterm([1 1 1 x],1,a39。,39。s39。)。 % LMI 1: x*a39。+a*xlmiterm([1 1 1 y],b,1,39。s39。)。 % LMI 1: b*y+y39。*b39。lmiterm([1 1 1 x],1,a139。,39。s39。)。 % LMI 1: x*a139。+a1*xlmiterm([1 2 2 ew4],.5*1,f239。*f2,39。s39。)。 % LMI 1: ew4*f239。*f2 (NON SYMMETRIC?)lmiterm([1 2 2 0],q)。 % LMI 1: qlmiterm([1 3 1 x],a,1)。 % LMI 1: a*xlmiterm([1 3 1 y],b,1)。 % LMI 1: b*ylmiterm([1 3 2 0],a1)。 % LMI 1: a1lmiterm([1 3 3 d1],.5*1,p1,39。s39。)。 % LMI 1: d1*p1 (NON SYMMETRIC?)lmiterm([1 3 3 ew5],.5*1,e39。*e,39。s39。)。 % LMI 1: ew5*e39。*e (NON SYMMETRIC?)lmiterm([1 4 1 0],a139。)。 % LMI 1: a139。lmiterm([1 4 4 d1],.5*1,inv(p1),39。s39。)。 % LMI 1: d1*inv(p1) (NON SYMMETRIC?)lmiterm([1 4 4 ew6],.5*1,f239。*f2,39。s39。)。 % LMI 1: ew6*f239。*f2 (NON SYMMETRIC?)lmiterm([1 5 1 0],e39。)。 % LMI 1: e39。lmiterm([1 5 5 0],inv(ew1)*I)。 % LMI 1: inv(ew1)*Ilmiterm([1 6 1 x],1,1)。 % LMI 1: xlmiterm([1 6 6 0],inv(q))。 % LMI 1: inv(q)lmiterm([1 7 1 x],f1,1)。 % LMI 1: f1*xlmiterm([1 7 7 0],inv(ew2)*I)。 % LMI 1: inv(ew2)*Ilmiterm([1 8 1 x],f2,1)。 % LMI 1: f2*xlmiterm([1 8 8 0],inv(ew3)*I)。 % LMI 1: inv(ew3)*Ilmiterm([2 1 1 x],1,1)。 % LMI 2: xlmiterm([3 1 1 ew1],1,1)。 % LMI 3: ew1lmiterm([4 1 1 ew2],1,1)。 % LMI 4: ew2lmiterm([5 1 1 ew3],1,1)。 % LMI 5: ew3lmiterm([6 1 1 ew4],1,1)。 % LMI 6: ew4lmiterm([7 1 1 ew5],1,1)。 % LMI 7: ew5lmiterm([8 1 1 ew6],1,1)。 % LMI 8: ew6ss3=getlmis。 [lp kp]=feasp(ss3)。x=dec2mat(ss3,kp,x)y=dec2mat(ss3,kp,y)ew1=dec2mat(ss3,kp,ew1)ew2=dec2mat(ss3,kp,ew2)ew3=dec2mat(ss3,kp,ew3)ew4=dec2mat(ss3,kp,ew4)ew5=dec2mat(ss3,kp,ew5)ew6=dec2mat(ss3,kp,ew6)d1=dec2mat(ss3,kp,d1)d=inv(d1)p=inv(x)k=y*inv(x)