关于罚函数法
为什么下面的例题(1)可以得到计算结果,例题(2)却不可以呢?例题(1)
(i)编写 M 文件 test.m
function g=test(x);
M=50000;
f=x(1)^2+x(2)^2+8;
g=f-M*min(x(1),0)-M*min(x(2),0)-M*min(x(1)^2-x(2),0)...
+M*abs(-x(1)-x(2)^2+2);
(ii)在Matlab命令窗口输入
[x,y]=fminunc('test',rand(2,1))
例题(2)
(i)编写 M 文件 fun.m
function g=fun(x)
M=1000;
f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x*(1)*x(2)+2*x(2)+1);
g=f-M*min(1.5+x(1)*x(2)-x(1)-x(2),0)-M*min(-x(1)*x(2)-10,0)+M*abs(x(1)+x(2));
(ii)在Matlab命令窗口输入
[x,val]=fminunc('fun',rand(2,1))