Mercurial > hg > ltpda
view m-toolbox/test/test_ltpda_xcorr.m @ 12:86aabb42dd84 database-connection-manager
Use utils.repository utilities
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% TEST_LTPDA_XCORR tests the cross-correlation function ltpda_xcorr. % % M Hewitson 19-06-07 % % $Id: test_ltpda_xcorr.m,v 1.13 2011/04/09 16:59:57 hewitson Exp $ % function test_ltpda_xcorr() % Make test AOs nsecs = 10; fs = 1000; pl = plist(); pl.append('nsecs', nsecs); pl.append('fs', fs); pl.append('tsfcn', 'sin(2*pi*7.433*t) + randn(size(t))'); a1 = ao(pl); % Filter one time-series % Make a filter f1 = miir(plist('type', 'highpass', 'fc', 20, 'fs', fs)); % Filter the input data [a2, f1] = filter(a1,plist(param('filter', f1))); % Make X-correlation from a1 to a2 a4 = xcorr(a1,a2, plist(param('Scale', 'coeff'))); % Plot results and history iplot(a4); plot(a4.hist); % Reproduce from history a_out = rebuild(a4); iplot(a_out) plot(a_out.hist) end