24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3167  |  回复: 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);

matlab计算线性矩阵不等式(LMI)时出现维数与同一列不一致的问题
1.JPG
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 312求调剂 +17 门路摸摸 2026-04-14 18/900 2026-04-15 13:11 by 西北望—风沙
[考研] 复试调剂 +21 积极向上; 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 297,工科调剂? +6 河南农业大学-能 2026-04-14 6/300 2026-04-14 16:04 by Art1977
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 314求调剂 +24 wakeluofu 2026-04-09 25/1250 2026-04-13 08:58 by lhj2009
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 一志愿东北大学控制工程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
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
信息提示
请填处理意见