%%  一橋大学「金融経済論I」2010年度、担当：塩路悦朗
%%  ソロー。スワンモデル %% 
%%  ファイル名：solowswan2.m
%%  作成：塩路悦朗
%%  バージョン1、2010年10月26日午前1時30分

%%  目的：「インパルス応答関数」の求め方の例として、
%%  最適化不要・ジャンプ変数なしの、ソロー・スワンモデルを使う。
%%  「本物」のソロー・スワンモデルを対数線形近似したものを用いる。
%%  このコードは、第0期にlnAがpermanentに1単位だけ上昇したとして、
%%  (つまりAがもとの値から1%増加したとして）
%%  その後内生変数がどのように変化するかを分析する。

clear all;

%% ------------------- [1] Parameter proc  ------------------------

% 以下のパラメーター値を適宜設定。
s=0.2 ; %貯蓄率
d=0.05 ; %減耗率
A=0.25 ; %全要素生産性＝TFPの元の定常値
alpha=0.3 ; %資本のシェア(資本分配率)

%% --------------------- [2] Steady State proc >> -----------------------
% 定常状態の値を計算する。
Kstar=(s*A/d)^(1/(1-alpha)) ; %定常状態における資本ストック
    
%% ------------------ [3]  SIMULATION proc  ----------------- %%    
    
% TIME&INITIAL VALUES 

     t = 200;			% 何期後の反応まで求めたいか
     S1=0 ;             % 初期において経済はもともとの定常状態にあったとする。
                        % k=ln(K)-ln(K*)と定義
     a=1 ;              % そこで突然lnAの値が1増加した。a=ln(A)-ln(A*)と定義
 
% SIMULATION 
      S = zeros(t,1) ;  %以下でkの反応を求める
      S(1)=S1;          %ショックがあった期のK
      Ss = S1;
        for i = 2:t ;
          q = (1-d+alpha*d)*Ss+d*a  ;
   	      S(i) = q;
          Ss = S(i);
     	end
      
% DRAWING FIGURES
   set(0,'DefaultLineLineWidth',1.5) ; %図の線の太さを変える
   figure(1) 
       plot(1:t,S) ;
       ylabel('変数')
       xlabel('期間');
       title('ln(K)のインパルス応答関数(元の定常値からの乖離で表す)','FontSize',14) ;
