
close all
clear all
clc

% routine che mi tira fuori i dati di validazione di altre stagioni

load dati_inverno

l=[3:size(dati,3)];
y_inv1 = squeeze(dati(:,3,:))';
y_inv1=y_inv1(l,:);
u_2 = squeeze(dati(:,1,1))'; %T esterna
u_3 = squeeze(dati(:,1,2))'; %Umidità
u_4 = squeeze(dati(:,3,1))'; %solare normale diretto kW/m^2
u_5 = squeeze(dati(:,4,1))'; %solare orizzontale diffuso kW/m^2

%compongo l'ingresso totale con quelli desiderati (mettere ;)
u_1= [u_2; u_4;u_5];

%% Detrend: rimozione media
%tolgo la media ai dati con la funzione centra oppure è possibile usare la
%funzione detrend() con medesimi risultati.
[my, y_inv]=centra(y_inv1);
[mu, u_inv]=centra(u_1);

frazione = 2/3; 
N = length(y_inv);

% validazione
y_valinv = y_inv(:,ceil(N*frazione):end);
u_valinv = u_inv(:,ceil(N*frazione):end);

figure
subplot(2,1,1)
plot(u_valinv');
grid on;
title('Ingressi validazione inverno');
ylabel('°C (scalati)');
legend('Temp Ext','Umidità Ext')
subplot(2,1,2)
plot(y_valinv');
grid on;
title('Uscite validazione inverno');
ylabel('°C (scalati)');

d_valinv = iddata(y_valinv',u_valinv');
  
save d_valinv d_valinv






