function [CPss,f] = testBartlett_wn(s)
%
% test di Bartlett (o del periodogramma cumulato): serve per verificare se una sequenza è rumore bianco
%
% Una sequenza x(k) generata da rumore bianco presenta una sequenza di autocorrelazione in cui tutti i
% campioni per shift > 0 sono nulli. Quindi, osservare valori della sequenza di auto-correlazione
% significativamente (in senso statistico) non-nulli corrisponde alla presenza di componenti
% non-random nella sequenza x(k). Per individuare componenti periodiche mischiate al rumore è 
% vantaggioso considerare il "periodogramma cumulato":
%
N = length(s);
%
[Pss, f] = psd_Welch(s);
est_mean_s = mean(s);
est_var_s = sum(s.^2 - est_mean_s)/(N-1);
CPss = zeros(length(f), 1);
for i=1:length(f)
  CPss(i) = sum(Pss(1:i)) / (N * est_var_s);
end;
