Introduction - If you have any usage issues, please Google them yourself
f=10*(x(1)+x(2)-5)^2+(x(1)-x(2))^2
主程序dfp11.m如下:
oldoptions=optimset( fminunc )
options=optimset(oldoptions, LargeScale , off ) 采用中型优化算法
options11=optimset(options, HessUpdate , dfp ) DFG法,混合二、三次插值
[x11,fval11,exitflag11,output11]=fminunc( fun ,[0,2.5],options11)