
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% questo file carica un file di dati di cui si vuole 
% fare la carta di controllo UWMA per diversi valori w
% del n. di elementi del gruppo. Si richiama la function
% uwma.m che calcola i punti da disegnare, la centerline
% e i limiti di controllo
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


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


w1=1; % prendo 4 diversi valori per w
w2=4;
w3=8;
w4=16;

%%% richiamo 4 volte la function uwma.m

[M1,CL,UCL1,LCL1]=uwma(T,w1,n);
[M2,CL,UCL2,LCL2]=uwma(T,w2,n);
[M3,CL,UCL3,LCL3]=uwma(T,w3,n);
[M4,CL,UCL4,LCL4]=uwma(T,w4,n);

% disegno le 4 carte di controllo UWMA

nsam=1:n;
t1=[0,125];
CL=[CL,CL];

figure(1)
% clf
subplot(411); plot(nsam,M1,'r',t1,CL,'g'); axis([0 125 335 365]);
hold
stairs(nsam,UCL1);
stairs(nsam,LCL1);
title(['UWMA(',int2str(w1),')']);
subplot(412); plot(nsam,M2,'r',t1,CL,'g'); axis([0 125 335 365]);
hold
stairs(nsam,UCL2);
stairs(nsam,LCL2);
title(['UWMA(',int2str(w2),')']);
subplot(413); plot(nsam,M3,'r',t1,CL,'g'); axis([0 125 335 365]);
hold
stairs(nsam,UCL3);
stairs(nsam,LCL3);
title(['UWMA(',int2str(w3),')']);
subplot(414); plot(nsam,M4,'r',t1,CL,'g'); axis([0 125 335 365]);
hold
stairs(nsam,UCL4);
stairs(nsam,LCL4);
title(['UWMA(',int2str(w4),')']);
xlabel('sample n.');