在 Matlab 函数中使用方程的变换？

Question

我正在尝试编写一个非常基本的平均似然估计器函数，该函数依赖于优化工具箱的fminsearch()函数。 我尝试将函数句柄作为输入传递并优化否定的总和。 函数的日志。 我的代码是：

function [ params, max ] = routine( fun )
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here

[filename,path] = uigetfile('*.'); #To grab the csv file location
name = strcat(path,filename);
data = csvread(name);

lh = @(x) sum(-log( fun(x) )); <-------- ERROR LINE

options = optimset('Display', 'off', 'MaxIter', 100000, 'TolX', 10^-20, 'TolFun', 10^-20);
[theta, max1] = fminsearch(lh, [0,1], options);

params = theta
max = max1

end

中间的那一行给了我这个错误：

Undefined function or variable 'data'.

Error in @(x)(1/(sqrt(2*pi)*x(2)))*exp((-(data-x(1)).^2)/(2*x(2)^2))

Error in @(x)sum(-log(fun(x)),data)

Error in fminsearch (line 200)
fv(:,1) = funfcn(x,varargin{:});

现在，在这段代码之外，这工作得非常好，令人气愤。

pan = @(x)sum(-log((1/(sqrt(2*pi)*x(2)))*exp((-(data-x(1)).^2)/(2*x(2)^2))));
options = optimset('Display', 'off', 'MaxIter', 100000, 'TolX', 10^-20, 'TolFun', 10^-20);
[theta, max1] = fminsearch(pan, [0,1], options);

出于某种原因，按照我在中间所做的方式转换方程会使data 。 看不见。 我应该如何正确地将转换应用到我的方程以使此代码有效？

Answer 1

查看匿名函数的文档。 这是相关的部分

表达式中的变量

函数句柄不仅可以存储表达式，还可以存储表达式求值所需的变量。

例如，为需要系数a 、 b和c的匿名函数创建函数句柄。

 a = 1.3; b = .2; c = 30; parabola = @(x) a*x.^2 + b*x + c;

由于a 、 b和c在您创建parabola时可用，因此函数句柄包含这些值。 即使您清除了变量，这些值仍然存在于函数句柄中：

 clear abc x = 1; y = parabola(x)

y = 31.5000

要为系数提供不同的值，您必须创建一个新的函数句柄

你的情况的值data您阅读将无法使用fun ，因为函数句柄之前定义的data 。 相反，我们可以让data成为fun的参数，然后我们一起避免这个问题。

---

调用代码

myfun = @(x,data) (1/(sqrt(2*pi)*x(2)))*exp((-(data-x(1)).^2)/(2*x(2)^2));
[params,max] = routine(myfun);

---

例程.m （需要更改的位）

function [ params, max ] = routine( fun )
...
lh = @(x) sum(-log( fun(x,data) ));
...
end