(完整word版)变尺度DFP算法MATLAB实例.docVIP

clear all;close all;clc; syms x1 x2t; f=4*(x1-5 )人 2+(x2-6 )人2; f_grad=[diff(f,x1);diff(f,x2)]; X0=[8;9]; n=10; epsonal=0.01; fx=inline(f); fx_grad=inline(f_grad); X=X0; Xk=X0; while 1 fx0=fx(X0(1),X0(2)); fx0_grad=fx_grad(X0(1),X0(2)); Hk=eye(2); Pk=-fx0_grad; k=0; while 1 param=Xk+t*Pk; fx1=fx(par