ERRORES POR TRUNCAMIENTO
Este tipo de error es debido a la limitación que existe en los cálculos, es decir
al número de operaciones y que en algunos casos existen calculo infinito como por
ejemplo la función:
Sen(X)=X-X3/3! +X5/5! -X7/7!..........
La función indicada anteriormente solo puede calcularse en forma aproximada, el
número de elementos que componen esta función es infinito, pero el computador
solo realiza operaciones finitas, entonces tomaremos 10,15 o 20 términos de la
función para calcular su valor, en este caso habremos truncado la función.
Note que el error de truncamiento, a diferencia del error de redondeo, no
depende directamente del sistema numérico que se emplee. El Teorema de
Taylor nos ayuda a encontrar el error por truncamiento, en la evaluación de
funciones, que no se encuentran en una librería estándar de una calculadora o
computadora.
TEOREMA DE TAYLOR
El polinomio de Taylor de grado “n” para la función f alrededor de xo.
Es
El residuo o error de truncamiento asociado con Pn(x).
Entonces el valor de la función exacta será:
f(x) = Pn(x) + Rn(x)
Consideramos Pn(x) como valor aproximado y Rn(x) error absoluto o de
truncamiento
En el caso específico de que xo = 0 el polinomio de Taylor se conoce como el
polinomio deMaclaurin y la serie de Taylor se conoce como la serie de
Maclaurin.
Ejemplos. -
1.- Sea la función f(x) = ⅇ−x
Encuentre un Polinomio de Taylor de Tercer Grado para f alrededor de X0=1,
aproxime a ⅇ−0.99 usando un polinomio de Taylor ¿Cuántas cifras decimales
de precisión se esperan?
Solución. -
X0=1 X=0.99 n=3
Si n=3 debemos derivar 4 veces para incluir el residuo
f(x0) =e-1 f’(x0) = -e-1 f “(x0) = e-1 f “’(x0) = -e-1
§(x) valor comprendido entre x y x
0
§(x)
Se asumirá = x0 = 1
Por lo tanto f IV(x0) =e-1
1 (−0.01)2 (−0.01)3
Luego P3(x)= (1 − (0.01) + − )=0.3715766907
𝑒 2¡ 3¡
1 (−0.01)4
R3(x)= =0.1532831005 x 10-9
𝑒 4¡
Lo que nos indica que el resultado es preciso hasta el 9 digito decimal.
2.- Use un polinomio de Taylor para la función f(x)=Ln(x) alrededor de “e”
para encontrar una aproximación de ln (3) y que sea exacta a 1x10-4
Solución. -
No sabemos cuál será el grado “n” para alcanzar la precisión solicitada, por lo
tanto, debemos tantear, asumiremos con n=2
X0=e X=3 n=2
(𝑥 − 𝑥0 )2
𝑃2 (𝑥) = 𝑓(𝑥0 ) + 𝑓 ′ (𝑥0 )(𝑥 − 𝑥0 ) + 𝑓′′(𝑥0 )
2¡
§ x ) (x−𝑥
𝑅2 (𝑥) = 𝑓 𝐼𝐼𝐼 ( ( )
3¡
0)
1 1
f(x0)=Ln(e)=1 f ’(x0)= f “(x0)=-
𝑒 𝑒2
1 1 (3−𝑒)2
P2(x)=1+ (3 − 𝑒) − =1.098267872
𝑒 𝑒2 2¡
§(x) 2
Se asumirá = x0 = e fIII(§(x))=
𝑒3
3 2 (3−𝑒)3
𝐼𝐼𝐼 §(x) (𝑥−𝑥0 )
R2(x)= 𝑓 ( ) = =0.3710563623x10-3
3¡ 𝑒3 6
Debe ser menor a 1 x 10-4
No cumple ¡¡¡¡¡¡¡¡¡¡¡¡ Calcular con n=3
2
f ”’(x0)=
𝑒3
1 1 (3 − 𝑒)2 2 (3 − 𝑒)3
𝑃3 (𝑥 ) = 1 + (3 − 𝑒) − 2 + 3
𝑒 𝑒 2¡ 𝑒 3¡
P3(x)=1.098638929
6
fIV(§(x))=− si §(x)=X0=e
𝑒4
6 (3−𝑒)4
R3(x)= −
𝑒4 24
=-0.4806957414 x10-4< 1x 10-4 OK
ALGORITMOS ESPECIALES
Antes de analizar Ecuaciones no lineales veremos, dos algoritmos que por sus
características serán de importancia, para entender mejor el tema que estudiaremos
(Ecuaciones no lineales).
EVALUACION DE FUNCIONES
Como se vio anteriormente el tratar de evitar excesivos cálculos nos beneficia
en la precisión de los resultados. Un ejemplo típico de este análisis es el polinomio
de Horner cuya estructura es:
Pn(x) =a0xn+a1xn-1+a2xn-2.....an-1x+an
a0,a1,a2….an=coefecicientes constantes.
x=variable independiente. y “n” grado del Polinomio
También este polinomio podemos derivarlo y tendrá la sgte estructura.
Qn-1(x) =P’n(x) = na0xn-1+ (n-1) a1xn-2+a2xn-2.....an-1
Podemos utilizar otro criterio para evaluar este polinomio o siguiendo el
ejemplo anterior tenemos lo sgte para un polinomio de 5to Grado:
P5(X)= a0x5 +a1x4+a2x3 +a3x2+a4x +a5
El número de o. aritméticas son: 15 Multiplicaciones
5 Sumas
20 Total operaciones Aritméticas
Volcando el Polinomio y factorizando sucesivamente tendremos:
P5(X)=a5+x(a4+a5x+a2x2+a1x3+a0x4)
P5(X)=a5+x(a4+(a3+a2x+a1x2+a0x3))
P5(X)=a5+x (a4+x (a3+x (a2+xa1 +a0x2)))
P5(X)=a5+x (a4+x (a3+x (a2+x (a1 +a0x))))
El número de o. aritméticas son: 5 Multiplicaciones
5 Sumas
10 Total operaciones Aritméticas
Como se sabe, menos operaciones aritméticas, mayor precisión
Por lo tanto, el algoritmo de Horner es muy eficiente con solo factorizar
adecuadamente el polinomio
Del ejemplo anterior podemos decir que el polinomio de Horner se reescribe:
Pn(X) =an+x (an-1+x (an-1+x (an-2…. +x (a1+a0x)….)))
P0=a0 p1=a1+p0x p2=a2+p1x p3=a2+p2x
Pi=ai+pix p = ai +p x
I=1 n
Po=0 Qo=Po
P1=P0 Q1=P1+Xqo
P2=P1+XP1 Q2=P1+XQ1
I=1 (N-1) Q2=Pi+Xx-1
q=p+xq
Diagrama de Flujo Evaluación Polinomio de Horner
Donde:
N = Grado de Polinomio
I = Contador
A(I) = Coeficiente del Polinomio
P = Polinomio Horner
Q = Derivada Polinomio Horner
X = Variable Independiente
Antro
Intro
Intro
Mostrar
RELACIONES DE RECURRENCIA
Existen ecuaciones lineales con dos, tres o más variables cuya única forma de
solución es asumir dos o más variables libre, tomaremos como ejemplo una ecuación
lineal con tres incógnitas.
axi+bxi-1+cxi-2=0
a,b,c=cte.
xi=2 n
ax2+bx1+cx0=0
x0
Conocidos se obtiene x2
xl
axi+bxi-1+cxi-2=0
xi=-(bxi-1+cxi-2)/a
I=2,3,4………
Podemos hacer un diagrama de flujo a partir de dos raíces iniciales de tal
manera que se puedan generar un conjunto de raíces solución.
Ejemplo
-3x2+4x1-x0=0
si x0=1
x1=2
x2=(-4x1+x0)3=-7/3
-3x3+4x2-x1=0
x3=(-4x2+x1) /3=34/9
3x4+4x3-x2=0
Podemos continuar de forma recurrente hasta
“n” soluciones-
Diagrama de Flujo Relaciones de Recurrencia de 3 Términos
Donde:
A, B, C = Constantes
N = Número de Soluciones >2
I = Contador
X(I) = Soluciones
X (0), X (1) = Soluciones Iniciales
Alternativa 1
Intro
Intro
Intro
Alternativa 2
Mostrar
