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

%%  目的：「インパルス応答関数」の求め方の例として、
%%  最適化不要・ジャンプ変数なしの、ソロー・スワンモデルを使う。
%%  ただし、「本物」のソロー・スワンモデルと違って、線形の生産関数を仮定。
%%  このコードは、第0期にbがpermanentに10%だけ上昇したとして、
%%  その後内生変数がどのように変化するかを分析する。
%%  （ショックに対する内生変数の反応を通時的に追ったものをインパルス応答関数という)
%%  これはループを組む練習問題にもなっている。

clear all;

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

% 以下のパラメーター値を適宜設定。
s=0.2 ; %貯蓄率
d=0.05 ; %減耗率
a=0.1 ; %生産関数の傾き
b=0.15 ; %生産関数の切片

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

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