Mercurial > hg > ltpda
view m-toolbox/test/utils/test_blwhitenoise.m @ 37:a4b7ceae0403 database-connection-manager
Show backtrace on unit test errors
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% 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)