- 最后登录
- 2009-5-17
- 在线时间
- 1 小时
- 阅读权限
- 5
- 注册时间
- 2009-4-26
- 积分
- 1
- 帖子
- 1
- 精华
- 0
- UID
- 89505
- 性别
- 保密
- 积分
- 1
- 帖子
- 1
- 精华
- 0
- UID
- 89505
- 性别
- 保密
|
我尝试用matlab编程求解。
一下是两个程序的代码,
myfun.m:
function f=myfun(x)
w=[46 23 13 8 35 17 10 6 28 14 8 5];
p=[1/20 1/15 23/60 1/2];
p=[p./3 p./3 p./3];
f=sum(sum(x))-sum(sum(w.*p*x));
pro.m:
clear;
format long
A=[];b=[];Aeq=[];beq=[];lb=zeros(12,1);ub=[];
options=optimset('largescale','off','MaxFunEvals',100000);
x0=ones(12,1);
[x,fval]=fmincon(@myfun,x0,A,b,Aeq,beq,lb,ub,[],options);
j=1;
for i=1:3
X(i, =x([j:j+3],1)';
j=j+4;
end
X
fval
pro是主程序了,运行后结果:
Optimization terminated: first-order optimality measure less than options.TolFun
and maximum constraint violation is less than options.TolCon.
Active inequalities (to within options.TolCon = 1e-006):
lower upper ineqlin ineqnonlin
1
X =
1.0e+015 *
0 0.12738413925575 2.94298697982253 1.49385216880625
0.19736933048247 0.00000000000000 1.24881998421207 0.00000000000000
0.08275009100300 0.08674104694720 0.09892429527400 0.14148336004783
fval =
-2.520702076217168e+015
由于解决整数非线性问题较难,就可以通过这个来得出近似结果。
结果表明: M1 M2 R1 是押宝的重点。(根据你给的条件,实际可能不一样)
%%%%%%%%%%%%%%%%%%%%%%%%
魔方新手,请多指点
[ 本帖最后由 linxin_wise 于 2009-4-27 09:43 编辑 ] |
|