1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - INRIA - Serge Steer
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
a = [0.8604043 , 0.0070020 , 0.6223373 , -1.356213 , -4.2301775
0.159714 , 0.0857791 , -0.2367751 , 0.4958580 , 0.6398817
-4.3054931 , -0.0365878 , 2.1784911 , 0.0314793 , 2.3728994
-3.7115779 , -0.5370809 , 2.4291716 , -0.6352663 , 2.6972781
6.3580966 , 0.1377712 , -4.0461243 , -0.5192899 , -4.0394083];
b = [-0.0532544 , -0.1494083 , -0.0098619
-0.0355030 , -0.0162722 , 0.1045365
-0.1301775 , 0.1736686 , -0.0611440
0.1834320 , 0.1757396 , -0.5956607
0.1775148 , -0.1186391 , 0.1439842];
c = [2 , 7 , -2 , 5 , 1
0 , -1 , 3 , 0 , 2];
d = [1 , 0 , 0
0 , 0 , 0];
sl=syslin('c',a,b,c,d);
slc=contrss(sl,%eps);
if size(slc.A,1)<>5 then bugmes();quit;end
slc=contrss(sl);
n=size(slc.A,1);
if n<>4 then bugmes();quit;end
[n,U,ind,V,Ac,Bc]=contr(a,b);
if norm(Ac(1:n,1:n)-slc.A)>1d-14 then bugmes();quit;end
slo=obsvss(sl);
if size(slo.A,1)<>5 then bugmes();quit;end
slo=obsvss(sl,0.0001);
n=size(slo.A,1);
if n<>3 then bugmes();quit;end
[n,U,ind,V,Ao,Co]=contr(a',c',0.0001);
if norm(Ao(1:n,1:n)'-slo.A)>1d-14 then bugmes();quit;end
slm=minss(sl);
if size(slm.A,1)<>5 then bugmes();quit;end
slm=minss(sl,0.0001);
if size(slm.A,1)<>2 then bugmes();quit;end
|
