Decomposição LU
Transcript of Decomposição LU
Eliminação de Gauss e Decomposição LU
Profa. Dra. Marli de Freitas Gomes HernandezCESET-UNICAMP
Histórico• Uma versão preliminar da eliminação de Gauss apareceu pela
primeira vez no livro chinês “Nove Capítulo de Artes Matemática”, em torno de 200 a.C. Até então o poder do método não tinha sido reconhecido.
• No ano de 1801 Carl Friedich Gauss utilizou o método para calcular a órbita do asteróide Ceres com pouquíssimas informações (anotações do astrônomo siciliano Giuseppe Piazzi quem batizou o asteróide com o nome ao observar-lo pela primeira vez).
• O trabalho de Gauss causou sensação quando Ceres reapareceu na constelação de virgem, local aproximado aos seus cálculos.
• Mais tarde o método foi popularizado quando Willian Jordan (engenheiro alemão) em 1888 publicou no seu livro de geodésica intitulado “Handbuch der Vermessungskund”.
• Embora as idéias tenham sido conhecidas antes, muitos vezes o credito pela popularização da decomposição LU é atribuída ao lógico e matemático britânico Alan Turing (precursor do computador), pelo seu trabalho de 1948 nesse assunto.
• Ao final dos anos 1970, a Fundação Nacional de Ciências e o Departamento de Energia dos EUA financiaram o desenvolvimento de rotinas de computacionais para inverter matrizes e resolver sistemas de equações lineares. Aquele pesquisa levou a um conjunto de programas Fortran chamada LINPAC que são uma referencia para muitos algoritmos computacionais de hoje. Inclusive o chamado MATLAB. As rotinas LIMPAC estão organizadas em torno de quatro fatorações de matrizes, uma das quais é a decomposição LU. C.B. Moler, J.J. Dongarra, G.W. Stewart e J.R. Brunch, os principais programadores do LINPAC, basearam muitas de suas idéias no trabalho de Jemes Boyle e Kenneth Dritz, do Laboratório Argonne (nos EUA).
Informações retiradas de [1]
Objetivo• Resolver um Sistema de equações lineares do tipo:
• onde aij ,i = 1,2,...,m e j=1,2,...,n coeficientes,
bi, i = 1,2,...,m constantes,
xj, j=1,2,...,n incógnitas.
)1.1(
...
...
...
2211
22222121
11212111
mnmnmm
nn
nn
bxaxaxa
bxaxaxabxaxaxa
• O sistema (1.1) pode ter:– Mais equações do que incógnitas (m > n);– Mais incógnitas do que equações (m < n);– O mesmo número de incógnitas e equações
(m = n).
• A solução de (1.1) podem ser:– Única;– Infinitas;– Não existente.
Operações elementares entre equações sem alterar o resultado
As operações elementares entre equações de um sistema linear do tipo (1.1) são:
1. Trocar as equações de posição
2. Multiplicar uma ou mais equações por constantes (chamamos múltiplos de equações):
3. Somar o múltiplo de uma equação por outra
Se aplicarmos qualquer operação elementar entre equações, em um sistema linear o resultado (x1,x2,...,xn) sempre será o mesmo como veremos a seguir sem demonstração.
Trocar as equações de posição:
mnmnmm
qnqnqq
pnpnpp
nn
bxaxaxa
bxaxaxa
bxaxaxa
bxaxaxa
...
...
...
...
2211
2211
2211
11212111
mnmnmm
pnpnpp
qnqnqq
nn
bxaxaxa
bxaxaxa
bxaxaxa
bxaxaxa
...
...
...
...
2211
2211
2211
11212111
1849
20232
zyxzyx
zyx
202329184
zyxzyxzyx
Sistema 1 Sistema 2
Exemplo: Dado o seguinte sistema:
Solução do Sistema 1 = Solução do Sistema 2, x=3 y=2 e z=4.
mnmnmm
qnqnqq
pnpnpp
nn
bxaxaxa
bxaxaxa
bxaxaxa
bxaxaxa
...
...
...
...
2211
2211
2211
11212111
mnmnmm
qnqnqq
pnpnpp
nn
bxaxaxa
bxaxaxa
bxaxaxa
bxaxaxa
...
...
...
...
2211
2211
2211
11212111
Multiplicar uma ou mais equações por constantes (chamamos múltiplos de equações):
Exemplo: Dado o Sistema 1:
1849
20232
zyxzyx
zyx
362289
20232
zyxzyx
zyx
Sistema 1 Sistema 3
2)184( zyx
Solução do Sistema 1 = Solução do Sistema 3, x=3 y=2 e z=4.
Somar o múltiplo de uma equação por outra:
mnmnmm
qnqnqq
pnpnpp
nn
bxaxaxa
bxaxaxa
bxaxaxa
bxaxaxa
...
...
...
...
2211
2211
2211
11212111
mnmnmm
qpnqnpnqpqp
qnqnqq
nn
bxaxaxa
bbxaaxaaxaa
bxaxaxa
bxaxaxa
...
)()(...)()(
...
...
2211
222111
2211
11212111
+
9
2)20232(
zyx
zyx
1844957520232
zyxzyxzyx
Exemplo: Dado o Sistema 1:
1849
20232
zyxzyx
zyx
Sistema 1
Sistema 4
Solução do Sistema 1 = Solução do Sistema 4, x=3 y=2 e z=4.
mnmnmm
nn
nn
bxaxaxa
bxaxaxabxaxaxa
...
...
...
2211
22222121
11212111
)3.1(
...
...
...
2
1
2
1
21
22221
11211
bvetor
m
Xvetor
n
AMatrix
mnmm
n
n
b
bb
x
xx
aaa
aaaaaa
Sistema na forma Matricial
)2.1(bAX
Sistema na forma de equações lineares
Colocar o sistema de equações lineares (1.1) na forma matricial
Podemos abreviar (1.3) escrevendo-o em forma de arranjo retangular de números denominado Matriz Aumentada do sistema.
Esse termo Matriz Aumentada foi introduzida pelo matemático norte americano Bôcher no seu livro “Introduction to Higher Algebra” em 1907. [1]
bvetor
m
Xvetor
n
AMatrix
mnmm
n
n
b
bb
x
xx
aaa
aaaaaa
2
1
2
1
21
22221
11211
...
...
...
)4.1(
...
...
...
2
1
21
22221
11211
AumentadaMatriz
mmnmm
n
n
b
bb
aaa
aaaaaa
Sistema na forma Matricial
Matriz Aumentada do sistema
AumentadaMatrizMatricial
Forma
bAbAX
O primeiro exemplo conhecido do uso de uma matriz aumentada para descrever sistemas lineares aparece no livro chinês “Nove Capítulos de Arte Matemática” publicado entre 200 a.C. e 100 a.C.durante a dinastia de Han.
• Problema proposto pelo manuscrito: Existem três tipos de milho, dos quais três montes do primeiro, dois do segundo e um do terceiro totalizam 39 medidas. Dois montes do primeiro, três do segundo e um do terceiro totalizam 34 medidas. Finalmente, um monte do primeiro, dois do segundo e três do terceiro totalizam 26 medidas. Quantas medidas de milho estão contidas em um monte de cada um dos tipos?
• O Problema leva a um sistema linear de três equações e três incógnitas, que o autor escreve como:
O arranjo do autor é colocado em colunas e e não em linhas com colchetes, como mostrado em (1.4). )5.1(
393426113232321
Informações retiradas de [1]
)6.1(391233413226321
321
321
321
xxxxxxxxx
393426
123132321
3
2
1
xxx
Aproveitando o sistema proposto em (1.5), vamos usá-lo como exemplo e colocá-lo em forma de sistema de equações (1.1), forma matricial (1.3), e na forma de matriz aumentada (1.4)
391233413226321
Forma de sistema de equações lineares
Matriz aumentada do sistema (1.6)Forma matricial do sistema (1.6)
Resolução de sistemas triangulares superiores da forma:Supondo que a matriz Anxn(quadrada) do sistema seja não singular, que
implica que os elementos da diagonal são não zero.
).7.1(......
22222
11212111
a
bxa
bxaxabxaxaxa
nnnn
nn
nn
).7.1(
...00
...0
...
2
1
2
1
222
11211
b
b
bb
x
xx
a
aaaaa
BAX
bvetor
n
Xvetor
n
AMatrix
nn
n
n
).7.1(
...00
...0
...
2
1
222
11211
c
b
bb
a
aaaaa
AumentadaMatriz
nnn
n
n
Forma de sistema de equações lineares
Forma matricial do sistema (1.7.a) Matriz aumentada do sistema (1.7.a)
nnnn
nnnnnnn
nn
nnnn
bxaxabxa
xaxabxaxaxaxabxa
)1()1()1()1)(1(
23232222
1)1()1(13121111
...
...
Solução de (1.7.a)
nnnn
nnnnnnn
nn
nnnn
abxaxabx
axaxabxaxaxaxabx
//)(
/)...(/)...(
)1)(1()1()1()1(
22232322
111)1()1(131211
Passo 1 - Explicitar aiixi i=1,2,....,n.
Passo 2 – Dividir a equação i por aii para obter xi, i=1,2,....,n.
Inversa de uma matriz M triangular superior com diagonal principal com elementos unitários, M-1.
).18.1(
1...00
...10
...1
2
112
bccc
C n
n
E considerando
veremos que a seguinte igualdade é satisfeita
).18.1(
1...00
...10
...1
2
112
aaaa
M n
n
Seja
IC
n
n
M
n
n
ccc
aaa
1...00
0...100...01
1...00
...10
...1
1...00
...10
...1
2
112
2
112
(1.18.b) e a inversa de (1.18.a).
Sabendo que MM-1=I , se M-1 = C a seguinte igualdade é satisfeita
Mais tarde será mostrado como calcular a inversa de (1.18.a), é mais fácil do que a inversa de A em (1.17.b) e em (1.3) sendo A quadrada(m=n) e não singular.
Veremos como resolver o sistema (1.7.a) na forma matricial.Aplicando a operação elementar, multiplicando em cada linha i a constante
1/aii, i = 1, 2,..,n, em (1.7.a) , as soluções dos dois sistemas serão a mesma.
nnnnnn
nn
nn
abxa
abxaxa
abxaxaxa
1)(
1)...(
1)...(
2222222
1111212111
nn
nn
nn
nn
abx
abx
aax
abx
aax
aax
222
2
22
22
11
1
11
12
11
121
...
...
Colocando na forma matricial
*
22
2
11
1
2
1
22
2
11
1
11
12
100
10
1
b
nn
n
X
n
E
n
n
ab
abab
x
xx
aaaa
aa
Como E = D-1A é da forma (1.18.a)
nnn
n
n
b
bb
be
a
aaaaa
A
sendo
2
1
222
11211
00
0
nnAdediagonal
nn
a
a
a
D
a
aa
D
100
010
001
00
0000
22
11
122
11
ADE
n
n
A
nn
n
n
D
nn
aaaa
aa
a
aaaaa
a
a
a
AD
11
100
10
1
00
0
100
010
001
22
2
11
1
11
12
222
11211
22
11
1
.
100
010
001
*1
22
2
11
1
2
1
22
11
1
b
nn
n
b
n
D
nn ab
abab
b
bb
a
a
a
bD
e
111
1)(
2
111
1
22
2
11
1
11
12
100
101
100
10
1
ADE
n
n
ADE
n
n
eee
aaaa
aa
bAXbDDAXbDDAX
bDADXoubEX
bEXbDAXDbAX
bAXbEXPova
I
1111111
111*1
*11
1*1
)(
)(
Para resolver 1.17.b, basta calcular:
bAX
bEX
beEED
1
*1
*11 ,,,,
Exemplo: Resolver o seguinte sistema linear:
944136322
3
32
321
xxxxxx
bA
xxx
946
400130322
3
2
1
940041306322
Forma de sistema de equações lineares
Forma matricial do sistema acimaMatriz aumentada do sistema
400030002
,946
,400130322
DbA
4934
3
493426
,
10031102311
440031
330
23
22
22
,
4100
0310
0021
*
111
bC
bDADD
*1
*
bEX
bEXbAX
4934
3
,
10031106711
,
10031102311
*1 bEE
XbE
49
1272423
4934
3
10031106711
*1
X
X
49
1272423
Solução
1
12
4923
32
23
31
2313
43
43
31
49
31
2313
12
23
23
2323
13
13
2313
010230
31
e
e
e
ee
e
e
ee
eeeee
• Como E é uma matriz 3x3, considerada pequena, ela será determinada algebricamente da forma rudimentar:
100010001
100
101
10031102311
100
101
,
10031102311
23
12111
23
12111 e
eeEEe
eeEE
100010001
10031102311
100010001
100
101
10031102311
2313
231312
23
13121
ee
eee
eee
EE
Inversão de E
Eliminação de Gauss
• Como visto, é muito mais fácil resolver sistemas lineares triangulares superiores em forma de sistemas de equações.
• E extremamente fácil na forma AX=B(matricial), A triangular superior.
• As mesmas operações elementares entre equações, são válidas para linhas da matriz aumentada (1.4).
• Usando operações elementares entre linhas na matriz aumentada ou equações no sistema de equações lineares, transformar um sistema linear qualquer em sistema linear triangular superior.
• Como visto anteriormente, usando as operações elementares entre equações no sistema de equações lineares ou entre linhas na matriz aumentada a solução do sistema permanece a mesma.
Eliminação de Gauss visa:
Eliminação de Gauss visa transformarusando operações elementares: Vamos representar elementos não nulos por ”*”
originalSistema
xxx
xxxxxx
n
n
n
**...**
**...****...**
21
21
21
dotransformaSistema
x
xxxxx
n
n
n
**
**...***...**
2
21
Operações elementaresentre equações
bvetorXvetor
n
AMatrix
x
xx
*
**
*...**
*...***...**
2
1
**
*
**
*...00
*...*0*...**
2
1
bvetorXvetor
n
AMatrix
x
xx
Operações elementares
Sistema original
Sistema transformado
*
**
*...**
*...***...**
*
**
*...00
*...*0*...**
Operaçõeselementaresentre linhas
Matriz aumentada do sistema original
Matriz aumentada do sistema Transformado
Como aplicar a eliminação de Gauss no sistema forma matriz aumentada - usando operações elementares entre linhas.
Aqui será adotado o seguinte:
• Operar o sistema na forma matriz aumentada, no meu ponto de vista, é mais claro, fácil e menos trabalhoso.
• Desejando, pode operar também na forma de sistemas de equações, ou até mesmo na forma matricial.
• Supor que a matriz A seja quadrada m=n e não singular.
• Pode aplicar eliminação de Gauss em matrizes singulares( se quadrada) ou com m ≠ n também. Esses casos serão tratados mais tarde .
• Adotado as seguintes notações; aij(k) e b i
(k), i = 1,2,...,m( i-ésima linhas) e j=1,2,...,n (j-ésima colunas) e k=1,...(k-ésima etapa da eliminação).
A eliminação (ou pivoteamento) se procede da esquerda para a direita, de cima para baixo, abaixo da diagonal
principal.
Pivô
Pivôs das fazes anteriores a k
AumentadaMatriz
kn
kk
k
k
kknk
kkkk
kkk
kkkk
b
b
bb
aa
aa
aaaaaaa
nn
kn
nk
nk
)(
)(
)(2
)(1
)()(
)()(
)()()(
)()()()(
00
00
02222
111211
Elementos a serem eliminados na faze k
Na fase k , escolhe-se o elemento pivô akk(k) (elemento referência) situado na posição da diagonal principal da coluna k e linha k.
O pivô akk (k) não será eliminado(zerado), somente os
elementos abaixo dele. Pivôs das fazes anteriores a k
pivô
AumentadaMatriz
kn
kk
k
k
kknk
kkkk
kkk
kkkk
b
b
bb
aa
aa
aaaaaaa
nn
kn
nk
nk
)(
)(
)(2
)(1
)()(
)()(
)()()(
)()()()(
00
00
02222
111211
Elementos a serem eliminados na faze k
• Caso o elemento akk(k) for zero ou próximo de zero,
escolher outro elemento abaixo da diagonal principal, na mesma coluna, apk
(k) ,não zero e p>k.• O pivô dessa coluna será apk
(k)
Posição do pivô, mas, a22
(2) = 0Pivô da faze 1
ap2(2) ≠ 0
AumentadaMatriz
n
p
nk
pk
b
b
bb
aaa
aaa
aaaaaaa
nnn
pnp
nk
nk
)2(
)2(
)2(2
)2(1
)2()2()2(
)2()2()2(
)2()2()2(
)2()2()2()2(
2
2
2222
111211
0
0
0
• Colocar a linha p na posição da linha k e vice versa e eliminar os elementos abaixo da posição do pivô.
• Exemplo: k=2.
Pivô da faze 1pivô da fase 2
AumentadaMatriz
n
p
nk
k
b
b
bb
aaa
aaa
aaaaaaa
nnn
n
pnpkp
nk
)2(
)2(2
)2(
)2(1
)2()2()2(
)2()2(2
)2(
)2()2()2(
)2()2()2()2(
2
222
2
111211
0
0
0
Elementos a serem eliminados na faze 2
• Continuar a eliminação (ou pivoteamento) até que k=n ou a posição do pivô seja ann
(n) e todo triangulo inferior à diagonal principal seja 0 (zero).
ultimo pivô
Eliminação de Gauss Terminada
Agora, é só terminar de resolver o sistema, basta usar o método já mostrado aqui, para sistemas triangulares superiores.
AumentadaMatriz
nn
nk
n
n
n
nnkk
nnn
nnnn
b
b
bb
a
aa
aaaaaaa
nn
kn
nk
nk
)(
)(
)2
)(1
)(
)()(
)()()(
)()()()(
000
00
02222
111211
Como fazer as operações elementares na eliminação (ou
pivotamento) de Gauss. Para cada fase k = 1,2,..,n, da eliminação (ou pivoteamento):
Determinar o pivô akk(k) ≠0 (ou não muito pequeno).
Aplicando operações elementares entre linhas.
Para cada elemento aik(k) que deverá ser eliminado (zerado), na i-ésima linha, i
= k+1,...,n, a abaixo da k-ésima da linha do pivô na mesma k-ésima coluna, determinar uma constante mik, de modo que, ao multiplicá-la pela k-ésima linha do pivô e somar com a i-ésima linha, esse elemento deverá ser zerado.
Valor do elemento aik(k)
na fase k
Valor do pivô akk(k) na
fase k
)(
)(
)1()()( 0
kkk
kik
ik
ikkkkik
kkk
kkkik
aam
aamaam
Exemplo: seja a11(1) o pivô. O objetivo, é zerar todos os elementos ai1
(1) i = 2,...,n, na coluna 1, abaixo da linha 1. Isso é:
)1(
1)1(
111
)1(1
)1(111 0
ii
ii
aam
aam
)1(1
)1(1
)1(13
)1(12
)1(11 baaaa n
)1(1
)1(1
)1(13
)1(12
)1(11 baaaa n
Pivô da fase 1 )1(11
)1(1
1 aam i
i
)1(
11
)1(31
aa
)1(
11
)1(21
aa
)1(
11
)1(1
aan )1(
1)1(
1)1(
13)1(
12)1(
11 baaaa n
)2(
)2(3
)2(2
)2(1
)2()2(3
)2(2
)2(3
)2(33
)2(32
)2(2
)2(23
)2(22
)2(1
)2(13
)2(12
)2(11
0
00
nnnnn
n
n
n
b
bbb
aaa
aaaaaaaaaa
)1(
)1(3
)1(2
)1(1
)1()1(3
)1(2
)1(1
)1(3
)1(33
)1(32
)1(31
)1(2
)1(23
)1(22
)1(21
)1(1
)1(13
)1(12
)1(11
nnnnnn
n
n
n
b
bbb
aaaa
aaaaaaaaaaaa
)1(1
)1(1
)1(13
)1(12
)1(11 baaaa n
Fase 1
)2(
2)2(
222
)2(2
)2(222 0
ii
ii
aam
aam
)2(2
)2(2
)2(23
)2(220 baaa n
pivô da fase 2
)2(22
)2(2
2 aam i
i
)2(
22
)2(32
aa
)2(
22
)2(2
aan
Fase 2
)2(2
)2(2
)2(23
)2(220 baaa n
)2(2
)2(2
)2(23
)2(220 baaa n
)3(
)3(3
)3(2
)3(1
)3()3(3
)3(3
)3(33
)3(2
)3(23
)3(22
)3(1
)3(13
)3(12
)3(11
00
000
nnnn
n
n
n
b
bbb
aa
aaaaaaaaa
)2(
)2(3
)2(2
)2(1
)2()2(3
)2(2
)2(3
)2(33
)2(32
)2(2
)2(23
)2(22
)2(1
)2(13
)2(12
)2(11
0
00
nnnnn
n
n
n
b
bbb
aaa
aaaaaaaaaa
)1(1
)1(1
)1(13
)1(12
)1(11 baaaa n
)3(3
)3(3
)3(3300 baa n
)3(
3)3(
333
)3(3
)3(333 0
ii
ii
aam
aam
)3(3
)3(3
)3(3300 baa n
pivô da fase 3
)3(33
)3(3
3 aam i
i
)3(
33
)3(3
aan
Fase 3
)2(2
)2(2
)2(23
)2(220 baaa n
)1(1
)1(1
)1(13
)1(12
)1(11 baaaa n
)3(
)2(3
)3(2
)3(1
)3()3(3
)3(3
)3(33
)3(2
)3(23
)3(22
)3(1
)3(13
)3(12
)3(11
00
000
nnnn
n
n
n
b
bbb
aa
aaaaaaaaa
)4(
)4(3
)4(2
)4(1
)4(
)4(3
)4(33
)4(2
)4(23
)4(22
)4(1
)4(13
)4(12
)4(11
000
000
nnn
n
n
n
b
bbb
a
aaaaaaaaa
Fase nParar
)3(3
)3(3
)3(3300 baa n
)2(2
)2(2
)2(23
)2(220 baaa n
)1(1
)1(1
)1(13
)1(12
)1(11 baaaa n
pivô da fase n (ultima)
)(
)(3
)(2
)(1
)(
)(3
)2(33
)(2
)(23
)(22
)(1
)(13
)(12
)(11
000
000
nn
n
n
n
nnn
nn
nn
nn
nn
nnn
b
bbb
a
aaaaaaaaa
Agora, basta terminar de resolver o sistema.
O sistema proposto em (1.5), no livro chinês, será usado como exemplo de eliminação de Gauss.
Pivô da fase 1 313,2
12
)1(11
)1(31
31)1(11
)1(21
21 aam
aamFase 1
391233413226321
26321
26321
12
13
3426)2(13)2(32)2(21)2(
3926)3(13)3(22)3(31)3(
3984018510
26321
Pivô da fase 2
4)1()4(
)2(22
)2(32
32
aam
Fase 2
18510
)1()4(
39)18()4(8)5()4(4)1()4(0
3984018510
26321
33120018510
26321
Agora é só terminar de resolver o sistema equivalente triangular superior.
Esse sistema é muito mais fácil de resolver, do que o original, tanto pelo método rudimentar como pelo de eliminação de Gauss.
Como já foi visto eliminação de Gauss, será usado eliminação de Gauss para terminá-lo, mas antes vamos colocá-lo em uma outra forma equivalente triangular inferior com diagonal principal unitária para facilitar ainda mais a resolução.
391233413226321
321
321
321
xxxxxxxxx
121)3312(
)1(1)1851(
11)26321(
3
32
321
x
xx
xxx
12331852632
3
32
321
x
xxxxx
Sistema equivalente Sistema original
Aplicando eliminação de Gauss
Colocar na forma de matriz aumentada com equações e icóginitas em ordem invertidas
Sistema equivalente triangular superior diagonal principal
unitária.
Dividir cada linha pelo respectivo elemento da
diagonal
33
Sistema equivalente forma matriz aumentada triangular inferior
diagonal principal unitária.
1
2
3
,26123180151233001
xxx
X
De agora em diante, para mostrar as operações, será colocado à frente da linha pivô i valor mij que irá multiplicar-la, na forma “(mij)” e uma seta desde esse valor até a linha a qual será somado.
+
+
)3()5(
26123180151233001
+
)2(
1221312012510101233001
sistemadosolução
x
x
x
X
4374
174
11
,
1211110012510101233001
1
2
3
• A eliminação de Gauss para esses tipos de sistema, continua sendo como já foi visto. Mas pode-se acontecer de:
• Caso obtenha equações(linhas) toda de zeros Basta colocá-las no final das equações (linhas).
Sistemas lineares com n≠m ou matriz A singular (determinante de A)=0.
Exemplo 1 Forma matriz elementar
++
+
+
+
+
000200000460332
~
00043023040
332)
21(
~
430460230000332
~
1094868
56299633)3()1()4()2(
6
2
pivô
pivô
0=2 significa (não existe) solução (obviamente 0≠2)
• Caso obtenha colunas de zeros desde a linha do pivô(inclusive), busque outro pivô na primeira coluna à direita na mesma linha.
Sistema fica com duas equações e três incógnitas. Significa que existe infinitas soluções.
Para cada α (constante) existe uma solução
2
2~
000042006321
~
21004300
6321
)21(
~
210042006321
~
4221168426321)2()1(
3
2
1
321
321
321
321
21
321
321
321
321
321
321
32
2
1
xxx
xxxxxxxxx
xxxxx
xxx
xxxxxxxxx
xxxxxxxx
x
x
pivô
pivô
+
+
+
Decomposição LUUma decomposição LU ou uma fatoração LU de uma matriz quadrada A e uma fatoração A=LU na qual L é triangular inferior e U triangular superior.
)(
Re
sin:
1
11
XdesoluçãoYUXUXYUX
bLYbLYbUXLbXALLbAX
solverComo
gularnãonnquadradaMatrizAseja
YU
1111
sin:
LULUAAinverter
gularnãonnquadradaMatrizLUAseja
)2(
)1(
)2(
21311
)2()2(3
)2(2
)2(3
)2(33
)2(32
)2(2
)2(13
)2(22
)2(1
)2(13
)2(12
)2(11
)1()1(3
)1(2
)1(1
)1(3
)1(33
)1(32
)1(31
)1(2
)1(13
)1(22
)1(21
)1(1
)1(13
)1(12
)1(11
)1(11
)1(1
)1(11
)1(31
)1(11
)1(21
)2()2(3
)2(2
)2(3
)2(33
)2(32
)2(2
)2(13
)2(22
)2(1
)2(13
)2(12
)2(11
)1()1(3
)1(2
)1(1
)1(3
)1(33
)1(32
)1(31
)1(2
)1(13
)1(22
)1(21
)1(1
)1(13
)1(12
)1(11
)1(11
)1(21
)1(11
)1(31
)1(11
)1(1
0
00
100
010
001
0001
0
00
~
A
nnnn
n
n
n
A
nnnnn
n
n
n
M
n
A
nnnn
n
n
n
A
nnnnn
n
n
n
mmm
n
aaa
aaaaaaaaaa
aaaa
aaaaaaaaaaaa
aa
aaaa
aaa
aaaaaaaaaa
aaaa
aaaaaaaaaaaa
aa
aa
aan
Decomposição LU é feita usando eliminação de Gauss, registrando em uma matriz diagonal unitária, os valores multiplicados pela linha pivô ii com o
objetivo de somar às linhas (k=i,i+1,...n) para eliminar(zerar) os elemento ki,
++
+
)3()2(
)2(
)3()2(
211
)3()3(3
)3(3
)3(33
)3(2
)3(13
)3(22
)3(1
)3(13
)3(12
)3(11
)2()2(3
)2(2
)2(3
)2(33
)2(32
)2(2
)2(13
)2(22
)2(1
)2(13
)2(12
)2(11
)2(22
)2(2
)2(22
)2(32
)3()3(3
)3(3
)3(33
)3(2
)3(13
)3(22
)3(1
)3(13
)3(12
)3(11
)2()2(3
)2(2
)2(3
)2(33
)2(32
)2(2
)2(13
)2(22
)2(1
)2(13
)2(12
)2(11
)2(22
)2(32
)1(11
)1(1
00
000
0
00
100
0100
00100001
00
000
~
0
00
A
nnn
n
n
n
A
nnnn
n
n
n
M
n
A
nnn
n
n
n
A
nnnn
n
n
nmm
n
aa
aaaaaaaaa
aaa
aaaaaaaaaa
aa
aa
aa
aaaaaaaaa
aaa
aaaaaaaaaa
aa
aan
O mesmo que
++
U
nn
n
n
n
L
nn
UA
nnn
nn
n
nn
nn
nn
nnn
A
nnnnn
n
n
n
M
n
M
n
M
nnn
nnn
u
uuuuuuuuu
UAMMMM
a
aaaaaaaaa
aaaa
aaaaaaaaaaaa
aa
aaaa
aa
aa
aa
n
n
000
000
...
000
000
100
010
001
0001
100
0100
00100001
100
01000000100001
333
21322
1131211
)1()2()2()1(
)(
)(3
)(33
)(2
)(13
)(22
)(1
)(13
)(12
)(11
)1()1(3
)1(2
)1(1
)1(3
)1(33
)1(32
)1(31
)1(2
)1(13
)1(22
)1(21
)1(1
)1(13
)1(12
)1(11
)1(11
)1(1
)1(11
)1(31
)1(11
)1(21
)2(22
)2(2
)2(22
)2(32
)1()1)(1(
)1(1(
1
)(
)1()2(
)1(
1)2()2(
1)1()1(
100
0100
00100001
100
0100
00100001
100
010
001
0001
100
010
001
0001
)2(22
)2(2
)2(22
)2(32
1
)2(22
)2(2
)2(22
)2(32
)1(11
)1(1
)1(11
)1(31
)1(11
)1(21
1
)1(11
)1(1
)1(11
)1(31
)1(11
)1(21
M
n
M
n
M
n
M
n
aa
aa
aa
aa
aa
aaaa
aa
aaaa
1)1()1(
100
01000000100001
100
01000000100001
)1()1)(1(
)1(1(
1
)1()1)(1(
)1(1(
nnM
nnn
nnn
M
nnn
nnn
aa
aa
AUMMMM
aaaa
aaaaaaaaaaaa
a
aaaaaaaaa
aa
aa
aa
aa
aaaa
L
nn
A
nnnnn
n
n
n
UA
nnn
nn
n
nn
nn
nn
nnn
M
nnn
nnn
M
n
M
n
n
n
1)1(1)2(1)2(1)1(
)1()1(3
)1(2
)1(1
)1(3
)1(33
)1(32
)1(31
)1(2
)1(13
)1(22
)1(21
)1(1
)1(13
)1(12
)1(11
)(
)(3
)(33
)(2
)(13
)(22
)(1
)(13
)(12
)(11
)1()1)(1(
)1(1(
)2(22
)2(2
)2(22
)2(32
)1(11
)1(1
)1(11
)1(31
)1(11
)1(21
...
000
000
100
01000000100001
100
0100
00100001
100
010
001
0001
)(1)1(
1)2(1)1(
LMMMM
nnnn
nn
LMMMM
nnn
nnnnn
nn
M
nnn
nnn
M
n
M
n
nn
nn
n
mmmmm
m
aa
aa
aa
aa
aa
aa
aa
aa
aa
aa
aaaa
1)1(1)2(1)2(1)1(
1)1(1)2(1)2(1)1(
1)1(1)2(1)1(
...
)1(21
2)1(1)1(
21
...
)1()1)(1(
)1(1(
)2(22
)2(2
)1(11
)1(1
)2(22
)2(2)1(
)1(11
)1(1)1(
)1(11
)1(21
)1()1)(1(
)1(1(
)2(22
)2(2
)2(22
)2(32
)1(11
)1(1
)1(11
)1(31
)1(11
)1(21
101000010001
1
01
00
001
0001
100
01000000100001
100
0100
00100001
100
010
001
0001
L
nnnn
nn
lllll
l
101000010001
)1(21
2)1(1)1(
21
Como se percebe, a matriz U é toda de zeros abaixo da diagonal principal e a matriz L é toda de zeros acima da unitária diagonal principal . Computacionalmente, para economizar memória, a matriz L e U são armazenadas em uma só matriz e um vetor K com registro das trocas de linhas feitas durante a decomposição LU.
U
nn
n
n
n
L
nnnn
nn
u
uuuuuuuuu
e
lllll
l
000
000
101000010001
333
21322
1131211
)1(21
2)1(1)1(
21
K
n
LeUmatrizesdasntoArmazename
nnnnn
n
n
n
k
kkk
e
ulll
uulluuuluuuu
3
2
1
321
3333231
2132221
1131211
Ki é o índice da k-ésima
linha original A.
Para exemplificar voltaremos ao exemplo do livro chinês (1.6).
Para que o exemplo seja completo usaremos a técnica do pivô sendo o maior elemento da coluna, desde a linha do pivô para baixo.
3121
)0()1()0()0(
31
32
321132123
123
123132321
321
mm
AKpivômaior
AK
+ +
L
L
UA
K
m
pivômaior
A
K
mmm
154
31
0132
001
101001
5120031
350
123
123
35/
34
38
340
31
350
123
123
3231
21
54
)3(
)3(
23
)2(
)2(
+
UeL
K
512
54
31
31
35
32
123,
123
Arm
azen
amen
to d
e L
e U
1
12500121
530
363
52
31
100010001
32
12500121
530
3650
31
100010
0321
31
51
12500
0530
0031
100511031
321
1255331
100010001
5120031
350
123
123
U
I
I
U
K
+
+
+
1
154
51
0132
001
100010001
54
1031
0132
001
1540
010001
31
32
100010001
154
31
0132
001
L
I
I
L
+ +
+
sistemadosoluçãoX
b
LU
XK
LULU
xxx
4114
174
37
263439
125
31
121
121
32
125
121
31
127
123
125
31
121
121
32
125
121
31
127
154
51
0132
001
12500121
530
121
52
31
11
1111
3
2
1
Bibliografia[1] ANTON, H. & BUSBY, R. Algebra Linear Contemporânea. Editora
Bookman. Porto Alegre. 2006.[2] RUGGIERO, M.G. & LOPES, V.L.R. Cálculo Numérico – Aspectos
Computacionais, Pearson Education. São Paulo. 1996.