%generazione dei dati sperimentali
%I "dati sperimentali" sono in questo caso dei dati prodotti da una
%simulazione con il flash virtuale (quindi e' la risposta al campione
%unitario), piu' del rumore aggiunto, ed un "h" prestabilito: e' da vedere
%se viene correttamente identificato dall'algoritmo di analisi dei dati.


corrosione=0;   %0= per generare dati campione sano
                %1= per generare dati campione corroso
dati_veri_dominio;
save L_sano Lungh;
calore_1d_fd; %=> T_vero_x0

%CAMPIONE NON CORROSO (es. n=50)
%da digitare per generare T_vero_x0:
%aggiunta del rumore per simulare  le misure sperimentali
std_err = 1e-3;
rumore = std_err*randn(1,N+1); %<=>  distribuzione N(0,var)
T_vero_x0=T_vero_x0+rumore; %<=>  distribuzione N(0,var)
save T_vero_x0 T_vero_x0;
disp('dati sperimentali relativi al campione non corroso: T_vero_x0.mat');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

corrosione=1;   %0= per generare dati campione sano
                %1= per generare dati campione corroso
dati_veri_dominio;
calore_1d_fd; %=> T_vero_x0_corr

%CAMPIONE CORROSO (quindi con un n minore: es. n=20)
%da digitare per generare T_vero_x0_corr:
%aggiunta del rumore per simulare  le misure sperimentali
std_err = 1e-3;
rumore = std_err*randn(1,N+1); %<=>  distribuzione N(0,var)
T_vero_x0_corr=T_vero_x0+rumore; %<=>  distribuzione N(0,var)
save T_vero_x0_corr T_vero_x0_corr;
disp('dati sperimentali relativi al campione non corroso: T_vero_x0_corr.mat');
