
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% questo file carica un file di dati di cui si vuole 
% fare la carta di controllo EWMA per diversi valori di 
% lambda. Si richiama la function ewma.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


lambda1=1; % prendo 4 diversi valori per lambda
lambda2=0.9;
lambda3=0.5;
lambda4=0.1;

L=3;

%%% richiamo 4 volte la function ewma.m

[Z1,CL,UCL1,LCL1]=ewma(T,lambda1,L,n);
[Z2,CL,UCL2,LCL2]=ewma(T,lambda2,L,n);
[Z3,CL,UCL3,LCL3]=ewma(T,lambda3,L,n);
[Z4,CL,UCL4,LCL4]=ewma(T,lambda4,L,n);

% disegno le 4 carte di controllo EWMA

nsam=1:n;
t1=[0,125];
CL=[CL,CL];
size(UCL1)

figure(1)
% clf
subplot(411); plot(nsam,Z1,'r',t1,CL,'g'); %axis([0 125 335 365]);
hold
stairs(nsam,UCL1);
stairs(nsam,LCL1);
title(['EWMA(',num2str(lambda1),')']);
subplot(412); plot(nsam,Z2,'r',t1,CL,'g'); %axis([0 125 335 365]);
hold
stairs(nsam,UCL2);
stairs(nsam,LCL2);
title(['EWMA(',num2str(lambda2),')']);
subplot(413); plot(nsam,Z3,'r',t1,CL,'g'); %axis([0 125 335 365]);
hold
stairs(nsam,UCL3);
stairs(nsam,LCL3);
title(['EWMA(',num2str(lambda3),')']);
subplot(414); plot(nsam,Z4,'r',t1,CL,'g'); %axis([0 125 335 365]);
hold
stairs(nsam,UCL4);
stairs(nsam,LCL4);
title(['EWMA(',num2str(lambda4),')']);
xlabel('sample n.');