Mercurial > hg > ltpda

diff m-toolbox/test/MDC1/test_ltpda_mdc1_x2acc.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/MDC1/test_ltpda_mdc1_x2acc.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,31 @@
+% Test the time-domain conversion to acceleration
+% 
+% M Hewitson 11-04-08
+% 
+% $Id: test_ltpda_mdc1_x2acc.m,v 1.1 2008/04/14 15:39:44 hewitson Exp $
+% 
+
+% Make a white-noise AO
+fs    = 10;
+nsecs = 100000;
+
+a = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs));
+
+% Apply filter
+b = ltpda_mdc1_x2acc(a, plist('Omega2', -1.3e-6));
+
+% Split off transients
+a = split(a, plist('split_type', 'times', 'times', [100 nsecs-100]));
+
+% Make spectrum
+bxx = ltpda_lpsd(b);
+
+iplot(bxx);
+
+% Make TF
+tf = ltpda_tfe(a,b);
+
+iplot(tf(1,2))
+
+
+% END