Mercurial > hg > ltpda

comparison m-toolbox/test/test_median.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_MEDIAN test median() operator for analysis objects.
2 %
3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 %
5 % DESCRIPTION: TEST_MEDIAN Test median() operator for analysis objects with
6 % tsdata, fsdata cdata and xydata
7 %
8 % CALL: test_median;
9 %
10 % HISTORY: 31-05-2007 Diepholz
11 % Creation.
12 %
13 % VERSION: $Id: test_median.m,v 1.5 2009/08/10 13:37:58 hewitson Exp $
14 %
15 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
16 function test_median()
17
18 % Create tsdata with a sinewave
19 pl = plist('fs', 100, 'nsecs', 5, 'waveform', 'sine wave', 'f', 1.23, 'phi', 30);
20
21 ts_data = ao(pl);
22
23 % Create cdata with a 3x3 Matrix
24 matr = [1 2 3; 4 5 6; 7 8 9];
25
26 c_data = cdata(matr);
27 c_data = ao(c_data);
28
29 % Create xydata with a exp function
30 x = 1:1000;
31 y = exp(x.^.5/2);
32
33 xy_data = xydata(x,y);
34 xy_data = ao(xy_data);
35
36 % Some tests
37 ao_out = median(ts_data);
38 ao_out = median(ts_data, ts_data, ts_data);
39 ao_out = median(ts_data, xy_data);
40
41 ao_out = median(ts_data, xy_data);
42 ao_out = median(ts_data, xy_data);
43
44
45 ao_out = median(c_data, plist('axis', 'y'));
46
47 % Create different parameter lists to control the median test
48 pl1 = plist('xdata', 't', 'ydata', 'x');
49 ao_out = median(ts_data, pl1);
50
51 pl2 = plist('xdata', 'x', 'ydata', 'y');
52 ao_out = median(xy_data, pl2);
53
54 pl3 = plist('xdata', 'x');
55 ao_out = median(xy_data, pl3);
56
57 pl4 = plist('xdata', 'x');
58 ao_out = median(xy_data, pl4);
59
60 pl5 = plist('ydata', 'y', 'xdata', 'x');
61 ao_out = median(xy_data, pl5);
62
63 end