| 查看: 3168 | 回复: 10 | ||
wangxingye新虫 (小有名气)
|
[求助]
matlab计算线性矩阵不等式（LMI）时出现维数与同一列不一致的问题 已有2人参与
|
|
遇到一个求解线性矩阵不等式的问题，利用MATLAB进行编程求解的过程中遇到“Error using lmiterm (line 293) lhs of LMI #2, block (4,2): term dimensions incompatible withother terms in same row”，但是我已经检查了许多次了，不论是pos(2,4)还是pos(4,2),不管是所在行还是所在列，矩阵的维数都和同一行和列的其他矩阵一致，不应该有这样的问题啊，实在是没办法了，求教大神帮忙看一下，我把不等式和程序都贴上来。 程序如下： clear all; clc; ms=320; mt=40; ks=18000; kt=200000; cs=1000; ct=10; g = 9.8; %控制器参数 rou = 1.6; ite = 10000; gama = 1; w1 = 8*pi;%4Hz w2 = 16*pi;%8Hz wc =(w1+w2)/2; xstmax = 0.1;%最大动行程100mm umax = 2500;%最大主动出力2000N A=[0 0 1 -1; 0 0 0 1; -ks/ms 0 -cs/ms cs/ms; ks/mt -kt/mt cs/mt -(cs+ct)/mt];%系统矩阵 B=[0; 0; 1/ms; -1/mt];%控制矩阵 B1=[0; -1; 0; ct/mt]; C1=[-ks/ms 0 -cs/ms cs/ms];%输出矩阵 C2=[1/xstmax 0 0 0; 0 kt/(ms+mt)/g 0 0]; D1=1/ms;%传递矩阵 %定义LMI setlmis([]); %首先定义决策变量 P = lmivar(1, [4 1]); P1= lmivar(1, [4 1]); Q = lmivar(1, [4 1]); F = lmivar(2, [4 4]); %4*4矩阵 K1= lmivar(2, [1 4]); %1*4矩阵 %2nd pos(1 1) lmiterm([2 1 1 Q], -1, 1); %2nd pos(1 2) lmiterm([2 1 2 P], 1, 1); lmiterm([2 1 2 F], -1, 1); %2nd pos(1 6) lmiterm([2 1 6 Q], wc, 1); %2nd pos(2 2) lmiterm([2 2 2 Q], w1*w2, -1); lmiterm([2 2 2 F], A, 1, 's'); lmiterm([2 2 2 K1], B, 1, 's'); %2nd pos(2 3) lmiterm([2 2 3 0], B1); %2nd pos(2 4) lmiterm([2 2 4 F'], 1, C1'); lmiterm([2 2 4 K1'], 1, D1'); %2nd pos(2 5) lmiterm([2 2 5 Q], wc, -1); %2nd pos(3 3) lmiterm([2 3 3 0],-1*gama^2); %2nd pos(4 4) lmiterm([2 4 4 0], -1); %2nd pos(5 5) lmiterm([2 5 5 Q], -1, 1); %2nd pos(5 6) lmiterm([2 5 6 P], 1, 1); lmiterm([2 5 6 F], -1, 1); %2nd pos(6 6)计算线性矩阵不等式（LMI）时 lmiterm([2 6 6 Q], w1*w2, -1); lmiterm([2 6 6 F], A, 1, 's'); lmiterm([2 6 6 K1], B, 1, 's'); %2nd pos(6 7) lmiterm([2 6 7 0], B1); %2nd pos(6 8) lmiterm([2 6 8 F'], 1, C1'); lmiterm([2 6 8 K1'], 1, D1'); %2nd pos(7 7) lmiterm([2 7 7 0],-1*gama^2); %2nd pos(8 8) lmiterm([2 8 8 0], -1);  1.JPG |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有165人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有4人回复
wangxingye
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1510.5
- 散金: 6
- 红花: 1
- 帖子: 97
- 在线: 49.7小时
- 虫号: 3908308
- 注册: 2015-06-04
- 性别: GG
- 专业: 机械工程
2楼2017-09-06 20:49:27
wangxingye
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1510.5
- 散金: 6
- 红花: 1
- 帖子: 97
- 在线: 49.7小时
- 虫号: 3908308
- 注册: 2015-06-04
- 性别: GG
- 专业: 机械工程
3楼2017-09-08 07:36:25
wangxingye
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1510.5
- 散金: 6
- 红花: 1
- 帖子: 97
- 在线: 49.7小时
- 虫号: 3908308
- 注册: 2015-06-04
- 性别: GG
- 专业: 机械工程
4楼2017-09-12 19:38:45
fxwangwei87
金虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 447.4
- 散金: 1086
- 红花: 38
- 帖子: 1117
- 在线: 497.9小时
- 虫号: 1100466
- 注册: 2010-09-16
- 性别: GG
- 专业: 工程热物理相关交叉领域
5楼2017-09-13 15:18:03
wangxingye
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1510.5
- 散金: 6
- 红花: 1
- 帖子: 97
- 在线: 49.7小时
- 虫号: 3908308
- 注册: 2015-06-04
- 性别: GG
- 专业: 机械工程
6楼2017-09-13 15:21:54
wangxingye
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1510.5
- 散金: 6
- 红花: 1
- 帖子: 97
- 在线: 49.7小时
- 虫号: 3908308
- 注册: 2015-06-04
- 性别: GG
- 专业: 机械工程
7楼2017-09-16 12:50:39
wangxingye
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1510.5
- 散金: 6
- 红花: 1
- 帖子: 97
- 在线: 49.7小时
- 虫号: 3908308
- 注册: 2015-06-04
- 性别: GG
- 专业: 机械工程
8楼2017-09-20 12:52:39
9楼2019-08-13 09:42:02
10楼2020-04-12 20:12:51
5