Nombre del alumno:
Juan Pablo Guzmán Arcos.
Catedrático:
Ángel Jiménez Vásquez.
Carrera:
Ing. Civil.
Semestre:
4* “C”
Materia:
Métodos numéricos.
Tema:
UNIDAD II: Solución de ecuaciones no lineales de una
variable.
INDICE.
2.1------------ Búsqueda de valores iniciales tabulación y
graficación.
2.2------------ Métodos cerrados y sus interpretaciones
geométricas (bisección y regla falsa).
2.3------------ Métodos abiertos y sus interpretaciones
geométricas como sus criterios de convergencia (Newton
Raphson, secante).
2.4------------ aplicación de la solución de ecuaciones no
lineales.
2.5------------ Uso de herramientas computacionales.
2.1- Búsqueda de valores iniciales, tabulación y
graficacion.
Este paso a mi parecer es de los mas importante, en la clase
lo obteníamos derivando la función dada y dando un valor a
esta para ver si cumple con los requisitos dicho valor (y así
tomando como valor inicial).
El uso de cualquier logaritmo numérico para encontrar las
raíces de f(x)=0 requiere uno o mas valores iniciales; además,
en métodos como la bisección y el de posición falsa, los dos
valores iniciales requeridos deben estar a los lados de la raíz
buscada, y sus valores funcionales correspondientes tienen
que ser de signos opuestos.
A continuación se dan algunos alineamientos generales para
obtener valores aproximados a las raíces de f(x)=0.
Por lo general, la ecuación cuya raíces se buscan tienen
algún significado físico; entonces, a partir de
consideraciones físicas pueden estimarse valores
aproximados a las raíces. Este razonamiento es
particular para cada ecuación.
Otra manera de conseguir información sobre la función,
que permita determinar valores iniciales “adecuados”,
consiste en obtener su grafica aproximada mediante un
análisis de f(x), a la manera clásica del cálculo
diferencial e integral, o bien como se ha venido
sugiriendo, con algún software comercial o, en el mejor
de los casos, empleando ambos.
Tabulación y graficación.
Tabulación se refiere al hecho de calcular valores parciales
para una función y compararlos en una tabla, de ahí el
nombre de tabular.
El método general para graficar cualquier función es el de
tabulación. Consiste en dar valores a la variable x y con ellos
calcular los correspondientes a la variable y, los cuales se
van anotando en una tabla.
Después se localiza en el plano cartesiano cada punto
tabulado así y se unen para obtener la forma de la gráfica
buscada.
Por ejemplo, para graficar y= -2x-1, dando valores a la x de
-2, -1, 0, 1, 2 y 3 se construye la siguiente tabla
X -2 -1 0 1 2 3
y -5 -3 -1 1 3 5
Y llevando esos puntos al plano cartesiano se obtiene la recta
de la figura siguiente:
De igual forma, para graficar y=x´2-
10x + 24, dando valores a la x, por
ejemplo de x=2 se obtiene para la
y=(2)´2-10 (2) + 24 y= 8.
Repitiendo el procedimiento para los valores de x de:
2,3,4,5,6,7 y 8 y concentrando los valores en una tabla:
X 2 3 4 5 6 7 8
y 8 3 0 -1 0 3 8
Estos puntos localizados en el plano cartesiano dan la
siguiente figura:
Y uniéndolos se llega a la parábola que se muestra en la
figura siguiente:
2.2- Métodos cerrados y su interpretación geométrica
(bisección y regla falsa).
Conocido también como de corte binario, de participación en
dos intervalos iguales o método Bolzano, es un método de
búsqueda incremental en el que el intervalo se divide siempre
en 2. Si la función cambia de signo sobre un intervalo, se
evalúa el valor de la función en el punto medio. La posición
de la raíz se determina situándola en el punto medio del
subintervalo dentro del cual ocurre un cambio de signo. El
proceso se repite hasta obtener una mejor aproximación.
Pasos, métodos de bisección.
Paso 1.- elegir los valores iniciales interior XI y superior Xu de
forma que la función cambia de signo sobre el intervalo. Esto
se puede verificar asegurándose de que f(XI) f(Xu) < 0.
Paso 2.- la primera aproximación a la raíz, se determina
como:
Paso3.- realice las siguientes evaluaciones para determinar
en que subintervalo cae la raíz.
a) Si f(XI) f(Xu) < 0 entonces la raíz se encuentra dentro del
subintervalo . por lo tanto, tome Xu=XI y continúe en el
paso 2.
b) Si f(XI) f(Xu) > 0 entonces la raíz se encuentra dentro
del subintervalo superior. Por lo tanto resuélvase XI=Xr y
continúe en el paso 2.
c) Si f(XI) f(Xu) = 0. La raíz es igual a Xr; determina el
calculo.
Ejemplo:
Use la aproximación grafica para determinar el coeficiente de
razonamiento c necesario para que una paracaidista de masa
m=68.1 kg tenga una velocidad de 40 m/s después de una
caída libre de t=10 s. La aceleración de la gravedad es de 9.8
m/s´2.
Varios valores de c pueden ser sustituidos en el lado derecho
de esta ecuación para calcular
c f(c)
4 34.115
8 17.653
12 6.067
14 1.569
16 -2.269
20 -8.401
La aproximación grafica para determinar las raíces de una
ecuación.
Solución:
12+16
Xr= =14
2
En el siguiente calculo el producto del valor de la función de
un limite inferior y un punto medio es:
F(12) f(14)= 6.067 (1.569)= 9.517
La raíz debe estar localizada entre 14 y 16
14+ 16
Xr= =15
2
∈t- 1.5%
F(14) f(16) – 1.569 (0.425) – 0.666
Por lo tanto la raíz esta localizada entre 14 y 15.
14+ 15
Xr= =14.5
2
∈t= 1.9 %
|Ea|= |15−14
15 |
∗100 %=6.667 %
Iteració Xf Xu Xr Ea % Es %
n
1 12 16 14 5.279
2 14 16 15 6.66 1.487
7
3 14 15 14.5 3.44 1.896
8
4 14.5 15 14.75 1.69 0.204
5
5 14.7 15 14.875 0.84 0.641
5 0
6 14.7 14.875 14.8125 0.42 0.219
5 2
-Regla falsa.
La falsa posición es una alternativa basada en una
visualización grafica.
Este método alternativo que aprovecha la visualización
grafica es la de unir f(x) y f(Xu) con una línea recta. La
intersección de esta línea en e eje de las x re presenta una
mejor estimación de la raíz. El echo de que se remplace la
curva por una línea recta de una “posición falsa” de la raíz.
También a este método se le conoce como interpolación
lineal.
Usando triángulos semejantes la intersección de la recta en el
eje x puede ser estimado como:
f ( XI ) f ( Xu)
=
Xr−XI Xr−Xu
Por lo tanto
f ( Xu)(XI −Xu)
Xr= Xu−
f ( XI ) −f ( Xu)
Esta es la formula de la falsa posición.
Grafica método de la regla falsa
2.2.- Métodos abiertos y sus interpretaciones geométricas
como sus criterios de convergencia (Newton Raphson,
secante).
Los métodos abiertos utilizan una formula para predecir la
raíz. Esta fórmula puede desarrollarse como una iteración de
un punto fijo (también llamada iteración de un punto o
sustitución sucesiva o método de punto fijo).
- Newton-Raphson.
Entre los métodos de aproximaciones sucesivas para
encontrar alguna de las raíces de una ecuación algebraica o
trascendente, el cual de Newton-Raphson es el que presenta
mejores características de eficiencia, debido a que casi
siempre converge a la solución y lo hace en un numero
reducido de iteraciones.
Este método es aplicable tanto en ecuaciones algebraicas
como trascendentes y con el es posible obtener raíces
complejas. Tal vez, de las formulas de Newton-Raphson sea
lamas ampliamente utilizada.
Si el valor inicial para la raíz xi, entonces se puede trzar una
tangente desde el punto [ xi , f (xi) ] de la curva. Por lo común, el
punto donde esta tangente cruza el eje x representa una
aproximación mejorada de la raíz.
El método de Newton-Raphson se deduce a partir de esta
interpretación geométrica.
El método de Newton-Raphson, como todos los de
aproximaciones sucesivas, parte de una primera
aproximación de una formula de recurrencia se acercara a la
raíz buscada, de tal manera que la nueva aproximación se
localiza en la intersección de la tangente a la curva de la
función en el punto y el eje de las abscisas.
De la figura se tiene que la primer derivada en x es
equivalente a la pendiente:
f ´ ( x 1 )= f ´ ¿ ¿
¿
Que se reordena para obtener :
f ( x 1)
x Y −1 =x1 −
f (x Y )
La cual se conoce como la formula de Newton-Raphson.
- Secante.
La recta secante es una recta que corta una circunferencia en
dos puntos.
Conforme estos puntos de corte se acercan, dicha recta se
aproximan a un punto y, cuando solo existe un punto que toca
la circunferencia, se le llama tangente.
Dados los puntos de intersección A y B puede calcularse la
ecuación de la recta secante empleando para saber la
respuesta de esta operación se emplea en matemáticas la
ecuación de la recta que pasa por dos puntos.
Ene análisis numérico el método de la secante es un método
para encontrar los ceros de una función de forma iterativa.
Uno de los objetivos de este método es eliminar el problema
de la derivada de la función, ya que existe funciones que
describen fenómenos físicos en la vida real, cuya derivada es
muy compleja. El método de la secante es muy similar a la de
Newton con la diferencia principal que en este método de la
secante no requiere de la segunda derivada.
El método se basa en obtener la ecuación de la recta que
pasa por los puntos (xn-1), f(xn-1)) y (xn, f(xn)). A dicha recta
se le llama secante por cortar la grafica de la función.
Posteriormente se escoge como siguiente elemento de la
relación de recurrencias, xn+1, la intersección de la recta
secante con el je de abscisas obteniendo la formula.
Este método, a diferencia del de bisección y regla falsa, casi
nunca falla ya que solo requiere de dos puntos al principio, y
después el mismo método se va retroalimentando. Lo que
hace básicamente es ir tirando rectas secantes a la curva de
la ecuación que se tiene originalmente, y va checando la
intersección de esas rectas con el eje de la x para ver si es la
raíz que se busca.
2.4.- Aplicación de la solución de ecuaciones no lineales.
Un rectángulo tiene 40 cm2 de área y su diagonal es 10 cm.
Halla las dimensiones del rectángulo.
Sean x e y las dimensiones (base y altura) del rectángulo.
Sabemos que el área A de un rectángulo, de base x y altura
y… viene dado por: A= x*y
Por tanto nuestra primera ecuaciones es: 40= xy
Por otro lado, sabemos que la diagonal de un rectángulo es la
hipotenusa de un triangulo rectángulo, como se indica en la
figura:
El teorema de pitágoras nos dice que nos dice que:
D= √ x 2 + y 2 , o lo que es lo mismo, D2=x 2 + y 2
Por tanto como D=10, tenemos la ecuación: 1002=x 2 + y 2
Las dimensiones X e Y del rectángulo serán la solución del
sistema de ecuaciones:
Resolvemos por sustitución despejando la incognita X de la
segunda ecuación:
xy=48 ⇒ x=¿ 48
y ¿
2
48 2304 2
( )
y
+ y 2=100 ⇒
y
+ y =100 ⇒
⇒ 2304 + y 4 =100 y 2 ⇒ y 4 −100 y 2+2304=0
Hacemos el cambio de variable
z= y 2 . z 2−100 z +2304=0
−(−100)± √(−100)2−4.1 .2304 100 ±28 z=64
z=
2 a2.1
=
2 z =36 {
Deshaciendo el cambio: z= y 2 ⇒ y=± √ z
z=64 ⇒ y=± √ 64 ⇒ y=± 8
z=36 ⇒± √ 36 ⇒ y=± 6
Calculamos los valores de la incognita x para cada valor de y.
48 48
y=8⇒ x= = =6
y 8
48 48
y=6⇒ x= = =8
y 6
48 48
y=−6 ⇒ x= = =−8
y −6
Como las dimensiones de un área tiene que ser positivas, la
solución del problema es:
y=8 x =6 o y=6 x =8
2.5- El uso de las herramientas computacionales.
Errores en el cómputo.
Errores en el Computo.
Vocabulario: Palabras, byte y Kilobytes,
N´umeros Binarios = 0; 1= expresados en base 2.
1 byte = 1B ´ 8Bits = Bits para almacenar un car´acter “x” o “y”
1 Palabra = Bits requeridos para almacenar un número ´ 8Bits
= 1B
1 Kilobyte = 1KB = 210 bytes = 1024 bytes
512 Kilobytes = 29 = 524;288Bits =
Enteros y Puntos Flotantes
Los números enteros se representan
xent = sign£binario
¡2N¡1 < xent < 2N¡1
¡2147483647 ' ¡231 < xent < 231 ' 2147483647
para sistemas con palabras de 32 bits, (4bytes)
En general los números reales o puntos flotantes, se
representan en el computadores siguiendo la
Convención
xfloat = sign (xfloat) £M £ 2p
Donde sign (xfloat) es el signo del número xfloat, la mantisa
es M con 0 < M · 1 y el entero positivo,
p; el exponente. Obviamente cualquier número real puede ser
expresado de esa manera. Con algunas
Convenciones adicionales. Signo, exponente y mantisa se
expresan en números binario. Para representaciones
Punto flotantes con palabras de 4 bytes (32 bits) se separan:
un bit para el signo, 8 bits para el
Exponente y los 23 bits restantes para la mantisa.
con [x] la parte entera de x: Dado que p viene representado
por un numero binario de 8 dıgitos el
numero entero máximo que puede representar es 28 = 256 =)
p 2 [0; 255] , para lo cual se establece
una correspondencia a un intervalo que permita exponentes
negativos. De modo que el nuevo exponente
3
P = p ¡ 127 =) P 2 [¡127; 128] : Queda entonces determinar
como expresar la mantisa en números
binarios, para lo cual se sigue la manera tradicional de
expresar números fraccionarios
0;6254 = 6 £ 10¡1 + 2 £ 10¡2 + 5 £ 10¡3 + 4 £ 10¡4
En general una expresión B-fracción
∞
xf =∑ akB−k
n=1
Donde B es la base (10 en el caso decimal, o dos en el caso
binario). Para encontrar la expresión de la mantisa en
números binarios con encontrar el conjunto coeficientes ak en
base B=2.
Vale decir
∞
−k
0.6254=∑ ak 2
n=1
∞
−k
2∗0.6254=1+ 0.2508=a 1+ ∑ ak +12 ⇒a1 =1
n=1
∞
−k
2∗0.2508=1+0.5016=a2 + ∑ ak +22 ⇒ a2=0
n=1
∞
−k
2∗0.5016=1+0.0032=a3 + ∑ ak+ 32 ⇒a3 =1
n=1
∞
−k
2∗0.0032=1+0.0064=a 4 + ∑ ak + 42 ⇒ a4 =0
n=1
∞
−k
2∗0.0064=1+ 0.0128=a 5+ ∑ ak +52 ⇒ a5=0
n=1
∞
−k
2∗0.0128=1+0.0512=a6 + ∑ ak+ 62 ⇒ a 6=0
n=1
∞
−k
2∗0.0512=1+0.1024=a7 + ∑ ak +72 ⇒ a7=0
n=1
∞
−k
2∗0.1024=1+ 0.1024=a8 + ∑ ak + 82 ⇒a8 =0
n=1
Y asi finalmente 0.6254≃10100000
Nótese que es una representación aproximada, limitada por el
numero de cifras permitidas por el tamaño de las palabras y
que tal, y como ocurre en los números decimales el
coeficiente de las potencias mayores va primero.
Underflows y overflows
En precisión simple se usan 4 bytes (32 bits) para representar
un numero, por lo tanto el máximo numero real que puede ser
representado es
2128 =3.4∗1038=⌊ 0 ⌋ ⌊ 11111111 ⌋ ⌊111111111111111111111111 ⌋
Equivalente al minimo será 2−128=2.9∗10−39 . Si se utiliza la doble
precisión, 2 palabras (8 bytes=64 bytes) representaran cada
numero. De los 64 bits, 11 representaran el exponente y los
restantes 52 para la mantisa. En doble precisión
10−322 ≤ doble ≤10 308
Los números menores que los mínimos (tanto en precisión
simple como en doble) son consideraos underflows mientras
que cifras mayores son overflous. El siguiente programa
pseudocodigos identifica estas anomalías
under <------- 1.0
over<------- 1.0
mientras N< Nfinal
haga
under = under/2.0
over = over*2.0
fin mientras
escriba Nfinal, under, over.
Hay compiladores que asignan los underflows y overflows a
cero y eso hay que chequearlo.
