MÉTODOS NUMÉRICOS

Ingeniería Civil

Docente
Dra. María Virginia QUINTANA

MÉTODOS NUMÉRICOS
Técnicas mediante las cuales los modelos matemáticos son
resueltos usando solamente operaciones aritméticas, … para
algunos tediosos cálculos aritméticos.
La computadora es una máquinas que sólo hace lo que se le ordena; y
particularmente a los cálculos numéricos los hace muy rápido y muy bien, sin
fastidiarse.

Por eso, para el ingeniero moderno, los métodos numéricos

implican:
El manejo de hojas de cálculo en PC (Ej. Planilla Excel)
La programación en lenguaje “C++”, “Basic” o “Fortran”
La utilización de software matemático simbólico: Maple, MatLab”, “MathCad”,
“Mathematica”.
La utilización de software especifico: “SAP 2000”, “Plaxis”, “ABAQUS”

02/03/2020 Métodos Numéricos- Ingeniería Civil 2

 Temario…
Algoritmos y
Errores
5 0

2
4

02/03/2020 Métodos Numéricos- Ingeniería Civil 3

TEMA 0
ALGORITMOS Y ERRORES

1. ALGORITMOS
Definición: Un algoritmo es un conjunto finito de instrucciones que especifican una
secuencia de operaciones a realizar en forma ordenada para resolver un problema
específico.

Características: Un algoritmo debe ser simple, claro y preciso para que se pueda
traducir en programas que serán fáciles de seguir y corregir, además debe ser
eficiente para hacer un buen uso de los recursos de memoria y de tiempo de
procesamiento.

02/03/2020 Métodos Numéricos- Ingeniería Civil 4

 Capturado de: [Link]

02/03/2020 Métodos Numéricos- Ingeniería Civil 5

TEMA 0
ALGORITMOS Y ERRORES

2. MÉTODOS NUMÉRICOS

Definición: Son algoritmos que permiten la resolución de problemas cuyas

soluciones consisten en uno o más valores numéricos. Los valores numéricos
obtenidos por este procedimiento constituyen la solución numérica del problema.

Características: Permiten la reformulación de modelos matemáticos, de manera que

puedan ser resueltos en forma aproximada, mediante una cantidad finita de
operaciones matemáticas

02/03/2020 Métodos Numéricos- Ingeniería Civil 6

 Capturado de: [Link]

02/03/2020 Métodos Numéricos- Ingeniería Civil 7

TEMA 0
ALGORITMOS Y ERRORES

2. SEUDOCÓDIGOS
Definición: Es una herramienta que permite describir con claridad los distintos
pasos que componen un algoritmo. Es decir que es un lenguaje intermedio entre el
código de programación de alto nivel y el lenguaje coloquial que empleamos.

Elementos de un seudocódigo:
Datos de entrada: son todos aquellos datos imprescindibles para el desarrollo de
los distintos pasos que componen al seudocódigo.
Proceso: Instrucciones esenciales del algoritmo; indican de una manera clara y
precisa todos los pasos que permiten el objetivo deseado.
Salida de resultados: se refiere a los valores numéricos obtenidos con la
aplicación del algoritmo.

02/03/2020 Métodos Numéricos- Ingeniería Civil 8

 DATOS

CUERPO
PRINCIPAL
(desarrollo
del algoritmo)

SALIDA

02/03/2020 Métodos Numéricos- Ingeniería Civil 9

TEMA 0
ALGORITMOS Y ERRORES

3. PROGRAMACIÓN
Definición: Se denomina programa a un conjunto de instrucciones secuenciales
(comandos o sentencias), interpretables y ejecutables por la computadora.

Implican, por parte del programador, una secuencia de etapas a cumplir en el

tiempo:
 Descripción y análisis del problema a resolver.
 Diseño del algoritmo con Instrucciones específicas
 Planteo del seudocódigo
 Armado del código fuente en algún lenguaje de computadora (ejemplo: C++,
Fortran, Java, etc)
 Ejecución del programa en una computadora.
 Corrección, pruebas y optimización del programa.
 Documentación y mantenimiento del programa.
02/03/2020 Métodos Numéricos- Ingeniería Civil 10
 TEMA 0
ALGORITMOS Y ERRORES

4. TEORIA DE ERRORES
Los errores numéricos se generan al hacer uso de aproximaciones en el
manejo y representación de números y en la resolución aproximada de
modelos matemáticos.

Errores de redondeo: surgen cuando se utilizan aproximaciones en lugar de

números exactos, ya que no es posible trabajar con infinitas cifras decimales
(las computadoras trabajan en promedio con 10 dígitos significativos). El
error es mayor cuando se divide por números cercanos a cero:

02/03/2020 Métodos Numéricos- Ingeniería Civil 11

Errores de truncamiento: aparece cuando un proceso de límites o una serie

infinita son detenidos antes de llegar al valor del buscado. Por ejemplo se
obtiene al considerar la Serie de Taylor con n términos:

Errores de discretización: tienen lugar cuando se reemplaza un modelo

matemático dado, por otro más simple mediante el uso de fórmulas que
discretizan el problema original. En el caso de la definición de derivada:

02/03/2020 Métodos Numéricos- Ingeniería Civil 12

 TEMA 1
RESOLUCIÓN NUMÉRICA DE
PROBLEMAS DE VALORES INICIALES

Clase

 dy
 x   f x , y 

1. Resolución de PVI de primer orden dx

  0

y x  y0


2. Generalización a PVI de orden superior

02/03/2020 Métodos Numéricos- Ingeniería Civil 13

¿Para que estudiamos las ecuaciones diferenciales?

Ley de enfriamiento

02/03/2020 Métodos Numéricos- Ingeniería Civil 14

 ¿Que estudiamos hasta ahora de los PVI o ecuaciones
diferenciales ordinarias de primer orden (EDOs)?

Una serie de herramientas que nos permiten resolver estas

ecuaciones y dar solución a los problemas que describen.

EDO a Variables Separables . . . . . . . . . . . . . . . . . . . . y ' t   f t, y   g t  h y 

Ecuaciones diferenciales homogéneas. . . . . . . . . . . y ' t   f t, y   g 

y
t

EDO Lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . y ' t   P t  y t   Q t 

Ecuación de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . y ' t   P t  y t   Q t  y  t 

EDO Exacta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . M t, y   N t, y  y ' t   0

02/03/2020 Métodos Numéricos- Ingeniería Civil 15

En el mejor de los casos solamente algunas ED se pueden solucionar

analíticamente en forma cerrada. En la mayoría de los problemas de la
ingeniería que implican ED se requiere el uso de métodos numéricos

02/03/2020 Métodos Numéricos- Ingeniería Civil 16

 PVI de 1er orden: Método de Euler
El método de EULER es un método que nos permite encontrar soluciones
numéricas aproximadas para un PVI de primer orden.

dy
 x   f x , y 

dx

y x 0   y 0




02/03/2020 Métodos Numéricos- Ingeniería Civil 17

PVI de 1er orden: Método de Euler

El método de Euler es un método de un solo paso

Objetivo: Solucionar un PVI en forma discreta obteniendo una aproximación a la

solución en un punto a partir del valor de la solución en el punto anterior

Sea un PVI de primer orden:

y ' x   f x , y 
(problema genérico)
y x 0   y 0

Podemos usar el concepto de cociente incremental para aproximar la solución de

la EDO haciendo:
y x  h   y x  y x  h   y x 
y ' x   lim   y x  h   y x   hf x , y 
h 0 h h

Si denotamos por y x  a la solución exacta y por Y x  a la solución aproximada

Y x  h   Y x   hf x ,Y 

02/03/2020 Métodos Numéricos- Ingeniería Civil 18

 PVI de 1er orden: Método de Euler

Curva Sol. Exacta y x 

error= error que se comete por aproximar la derivada con el cociente incremental

02/03/2020 Métodos Numéricos- Ingeniería Civil 19

PVI de 1er orden: Método de Euler

Obtención del algoritmo
Primero dividimos el intervalo x 0, x N  en N subintervalos de ancho h ; o sea:
x  x0
h N (tamaño del paso o paso de integración)
N
Así se obtiene un conjunto discreto de N  1 puntos que cumplen:

Sol. Exacta
x j  x 0  jh, j  0,1,..., N
Sol. Aprox.
Se suele denotar a
x0  a
Cond. Inicial.
xN  b

02/03/2020 Métodos Numéricos- Ingeniería Civil 20

 PVI de 1er orden: Método de Euler
Partimos de la aproximación obtenida por el cociente incremental

02/03/2020 Métodos Numéricos- Ingeniería Civil 21

PVI de 1er orden: Método de Euler

02/03/2020 Métodos Numéricos- Ingeniería Civil 22

 PVI de 1er orden: Método de Euler

02/03/2020 Métodos Numéricos- Ingeniería Civil 23

PVI de 1er orden: Método de Euler

02/03/2020 Métodos Numéricos- Ingeniería Civil 24

Ejemplo
Consideremos el siguiente problema de valor inicial


y '  xy 1  x  2 f x, y   xy

  

y 1 2

La fórmula de Euler para este problema en el intervalo [1, 2], con un paso h = 0,2,
resulta N = 5, entonces es:

 
Y x j 1   Y x j   0.2 f x j ,Y x j  , con 0  j  N  1

La tabla de valores obtenida con la fórmula dada resulta:

j xj Yj
0 1,00 2,0000
1 1,20 2,4000
2 1,40 2,9760
3 1,60 3,8093
4 1,80 5,0282
5 2,00 6,8384

02/03/2020 Métodos Numéricos- Ingeniería Civil 25

Ahora, aplicamos la fórmula el método de Euler con N = 20 y N = 50.

Representando gráficamente los puntos obtenidos y comparándolos con la

solución exacta, dada por la función
t 2 1
y(t )  2e 2

02/03/2020 Métodos Numéricos- Ingeniería Civil 26

 COMENTARIOS
• La solución Y(x) esta tan solo definida en el conjunto discreto de puntos
x , x 
 0 N 

• Un cambio del tamaño del paso h, cambia los puntos de la red (nos sirve
para obtener mas precisión).

• Puede prefijarse el valor del paso h entonces resulta N una función de h.

• El método de Euler rara vez se utiliza en la práctica para obtener la

solución aproximada de un PVI, pero se estudia por su simplicidad en la
derivación de la fórmula y de la determinación del error.

• Los métodos de orden superior utilizan las mismas técnicas, pero el

álgebra que requieren es mucho más complicada

02/03/2020 Métodos Numéricos- Ingeniería Civil 27

Método de Euler – Análisis del error

02/03/2020 Métodos Numéricos- Ingeniería Civil 28

 Errores de los métodos de aproximación

02/03/2020 Métodos Numéricos- Ingeniería Civil 29

Método de Euler – Análisis del error Local

02/03/2020 Métodos Numéricos- Ingeniería Civil 30

 Método de Euler – Análisis del error Local

Método de Euler – Análisis del error Global

E x j ; h 

02/03/2020 Métodos Numéricos- Ingeniería Civil 32

 Consistencia, Estabilidad y Convergencia

Consistencia, Estabilidad y Convergencia

02/03/2020 Métodos Numéricos- Ingeniería Civil 34

 Métodos de Runge-Kutta
Los denominados Métodos de Runge-Kutta (R-K) se caracterizan por tener un
error local de discretización de orden elevado, hecho que les permite tener una
precisión muy alta en los resultados numéricos.

Métodos de Runge-Kutta de segundo orden

02/03/2020 Métodos Numéricos- Ingeniería Civil 36

 Métodos de Runge-Kutta de segundo orden

Solución del sistema

Nota: La elección de distintas combinaciones de los valores de las constantes, que

sean solución, conducen a la ecuación en diferencias de los distintos métodos de
segundo orden (casos particulares)

Métodos de Runge-Kutta de segundo orden

02/03/2020 Métodos Numéricos- Ingeniería Civil 38

 Métodos de Runge-Kutta de segundo orden

Métodos de Runge-Kutta de segundo orden

Interpretación grafica de los métodos de segundo orden
Euler Mejorado: Promedio de las pendientes

02/03/2020 Métodos Numéricos- Ingeniería Civil 40

 Métodos de Runge-Kutta de segundo orden
Interpretación grafica de los métodos
Euler Modificado o del Punto Medio: Pendiente del punto medio del intervalo

02/03/2020 Métodos Numéricos- Ingeniería Civil 41

Métodos de Runge-Kutta de orden superior

02/03/2020 Métodos Numéricos- Ingeniería Civil 42

 Métodos de Runge-Kutta de orden superior

donde

02/03/2020 Métodos Numéricos- Ingeniería Civil 43

Métodos de Runge-Kutta de orden superior

02/03/2020 Métodos Numéricos- Ingeniería Civil 44

 Conclusiones

02/03/2020 Métodos Numéricos- Ingeniería Civil 45

Ejemplos
Obtener una aproximación a la solución para x=0,5 utilizando el Método de
Euler Mejorado, con h=0,1 del siguiente PVI:



y '  2xy

  

y 0 1

Solución:

02/03/2020 Métodos Numéricos- Ingeniería Civil 46