clc
close all
clear all

%%
N = 1e7;
Fs = 1000;
%%

data_white.rate = Fs;
data_white.freq = noise(N, 1, 'white');
[ad_white, S, err_white, tau] = allan(data_white, 2.^[0:nextpow2(length(data_white.freq))]./Fs, 'white noise',0);
[p_white, f_white] = pwelch(data_white.freq, [], [], [], data_white.rate, 'onesided');

data_pink.rate = Fs;
data_pink.freq = noise(N, 1, 'pink');
[ad_pink, S, err_pink, tau] = allan(data_pink, 2.^[0:nextpow2(length(data_pink.freq))]./Fs, 'pink noise',0);
[p_pink, f_pink] = pwelch(data_pink.freq, [], [], [], data_pink.rate, 'onesided');


data_brown.rate = Fs;
data_brown.freq = noise(N, 1, 'brown');
[ad_brown, S, err_brown, tau] = allan(data_brown, 2.^[0:nextpow2(length(data_brown.freq))]./Fs, 'brown noise',0);
[p_brown, f_brown] = pwelch(data_brown.freq, [], [], [], data_brown.rate, 'onesided');

figure
hold on
loglogerr(tau, ad_white, err_white, 'b')
loglogerr(tau, ad_pink, err_pink, 'm')
loglogerr(tau, ad_brown, err_brown, 'k')
grid on

figure
hold on
semilogx(f_white, 10*log10(p_white), 'b')
semilogx(f_pink, 10*log10(p_pink), 'm')
semilogx(f_brown, 10*log10(p_brown), 'k')
grid on