5
[ ] [ b ] [ ] [ ] ] [ g 11 g 21 g 31 g 41 0 g 22 g 32 g 33 0 0 g 42 g 43 0 0 0 g 44 [ ] G t = MÉTODO DE CHOLESKY Dada uma matriz simétrica, ou seja, na qual os elementos acima e abaixo da diagonal principal são coincidentes, fazendo com que a sua transposta seja igual à matriz original, é possível aplicar o Método de Cholesky, desde que esta matriz também seja positiva e definida. Considere o sistema: 10X 1 +7X 2 +8X 3 +7X 4 = 32 7X 1 +5X 2 +6X 3 +5X 4 = 23 8X 1 +6X 2 +10X 3 +9X 4 = 33 7X 1 +5X 2 +9X 3 +10X 4 = 31 Sabe-se então que: 10 7 8 7 32 7 5 6 5 X = 23 8 6 10 9 33 7 5 9 10 31 Assume-se a primeira matriz (A) como da forma: a 11 a 12 a 13 a 14 a 21 a 22 a 23 a 24 a 31 a 32 a 33 a 34 a 41 a 42 a 43 a 44 Como a matriz inicial é simétrica, positiva e definida, pode-se simplificar os cálculos do Método de Decomposição por LU, calculando-se GG t , sendo G uma matriz triangular inferior com diagonal positiva. Logo: g 11 0 0 0 G = g 21 g 22 0 0

Método de Cholesky

Embed Size (px)

DESCRIPTION

Breve explicação sobre o método numérico de Cholesky.

Citation preview

MTODO DE CHOLESKYDada uma matriz simtrica, ou seja, na qual os elementos acima e abaixo da diagonal principal so coincidentes, fazendo com que a sua transposta seja igual matriz original, possvel aplicar o Mtodo de Cholesky, desde que esta matriz tambm seja positiva e definida.Considere o sistema:

10X1+7X2+8X3+7X4 = 327X1+5X2+6X3+5X4 = 238X1+6X2+10X3+9X4 = 337X1+5X2+9X3+10X4 = 31Sabe-se ento que:

10787

32

7565 X

= 23

86109

33

75910

31

Assume-se a primeira matriz (A) como da forma:

a11a12a13a14

a21a22a23a24

a31a32a33a34

a41a42a43a44

Como a matriz inicial simtrica, positiva e definida, pode-se simplificar os clculos do Mtodo de Decomposio por LU, calculando-se GGt, sendo G uma matriz triangular inferior com diagonal positiva. Logo:g11000

G = g21g2200

g31g32g330

g41g42g43g44

Os termos de G so dados como:g11 = = g21 = a21/g11 = g22 = = = g31 = a31/g11 = 8/ = 4g32 = (g21xg31)] = (x)] = = g33 = = = = g41 = = 7/ = 7g42 = [a42 (g41xg21)] = [5 (x)] = )] = g43 = [a43 (g31xg41)- (g32xg42)] = [9 (x) - (x)] = [9 ] = 3g44 = = = = = Ento, fica:

000 32

00 23

0 33

3

31 Assim:Y1 = = Y2 = 23 - x = 23 - = = Y3 = 33 x x = 33 - - = = Y4 = 31 - x - x - 3 x = 31 - - - = = E a Gt ser:

X4 = = 1; X3 = - 3 x 1 + = = 1

X2 = - x 1 - x 1 + = = 1; X1 = - x 1 - x 1 - x1 + = = 1Ento:

10X1+7X2+8X3+7X4 = 32

7X1+5X2+6X3+5X4 = 23

8X1+6X2+10X3+9X4 = 33

7X1+5X2+9X3+10X4 = 31O Mtodo de Gauss com Pivoteamento Parcial

1) O elemento akk(k) o pivot do K passo.

2) Se em algum passo K encontrarmos akk(k) = 0, isso significa que det (Ak) = 0.

Nesse caso, o sistema ainda pode ter soluo determinada (basta que det (A) 0).

O mtodo pode ser continuado simplesmente permutando a K equao com qualquer outra abaixo cujo coeficiente da K incgnita seja 0.

3) Anlise de propagao de erros de arredondamento para o algortmo de Gauss indicam a convenincia de serem todos multiplicadores (as constantes aik(k)/akk(k) do K passo) menores que 1 em mdulo; ou seja o pivot deve ser o elemento de maior valor absoluto da coluna, da diagonal (inclusive) para baixo.

Podemos ento em cada passo, escolher na coluna correspondente o elemento de maior valor absoluto, da diagonal (inclusive) para baixo, e fazer uma permutao nas equaes do sistema, de modo que esse elemento venha a ocupar a posio diagonal.]

]

]

[

[

[

b

]

[

]

]

[

g11g21g31g41 0g22g32g33

00g42g43

000g44

[

Gt =

[

]

]

[

=

G=

]

]

[

[

QUOTE QUOTE QUOTE QUOTE 0 QUOTE QUOTE QUOTE

00 QUOTE 3 QUOTE

000 QUOTE

Gt =

=

10 x 1 + 7 x 1 + 8 x 1 + 7 x 1 = 32

7 x 1+5 x 1 + 6 x 1 + 5 x 1 = 23

8 x 1 + 6 x 1 + 10 x 1 + 9 x 1 = 33

7 x 1+5 x 1 + 9 x 1 +10 x 1 = 31