Ing. Luis Manfredo Reyes
En aquellos casos en que no es posible o no es práctico encontrar la solución de una ecuación de una variable por métodos analíticos, se recurre a los métodos numéricos.
Existen excelentes libros que abordan el tema, en especial el de Burden y Faires. El objetivo del presente documento es mostrar la resolución de una ecuación utilizando Microsoft Excel.
METODO DE LA SECANTE
Este método
pretende reducir la cantidad de iteraciones necesarias para lograr la
convergencia de la solución. Se basa en el hecho de que el método de bisección
siempre utiliza la mitad del intervalo, pero no toma en cuenta que la solución
puede estar más cerca de uno de los valores (x0 o x1).
Para ello, se traza una línea entre
f(x0) y f(x1) y se calcula el x2 como el punto en que la línea intersecta
al eje x. (ver gráfica)
El nuevo valor
de x2 se calcula así:
VENTAJAS:
· Converge más rápidamente que bisección
· No necesita derivadas para el cálculo de aproximaciones (comparado con el método de Newton)
DESVENTAJAS:
· En ocasiones, aunque existe la solución, el método no converge
ALGORITMO DEL
METODO DE LA SECANTE:
1. Pedir los
datos de entrada: aproximaciones iniciales P0, P1,
tolerancia T y número máximo de iteraciones N0
2.
Definir: i=2; q0=f(p0); q1=f(p1)
3. Mientras
que I<=N0 seguir los pasos 4-7
4. Calcular p=p1-q1*(p1-p0)/(q1-q0)
5. Si |p-p1|<=T entonces
presentar resultado (p) y fin del proceso, si no:
6 i=i+1
7
Calcular: p0=p1;
q0=q1; p1=p; q1=f(p)
1.
Mensaje de error “El método fracasó despues de N0
iteraciones”
2.
Fin.
IMPLEMENTACION
DEL ALGORITMO:
Si se aplica
el método a una ecuación comparando con el de bisección, se observará que se
necesitan menos iteraciones para obtener la solución final.
El archivo de Excel donde se implementa el método es:
Ejemplo:
Resuelva la
ecuación x5+x-1=0
en la celda E7 se ingresa la tolerancia, que usualmente puede ser 0.00001 o cualquier otro valor
Se deben ingresar dos aproximaciones iniciales para el cálculo, a criterio del usuario
en este caso por ejemplo: 0 y 1
en la celda C10 se ingresa la fórmula, usando en lugar de x el valor B10, es decir
+B10^5+B10-1
Esta celda se copia al resto de la columna.
Excel realiza los cálculos y presenta la solución
RESUELVA LA
ECUACION: x5+X-1=0
|
TOLERANCIA=
|
0.00001
|
||
ITERACION
|
X
|
(FX)
|
MENOR A TOL
|
|
0
|
0
|
-1
|
TODAVIA NO
|
|
1
|
1
|
1
|
TODAVIA NO
|
|
2
|
0.5
|
-0.46875
|
TODAVIA NO
|
|
3
|
0.65957447
|
-0.21559547
|
TODAVIA NO
|
|
4
|
0.79547381
|
0.11398846
|
TODAVIA NO
|
|
5
|
0.74847225
|
-0.01663017
|
TODAVIA NO
|
|
6
|
0.75445642
|
-0.00110441
|
TODAVIA NO
|
|
7
|
0.7548821
|
1.1631E-05
|
TODAVIA NO
|
|
8
|
0.75487766
|
-8.0344E-09
|
SOLUCION
|
|
9
|
0.75487767
|
-5.8398E-14
|
SOLUCION
|
|
10
|
0.75487767
|
0
|
SOLUCION
|
|
11
|
0.75487767
|
0
|
SOLUCION
|
|
12
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
13
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
14
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
15
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
16
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
17
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
18
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
19
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
|
20
|
#¡DIV/0!
|
#¡DIV/0!
|
#¡DIV/0!
|
La solución se encuentra en 8 iteraciones
En forma
comparativa, el método de bisección necesita 21 iteraciones para llegar a la
solución.
No hay comentarios:
Publicar un comentario