view src/ltpda_polyreg/test_ltpda_polyreg.m @ 36:5eb86f6881ef
database-connection-manager
Remove commented-out code
author |
Daniele Nicolodi <nicolodi@science.unitn.it> |
date |
Mon, 05 Dec 2011 16:20:06 +0100 (2011-12-05) |
parents |
f0afece42f48 |
children |
|
line source
% function test_ltpda_polydetrend()
% A test script for the AO implementation of detrending.
%
% M Hewitson 02-02-07
%
% $Id: test_ltpda_polyreg.m,v 1.2 2009/12/17 08:05:13 mauro Exp $
%
clear all;
for N = -1:3
Nsecs = 100000;
%% Make test AOs
pl1 = plist('waveform', 'sine wave', 'f', 0.1/Nsecs, 'nsecs', Nsecs, 'fs', 10);
pl2 = plist('waveform', 'noise', 'nsecs', Nsecs, 'fs', 10);
a = ao(pl1) + 0.01.*ao(pl2);
iplot(a)
%% Detrend with MATLAB
pl = plist('N', N);
tic
c = detrend(a, pl);
toc
c.setName('matlab');
%% Detrend with mex
tic
[y, a] = ltpda_polyreg(a.data.y, N);
toc
ts = c.data;
ts.setY(y);
d = ao(ts);
d.setName('mex');
%% Plot
iplot(c,d);
iplot(c./d)
end