分析
这是一个约束优化问题,可用fmincon求解。
参考代码
f = @(x)exp(8.864-267.893/abs(x(1)))+exp(11.271-1227.294/x(2))+exp(13.223-12277.294/x(3))+exp(11.925-3632.467/x(4));s = 10000;x = fmincon(@(x)-f(x),ones(1,4)*s/4,[],[],ones(1,4),s,zeros(1,4),ones(1,4)*inf,[]);x = num2str(x)f = f(x)
得到的最优解为:
x =
-3.3140659e-018 3088.6297 1.1046886e-018 6911.3703
f =
1.4204e+005
说明
1、最优解的x1、x3在浮点数误差意义上为0;
2、由于x1的结果是很小的负值,导致目标函数的y1为无穷大,所以该项的x1应加绝对值符号(x1-x4都应该是正数吧?不然任取哪一个为一个很小的正数,都会导致目标函数为无穷大)。
温馨提示:答案为网友推荐,仅供参考