Mercurial > hg > ltpda

diff 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
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/test_ao_psd_var.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,20 @@
+
+%%
+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)')