%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% file per la costruzione di una carta di Shewart
% individuale. Si costruiscono la carta X e la
% carta MR (moving range), dopo aver controllato
% che i dati sono gaussiani e indipendenti
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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

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

%%%%%%%%%%%%%% CARTA MR %%%%%%%%%%%%%%%%%%%%%

samplen1=1:124;
n=length(T);
for i=1:n-1
    MR(i)=abs(T(i+1)-T(i));
end
MRbar=mean(MR)
sigMR=std(MR);
CLMR=[MRbar,MRbar];
UCLMR=3.267*MRbar
UCLMR=[UCLMR,UCLMR];
LCLMR=0
LCLMR=[LCLMR,LCLMR];


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

samplen=1:125;
CLm=mean(T)
CL=[CLm CLm]; % centerline della carta \bar X
UCLm=CLm+3*MRbar/1.128 % oppure UCLm=CLm+3*std(T)
UCLm=[UCLm UCLm]; 
LCLm=CLm-3*MRbar/1.128 % oppure LCLm=CLm-3*std(T)
LCLm=[LCLm LCLm];

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


tt=[0 125];

figure(2)
clf

subplot(211)
plot(samplen,T,'dk-',tt,UCLm,'k',tt,LCLm,'k',tt,CL,'k--')
xlabel('sample n.'),ylabel('Temperature');title('X chart')

subplot(212)
plot(samplen1,MR,'dk-',tt,UCLMR,'k',tt,LCLMR,'k',tt,CLMR,'k--')
xlabel('sample n.'),ylabel('Moving range');title('MR chart')

