| 查看: 3058 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
小_蚂蚁新虫 (小有名气)
|
[求助]
LMI
|
||
|
这样的LMI如何求出系统延时的上界？？跪求大神指导。我自己编的程序是这样的下面那个样子，但是是不对的，我想听听大神的方法。非常感谢各位大神 clc; clear; A=[-18/7 9 0; 1 -1 1; 0 -14.28 0];% n = length(A); B=[27/7 0 0]';% C=[1 0 0]; H=[1 0 0]; k=1; setlmis([]); % Define Variables %P Q is a symmetric matrix, has a block size of 2 and this block is symmetric Q = lmivar(1,[3 1]); P = lmivar(1,[3 1]); Y = lmivar(2,[3 1]); M1=lmivar(1,[3 1]); M2=lmivar(1,[3 1]); V =lmivar(1,[1 0]); % tau_max=lmivar(1,[1 0]); % Q, P, V > 0 lmiterm([-1 1 1 Q], 1, 1); lmiterm([-2 1 1 P], 1, 1); lmiterm([-3 1 1 V], 1, 1); %lmiterm([5 1 1 tau_max], 1, 1); % pos in (1, 1) lmiterm([-4 1 1 P], -1, A, 's'); lmiterm([-4 1 1 Q], -1, 1); lmiterm([-4 1 1 M1], -1, 1); lmiterm([-4 1 1 -M1], -1, 1); % pos (1, 2) lmiterm([-4 1 2 Y], 1, H); lmiterm([-4 1 2 -M1], 1, 1); lmiterm([-4 1 2 M2], -1, 1); % pos(1, 3) lmiterm([-4 1 3 P], -1, B); lmiterm([-4 1 3 V], -k*C',1); % pos(1, 4) lmiterm([4 1 4 P], -A', 1); % pos(1, 5) lmiterm([4 1 5 -M1], -1, 1); % pos(2, 2) lmiterm([-4 2 2 Q], 1, 1); lmiterm([-4 2 2 M2], 1, 1); lmiterm([-4 2 2 -M2], 1, 1); % pos(2, 4) lmiterm([4 2 4 -Y], -H', -1); % pos(2, 5) lmiterm([4 2 5 -M2], 1, -1); % pos(3, 3) lmiterm([-4 3 3 V], 2, 1); % pos(3, 4) lmiterm([4 3 4 P], -B', 1); % pos(4, 4) lmiterm([-4 4 4 P], 1, 1); % pos(5, 5) lmiterm([-4 5 5 P], 1, 1); lmis = getlmis; [t,xopt]=gevp(lmis,1); t tt=1/t QQ = dec2mat(lmis,xopt,Q) PP = dec2mat(lmis,xopt,P) YY = dec2mat(lmis,xopt,Y) vv=dec2mat(lmis,xopt,V) K=inv(PP)*YY    发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
小_蚂蚁
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 601.7
- 帖子: 113
- 在线: 17.6小时
- 虫号: 5424975
- 注册: 2016-12-30
- 性别: MM
- 专业: 导航、制导与传感技术
|
你上面那个例子，我之前看到过，也套用了，但是有一个问题，你的例子的矩阵，它的对角线上都有a，但我那个LMI对角线上的不全是tau，我下面做了两种变换，但都求不出来。我真不知道要怎么办，求您指导啊。第一张图，是根据定理里面的lmi变换出来的，第二张图，是根据证明里面的不等式变换出来的，这里出现了一个问题 ：M1，P，M2都是变量，矩阵中M1P^(-1)M2这一项怎么在LMI-tool表示，还有H'Y'P^(-1)PY这里面Y和P也都是变量，这要怎么办。   发自小木虫Android客户端 |
6楼2017-07-31 10:10:45
youandiandhe
铁杆木虫 (知名作家)
- 应助: 25 (小学生)
- 金币: 13505.9
- 散金: 5605
- 红花: 29
- 帖子: 7999
- 在线: 2667.4小时
- 虫号: 3483828
- 注册: 2014-10-18
- 专业: 控制理论与方法
2楼2017-07-29 22:28:14
youandiandhe
铁杆木虫 (知名作家)
- 应助: 25 (小学生)
- 金币: 13505.9
- 散金: 5605
- 红花: 29
- 帖子: 7999
- 在线: 2667.4小时
- 虫号: 3483828
- 注册: 2014-10-18
- 专业: 控制理论与方法
3楼2017-07-29 22:28:35
youandiandhe
铁杆木虫 (知名作家)
- 应助: 25 (小学生)
- 金币: 13505.9
- 散金: 5605
- 红花: 29
- 帖子: 7999
- 在线: 2667.4小时
- 虫号: 3483828
- 注册: 2014-10-18
- 专业: 控制理论与方法
4楼2017-07-29 22:28:52
10