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 xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns4="http://www.w3.org/1999/xhtml" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:id="rankqr" xml:lang="en">
15
<refname>rankqr</refname>
16
<refpurpose>fatora��o QR com revela��o do posto</refpurpose>
19
<title> Seq��ncia de Chamamento </title>
20
<synopsis>[Q,R,JPVT,RANK,SVAL]=rankqr(A, [RCOND,JPVT])</synopsis>
23
<title>Par�metros</title>
28
<para>matriz de reais ou complexos</para>
34
<para> n�mero real usado para determinar o posto efetivo de
35
<literal>A</literal>, que � definido como sendo a ordem da maior
36
submatriz regente triangular <literal>R11</literal> na fatora��o QR
37
com pivoteamento de <literal>A</literal>, cujo n�mero de
38
condicionamento estimado � < <literal>1/RCOND</literal>.
46
vetor de inteiros nas entradas, se <literal>JPVT(i)</literal>
47
n�o � 0, a <literal>i</literal>-�simo coluna de <literal>
50
�permtutada para a frente de <literal>AP</literal>,
51
sen�o, a coluna <literal>i</literal> � uma coluna livre. Na sa�da,
52
se <literal>JPVT(i) = k</literal>, ent�o a
53
<literal>i</literal>-�sima coluna de <literal>A*P</literal> era a
54
<literal>k</literal>-�sima coluna de <literal>A</literal>.
62
posto efetivo de <literal>A</literal>, i.e., a ordem da
63
submatriz <literal>R11</literal>. � o mesmo que a ordem da submatriz
64
<literal>T1</literal> na fatora��o ortogonal completa de
72
<para>vetor de reais com 3 componentes; as estimativas de alguns dos
73
valores singulares do fator triangular <literal>R</literal>.
76
<literal>SVAL(1)</literal> � o maior valor singular de
77
<literal>R(1:RANK,1:RANK)</literal>;
80
<literal>SVAL(2)</literal> � o menor valor singular de
81
<literal>R(1:RANK,1:RANK)</literal>;
84
<literal>SVAL(3) </literal>� o menor valor singular de
85
<literal>R(1:RANK+1,1:RANK+1)</literal>, se <literal>RANK</literal>
86
< <literal>MIN(M,N)</literal>, ou de
87
<literal>R(1:RANK,1:RANK)</literal>, caso contr�rio.
94
<title>Descri��o</title>
95
<para>Computa (opcionalmente) uma fatora��o QR com revela��o do posto de
96
uma matriz de reais geral M-por-N, ou de complexos <literal>A</literal>,
97
que pode ser deficiente de posto, e estima seu posto efetivo usando
98
estimativa de condi��o incremental.
100
<para>A rotina usa uma fatora��o QR com pivoteamento de colunas:</para>
101
<programlisting role=""><![CDATA[
102
A * P = Q * R, onde R = [ R11 R12 ],
106
com <literal>R11</literal> definida como a maior submatriz regente
107
cujo n�mero de condi��o estimado � menor que <literal>1/RCOND</literal>. A
108
ordem de <literal>R11</literal>, <literal>RANK</literal>, � o posto
109
efetivo de<literal>A</literal>.
111
<para>Se a fatora��o triangular revela o posto (que ser� o caso se as
112
colunas regentes forem bem condicionadas), ent�o
113
<literal>SVAL(1)</literal> tamb�m ser� uma estimativa para o maior valor
114
singular de <literal>A</literal>, e <literal>SVAL(2)</literal> e
115
<literal>SVAL(3)</literal> ser�o estimativas para o
116
<literal>RANK</literal>-�simo e <literal>(RANK+1)</literal>-�simo valores
117
singulares de <literal>A</literal>, respectivamente.
119
<para>Examinando-se estes valores, pode-se confirmar que o posto � bem
120
definido a respeito do valor escolhido de <literal>RCOND</literal>. A
121
raz�o <literal>SVAL(1)/SVAL(2)</literal> � uma estimativa do n�mero de
122
condicionamento de <literal>R(1:RANK,1:RANK)</literal>.
126
<title>Exemplos</title>
127
<programlisting role="example"><![CDATA[
128
A=rand(5,3)*rand(3,7);
129
[Q,R,JPVT,RANK,SVAL]=rankqr(A,%eps)
133
<title> Ver Tamb�m </title>
134
<simplelist type="inline">
136
<link linkend="qr">qr</link>
139
<link linkend="rank">rank</link>
144
<title> Fun��es Usadas</title>
145
<para>Rotinas da biblioteca Slicot MB03OD, ZB03OD.</para>