load dati_autoclave.txt
T=dati_autoclave;
% for i=101:125
%     T(i)=T(i)+6;
% end


%%%%%%%%%%%%%%%%%% ISTOGRAMMA DATI %%%%%%%%%%%%%%
figure(1)
clf
hist(T,336:2:364) % faccio un istogramma dei dati
xlabel('temperature')
ylabel('frequency')
title('Histogram of temperature');
grid
pause

%%%%%%%%%%%%%%%%% VERIFICA GAUSSIANITA' E INDIPENDENZA %%%%%%%%%%%%%

figure(1)
clf
qqplot(T); % se i punti del qqplot stanno su una retta allora i dati
pause      % sono approssimativamente gaussiani

figure(1)
clf
plot(T); % per vedere se i dati non seguono trend
xlabel('observation n.');
ylabel('temperature');
pause

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tsub=(reshape(T,5,25)); % raggruppo i dati in una matrice 5x25
tmean=mean(tsub); % faccio la media per colonne
range=max(tsub)-min(tsub); % faccio il range per colonne

%%%%%%%%%%%%%%%% CARTA \bar X %%%%%%%%%%%%%%%%%

samplen=1:25;
rbar=mean(range);
CLm=mean(tmean);
CL=[CLm CLm]; % centerline della carta \bar X
UCLm=CLm+0.577*rbar; % limite di controllo sup.
UCLm=[UCLm UCLm]; 
LCLm=CLm-0.577*rbar; % limite di controllo inf.
LCLm=[LCLm LCLm];

%%%%%%%%%%%%%%%% CARTA R %%%%%%%%%%%%%%%%%%%%%

CLr=rbar;
CLr=[CLr CLr]; % centerline della carta R
UCLr=2.114*rbar;    % limite di controllo sup.
UCLr=[UCLr UCLr];
LCLr=0;             % limite di controllo inf.
LCLr=[LCLr LCLr];



tt=[0 25];

figure(2)
clf

subplot(211)
plot(samplen,tmean,'k',samplen,tmean,'dk',tt,UCLm,'k',tt,LCLm,'k',tt,CL,'k--')
xlabel('sample n.'),ylabel('Mean of temperature'),title('Mean of temperature, LCL = 344.57, UCL = 355.41')

subplot(212)
plot(samplen,range,'dk-',tt,UCLr,'k',tt,LCLr,'k',tt,CLr,'k--')
xlabel('sample n.'),ylabel('Range of temperature'),title('Range of temperature, LCL = 0, UCL = 19.9')

