Solow模型推导模拟

Solow模型推导模拟

1 Solow模型推导

在存在资本折旧、技术进步和人口增长条件下，有效劳动人均资本为
$$\begin{array}{rl}\dot{k}(t)& =\frac{\dot{K}(t)}{A(t)L(t)}-\frac{K(t)}{[A(t)L(t){]}^2}[A(t)\dot{L}(t)+L(t)\dot{A}(t)]\\ & =\frac{\dot{K}(t)}{A(t)L(t)}-\frac{K(t)}{A(t)L(t)}\frac{\dot{L}(t)}{L(t)}-\frac{K(t)}{A(t)L(t)}\frac{\dot{A}(t)}{A(t)}\end{array}$$
其中$\dot{K}(t)=sY(t)-\delta K(t)$,$k(t)=K(t)/A(t)L(t)$。假设技术进步和人口增长率不变，分别为$g、n$，即
$$\dot{L}(t)/L(t)=n;\dot{A}(t)/L(t)=g$$
代入上式得到
$$\begin{array}{rl}\dot{k}(t)& =\frac{sY(t)-\delta K(t)}{A(t)L(t)}-k(t)n-k(t)g\\ & =s\frac{Y(t)}{A(t)L(t)}-\delta k(t)-nk(t)-gk(t)\end{array}$$
其中$f(k)=Y(t)/A(t)L(t)$，于是得到
$$\dot{k}(t)=sf(k(t))-(n+g+\delta )k(t)$$
上述方程是solow的关键方程。假设生产函数满足如下形式
$$Y=K^{\alpha }(AL{)}^{1-\alpha }$$
则有效劳动的平均收入$y(t)=Y/AL=k(t{)}^{\alpha }$,solow方程变为
$$\dot{k}(t)=sk(t{)}^{\alpha }-(n+g+\delta )k(t).$$

---

2 Solow模型模拟

假设一国家(地区)储蓄率$s=0.4$,资本贡献的份额$\alpha =0.4$，技术进步率$g$和人口增长率$n$均为0.01，资本折旧率为$\delta =0.03$。使用matlab求解微分方程$\dot{k}(t)=s k(t)^{\alpha}-(n+g+\delta) k(t) $的积分曲线。这里令初始$k(0)$从5到50取值不等。

s=0.4;
alpha=0.4;
n=0.01;
g=0.01;
delta=0.03;
F = @(t,k) s*k(1)^alpha -(n+g+delta)*k(1);
for k0 = 0:5:50
    [t,k]=ode45(F,[0 200],k0);
    plot(t,k,'k')
    grid on;
    hold on;
    xlabel("t");ylabel("k")
end

可以看出无论初始$k>0$为多少，最终都会都会收敛到稳态$k=31.984$，此时稳态$k$增长率为0。此外，也可以基于储蓄率变化如何影响稳态的有效劳动人均资本$k$。

% -----------------持平投资曲线------------------
plot(k,(n+g+delta)*k,"-k")
hold on
M1=find(abs(0.4*k.^alpha-(n+g+delta)*k)<0.001);
N1=(n+g+delta)*k(M1(2));
plot([k(M1(2)),k(M1(2))],[0,N1],"--r")
hold on
text(k(M1(2))+0.03,N1-0.03,'A')
M2=find(abs(0.5*k.^alpha-(n+g+delta)*k)<0.001);
N2=(n+g+delta)*k(M2(2));
plot([k(M2(2)),k(M2(2))],[0,N2],"--r")
text(k(M2(2))+0.05,N2-0.05,'B')
grid on;
xlabel('k');ylabel('sf(k),(n+\delta+g)k')
legend("s=0.4","s=0.5")

当储蓄率$s$提升，稳态$k$提升但后期也会收敛达到新的稳态。

---

-END-