嗨,劳拉,
你在正确的道路,但我就会相应的给你一些建议:
1)。当使用函数“fcontour”情节的轮廓一个函数,MATLAB试图提供函数的输入(x, y)向量。这是可能的,底层函数定义内的表达式应该支持向量操作。万博1manbetx在上面你的定义中,您使用表达式“x。^ 2的支持基础是一万博1manbetx个矢量,从而提高向量的每个组件x ' 2,返回一个包含结果的向量。然而,表达“x * y”不能应用于两个向量相同的形式(两行或两列向量)。允许操作“x * y”应用于两个向量elementwise,请考虑增加一个句号”。在相应的操作,即x。* y '。这将使底层的计算效率更高。在你的情况中,你可以改变你的函数定义如下,
fcontour (@ (x, y) 0.9097 + (-0.2332 . * x) + 0.2645。* y + (-1.07) * x ^。2 + (-0.02072)。* x。* y + 0.5786。* y。^ 2);
fcontour (@ (x, y) 0.8906 - 9.46 e + 12。* x + 9.46 e12汽油。* y + 5.20 e + 13。* x。^ 2 + 0.03807。* y。* x + 5.208 e + 13。* y ^ 2);
2)。你的目标函数无约束,表情如“x ^ 2”。后者的表达实际上是等于“x ^ 0.2”,导致的结果如果“x”是负的。因此,它可能是可取的0定义为下界设计变量,目标函数没有实现虚值,当然取决于您的应用程序和您的需要。我已经习惯你的脚本为0账户也下界,即
% %清晰的输出窗口和工作空间
clc,清晰;
% %目标函数的等值线
fcontour (@ (x, y) 0.9097 + (-0.2332 . * x) + 0.2645。* y + (-1.07) * x ^。2 + (-0.02072)。* x。* y + 0.5786。* y。^ 2);
fcontour (@ (x, y) 0.8906 - 9.46 e + 12。* x + 9.46 e12汽油。* y + 5.20 e + 13。* x。^ 2 + 0.03807。* y。* x + 5.208 e + 13。* y ^ 2);
% %定义帕累托优化选项和解决面前的问题
%选择= optimoptions (gamultiobj) ';
选择= optimoptions (“gamultiobj”,“PopulationSize”,61,“ParetoFraction”,0.71);
健身= @simpleMultiObjective;
据nvar = 2;
磅= 0 (2,1);
[x, fval] = gamultiobj(健身、据nvar、[]、[][],[],磅,[],[],选择);
图
情节(fval (: 1) fval (:, 2),“o”)
% %定义多目标函数
函数z = simpleMultiObjective (xy)
x = xy (1);y = xy (2);
z (1) = 0.9097 + (-0.2332 . * x) + 0.2645。* y + (-1.07) * x ^。2 + (-0.02072)。* x。* y + 0.5786。* y。^ 2;
z (2) = 0.8906 - 9.46 e + 12。* x + 9.46 e12汽油。* y + 5.20 e + 13。* x。^ 2 + 0.03807。* y。* x + 5.208 e + 13。* y ^ 2;
结束
请注意,我有增加的值“PopulationSize”和“ParetoFraction”优化选项时使用“gamultiobj”函数以获得一个密集的帕累托。还的更多信息,请访问下面的文档页面,
这样,下面的帕累托前后者运行脚本时,就会产生
我希望这可以帮助。
亲切的问候,
安德烈亚斯
