0
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
1 % TEST.M
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
2 %
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
3 %
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
4 % written by ao2m / $Id: test.m,v 1.1 2007/03/19 15:43:09 hewitson Exp $%
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
5 % based on analysis object:
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
6 % name: lpsd(x12ns)
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
7 % created: 19-Mar-2007 11:38:01
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
8 % original m-file:
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
9 %
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
10 %
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
11
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
12
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
13 a23 = ao(plist([param('filename', 'eta12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
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]) ]))) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
15 a21 = split(a22, plist([param('times', [1000 13000]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
16 a20 = ao(plist([param('filename', 'phi12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
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]) ]))) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
18 a18 = split(a19, plist([param('times', [1000 13000]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
19 a17 = ao(plist([param('filename', 'eta1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
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]) ]))) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
21 a15 = split(a16, plist([param('times', [1000 13000]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
22 a14 = ao(plist([param('filename', 'phi1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
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]) ]))) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
24 a12 = split(a13, plist([param('times', [1000 13000]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
25 a11 = ao(plist([param('filename', 'x12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
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]) ]))) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
27 a9 = split(a10, plist([param('times', [1000 13000]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
28 a8 = ltpda_timedomainfit(a9, a12, a15, a18, a21);
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
29 a7 = ao(plist([param('filename', 'eta12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
30 a6 = ao(plist([param('filename', 'phi12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
31 a5 = ao(plist([param('filename', 'eta1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
32 a4 = ao(plist([param('filename', 'phi1.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
33 a3 = ao(plist([param('filename', 'x12.txt') param('filepath', '/Users/hewitson/working/ltp/ltpda/software/m-toolbox/test/tdfit') ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
34 a2 = ltpda_lincom(a3, a4, a5, a6, a7, a8, plist([param('coeffs', [1 -0.21658 0.079907 -0.064982 -0.014581]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
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]) ]));
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
36
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
37
|
Daniele Nicolodi <nicolodi@science.unitn.it>
parents:
diff
changeset
|
38 % END
|