|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 giftdreamer 于 2011-7-4 23:13 编辑
LMI形式如(44)、(45)所示 ,论文中的结果如另一幅图所示,但我计算的结果和他们的不一致
程序如下 不知道什么原因?
望大家指点
- A1 = [0 1;17.294 0];
- B1 = [0;-0.1765];
- A2 = [0 1;9.365 0];
- B2 = [0;-0.0349];
- % LMI Description
- setlmis([]); % 定义一个LMI系统
- X=lmivar(1,[2 1]); % inv(P) 定义LMI中的矩阵变量
- N1=lmivar(2,[1 2]);
- N2=lmivar(2,[1 2]);
- % 1 LMI
- lmiterm([-1 1 1 X],1,1); % X > 0
- % 2 LMI:
- lmiterm([2 1 1 X],A1,1,'s');
- lmiterm([2 1 1 N1],-B1,1,'s');
- lmiterm([2 1 1 X],10,1);
- % 3 LMI
- lmiterm([3 1 1 X],A2,1,'s');
- lmiterm([3 1 1 N2],-B2,1,'s');
- lmiterm([3 1 1 X],10,1);
- % 4 LMI
- lmiterm([4 1 1 X],A1,1,'s');
- lmiterm([4 1 1 X],A2,1,'s');
- lmiterm([4 1 1 N2],-B1,1,'s');
- lmiterm([4 1 1 N1],-B2,1,'s');
- lmiterm([4 1 1 X],20,1);
- LMIs=getlmis;
- [alpha,popt]=feasp(LMIs);
- X=dec2mat(LMIs,popt,1);
- N1=dec2mat(LMIs,popt,2);
- N2=dec2mat(LMIs,popt,3);
-
- P=inv(X);
- K1=N1*P;
- K2=N2*P;
-
- % show results
- P=inv(X)
- K1
- K2
复制代码
计算结果
- P =
- 345.1201 36.3569
- 36.3569 5.5142
- K1 =
- 1.0e+003 *
- -1.3727 -0.1959
- K2 =
- 1.0e+003 *
- -4.0205 -0.5748
复制代码
|
|