Mercurial > hg > ltpda

comparison m-toolbox/test/tdfit/test.m @ 0:f0afece42f48

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:f0afece42f48
1 % TEST.M
2 %
3 %
4 % written by ao2m / $Id: test.m,v 1.1 2007/03/19 15:43:09 hewitson Exp $%
5 % based on analysis object:
6 % name: lpsd(x12ns)
7 % created: 19-Mar-2007 11:38:01
8 % original m-file:
9 %
10 %
11
12
13 a23 = ao(plist([param('filename', 'eta12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
14 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]) ]))) ]));
15 a21 = split(a22, plist([param('times', [1000 13000]) ]));
16 a20 = ao(plist([param('filename', 'phi12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
17 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]) ]))) ]));
18 a18 = split(a19, plist([param('times', [1000 13000]) ]));
19 a17 = ao(plist([param('filename', 'eta1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
20 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]) ]))) ]));
21 a15 = split(a16, plist([param('times', [1000 13000]) ]));
22 a14 = ao(plist([param('filename', 'phi1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
23 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]) ]))) ]));
24 a12 = split(a13, plist([param('times', [1000 13000]) ]));
25 a11 = ao(plist([param('filename', 'x12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
26 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]) ]))) ]));
27 a9 = split(a10, plist([param('times', [1000 13000]) ]));
28 a8 = ltpda_timedomainfit(a9, a12, a15, a18, a21);
29 a7 = ao(plist([param('filename', 'eta12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
30 a6 = ao(plist([param('filename', 'phi12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
31 a5 = ao(plist([param('filename', 'eta1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
32 a4 = ao(plist([param('filename', 'phi1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
33 a3 = ao(plist([param('filename', 'x12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
34 a2 = ltpda_lincom(a3, a4, a5, a6, a7, a8, plist([param('coeffs', [1 -0.21658 0.079907 -0.064982 -0.014581]) ]));
35 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]) ]));
36
37
38 % END