SOLUCIONES DE ECUACIONES DE UNA VARIABLE
Existen diversos métodos para obtener la raíz o solución aproximada de una
ecuación de una variable. Los métodos que estudiaremos, en términos generales,
solamente proporcionan una raíz a la vez y en nuestro caso también sólo nos
interesan las raíces reales. A continuación, estudiaremos algunos métodos.
• MÉTODO DE BISECCIÓN
Este método se utiliza para obtener la raíz aproximada de una ecuación de la
forma f(x) = 0. Si la función f(x) es continua en el intervalo [a, b] y si además f(a) y
f(b) poseen signos diferentes, es decir f(a) * f(b) < 0, entonces dicha función poseerá
un cero o raíz en dicho intervalo. Este método se basa en el teorema del valor
intermedio, el cual establece que si existe un número c ]a, b[ y si además
existe un número k de tal manera que se cumpla: f(a) < k < f(b) , entonces se
cumple: f(c) = k.
Por lo tanto, para aplicar el método de bisección deberemos verificar que se
cumpla: f(a) * f(b) < 0.
Luego para obtener el primer valor de aproximación se utilizará la siguiente
fórmula de aproximación:
𝒃𝒊 − 𝒂𝒊 𝒂𝒊 + 𝒃𝒊
𝑷𝒊 = 𝒂𝒊 + =
𝟐 𝟐
Si a1 = a ^ b1 = b
𝒃𝟏 − 𝒂𝟏 𝒂+𝒃
𝑷𝟏 = 𝒂𝟏 + =
𝟐 𝟐
Para saber si dicho valor obtenido es el deseado, deberemos averiguar el error de
aproximación de la siguiente manera:
𝑬𝒓𝒓𝒐𝒓 = |𝒇(𝑷𝟏) | < 𝜖
Si el criterio de comparación no se satisface, entonces deberemos averiguar un
nuevo valor de aproximación, para ello es necesario establecer un nuevo
intervalo de la siguiente manera:
𝑺𝒊 𝒔𝒆 𝒄𝒖𝒎𝒑𝒍𝒆: 𝒇(𝑷𝟏) ∗ 𝒇(𝒂) < 0 , 𝒍𝒂 𝒓𝒂í𝒛 𝒆𝒔𝒕á 𝒆𝒏 [𝒂𝟐 , 𝑷𝟏 ]
a2 = a1 ^ b2 = P 1
𝒃𝟐 − 𝒂𝟐 𝒂 + 𝑷𝟏
𝑷𝟐 = 𝒂𝟐 + =
𝟐 𝟐
𝑬𝒓𝒓𝒐𝒓 = |𝑷𝟐 − 𝑷𝟏 | < 𝜖
𝑺𝒊 𝒔𝒆 𝒄𝒖𝒎𝒑𝒍𝒆: 𝒇(𝑷𝟏) ∗ 𝒇(𝒃) < 0 , 𝒍𝒂 𝒓𝒂í𝒛 𝒆𝒔𝒕á 𝒆𝒏 [𝑷𝟏 , 𝒃𝟐 ]
a2 = P1 ^ b2 = b 1
𝒃𝟐 − 𝒂𝟐 𝑷𝟏 + 𝒃
𝑷𝟐 = 𝒂𝟐 + =
𝟐 𝟐
𝑬𝒓𝒓𝒐𝒓 = |𝑷𝟐 − 𝑷𝟏 | < 𝜖
Esto deberá realizarse sucesivamente hasta que el criterio de paro se cumpla. Es
conveniente elegir un intervalo inicial pequeño.
Ejemplo:
La velocidad vertical de un cohete se calcula mediante la siguiente fórmula:
𝐦𝐨
𝐯 = 𝐮𝐥𝐧 ( ) − 𝐠𝐭
𝐦𝐨 − 𝐪𝐭
donde:
g = 9.81 m/seg2
q = tasa de consumo de combustible = 2630 kg/seg
u= velocidad con la que se expele el combustible = 7920km/h
mo = masa inicial del cohete = 145000 kg
Emplee el MÉTODO DE BISECCIÓN para determinar el tiempo “t”, para el cual el
cohete alcanza una velocidad de 985 m/s, con una precisión de 10 -12. MUESTRE
EN FORMA DE TABLA: #ITERACIONES, t0, t1, VALOR APROXIMADO, ERROR.
EMPLEE QUINCE DECIMALES.
Solución:
Reescribir la función de trabajo de la forma f(t) = 0
𝐦𝐨
𝐯 = 𝐮𝐥𝐧 ( ) − 𝐠𝐭
𝐦𝐨 − 𝐪𝐭
𝐦𝐨
𝐮𝐥𝐧 ( ) − 𝐠𝐭 − 𝐯 = 𝟎
𝐦𝐨 − 𝐪𝐭
𝐦𝐨
𝒇(𝒕) = 𝐮𝐥𝐧 ( ) − 𝐠𝐭 − 𝐯
𝐦𝐨 − 𝐪𝐭
Para elegir un intervalo de análisis podemos graficar la función f(t), sustituyendo
previamente los valores de ‘u’, ‘q’, ‘mo’, ‘g’, ‘v’, y luego elegimos un intervalo que
contenga la raíz buscada. Los resultados se muestran a continuación:
Método de la bisección
introduzca el valor de la velocidad v: 985
introduzca el valor de la gravedad g: 9.81
introduzca el valor de masa mo: 145000
introduzca el valor de la velocidad del combustible u: 2200
introduzca el valor de la tasa de consumo de combustible q: 2630
introduzca el valor de to: 23.38816058
introduzca el valor de t1: 23.38816059
introduzca la función: u*log(mo/(mo-q*t))-g*t-v
introduzca el valor de precisión:1e-12
n to t1 t error
1 23.388160580000001 23.388160589999998 23.388160585000001 1.31e-07
2 23.388160585000001 23.388160589999998 23.388160587500000 2.50e-09
3 23.388160585000001 23.388160587500000 23.388160586250002 1.25e-09
4 23.388160586250002 23.388160587500000 23.388160586875003 6.25e-10
5 23.388160586875003 23.388160587500000 23.388160587187500 3.12e-10
6 23.388160587187500 23.388160587500000 23.388160587343748 1.56e-10
7 23.388160587187500 23.388160587343748 23.388160587265624 7.81e-11
8 23.388160587187500 23.388160587265624 23.388160587226562 3.91e-11
9 23.388160587187500 23.388160587226562 23.388160587207032 1.95e-11
10 23.388160587187500 23.388160587207032 23.388160587197266 9.77e-12
11 23.388160587197266 23.388160587207032 23.388160587202151 4.88e-12
12 23.388160587197266 23.388160587202151 23.388160587199707 2.44e-12
13 23.388160587197266 23.388160587199707 23.388160587198485 1.22e-12
14 23.388160587197266 23.388160587198485 23.388160587197873 6.11e-13
El valor aproximado del tiempo es: t= 23.388160587197873
• MÉTODO DE ITERACIÓN DE PUNTO FIJO
Este método se utiliza para obtener la raíz aproximada de una ecuación de la
forma f(x) = 0, pero en el proceso de cálculo no se utiliza directamente f(x).
Una función g(x) poseerá un punto fijo en x = p, si se cumple que: g(p) = p.
Si la función g(x) posee un punto fijo en x = p, entonces la función f(x) de la forma:
f(x) = x - g(x) poseerá un cero o raíz en x = p.
Para aplicar el método de iteración de punto fijo, deberá utilizarse una función
g(x).
A partir de: f(x) = 0, se establece una función g(x) de tal manera que se cumpla:
x = g(x)
Dicha función g(x) debe satisfacer lo siguiente:
a) Si la función g es continua en [a , b] y además g(x) [a , b], x [a , b],
entonces poseerá un punto fijo en [a , b]
b) Si además g’(x) está definida en ]a , b[ y si existe una constante positiva K <1,
tal que se cumpla | g’(x) | K, x ]a , b[, entonces el punto fijo es único.
De las condiciones anteriores, solamente la primera nos interesa aplicar.
Para iniciar el proceso de obtención de la raíz aproximada, una vez establecida la
función g(x), deberemos seleccionar un valor inicial.
Si x = p, entonces la fórmula de aproximación viene dada así:
p = g(p)
Si po es el punto de inicio, el primer valor aproximado se obtendrá así:
𝑷 = 𝒈(𝑷𝒐)
Para saber si dicho valor obtenido es el deseado, deberemos averiguar el error de
aproximación de la siguiente manera:
𝑬𝒓𝒓𝒐𝒓 = |𝑷 − 𝑷𝒐 | < 𝜺
Si el criterio de comparación no se satisface, entonces deberemos averiguar un
nuevo valor de aproximación, para ello es necesario un nuevo valor inicial.
Luego para la siguiente iteración, el nuevo valor inicial se obtendrá así:
Po = P
Es decir, la raíz obtenida en la iteración previa se convierte en valor inicial para
la siguiente y así sucesivamente hasta que el criterio de paro se cumpla.
Con el nuevo valor inicial se procede a averiguar la nueva raíz y el nuevo error.
Ejemplo:
Emplee el método de iteración de punto fijo para obtener el valor de x que
satisfaga la siguiente ecuación: x4 – 3x2 – 3 = 0, con una precisión: є = 10-5, en el
intervalo [ 1, 2 ]. MUESTRE EN FORMA DE TABLA: #ITERACIONES, p0, VALOR
APROXIMADO, ERROR. EMPLEE QUINCE DECIMALES.
Solución:
A partir de la ecuación: x4 – 3x2 – 3 = 0, deberemos generar una función g(x) que
cumpla las condiciones del punto fijo.
Estableciendo la siguiente función:
𝟒
𝒈(𝒙) = √𝟑 + 𝟑𝒙𝟐
Para averiguar si dicha función cumple la condición a) de punto fijo, graficamos
dicha función g(x) en el intervalo dado, si no cumpliese, entonces deberá
establecerse otra función hasta que cumpla la condición a) de punto fijo.
Una vez establecida la función g(x) adecuada, debemos elegir un punto fijo dentro
del intervalo y luego comenzamos el proceso de cálculo del valor aproximado. Los
resultados se muestran a continuación:
Ejercicios
La distancia D(t) recorrida por un automóvil se establece mediante la siguiente
ecuación: D(t) = –70 + 7t + 70e-t/10. Aproxime el valor de “t” para el cual la
distancia es igual a: 37.855364040665670.
Emplee el MÉTODO DE BISECCIÓN para determinar el tiempo, con una exactitud
-10
de 10 . MUESTRE EN FORMA DE TABLA: #ITERACIONES, t0, t1, VALOR
APROXIMADO, ERROR. EMPLEE QUINCE DECIMALES.
Emplee el método de iteración de punto fijo para obtener el valor de x que
satisfaga la siguiente ecuación: x4 + 2x2 – x – 3 = 0, con una precisión: є = 10-5,
en el intervalo [ –0.88, –0.87]. MUESTRE EN FORMA DE TABLA: #ITERACIONES,
p0, VALOR APROXIMADO, ERROR. EMPLEE QUINCE DECIMALES.
