在函数体内调用 Octave 插值函数

Question

我正在尝试将 Octave 插值函数包装在函数体中，

function FUN = inter(p);
    FUN = interpn (x1, x2, x3, x4, x5, A, p(1), p(2), p(3), p(4), p(5), "spline");
end

我这样做的原因是因为我正在使用一个包，该函数需要一个字符串名称函数，在这种情况下是 packageFunction("inter", argument1);

问题现在正在调用，例如，

  disp("value = "), inter([10 2 4 3 4])

不起作用；没有看到向量错误：'x1' undefined ,

当然，向量 xi 和矩阵 A 是在函数体之上定义的。非常感谢这方面的建议，

谢谢，达米尔

---------- 在文件 example1.m 中

[a b c] = fminuit('gaussian','mnplot',[10 166 33],[x;y;dy])

---------- 在文件 gaussian.m 中

function f = gaussian(par,data);
%theoretical function
f = par(1)/(sqrt(2*pi)*par(3)) * exp(-.5*((data(1,:)- 
par(2))./par(3)).^2); 
if (size(data,1)==2),    %chi-square, error = 1
  f = sum((data(2,:) - f).^2);
 elseif (size(data,1)>2), %chi-square, error = 3rd row of data
 f = sum(((data(2,:) - f)./data(3,:)).^2);
end

Answer 1

鉴于您正在使用需要字符串作为函数的旧函数，下面的第一个解决方案将不起作用。然而，这是正确的做法。将旧函数更改为使用函数句柄而不是字符串将是我的首选解决方案。但是，您也可以使用下面的替代解决方案，它使用全局变量。这不是推荐的方法 (we should strive to avoid globals)，但可以解决您近期的问题。

正确方法：使用匿名函数

您应该使用匿名函数，这些函数可以在定义变量时捕获它们：

inter = @(p)interpn (x1, x2, x3, x4, x5, A, p(1), p(2), p(3), p(4), p(5), "spline");

现在inter(p) 就像inter 被声明为一个普通函数一样工作。但是x1、x2等在定义inter时定义的值将存储在inter中。

如上所述，您传递inter 的函数必须编写为接受函数句柄。

糟糕的快速解决方案：使用全局变量

首先，创建一个文件inter.m，内容如下：

function FUN = inter(p);
    global x1 x2 x3 x4 x5 A
    FUN = interpn (x1, x2, x3, x4, x5, A, p(1), p(2), p(3), p(4), p(5), "spline");
end

接下来，在您调用 inter 的脚本函数中，再次声明全局变量（目前 MATLAB 警告您应该在给它们赋值之前将它们声明为全局变量，在未来的版本中将需要这样做）：

global x1 x2 x3 x4 x5 A
x1 = ...
x2 = ...
% etc

inter([10 2 4 3 4])
% or:
fminuit('inter',...)