Escuela Superior Politécnica de Chimborazo Facultad de Mecánica Carrera de Mantenimiento Industrial
Escuela Superior Politécnica de Chimborazo Facultad de Mecánica Carrera de Mantenimiento Industrial
CHIMBORAZO
FACULTAD DE MECÁNICA
MATERIA:
MÉTODOS NUMÉRICOS
INTEGRANTES:
JOEL CHOCO 2710
DAVID GAIBOR 2688
ANELLY RUIZ 2671
AÑO LECTIVO:
2022-2023
CONTENIDO
ÍNDICE DE TABLAS.....................................................................................................................2
ÍNDICE DE GRÁFICOS................................................................................................................3
CAPÍTULO 1...................................................................................................................................4
1. ¿Qué es la Interpolación?....................................................................................................4
EJERCICIOS RESUELTOS........................................................................................................14
EJERICICIOS PROPUESTOS....................................................................................................40
BIBLIOGRAFÍA...........................................................................................................................43
ÍNDICE DE TABLAS
ÍNDICE DE GRÁFICOS
CAPÍTULO 1
1. ¿Qué es la Interpolación?
La interpolación consiste en que se debe encontrar con que se tiene que estimar valores
intermedios entre puntos asociados con datos. El método más común que se usa para este
propósito es la interpolación polinomial. Recuerde que la fórmula general para un
polinomio de enésimo grado es
2 n
f ( x )=a0 +a1 x+ a2 x +… … … … … ..+ an x
1.1 Tipos de interpolación
Interpolación por Lagrange
Interpolación lineal
Interpolación polinomial
Interpolación de newton
Interpolación con splines cúbicos
1.1.1. Polinomios de Lagrange
El método de aproximación polinomial estudiado en la sección anterior requiere la
solución de un sistema de ecuaciones algebraicas lineales que, cuando el grado del
polinomio es alto, puede presentar inconvenientes. Existen otros métodos de aproximación
polinomial en que no se requiere resolver un sistema de ecuaciones lineales y los cálculos
se realizan directamente; entre éstos se encuentra el de aproximación polinomial de
Lagrange.
En éste nuevamente se parte de una función desconocida f (x)dada en forma tabular, y se
asume que un polinomio de primer grado (ecuación de una línea recta) puede escribirse:
p ( x ) =a0 ( x – x 1) + a1 ( x – x 0 ) (5.12)
donde x 1 y x 0 son los argumentos de los puntos conocidos [ x0 , f (x0 )],[ x 1 , f ( x1 )], y a 0 y a 1
son dos coeficientes por determinar. Para encontrar el valor de a 0, se hace x=x 0 en la
ecuación 5.12, que al despejar da
p (x 0) f (x 0)
a 0= = (5.13)
x0 −x1 x 0−x 1
Y para hallar el valor de a 1, se sustituye el valor de x con el de x 1, con lo que resulta
p (x1 ) f (x 1)
a 1= = (5.14)
x 1−x 0 x1 −x0
de tal modo que, al sustituir las ecuaciones 5.13 y 5.14 en la 5.12, queda
f (x 0 ) f (x1 )
p ( x) =
x 0−x 1
( x – x1)+
x 1−x 0
( x – x 0 ) (5.15)
Donde
( x – x 1 ) ( x−x 2)
L0 ( x ) =
(x ¿ ¿ 0−x 1)( x 0− x2 )¿
( x – x 0 ) (x−x 2 )
L1 ( x ) =
( x ¿ ¿ 1−x 0 )(x 1−x 2) ¿
( x – x 0 ) (x−x 1)
L2 ( x ) =
( x ¿ ¿ 2−x0 )( x 2−x 1)(5.21)¿
Pn ( x ) =L0 ( x ) f ( x 0 ) + L1 ( x ) f ( x 1 ) +… … … .+ Ln ( x ) f ( x n )
( x – x 1 ) ( x−x 2)
L0 ( x ) =
(x ¿ ¿ 0−x 1)( x 0− x2 )¿
( x – x 0 ) (x−x 2 )
L1 ( x ) =
( x ¿ ¿ 1−x 0 )(x 1−x 2) ¿
.
.
.
n
( x – x j)
Li ( x )= ⏞
⏟
ℿ
( x ¿ ¿ i−x j )(5.23)¿
⏟
j=0
j≠ i
Ejemplo
Obtener el polinomio de interpolación usando la forma de interpolación de
LaGrange con la siguiente tabla de valores, e interpolar en el punto x=4
XK −6 6 −1
YK 8 −16 2
Table 1 Datos para interpolación
n n n
x−x 1
p( x )=∑ Y k ∏ =∑ 4 k L k ( x ) ¿ ¿
k=0 i=0 ,i ≠ k x k −x 1 k=0
(x ¿ ¿ 0 ; y 0)=(−6,8)¿
¿ ¿)¿(6,16)
(x ¿ ¿ 2; y 2 )¿ =(−1 ,−2)
( x −x1 ) ( x−x 2 ) 1 1 1 1
L0 ( x ) = = ( x−6 )( x +1 )= x 2− x−
( 0 1) ( 0 2)
x −x x −x 60 60 12 10
( x−x 0 ) ( x−x 2 ) 1 1 1 1
L1 ( x ) = = ( x +6 ) ( x+1 )= x 2 + x+
( x 1−x 0 ) ( x 1−x 2) 84 84 12 14
( x−x 0 ) ( x−x 1 ) 1 −1 2 36
L2 ( x )= = ( x+6 )( x−6 )= x +
( x 2−x 0 ) ( x 2−x 1 ) 35 35 35
2
p( x )=∑ y k Lk ( x ) =8 L0 ( x )−16 L1 ( x )−2 L2 ( x ) =−4−2\
k=0
Gráfico 1 Interpolación por el método de Lagrange
1
L0 (−4 )=
2
−1
L1 (−4 )=
14
4
L2 (−4 )=
7
2
p(−4)=∑ y k Lk (−4 )=8 L0 (−4 )−16 L1 (−4 )−2 L2 (−4 )=4
k =0
p( x )=−4−2 x
−4 −2 −4
−2 4
Table 2 Datos para interpolación por Lagrange
P(−4)=−2(−4)−4=4
1.1.2. Interpolación polinomial de newton en diferencias divididas
Forma general de los polinomios de interpolación de Newton
El análisis anterior puede generalizarse para ajustar un polinomio de enésimo grado a n + 1
puntos asociados con datos. El polinomio de enésimo grado es
f n ( x )=b0 +b1 ( x−x 0 ) +…
+b n ( x −x0 ) ( x−x 1 ) … .(x−x n−1)
Como se hizo antes con las interpolaciones lineales y cuadráticas, los puntos asociados con
datos se utilizan para evaluar los coeficientes b 0 , b 1 ,... , bn. Para un polinomio de n-ésimo
grado se requieren n + 1 puntos asociados con dato s :
[ x 0 , f ( x 0)],[ x 1 , f ( x 1)], ... ,[xn , f (xn)]. Usamos estos puntos asociados con datos y las
siguientes ecuaciones para evaluar los coeficientes:
b 0=f ( x 0)
b 1=f [x 1 , x 0 ]
b 2=f [x 2 , x 1 , x 0 ]
.
.
.
b n=f [x n , xn−1 , … . , x 1 , x 0 ]
donde las evaluaciones de la función colocadas entre paréntesis son diferencias divididas
finitas.
Como se dijo antes, existe una gran variedad de formas alternativas para expresar una
interpolación polinomial. El polinomio de interpolación de Newton en diferencias
divididas es una de las formas más populares y útiles. Antes de presentar la ecuación
general, estudiaremos las versiones de primero y segundos grados por su sencilla
interpretación visual
Partiendo de n puntos (x, y), podemos obtener un polinomio de grado n − 1. El método que
se utilizaría es el de las diferencias divididas para obtener los coeficientes, el cual facilita la
tarea de resolver un sistema de ecuaciones usando el cociente de sumas y restas. Dada una
colección de n puntos de x y sus imágenes f(x), se pueden calcular los coeficientes del
polinomio interpolante utilizando las siguientes expresiones:
f |x k|=f (x k )
f ( x k +1 )−f ( x k )
f |x k , x k + 1|=
x k + 1−x k
f ( x k +1 , x k + 2 )−f (x k , x k +1)
f |x k , x k + 1, x k +2|=
x k +2−x k
Finalmente, a partir de los valores obtenidos, se pueden obtener dos formas de representar
el polinomio:
Progresivo (desde 0 hasta n − 1):
Pn−1=f ( x 0 )+ f ( x 0 , x 1 )∗( x−x 0 ) +f ( x 0 , x 1 , x2 )∗( x−x 0 ) ( x−x 1 ) ++++… … … … ..+ f ( x 0 , x1 , … … , xn )∗( x−x 0 )∗(
Sean dos puntos( x 0 , y 0 ) ,(x1 , y 1) la interpolación lineal consiste en hallar una estimación
del valor y, para un valor x tal que x0<x<x1.
f ( x1 ) + f ( x 0 )
f 1 ( x )=f ( x 0) + ( x −x 0)
x1 −x0
En cada uno de estos intervalos, debemos definir una función polinomial de grado 2,
como sigue:
{ }
a1 x 2 +b1 x+ c 1
s ( x )= a x 2 +b x+ c si x ∈ [ 3,4,5 ] , [ 4,5,7 ] , [ 7,9 ]
2 2 2
2
a2 x +b 2 x+ c 2
Primero, hacemos que el splines pase por los puntos de la tabla de datos. Es decir, se
debe cumplir que:
s ( 9 )=0.5→ 81 a3 +9 b 3+ c 3=0.5
Hasta aquí, tenemos un total de 6 ecuaciones vs. 9 incógnitas.
El siguiente paso es manejar la existencia de las derivadas continuas. En el caso de las
splines de grado 2, necesitamos que la spline tenga derivada continua de orden k −1=1, es
decir, primera derivada continua.
Calculamos primero la primera derivada:
{ }
2 a1 x+ b1
'
s ( x )= 2 a2 x+ b2 si x ∈ [ 3,4,5 ] , [ 4,5,7 ] , [ 7,9 ]
2 a3 x+ b3
Vemos que esta derivada está formada por segmentos de rectas, que pudieran presentar
discontinuidad en los cambios de intervalo. Es decir, las posibles discontinuidades son
x=4.5 y x=7 . Por lo tanto para que sx sea continua, se debe cumplir que:
o lo que es lo mismo,
9 a 1+ b1=9 a 2+ b2
También debe cumplirse que:
2 a2 ( 7 ) +b2 =2 a3 x +b3
o lo que es lo mismo,
14 a 2+ b2=14 a3 x +b 3
Así, tenemos un total de 8 ecuaciones vs. 9 incógnitas; esto nos da un grado de libertad
para elegir alguna de las incógnitas. Elegimos por simple conveniencia a 1=0.
De esta forma, tenemos un total de 8 ecuaciones vs. 8 incógnitas. Estas son las
siguientes:
3 b1 +c 1
4.5 b 1+ c 1=1
49 a2+ 7 b2 +c 2=2.5
49 a3+ 7 b3 +c 3=2.5
81 a3 +9 b 3+ c 3=0.5
b 1=9 a2 +b 2
14 a 2+ b2=14 a3 +b3
Sustituyendo estos valores (junto con a 1=0 ), obtenemos la función spline cuadrática que
interpola la tabla de datos dada:
{ }
−x +5.5
s ( x )= 0.64 x −6.76 x+18.46 si x ∈ [ 3,4,5 ] , [ 4,5,7 ] , [ 7,9 ]
2
La gráfica que se muestra a continuación contiene tanto los puntos iniciales de la tabla de
datos, así como la spline cuadrática. Esta gráfica se generó usando Matemática.
(CHAPRA, 2015)
1.1.5. Interpolación cubica a trozos: splines cúbicos
La idea básica del spline cúbica es representar la función mediante una función cúbica
diferente en cada intervalo entre los puntos de datos. Es decir, si hay N puntos de datos,
entonces el spline f(x) es la función.
En este caso, cada polinomio P(x) a través del cual construimos los splines en [m, n] es de
grado 3. Esto significa que tendrá la forma P ( x ) =a x 3+ b x2 +cx + d . Se tendrá entonces
cuatro variables para cada intervalo (a, b, c, d) y una nueva condición para cada punto
común en dos intervalos, con respecto a la segunda derivada.
Al igual que en los splines cuadráticos, existen condiciones, tales como asegurar que la
derivada en un punto siempre coincide para ambos "lados" de la función definida en
tramos que pasa a través de un punto tan común. Asegurar que la segunda derivada en un
punto siempre coincide para ambos "lados" de la función definida en los tramos que pasan
por ese punto común.
Consideremos n + 1 puntos ( x 0 , y 0 ) ,( x 1 , y 1 ),....(x n ,Y n) cuyas abscisas están ordenadas de
manera Creciente a=x 0 < x1 < x 2 <..< x n−1< x n =b (nodos).
La función S: [a, b]→→ R es un spline cúbico si existen n polinomios S K (x ), cada uno
definido sobre un intervalo, los cuales lo podemos escribir en la forma:
S K ( x ) =a K ¿
Para x ∈[ x k , x k +1]y k= 0, 1, 2, …, n-1
. S ( x k )= y k donde k= 0, 1, 2, …, n condición de interpolación
. Sk ( x k+1 ) =S k+1 (x k +1) donde k =0 , 1 ,2, … , n−3 , n−2 continuidad de los splines
. S ' k ( x k+1 ) =S ' k+1 ( x k +1) donde k= 0, 1, 2, … , n−3 , n−2continuidad de la derivada
. S ' ' k ( x k+1 ) =S ' ' k +1(x k +1) donde k= 0, 1, 2, … , n−3 , n−2 continuidad de la segunda
derivada
Una de las siguientes condiciones de frontera se satisface
'' ''
S ( x 0 ) =S ( x n )=0 frontera libre o natural
' ' '
S ( x0 ) =f ( x 0 ) y S ( x n )=f ' ( x n) frontera sujeta
3 2
s(x )=a1 x + b1 x +c 1 x +d 1 si x ∈[ x 0 , x 1]
3 2
s(x )=a n−1 x + bn−1 x +c n−1 x +d n−1 si x ∈[ x n−1 , xn ]
(HURTADO, 2015)
EJERCICIOS RESUELTOS
1. Determinar la función de interpolación que pasa por los puntos (−1,0 ) , ( 4,2 )
tenemos los putos
f ( x 0 , y 0 )=(−1,2 )
Q ( x 1, y1 ) =( 4,2 )
y 1− y 0
f ( g )= y 0 +
x 1 , −x 0
( x −( x 0 ) )
2−0
f ( g )=0+ ( x −(−1 ) )
4−(−1 )
2 2
¿ x+
3 5
2 2 4
f ( X )= ( 1 ) + =
5 5 5
Interpolacion lagrange
3. Aplicar la interpolacion de lagrange para encontrar el polinomio que interpole los
datos siguientes:
( 0,1 ) , (1,3 ) y ( 2,0 )
F (xi) L(xi)
n 3 f ( x 0 ) =1 L¿
Xi 012 f ( x 1 ) =3 2
L ( x1 ) =−(x −2 x )
f ( x 2 ) =0 L ( x2 ) =1/2(x 2−x)
Iteracion 2
( x−x 0 ) ( x−x 2 )
L1 ( x )=
( x 1−x 0 ) ( x 1−x 2)
( x−0 )( x−2 )
L1 ( x ) =
( 1−0 )( 1−2 )
¿−( x −2 x )
2
Iteracion 3
( x−x 0 ) ( x−x 1 )
L2 ( x )=
( x 2−x 0 ) ( x 2−x 1 )
( x−0 )( x−1 )
L2 ( x ) =
( 2−0 )( 2−1 )
¿ 1/2(x 2−x)
P ( x ) =f ( x 0 ) L0 ( x ) + f ( x 1 ) L1 ( x ) +f ( x2 )L2 ( x )
x p(x )
−2 −18
−1 −6
0 1
1 3
2 0
3 −8
Table 6 Datos obtenidos del reemplazo en la ecuación de interpolación
4. SPLINE CUBICO
X −2 −1 1 3
Y 3 1 2 −1
Table 7 Datos para interpolar por spline cúbico
Interpolar el valor correspondiente a x=2
{
a1 x 3 +b 1 x 2+ c 1 x+ d 1 para−2≤ x ←1
S ( x )= a x 3 +b x 2+ c x+ d para−1≤ x <1
2 2 2 2
3 2
a3 x +b3 x + c 3 x +d 3 para 1 ≤ x ≤3
s1 (−2 )=3
s1 (−1 )=1
s2 (−1 )=1
s2 ( 1 )=2
s3 ( 1 )=2
s3 ( 3 )=−1
a 1(−2)3 +b 1(−2)2 +c 1 (−2)+ d 1=3
3 2
a1 (−1) +b 1 (−1) +c 1 (−1)+d 1=1
a2 (−1)3 +b 2 (−1)2 +c 2 (−1)+d 2=1
3 2
a2 (1) +b 2 (1) +c 2 (1)+ d 2=2
a3 (1)3 +b 3 (1)2 +c 3 (1)+ d 3=2
3 2
a3 (3) +b3 (3) +c 3 (3)+d 3=−1
−8 a1 + 4 b1−2c 1 +d 1=3
−a 1+ b1−c 1 +d 1=1
−a 2+ b2−c 2 +d 2=1
a 2+ b2 +c 2+ d 2=2
a3 + b3 +c 3+ d 3=2
27 a 3+ 9 b3 +3 c 3+ d 3=−1
Se observa
' 2
s k ( x ) =3 ak x +2 bk x +c k
Por lo tanto
s' 1 (−1 ) =s ' 2 (−1 ) y s' 2 (1 ) ¿ s ' 3 ( 1 )
Corresponden a:
2 2
3 (−1 ) a1+ 2 (−1 ) b1 +c 1=3 (−1 ) a2 +2 (−1 ) b 2+ c 2
2 2
3 ( 1 ) a2+ 2 ( 1 ) b2 +c 2=3 ( 1 ) a3 +2 (1 ) b3 + c3
3 a1 −2b 1+ c 1−3 a2 +2 b2−c2=0
3 a2 −2b 2+ c 2−3 a3 +2 b3−c 3=0
Se observa
'
s k ( x ) =6 ak x+ 2b k
Por lo tanto
' ' ' '
s 1 (−1 ) =s 2 (−1 ) +1 y s 2 ( 1 ) ¿ s 3 ( 1 )
2
6 (−1 ) a1 +2 b1=6 (−1)a 2+2 b2
6( 1)a2 +2 b2=6(1)a 3+ 2b 3
−6 a 1+2 b 1+6 a2−2 b2=0
6 a 2+2 b 2−6 a3 +2 b3 =0
Se observa
s' k ( x ) =6 ak x+ 2b k
Por lo tanto
' '
s ' 1 (−2 )=0 y s ' 3 ( 3 )
6 (−2 )2 a1 +2 b1=0
6 ( 3 ) a2 +2 b2=0
−12 a1 +2 b1=0
18 a3 +2 b3 =0
Sistema de ecuaciones
−8 a1 + 4 b1−2c 1 +d 1=3
−a 1+ b1−c 1 +d 1=1
−a 2+ b2−c 2 +d 2=1
a 2+ b2 +c 2+ d 2=2
a3 + b3 +c 3+ d 3=2
27 a 3+ 9 b3 +3 c 3+ d 3=−1
3 a1 −2b 1+ c 1−3 a2 +2 b2−c2=0
3 a2 −2b 2+ c 2−3 a3 +2 b3−c 3=0
−6 a 1+2 b 1+6 a2−2 b2=0
6 a 2+2 b 2−6 a3 +2 b3 =0
−12 a1 +2 b1=0
18 a3 +2 b3 =0
Xk 1 −4 −7
Yk 10 10 34
Table 10 Datos para interpolación por Lagrange
2
6 X +3 X
n n n
x−x 1
p( x )=∑ 4 k ∏ =∑ 4 kLk ( x )
k=0 i=0 ,i ≠ k x k −x 1 k=0
( x ¿ ¿ 0 , y 0)=( 1,10 ) , ( x 1 , y 1 ) =(−4,10 ) , ( x 2 , y 2 )=(7,34) ¿
( x−x1 ) ( x−x 2 ) 1 1 2 1 7
L0 (x)= = ( x+ 4 ) ( x+ 7 )= x + +
( x0 −x1 ) ( x 0−x 2 ) 40 40 40 10
( x−x 0 ) ( x−x 2 ) −1 −1 2 2 7
L1 ( x )= = ( x −1 )( x +7 )= x− +
( x 1−x 0 ) ( x 1−x 2) 15 15 5 15
( x−x 0 ) ( x−x 1 ) 1 1 2 1 1
L2 ( x )= = ( x−1 ) ( x+ 4 )= x + x−
( x 2−x 0 ) ( x 2−x 1 ) 24 24 8 6
1
L0 (−3)=
10
16
L1 (−3)=
15
−1
L2 (−3 )=
16
2
p(−3)=∑ y k L k (−3 )=10 L0 (−3 )+10 L1 (−3 ) +34 L2 (−3 )=6
k=0
−3 1 3 6
1 0 6
Table 11 Datos obtenidos para interpolar en la ecuación correspondiente
p(−3)=(1(−3)+3)∙(−3)+6=6
7. Calcular el polinomio de Lagrange usando los siguientes datos
X 1 −3 5 7
Y −2 1 2 −3
Table 12 Datos para interpolar por Lagrange
n n
(x −xi )
p(x )=∑ f ( x 1) L1 (x )Li ( x )= ∏ ¿
i=0 i=0 j ≠1 ( x ¿ ¿ i−x j )
−( x−3 ) (x−5)(x−7) ( x−1 )( x−5 )( x−7 ) ( x−1 ) ( x+3 ) ( x−7) ( x +1 ) ( x+3 ) (x−5)
f ( x)= - − −
48 320 40 40
f (x)=−0.0739 x 3+0.3906 x 2+ 0.624 x−2.978
(x ¿ ¿ 0 ; y 0)=(−6,8)¿
(x ¿ ¿ 1; y 1 )=(6,16)¿
(x ¿ ¿ 2; y 2 )=(−1 ,−2)¿
( x −x1 ) ( x−x 2 ) 1 1 2 1 1
L0 ( x ) = = ( x−6 )( x +1 )= x − x−
( x 0 −x1 ) ( x 0−x 2 ) 60 60 12 10
( x−x 0 ) ( x−x 2 ) 1 1 2 1 1
L1 ( x ) = = ( x +6 ) ( x+1 )= x + x+
( x 1−x 0 ) ( x 1−x 2) 84 84 12 14
( x−x 0 ) ( x−x 1 ) 1 −1 2 36
L2 ( x )= = ( x+6 )( x−6 )= x +
( x 2−x 0 ) ( x 2−x 1 ) 35 35 35
2
p( x )=∑ y k Lk ( x ) =8 L0 ( x )−16 L1 ( x )−2 L2 ( x ) =−4−2\
k=0
p( x )=−4−2 x
−4 −2 −4
−2 4
Table 14 Datos obtenidos para interpolar en la ecuación correspondiente
P(-4)=-2(-4)-4=4
9. Obtener el polinomio de interpolación usando la fórmula de interpolación de
newton en diferencia divididas con los datos de la tabla que aparece a
continuación, interpolar en x=3
XK −4 6 2
YK −20 −30 −2
Table 15 Datos para interpolar por diferencias divididas
n i=1
pn ( x ) =f ( x 0 ) + ∑ f [ x 0 , x 1 , … x i ] ∏ ( x−xj )
i=1 j=0
−4 −20
6 −30 −1
2 −2 −7 −1
Table 16 Datos obtenidos para reemplazar en la ecuación correspondiente
P(0)(x )=−20
P 1(x)=−1( x+ 4)+ p 0( x)=−x−24
2
P(2)(x)=−( x +4 )( x−6)+ p1( x)=x−x
P( x )=−20−1( x+ 4)−1(x + 4)( x +6)=x−x 2
2
P( x )=x−x
( x− x0 ) ( x−x 2 ) ( x−x 3 ) −1 −1 3 7 2 7 3
L1 ( x ) = = ( x−5 ) ( x−3 ) ( x+ 1 )= x+ x− x−
( x 1−x 0 ) ( x 1−x 2) ( x 1−x 3 ) 320 320 320 320 64
2
p( x )=∑ y k Lk ( x ) =−387 L0 ( x ) +333 L1 ( x )−83 L2 ( x )−363 ( x )=2+3 x−x 2−3 x 3
k=0
Gráfico 9 Interpolación por Lagrange
11. Interpolación lineal
Km 50 60 70 ?
$ 3 3.5 ? 90
{3=m ⋅50+n−1
3,5=m⋅60+ n
3.5=60(0.05)+n
3.5=3+n
−3=−50 m−n
3.5=60 m+ n n=3.5−3
0,5=10 m
. n=0.5
0.5
m=
10
m=0.05
y=0.05 x +0.5
y=0.05 ( 70 )+ 0.5
y=4
5=0.05 x +0.5
5−0.5=0.05 x
5−0.5
x=
0.05
x=90
12. Interpolación de Lagrange
Aplicar la interpolación de Lagrange para encontrar el polinomio que interpole los
siguientes datos(0,1)(1,3)(2,0)
Datos
n−1
n=3 Pn ( x ) =∑ f ( x ) Li ( x )
i=0
n−1
x −xj
i=0,1,2 Li ( x )=∏
i=0 xi−xj
j=0,1,2
Interacción 1
i=0 j=1,2
(x−x 1 )( x−x 2)
Lo ( x )=
( x o−x 1 )( x o−x 2)
( x−1)( x−2)
¿
( 0−1)(0−2)
2
x −2 x−x +2
¿
(−1 ) (−2 )
2
x −3 x +2
¿
2
1 2
¿ x −3 x+2
2
Interacción 2
i=0 j=0,2
( x− xo )(x−x 2 )
L1 ( x)=
( x1 −x0 )(x 1−x 2)
(x−0)( x−2)
¿
(1−0)(1−2)
2
x −2 x−0+ 2
¿
1 (−1 )
x 2−2 x
¿
−1
Interacción 3
¿−(x 2−2 x)
i=2 j=0,1
( x−x o)(x−x 1 )
L2 ( x)=
( x2 −x o)(x 2−x 1)
(x−0)( x−1)
¿
(2−0)(2−1)
2
x −x
¿
2( 1)
1
¿ (x ¿¿ 2−x )¿
2
1 2
¿ ( x −3 x+ 2 )−3 ( x −2 x )
2
2
1 2 3 2
¿ x − x+1−3 x +6 x
2 2
5 2 9
¿− x + x +1
2 2
P( x ) −18 −6 1 3 0
X −2 −1 0 1 2
f (x1 )−f (x o )
1 ra=f (x o , x 1)=
x 1−x 0
0−1
¿
−1−2
¿−1
2da
f (x 1 , x 2)−f ( xo , x 1)
f ( x o , x 1 , x 2)=
x 2−x 0
1−(−1)
¿
0−(−1)
=1
x Y
0−1
=1
-2 1 −1−2
1−(−1)
=1
-1 0 0−(−2)
1−0
=1
0−(−1)
0 1
x 0=−2
x 1=−1
x 2=0
P ( x ) =ao + a1 ( x −x o ) +a2 ( x−x o ) ( x−x 1 )
¿ 1±1(x−(−2))+1( x−2)(x−1)¿ 1−x−2+ x 2 +3 x+ 2
2
¿ x + 2 x +1
0 -1
1 6
2 31
3 18
( x−1)( x−2)(x−3)
¿
(0−1)(0−2)(0−3)
x 3 −6 x2 +11 x−6
¿
−6 ( x−x )( x −x )(x−x )
0 2 3
L1 ( x)=
( x1 −x0 )( x 1−x 2)(x 1−x 3 )
( x−0)( x−2)(x−3)
¿
(1−0)(1−2)(1−3)
3 2
x −5 x +6 x
¿
2
( x−0)( x−1)(x−3)
¿
(1−0)(2−0)(2−3)
x 3 −4 x 2 +3 x
¿
−2
( x−x 0 )( x−x 1)( x−x 2 )
L3 ( x)=
( x 3−x 0)(x 3−x 1 )( x3 −x2 )
( x−0)(x−1)( x−2)
¿
(3−0)(3−1)(3−2)
3 2
x −3 x +2 x
¿
6
1 3 6 3 31 3 18 3
P3 ( x ) = ( x −6 x +11 x−6 ) + ( x −5 x +6 x )− ( x −4 x +3 x ) + ( x −3 x + 2 x )
2 2 2 2
6 2 2 6
x3 2 11 5 2 31 3 2 93 18 3 2
¿ −x + x −1+ 3 x −15 x +18 x − x +62 x + x + x −9 x +6 x
6 6 2 2 6
28 3 2 62
¿− x +37 x − x−1
3 3
-2 -4
0 0
2 4
3 21
a¿
a¿
a¿
a¿
−8 a+ 4 b−2 c +d=4
d=0
8 a+ 4 b+ 2 c+ d=4
27 a+ 9 b+3 c +d=21
−8 a+ 4 b−2 c=−4
8 a+4 b+2 c=4
8 b=0
.
b=0
8 a+ 2c =4
27 a+ 3 c=21
27 a+ 3 ( 2−4 a )=21
15 a=15
a=1
4−8 a
c=
2
¿ 2−4 a
c=−2
Gráfico 12 Interpolación polinómica
16. Empleé la porción de la tabla de vapor que se da para el H2O supe calentada a 200
MPa, para a) encontrar la entropía correspondiente s para un volumen específico v
de 0.108 m3/kg con interpolación lineal, b) encontrar la misma entropía
correspondiente con el uso de interpolación cuadrática, y c) hallar el
Volumen correspondiente a una entropía de 6.6 con el empleo de interpolación
inversa.
m3 0.10377 0.11144 0.1252
V( )
kg
kJ 6.4147 6.5453 6.7664
s( )
Kg∗K
Table 25 Datos para hallar el volumen y la entropía
Para la interpolación lineal tomamos los 2 primeros valores de la tabla
x 0.10377 0.111444
y 6.4147 6.5453
Table 26 Datos para interpolación lineal
aplicamos la siguiente ecuación
P1 ( x )=a0 ( x−x 1 )( x 0 −x1 ) + a1 ( x −x0 ) ( x 1−x 0 )
P1 ( x )=6.4147 ( x −0.11144 ) ( 0.10377−0.11144 ) +6.5453 ( x−0.10377 )( 0.11144−0.10377 )
¿ 17.02738 x−8.31168
Al evaluarla en el punto 0,108 nos queda que la entropía es de 6,472 kj/kg
17. i=corriente, v=voltaje aproxime el valor de v para i=3.5en polinomio de newton
con diferencias divididas
X =i 1 2 3 4
Y =v 120 94 75 62
2 3
P3 ( x ) =f ( x )+ ∆ f ( x ) ( x−x 0 ) +∆ f ( x ) ( x−x 0 )( x −x1 ) + ∆ f ( x ) ( x−x 0 ) ( x−x 1 ) ( x−x 2 )
f ( x 1 )−f ( x 0 )
f ( x1 , x0)=
x 1−x 0
94−120
¿ =−26
2−1
f ( x2 ) −f ( x 1 )
f ( x2 , x1)=
x 2−x 1
75−94
¿ =−19
3−2
f ( x 3 )−f ( x2 )
f ( x 3 , x2 ) =
x 3−x 2
62−75
¿ =−13
4−3
f ( x 2 , x 1 )−f ( x1 , x0 )
f ( x 2 , x 1 , x0 ) =
x 2−x 0
−19+26
¿ =3.5
31
f ( x 3 , x 2 )−f ( x2 , x1 )
f ( x 3 , x2 , x1 ) =
x 3−x 1
−13+19
¿ =3
4−2
f ( x 3 , x 2 , x 1 )−f ( x2 , x1 , x 0 )
f ( x 3 , x2 , x1 ) =
x 3−x 0
3−3.5
¿ −0.167
4−1
x f ( x) ∆ f ( x) ∆ 2 f (x) ∆ 3 f ( x)
1 120 -26 3.5 -0.167
2 94 -19 3
3 75 -13
4 62
P3 ( x ) =120−26 ( x−1 ) +3.5 ( x−1 ) ( x−2 )−0.167 ( x−1 ) ( x−2 ) ( x −3 )
i=x=3.5
P3 ( x =3.5 )=67.8125
18. Aplicando la interpolacion para obtener un polinomio que pasa por los puntos
siguientes (-2,1) (-1,0) (0,1)
P ( x ) =ao + a1 ( x −x o ) +a2 ( x−x o ) ( x−x 1 )
Datos
n=3
i xi f ( x 1) 1 ra 2 da
0 −2 1 −1 1
1 −1 0 1
2 0 1
f (x 1)−f (x o ) f ( x 1 , x 2)−f ( xo , x 1)
1ra f ( x o , x 1)= 2da f (x o , x 1 , x 2)=
x 1−x 0 x 2−x 0
0−1 1−(−1)
¿
−1−2 ¿
0−(−1)
=-1
=1
x Y
-2 1
0−1
=1
-1 0 −1−2
1−(−1)
=1
1−0 0−(−2)
0 1 =1
0−(−1)
x 0=−2
a 0=1
x 1=−1
a 1=−1
a 2=1 x 2=0
P ( x ) =ao + a1 ( x −x o ) +a2 ( x−x o ) ( x−x 1 )
¿ 1±1( x−(−2))+1(x−2)(x−1)
2
¿ 1−x−2+ x +3 x+ 2
2
¿ x + 2 x +1
19. Dada una funcion f ∈¿ 3(l]) y los nodos x 0=a−¿ 1, x 1=a , x=a+ ¿con (1>0, i=1,2 y a
∈ I ¿,se pide:
a.- construir el polinomio de interpolación p(x) de f(x) con los nodos dados
b.= utilizando p(x) obtener la formula
f ( fP bd )−f (f ) f ( f )−f ( f −b c )
f ( a ) a bc P b d e − g
bd bc
a) Construir nuestra tabla
X f(x)
x 0−a−¿ f(a-(1)
f ( f )−f (f −bc )
x 1=a f(a)
f −(f −bc )
f ( fP b d )−f (f )
X=a+( f(a+() f ¿¿¿
fP bd −f
m(x)=f(a-(1) + f ( a )−f ¿¿
20. Con los mismos datos del ejercicio anterior , calcular la expresión de error de la
formula anterior:
njc
f (ε )
un ( x ).∨f ( x)−m n (x )∨. y ( x );
( nP 1 ) ! n
Con y n ( x ) . maxz Π ni=0 ¿
f nP 1( ε ). max∨f nP 1(x )∨, ∀ x ∈[a , b]
En nuestro caso, el error es;
f 3(ε )
u ( x ) .|f ( x )−m ( x )| . y(x);
3!
Con y (x).max|( x−x 0 ) ( x−x 1 ) ( x−x )| , ∀ x ∈ [ x 0 , x ]
f ( ε ) . max|f ( x )|, ∀ x ∈ [ x 0 x ]
3 3
8. Calcular la recta que pasa por los puntos a (−3 ,−2 ) y b ( 3,4 ) .interpolar el valor de la
función para x=2
10. calcular la parábola que pasa por los puntos a (−3,1 ) , b (−1 ,−3 ) y c ( 2,5 ) . Interpolar
el valor de la función x=−2
11. Calcular el polinomio que interpola al conjunto de datos
{(−1,2); (2,8);( 5 ,−3); (8,10)}
12. Calcular usando el polinomio p(x) que interpola a la función f (x)=ln(x ) en los
puntos que resulta de dividir (1,4) en 6 partes iguales
13. Los datos correspondientes al censo de una población se recogen en la siguiente
tabla
Año 1950 1960 1970 1980 1990 2000
¿ habitantes 123.5 131.2 150.7 141.3 203.2 240.5
14. Se ha observado que la vida media, en minutos, de una bacteria varía en función de
la temperatura del medio en el que vive según la siguiente tabla
Temperatura 6° 9° 12 ° 15 ° 16 °
15. Se tiene los siguientes datos sobre la evaluación del índice de precios al consumo
(IPC) en los últimos años