matlab黄金分割法求解.doc

《matlab黄金分割法求解.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第4次作业用matlab求解f(t)=(t^2-1)^2+(t-1)^2+3[-10,10]最小植（黄金分割法）。x为最优的t值，minf为最优的函数值，f为函数，t为函数f的变量，a,b为下单峰区间[a,b],eps为精确度，M函数文件：function[x,minf]=minHJ(f,a,b,eps)formatlong;ifnargin==3eps=1.0e-6;Endl=a+0.382*(b-a);u=a+0.618*(b-a);a=-10k=1;tol=b-a;whiletol>eps&&k<100000fl=subs(f,f

2、indsym(f),l);%subs是置换函数iffl>fua=l;l=u;u=a+0.618*(b-a);elseb=u;u=l;l=a+0.382*(b-a);endk=k+1;tol=abs(b-a);endifk==100000disp('找不到最小值！');x=NaN;minf=NaN;return;endx=(a+b)/2;minf=subs(f,findsym(f),x);formatshort;窗口M文件调用：>>symst>>f=(t.^2-1).^2+(t-1).^2+3f=(t^2-1)^2+(t-1)^2+3>>

3、[x,minf]=minHJ(f,-10,10,0.000001)a=-10x=1.0000minf=3.0000