diff options
Diffstat (limited to 'modules/cacsd/macros/armac.sci')
| -rwxr-xr-x | modules/cacsd/macros/armac.sci | 47 |
1 files changed, 47 insertions, 0 deletions
diff --git a/modules/cacsd/macros/armac.sci b/modules/cacsd/macros/armac.sci new file mode 100755 index 000000000..bf71e6a34 --- /dev/null +++ b/modules/cacsd/macros/armac.sci @@ -0,0 +1,47 @@ +// 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 [ar]=armac(a,b,d,ny,nu,sig) + // just build a tlist for storing armacx coefficients + // A(z^-1)y= B(z^-1)u + D(z^-1)sig*e(t) + // a=<Id,a1,..,a_r>; matrix (ny,r*ny) + // b=<b0,.....,b_s>; matrix (ny,(s+1)*nu) + // d=<Id,d1,..,d_p>; matrix (ny,p*ny); + // ny : dim of observation y + // nu : dim of control u + // sig : standard deviation (ny,ny); + // + //! + + [na,la]=size(a); + if na<>ny then + error(msprintf(gettext("%s: Wrong size for input argument #%d: row dimension must be equal to %d.\n"),.. + "armac",1, ny)); + + end + [nb,lb]=size(b); + if nb<>0 & nb<>ny then + error(msprintf(gettext("%s: Wrong size for input argument #%d: row dimension must be equal to %d.\n"),.. + "armac",2, ny)); + end; + if lb<>0 & nu<>0 then + if modulo(lb,nu)<>0 then + error(msprintf(gettext("%s: Wrong size of input argument #%d: Number of columns are incompatible with %s.\n"),.. + "armac",2,"nu")); + end; + end + [nd,ld]=size(d); + if nd<>ny then + error(msprintf(gettext("%s: Wrong size for input argument #%d: row dimension must be equal to %d.\n"),.. + "armac",3, ny)); + + end + ar=tlist(["ar","a","b","d","ny","nu","sig"],a,b,d,ny,nu,sig); +endfunction |