%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% file per costruire la carta di controllo u
% i limiti di controllo sono ottenuti tramite
% l'approssimazione della Poisson con la normale
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

load computer.txt
x=computer;
n=length(x);
r=5;
u=computer/r;
alpha=0.05; % livello di significativita' per il
            % test di dispersione

%%%%%%% verifica della poissonianita' %%%%%%%%%%
%%%%%%% test di dispersione  %%%%%%%%%%%%%%%%%%%

csi=(n-1)*var(x)/mean(x) % devo confrontare questo
                         % valore con quelli delle
                         % tavole
liminf=chi2inv(alpha,n-1) % calcolo i limiti per la
limsup=chi2inv(1-alpha,n-1) % verifica della poissonianita'
                            % con livello di significativita'
                            % alpha e n-1 gradi di liberta'
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


cbar=sum(x)/n;
ubar=cbar/r

UCL=ubar+3*sqrt(ubar/r)
LCL=ubar-3*sqrt(ubar/r)

t=[0,20];
t1=1:20;
lsup=[UCL,UCL];
linf=[LCL,LCL];
cl=[ubar,ubar];

plot(t1,u,'-o',t,lsup,'b',t,linf,'b',t,cl,'b');
xlabel('N. campione');
ylabel('Non conformità per unità');