a211046578f60be2cfc624d358b0cfce11710f93
[scilab.git] / scilab / modules / cacsd / tests / nonreg_tests / bug_7794.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2016 - Scilab Enterprises - Adeline CARNIS
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 //<-- CLI SHELL MODE -->
9 // <-- Non-regression test for bug 7794 -->
10 //
11 // <-- Bugzilla URL -->
12 // http://bugzilla.scilab.org/show_bug.cgi?id=7794
13 //
14 // <-- Short Description -->
15 //  Example of findABCD help page failed.
16
17 A = [ 0.5, 0.1,-0.1, 0.2;
18       0.1, 0,  -0.1,-0.1;      
19      -0.4,-0.6,-0.7,-0.1;  
20       0.8, 0,  -0.6,-0.6];      
21 B = [0.8;0.1;1;-1];
22 C = [1 2 -1 0];
23 SYS=syslin(0.1,A,B,C);
24 nsmp=100;
25 U=prbs_a(nsmp,nsmp/5);
26 Y=(flts(U,SYS)+0.3*rand(1,nsmp,'normal'));
27
28 // Compute R
29 S=15;
30 [R,N1,SVAL] = findR(S,Y',U');
31 N=3;
32 SYS1 = findABCD(S,N,1,R) ;SYS1.dt=0.1;
33
34 SYS1.X0 = inistate(SYS1,Y',U');
35 expectedX0 = [0.180418869; 0.18003251391; -0.06993989023];
36 assert_checkalmostequal(SYS1.X0, expectedX0, 1.e-8);
37
38