Mercurial > hg > ltpda

comparison m-toolbox/test/test_ao_psd_var.m @ 0:f0afece42f48

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:f0afece42f48
1
2 %%
3 a = ao(plist('waveform', 'noise', 'nsecs', 1e6, 'sigma', 25, 'yunits', 'm'));
4
5 Navs = [2:10 20:20:100 200:200:1000 5000 10000];
6
7 for jj=1:length(Navs)
8 S = a.psd(plist('navs',Navs(jj),'win','BH92'));
9 rat(jj) = mean((S.dy)./S.y);
10 N(jj) = S.data.navs;
11 end
12
13 figure
14 loglog(Navs, rat, 'b+');
15 hold on
16 grid on
17 loglog(Navs, 1./sqrt(N),'k');
18 xlabel('N');
19 ylabel('(S.dy)./S.y');
20 legend('','1/sqrt(N)')