METODOS NUMERICOS
CAPITULO 1
SOLUCION DE
ECUACIONES IMPLICITAS
Ing. [Link]. Roberto Parra Zeballos
TEMARIO
Métodos cerrados:
Métodos gráficos
Método de bisección
Método de la posición falsa
Métodos abiertos
Iteración simple de punto fijo
Método de Newton-Raphson
Método de la secante
Raíces de polinomios
Método de Müller
Método de Bairstow
MÉTODOS GRÁFICOS
Los métodos gráficos consisten en graficar la
función f(x) y observar donde la función cruza el
eje x.
Ilustración de las formas generales en que puede
ocurrir una raíz en un intervalo pre-escrito por los
límites inferior xl y superior xu. Las figuras a) y c)
muestran que si f(xl) y f(xu) tienen el mismo signo,
entonces no habrá raíces dentro del intervalo o
habrá un número par de ellas. Las figuras b) y d)
muestran que si la función tiene signos diferentes
en los puntos extremos, entonces habrá un número
impar de raíces dentro del intervalo.
Ilustración de algunas excepciones a
los casos generales. a) Pueden
ocurrir raíces múltiples cuando la
función es tangencial el eje y. En
este caso, aunque los puntos
extremos son de signos opuestos,
hay un número par de intersecciones
con el eje x en el intervalo.
b) Función discontinua donde los
puntos extremos de signo opuesto
contienen un número par de raíces.
Se requiere de estrategias especiales
para determinar las raíces en estos
casos.
EJEMPLO 1
Encontrar la raíz de:
EJEMPLO 2 (CONT.)
Grafica de: f(x) = sen 10x + cos 3x
MÉTODO DE LA BISECCIÓN
Se trata de encontrar los ceros de
f(x) = 0
Donde f es una función continua en [a,b] con f(a) y f(b) con
signos diferentes.
y
f(a)
y = f(x)
b
x
a f(b)
MÉTODO DE LA BISECCIÓN
De acuerdo con el teorema del valor medio, existe p [a,b] tal
que f(p) = 0.
El método consiste en dividir a la mitad el intervalo y localizar
la mitad que contiene a p.
El procesos se repite hasta la lograr la precisión deseada.
MÉTODO DE LA BISECCIÓN
Primera iteración del algoritmo
y Mitad del intervalo que
contiene a p
f(a)
y = f(x)
f(p1)
b
x
a f(b)
p
p1=(a+b)/2
MÉTODO DE LA BISECCIÓN
Segunda iteración del algoritmo
y Mitad del intervalo que
contiene a p
y = f(x)
f(a)
b
x
a =p1 f(b)
f(p2)
p
p2=(a+b)/2
y y
f(a) p=(a+b)/2 f(a) p=(a+b)/2
b b
x
a a
f(b) f(b)
f(a)*f(p) > 0 a = p f(a)*f(p) < 0 b = p
ERROR EN EL MÉTODO DE BISECCIÓN
Criterio de paro, Error aproximado Termina cuando :
RESOLVER LA ECUACIÓN: F(X) = X4-2X3-4X2+4X+4 EN EL INTERVALO [1, 2]
function v=medianas(f,a,b,Tol,Nit)
else
format short
f=inline(f); if f(a)*f(xr)<0
v=[]; b=xr;
xr=(a+b)/2;
i=0;
else
fprintf('\n') a=xr;
fprintf('i a_i b_i xr_i
f(xr_i)\n ') end
v=[v; i a b xr f(xr)]; v=[v; i a b xr
while (i<=Nit) f(xr)];
xr=(a+b)/2;
if(abs(f(xr))<= Tol)
i=i+1;
fprintf('El método tuvo exito end
despues de%3.0f iteraciones',i)
fprintf('\n')
end
fprintf('Resultados obtenidos') if (i == Nit)
fprintf('\n')
fprintf('i a_i b_i
fprintf('El método
xr_i(SOL) f(xr_i)\n ') fracasó despues de %3.0f
v=[v; i a b xr f(xr)]; iteraciones',i)
break end
RESULTADOS DEL PROGRAMA
>> syms x
>> f = x^4-2*x^3-4*x^2+4*x+4;
>> v = medianas(f,1, 2, 0.00001, 50)
i a_i b_i c_i f(c_i)
El método tuvo exito despues de 18 iteraciones
Resultados obtenidos
i a_i b_i c_i(SOL) f(c_i)
v=
0 1.0000 2.0000 1.5000 -0.6875 10.0000 1.4141 1.4146 1.4146 -0.0027
0 1.0000 1.5000 1.5000 -0.6875
11.0000 1.4141 1.4143 1.4143 -0.0007
12.0000 1.4142 1.4143 1.4142 0.0002
1.0000 1.2500 1.5000 1.2500 1.2852
13.0000 1.4142 1.4142 1.4142 -0.0003
2.0000 1.3750 1.5000 1.3750 0.3127 14.0000 1.4142 1.4142 1.4142 -0.0000
3.0000 1.3750 1.4375 1.4375 -0.1865 15.0000 1.4142 1.4142 1.4142 0.0001
4.0000 1.4063 1.4375 1.4063 0.0637 16.0000 1.4142 1.4142 1.4142 0.0000
5.0000 1.4063 1.4219 1.4219 -0.0613 17.0000 1.4142 1.4142 1.4142 0.0000
6.0000 1.4141 1.4219 1.4141 0.0012 18.0000 1.4142 1.4142 1.4142 0.0000
7.0000 1.4141 1.4180 1.4180 -0.0300
8.0000 1.4141 1.4160 1.4160 -0.0144
9.0000 1.4141 1.4150 1.4150 -0.0066
NÚMERO DE ITERACIONES
El error absoluto en la primera iteración es:
El error absoluto en la iteración n-ésima es:
Si el error deseado es Ead, el número de iteraciones
será:
MÉTODO DE FALSA POSICIÓN
Este método considera cual límite del
intervalo está más próximo a la raíz.
f(xu)
De la figura
Despejando
xl xr
xu
f(xl) f(xr)
DESARROLLO DE FALSA POSICION
METODOS ABIERTOS
a) Método cerrado y b) y c)
métodos abiertos.
El método b) diverge y el c)
converge.
FIN 1RA. PARTE
