标题:关于罚函数法
只看楼主
月下寒心
Rank: 1
等 级:新手上路
帖 子:98
专家分:0
注 册:2009-3-21
结帖率:97.5%
 问题点数:0 回复次数:0 
关于罚函数法
为什么下面的例题(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))
搜索更多相关主题的帖子: 罚函数法 
2009-08-12 15:15



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-282386-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 1.089072 second(s), 8 queries.
Copyright©2004-2025, BCCN.NET, All Rights Reserved