此 Matlab 函数的主体有效，但函数本身无效（interp1 错误）

Question

我为 Euler 策略迭代算法编写了以下子代码（带有注释的参数）。当我尝试运行函数的主体（全局以下的所有内容）时，比如 a1 = 1，它起作用，并且 returns 一个标量。但是，当我将函数调用为 euler_diff_test(1) 时，出现错误。（粘贴在下方）

function diff = euler_diff_test(a1)
%the following comments are example parameters. They are in the global line originally.
% r = 0.2, a = 0.5, y = 1.1, a_grid = linspace(0.5,7,100)
%policy_guess = zeros(2,N);
%policy_guess(1,:) = 0.3*a_grid;
%policy_guess(2,:) = 0.3*a_grid;
% M = zeros(2,2); %M for markov transition kernel
% M(1,1) = p;
% M(2,2) = p;
% M(2,1) = 1-p;
% M(1,2) = 1-p;
% j = 1
global r a y a_grid policy_guess M j;

c = (1+r)*a + y - a1; %consumption formula

if c<=1e-02  %don't care about consumption being negative
    
    diff = 888888888888888888888;
else
    policy_func = interp1(a_grid', policy_guess', a1, 'linear');

    diff = 1/c - beta*(1+r)*(1 ./ policy_func)*M(j,:)'; 
end

end

错误读取：非常感谢任何帮助！

Answer 1

问题是您不了解全局变量，也不了解它们的工作原理！

您似乎在做类似的事情：

N=100; p=0.1;
r = 0.2, a = 0.5, y = 1.1, a_grid = linspace(0.5,7,100)
policy_guess = zeros(2,N);
policy_guess(1,:) = 0.3*a_grid;
policy_guess(2,:) = 0.3*a_grid;
M = zeros(2,2); %M for markov transition kernel
M(1,1) = p;
M(2,2) = p;
M(2,1) = 1-p;
M(1,2) = 1-p;
euler_diff_test(1)

这是导致您显示的错误的原因。当然是！

首先，您需要了解什么是全局以及什么是工作空间。每个功能都有自己的工作空间或“范围”。这意味着函数本身只能看到工作区内定义的变量。

全局变量是所有工作空间都存在的，任何人都可以修改。您似乎希望在函数外部定义所有这些变量，在您的函数内部。但是意识到！定义变量时，它们不是全局的。函数开始，在它的第一行，它只知道 a1 的存在。然后，稍后，您将一堆变量定义为全局变量，函数不知道。那么这个函数有什么作用呢？只为您将它们创建为空。

如果您希望在主脚本范围内创建的变量是全局变量，则需要将它们声明为全局变量，而不是在函数内部。因此，从函数中删除行 global ...，并将其放在声明所有变量的脚本之上，即在

之上

% here! 
N=100; p=0.1;
...

在我的例子中。

现在，重要的是：全局变量不好。当你有全局变量时，你不知道谁修改了，而且很容易忘记它们发生了什么，因为每个使用变量 a 的函数都会修改全局变量 a，所以调试起来很痛苦。几乎没有人因此而使用全局变量。最好的方法是将它们作为输入传递给函数，即将函数定义为：

function diff = euler_diff_test(a1,r, a, y, a_grid, policy_guess, M, j)

此 Matlab 函数的主体有效，但函数本身无效（interp1 错误）

Body of this Matlab function works, but not the function itself (interp1 error)

matlab

interpolation

numerical-methods