| 查看: 1798 | 回复: 3 | ||
su2gges2tion木虫 (小有名气)
|
[求助]
使用matlab中LMI工具箱求解线性矩阵不等式程序
|
|
求解矩阵不等式A1'P+PA1+F'F+PEE'P-K1'B'P-PBK1<0;A1'P+PA1+F'F+PEE'P-K2'B'P-PBK2<0;A2'P+PA2+F'F+PEE'P-K1'B'P-PBK1<0;A2'P+PA2+F'F+PEE'P-K2'B'P-PBK2<0;P>0.其中B为三阶单位矩阵，K1，K2，P为所求的未知矩阵，为什么我写的程序求出的tmin是大于零的，求出的P矩阵不是正定的，并且求得K1，K2均是零矩阵。我的程序哪里出错了，求指教，非常感谢！ %这个问题是二次稳定性问题中提取出来的，首先确定已知的LMI，在调用函数feasp %常数 A1=[-1 -20 0;20 -1 20;0 5.45 -5.45]; A2=[-1 20 0;-20 -1 20;0 5.45 -5.45]; E=[0 0 0 0 0 0 0 0 0;0 0 0 0 0 1 0 0 0;0 0 0 0 0 0 0 sqrt(0.2) sqrt(0.2)]; F=[0 0 0;0 0 0;0 0 0;0 0 0;0 0 0;0 0 sqrt(1);0 0 0;0 sqrt(0.2) 0;0 0 sqrt(0.2)]; %初始化LMI setlmis([]); %定义变量 P=lmivar(1,[3 1]); K1=lmivar(2,[3 3]); K2=lmivar(2,[3 3]); %添加项 lmiterm([1 1 1 P],1,A1,'s'); lmiterm([1 1 1 0],F'*F); lmiterm([1 1 1 P],-1,K1); lmiterm([1 1 1 P],-K1,-1); lmiterm([1 1 2 P],1,E); lmiterm([1 2 2 0],-1); lmiterm([2 1 1 P],1,A1,'s'); lmiterm([2 1 1 0],F'*F); lmiterm([2 1 1 P],-1,K2); lmiterm([2 1 1 P],-K2,-1); lmiterm([2 1 2 P],1,E); lmiterm([2 2 2 0],-1); lmiterm([3 1 1 P],1,A2,'s'); lmiterm([3 1 1 0],F'*F); lmiterm([3 1 1 P],-1,K1); lmiterm([3 1 1 P],-K1,-1); lmiterm([3 1 2 P],1,E); lmiterm([3 2 2 0],-1); lmiterm([4 1 1 P],1,A2,'s'); lmiterm([4 1 1 0],F'*F); lmiterm([4 1 1 P],-1,K2); lmiterm([4 1 1 P],-K2,-1); lmiterm([4 1 2 P],1,E); lmiterm([4 2 2 0],-1); lmiterm([-5 1 1 P],1,1); %获取LMI系统描述 lmisys=getlmis; [tmin,xfeas]=feasp(lmisys)%options参数可以自己设置，也可以默认 |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 矩阵及其子矩阵的奇异值之间的不等式或者等式关系
已经有9人回复
- 关于matlab的LMI求助
已经有8人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- LMI中如何定义一个非零的矩阵变量
已经有8人回复
- .dat文件在matlab中直接应用吗，以矩阵形式
已经有12人回复
- matlab中如何随机不重复选取矩阵中的元素
已经有4人回复
- matlab线性矩阵不等式求解做优化问题
已经有10人回复
su2gges2tion
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2660.3
- 散金: 46
- 红花: 1
- 帖子: 228
- 在线: 100.8小时
- 虫号: 2574753
- 注册: 2013-07-31
- 性别: GG
- 专业: 高电压与绝缘
2楼2013-08-21 10:13:48
3楼2020-04-22 07:20:17
|
本帖内容被屏蔽 |
4楼2022-06-07 21:24:54
