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
40
41
42
43
44
45
46
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
|
