嵌套输出函数的示例
创建和使用输出函数中的示例不需要输出函数在每次迭代后保留数据。如果不需要保存各次迭代之间的数据,可以将输出函数编写为函数文件,并从命令行直接调用优化函数。但是,要使输出函数在每次迭代后记录数据,请编写一个用于实现以下目的的文件:
以嵌套函数的形式包括输出函数 - 有关详细信息,请参阅 MATLAB 编程基础中的嵌套函数。
调用优化函数。
在以下示例中,函数文件还包含目标函数作为局部函数。您也可以将目标函数编写为单独的文件或匿名函数。
嵌套函数可以访问其所在的文件中的变量。因此,此方法使输出函数能够在每次迭代后保留变量。
以下示例使用输出函数记录以下求解中的 fminsearch 迭代:
minxf(x)=minxex1(4x12+2x22+x1x2+2x2).
输出函数以矩阵(称为 history)的形式返回点序列。
要运行此示例,请执行下列步骤:
使用 MATLAB 编辑器打开一个新文件。
将以下代码复制并粘贴到此文件。
function [x fval history] = myproblem(x0)
history = [];
options = optimset('OutputFcn', @myoutput);
[x fval] = fminsearch(@objfun, x0,options);
function stop = myoutput(x,optimvalues,state);
stop = false;
if isequal(state,'iter')
history = [history; x];
end
end
function z = objfun(x)
z = exp(x(1))*(4*x(1)^2+2*x(2)^2+x(1)*x(2)+2*x(2));
end
end
将文件作为 myproblem.m 保存到 MATLAB 路径上的某个文件夹中。
在 MATLAB 提示符下,输入
[x fval history] = myproblem([-1 1]);
函数 fminsearch 返回最佳点 x 及 x 处的目标函数值 fval。
x,fval
x =
0.1290 -0.5323
fval =
-0.5689
此外,输出函数 myoutput 向 MATLAB 工作区返回矩阵 history,该矩阵包含算法在每次迭代中生成的点。history 的前四行为
history(1:4,:)
ans =
-1.0000 1.0000
-1.0000 1.0000
-1.0750 0.9000
-1.0125 0.8500
history 最后一行的点与最佳点 x 相同。
history(end,:)
ans =
0.1290 -0.5323
objfun(history(end,:))
ans =
-0.5689