查看: 1102|回复: 3

matlab优化问题

[复制链接]

字体大小: 正常放大

simonsnow

2 主题	7 听众	380 积分

234

升级 26.67%

TA的每日心情

	擦汗 2015-10-26 18:13

签到天数: 8 天

[LV.3]偶尔看看II

电梯直达

1^#

发表于 2015-7-23 20:43 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

10体力

本帖最后由 simonsnow 于 2015-7-23 22:28 编辑

这个问题用matlab怎么解啊，大神路过麻烦帮下忙，提供一下思路也好，谢谢啦
360反馈意见截图16341022282553.png

最佳答案

shengshengchina 查看完整内容

再发一遍 lingo程序： sets: a/1..4/:x,c; aa(a,a): Q; endsets data: c=6 8 4 2; Q=-1 2 0 0 2 -1 2 0 0 2 -1 2 0 0 2 -1; enddata max=@sum(a:c*x)+1/2*@sum(a(j):x(j)*@sum(a(i):x(i)*Q(i,j))); x(1)*x(2)+x(3)*x(4)>-1; x(1)*x(2)+x(3)*x(4)-3; @sum(a:x)

zan

lingo, 优化

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

相关帖子

234

使用道具举报

shengshengchina

1 主题	3 听众	309 积分

升级 3%

TA的每日心情

	开心 2016-6-7 09:46

签到天数: 53 天

[LV.5]常住居民I

自我介绍: 200 字节以内<br />
不支持自定义 Discuz! 代码

群组: 数学建模

群组: 2011年第一期数学建模

2^#

发表于 2015-7-23 20:43 |只看该作者 |招呼Ta 关注Ta

再发一遍

lingo程序：

sets:
a/1..4/:x,c;
aa(a,a):  Q;
endsets
data:
c=6 8 4 2;
Q=-1 2 0 0
2 -1 2 0
0 2 -1 2
0 0 2 -1;
enddata
max=@sum(a:c*x)+1/2*@sum(a(j):x(j)*@sum(a(i):x(i)*Q(i,j)));
x(1)*x(2)+x(3)*x(4)>-1;
x(1)*x(2)+x(3)*x(4)<1;
@sum(a:x)>-3;
@sum(a:x)<2;
@for(a :  @free(x);@bnd(-1,x,1));

结果：
  Global optimal solution found.
Objective value:                            16.33333
Objective bound:                            16.33334
Infeasibilities:                            0.000000
Extended solver steps:                            11
Total solver iterations:                         706

                     Variable          Value       Reduced Cost
                        X( 1)       1.000000          -3.000000
                        X( 2)       1.000000          -6.333333
                        X( 3)    0.6666667          0.000000
                        X( 4)    -0.6666667          0.000000

我也说一句

发表

使用道具举报

shengshengchina

1 主题	3 听众	309 积分

升级 3%

TA的每日心情

	开心 2016-6-7 09:46

签到天数: 53 天

[LV.5]常住居民I

自我介绍: 200 字节以内<br />
不支持自定义 Discuz! 代码

群组: 数学建模

群组: 2011年第一期数学建模

3^#

发表于 2015-7-24 16:34 |只看该作者 |招呼Ta 关注Ta

lingo程序：

sets:
a/1..4/:x,c;
aa(a,a)

;
endsets
data:
c=6 8 4 2;
Q=-1 2 0 0
2 -1 2 0
0 2 -1 2
0 0 2 -1;
enddata
max=@sum(a:c*x)+1/2*@sum(a(j):x(j)*@sum(a(i):x(i)*Q(i,j)));
x(1)*x(2)+x(3)*x(4)>-1;
x(1)*x(2)+x(3)*x(4)<1;
@sum(a:x)>-3;
@sum(a:x)<2;
@for(a

free(x);@bnd(-1,x,1));

结果：
  Global optimal solution found.
  Objective value:                            16.33333
  Objective bound:                            16.33334
  Infeasibilities:                            0.000000
  Extended solver steps:                            11
  Total solver iterations:                         706

                     Variable          Value       Reduced Cost
                        X( 1)       1.000000          -3.000000
                        X( 2)       1.000000          -6.333333
                        X( 3)    0.6666667          0.000000
                        X( 4)    -0.6666667          0.000000