Sesion5 Teo MetodosNumericos 2023 2
Sesion5 Teo MetodosNumericos 2023 2
Sesion5 Teo MetodosNumericos 2023 2
Numéricos
M. Iterativos S.E.L - S5
Hermes Pantoja Carhuavilca
([email protected])
Rósulo Perez Cupe
([email protected])
Jimmy Mendoza Montalvo
([email protected])
Máximo Obregón Ramos
([email protected])
Jose Mantari Laureano
([email protected])
Profesores: Utec-Ciencias
Índice
1 Método de Jacobi
2 Método de Gauss Seidel
3 Convergencia de los
Métodos Iterativos
Logros de Aprendizaje
{x (k ) }, k = 0, 1, 2, . . .
lim x (k ) = x.
k →∞
x = Tx + c
x (k +1) = Tx (k ) + c, para k = 0, 1, 2, . . .
Si logramos despejar:
1
(b1 − a12 x2 − a13 x3 − · · · − a1n xn )
x1 =
a11
1
x2 = (b2 − a21 x1 − a23 x3 − · · · − a2n xn )
a22
.. .. ..
. . .
1
xn = (bn − an1 x1 − an2 x2 − · · · − ann−1 xn−1 )
ann
entonces podemos definir un proceso iterativo:
Universidad de Ingeniería y Tecnología Métodos Numéricos September 7, 2023 6 / 40
Forma algebraica del método de Jacobi
(k +1) 1
(k ) (k ) (k )
x1 = b1 − a12 x2 − a13 x3 − · · · − a1n xn
a11
(k +1) 1 (k ) (k ) (k )
x2 = b2 − a21 x1 − a23 x3 − · · · − a2n xn
a22
.. .. ..
. . .
(0)
x1
(0)
x2
(k +1) 1 (k ) (k ) (k )
xn = bn − an1 x1 − an2 x2 − · · · − ann−1 xn−1 , con x (0) = ..
ann
.
(0)
xn
un punto semilla dado, en general:
n
(k +1) 1 bi −
X (k )
xi = aij xj , para i = 1, 2, . . . , n siempre que aii ̸= 0
aii
j=1,j̸=i
−1 3 1 8 210
(D − L − U)x = b
Dx = (L + U)x + b
x = D −1 (L + U)x + D −1 b
Entonces queda definido el método iterativo:
x (k +1) = D −1 (L + U)x (k ) + D −1 b
(k +1) 90 (k ) 2 (k ) 1 (k ) 1 (k )
x1 = + 0x1 + x2 + x3 + x4
6 6 6 6
(k +1) −150 4 (k ) (k ) 3 (k ) 2 (k )
x2 = + x1 + 0x2 + x3 − x4
8 8 8 8
(k +1) 300 1 (k ) 2 (k ) (k ) 0 (k )
x3 = − x1 − x2 + 0x3 + x4
5 5 5 5
(k +1) 210 1 (k ) 3 (k ) 1 (k ) (k )
x4 = + x1 − x2 − x3 + 0x4
8 8 8 8
(0) (0) (0) (0)
Iniciando en el punto x1 = 20 ; x2 =5 ; x3 = 30 ; x4 = 10
(1) 90 2 1 1
x1 = + 0(20) + (5) + (30) + (10) = 23.333
6 6 6 6
(1) −150 4 3 2
x2 = + (20) + 0(5) + (30) − (10) = 0.000
8 8 8 8
(1) 300 1 2 0
x3 = − (20) − (5) + 0(30) + (10) = 54.000
5 5 5 5
(1) 210 1 3 1
x4 = + (20) − (5) − (30) + 0(10) = 23.125
8 8 8 8
(2)90 2 1 1
x1 + 0(
= )+ ( )+ ( )+ ( ) =
6 6 6 6
(2) −150 4 3 2
x2 = + ( ) + 0( )+ ( )− ( ) =
8 8 8 8
(2) 300 1 2 0
x3 = − ( )− ( ) + 0( )+ ( ) =
5 5 5 5
(2) 210 1 3 1
x4 = + ( )− ( )− ( ) + 0( ) =
8 8 8 8
asi sucesivamente.
k =2
k =3
x = Tx + c ⇐⇒ x k +1 = Tj x k + cj
Es decir:
2 1 1 90
0
(k +1) (k )
x1 6 6 6 x1 6
(k +1)
4 3 2
(k )
−150
x2
0 − x2
8 8 8 8
= +
(k +1) 1 2 0 (k ) 300
x3 − − 0 x3
5 5 5 5
(k +1)
1 3 1
(k )
210
x4 − − 0 x4
| 8 8 {z 8 } | 8
{z }
Tj =D −1 (L+U) cj =D −1 b
Universidad de Ingeniería y Tecnología Métodos Numéricos September 7, 2023 16 / 40
MÉTODO DE GAUSS
SEIDEL
2
Método de Gauss Seidel
Observe que en el método de Jacobi, cuando se calculan las componentes del
vector x (k +1) , sólo se usan las componentes del vector x (k ) , sin embargo, note que
(k +1) (k +1) (k +1)
para obtener xi , se podrían haber usado las componentes x1 hasta xi−1
porque ellas ya han sido calculadas.
n
(k +1) (k )
X
Jacobi: xi = (bi − aij xj )/aii , para i = 1, 2, . . . , n
j=1,j̸=i
i−1 n
(k ) (k )
X X
= (bi − aij xj − aij xj )/aii , para i = 1, 2, . . . , n
j=1 j=i+1
i−1 n
(k +1) (k +1) (k )
X X
Gauss-Seidel: xi = (bi − aij xj − aij xj )/aii , para i = 1, 2, . . . , n
j=1 j=i+1
90 2 1 1
x1 = + 0x1 + x2 + x3 + x4
6 6 6 6
−150 4 3 2
x2 = + x1 + 0x2 + x3 − x4
8 8 8 8
300 1 2 0
x3 = − x1 − x2 + 0x3 + x4
5 5 5 5
210 1 3 1
x4 = + x1 − x2 − x3 + 0x4
8 8 8 8
(k +1) 90 (k ) 2 (k ) 1 (k ) 1 (k )
x1 = + 0x1 + x2 + x3 + x4
6 6 6 6
(k +1) −150 4 (k +1) (k ) 3 (k ) 2 (k )
x2 = + x1 + 0x2 + x3 − x4
8 8 8 8
(k +1) 300 1 (k +1) 2 (k +1) (k ) 0 (k )
x3 = − x1 − x2 + 0x3 + x4
5 5 5 5
(k +1) 210 1 (k +1) 3 (k +1) 1 (k +1) (k )
x4 = + x1 − x2 − x3 + 0x4
8 8 8 8
(0) (0) (0) (0)
Iniciando en el punto x1 = 20 ; x2 =5 ; x3 = 30 ; x4 = 10
(1) 90 2 1 1
x1 = + 0(20) + (5) + (30) + (10) = 23.333
6 6 6 6
(1) −150 4 3 2
x2 = + (23.333) + 0(5) + (30) − (10) = 1.666
8 8 8 8
(1) 300 1 2 0
x3 = − (23.333) − (1.666) + 0(30) + (10) = 54.666
5 5 5 5
(1) 210 1 3 1
x4 = + (23.333) − (1.666) − (54.666) + 0(10) = 21.708
8 8 8 8
(2)90 2 1 1
x1 + 0(
= )+ ( )+ ( )+ ( ) =
6 6 6 6
(2) −150 4 3 2
x2 = + ( ) + 0( )+ ( )− ( ) =
8 8 8 8
(2) 300 1 2 0
x3 = − ( )− ( ) + 0( )+ ( ) =
5 5 5 5
(2) 210 1 3 1
x4 = + ( )− ( )− ( ) + 0( ) =
8 8 8 8
asi sucesivamente.
k =2
k =3
(D − L − U)x = b
(D − L)x = Ux + b
x = (D − L)−1 Ux + (D − L)−1 b
Se define el método iterativo como:
1 1 1
15
0
(k +1) (k )
x 3 6 6 x
1 1
45
1 11 1
(k +1) 0 (k ) −
x2 − x2 4
6 24 6
= +
123
(k +1)
x3 0 − 2 13 1 (k )
x3
−
2
15 60 30
789
(k +1)
1 119 19
(k )
x4 0 − − x4
| 240 {z 960 240 } | 32 {z }
Tgs =(D−L)−1 U cgs =(D−L)−1 b
A=D−L−U
x (k +1) = Tx (k ) + c
Método T c
Jacobi −1
D (L + U) −1
D b
Gauss-Seidel (D − L)−1 U (D − L)−I b
Teorema
La sucesión x (k +1) = Tx (k ) +c, para k ≥ 0 converge a la solución única x = Tx +c
si y sólo si ρ(T ) < 1.
2 2 2
λ1 = 0, λ2 = =⇒ ρ(Tgs ) = Max{|0|, | |} =
15 15 15
2
4 Dado que el radio espectral ρ(Tgs ) = 15 < 1 entonces el Método de
Gauss-Seidel es convergente.
Nota: El método también es convergente por ser la matriz de coeficientes
diagonal estrictamente dominante.
Ea = ||x (k +1) − x (k ) ||
El estimado del error relativo en la iteración k + 1 es:
||x (k +1) − x (k ) ||
Er = .
||x (k +1) ||