matlab - MATLAB fsolve 超过 100 个时间段

Question

我正在尝试解决超过 100 个时间段的方程组。我需要找到 k1, k2, ... 和 c1, c2, ... y1, y2 等等。我有下面时间零周期的方程。我也知道 k100 = k_ss 由代码解决。我如何使用 fsolve 获取周期 1 到 100 之间的 k、c 和 y 值。谢谢！

% Parameters
alpha=1/3;          
beta=0.9;      
delta=.05;      
A=100;
k0=5;


k_ss = A*(alpha/((1/beta) - (1 - delta))).^(1/(1-alpha)); % this is also k100

y0 = A.^(1- alpha)*k0.^(alpha);
k1= @(c0, k1) y0+((1-delta)*k0)-c0
c1 = @(c0, k1) c0*beta*((1 - delta) + alpha*A.^(1-alpha)*k0.^(alpha-1));

标签： matlabsolver