diff options
Diffstat (limited to 'modules/cacsd/macros/minreal.sci')
| -rwxr-xr-x | modules/cacsd/macros/minreal.sci | 67 |
1 files changed, 67 insertions, 0 deletions
diff --git a/modules/cacsd/macros/minreal.sci b/modules/cacsd/macros/minreal.sci new file mode 100755 index 000000000..c8e15aead --- /dev/null +++ b/modules/cacsd/macros/minreal.sci @@ -0,0 +1,67 @@ +// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab +// Copyright (C) INRIA - +// +// This file must be used under the terms of the CeCILL. +// This source file is licensed as described in the file COPYING, which +// you should have received as part of this distribution. The terms +// are also available at +// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt + +function [a,b,c]=minreal(a,b,c,domaine,tol) + // + + [lhs,rhs]=argn(0) + select typeof(a) + case "state-space" then + if lhs<>1 then + error(msprintf(gettext("%s: Wrong number of output arguments: %d expected.\n"),"minreal",1)), + end; + select rhs + case 1 then + istol = %f; + case 2 then + istol = %t, + tol = b, + else + error(msprintf(gettext("%s: Wrong number of input arguments: %d or %d expected.\n"),"minreal",1,2)), + end; + [a,b,c,d,x0,dom] = a(2:7); + if dom == [] then + error(96,1); + end + domaine="c"; + if dom<>"c" then + domaine="d"; + end + case "constant" then + if lhs<>3 then + error(msprintf(gettext("%s: Wrong number of output arguments: %d expected.\n"),"minreal",3)), + end; + select rhs + case 4 then istol = %f + case 5 then istol = %t, + else + error(msprintf(gettext("%s: Wrong number of input arguments: %d or %d expected.\n"),"minreal",4,5)); + end; + else + error(91,1); + end; + // + wc = lyap(a', -b*b', domaine); + wo = lyap(a, -c'*c, domaine); + if ~istol then + [r,n] = equil1(wc,wo); + else + [r,n] = equil1(wc,wo,tol); + end; + n1 = n(1); + ri = inv(r); + r = r(1:n1,:); + ri = ri(:,1:n1); + a = r*a*ri; + b = r*b; + c = c*ri; + if lhs == 1 then + a = syslin(dom,a,b,c,d,r*x0); + end +endfunction |