example2_18.m

来自「Matlab 6.X辅助优化计算与设计----源代码」· M 代码 · 共 15 行

options=optimset('LargeScale','off');
options=optimset(options,'MaxFunEvals', 300);
options=optimset(options,'HessUpdate','dft');
%options=optimset(options,'LineSearch','cubicpoly');
% 采用本行参数设置时，表示计算方法为三次样条内插搜索法
% 否则为混合多项式内插搜索法
%以下代码同上
x0=[-1.9,2]
GRAD=[100*(4*x(1)^3-4*x(1)*x(2))+2*x(1)-2; …
100*(2*x(2)-2*x(1)^2); banplot2(x)]';
f='100*(x(2)-x(1)^2)^2+(1-x(1))^2';
options = optimset(options,'gradobj','on');
disp('[x,fval,exitflag,output] = fminunc({f,GRAD},x,options);');
[x,fval,exitflag,output] = fminunc({f,GRAD},x0,options);