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:ns5="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="bdiag" xml:lang="en">
15
<refname>bdiag</refname>
16
<refpurpose>diagonaliza��o em blocos, autovetores
21
<title> Seq��ncia de Chamamento </title>
22
<synopsis>[Ab [,X [,bs]]]=bdiag(A [,rmax])</synopsis>
25
<title>Par�metros</title>
30
<para>matriz quadrada de reais ou complexos </para>
36
<para>n�mero real</para>
42
<para>matriz quadrada de reais ou complexos</para>
48
<para>matriz de reais ou complexos n�o-singular</para>
54
<para>vetor de inteiros</para>
60
<title>Descri��o</title>
61
<programlisting role=""><![CDATA[
62
[Ab [,X [,bs]]]=bdiag(A [,rmax])
65
realiza a diagonaliza��o em blocos da matriz <literal>A</literal>.
66
bs fornece a estrutura dos blocos (respectivos tamanhos dos blocos).
67
<literal>X</literal> � a matriz mudan�a de base i.e <literal>Ab =
73
<literal>rmax</literal> controla o condicionamento de
74
<literal>X</literal>; o valor padr�o � a norma L1 de
77
<para>Para encontrar a forma diagonal (se existir) escolha um valor
78
suficientemente grande para <literal>rmax</literal>
79
(<literal>rmax=1/%eps</literal> , por exemplo). Genericamente, (para uma
80
matriz A de reais aleat�ria) os blocos s�o (1x1) e (2x2) e
81
<literal>X</literal> � a matriz de autovetores.
85
<title>Exemplos</title>
86
<programlisting role="example"><![CDATA[
87
//No caso de reais: blocos 1x1 e 2x2
88
a=rand(5,5);[ab,x,bs]=bdiag(a);ab
90
//No caso de complexos: blocos 1x1 complexos
91
[ab,x,bs]=bdiag(a+%i*0);ab
95
<title> Ver Tamb�m </title>
96
<simplelist type="inline">
98
<link linkend="schur">schur</link>
101
<link linkend="sylv">sylv</link>
104
<link linkend="spec">spec</link>