5f5e4f802d9d94aafc62527c9df1bf7216235aba
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / syredi.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 [fact,b2,b1,b0,c1,c0,zzeros,zpoles]=syredi(1,4,[2*%pi/10,4*%pi/10,0,0],0.02,0.001);
13 ref_fact = 0.005976;
14 if norm(fact-ref_fact) > 0.0000037 then pause,end;
15 if or(b2<>[1 1 1]) then pause,end;
16 if or(size(b1)<>[1 3]) then pause,end;
17 if b1(3)<> 1 then pause,end;
18 if or(b0<>[1 1 0]) then pause,end;
19 if or(size(c1)<>[1 3]) then pause,end;
20 if or(size(c0)<>[1 3]) then pause,end;
21 if c0(3)<> 0 then pause,end;
22 if or(size(zzeros)<>[1 5]) then pause,end;
23 if or(size(zpoles)<>[1 5]) then pause,end;
24
25 import_from_hdf5(test_path+"syredi.h5");
26 [fact1, b21, b11, b01, c11, c01, zzeros1, zpoles1] = syredi(1, 4, [1, 2, 0, 0], 0.02, 0.001);
27 if fact1 <> fact1_ref then pause, end
28 if b21 <> b21_ref then pause, end
29 if b11 <> b11_ref then pause, end
30 if b01 <> b01_ref then pause, end
31 if c11 <> c11_ref then pause, end
32 if c01 <> c01_ref then pause, end
33 if zzeros1 <> zzeros1_ref then pause, end
34 if zpoles1 <> zpoles1_ref then pause, end
35
36 [fact2, b22, b12, b02, c12, c02, zzeros2, zpoles2] = syredi(3, 4, [1, 2, 3, 3.01], 0.1, 0.001);
37 if fact2 <> fact2_ref then pause, end
38 if b22 <> b22_ref then pause, end
39 if b12 <> b12_ref then pause, end
40 if b02 <> b02_ref then pause, end
41 if c12 <> c12_ref then pause, end
42 if c02 <> c02_ref then pause, end
43 if zzeros2 <> zzeros2_ref then pause, end
44 if zpoles2 <> zpoles2_ref then pause, end
45