License Header change: Removed the LICENSE_END before beta
[scilab.git] / scilab / modules / cacsd / macros / cont_mat.sci
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) INRIA -
3 //
4 // Copyright (C) 2012 - 2016 - Scilab Enterprises
5 //
6 // This file is hereby licensed under the terms of the GNU GPL v2.0,
7 // pursuant to article 5.3.4 of the CeCILL v.2.1.
8 // This file was originally licensed under the terms of the CeCILL v2.1,
9 // and continues to be available under such terms.
10 // For more information, see the COPYING file which you should have received
11 // along with this program.
12
13 function c=cont_mat(a,b)
14     //c=cont_mat(a,b) or c=cont_mat(sl) is the controllability matrix.
15     // of the pair a,b or of the system sl=[a,b,c,d] (syslin list)
16     //                 2       n
17     //i.e. c=[b, ab, ab,...; ab ]
18     //!
19     [lhs,rhs]=argn(0)
20     select typeof(a)
21     case "constant"  then
22         if rhs==1 then
23             error(msprintf(gettext("%s: Wrong number of input arguments: %d expected"),"cont_mat",2)),
24         end
25         [m,n]=size(a)
26         if m<>n then error(20,1),end
27         [mb,nb]=size(b);
28         if mb<>n then error(60),end
29
30     case "state-space" then
31         if rhs==2 then
32             error(msprintf(gettext("%s: Wrong number of input arguments: %d expected"),"cont_mat",1)),
33         end
34         [a,b]=a([2,3])
35         [n,n]=size(a)
36     else
37         if rhs==1 then
38             error(msprintf(gettext("%s: Wrong type for input argument #%d: Linear dynamical system expected.\n"),"cont_mat",1))
39         else
40             error(msprintf(gettext("%s: Wrong type of input argument #%d: Array of floating point numbers expected.\n"),"cont_mat",1))
41         end
42
43     end;
44     c=b;for k=1:n-1, c=[b,a*c],end
45 endfunction