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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 关于matlab的LMI求助
91/1返回列表
查看: 2592  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fxlily2008

木虫 (小有名气)

[求助] 关于matlab的LMI求助

各位学长学姐，我从来没有用过matlab的LMI工具箱编写过程序，现在我想用它的mincx实现一个极简单的功能:
target:  min(0.1*x1+0.2*x2+0.3*x3+0.4*x4)
   s.t. :  x1,x2,x3,x4>=0
            x1+x2+x3+x4>10

我试着编了一个小程序，但是结果好像不对，请高手帮助一下吧，不胜感激，小弟在线等。

最好回复一个完整的程序吧，谢谢！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-05-15 19:35:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-05-26 12:36:34
你这个条件就有问题，没有解！
赞一下 回复此楼
2楼2013-05-15 20:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxlily2008

木虫 (小有名气)
引用回帖:
2楼: Originally posted by lijie169 at 2013-05-15 20:07:30
你这个条件就有问题，没有解！

不是啊，虽然x1,x2,x3,x4可以无限大，但min(……)本身就是约束啊，理论上正确的结果应该是x1=10.0000001，x2=x3=x4=0对吧。但是用matlab的LMI怎么写程序啊，求帮助
赞一下 回复此楼
3楼2013-05-15 21:09:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxlily2008

木虫 (小有名气)
另外，其中的等于号可以不考虑，没关系的
赞一下 回复此楼
4楼2013-05-15 21:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

idnk

金虫 (小有名气)
★ ★
感谢参与，应助指数 +1
1592203609: 金币+2, 应助指数-1, 非应助帖不要选择应助 2013-05-15 23:00:23
把你程序贴出来，大家一起看看什么问题
赞一下 回复此楼
5楼2013-05-15 22:36:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elvinsu

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fxlily2008: 金币+5, ★★★很有帮助 2013-05-18 14:45:27
xiegangmai: 金币+1, 谢谢参与！ 2013-05-19 08:43:35
LMI是可以处理， 但是用线性规划不是更直接吗？ x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)
赞一下 回复此楼
科学的历程就是在假设求证中不断进步--巴瑞•马歇尔教授
6楼2013-05-18 11:56:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxlily2008

木虫 (小有名气)
引用回帖:
6楼: Originally posted by elvinsu at 2013-05-18 11:56:08
LMI是可以处理， 但是用线性规划不是更直接吗？ x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)

哦，这样啊，这个函数我不熟，我先试试，谢谢前辈指导，我先给你5分吧
赞一下 回复此楼
7楼2013-05-18 14:44:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elvinsu

金虫 (小有名气)
引用回帖:
7楼: Originally posted by fxlily2008 at 2013-05-18 14:44:53
哦，这样啊，这个函数我不熟，我先试试，谢谢前辈指导，我先给你5分吧...

不敢当，谢谢加分！
回复此楼
科学的历程就是在假设求证中不断进步--巴瑞•马歇尔教授
8楼2013-05-19 21:00:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

EVRYIBISC

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
fxlily2008: 金币+5, ★★★很有帮助 2013-06-07 09:04:54
xiegangmai: 金币+2, 谢谢参与！ 2013-06-07 20:36:25
clear all;
close all;
clc

setlmis([]);

x1 = lmivar(1,[1,1]);
x2 = lmivar(1,[1,1]);
x3 = lmivar(1,[1,1]);
x4 = lmivar(1,[1,1]);

br1 = newlmi;

lmiterm([-br1 1 1 x1],1,1);

br2 = newlmi;

lmiterm([-br2 1 1 x2],1,1);

br3 = newlmi;

lmiterm([-br3 1 1 x3],1,1);

br4 = newlmi;

lmiterm([-br4 1 1 x4],1,1);

br5 = newlmi;

lmiterm([-br5 1 1 x1],1,1);
lmiterm([-br5 1 1 x2],1,1);
lmiterm([-br5 1 1 x3],1,1);
lmiterm([-br5 1 1 x4],1,1);
lmiterm([br5 1 1 0],10);

lmisys = getlmis;

n = decnbr(lmisys);
c = zeros(n,1);
c(n,1)=1;

for jj = 1:n
    [x1j,x2j,x3j,x4j] = defcx(lmisys,jj,x1,x2,x3,x4);
    c(jj) = 0.1*x1j+0.2*x2j+0.3*x3j+0.4*x4j
end
%c = mat2dec(lmisys,zeros(2,2),zeros(1,2),1);

options = [1e-5,0,0,0,0];

[copt,xopt] = mincx(lmisys,c,options);

x1opt = dec2mat(lmisys,xopt,x1)
x2opt = dec2mat(lmisys,xopt,x2)
x3opt = dec2mat(lmisys,xopt,x3)
x4opt = dec2mat(lmisys,xopt,x4)

运行结果:
c =

    0.1000
         0
         0
    1.0000


c =

    0.1000
    0.2000
         0
    1.0000


c =

    0.1000
    0.2000
    0.3000
    1.0000


c =

    0.1000
    0.2000
    0.3000
    0.4000


Solver for linear objective minimization under LMI constraints

Iterations   :    Best objective value so far

     1                   1.734323
     2                   1.133169
***                 new lower bound:     0.563924
     3                   1.022152
***                 new lower bound:     0.834201
     4                   1.013815
***                 new lower bound:     0.977789
     5                   1.000993
***                 new lower bound:     0.997951
     6                   1.000010
***                 new lower bound:     0.998732
     7                   1.000010
***                 new lower bound:     0.999878
     8                   1.000001
***                 new lower bound:     0.999973
     9                   1.000001
***                 new lower bound:     0.999994

Result:  feasible solution of required accuracy
          best objective value:     1.000001
          guaranteed relative accuracy:  6.36e-06
          f-radius saturation:  0.000% of R =  1.00e+09


x1opt =

   10.0000


x2opt =

   1.2620e-06


x3opt =

   6.3082e-07


x4opt =

   4.2051e-07

>>
赞一下(3人) 回复此楼
9楼2013-06-06 21:35:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fxlily2008 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +4 shxchem 2026-03-20 6/300 2026-03-21 11:54 by 无际的草原
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿武理材料工程348求调剂 +3 ￣^￣゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by ！本暗一次！
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
信息提示
关闭
请填处理意见
关闭