MAT 373

MAT 205

METODOS NUMERICOS

PROGRAMACION EN VBA

Con aplicaciones a la ingenieria

Ing. Leonardo Lopez T.

Prof. [Link] e Industrial
U.A.G.R.M.

SANTA CRUZ 2019

PROLOGO

El presente texto ha sido desarrollado para estudiantes de ingenieria que tengan un

conocimiento inicial de solucionar en forma numerica un problema en computadora. Esto
particularmente como continuacion de la materia Introduccion a La Informatica Basica
(MAT204).

Se ha tenido cuidado particular en la eleccion de los metodos que se describen , que tienen
un alto grado de utilizacion en el campo de la Ingenieria, que muy bien pueden ser
desarrollados con la ayuda de la computacion moderna.

Debe indicarse que los metodos descritos , deben previamente efectuarse un Analisis
Numerico , antes de ser aplicados lo que da como consecuencia una vision clara de cada
metodo , y de esta forma poder comprenderlos en su totalidad.

La combinacion de los metodos y el lenguaje de alto nivel VISUAL BASIC 13 el VBA de

Excel con Macros es con el proposito satisfacer las necesidades de los metodos numericos
en cuanto a calculo ,especialmente cuando son iterativos.
En principio se explica un tema muy importante , como es el de la teoriade errores , por
redondeo en punto flotante , donde se presentara el desarrollo ,para mayor conocimiento
de la evolucion de los mencionados errores.
En general el texto abarca 4 temas importantes que son:

1.-Teoria de errores
2.-Ecuaciones algebraicas no lineales, Sistemas de Ecuaciones Lineales
3.-Interpolacion Derivacion e Integracion
4.-Ecuaciones Diferenciales Ordinarias.

Si bien este texto esta orientado para los universitarios ,puede tambien ser una herramienta
muy util para Ingenieros que se les presente problemas muy puntuales en el ejercicio de su
profesion , y que bajo un previo analisis aplicar un metodo apropiado que le conducira a una
solucion optimizada .
Introduccion al Calculo Numerico
Uno de los principales propositos de las matematicas es el de proporcionar soluciones a los
problemas que se plantean dentro del areas del conocimiento humano.

Particularmente en el campo de la ingenieria la Matematica Tradicional contribuye en gran

medida a solucionar muchos de los problemas planteados , pero no son todos , es entonces
en que se recurre al Calculo Numerico. Se debe tener presente que ambas diciplinas tienen
mucho en comun al dar soluciones , pero por diferentes razones.

Las diferencias que se puedan enumerar entre la Matematica Tradicional y el Calculo

Numerico son:

1.-Grado de precision de las soluciones

2.-Eficacia de la solucion y del metodo.

La solucion exacta solo puede concebirse a travez de la Matematica Tradicional , puesto que
prevee un numero infinito de [Link] debe indicarse que el metodo de resolucion
esta orientado mas que todo a la estructura a seguir , que a conseguir eficacia y fidelidad.

El calculo numerico solo calcula soluciones aproximadas puesto que ejecuta un numero
finito de pasos , a la vez que dicha solucion debe satisfacer directamente al problema
planteado , pero en forma real.

El Calculo Numerico nacio casi junto con la Matematica Tradicional , pero su evolucion y
actualizacion es reciente , debido al Chip de Silicio , puesto que muchos Metodos se han
desarrollado gracias al impulso de la Informatica, los metodos antiguos pueden ahora ser
optimizados.

ANALISIS NUMERICO

CALCULO NUMERICO
METODOS NUMERICOS

METODOS NUMERICOS Y ANALISIS NUMERICO.

Cuando se habla de Analisis Numerico nos referimos a estructuras matematicas que
aparentemente solo se ocupan los Cientificos e Investigadores que necesitan conceptos
avanzados.
Los metodos Numericos tambien podriamos decir que los utilizan Ingenieros , Universitarios
de Tecnologia , profesionales considerados como gente practica

Sin embargo analizando los conceptos anteriormente mensionados se puede cometer un

grave error si se aplica un Metodo Numerico a un determinado problema sin antes haber
hecho un Analisis Numerico .

Los Metodos Numericos sirven para encontrar una solucion inmediata , el Analisis Numerico
busca el mejor metodo para aplicarlo , indicando el grado de precision de la respuesta o
solucion.
Entonces idealmente un metodo debe ser aplicado para una solucion que tenga el menor
porcentaje de error.

Los metodos numericos se analizaran y luego se diseñaran algoritmos con diagramas de flujo,
y luego se codificran en VB_13 y VBA Excel. Los algoritmos en su mayoria son del tipo
ITERATIVO ademas los calculos aritmeticos internos se relizaran en el sistema de numeros
binarios , los errores por redondeo pueden acumularse en porcentajes considerables , es por
esta razon que debe realizarse el menor numero de operaciones , y si un determinado metodo
cumple con esta condicion ; estamos hablando entonces de un Metodo Numerico eficiente.

TEORIA DE ERRORES EN CALCULO NUMERICO

Los datos que son el objeto principal en los procesos de calculo,pueden tener diferentes
tipos de errores al igual que los resultados.

Estos errores podemos clasificar de acuerdo a su naturaleza, es decir

-Errores por sus Caracteristicas

-Errores por su procedencia.

Errores por sus Caracteristicas.-Estos son:

Sea X numero exacto
X* numero aproximado
|
Ea= |X - X* |=Error Absoluto**
Er= |1 - X* / X | =Error Relativo
E%= |1 - X* / X|*100= Error Porcentual

Debe indicarse que el error relativo se calcula en base al valor exacto, valor que casi
nunca tendremos, pero al utilizarse tecnicas iterativas que , en base a cada repeticion
se mejora la aproximacion al valor exacto se toma como referencia , esa ultima
aproximacion como valor exacto.
Errores por su Procedencia.-

-Inherentes
-Redondeo
-Truncamiento

Los errores inherentes, son aquellos, que los datos, antes de ser utilizados ya lo
tienen, estos pueden ser por errores humanos,deficiencia de los instrumentos de
medicion, resultado de experimentos de laboratorio, o finalmente producto de un
proceso de datos anterior.

Los errores por Redondeo y Truncamiento se producen en el momento del

tratamiento de datos, tiempo de [Link] entender la evolucion de este tipo de
errores , es conveniente analizar las caractristicas del medio de calculo , esn este
caso la forma de almacenar los numeros de un computador, y esto se realiza a travez
de la aritmetica de punto flotante.(APF).

ARITMETICA DE PUNTO FLOTANTE

La capacidad de almacenamiento de los computadores es limitada, por lo tanto a veces

es necesario trabajar con un numero fijo de cifras decimales, o bien con un numero
fijo de cifras significativas esto hace que en muchos momentos del calculo los
numeros tengan que ser redondeados eso proboca que muchos digitos puedan ser
eliminados produciendo un error llamado un error por redondeo.

El e rror por redondeo se ve en toda operación aritmetica, de punto flotante,

este significa que constantemente en cada operación el punto decimal se va moviendo.

Ej: Sean los numeros.

Formato Usual Formato APF Notacion Cientifica
73.465 0.73465x 102 0.73465 E+2
495.39 0.49539 x 103 0.49539 E+3
0.00657 0.657 x 10-2 0.657 E-2

Recordemos que cada número esta representado por una fracción al que lo
llamaremos mantiza, la cual esta multiplicada por una potencia del numero de base 10
y al que llamaremos exponente.

Entonces cualquier numero X podemos representarlo de la sgt. Manera

X=ƒ*10e
 donde “e” significa el exponente que esta en función del numero de digitos
recorrido por el punto flotante, si es hacia la izquierda es positivo, si es a la derecha
negativo.
El resultado de una operación aritmetica consta en general de dos partes.

1. La ubicación del punto flotante

2. El redondeo de acuerdo al numero de digitos significativos.

Cuando se suman por ejemplo dos numeros el de manera exponente debe

acomodarse de acuerdo al mayor.

Ej:
0.1964 x 103
+ 0.1263 x 101
-------------------
0.1964 x 103
+ 0.001263 x 103
------------------
0.197663 x 10 3=0.1976 x 103 + 0.63 x 10-1

Del ultimo ejemplo podemos representar una forma general en el cual el

resultado de una operación nos dice lo siguiente.

Y=ƒy x 10e + gy x 10 e-t

ƒy debe ser menor que uno y mayor o igual a 1/10

|ƒy| < 1 y |ƒy| >=1/10

El caso de gy es diferente puesto que no podemos garantizar que este

normalizado, pero podemos dar los sgtes. Limites.

0<=gy< 1
 Existen 2 tipos de Redondeo:
Truncado
Simetrico

El redondeo truncado no tiene en cuenta gy es decir lo elimina completamente sin

importar su valor, en este caso para evaluar el maximo error relativo debido al redondeo
truncado se debera tomar como valor maximo gy= 1 y fy=0.1 .
Teniendo entonces la sgt. Operación.

|ey |/|y|=gy x 10e-t/fy x 10e = 1 x 10e-t / 0.1 x 10e = 1 x 10-t+1

Si recordamos que “t” es el numero de digitos de la mantiza de cualquier
numero de punto flotante, obtenemos un resultado interesante, el maximo error
relativo por redondeo de una operación aritmetica no depende del tamaño de las
cantidades, sino del numero de digitos significativos.

El tipo mas conocido de redondeo se lo denomina redondeo simetrico y que

tiene tres reglas importantes.

REGLAS PARA EL REDONDEO SIMETRICO

1. Todos los digitios despues del ultimo digito significativo permanecen iguales
si el proximo digito es menor que 5.
2. Se añade un digito al ultimo digito significativo si el proximo digito es mayor
que 5
3. Si el primer digito suprimido es igual a 5 se aumenta una unidad a la ultima
cifra conservada si esta es impar y no se modifica si es par.

En caso de las computadoras estandar la regla 3 se agrega a

la regla 2.

Al hablar de aritmetica de punto flotante nos referimos basicamente al punto

decimal que constantemente recorre hacia la izquierda haciendo siempre fluctua entre
cero y uno, esto tiene una norma ventaja a uno esto puede ser posible utilizando
tambien notación cientifica.

El redondeo simetrico es mas facil de encontrar el error absoluto, por ejemplo.

Si gy <1/2 |ey|=gy x 10e-t

 Si gy >1/2 |ey|=|1 - gy| x 10e-t
De cualquier manera en ambos casos siempre se multiplicará un factor que
nunca es mayor que ½ que por lo tanto el maximo error relativo absoluto sera:

max y  0.5 • 10e −t

Si calculamos el maximo error relativo absoluto tendremos

Abs(ey/y)=(g(x)*10e-t)/(f(x)*10e)=(0.5*10e-t/0.1*10e)=5*10-t = 0.5*10-t+1

El ultimo resultado sera el maximo error relativo absoluto producido por el

redondeo Simetrico y que es la mitad del Redondeo Truncado , por lo tanto es mas
eficiente el redondeo simetrico.

PROPAGACION DEL ERROR POR REDONDEO

Es de mucha importancia especialmente en analisis numerico la forma en como

se propaga un error a traves de la computadora, es decir si su efecto aumenta o
disminuye al realizarse las diferentes operaciones aritmeticas.

Como primera medida debemos encontrar expresiones para el error absoluto y

el error relativo en cada una de las 4 operaciones aritmeticas, en función de sus
operadores y sus errores. Posteriormente desarrollaremos un metodo que nos permita
determinar un limite a los errores en el total de los calculos.

SUMA

X,Y son numeros exactos, son los numeros aproximados.

x = x* + ex
y = y* + ey
x + y = x* + ex + y* + ey
x + y = (ex + ey) + x* + y*
ex+y = ex + ey
Resta
ex-y = ex - ey

Multiplicación
ex*y = x* ey + y* ex

División
ex/y = ex /y* - x* ey/y*2

Con estas consideraciones podemos formular una tabla para las cuatro
operaciones aritmetica, esta tabla calculará los errores relativos, con todas las
consideraciones que se tomo al principio, tendremos los sgt.

Suma
e(x+y/x*+y*) = (x*/(x*+y*))(ex /x*) + (y*/(x*+y*))(ey/y*/)

Resta
e(x-y/x*-y*) = (x*/(x*+y*))(ex /x*) - (y*/(x*+y*))(ey/y*/)

Multiplicación
ex.y/x*.y* = ex /x* + ey/y*
División
ex/y/x*/y* = ex /x* - ey/y*

Es importante comprender claramente el significado de cada una de estas

formulas, los errores pueden ser de cualquier tipo, es mas los valores aproximados X
y Y pueden ser resultado de datos experimentales que contienen errores inherentes,
las formulas nos permiten calcular los errores por redondeo a partir de los datos
iniciales, por ejemplo podemos averiguar la evolución de los errores en el sgte calculo.

U = ( x+y ) * z

e x + y =  x +x y  *  exx  +  x +y y  *  eyx  + r1
x+ y

       

1
e x + y = 5 *10
x+ y −t
=
2
*10− t +1

Cuando se trata de la segunda operación el error relativo de la multiplicación

debe venir acompañado por un error llamado de redondeo, es decir:

 ey  e
e u = x +x y *  exx  + x +y y * 
u + z +r
y  z m
   

Podemos hacer comparaciones entre cada una de las operaciones pero lo mas
importante es que podemos evaluar el maximo error relativo.

e 
e u =  xx++ yy  + ezz + rm
u

 

e u = 5 *10
u −t
+ 5 * 10 − t

e u  1*10
u − t +1

Ejemplo Nro 1: Sumar las siguientes cantidades manualmente .

a) De menor a mayor
b) De mayor a menor

1. 0.2897*100
2. 0.4976*100
3. 0.2488*101
4. 0.7259*101
5. 0.1638*102
 6. 0.6249*102
7. 0.2162*103
8. 0.5233*103
9. 0.1403*104
10. 0.5291*104

a) Sumando de menor a mayor

0.1053x102
0.2897x100 +0.1638x102
+ 0.4976x100 0.2691x102
0.7873x100 +0.6249x102
0.07873x101 0.8940x102
+0.2488 x101 +0.2162x103
0.32753x101 0.08940 x103
0.3275 x 101 +0.2162 x103
+0.7259 x 101 0.3056 x103
1.0534x101 +0.5233 x103
0.10534x102 0.8289 x103
0.1053 x102 + 0.1403 x104
0.08289 x104
+ 0.1403 x104
0.22319 x104
0.2232 x104
0.2232 x 104
+ 0.5291 x 104
0.7523 x 104
b) Sumando de Mayor a Menor(Desarrollo Idem Inciso a))

0.6694 x 104
0.7217 x 104
0.7473 x 104
0.7495 x104
0.7511 x104
0.7518 x104
 12

0.7520 x104
0.7520 x104
0.7520 x104
Valor exacto=0.7522904 x104
Del ejemplo anterior podemos indicar que el porcentaje de arrastre entre una
suma ascendente o descendente es de gran diferencia, notandose en comparación con
la cifra exacta el porcentaje es 30 veces mayor, por lo tanto se recomienda sumar
siempre en forma ascendente.

Ejemplo Nro 2: sumar 4 numeros muy proximamente iguales, de acuerdo a

las sgts. Condiciones

x1,x2,x3,x4
Si i=1,2,3,4
xi = x0 +i
Si  i x0 x 10
y=x1+x2+x3+x4
Si X1=0.5243 x 100: X2=0.5262 x 100: X3=0.5226 x100: X4=0.5278 x 100

Sumando de forma secuencial se obtiene redondeando despues de cad

operación(Redondeo Simetrico), se obtiene :

Y=0.2102 x101

Si agrupamos Y=(X1+X2)+(X3+X4)

Y=0.2101 x101
El valor exacto es Y=0.21009 x 101 el valor mas proximo es sumando agrupando.

Regla.-Cuando se trata de sumas de “n” numeros muy proximamente iguales, el error

por redondeo se reduce si se tiene n2 numeros, agrupando en “n” grupos de “n”
elementos cada uno,de tal manera que despues se suman las “n” sumas parciales.
 13

Ejemplo Nro3.: En este ejemplo se ejecutará la resta de dos numeros, muy

proximamente iguales, de acuerdo al sgte esquema

Sea x,y muy proximo iguales

Z= x-y
𝑒𝑧 𝑋 𝑌
=( )𝑖𝑥 − ( ) 𝑖𝑦 + 𝑟𝑚
𝑧 𝑋−𝑌 𝑋−𝑌

Si los errores inherentes son por ejemplo numeros perfectamente redondeados

entonces:

ex <= 0.5 x 10−𝑡+1 ey<=0.5 x 10−𝑡+1

𝑔𝑥 = 0.5 ∗ 10𝑒−𝑡 𝑔𝑦 = 0.5 ∗ 10𝑒−𝑡

Consideramos el problema utilizando numeros reales por ejemplo:

4
x = 0.5628*10
4
y = 0.5631*10

ex 0.5 ∗ 10e−t
=
x fx ∗ 10e
ex 0.5∗100
= ≤ 0.01%
x 0.5628∗104

ey 0.5 ∗ 10e−t
=
y fy ∗ 10e
ey 0.5 ∗ 100
= ≤ 0.01%
y 0.5631 ∗ 104

Z = -0.0003*104= -3
 14

ez/z= 𝑎𝑏𝑠 ((0.5(+3−)0.5)) = 0.3333 = 33%

Los errores relativos individuales son pequeños como se ve
en lineas arriba, pero el error producido como resultado de la
operación resta es my grande casi el 33%.

LISTA DE RECOMENDACIONES PARA LOGRAR MAYOR PRECISION

Podemos hacer un resumen de los ejemplos anteriores de acuerdo a lo sgte.:

1. Cuando se van a sumar o restar tarbajar primero con los numeros mas
pequeños
2. De ser posible evitar siempre las sutracción de dos numeros muy
proximamente iguales.
3. Una expresión del tipo a (b-c) puede ser escritas = ab-ac tambien la expresión
(a-b)/c=a/c-b/c.

Si hay numeros muy aproximadamente iguales en el parentesis se debe

ejecutar la resta antes que la multiplicación esto evitará errores por redondeo
adicionales.
4. Si no se puede cumplir las anteriores reglas debemos de tratar de simplificar
el maximo las operaciones aritmeticas utilizando criterios que permitan en el menor
numero de calculos reslutados.
15