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

load circuiti2.txt
x=circuiti2;
n=length(x);
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  % centerline

UCL=cbar+3*sqrt(cbar); % limiti di controllo
LCL=cbar-3*sqrt(cbar);

t=[0,26];
t1=1:26;
lsup=[UCL,UCL];
linf=[LCL,LCL];
cl=[cbar,cbar];

plot(t1,x,'-o',t,lsup,'b',t,linf,'b',t,cl,'b');
xlabel('N. campione');
ylabel('N. di non conformità');