掌握线性矩阵不等式的求解方法及Matlab工具箱

1、精选优质文档-倾情为你奉上年级班号学号专业姓名实验名称数学软件实验课程考察报告线性矩阵不等式问题的求解方法实验类型设计型综合型创新型实验目的或要求1. 掌握线性矩阵不等式的求解方法及Matlab工具箱。2. 求解两个矩阵不等式问题。3. 验证结果的正确性。实验程序的源代码1.A=0 1 0 0;1 0 0 0;0 0 0 1;0 0 1 0;B=0 0;1 0;0 0;0 1;H=0 0 0 0;-1 0 1 0;0 0 0 0;1 0 -1 0;I=eye(4);setlmis();Y=lmivar(1,2 0;2 0);L=lmivar(2,2 4);gamma=lmivar(1,1 0)

2、;beta=lmivar(1,1 0);lmiterm(1 1 1 Y,1,A','s');lmiterm(1 1 1 L,B,1,'S');lmiterm(1 2 1 0,I);lmiterm(1 3 1 Y,H,1);lmiterm(-2 1 1 Y,1,1);lmiterm(1 2 2 gamma,5.*1,-I,'s');lmiterm(1 3 3 beta,5.*1,-I,'s');lmisys=getlmis;tmin,xfeas=feasp(lmisys);Y_ans=dec2mat(lmisys,xfeas

3、,Y)L_ans=dec2mat(lmisys,xfeas,L)gamma_ans=dec2mat(lmisys,xfeas,gamma)beta_ans=dec2mat(lmisys,xfeas,beta)q=Y_ans*A'+A*Y_ans+L_ans*B'+B*L_ans;z=zeros(4,4);r=q I Y_ans.*H'I gamma_ans.*I. z;H.*Y_ans. z beta_ans.*I;v=eig(r);display(v);2.A=0 1;1 0;B=0;1;I=eye(2);setlmis();Y=lmivar(1,1 0;1 0);L

4、=lmivar(2,1 2);gamma=lmivar(1,1 0);beta=lmivar(1,1 0);lmiterm(1 1 1 Y,1,A','s');lmiterm(1 1 1 L,B,1,'S');lmiterm(1 2 1 0,I);lmiterm(1 3 1 Y,1,1);lmiterm(-2 1 1 Y,1,1);lmiterm(1 2 2 0,-gamma.*1);lmiterm(1 3 3 beta,5.*1,-I,'s');lmisys=getlmis;tmin,xfeas=feasp(lmisys);Y_ans=dec2mat(lmisys,xfeas,Y)L_ans=dec2mat(lmisys,xfeas,L)gamma_ans=dec2mat(lmisys,xfeas,gamma)beta_ans=dec2mat(lmisys,xfeas,beta)q=Y_ans*A'+A*Y_ans+L_ans*B'+B*L_ans;z=zeros(2,2);r=q I Y_ans;I -gamma_ans*I z;Y_ans z -beta_ans*I;