查看: 2417|回复: 0

菜鸟求教怎么调对这个程序。。。内点法求极小值x(1)-1)^2+x(2)^2 x1>=0,x2>=0

字体大小: 正常放大

7 主题	8 听众	317 积分

升级 5.67%

TA的每日心情

	慵懒 2016-10-11 09:39

签到天数: 173 天

[LV.7]常住居民III

自我介绍: 想学习数学建模

群组: 数学建模培训课堂2

群组: 2011年第一期数学建模

群组: 数学建模培训课堂1

电梯直达

1^#

发表于 2012-12-9 22:40 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

function mytest
global x
r(1)=10;
c=0.1;
x0=[3,3];
e=0.0001;
for k=1:20
[x,feval,exitflag,output]=fminunc(myfun1,x0);
if r(k)*(-log(x(1,1))-log(x(2,1)))<= e
      x
      k
      break
else
      r(k+1)=c*r(k);
end
end
function f=myfun1(x)
global x r
f=(x(1)-1)^2+x(2)^2-r(k)*log(x(1))-r(k)*log(x(2));

zan

程序