Metodo de Cholesky

9.2.
METODO DE CHOLESKY
Primera versin
Tambin para resolver el sistema
Ax=b
para aplicar cholesky se debe
cumplir lo siguiente:
1)
es simtrico, es decir debe cumplir
2)
sea definida positiva.
A= A t
Ejemplo 9.2.1.
Ax=b
Desarrolle:
[ ][ ] [ ]
4 2 1 x1
1
=
2 2 0 x2
2
1 0 3 x3
3
Dado :
[ ]
4 2 1
A= 2 2 0
1 0 3
entonces
A= A t
Viendo si es definida positiva:
A 1=[ 4 ]
[ ]
A 2= 4 2
2 2
[ ]
4 2 1
A 3= 2 2 0
1 0 3
entonces
det A 1=4 >0
entonces
entonces
det A 2=4 >0
det A 3=10>0
Por lo tanto A es positiva, entonces se puede aplicar cholesky.

Se factoriza como el producto de dos matrices triangulares.
A=LLt
[ ][
][
l 11 0 0
l 11 l21 l 31
4 2 1
2 2 0 = l 21 l 22 0 0 l22 l 32
1 0 3
l 31 l 32 l33
0 0 l 33
Primera columna:
Segunda columna:
4=( l 11 )2
Entonces
l 11=2
2=l 11l 21
2=l 21l 11
Entonces
l 21=1
2=( l 21 )2 + ( l 22 )2
1=l 31l 11
Entonces
l 31=
Entonces
0=l 21l31 +l 22l 32
1
2
l 32 =
l 22=1
Entonces
1
2
Tercera columna:
1=l 11l 31
0=l 21u13 +l 22u 23
3=( l31 )2 + ( l 32 )2+ ( l 33) 2
Entonces
[ ][
2
0
4 2 1
1
2 2 0= 1
1 0 3
1/2 1/2
Calculando
[]
x1
x= x 2
x3
[]
1
b= 2
3
Ax=b
Se sabe
A=LLt
][
0
2 1
0 0 1
10/ 2 0 0
tal que
10
l 33 =
2
1/ 2
1/2
10/2
Ax=b
]
con
[ ]
4 2 1
A= 2 2 0
1 0 3
( LLt )x=b
L( Ltx )=b
Sea
L x= y
la cual
Ly =b
[]
y1
y= y 2
y3
][ ] []
y1
0
1
0 y2 = 2
10 /2 y 3 3
2
0
1
1
1 /2 1 /2
2 y1 =1
Entonces
y 1+ y 2=2
Entonces
y 1=
1
2
y 2=
3
2
1
1
10
y 1 y 2+ y 3=2 Entonces
2
2
3
y 3=
Ltx= y
7 10
10
2 1
0 1
0 0
][ ] [ ]
x1
1/2
1/2
1/2 x 2 = 3 /2
10/2 x3 7 10/10
10 x = 7 10
3
2
x 3=
10
7
5
1
3
x 2 x3 =
2
2
x 2=
Entonces
11
5
Entonces
1
1
2 x 1 + x 2+ x 3=
2
2
x 1=
6
5
Entonces
[ ]
6/5
x= 11/5
7 /5
1.
MTODO DE NEWTON
RAPHSON
El Mtodo de Newton - Raphson es ampliamente utilizado para encontrar
f ( x )=0 , ya que converge rpidamente, la
las races de la ecuacin
contra es que uno debe conocer la derivada de
f (x)
y se necesita
una aproximacin inicial muy cercana a la raz.

Se requiere que
f (x)
sea doblemente continua y diferenciable en
[ a , b] .
Algoritmo:
P-1.-
f ( x )=0
Dado la ecuacin
tal que existe la raz
x [ a , b]
por T.B.
P-2.- Generar la
{ x i } x
mediante la relacin
x i+1=x i
f ( xi )
1
f ( xi )
i=0, 1,2, 3,
Convergencia de N-R.
Existe
{ x i } x
Esto significa que
|[
(2)
f ( x0) f ( x0)
si
xi
f ( x0)]
1
<1
f ( x 0 ) f (2) ( x0 ) > 0
est muy cercano a la raz.
P-3. -
Dejar de iterarse
m n +1
=0.5 10
para el caso
de cifras
exactas (
Significativas
)
|x ix i1| E
=0.5 10k
para el caso de
cifras
Decimales exactas
Caso contrario ir al P-2 en
Ejemplo 7.1.
2 x
f ( x )=x e 1=0 ( 1 )
a) Por Teorema de Bolzano localizar el intervalo donde exista la raz.
x 0 ,1
m=1
b) Verificar su convergencia por N.R.
f ( x )=x 2 e x 1
Si
entonces
f 1 ( x )=x ( 2+ x ) e x , f (2) ( x )=[ ( x +2 )22 ] e x

x 0=0
Veremos si
es vlido o no.
f ( 0 ) f (2) ( 0 ) <0
entonces
x 0=0
no es vlido para iterar.
x 0=1 es vlido o no
Veremos si
f ( 1 ) f (2) ( 1 ) >0
entonces
x 0=1
es vlido para iterar.
|[ ] |
(2 )
f ( 1) f ( 1)
1
f ( 1)
=0.491649373 1(cumple)
{ xi } x
por N.R.
c) En la iteracin Cuantas cifras significativas exactas tiene la solucin ?

Con
x 0=1
y con
f ( x )=x 2 e x 1
f 1 ( x )=x ( 2+ x ) e x
x i+1=x i
f ( xi )
1
f ( xi )
2
xi
( x i ) e 1
x i+1=x i
xi ( 2+ x i ) e x
Escogeremos para dos cifras significativas exactas ( n = 2 )
|x ix i1| 0.5 x 10m n+1
Se deja de iterar si
Con
m=1 ,n=2
|x ix i1| 0.5 x 102
i=0
Iteracin inicial
x 1=x 0
f ( x0 )
1
f ( x0 )
=1
f (1)
f 1 ( 1)
Entonces
|x 1x 0|=0.210706853=0.2=0.2 x 10 0 0.5 x 102 ( Falso)

i=1
Primera Iteracin
x 2=x 1
f ( x1 )
1
f ( x1 )
=0.789293147
f ( 0.789293147 )
f 1 ( 0.789293147 )
x 1=0.789293147
Entonces
x 2=0.71261243
|x 2x 1|=0.076680717=0.1=0.1 x 10 0 0.5 x 102 ( Falso)

i=2
Segunda Iteracin
x 3=x 2
f ( x2 )
1
f ( x2 )
=0.71261243
f ( 0.71261243 )
f 1 ( 0.71261243 )
Entonces
x 3=0.703582624
|x 3x 2|=0.009029806=0.01=0.1 x 101 0.5 x 102 ( Falso)

i=3
Segunda Iteracin
x 4=x 3
f ( x3 )
1
f ( x3 )
=0.703582624
f ( 0.703582624 )
f 1 ( 0.703582624 )
Entonces
x 4=0.703467441
|x 4x 3|=0.000115183=0.0001=0.1 x 103 0.5 x 102 (cumple)

Entonces
x 4=0.703467441
es solucin con dos cifras significativas
exactas.
1.1. Algoritmo del Punto Fijo en dos variables :
P-1.- De ( 1 ) despejar
respectivamente para obtener una
relacin de la siguiente forma.

x=f ( x , y )=f
(2)
y=g ( x , y )=g
P-2.- De ( 2 ) generar la sucesin:
{ x i } x , { y i } y
Mediante la siguiente relacin de recurrencia:
x i+1=f ( xi , y i )
y i+1=g ( x i , y i )
Donde
P-3. -
i=0, 1,2,
Dejar de iterarse
=0.5 10m n +1
para el caso de
cifras
exactas (
Significativas
)
|x ix i1| E
=0.5 10k
para el caso de
cifras
Decimales exactas
=0.5 10m n +1
cifras
para el caso de
exactas (
Significativas
| y i y i1| E
=0.5 10k
para el caso de
cifras
Decimales exactas
Caso contrario ir al P-2 en
Condicin de convergencia del punto fijo:
{ xi } x , { yi } y
Si se cumple lo siguiente:
|f x|( x 0 , y 0 ) +|f y|( x 0 , y 0 ) L<1
|g x|( x 0 , y 0 ) +|g y|( x 0 , y 0 ) L< 1
Donde
f x=
f (x, y) g(x , y)
=
x
x gx
Ejemplo 8.1.1.
F=F ( x , y ) =xe 2 y =0 ( 1 )
G=G ( x , y )=ex + y=0 ( 2 )
a) Por Teorema de Bolzano localizar el punto inicial
(i)
Primero formamos de (1) la forma

De (2)
De
( )
se obtiene
( x0 , y0)
F(x)
y=ex ( )
x
en (1)
se obtiene
F ( x )=xe2 e
(ii)
Aplico Teorema de Bolzano para intervalos de longitud 1.
x 0 ,1
m=1
(iii)
De (ii) igual para intervalo de longitud 0.1
x 0.1 ,0.2
Sea
x 0=0.1
Luego de
( )
y 0=ex =e0.1
0
entonces
y 0=0.9
b) Verificar convergencia
{ x i } x , { y i } y ?
F ( x , y )=x e 2 y =0
2y
se tiene
x=e
se tiene
y=e
entonces
f =f ( x , y ) =e 2 y
x
G ( x , y )=e + y=0
x
g=g ( x , y ) =e
entonces
f x =0
f y =2 e2 y
g y =0
gx =ex
|f x|( x , y ) +|f y|( x , y ) <1?

0
|f x|( 0.1,0.9 )+|f y|(0.1,0.9 ) <1

0+|2 e 2(0.9)|<1
0+0.3<1( cumple)
|g x|( x
, y0 )
+|g y|( x , y ) <1?

0
|g x|(0.1 ,0.9 ) +|g y|( 0.1,0.9 )< 1

|ex|+0<1
0.9+0<1( cumple)
Por lo tanto f y g cumplen la condicin de convergencia.
c) Por punto fijo obtener una solucin con dos cifras significativas exactas (
n=2 ).
Se dejar de iterar si:
|x ix i1| 0.5 x 10m n+1
con
|x ix i1| 0.5 x 1012+1
n=2 Entonces
| y i y i1| 0.5 x 10mn+ 1

con
| y i y i1| 0.5 x 1012+1
|x ix i1| 0.5 x 102
i=1, 2,3,
con
| y i y i1| 0.5 x 102
Si es verdadera entonces (
xi , yi
es solucin con dos cifras
significativas exactas.
Teniendo :
f =x=e 2 y
g= y =ex
entonces
entonces
x i+1=e 2 y
y i+1=ex
con
i=0
i=0, 1,2, 3,
( x 0 , y 0 ) =( 0.1 ,0.9 )
x 1=e 2 y =e 2(0.9)
Entonces
x 1=0.165299
y 1=ex =e(0.1)
Entonces
y 1=0.904837
|x 1x 0|=0.0652989=0.1=0.1 x 10 0 0.5 x 102 (no cumple)

| y 1 y 0|=0.00483742=0.004=0.4 x 102 0.5 x 102 (cumple)
i=1
( x 1 , y 1 )=( 0.165299,0.904837 )
x 2=e 2 y =e 2(0.904837)
Entonces
x 2=0.163707
y 2=ex =e(0.165299)
Entonces
y 2=0.84764
|x 1x 0|=0.00159153=0.001=0.1 x 102 0.5 x 102 (cumple)

| y 1 y 0|=0.0571971=0.1=0.1 x 100 0.5 x 102(no cumple)
i=2
( x 2 , y 2 )= ( 0.163707 ,0.84764 )
2 y1
x 2=e =e
2(0.84764 )
y 2=ex =e(0.163707)
1
Entonces
x 3=0.183548
Entonces
y 3=0.84899
|x 2x 1|=0.0198404=0.01=0.1 x 101 0.5 x 102 (no cumple)

| y 2 y 1|=0.00135012=0.001=0.1 x 102 0.5 x 102 (cumple )
i=3
( x 3 , y 3 ) =( 0.183548 ,0.84899 )
x 4=e2 y =e2(0.84899)
Entonces
x 4=0.183053
y 4 =ex =e(0.183548)
Entonces
y 4 =0.832312
|x 4x 3|=0.00159153=0.001=0.1 x 102 0.5 x 102( cumple)

| y 4 y 3|=0.0166783=0.02=0.2 x 101 0.5 x 102 (no cumple )
i=4
( x 4 , y 4 ) =( 0.183053 ,0.832312 )
2 y4
x 5=e =e
2(0.832312 )
y 5=ex =e(0.183053 )
4
Entonces
x 5=0.189262
Entonces
y 5=0.832724
|x 5x 4|=0.00620899=0.01=0.1 x 101 0.5 x 102 (no cumple)
| y 5 y 4|=0.000412058=0.0004=0.4 x 103 0.5 x 102 (cumple )
i=5
( x 5 , y 5 ) =( 0.189262,0.832724 )
x 6=e 2 y =e 2(0.832724 )
5
y 6=ex =e(0.189262)
5
Entonces
Entonces
x 6=0.189106
y 6=0.82757
|x 6x 5|=0.000155909=0.0001=0.1 x 103 0.5 x 102 (cumple )

| y 6 y5|=0.00515436=0.005=0.5 x 102 0.5 x 102 (cumple )
Por lo tanto
( x 6 , y 6 ) =( 0.189106 ,0.82757 )=( x , y )
significativas exacta.
es raz con dos cifras

Metodo de Cholesky

Cargado por

Copyright:

Formatos disponibles

Metodo de Cholesky

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodo de Cholesky

Cargado por

Copyright:

Formatos disponibles

9.2.

para aplicar cholesky se debe

es simtrico, es decir debe cumplir

sea definida positiva.

Viendo si es definida positiva:

det A 1=4 >0

det A 2=4 >0

Por lo tanto A es positiva, entonces se puede aplicar cholesky.

0=l 21l31 +l 22l 32

una aproximacin inicial muy cercana a la raz.

sea doblemente continua y diferenciable en

tal que existe la raz

Esto significa que

est muy cercano a la raz.

f 1 ( x )=x ( 2+ x ) e x , f (2) ( x )=[ ( x +2 )22 ] e x

no es vlido para iterar.

es vlido para iterar.

c) En la iteracin Cuantas cifras significativas exactas tiene la solucin ?

Escogeremos para dos cifras significativas exactas ( n = 2 )

|x ix i1| 0.5 x 10m n+1

|x ix i1| 0.5 x 102

|x 1x 0|=0.210706853=0.2=0.2 x 10 0 0.5 x 102 ( Falso)

|x 2x 1|=0.076680717=0.1=0.1 x 10 0 0.5 x 102 ( Falso)

|x 3x 2|=0.009029806=0.01=0.1 x 101 0.5 x 102 ( Falso)

|x 4x 3|=0.000115183=0.0001=0.1 x 103 0.5 x 102 (cumple)

es solucin con dos cifras significativas

1.1. Algoritmo del Punto Fijo en dos variables :

respectivamente para obtener una

relacin de la siguiente forma.

Mediante la siguiente relacin de recurrencia:

Condicin de convergencia del punto fijo:

|f x|( x 0 , y 0 ) +|f y|( x 0 , y 0 ) L<1

|g x|( x 0 , y 0 ) +|g y|( x 0 , y 0 ) L< 1

a) Por Teorema de Bolzano localizar el punto inicial

Primero formamos de (1) la forma

Aplico Teorema de Bolzano para intervalos de longitud 1.

De (ii) igual para intervalo de longitud 0.1

|f x|( x , y ) +|f y|( x , y ) <1?

|f x|( 0.1,0.9 )+|f y|(0.1,0.9 ) <1

+|g y|( x , y ) <1?

|g x|(0.1 ,0.9 ) +|g y|( 0.1,0.9 )< 1

|x ix i1| 0.5 x 10m n+1

|x ix i1| 0.5 x 1012+1

| y i y i1| 0.5 x 10mn+ 1

| y i y i1| 0.5 x 1012+1

|x ix i1| 0.5 x 102

| y i y i1| 0.5 x 102

es solucin con dos cifras

|x 1x 0|=0.0652989=0.1=0.1 x 10 0 0.5 x 102 (no cumple)

|x 1x 0|=0.00159153=0.001=0.1 x 102 0.5 x 102 (cumple)

|x 2x 1|=0.0198404=0.01=0.1 x 101 0.5 x 102 (no cumple)

|x 4x 3|=0.00159153=0.001=0.1 x 102 0.5 x 102( cumple)

|x 5x 4|=0.00620899=0.01=0.1 x 101 0.5 x 102 (no cumple)

| y 5 y 4|=0.000412058=0.0004=0.4 x 103 0.5 x 102 (cumple )

|x 6x 5|=0.000155909=0.0001=0.1 x 103 0.5 x 102 (cumple )

( x 6 , y 6 ) =( 0.189106 ,0.82757 )=( x , y )

es raz con dos cifras

También podría gustarte