Mercurial > hg > ltpda

diff 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
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/test/tdfit/test.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,38 @@
+% 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