| 查看: 849 | 回复: 0 | ||
[求助]
线性矩阵不等式,feasp求解，原文算例无法重现
|
|
求助~有大神可以帮忙看一下程序吗？ 鲁棒控制+时滞问题，想复现一下文章的例子2. 这个lmi程序在matlab里可以求解，tmin<0，有值，但是得到的结果不对，得到的反馈律K接近0，跟文中的算例的结果不一样，不知道是哪里错了，可以有人帮忙解决一下吗？多谢了~ %% lmi 原文算例 2009Theorem4 clear E=[1 0; 0 0]; A=[0 0.46;-1 -2]; B1=[1 1.5]'; D=[0.2 0; 0 0.2]; D1=D; H=eye(2); H1=[1 1]'; [~,na]=size(A); taomax0=1.5; %% LMI定义 setlmis([]) %都是任意阶数为A的矩阵 Q1=lmivar(2,[na,na]);%是na*na的非对称阵 Z1=lmivar(2,[na,na]); % Q1=lmivar(1,[na,1]);%是na*na的对称阵 % Z1=lmivar(1,[na,1]); N=lmivar(2,[1,na]); %N的阶数（x,y），x是B的列数，是u的阶数,y已更正为na阶 X=lmivar(2,[na,na]); %X非奇异 X1=lmivar(2,[na,na]); Y1=lmivar(2,[na,na]); eps0=lmivar(1,[1,0]); %定义标量的时候，后面那个应该是写0 但是写0写1都不影响 eps1=lmivar(1,[1,0]); %% 不等式描述 % lmiterm([-1 1 1 Q1],1,1); lmiterm([-2 1 1 Z1],1,1); % 2nd LMI %sum11 lmiterm([3 1 1 X],A,1,'s'); lmiterm([3 1 1 X1],taomax0,1); lmiterm([3 1 1 Y1],1,1,'s'); lmiterm([3 1 1 Q1],1,1); lmiterm([3 1 1 eps0],1,D*D'); lmiterm([3 1 1 eps1],1,D1*D1'); %sum12 lmiterm([3 1 2 N],B1,1); lmiterm([3 1 2 Y1],-1,1); %sum22 lmiterm([3 2 2 Q1],-1,1); %(theta12) lmiterm([3 1 3 -X],taomax0,A'); lmiterm([3 1 3 eps0],1,taomax0*(D*D')); lmiterm([3 1 3 eps1],1,taomax0*(D1*D1')); lmiterm([3 2 3 -N],taomax0,B1'); %(theta22) lmiterm([3 3 3 Z1],-taomax0,1); lmiterm([3 3 3 eps0],1,taomax0*taomax0*(D*D')); lmiterm([3 3 3 eps1],1,taomax0*taomax0*(D1*D1')); %(theta13) lmiterm([3 1 4 -X],1,H'); lmiterm([3 2 5 -N],1,H1'); %(theta33) lmiterm([3 4 4 eps0],-eye(2),1); lmiterm([3 5 5 eps1],-eye(2),1); % lmiterm([3 4 4 eps0],-1,1); % lmiterm([3 5 5 eps1],-1,1); lmis=getlmis; options(5)=1; [tmin,feas]=feasp(lmis,options) %% 求解控制参数 N=dec2mat(lmis,feas,N) X=dec2mat(lmis,feas,X) Y1=dec2mat(lmis,feas,Y1) X1=dec2mat(lmis,feas,X1) Q1=dec2mat(lmis,feas,Q1) Z1=dec2mat(lmis,feas,Z1) eps0=dec2mat(lmis,feas,eps0) eps1=dec2mat(lmis,feas,eps1) K=N/X |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 2009.pdf
2015-11-17 09:58:47, 444.21 K
