* Bugs 16090 16163 fixed: prbs_a(n,nc) infringed nc & slow. Page overhauled
[scilab.git] / scilab / modules / cacsd / tests / nonreg_tests / bug_13581.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2017 - 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 // <-- NO CHECK REF -->
10
11 // <-- Non-regression test for bug 13581 -->
12 //
13 // <-- Bugzilla URL -->
14 // http://bugzilla.scilab.org/13581
15 //
16 // <-- Short Description -->
17 // The example of inistate help page did not work on Windows 64bit
18
19 //generate data from a given linear system
20 A = [ 0.5, 0.1,-0.1, 0.2;
21 0.1, 0,  -0.1,-0.1;
22 -0.4,-0.6,-0.7,-0.1;
23 0.8, 0,  -0.6,-0.6];
24 B = [0.8;0.1;1;-1];
25 C = [1 2 -1 0];
26 SYS=syslin(0.1,A,B,C);
27 nsmp=100;
28 U=prbs_a(nsmp,nsmp/5);
29 Y=(flts(U,SYS)+0.3*rand(1,nsmp,"normal"));
30 // Compute R
31 S=15;
32 [R,N1,SVAL] = findR(S,Y',U');
33 N=3;
34 SYS1 = findABCD(S,N,1,R) ;
35 SYS1.dt=0.1;
36 x0 = inistate(SYS1,Y',U');
37 expected = [1.358882940201 ; -0.171852413543 ; 1.076698107654];
38 assert_checkalmostequal(x0, expected);