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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » matlab计算线性矩阵不等式（LMI）时出现维数与同一列不一致的问题
111/2返回列表12下一页
查看: 2869  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭