大连理工优化方法大作业MATLAB编程15页.docx

function x,dk,k=fjqx(x,s)flag=0;a=0;b=0;k=0;d=1;while(flag=0)p,q=getpq(x,d,s);if (p=0)&(q=0)dk=d;x=x+d*s;flag=1;end k=k+1; if(p=0)&(q0)a=d;d=min2*d,(d+b)/2;endend%定义求函数值的函数fun，当输入为x0=（x1，x2）时，输出为ffunction f=fun(x) f=(x(2)-x(1)2)2+(1-x(1)2;function gf=gfun(x) gf=-4*x(1)*(x(2)-x(1)2)+2*(x(1)-1),2*(x(2)-x(1)2);function p,q=getpq(x,d,s) p=fun(x)-fun(x+d*s)+0.20*d*gfun(x)*s