13
13
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:lang="en" xml:id="findAC">
15
<refname>findAC</refname>
16
<refpurpose> discrete-time system subspace identification</refpurpose>
19
<title>Calling Sequence</title>
20
<synopsis>[A,C] = findAC(S,N,L,R,METH,TOL,PRINTW)
21
[A,C,RCND] = findAC(S,N,L,R,METH,TOL,PRINTW)
25
<title>Arguments</title>
30
<para>integer, the number of block rows in the block-Hankel matrices</para>
48
<para>matrix, relevant part of the R factor of the concatenated block-Hankel matrices computed by a call to findr.</para>
54
<para>integer, an option for the method to use</para>
59
<para> MOESP method with past inputs and outputs;</para>
65
<para> N4SID method;</para>
77
<para>the tolerance used for estimating the rank of matrices. If TOL > 0, then the given value of TOL is used as a lower bound for the reciprocal condition number. Default: prod(size(matrix))*epsilon_machine where epsilon_machine is the relative machine precision.</para>
83
<para>integer, switch for printing the warning messages.</para>
88
<para>= 1: print warning messages;</para>
94
<para>do not print warning messages.</para>
106
<para>matrix, state system matrix</para>
112
<para>matrix, output system matrix</para>
118
<para>vector of length 4, condition numbers of the matrices involved in rank decision</para>
124
<title>Description</title>
126
finds the system matrices A and C of a discrete-time system, given the
127
system order and the relevant part of the R factor of the concatenated
128
block-Hankel matrices, using subspace identification techniques (MOESP
133
<para>[A,C] = findAC(S,N,L,R,METH,TOL,PRINTW) computes the system matrices A and C. The model structure is: x(k+1) = Ax(k) + Bu(k) + Ke(k), k >= 1, y(k) = Cx(k) + Du(k) + e(k), where x(k) and y(k) are vectors of length N and L, respectively.</para>
136
<para>[A,C,RCND] = findAC(S,N,L,R,METH,TOL,PRINTW) also returns the vector RCND of length 4 containing the condition numbers of the matrices involved in rank decisions.</para>
140
Matrix R, computed by findR, should be determined with suitable arguments
145
<title>Examples</title>
146
<programlisting role="example"><![CDATA[
15
<refname>findAC</refname>
16
<refpurpose> discrete-time system subspace identification</refpurpose>
19
<title>Calling Sequence</title>
20
<synopsis>[A,C] = findAC(S,N,L,R,METH,TOL,PRINTW)
21
[A,C,RCND] = findAC(S,N,L,R,METH,TOL,PRINTW)
25
<title>Arguments</title>
30
<para>integer, the number of block rows in the block-Hankel matrices</para>
48
<para>matrix, relevant part of the R factor of the concatenated block-Hankel matrices computed by a call to findr.</para>
54
<para>integer, an option for the method to use</para>
59
<para> MOESP method with past inputs and outputs;</para>
65
<para> N4SID method;</para>
77
<para>the tolerance used for estimating the rank of matrices. If TOL > 0, then the given value of TOL is used as a lower bound for the reciprocal condition number. Default: prod(size(matrix))*epsilon_machine where epsilon_machine is the relative machine precision.</para>
83
<para>integer, switch for printing the warning messages.</para>
88
<para>= 1: print warning messages;</para>
94
<para>do not print warning messages.</para>
106
<para>matrix, state system matrix</para>
112
<para>matrix, output system matrix</para>
118
<para>vector of length 4, condition numbers of the matrices involved in rank decision</para>
124
<title>Description</title>
126
finds the system matrices A and C of a discrete-time system, given the
127
system order and the relevant part of the R factor of the concatenated
128
block-Hankel matrices, using subspace identification techniques (MOESP
133
<para>[A,C] = findAC(S,N,L,R,METH,TOL,PRINTW) computes the system matrices A and C. The model structure is: x(k+1) = Ax(k) + Bu(k) + Ke(k), k >= 1, y(k) = Cx(k) + Du(k) + e(k), where x(k) and y(k) are vectors of length N and L, respectively.</para>
136
<para>[A,C,RCND] = findAC(S,N,L,R,METH,TOL,PRINTW) also returns the vector RCND of length 4 containing the condition numbers of the matrices involved in rank decisions.</para>
140
Matrix R, computed by findR, should be determined with suitable arguments
145
<title>Examples</title>
146
<programlisting role="example"><![CDATA[
147
147
//generate data from a given linear system
148
148
A = [ 0.5, 0.1,-0.1, 0.2;
149
149
0.1, 0, -0.1,-0.1;
165
165
[A,C] = findAC(S,N,L,R,METH,TOL);
166
166
]]></programlisting>
168
<refsection role="see also">
169
<title>See Also</title>
170
<simplelist type="inline">
172
<link linkend="findABCD">findABCD</link>
175
<link linkend="findBD">findBD</link>
178
<link linkend="findBDK">findBDK</link>
181
<link linkend="findR">findR</link>
184
<link linkend="sorder">sorder</link>
187
<link linkend="sident">sident</link>
168
<refsection role="see also">
169
<title>See Also</title>
170
<simplelist type="inline">
172
<link linkend="findABCD">findABCD</link>
175
<link linkend="findBD">findBD</link>
178
<link linkend="findBDK">findBDK</link>
181
<link linkend="findR">findR</link>
184
<link linkend="sorder">sorder</link>
187
<link linkend="sident">sident</link>