Mercurial > hg > ltpda

view m-toolbox/test/test_ao_psd_var.m @ 52:daf4eab1a51e database-connection-manager tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix. Default password should be [] not an empty string
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 07 Dec 2011 17:29:47 +0100
parents f0afece42f48
children
line wrap: on
line source


%%
a = ao(plist('waveform', 'noise', 'nsecs', 1e6, 'sigma', 25, 'yunits', 'm'));

Navs = [2:10 20:20:100 200:200:1000 5000 10000];

for jj=1:length(Navs)
  S = a.psd(plist('navs',Navs(jj),'win','BH92'));
  rat(jj) = mean((S.dy)./S.y);
  N(jj) = S.data.navs;
end

figure
loglog(Navs, rat, 'b+');
hold on
grid on
loglog(Navs, 1./sqrt(N),'k');
xlabel('N');
ylabel('(S.dy)./S.y');
legend('','1/sqrt(N)')