Mercurial > hg > ltpda

comparison m-toolbox/test/make_ltpda_objects.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 % make some test objects
2
3
4
5 fs = 10;
6
7 pl1 = plist();
8 pl1 = append(pl1, param('fs', fs));
9 pl1 = append(pl1, param('nsecs', 1));
10 pl1 = append(pl1, param('waveform', 'sine wave'));
11 pl1 = append(pl1, param('f', 1.23));
12 pl1 = append(pl1, param('phi', 30));
13
14 % tsdata
15 a1 = ao(pl1);
16
17 % fsdata
18 f = 1:10;
19 xx = randn(size(f));
20 a2 = ao(fsdata(f,xx,fs));
21
22 % xydata
23 x = 1:10;
24 y = exp(x.^.5/2);
25
26 a3 = ao(xydata(x,y));
27
28 % cdata
29 a4 = ao([1 -2 3; 4 5 -6; -7 8 9]);
30
31 a5 = ao(plist('tsfcn', 'randn(size(t))', 'nsecs', 1, 'fs', fs));
32 a6 = ao(plist('tsfcn', 'randn(size(t))', 'nsecs', 1, 'fs', fs));
33
34
35 % specwin
36 s = specwin('Hanning', 100);
37
38 % pzmodel
39 pzm = pzmodel(1, pz(1), pz(10));
40
41 % mfir
42 firf = mfir(plist('type', 'highpass', 'gain', 2, 'fs', 100, 'fc', 1));
43
44 % miir
45 iirf = miir(plist('type', 'highpass', 'gain', 2, 'fs', 100, 'fc', 1));
46
47 % time
48 t = time('2007-06-06 14:00:05.000');
49
50 % timespan
51 ts = timespan('14:00:00.000', '14:00:05.000');