Hidrico

Descargar como pptx, pdf o txt
Descargar como pptx, pdf o txt
Está en la página 1de 11

MÉTODOS HÍBRIDOS PARA LA SOLUCIÓN

DE RAÍCES DE ECUACIONES
MÉTODOS HÍBRIDOS
Usualmente los métodos iterativos para resolver ecuaciones no lineales (tales como el método
de Newton y el método de la secante), se presentan de manera aislada mientras que en
la práctica computacional, estos métodos nunca se usan de forma aislada (conocidas sus
debilidades), sino acompañados de otros métodos que les agregan robustez, como el método
de bisección. Los paquetes más conocidos implementan, para resolver ecuaciones no lineales
generales, un híbrido conocido como el método de Dekker-Brent que usa secante,
interpolación cuadrática inversa y bisección. En este trabajo se construyen híbridos más
elementales pero robustos: Newton-Bisección y Secante-Bisección y sus respectiva
MÉTODO: SECANTE-BISECCIÓN
Uno de los problemas del método de la secante es que aunque un cero x* de f esté entre las
aproximaciones iniciales x0 y x1, no se garantiza que el método converja (si converge) a x∗,
Para evitar comportamientos no deseables y asegurar que siempre la nueva iteración esté en un
intervalo que contenga al cero que se quiere aproximar, combinamos el método de la secante
con el método de bisección. Para describir el procedimiento usamos int{b,c} para denotar el más
pequeño intervalo cerrado que contiene a b y c, es decir int{b,c} es [c,b] o [b,c].
ALGORITMO: MÉTODO SECANTE-BISECCIÓN
1) Iniciamos con un intervalo [a,b] en el que f(a) y f(b)
tengan signos opuestos. Inicialmente el intervalo de
bisección es [c,b] pues iniciamos con c = a.
2) En el proceso se actualizan a, b y c de la siguiente
manera:
i. Si f(a) y f(b) tienen signos opuestos, se aplica una iteración
del método de la secante.
ii. Si f(a) y f(b) tienen signos iguales, se aplicará una iteración
del método de la secante solo si se conoce que xn+1 estará en
int{c,b}. Sino, se aplica bisección tomando como b el punto
medio del intervalo int{c,b}.
iii. Se actualiza el intervalo int{c,b} actualizando “c” por
comparación del signo de f(c) y f(b).
ALGORITMO: MÉTODO SECANTE-BISECCIÓN
El algoritmo requiere saber a priori si xn+1 ∈ int{c,b} pero sin aplicar la iteración.
Digamos que int{c,b} = [p,q], entonces
MÉTODO: NEWTON-BISECCIÓN
Si el método de Newton falla en una iteración, podemos usar
bisección para dar un pequeño salto y regresar al método de
Newton lo más pronto posible. Supongamos que:
f(a) f(b) < 0. Sea x0 = a o x0 = b. En cada iteración una nueva
aproximación x′ es calculada y a y b son actualizados como sigue:
ALGORITMO:
MÉTODO
NEWTON-
BISECCIÓN
MÉTODO DE BRENT
Es un algoritmo de búsqueda de raíz que combina el método de bisección, el método
de la secante y la interpolación cuadrática inversa.

A) Método de la Bisección:
Supongamos que deseamos resolver la función
f=0, por el método de bisección si recordamos
obtenemos dos valores de a0 y b0, de tal forma
que f tiene signos distintos.
Entonces, si f es continua los valores intermedios
del teorema indica la existencia de una solución
entre a0 y b0.
MÉTODO DE BRENT
B) INTERPOLACIÓN CUADRÁTICA INVERSA:
Cada iteración, tiene tres puntos involucrados:
 b : es la iteración actual -aproximación actual de la raíz
 a : contra punto donde la función tiene signos opuestos.
 b anterior : es la iteración anterior.

Posteriormente se elige valor del nuevo contra punto si la


función tiene signos opuestos a posterior=a de lo contrario a
posterior=b.
MÉTODO DE EULER-LEAL
Emplea dos métodos más simples, intuitivos y fáciles de analizar, los cuales son el
método de Euler, para ecuaciones diferenciales, y el método de la bisección.

A) Ecuación general sin despejar:

Para la cual calcularemos los distintos valores


aproximados de x’ en un intervalo t0 ≤ t ≤ tf , con
la condición inicial x(t0) = x0; la cual, al sustituirla
en la función f:
MÉTODO DE EULER-LEAL
Entonces, dado que no tenemos despejada de manera explícita la variable x’ , para encontrar su
valor, de manera aproximada, usaremos el método de la bisección, donde tomamos dos valores de
arranque inamovibles αi y αd,
αi < αd, tales que

B) ITERACIONES DEL NÚCLEO DEL ALGORITMO:

También podría gustarte