clear all; clc


% load data_impluse_as_eA;
% load data_impluse_as_eP;
% load data_impluse_as_eO;

dynare IRF_as_no_policy_old_calibration_eA;
dynare IRF_as_no_policy_old_calibration_eP;
dynare IRF_as_no_policy_old_calibration_eO;

load data_impluse_as_eA_haoning;
load data_impluse_as_eP_haoning;
load data_impluse_as_eO_haoning;
load data_steady_state_as_haoning
Rb_ss = steady_state_as(7);


% h = legend('$R^{B}\,\mbox{(high \ensuremath{\overline{\omega}})}$','$R^{B}\,\mbox{(low \ensuremath{\overline{\omega}})}$',2);
% set(h,'Interpreter','latex','fontsize',14)



% impluse_as_eA_unconventional =[I_eA K_eA Y_eA C_eA qk_eA eps1_eA eps2_eA Rb_eA omegat_eA spread_eA S_eA haircut_eA N_eA Phi_k_tilde_eA l_ppi_eA l_omega_eA];
% %                               1    2    3    4    5     6        7      8     9          10      11     12        13  
% save data_impluse_as_eA_unconventional impluse_as_eA_unconventional;



figure(1)
subplot(3,4,1)
plot(1:20,impluse_as_eA(:,3)*111,'-b','LineWidth',2);
hold on
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$Y$','$I$','$C$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
ylabel('$Y$','Interpreter','latex','fontsize',14);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;
%set(gca,'xtick',[])
%NumTicks = 4;
%L = get(gca,'YLim');
%set(gca,'YTick',linspace(L(1),L(2),NumTicks))



subplot(3,4,3)
plot(1:20,impluse_as_eA(:,8)*105,'-b','LineWidth',2);
ylabel('$R^{B}$','Interpreter','latex','fontsize',14);
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$P^{K}$','$R^{B}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;


subplot(3,4,2)
plot(1:20,impluse_as_eA(:,5)*123,'-b','LineWidth',2);
ylabel('$P^{K}$','Interpreter','latex','fontsize',14);
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$P^{K}$','$R^{B}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;



subplot(3,4,4)
plot(1:20,impluse_as_eA(:,16)*1057,'-b',1:20,impluse_as_eA(:,15)*117,'--r','LineWidth',2);
hold on
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$R^{K}-R^{B}$','$r^{\pi}$','$r^{\omega}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;

% h2 = legend('$r^{\omega}$','$r^{\pi}$');
% set(h2, 'FontSize', 12,  'Interpreter','latex','Location', 'NorthEast', 'Box', 'on');

%set(gca,'xtick',[])
%NumTicks = 4;
%L = get(gca,'YLim');
%set(gca,'YTick',linspace(L(1),L(2),NumTicks))

subplot(3,4,5)
plot(1:20,impluse_as_eP(:,3)*119,'-b','LineWidth',2);
hold on
%plot(1:20,C_eP*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$Y$','$I$','$C$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
ylabel('$Y$','Interpreter','latex','fontsize',14);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;
%set(gca,'xtick',[])
%NumTicks = 4;
%L = get(gca,'YLim');
%set(gca,'YTick',linspace(L(1),L(2),NumTicks))



subplot(3,4,6)
plot(1:20,impluse_as_eP(:,5)*112,'-b','LineWidth',2);
ylabel('$P^{K}$','Interpreter','latex','fontsize',14);
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$P^{K}$','$R^{B}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;


subplot(3,4,7)
plot(1:20,impluse_as_eP(:,8)*117,'-b','LineWidth',2);
ylabel('$R^{B}$','Interpreter','latex','fontsize',14);
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$P^{K}$','$R^{B}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;




subplot(3,4,8)
plot(1:20,impluse_as_eP(:,16)*117,'-b',1:20,impluse_as_eP(:,15)*117,'--r','LineWidth',2);
hold on
%plot(1:20,C_eP*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;

h2 = legend('$r^{\omega}$','$r^{\pi}$');
set(h2, 'FontSize', 18,  'Interpreter','latex','Location', 'NorthEast', 'Box', 'on');

%h = legend('$r^{\omega}$','$r^{\pi}$',2);
%set(h,'Interpreter','latex','fontsize',12)
xlim([1,20]);
%set(gca,'xtick',[])
%NumTicks = 4;
%L = get(gca,'YLim');
%set(gca,'YTick',linspace(L(1),L(2),NumTicks))
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;




subplot(3,4,9)
plot(1:20,impluse_as_eO(:,3)*126,'-b','LineWidth',2);
hold on
%plot(1:20,C_eO*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$Y$','$I$','$C$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
ylabel('$Y$','Interpreter','latex','fontsize',14);
%set(gca,'xtick',[])
%NumTicks = 4;
%L = get(gca,'YLim');
%set(gca,'YTick',linspace(L(1),L(2),NumTicks))
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;



subplot(3,4,10)
plot(1:20,impluse_as_eO(:,5)*117,'-b','LineWidth',2);
ylabel('$P^{K}$','Interpreter','latex','fontsize',14);
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$P^{K}$','$R^{B}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;


subplot(3,4,11)
plot(1:20,impluse_as_eO(:,8)*117,'-b','LineWidth',2);
ylabel('$R^{B}$','Interpreter','latex','fontsize',14);
%plot(1:20,C_eA*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$P^{K}$','$R^{B}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;



subplot(3,4,12)
plot(1:20,impluse_as_eO(:,16)*117,'-b',1:20,impluse_as_eO(:,15)*117,'--r','LineWidth',2);
hold on
%plot(1:20,C_eO*117,'-.', 'LineWidth',2);
%title('$I_{t}$','Interpreter','latex','fontsize',14);
%grid on;
% h = legend('$R^{K}-R^{B}$','$r^{\pi}$','$r^{\omega}$',2);
% set(h,'Interpreter','latex','fontsize',14)
xlim([1,20]);
%set(gca,'xtick',[])
%NumTicks = 4;
%L = get(gca,'YLim');
%set(gca,'YTick',linspace(L(1),L(2),NumTicks))
set(gca,  'FontSize', 14, 'LineWidth', 1, 'Box', 'on','fontname','helvetica');
grid on;

% h2 = legend('$r^{\omega}$','$r^{\pi}$');
% set(h2, 'FontSize', 12,  'Interpreter','latex','Location', 'NorthEast', 'Box', 'on');



print -depsc2 irf_no_policy_all_new.eps
















% subplot(3,4,1)
% plot(1:20,impluse_as_eA(:,3)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eA*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% ylabel('$Y$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$Y$','$I$','$C$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,2)
% plot(1:20,impluse_as_eA(:,5)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eA*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% ylabel('$P_{K}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$P^{K}$','$R^{B}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% 
% subplot(3,4,3)
% plot(1:20,impluse_as_eA(:,8)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eA*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$P^{K}$','$R^{B}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% ylabel('$R_{B}$','Interpreter','latex','fontsize',14);
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,4)
% plot(1:20,impluse_as_eA(:,10)*117,'-k',1:20,impluse_as_eA(:,15)*117,'--b',1:20,impluse_as_eA(:,16)*117,'-.r','LineWidth',2);
% hold on
% %plot(1:20,C_eA*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$R^{K}-R^{B}$','$r^{\pi}$','$r^{\omega}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% ylabel('Spread','Interpreter','latex','fontsize',14);
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,5)
% plot(1:20,impluse_as_eP(:,3)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eP*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% ylabel('$Y$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$Y$','$I$','$C$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,6)
% plot(1:20,impluse_as_eP(:,5)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eP*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% ylabel('$P_{K}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$P^{K}$','$R^{B}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% 
% subplot(3,4,7)
% plot(1:20,impluse_as_eP(:,8)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eP*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$P^{K}$','$R^{B}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% ylabel('$R_{B}$','Interpreter','latex','fontsize',14);
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,8)
% plot(1:20,impluse_as_eP(:,10)*117,'-k',1:20,impluse_as_eP(:,15)*117,'--b',1:20,impluse_as_eP(:,16)*117,'-.r','LineWidth',2);
% hold on
% %plot(1:20,C_eP*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% %grid on;
% h = legend('$R^{K}-R^{B}$','$r^{\pi}$','$r^{\omega}$');
% set(h,'Interpreter','latex','fontsize',14)
% ylabel('Spread','Interpreter','latex','fontsize',14);
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% 
% subplot(3,4,9)
% plot(1:20,impluse_as_eO(:,3)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eO*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% ylabel('$Y$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$Y$','$I$','$C$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,10)
% plot(1:20,impluse_as_eO(:,5)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eO*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% ylabel('$P_{K}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$P^{K}$','$R^{B}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% 
% subplot(3,4,11)
% plot(1:20,impluse_as_eO(:,8)*117,'-b','LineWidth',2);
% hold on
% %plot(1:20,C_eO*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$P^{K}$','$R^{B}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% ylabel('$R_{B}$','Interpreter','latex','fontsize',14);
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
% subplot(3,4,12)
% plot(1:20,impluse_as_eO(:,10)*117,'-k',1:20,impluse_as_eO(:,15)*117,'--b',1:20,impluse_as_eO(:,16)*117,'-.r','LineWidth',2);
% hold on
% %plot(1:20,C_eO*117,'-.', 'LineWidth',2);
% %title('$I_{t}$','Interpreter','latex','fontsize',14);
% %grid on;
% % h = legend('$R^{K}-R^{B}$','$r^{\pi}$','$r^{\omega}$',2);
% % set(h,'Interpreter','latex','fontsize',14)
% ylabel('Spread','Interpreter','latex','fontsize',14);
% xlim([1,20]);
% set(gca,'xtick',[])
% %NumTicks = 4;
% %L = get(gca,'YLim');
% %set(gca,'YTick',linspace(L(1),L(2),NumTicks))
% 
%  print -depsc2 irf_no_policy_all_2.eps

