signal_processing plugged.
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / rpem.dia.ref
1
2 // -*- scilab -*-
3
4 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
5
6 // Copyright (C) 2009 - Digiteo - Vincent Liard
7
8 //
9
10 // This file must be used under the terms of the CeCILL.
11
12 // This source file is licensed as described in the file COPYING, which
13
14 // you should have received as part of this distribution.  The terms
15
16 // are also available at    
17
18 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
19
20 // 'rpem' unit tests
21
22 test_path = pathconvert(SCI + '/modules/signal_processing/tests/unit_tests');
23
24 // requires test_utilities.sci
25
26 exec(test_path + 'test_utilities.sci');
27
28 n = 2;
29
30 theta = zeros(1,3*n);
31
32 p = eye(3*n,3*n);
33
34 phi = theta;
35
36 psi = theta;
37
38 l = theta;
39
40 w0 = list(theta, p, l, phi, psi);
41
42 u0 = [1 2 3 4];
43
44 y0 = [5 2 7 3];
45
46 lambda = [1 1 1];
47
48 kappa = [1 1 1];
49
50 c = 500;
51
52 // test 1 ****************************************************************
53
54 //load(test_path + 'rpem_ref1.dat', 'w1_ref1', 'v_ref1');
55
56 import_from_hdf5(test_path+"rpem_ref1.h5");
57
58 check(rpem(w0, u0, y0), w1_ref1);
59
60 [w1,v] = rpem(w0, u0, y0);
61
62 check(w1, w1_ref1);
63
64 check(v, v_ref1);
65
66 // test 2 ****************************************************************
67
68 //load(test_path + 'rpem_ref2.dat', 'w1_ref2', 'v_ref2');
69
70 import_from_hdf5(test_path+"rpem_ref2.h5");
71
72 check(rpem(w0, u0, y0, lambda), w1_ref2);
73
74 [w1,v] = rpem(w0, u0, y0, lambda);
75
76 check(w1, w1_ref2);
77
78 check(v, v_ref2);
79
80 // test 3 ****************************************************************
81
82 //load(test_path + 'rpem_ref3.dat', 'w1_ref3', 'v_ref3');
83
84 import_from_hdf5(test_path+"rpem_ref3.h5");
85
86 check(rpem(w0, u0, y0, lambda, kappa), w1_ref3);
87
88 [w1,v] = rpem(w0, u0, y0, lambda, kappa);
89
90 check(w1, w1_ref3);
91
92 check(v, v_ref3);
93
94 // test 4 ****************************************************************
95
96 //load(test_path + 'rpem_ref4.dat', 'w1_ref4', 'v_ref4');
97
98 import_from_hdf5(test_path+"rpem_ref4.h5");
99
100 check(rpem(w0, u0, y0, lambda, kappa, c), w1_ref4);
101
102 [w1,v] = rpem(w0, u0, y0, lambda, kappa, c);
103
104 check(w1, w1_ref4);
105
106 check(v, v_ref4);
107
108 //****************************************************************
109