CURSO DE MÉTODOS NUMÉRICOS PARA INGENIERÍA
INGENIERÍA EN ENERGÍA
SEMESTRE 2022 - II
DOCENTE: MSc. LIZANDRO BALDOMERO REYNA ZEGARRA
METODO DE PUNTO FIJO
Definición de Punto Fijo: Un punto fijo de una función g es un numero real talque
Geometricamente los puntos fijos de una función son los puntos de intersección de y la recta
.
𝑦
𝑦 =𝑥
𝑦 =𝑔 ( 𝑥 )
( 𝑥 , 𝑔 (𝑥 ))
𝑥
Definición de iteración de punto fijo: La iteración para es llamado iteración de punto fijo.
Asi:
Si
Si
Si
.
.
.
TEOREMA 1: Asuma que es una función continua y que es una sucesión generada por iteración de
punto fijo. Si , entonces es un punto fijo de .
INTERPRETACIÓN GRAFICA DE LA ITERACIÓN DE PUNTO FIJO
En las siguientes figuras se muestran los casos de convergencia y divergencia
CONVERGENCIA OSCILANTE:
𝑦
𝑦 =𝑥
( 𝑝0 , 𝑔 ( 𝑝0 ) )
( 𝑝 1 , 𝑝1 )
( 𝑝 2 , 𝑔 ( 𝑝 2) )
( 𝑝2 , 𝑝2 ) ( 𝑝1 , 𝑔 ( 𝑝 1 ) )
𝑦 =𝑔 ( 𝑥 )
𝑝 0𝑝 2𝒑𝑝𝑝
3 1 𝑥
CONVERGENCIA MONOTONA
𝑦
𝑦 =𝑥
( 𝑝 1 , 𝑝1 )
( 𝑝0 , 𝑔 ( 𝑝0 ) )
(𝑝2 , 𝑝2 )
( 𝑝1 , 𝑔 ( 𝑝 1 ) )
𝑦 =𝑔 ( 𝑥 )
( 𝑝 2 , 𝑔 ( 𝑝 2) )
𝒑 p2 p1 p0 𝑥
DIVERGENCIA OSCILANTE:
𝑦 =𝑔 ( 𝑥 )
𝑦 𝑦 =𝑥
( 𝑝 2 , 𝑔 ( 𝑝 2) ) ( 𝑝3 , 𝑝3 )
( 𝑝 0 , 𝑔 ( 𝑝 0 ) ) ( 𝑝 1 , 𝑝1 )
( 𝑝2 , 𝑝2 ) ( 𝑝1 , 𝑔 ( 𝑝 1 ) )
( 𝑝3 , 𝑔 ( 𝑝3 ))
2𝐩
pp 0 p1 p3 𝑥
DIVERGENCIA MONOTONA:
𝑦 𝑦 =𝑔 ( 𝑥 )
( 𝑝 2 , 𝑔 ( 𝑝 2) ) 𝑦 =𝑥
( 𝑝1 , 𝑔 ( 𝑝 1 ) )
( 𝑝2 , 𝑝2 )
( 𝑝0 , 𝑔 ( 𝑝0 ) ) ( 𝑝 , 𝑝 )
1 1
𝐩𝑝𝑝
0 1𝑝 2 𝑥
Aquí, el propósito del desarrollo de Punto Fijo es usarlo como una técnica para aproximar raíces de
ecuaciones no lineales. La idea es resolver la ecuación no lineal,
… (1)
Para llevar a cabo esto se escribe en la forma , y se tiene la ecuación
… (2)
En lugar de la ecuación
Luego, de se tiene que
… (3)
Luego, si la solución existe, es claro que esta solución satisface (1) si y solo si satisface (3). El mecanismo
usado fundamentalmente es la iteración de punto fijo junto con el teorema 1.
EJEMPLO
Aproximar la raíz de la ecuación
SOLUCION
Paso 1: De se tiene que . Luego,
Paso 2: Usando el método grafico determinar el valor inicial
1.5
Paso 3: Verificar la convergencia
La derivada evaluada en es 0.137. Es decir, . Luego la iteración converge.
Paso 4: Iterar
𝑝 1=𝑔 ( 𝑝 0 ) =√ √ 1.5+1=1.49156
𝑝 2=𝑔 ( 𝑝 1 ) =√ √ 1.49156+1=1.49040
𝑝 3=𝑔 ( 𝑝 2 ) =√ √ 1.49040+1=1.49024
𝑝 4 =𝑔 ( 𝑝 3 ) =√ √ 1.49024+1=1.49022
𝑝 5=𝑔 ( 𝑝 4 )= √ √1.49022+1=1.49022
Resumiendo en una tabla
OBSERVACION: La convergencia de depende de como se obtenga
EJEMPLO
Hallar una raíz positiva de la ecuación que genera la intersección de las curvas y
SOLUCION
Paso1: igualar los miembros derechos de las ecuaciones de las curvas
= de donde 𝑥 2 −2 𝑥 − √ 9− 𝑥2 =0
Paso2: Escribir la ecuación del paso 1, en la forma , de donde
Paso3: Determinar usando el método grafico
𝑝 0=2.6
Paso 4: verificar la convergencia
evaluada en es 0.05078, es decir, . Hay convergencia.
Luego, la iteración converge
Paso 5: Iterar
=2.58779
=2.58713
=2.58710
=2.58709
=2.58709
Recogiendo los valores en una tabla, se tiene
OBSERVACION: muestre que si la ecuación del paso 1 lo escribimos en la forma . , genera una iteración que
diverge. Compruebe haciendo 5 iteraciones.
EJEMPLO: Hallar una raíz de con 4 iteraciones.
[ ]
𝑥=0 1 2 0 .5 − 1.5 1
𝐿 ( 1 , :) = 0 0 0
y
0 0 0
𝑤=3
𝑓𝑜𝑟 𝑘=2
[ ]
0 0 0
𝐿= 0 0 0 𝑉 =1
0 0 0
𝑗=1
1 𝑆𝐼
𝑓𝑜𝑟 𝑘=1 ( 𝑥 −0 )
( 𝑥 − 𝑥 ( 1) ) 𝑉 =1.
𝑉 =1 𝑉 =1. 1− 0
𝑥 ( 2) − 𝑥 ( 1)
𝑗=1
𝑗=2
1 𝑁0
𝑗=2 𝑖𝑓 2 ≠ 2 𝑁 0
𝑗=3
𝑖𝑓 1≠ 2 𝑆𝐼
( 𝑥 −1 ) 𝑖𝑓 2 ≠ 3 𝑆𝐼
( 𝑥 − 𝑥 ( 2) ) 𝑉 =1.
𝑉 =1. 0 −1 ( 𝑥 − 𝑥 ( 1) ) ( 𝑥 − 𝑥 (3 ))
𝑥 ( 1) − 𝑥 ( 2) ( 𝑥 − 0) ( 𝑥 −2)
𝑉= 𝑉=
𝑥 ( 2) − 𝑥 ( 1) 𝑥 ( 2) − 𝑥 ( 3 ) 1 −0 1 −2
𝑗=3
𝑖𝑓 1≠ 3 𝑆𝐼 2
−1 𝑥 + 2 𝑥
[ ]
( 𝑥 − 𝑥 ( 2 ) ) ( 𝑥 − 𝑥 ( 3 ) ) 𝑉 = ( 𝑥 −1 ) ( 𝑥 − 2 ) 0 .5 − 1.5 1
𝑉= 0 −1 0 − 2
𝐿 ( 2 , : )= − 1 2 0
𝑥 ( 1) − 𝑥 ( 2) 𝑥 ( 1) − 𝑥 ( 3) 0 0 0
2
𝑉 = 0.5 𝑥 −1.5 𝑥+1
𝑓𝑜𝑟 𝑘=3
𝑉 =1
𝑗=1
1 𝑆𝐼
( 𝑥 − 𝑥 ( 1) ) ( 𝑥 −0 )
𝑉 =1.
𝑉 =1. 2− 0 𝐶 =[ 11 1 ]
𝑥 ( 3) − 𝑥 ( 1)
𝑗=2 𝑖𝑓 3 ≠ 2 𝑆𝐼 2
𝑃 ( 𝑥 ) = 1 𝑥 + 1 𝑥 +1
( 𝑥 − 𝑥 (1 ) ) ( 𝑥 − 𝑥 ( 2 ) ) 𝑉 = ( 𝑥 − 0 ) ( 𝑥 −1 )
𝑉=
𝑥 ( 3 ) − 𝑥 (1 ) 𝑥 ( 3 ) − 𝑥 ( 2 ) 2−0 2 −1
2
0.5 𝑥 − 0.5 𝑥
[ ]
0 .5 − 1.5 1
𝐿 ( 3 ,: )= − 0.5 1 0
0.5 − 0.5 0
