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

PAGE PAGE 5 年级 班号 学号 专业 姓名 实验 名称 数学软件实验课程考察报告 线性矩阵不等式问题的求解方法 实验类型 设计型 综合型 创新型 √ 实验目的或要求 掌握线性矩阵不等式的求解方法及Matlab工具箱。 求解两个矩阵不等式问题。 验证结果的正确性。 实 验 程 序 的 源 代 码 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]); 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,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=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]; v=eig(r); display(v); 实 验 结 果 及 分 析 一： 二： 心 得 体 会 成 绩 评 定 教师签名： 年 月 日 已知： 求 满足 问题2 已知： 求 满足