| 查看: 2595 | 回复: 8 | |||
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 我试着编了一个小程序，但是结果好像不对，请高手帮助一下吧，不胜感激，小弟在线等。 最好回复一个完整的程序吧，谢谢！！ |
回复此楼» 猜你喜欢
298求调剂
已经有3人回复
-
336求调剂
已经有4人回复
-
296求调剂
已经有3人回复
-
266求调剂
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有10人回复
-
初试 317
已经有4人回复
-
化学调剂
已经有3人回复
-
286求调剂
已经有7人回复
-
0703化学调剂
已经有3人回复
-
一志愿重庆大学085700资源与环境总分308求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助一个LMI问题
已经有10人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- 求助关于Matlab中字符串的转化问题
已经有4人回复
- 纯新手，全方面求助关于MATLAB........
已经有9人回复
- matlab画图求助
已经有11人回复
- 继续求助Matlab绘图！
已经有8人回复
- 求助matlab编程
已经有13人回复
- 求助svm在matlab中的应用
已经有7人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 求助！！关于MATLAB
已经有6人回复
- 【求助】会用Matlab的LMI toolbox的请帮帮忙~~
已经有5人回复
- 【求助】时间同步的MATLAB仿真
已经有5人回复
- 【求助】matlab如何只读取图像的一部分？
已经有11人回复
- 【求助】matlab中如何将程序的结果输出到一个文件中
已经有11人回复
- 求助一个差分方程的MATLAB实现
已经有3人回复
- 【求助】关于matlab的高维数组
已经有6人回复
- 【求助】LMI求解-可行性解的问题
已经有18人回复
- 【求助】初学matlab
已经有11人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
- 【求助】用matlab编过有限差分法的请进
已经有9人回复
2楼2013-05-15 20:07:30
fxlily2008
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3438.7
- 散金: 10
- 帖子: 141
- 在线: 210.9小时
- 虫号: 1680266
- 注册: 2012-03-10
- 专业: 控制理论与方法
3楼2013-05-15 21:09:02
fxlily2008
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3438.7
- 散金: 10
- 帖子: 141
- 在线: 210.9小时
- 虫号: 1680266
- 注册: 2012-03-10
- 专业: 控制理论与方法
4楼2013-05-15 21:10:13
idnk
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1158
- 散金: 11
- 帖子: 240
- 在线: 100.8小时
- 虫号: 150845
- 注册: 2005-12-29
- 性别: GG
- 专业: 控制理论与方法
5楼2013-05-15 22:36:49
elvinsu
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1000.1
- 帖子: 118
- 在线: 218小时
- 虫号: 2466860
- 注册: 2013-05-16
- 性别: GG
- 专业: 控制理论与方法
6楼2013-05-18 11:56:08
fxlily2008
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3438.7
- 散金: 10
- 帖子: 141
- 在线: 210.9小时
- 虫号: 1680266
- 注册: 2012-03-10
- 专业: 控制理论与方法
7楼2013-05-18 14:44:53
elvinsu
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1000.1
- 帖子: 118
- 在线: 218小时
- 虫号: 2466860
- 注册: 2013-05-16
- 性别: GG
- 专业: 控制理论与方法
8楼2013-05-19 21:00:21
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
fxlily2008: 金币+5, ★★★很有帮助 2013-06-07 09:04:54
xiegangmai: 金币+2, 谢谢参与！ 2013-06-07 20:36:25
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 >> |
9楼2013-06-06 21:35:52
10