clear all
load('dadoscompletos1.mat');
data=dadoscompleto;
%i_c_a, i_c_f, i_c_t, fbk, m, x, c, s,j_r, e_r, i;
data=data(:,1:14);
data_C=data(:,1);
data_I=data(:,2);
data_Q=data(:,3);
data_e_r=data(:,4);
data_j_B=data(:,5);
data_r_ib=data(:,6);
data_r_bw=data(:,7);
data_pie_star=data(:,8);
data_pie=data(:,9);
data_b_i=data(:,10);
data_b_iw=data(:,11);
data_Y=data(:,12);
data_Bstar=data(:,13);
data_spr_b=data(:,14);

time=2007.02:0.25:2016.75;
time=time.';
T=size(time,1);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure;
hold on
plot (data_r_bw);
title('Juros consignados pessoal');
xlabel('Trimestral');
ylabel('Juros consigando pessoal (a.a.)');
hp_r_bw=hpfilter(data_r_bw, 1600);
plot(hp_r_bw);
hold off
%data_hp_r_bw=destrend(ts,'linear', BP);
%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_r_ib);
title('Formacao bruta de capital');
xlabel('Trimestral');
ylabel('Em milhões');
hp_r_bi=hpfilter(data_r_ib, 1600);
plot(hp_r_ib);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_Bstar);
title('Exportações');
xlabel('Trimestral');
ylabel('Em milhões');
hp_Bstar=hpfilter(data_Bstar, 1600);
plot(hp_Bstar);


%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_I);
title('Crédito pessoal não consignado');
xlabel('Trimestral');
ylabel('Crédito pessoal não consignado R$(Milhões)');
hp_I=hpfilter(data_I, 1600);
plot(hp_I);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_b_iw);
title('Crédito pessoal consignado total');
xlabel('Trimestral');
ylabel('Crédito pessoal consignados total R$(Milhões)');
hp_b_iw=hpfilter(data_b_iw, 1600);
plot(hp_b_iw);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_Q);
title('Lucro liquido dos bancos');
xlabel('Trimestral');
ylabel('Lucros liquido dos bancos em milhoes');
hp_Q=hpfilter(data_Q, 1600);
plot(hp_Q);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_C);
title('Divida Externa Bruta');
xlabel('Trimestral');
ylabel('Divida Externa Bruta em milhoes');
hp_C=hpfilter(data_C, 1600);
plot(hp_C);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
%plot(data_y_exp);
%title('Taxa de juros SELIC');
%xlabel('Trimestral');
%ylabel('Taxa de juros SELIC');
%hp_y_exp=hpfilter(data_y_exp, 1600);
%plot(hp_y_exp);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_b_i);
title('Divida Externa Bruta');
xlabel('Trimestral');
ylabel('Divida Externa Bruta em milhoes');
hp_b_i=hpfilter(data_b_i, 1600);
plot(hp_b_i);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_j_B);
title('Taxa de cambio real');
xlabel('Trimestral');
ylabel('Taxa de cambio real');
hp_j_B=hpfilter(data_j_B, 1600);
plot(hp_j_B);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_e_r);
title('Taxa de cambio real');
xlabel('Trimestral');
ylabel('Taxa de cambio real');
hp_e_r=hpfilter(data_e_r, 1600);
plot(hp_e_r);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
%plot(data_b_iw);
%title('Taxa de cambio real');
%xlabel('Trimestral');
%ylabel('Taxa de cambio real');
%hp_b_iw=hpfilter(data_b_iw, 1600);
%plot(hp_b_iw);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_pie_star);
title('Taxa de cambio real');
xlabel('Trimestral');
ylabel('Taxa de cambio real');
hp_pie_star=hpfilter(data_pie_star, 1600);
plot(hp_pie_star);

%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_pie);
title('Taxa de cambio real');
xlabel('Trimestral');
ylabel('Taxa de cambio real');
hp_pie=hpfilter(data_pie, 1600);
plot(hp_pie);


%Primeira parte: Decomposicao da tendência e componente ciclico
%figure
plot(data_Y);
title('Taxa de cambio real');
xlabel('Trimestral');
ylabel('Taxa de cambio real');
hp_Y=hpfilter(data_Y, 1600);
plot(hp_Y);


%filtro HP
plot(data_spr_b);
title('Taxa do spread bancário');
xlabel('Trimestral');
ylabel('Taxa do spread bancário');
hp_spr_b=hpfilter(data_spr_b);
plot(hp_spr_b);

save ('data_hp.xlsx','hp_r_bw','hp_Bstar','hp_C','hp_I','hp_e_r','hp_Q','hp_j_B','hp_pie_star','hp_pie','hp_bi','hp_b_iw','hp_Y','hp_spr_b')

%Criando gráfico com todas as variáveis
subplot(3,3,1)
plot(hp_r_bw)
title('Juros consignados pessoal (a.a)')

subplot(3,3,2)
plot(hp_r_bi)
title('Formacao bruta de capital')

subplot(3,3,3)
plot(hp_Bstar)
title('Exportações')


subplot(3,3,4)
plot(hp_I)
title('Crédito pessoal não consignado R$(Milhões)')

subplot(3,3,5)
plot(hp_Q)
title('Crédito pessoal consignados total R$ (Milhões)')

subplot(3,3,6)
plot(hp_I)
title('Lucros liquido dos bancos em milhoes')

subplot(3,3,7)
plot(hp_C)
title('Divida Externa Bruta em milhoes')

subplot(3,3,8)
plot(hp_e_r)
title('Taxa de cambio real')

subplot(3,3,9)
plot(hp_j_B)
title('Taxa de juros SELIC')


suplot(2,2,1)
plot(hp_pie_star)
title('')

suplot(2,2,2)
plot(hp_pie)
title('')

suplot(2,2,3)
plot(hp_bi)
title('')

suplot(2,2,4)
plot(hp_b_iw)
title('')

suplot(2,2,5)
plot(hp_Y)
title('')

suplot(2,2,6)
plot(hp_spr_b)
title('')
