应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » matlab计算线性矩阵不等式（LMI）时出现维数与同一列不一致的问题
汕头大学海洋科学接受调剂
111/2返回列表12下一页
查看: 3170  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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);

matlab计算线性矩阵不等式（LMI）时出现维数与同一列不一致的问题
1.JPG
回复此楼

» 猜你喜欢
1楼 2017-09-06 20:38:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
刚开始接触LMI不久，还请大家多多指教，我自己先顶一顶！
赞一下 回复此楼
2楼2017-09-06 20:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
内容已删除
回复此楼
3楼2017-09-08 07:36:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
求大神指导，帖子都要沉了

发自小木虫Android客户端
赞一下 回复此楼
4楼2017-09-12 19:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxwangwei87

金虫 (著名写手)
帮助顶一下，我现在需要求解非线性不等式组，能给我看看怎么解答吗？
赞一下 回复此楼
做一只快乐的小金虫，专心科研，面朝大海，春暖花开~~~~
5楼2017-09-13 15:18:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
我也是才接触不久，现在连最简单的不等式都解不出来，怕你帮不上你忙啊

发自小木虫Android客户端
赞一下 回复此楼
6楼2017-09-13 15:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
再顶一次，真是没想到好的解决方案，期待有虫友来帮助解惑

发自小木虫Android客户端
赞一下 回复此楼
7楼2017-09-16 12:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxingye

新虫 (小有名气)
最后顶一次，抱着最后一点希望，恳请有经验的同仁能够给予指点！

发自小木虫Android客户端
赞一下 回复此楼
8楼2017-09-20 12:52:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小郝a学习

银虫 (小有名气)
帮顶，目前还没接触到这
赞一下 回复此楼
9楼2019-08-13 09:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Verney

新虫 (初入文坛)

【答案】应助回帖

三年前的帖子了...
原代码中有关于矩阵变量转置的语法错误
比如
%2nd pos(2 4)
lmiterm([2 2 4 F'], 1, C1');
lmiterm([2 2 4 K1'], 1, D1');
这里如果想做矩阵变量K1的转置应该是lmiterm([2 2 4 -K1], 1, D1');
具体参见lmiterm help
赞一下 回复此楼
10楼2020-04-12 20:12:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangxingye 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 材料工程281还有调剂机会吗 +43 xaw. 2026-04-11 44/2200 2026-04-15 12:46 by 西北望—风沙
[考研] 289 分105500药学专硕求调剂(找B区学校) +3 白云123456789 2026-04-13 3/150 2026-04-15 12:42 by 西北望—风沙
[考研] 307中医考研调剂 +3 于以采蘩 2026-04-14 3/150 2026-04-15 11:16 by xtt..
[考研] 材料相关专业344求调剂双非工科学校或课题组 +21 hualkop 2026-04-12 22/1100 2026-04-14 17:44 by lhj2009
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +8 慕绝cc 2026-04-09 8/400 2026-04-13 14:08 by 张zhihao
[考研] 2本，初试303，0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 一志愿东北大学控制工程085406数二英二385，求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 一志愿矿大，材料工程专硕314分，0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
信息提示
关闭
请填处理意见
关闭