1
<?xml version="1.0" encoding="ISO-8859-1"?>
3
* Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
4
* Copyright (C) 2008 - INRIA
6
* This file must be used under the terms of the CeCILL.
7
* This source file is licensed as described in the file COPYING, which
8
* you should have received as part of this distribution. The terms
9
* are also available at
10
* http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
13
<refentry version="5.0-subset Scilab" xml:id="spantwo" xml:lang="en"
14
xmlns="http://docbook.org/ns/docbook"
15
xmlns:xlink="http://www.w3.org/1999/xlink"
16
xmlns:svg="http://www.w3.org/2000/svg"
17
xmlns:ns4="http://www.w3.org/1999/xhtml"
18
xmlns:mml="http://www.w3.org/1998/Math/MathML"
19
xmlns:db="http://docbook.org/ns/docbook">
21
<pubdate>$LastChangedDate$</pubdate>
25
<refname>spantwo</refname>
27
<refpurpose>soma e interse��o de subespa�os</refpurpose>
31
<title> Seq��ncia de Chamamento </title>
33
<synopsis>[Xp,dima,dimb,dim]=spantwo(A,B, [tol])</synopsis>
37
<title>Par�metros</title>
44
<para>duas matrizes de reais ou complexos com igual n�mero de linhas
53
<para>matriz quadrada n�o-singular </para>
58
<term>dima, dimb, dim</term>
61
<para>inteiros, dimens�es dos subespa�os </para>
69
<para>n�mero real n�o-negativo</para>
76
<title>Descri��o</title>
78
<para>Dadas duas matrizes <literal>A</literal> e <literal>B</literal> com
79
o mesmo n�mero de linhas, retorna uma matriz quadrada
80
<literal>Xp</literal> (n�o-singular, mas n�o necessariamente ortogonal)
85
[A1, 0] (dim-dimb linhas)
86
Xp*[A,B]=[A2,B2] (dima+dimb-dim linhas)
87
[0, B3] (dim-dima linha)
92
<para>As primeiras <literal>dima</literal> colunas de
93
<literal>inv(Xp)</literal> geram Im(<literal>A</literal>).</para>
95
<para>As colunas de <literal>dim-dimb+1</literal> at�
96
<literal>dima</literal> de <literal>inv(Xp)</literal> geram a interse��o
97
de Im(A) e Im(B).</para>
99
<para>As primeiras <literal>dim</literal> colunas de
100
<literal>inv(Xp)</literal> geram
101
Im(<literal>A</literal>)+Im(<literal>B</literal>).</para>
103
<para>As colunas de <literal>dim-dimb+1</literal> at�
104
<literal>dim</literal> de <literal>inv(Xp)</literal> geram
105
Im(<literal>B</literal>).</para>
107
<para>A matrix <literal>[A1;A2]</literal> tem posto-linha cheio (=
108
posto(A)), a matrix <literal>[B2;B3]</literal> tem posto-linha cheio
109
(=posto(B)), a matriz <literal>[A2,B2]</literal> tem posto-linha (=posto(A
110
inter B)) e a matriz <literal>[A1,0;A2,B2;0,B3]</literal> tem posto-linha
111
cheio (=posto(A+B)).</para>
115
<title>Exemplos</title>
117
<programlisting role="example">
123
B=[1,2,0,0]';C=[4,0,0,1];
124
Sl=ss2ss(syslin('c',A,B,C),rand(A));
125
[no,X]=contr(Sl('A'),Sl('B'));CO=X(:,1:no); //Parte control�vel
126
[uo,Y]=unobs(Sl('A'),Sl('C'));UO=Y(:,1:uo); //Parte inobserv�vel
127
[Xp,dimc,dimu,dim]=spantwo(CO,UO); //Decomposi��o de Kalman
128
Slcan=ss2ss(Sl,inv(Xp));
134
<title> Ver Tamb�m </title>
136
<simplelist type="inline">
137
<member><link linkend="spanplus">spanplus</link></member>
139
<member><link linkend="spaninter">spaninter</link></member>