Universidad Autónoma del Estado de México
UAEM
“Facultad de Ingeniería”
Formulario:
“Métodos Numéricos”
Elaborado por:
David Isaías Jaimes Reyes
Estudiante de Ingeniería en Electrónica
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
Número Máquina
…
0 1 1 0 0 0 0 1 0 1 1 0 1 0 0 0 0 1 0 1 0 1 0 1 0 1 0 0 1 0 1 0
Exponente o Carácter Mantisa (parte decimal)
Signo Signo del exponente
del
número
Numero Real = signo(mantisa)* y
Tipos de Errores
Error absoluto Error Relativo Error Porcentual
| | | |
Algoritmo del Método de la Bisección donde es una tolerancia
preestablecida
1. Verificar que la función sea y continua
en[ ] Método de Newton – Raphson
( ) ( ) ( )
[ ]
( )
2. Obtener la aproximación de la raíz mediante
Método de la Secante
( )( )
[ ]
3. Determinar en que subintervalo se encuentra la ( ) ( )
raíz
( ) ( ) [ ] Ajuste de una Línea Recta (Regresión Lineal)
( ) ( ) [ ]
( ) ( ) es la solución
- Ecuaciones Normales
4. Estimar el error ∑ ∑ ∑
∑ (∑ )
donde es una tolerancia
̅ ̅ (promedio)
preestablecida
- Error estándar del estimado (Desviación
Algoritmo del Método de la Regla Falsa (falsa
estándar)
posición, regula falsi)
1. Verificar que la función sea y continua ⁄ √
en[ ]
( ) ( )
∑( )
2. Obtener la aproximación de la raíz mediante
( )( )
[ ] ∑( ̅)
( ) ( )
- Coeficiente de Determinación
3. Determinar en que subintervalo se encuentra la
raíz
( ) ( ) [ ]
( ) ( ) [ ]
( ) ( ) es la solución - Coeficiente de Correlación
4. Estimar el error ( √ )
David Isaías Jaimes Reyes 1
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
Regresión Polinomial
∑ ∑( )
- Ecuaciones Normales
( ) (∑ ) (∑ ) ∑
(∑ ) (∑ ) (∑ ) ∑
(∑ ) (∑ ) (∑ ) ∑
- Error estándar
⁄ √
( )
∑( ̅)
Regresión Lineal Múltiple
∑( )
- Ecuaciones Normales (notación matricial)
∑ ∑ ∑
∑ ∑ ∑ ∑
[∑ ∑ ∑ ]{ } {∑ }
- Error Estandar
⁄ √
( )
Linealización de funciones no lineales
Ajuste Exponencial
- Ecuaciones Normales
∑ ∑ (∑ )
(∑ ) (∑ )
(∑ )(∑ ) (∑ )(∑ )
(∑ ) (∑ )
David Isaías Jaimes Reyes 2
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
Ajuste de Potencia Simple
- Ecuaciones Normales
∑( ) ∑( ) ∑( )
∑( ) ∑( ) ∑( )( )
Interpolación
Polinomio de Taylor
( ) ( )
( ) ( ) ( )( ) ( ) ( )
Interpolación de LaGrange
Este polinomio está dado por:
( ) ( ) ( ) ( ) ( ) ( ) ( )
donde:
( )( )( ) ( )
( )
( )( )( ) ( )
( )
∏
( )
Interpolación de Hermite
( ) ∑ ( ) ( ) ∑ ( )̂ ( )
Donde
( ) [ ( ) ( )] ( )
Y
̂ ( ) ( ) ( )
Interpolación de Newton (Diferencias divididas)
( ) ( ) ( )( ) ( )( ) ( )
Para las constantes apropiadas
( ) ( )
( ) ( )
[ ]
( ) ( ) ( ) ( )
[ ] [ ] [ ] [ ]
[ ]
David Isaías Jaimes Reyes 3
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
( ) Primeras diferencias Segundas diferencias divididas Terceras diferencias divididas
divididas
[ ]
( ) ( )
[ ]
[ ] [ ] [ ]
[ ]
( ) ( ) [ ] [ ]
[ ] [ ]
[ ] [ ] [ ]
[ ]
( ) ( ) [ ] [ ]
[ ] [ ]
[ ] [ ] [ ]
[ ]
( ) ( ) [ ] [ ]
[ ] [ ]
[ ] [ ] [ ]
[ ]
( ) ( )
[ ]
[ ]
Interpolación Iterada (Método de Neville)
( ) ( )
Criterio de paro:
| |
David Isaías Jaimes Reyes 4
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
Integración Numérica
Fórmulas de Newton – Cotes cerradas simples
Nombre Fórmula
Regla del
trapecio simple ∫ ( ) [ ( ) ( )] ( )
Regla de
∫ ( ) [ ( ) ( ) ( )] ( )
Simpson ⁄
( )
Regla de ( ) [ ( ) ( ) ( ) ( )]
∫ ( )
Simpson ⁄
Regla de Boole ( ) [ ( ) ( ) ( ) ( ) ( )]
∫
( )
Fórmulas de integración compuestas
Nombre Características Fórmula
Punto medio ⁄
Sea [ ], par y y
los nodos ( )( ) para ∫ ( ) ∑ ( ) ( )
cada
Regla del trapecio Sea [ ], y los nodos
∫ ( ) [ ( ) ∑ ( ) ( )]
para cada
( )
⁄
Simpson ⁄ Sea [ ], par y y los
nodos para cada ∫ ( ) [ ( ) ∑ ( )
∑ ( ) ( )]
( )
( )
David Isaías Jaimes Reyes 5
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
Cuadratura Gaussiana
( ) ( ) ( )
∫ ( ) ∫ ( )
∫ ( ) ∑ ( )
Raíces Coeficientes
Diferenciación Numérica
Fórmulas de primeras derivadas
Nombre Fórmula
2 puntos hacia adelante ( ) [ ( ) ( )] ( )
3 puntos hacia adelante ( ) [ ( ) ( ) ( )] ( )
3 puntos centrada ( ) [ ( ) ( )] ( )
5 puntos centrada ( ) [ ( ) ( ) ( ) ( )] ( )(
)
5 puntos hacia adelante ( ) [ ( ) ( ) ( ) ( ) ( )] ( )
( )
3 puntos hacia atrás ( ) [ ( ) ( ) ( )]
5 puntos hacia atrás ( ) [ ( ) ( ) ( ) ( ) ( )]
David Isaías Jaimes Reyes 6
UAEM Facultad de Ingeniería Formulario: Métodos Numéricos
Derivadas de orden superior ( )
Nombre Fórmula
3 puntos centrada ( ) ( ) [ ( ) ( ) ( )] [ ]
3 puntos hacia adelante ( ) ( ) [ ( ) ( ) ( )] [ ]
4 puntos hacia adelante ( ) ( ) [ ( ) ( ) ( ) ( )] [ ]
5 puntos centrada ( ) ( ) [ ( ) ( ) ( ) ( ) ( )]
[ ]
4 puntos hacia adelante( ) ( ) [ ( ) ( ) ( ) ( )] [ ]
5 puntos centrada ( ) ( ) [ ( ) ( ) ( ) ( )] [ ]
5 puntos hacia adelante ( ) ( ) [ ( ) ( ) ( ) ( ) ( )] [ ]
5 puntos centrada ( ) ( ) [ ( ) ( ) ( ) ( ) ( )] [ ]
Extrapolación de Richardson
( ) ( ) ( ) ( )
( )
( )
( )
( )
( ) ( )
( )
( ) ( ) ( )
( ) ( )
( )
David Isaías Jaimes Reyes 7
FORMULARIO DE METODOS Para la solución se tomara el METODO DE ILLINOIS
NUMERICOS siguiente algoritmo:
Los valores
METODO PARA LA SOLUCION DE (𝑥𝐷 − 𝑥𝐼 ) (𝑥𝐼 , 𝑓(𝑥𝐼 )); (𝑥𝐷 , 𝑓(𝑥𝐷 ))
𝑥𝑀 = 𝑥𝐷 − 𝑓(𝑥𝐷 )
ECUACIONES NO LINEALES 𝑓(𝑥𝐷 ) − 𝑓(𝑥𝐼 )
Si:
METODO DE PUNTO FIJO Si:
𝑓(𝑥𝐷 ) ∗ 𝑓(𝑥𝑀 ) > 0
𝑥 = 𝑔(𝑥) 𝑓(𝑥𝑀 ) > 0; 𝑥𝑀 = 𝑥𝐷
Entonces
Converge hacia la raíz, caso Else
contrario buscar otra función 𝑥𝐷 = 𝑥𝐼 → 𝑓(𝑥𝐷 ) = 𝑓(𝑥𝐼 )
𝑓(𝑥𝑀 ) < 0; 𝑥𝑀 = 𝑥𝐼
𝑔(𝑥) 𝑥𝐼 = 𝑥𝑀 → 𝑓(𝑥𝐼 ) = 𝑓(𝑥𝑀 )
METODO DE LA BISECCION
METODO DE NEWTON RAPSON Si:
Se basa en el promedio que esta
𝑓(𝑥) 𝑓(𝑥𝐷 ) ∗ 𝑓(𝑥𝑀 ) < 0
𝑥𝑖+1 = 𝑥𝑖 − se pueda generar
𝑓 ′ (𝑥)
𝑥𝐷 + 𝑥𝐼 𝑓(𝑥𝐷 )
METODO DE LA SECANTE 𝑥𝑀 = 𝑥𝐷 = 𝑥𝐷 ; 𝑝𝑒𝑟𝑜 𝑓(𝑥𝐷 ) =
2 2
Se asemeja al NEWTON RAPSON 𝑓(𝑥𝐷 ) > 0 𝑥𝐼 = 𝑥𝑀 → 𝑓(𝑥𝐼 ) = 𝑓(𝑥𝑀 )
debido a que la derivada se puede
𝑓(𝑥𝐼 ) < 0 POLINOMIO Y SUS ITERACIONES
expresar de la siguiente forma
Si METODO DE HORNER ITERADO
𝑓(𝑥𝑖 )
𝑥𝑖+1 = 𝑥𝑖 − ′
𝑓 (𝑥𝑖 ) 𝑓(𝑥𝑀 ) > 0; 𝑥𝑀 = 𝑥𝐷 Para un polinomio
𝑓(𝑥𝑖 ) − 𝑓(𝑥𝑖−1 ) Else 𝑝(𝑥) = (𝑥 − 𝑡) ∗ 𝑞(𝑥) + 𝑅
𝑓 ′ (𝑥) =
(𝑥𝑖 − 𝑥𝑖−1 )
𝑓(𝑥𝑀 ) < 0; 𝑥𝑀 = 𝑥𝐼 Haciendo 𝑥 = 𝑡
𝑓(𝑥𝑖 )
𝑥𝑖+1 = 𝑥𝑖 − ACELERACION DE 𝑝(𝑡) = 𝑏0 = 𝑅
𝑓(𝑥𝑖 ) − 𝑓(𝑥𝑖−1 )
(𝑥𝑖 − 𝑥𝑖−1 ) CONVERGENCIA Derivando la función
𝑥𝑖+1 METODO DE AITKEN 𝑝′ (𝑥) = 𝑞(𝑥) + (𝑥 − 𝑡) ∗ 𝑞′(𝑥)
(𝑥𝑖 − 𝑥𝑖−1 )
= 𝑥𝑖 − 𝑓(𝑥𝑖 ) Donde Haciendo 𝑥 = 𝑡
𝑓(𝑥𝑖 ) − 𝑓(𝑥𝑖−1 )
𝑥 = 𝑔(𝑥) 𝑝(𝑡) = 𝑞(𝑡) = 𝑐1
METODO DE LA FALSA POSICION
(𝑥𝑖+1 − 𝑥𝑖 )2 Aplicando NEWTON RAPSON
𝑥𝑖′ = 𝑥𝑖 −
𝑥𝑖+2 − 2𝑥𝑖+1 + 𝑥𝑖
𝑓(𝑡𝑖 ) 𝑏0
Se generan columnas de 𝑥′ y la 𝑡𝑖+1 = 𝑡𝑖 − ′
= 𝑡0 −
𝑓 (𝑡𝑖 ) 𝑐1
segunda evaluación es con la
segunda columna y así con la METODO DE LIN
tercera 𝑓(𝑡𝑖 ) 𝑏0
𝑡𝑖+1 = 𝑡𝑖 − ′
= 𝑡0 −
METODO DE STIFFENSEN 𝑓 (𝑡𝑖 ) 𝑏1
Donde Se estabiliza:
𝑥 = 𝑔(𝑥) 𝑏0
𝑡𝑖+1 = 𝑡0 − 𝜆 ∗
𝑏1
(𝑥1 − 𝑥0 )2
𝑥0′ = 𝑥0 − 𝑓(0)
𝑥2 − 2𝑥1 + 𝑥0 𝜆=−
𝑡0 ∗ 𝑓(𝑡0 )
Donde para la segunda iteración
𝑥3 = 𝑥0
𝑥4 = 𝑥1 ; 𝑥5 = 𝑥2
RAICES COMPLEJAS
METODO DE MULLER
Valores iniciales:
𝑥0 , 𝑥1 𝑦 𝑥2
Evaluar en la función
𝑓(𝑥0 ), 𝑓(𝑥1 ) 𝑦 𝑓(𝑥2 )
Calcular las =ℎ
ℎ0 = 𝑥1 − 𝑥0
ℎ1 = 𝑥2 − 𝑥1
Calcular 𝛿1 𝑦 𝛿2 ,
𝑓(𝑥1 ) − 𝑓(𝑥0 )
𝛿0 =
𝑥1 − 𝑥0
𝑓(𝑥2 ) − 𝑓(𝑥1 )
𝛿1 =
𝑥2 − 𝑥1
Calculo de 𝑎, 𝑏 𝑦 𝑐
𝛿1 − 𝛿0
𝑎=
ℎ1 − ℎ0
𝑏 = 𝑎 ∗ ℎ1 + 𝛿1
𝑐 = 𝑓(𝑥2 )
Cálculo del valor 𝑥3
−2𝑐
𝑥3 = 𝑥2 +
𝑏 ± √𝑏 2 − 4𝑎𝑐
Se evalúa el denominador
|𝑏 + √𝑏 2 − 4𝑎𝑐|
|𝑏 − √𝑏 2 − 4𝑎𝑐|
Se escoge el mayor valor
absoluto