3890c5f22751fa998a9d7b8eeec11973fc6eb890
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / rpem.tst
1 // -*- scilab -*-
2
3 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
4 // Copyright (C) 2009 - Digiteo - Vincent Liard
5 //
6 // This file must be used under the terms of the CeCILL.
7 // This source file is licensed as described in the file COPYING, which
8 // you should have received as part of this distribution.  The terms
9 // are also available at    
10 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
11
12 // 'rpem' unit tests
13
14 test_path = pathconvert(SCI + '/modules/signal_processing/tests/unit_tests');
15
16 // requires test_utilities.sci
17 exec(test_path + 'test_utilities.sci');
18
19 n = 2;
20 theta = zeros(1,3*n);
21 p = eye(3*n,3*n);
22 phi = theta;
23 psi = theta;
24 l = theta;
25 w0 = list(theta, p, l, phi, psi);
26 u0 = [1 2 3 4];
27 y0 = [5 2 7 3];
28 lambda = [1 1 1];
29 kappa = [1 1 1];
30 c = 500;
31
32 // test 1 ****************************************************************
33 load(test_path + 'rpem_ref1.dat', 'w1_ref1', 'v_ref1');
34
35 check(rpem(w0, u0, y0), w1_ref1);
36
37 [w1,v] = rpem(w0, u0, y0);
38 check(w1, w1_ref1);
39 check(v, v_ref1);
40 // test 2 ****************************************************************
41 load(test_path + 'rpem_ref2.dat', 'w1_ref2', 'v_ref2');
42
43 check(rpem(w0, u0, y0, lambda), w1_ref2);
44
45 [w1,v] = rpem(w0, u0, y0, lambda);
46 check(w1, w1_ref2);
47 check(v, v_ref2);
48 // test 3 ****************************************************************
49 load(test_path + 'rpem_ref3.dat', 'w1_ref3', 'v_ref3');
50
51 check(rpem(w0, u0, y0, lambda, kappa), w1_ref3);
52
53 [w1,v] = rpem(w0, u0, y0, lambda, kappa);
54 check(w1, w1_ref3);
55 check(v, v_ref3);
56 // test 4 ****************************************************************
57 load(test_path + 'rpem_ref4.dat', 'w1_ref4', 'v_ref4');
58
59 check(rpem(w0, u0, y0, lambda, kappa, c), w1_ref4);
60
61 [w1,v] = rpem(w0, u0, y0, lambda, kappa, c);
62 check(w1, w1_ref4);
63 check(v, v_ref4);
64 //****************************************************************