黄金分割法、斐波那契法求极值.doc

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

1、５、用黄金分割法求极值：；%建立原方程函数；functiony=fx(x)ifnargin==1y=x+20/x;endend%a为区间下限，b为区间上限，e为精度；%fx（x）为原方程函数；function[xj,yj]=huangjin(a,b,e)a=input('Pleaseenterthevalueofa:');b=input('Pleaseenterthevalueofb:');e=input('Pleaseenterthevalueofe:');whileb-a>ex1=a+0.382*(b-

2、a);x2=a+0.618*(b-a);iffx(x1)

3、rgin==1y=[0];y(1)=[1];y(2)=[1];fork=3:1:ny(k)=[y(k-1)+y(k-2)];endfn=y(n);endend%求解应计算次数的函数；%s为（b-a）/e的值，其中（a，b）为单峰区间，e为精度；functionn=cishu(s)ifnargin==1n=1;whileF(n)

4、evalueofa:');b=input('Pleaseenterthevalueofb:');e=input('Pleaseenterthevaludofe:');s=(b-a)/e;n=cishu(s);x2=a+F(n-1)/F(n)*(b-a);x1=b-F(n-1)/F(n)*(b-a);while(b-a)>eiffx(x2)>=fx(x1)b=x2;x2=x1;x1=a+b-x2;elsea=x1;x1=x2;x2=a+b-x1;endendxj=(a+b)/2;yj=fx(xj);end此题

5、中，a=-10，b=10，e=0.01，程序运行结果为：xj=1.9930，yj=2.8845e-004，若原方程改变，只需改变原方程函数即可。