-
Notifications
You must be signed in to change notification settings - Fork 2
/
Insert_App_2_loop_PLT.m
33 lines (26 loc) · 1.01 KB
/
Insert_App_2_loop_PLT.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
% Application 2 - New Keynesian model in Brendon et al. (2013). For analysis of the M matrix
% in the case of a PLT interest rate rule. We 'comment out' certain parameters which are set in the 'loop' file.
%Calibration
betta = 0.99;
X1_min = betta-1;
%sigma = 1;
kappa = (1-0.85)*(1-0.85*betta)*(2+sigma)/0.85;
rho_i = 0.4;
%thetta_pi = 1.5;
thetta_p = thetta_pi;
thetta_y = 0;
%thetta_dy = 1.6;
%Matrices for xstar
F = [0 1 zeros(1,13)]; G = 0; H = 0;
%Reference regime
B1 = [1 -1 0 0 0; 0 1 -(1-rho_i)*(thetta_y+thetta_dy) 0 -(1-rho_i)*thetta_p; 1/sigma 0 1 0 0; 0 0 -kappa 1 0; 0 0 0 -1 1];
B2 = [zeros(2,5); 0 0 1 1/sigma 0; 0 0 0 betta 0; zeros(1,5)];
B3 = [0 0 0 0 0; 0 rho_i -(1-rho_i)*thetta_dy 0 0; zeros(2,5); 0 0 0 0 1];
B4 = [0;0;1;0;0];
B5 = [0;0;0;0;0];
%Alternative regime
B1_tild = [1 0 0 0 0; 0 1 -(1-rho_i)*(thetta_y+thetta_dy) 0 -(1-rho_i)*thetta_p; 1/sigma 0 1 0 0; 0 0 -kappa 1 0; 0 0 0 -1 1];
B2_tild = B2;
B3_tild = B3;
B4_tild = B4;
B5_tild = [X1_min; 0; 0; 0; 0];