Mercurial > hg > ltpda

comparison 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
comparison
equal deleted inserted replaced
-1:000000000000 0:f0afece42f48
1 % Test the time-domain conversion to acceleration
2 %
3 % M Hewitson 11-04-08
4 %
5 % $Id: test_ltpda_mdc1_x2acc.m,v 1.1 2008/04/14 15:39:44 hewitson Exp $
6 %
7
8 % Make a white-noise AO
9 fs = 10;
10 nsecs = 100000;
11
12 a = ao(plist('tsfcn', 'randn(size(t))', 'fs', fs, 'nsecs', nsecs));
13
14 % Apply filter
15 b = ltpda_mdc1_x2acc(a, plist('Omega2', -1.3e-6));
16
17 % Split off transients
18 a = split(a, plist('split_type', 'times', 'times', [100 nsecs-100]));
19
20 % Make spectrum
21 bxx = ltpda_lpsd(b);
22
23 iplot(bxx);
24
25 % Make TF
26 tf = ltpda_tfe(a,b);
27
28 iplot(tf(1,2))
29
30
31 % END