
function [M,CL,UCL,LCL] = uwma(T,w,n1)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% questa function calcola i punti da disegnare sulla
% carta UWMA, la centerline e i limiti di controllo
% partendo dal vettore di dati T, dal n. w
% di elementi per ogni gruppo e dalla lunghezza n1
% del vettore T
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

CL=mean(T); % centerline
SD=std(T);    % SD del vettore T

for i=1:n1
    if i<w  % quando i<w non ho abbastanza elem. per il gruppo
        M(i)=sum(T(1:i))/i; 
        UCL(i)=CL+3*SD/sqrt(i);
        LCL(i)=CL-3*SD/sqrt(i);
    else   % quando i>w posso fare i gruppi di dim. w
        M(i)=sum(T(i-w+1:i))/w;
        UCL(i)=CL+3*SD/sqrt(w);
        LCL(i)=CL-3*SD/sqrt(w);
    end
end

return