matlab进退法-切线法求极值.docx

《matlab进退法-切线法求极值.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、7、用进退法求极值：；a=0.5，h=0.5,ε=0.01。%建立函数f（x）。functiony=f(x)ifnargin==1y=x^2-10*x-36;endend%建立用进退法求解极值の程序。functiony=jintui(a,h,e)a=input('enterthevalueofa:');h=input('enterthevalueofh:');e=input('enterthevalueofe:');m=1;n=1;x0=a;x1=x0+h;whileabs(x1-x0)>eiff(x0)>f(x1)x0=x1;x1=x0+h*2^m;m=m+1;elsex1=x0;x

2、0=x1-h/2^n;n=n+1;endendx=(x0+x1)/2;y=f(x);end运行结果如下：enterthevalueofa:0.5enterthevalueofh:0.5enterthevalueofe:0.01ans=-59.99228、用切线法求极值：；%建立原函数；functiony=f(x)ifnargin==1y=x^2-16*x+35;endend%建立一阶导函数；functionyd=fd(x)ifnargin==1yd=2*x-16;endend%建立二阶导函数；functionydd=fdd(x)ifnargin==1ydd=2;endend%切线法求极

3、值函数；%x0为给定の初值，n为计算次数；%f（x）、fd（x）、fdd（x）分别为原函数，一阶导函数，二阶导函数；function[xj,yj]=qiexian(x0,n)x0=input('Pleaseenterthevalueofx0:');n=input('Pleaseenterthevalueofn:');fork=1:1:nx1=x0-fd(x0)/fdd(x0);x0=x1;endxj=x0;yj=f(xj);end此题中，自定义初值x0=2，n=30；xj=8，yj=-29；