matlab - 此 Matlab 函数的主体有效，但函数本身无效(interp1 错误)

Question

我已经为 Euler 策略迭代算法编写了以下子代码(带有注释的参数)。当我尝试运行函数体(全局以下的所有内容)时，例如 a1 = 1，它可以工作，并返回一个标量。但是，当我将该函数调用为 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)