--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/utils/test_blwhitenoise.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,72 @@
+% A test script for utils.math.blwhitenoise
+% 
+% L. Ferraioli 18-11-08
+% 
+% $Id: test_eigcsd.m,v 1.2 2009/06/10 16:12:06 luigi Exp $
+% 
+
+%%
+
+clear all
+
+%%
+
+% sampling frequency
+fs = 10;
+% higher bandwhidth frequency
+fh = fs/2;
+% lower bandwidth frequency
+fl = 1e-3;
+% number of data points
+npts = 1e6;
+
+% noise generation
+wn = utils.math.blwhitenoise(npts,fs,fl,fh);
+
+wnao = ao(tsdata(wn,fs));
+% wnao.setNane('BLWN');
+
+%% plot and statistics
+
+
+iplot(wnao)
+
+% Check histogram
+wnao_hist = hist(wnao,plist('N',500,'Norm','true'));
+wnao_norm = normdist(wnao,plist('N',500,'Norm','true'));
+
+plpl = plist( ...
+  'Yscales', {'All','log'}, ...
+  'Legends', {'Histogram', 'Gaussian'});
+
+iplot(wnao_hist,wnao_norm,plpl)
+
+%% fft
+
+WNf = fft(wnao);
+iplot(WNf)
+
+%% psd
+
+% expected white noise at fs
+f = [logspace(-5,fs/2,300)]';
+expec = ao(plist('fsfcn', 'ones(size(f))', 'f', f));
+expec = expec.*(2/fs);
+WNxx = psd(wnao,plist('order',-1,'nfft',1e4));
+iplot(WNxx,expec)
+
+%% compare randn
+
+aa = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', 1e5));
+
+iplot(wnao)
+
+% Check histogram
+aa_hist = hist(aa,plist('N',500,'Norm','true'));
+aa_norm = normdist(aa,plist('N',500,'Norm','true'));
+
+plpl = plist( ...
+  'Yscales', {'All','log'}, ...
+  'Legends', {'Histogram', 'Gaussian'});
+
+iplot(aa_hist,aa_norm,plpl)