1

我正在解决我的宏观经济学课上的一个问题。考虑以下等式：

在此处输入图像描述

这里，k 是固定的，c(k) 是通过 Matlab 中的“interp1”函数定义的。这是我的代码：

beta = 0.98;
delta = 0.13;
A = 2;
alpha = 1/3;
n_grid = 1000;   % Number of points for capital
k_grid = linspace(5, 15, n_grid)';
tol = 1e-5;
max_it = 1000;

c0 = ones(n_grid, 1);
new_k = zeros(n_grid, 1);
dist_c = tol + 1;
it_c = 0;

while dist_c > tol && it_c < max_it
   c_handle = @(k_tomorrow) interp1(k_grid, c0, k_tomorrow, 'linear', 'extrap');

   for i=1:n_grid
       % Solve for k'
       euler = @(k_tomorrow) (1/((1-delta)* k_grid(i) + A * k_grid(i)^alpha - k_tomorrow)) - beta*(1-delta + alpha*A*k_tomorrow^(alpha - 1))/c_handle(k_prime);
       new_k(i) = fzero(euler, k_grid(i));    % What's a good guess for fzero?    
   end

   % Compute new values for consumption
   new_c = A*k_grid.^alpha + (1-delta)*k_grid - new_k;

   % Check convergence
   dist_c = norm(new_c - c0);
   c0 = new_c;
   it_c = it_c + 1;

end

当我运行此代码时，对于某些索引 $i$，它运行良好并且 fzero 可以找到解决方案。但是对于索引，它只返回 NaN 并退出而不找到根。这是经济学中一个表现良好的问题，我们正在寻找的解决方案确实存在，并且我尝试实现的算法保证可以工作。但是我在 MATLAB 中解决这个问题没有太多经验，我想我在某个地方犯了一个愚蠢的错误。关于如何进行的任何想法？

这是典型的错误信息：

退出 fzero：中止对包含符号更改的区间的搜索，因为在搜索期间遇到了复杂的函数值。（-2.61092 处的函数值为 0.74278-0.30449i。）检查函数或使用不同的起始值重试。

提前非常感谢！

4

1 回答 1

1

唯一能产生复数的项是

k'^(alpha - 1) = k'^(-2/3)

您可能希望根据立方根的实际变体得到结果，您可以得到

sign(k') * abs(k')^(-2/3)

或更一般地，避免除以零

k' * (1e-16+abs(k'))^(alpha - 2)
于 2020-04-17T05:48:34.480 回答