[raw_data, text] = xlsread('data.xlsx');

consumption=raw_data(:,1);

investment=raw_data(:,2);

output=raw_data(:,3);

hours=raw_data(:,4);

[trend_c,c]=one_sided_hp_filter_kalman(consumption,1600);

[trend_i,i]=one_sided_hp_filter_kalman(investment,1600);

[trend_y,y]=one_sided_hp_filter_kalman(output,1600);

[trend_n,n]=one_sided_hp_filter_kalman(hours,1600);


save data c i y n;