% TEST.M 
% 
% 
% written by ao2m / $Id: test.m,v 1.1 2007/03/19 15:43:09 hewitson Exp $% 
% based on analysis object:
%            name: lpsd(x12ns)
%         created: 19-Mar-2007 11:38:01
% original m-file: 
% 
% 
 
 
a23 = ao(plist([param('filename', 'eta12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a22 = filter(a23, plist([param('filter', miir(plist([param('type', 'bandpass') param('gain', [1]) param('fs', [32.46]) param('order', [2]) param('fc', [0.003        0.03]) param('ripple', [0.5]) ]))) ]));
a21 = split(a22, plist([param('times', [1000  13000]) ]));
a20 = ao(plist([param('filename', 'phi12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a19 = filter(a20, plist([param('filter', miir(plist([param('type', 'bandpass') param('gain', [1]) param('fs', [32.46]) param('order', [2]) param('fc', [0.003        0.03]) param('ripple', [0.5]) ]))) ]));
a18 = split(a19, plist([param('times', [1000  13000]) ]));
a17 = ao(plist([param('filename', 'eta1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a16 = filter(a17, plist([param('filter', miir(plist([param('type', 'bandpass') param('gain', [1]) param('fs', [32.46]) param('order', [2]) param('fc', [0.003        0.03]) param('ripple', [0.5]) ]))) ]));
a15 = split(a16, plist([param('times', [1000  13000]) ]));
a14 = ao(plist([param('filename', 'phi1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a13 = filter(a14, plist([param('filter', miir(plist([param('type', 'bandpass') param('gain', [1]) param('fs', [32.46]) param('order', [2]) param('fc', [0.003        0.03]) param('ripple', [0.5]) ]))) ]));
a12 = split(a13, plist([param('times', [1000  13000]) ]));
a11 = ao(plist([param('filename', 'x12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a10 = filter(a11, plist([param('filter', miir(plist([param('type', 'bandpass') param('gain', [1]) param('fs', [32.46]) param('order', [2]) param('fc', [0.003        0.03]) param('ripple', [0.5]) ]))) ]));
a9 = split(a10, plist([param('times', [1000  13000]) ]));
a8 = ltpda_timedomainfit(a9, a12, a15, a18, a21);
a7 = ao(plist([param('filename', 'eta12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a6 = ao(plist([param('filename', 'phi12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a5 = ao(plist([param('filename', 'eta1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a4 = ao(plist([param('filename', 'phi1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a3 = ao(plist([param('filename', 'x12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
a2 = ltpda_lincom(a3, a4, a5, a6, a7, a8, plist([param('coeffs', [1    -0.21658    0.079907   -0.064982   -0.014581]) ]));
a1 = ltpda_lpsd(a2, plist([param('Kdes', [100]) param('Kmin', [1]) param('Jdes', [1000]) param('Win', specwin('Kaiser', 1000, 250)) param('Nfft', [1000]) param('Olap', [0.79592]) ]));
 
 
% END
