signal_processing plugged.
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / rpem.tst
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) 2011 - DIGITEO
3 //
4 // This file must be used under the terms of the CeCILL.
5 // This source file is licensed as described in the file COPYING, which
6 // you should have received as part of this distribution.  The terms
7 // are also available at    
8 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
9
10 test_path = SCI + '/modules/signal_processing/tests/unit_tests/';
11
12 n = 2;
13 theta = zeros(1,3*n);
14 p = eye(3*n,3*n);
15 phi = theta;
16 psi = theta;
17 l = theta;
18 w0 = list(theta, p, l, phi, psi);
19 u0 = [1 2 3 4];
20 y0 = [5 2 7 3];
21 lambda = [1 1 1];
22 kappa = [1 1 1];
23 c = 500;
24
25 import_from_hdf5(test_path + "rpem_ref1.h5");
26 w1 = rpem(w0, u0, y0);
27 if w1 <> w1_ref then pause, end
28 [w1,v] = rpem(w0, u0, y0);
29 if w1 <> w1_ref1 then pause, end
30 if v <> v_ref1 then pause, end
31
32 import_from_hdf5(test_path + "rpem_ref2.h5");
33 w1 = rpem(w0, u0, y0, lambda);
34 if w1 <> w1_ref2 then pause, end
35
36 [w1,v] = rpem(w0, u0, y0, lambda);
37 if w1 <> w1_ref2 then pause, end
38 if v <> v_ref2 then pause, end
39
40 import_from_hdf5(test_path + "rpem_ref3.h5");
41
42 w1 = rpem(w0, u0, y0, lambda, kappa);
43 if w1 <> w1_ref3 then pause, end
44
45 [w1,v] = rpem(w0, u0, y0, lambda, kappa);
46 if w1 <> w1_ref3 then pause, end
47 if v <> v_ref3 then pause, end
48
49 import_from_hdf5(test_path + "rpem_ref4.h5");
50 w1 = rpem(w0, u0, y0, lambda, kappa, c);
51 if w1 <> w1_ref4 then pause, end
52
53 [w1,v] = rpem(w0, u0, y0, lambda, kappa, c);
54 if w1 <> w1_ref4 then pause, end
55 if v <> v_ref4 then pause, end