多维搜索算法程序.docVIP

拟牛顿算法 算法流程图/算法步骤 1、DFP算法： （1） BFGS算法 （1） 模式搜索算法 （1） 程序代码 1、DFP算法 %拟牛顿法中DFP算法求解f=x1*x1+2*x2*x2-2*x1*x2-4*x1的最小值 tic format long syms x1 x2 r; f= x1*x1+2*x2*x2-2*x1*x2-4*x1; %要最小化的函数 x=[x1,x2]; df=jacobian(f,x); %函数f的偏导 df=df.; x0=[0,0]; g1=subs(df,x,x0); %起始点的梯度 epsilon=1e-6; %0.000001为搜索精度 k=0; H0=[1 0;0 1]; %初始矩阵为二阶单位阵 while(norm(g1)epsilon) %迭代终止条件||g1||=epsilon if k==0 d=-H0*g1; %负梯度方向 else H1=H0+pk*pk/(qk*pk)-H0*qk*qk*H0/(qk*H0*qk); d=-H1*g1;H0=H1; end z=subs(f,x,x0+r*d); %关于r的函数 result1=jintuifa(z,r); %进退法求下单峰区间 result2= gol