view m-toolbox/test/MDC1/test_ltpda_mdc1_x2acc.m @ 0:f0afece42f48
Import.
| author |
Daniele Nicolodi <nicolodi@science.unitn.it> |
| date |
Wed, 23 Nov 2011 19:22:13 +0100 (2011-11-23) |
| parents |
|
| children |
|
line source
% 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
