clear all
close all

x = load('cz.txt'); % load the data
t = x(:,1);
usd = x(:,3);
cpi = x(:,4);
cpi_gr = x(:,5);
dem = x(:,6);
pribor = x(:,7);

figure
plot(t,dem)
title('kuzr czk/dem')

kos = dem .^0.65 .* usd .^0.35;

figure
plot(t,kos)
title('efektivni smenny kurz')
grid

q = dem .* cpi_gr ./ cpi;

figure
plot(t,q)
title('realny smenny kurz')
grid

pie12 = cpi(13:end) ./ cpi(1:end-12) - 1;
pie12 = [ones(12,1)*nan;pie12];

lcpi = log(cpi);

pie12L = lcpi(13:end) - lcpi(1:end-12);
pie12L = [nan(12,1);pie12L];

figure
plot(t,[pie12 pie12L]*100)
legend('inflace CPI','inflace CPI - log aproximace')
ylabel('%')
grid

pie1 = cpi(2:end) ./ cpi(1:end-1) - 1;
pie1 = [nan;pie1]*12;

figure
plot(t,[pie12 pie1]*100)
legend('mezirocni inflace','mezimesicni inflace  - anualizovana')
grid

pribor = pribor/100;
rir = pribor(1:end-12) - pie12(13:end);
rir = [rir;nan(12,1)];

figure
plot(t,rir*100)
title('realna urokova mira')
ylabel('%')
grid






















    












