view m-toolbox/test/test_ltpda_ltfe.m @ 23:a71a40911c27
database-connection-manager
Update check for repository connection parameter in constructors
author |
Daniele Nicolodi <nicolodi@science.unitn.it> |
date |
Mon, 05 Dec 2011 16:20:06 +0100 (2011-12-05) |
parents |
f0afece42f48 |
children |
|
line source
% TEST_LTPDA_LTFE test the ltpde_ltfe function.
%
% M Hewitson 23-05-07
%
% $Id: test_ltpda_ltfe.m,v 1.19 2009/04/07 17:45:22 ingo Exp $
%
function test_ltpda_ltfe()
% Make test AO
nsecs = 10000;
fs = 10;
pl = plist('nsecs', nsecs, 'fs', fs, 'tsfcn', 'randn(size(t))');
a1 = ao(pl);
% Filter one time-series
% Make a filter
f1 = miir(plist('type', 'highpass', 'fc', 0.1));
% Filter the input data
[a2, f1] = filter(a1,plist('filter', f1));
% Make a filter
f1 = miir(plist('type', 'lowpass', 'fc', .03));
% Filter the input data
[a3, f1] = filter(a1,plist('filter', f1));
% Split data
pl = plist('times', [10 a1.data.nsecs]);
aoin = split(a1, pl);
aoout = split(a2, pl);
aoout2 = split(a3, pl);
% Make TF from a1 to a2
pl = plist('Kdes', 1000, 'Kmin', 2, 'Jdes', 1000, 'Order', 1);
% log TF estiamtes
tfs = ltfe(aoin, aoout, pl);
% Plot
iplot(tfs)
% Recreate from History
a_out = rebuild(tfs);
iplot(a_out)
b = a_out;
plot(b.hist)
end
% END