Polinomios de Lagrange

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 5

TECNOLOGICO DE ESTUDIOS SUPERIORES DE ECATEPEC

DIVISION DE INGENIERIA QUIMICA Y BIOQUIMICAPROGRAMACION Y METODOS NUMERICOS


GRUPO: 3401. HERNANDEZ-DIAZ-LUCIA-GUADALUPE.

LOS POLINOMIOS DE LAGRANGE.


Polinomio interpolante de Lagrange de primer grado.
El método de los polinomios de Lagrange utiliza una estructura de construcción
para el polinomio de interpolación, basada en la combinación creciente de
acumuladores aritméticos y geométricos. El polinomio de interpolación de grado n
es el siguiente:
n
f n ( x )=∑ f ( x j ) Lnj ( x )
j=0

En la relación anterior se tiene la presencia de un acumulador geométrico Lnj ( x )


para el cual:
n
( x−x k )
Lnj ( x ) = ∏
k=0 ,k ≠ j ( x j −x k )

Para construir un polinomio de Lagrange de grado 1 se elabora el siguiente


desarrollo:

f 1 ( x )=f ( x0 ) L10 ( x )+ f ( x 1) L11 ( x )


1 1
( x−x k ) ( x−x )
f 1 ( x )=f (x 0) ∏ ( x −x ) + f ( x 1 ) ∏ x −xk
k=0 , k≠ 0 0 k k=0 ,k ≠ 1 ( 1 k)

( x −x1 ) (x −x0 )
f 1 ( x )=f ( x0 ) + f ( x ¿¿ 1) ¿
( x 0−x 1 ) (x 1−x 0)

Ejemplo 1. Calcular el polinomio interpolante de Lagrange de grado 1 para los dos


puntos P0(x0, y0) y P1(x1, y1), donde:
x0=2, y0=f(x0) =4
x1=4, y1=f(x1) =6
solución analítica.
Al sustituir los datos en la fórmula, se elabora el siguiente desarrollo:

( x−x 1 )
f 1 ( x )=f (x ¿¿ 0) ¿
( x− x0 )
( x ¿ ¿ 0−x1 )+f ( x¿¿ 1) ¿¿
( x 1−x 0 )
( x−4 ) ( x−2 ) ( 4 ) (6)
f 1 ( x ) =( 4 ) +( 6 ) = ( x−4 ) + ( x−2 ) =−2 x +8+ 3 x−6=x +2
( 2−4 ) ( 4−2 ) (−2 ) (2)

Polinomios de Lagrange Página 1 de 5 M en C Jose Luis Parra.


TECNOLOGICO DE ESTUDIOS SUPERIORES DE ECATEPEC
DIVISION DE INGENIERIA QUIMICA Y BIOQUIMICAPROGRAMACION Y METODOS NUMERICOS
GRUPO: 3401. HERNANDEZ-DIAZ-LUCIA-GUADALUPE.
Como resultado se tiene la función reconstruida de primer grado: f 1 ( x )=x +2.
Polinomio interpolante de Lagrange de segundo grado.
Para construir un polinomio de Lagrange de grado 2 se desarrolla la siguiente
sumatoria:
2
f 2 ( x )=∑ f ( x ¿¿ j)L2j (x )¿
j=0

f 2 ( x )=f ( x 0) L20 ( x )+ f ( x 1) L21 ( x ) + f ( x 2) L22 (x)

Al incluir los acumuladores geométricos L2j se tiene la siguiente expresión:


2 2 2
( x−x k ) (x−x k ) (x−x k )
f 2 ( x )=f (x ¿¿ 0) ∏ + f (x ¿¿ 1) ∏ +f (x¿¿ 2) ∏ ¿¿¿¿
k=0 , k ≠0 (x 0−x k ) k=0 , k≠ 1 (x 1−x k ) k=0 ,k ≠2 ( x ¿ ¿ 2−x k )

( x−x 1 )( x− x2 )
f 2 ( x )=f (x ¿¿ 0)= ¿
( x−x 0 ) ( x−x 2 ) ( x−x 0 ) ( x−x 1 )
( x ¿ ¿ 0−x 1 )(x 0−x 2)+ f (x¿ ¿1) + f (x 2) ¿¿
( x 1−x 0 ) ( x 1−x 2) ( x 2−x 0 ) ( x 2−x 1 )
Ejemplo 2. Calcular el polinomio interpolante de Lagrange de segundo grado para
los tres puntos P0 (x0, y0), P1 (x1, y1) y P2 (x2, y2) donde se tienen los siguientes
valores:
x0=1, y0=f(x0) =4
x1=2, y1=f(x1) =-1
x2=3, y2=f(x2) =2
solución analítica.
Al sustituir los datos en la formula, se elabora el siguiente desarrollo:

( x−x 1 )( x−x 2)
f 2 ( x )=f (x ¿¿ 0) ¿
( x−x 0 ) ( x−x 2 ) ( x−x 0) ( x− x1 )
( x ¿ ¿ 0−x1 )(x 0−x 2)+ f (x ¿¿ 1) + f (x 2) ¿¿
( x 1−x 0 ) ( x 1−x 2 ) ( x2 −x 0) ( x 2−x 1 )
( x −2 )( x−3 ) ( x −1 )( x−3 ) ( x−1 ) ( x−2 )
f 2 ( x )=(4) +(−1) +( 2)
( 1−2 )( 1−3 ) ( 2−1 )( 2−3 ) ( 3−1 ) ( 3−2 )

(x ¿¿ 2−5 x+ 6) ( x2 −4 x +3 )
f 2 ( x )=(4) +(−1) ¿
(−1 ) (−2 ) ¿¿

f 2 ( x ) ¿(2)( x¿¿ 2−5 x +6)+(1)( x ¿¿ 2−4 x+ 3)+(1)( x¿¿ 2−3 x +2)¿ ¿ ¿

Polinomios de Lagrange Página 2 de 5 M en C Jose Luis Parra.


TECNOLOGICO DE ESTUDIOS SUPERIORES DE ECATEPEC
DIVISION DE INGENIERIA QUIMICA Y BIOQUIMICAPROGRAMACION Y METODOS NUMERICOS
GRUPO: 3401. HERNANDEZ-DIAZ-LUCIA-GUADALUPE.
f 2 ( x )=2 x 2−10 x+12+ x 2 −4 x +3+ x 2−3 x +2

Como resultado se tiene la función reconstruida de segundo grado:


f 2 ( x )=4 x 2−17 x+ 17

Polinomio interpolante de Lagrange de tercer grado.


Para construir un polinomio de Lagrange de tercer grado se desarrolla la siguiente
sumatoria:
3
f 3 ( x )=∑ f ( x j ) L3j ( x )
j=0

f 3 ( x )=f (x ¿¿ 0)L30 (x )+ f ( x ¿¿ 1) L31 ( x)+ f ( x ¿¿ 2) L32 ( x)+ f ( x¿ ¿3) L33 (x)¿ ¿ ¿¿


3
(x−x k )
f 3 ( x )=f (x ¿¿ 0) ∏ 3
k =0 , k ≠0 ( x−x k )
(x ¿ ¿ 0−x k )+ f (x ¿¿ 1) ∏ 3
k=0 ,k ≠ 1 (x−x k )
(x ¿ ¿ 1−x k )+ f ( x¿ ¿2) ∏ 3
k=0 ,k ≠2
( x ¿ ¿ 2−x k )+ f (x¿ ¿3) ∏
k=0 ,k ≠ 3 (

( x−x 1)( x−x 2 )(x−x 3) ( x−x 0)( x −x2 )(x−x 3)


f 3 ( x )=f (x ¿¿ 0) + f (x ¿¿1) + f (x ¿¿ 2) ¿ ¿ ¿ ¿¿
( x 0−x 1 )(x 0−x 2)(x 0−x 3 ) ( x1 −x0 )( x 1−x 2)(x 1−x 3 )

El desarrollo anterior se puede resolver por el método tradicional de sustitución de


valores y simplificación algebraica, o adaptando la solución a la hoja de cálculo de
Excel. Los cuatro denominadores se asignan a cuatro términos Den0, Den1, Den2 y
Den3:
Den0=( x0 −x 1)(x 0−x 2)( x 0− x3 ) Den1=(x 1−x 0 )(x 1−x 2)( x 1−x 3)

Den2=(x 2−x 0 )(x 2−x 1)( x 2−x 3) Den3=( x3 −x0 )( x 3−x 1)(x 3−x 2 )

Estos términos se incluyen en el polinomio de Lagrange como se muestra a


continuación:
f ( x0 ) (x¿¿ 1) (x ¿¿ 2) f ( x ¿¿ 3)
f 3 ( x )= ( x−x 1 )( x−x 2) ( x−x 3 ) + f ( x−x 0 ) ( x−x 2 ) ( x−x 3 ) +f ( x−x 0 ) ( x−x 1 ) ( x−x 3 ) + (
Den0 Den1 Den 2 Den3

Al desarrollar los elementos dentro de los paréntesis se obtiene la siguiente


expresión:

(x¿ ¿0) 3 f ( x ¿¿ 1) 3
f 3 ( x )=f [ x + ( −x1 −x2− x3 ) x 2+ ( x1 x 2+ x 1 x3 + x 2 x 3 ) x+ x1 x 2 x 3 ] + x + ( −x 0−x 2−x 3 ) x 2 + ( x 0 x 2 +
Den0 Den1 [

Polinomios de Lagrange Página 3 de 5 M en C Jose Luis Parra.


TECNOLOGICO DE ESTUDIOS SUPERIORES DE ECATEPEC
DIVISION DE INGENIERIA QUIMICA Y BIOQUIMICAPROGRAMACION Y METODOS NUMERICOS
GRUPO: 3401. HERNANDEZ-DIAZ-LUCIA-GUADALUPE.

Se establecen las siguientes equivalencias:


b 0=−x 1−x 2 −x3 c 0=x 1 x2 + x 1 x 3 + x 2 x 3 d 0=x 1 x2 x 3

b 1=−x 0−x 2 −x3 c 1=x 0 x2 + x 0 x 3 + x 2 x 3 d 1=x 0 x2 x 3

b 2=−x 0−x 1−x 3 c 2=x 0 x1 + x 0 x 3 + x 1 x 3 d 2=x 0 x1 x 3

b 3=x 0−b x 1−x 2 c 3=x 0 x1 + x 0 x2 + x 1 x 2 d 3=x 0 x 1 x 2

Se puede elaborar el siguiente desarrollo con las relaciones anteriores:

f (x¿ ¿0) 3 f (x ¿¿1) 3 f (x¿ ¿2) 3 f ( x ¿¿ 3)


f 3 ( x )= [ x +b 0 x 2 +c 0 x + d0 ]+ [ x +b1 x 2+ c1 x + d1 ] + [ x +b2 x 2+ c2 x+ d 2 ]+
Den0 Den1 Den2 Den3

Reacomodando términos comunes se obtiene la siguiente expresión polinomial:

f 3 ( x )=¿

Adaptando la expresión anterior se tiene la siguiente tabla de equivalencias para


los parámetros a, b, c y d; del polinomio interpolante de tercer grado
f3(x)=ax3+bx2+cx+d:
a =f(x0) /Den0+f(x1) /Den1+f(x2) /Den2+f(x3) /Den3
b =f(x0) b0/Den0+f(x1) b1/Den1+f(x2) b2/Den2+f(x3) b3/Den3
c =f(x0) c0/Den0+f(x1) c1/Den1+f(x2) c2/Den2+f(x3) c3/Den3

Polinomios de Lagrange Página 4 de 5 M en C Jose Luis Parra.


TECNOLOGICO DE ESTUDIOS SUPERIORES DE ECATEPEC
DIVISION DE INGENIERIA QUIMICA Y BIOQUIMICAPROGRAMACION Y METODOS NUMERICOS
GRUPO: 3401. HERNANDEZ-DIAZ-LUCIA-GUADALUPE.
d =f(x0) d0/Den0+f(x1) d1/Den1+f(x2) d2/Den2+f(x3) d3/Den3

Ejemplo 3. Calcular el polinomio interpolante de Lagrange de tercer grado para


los cuatro puntos P0 (x0, y0), P1 (x1, y1), P2(x2, y2) y P3 (x3, y3) donde se tienen los
siguientes valores:
X0=1; y0=f(x0) =-1 X2=5 y2=f(x2) =-2

X1=3; y1=f(x1) =2 X3=7 y3=f(x3) =3


Solución analítica:
El polinomio interpolante de Lagrange de tercer grado se obtiene con el método
tradicional, mediante el siguiente desarrollo:
( x−x 1)( x−x 2 )(x−x 3) ( x−x 0)( x −x2 )(x−x 3)
f 3 ( x )=f (x ¿¿ 0) + f (x ¿¿1) + f (x ¿¿ 2) ¿ ¿ ¿ ¿¿
( x 0−x 1 )(x 0−x 2)(x 0−x 3 ) ( x1 −x0 )( x 1−x 2)(x 1−x 3 )

( 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 3 ( x )=(−1 ) +( 2 ) + (−2 ) +(3)
( 1−3 ) ( 1−5 ) ( 1−7 ) ( 3−1 )( 3−5 ) ( 3−7 ) ( 5−1 )( 5−3 ) (5−7 ) (7−1)(7−3)(7−5)

( x2 −8 x+15 ) ( x−7 ) ( x2 −6 x+5 ) ( x−7 ) ( x 2−4 x +3 ) ( x−7 ) ( x2 −4 x +3 ) ( x−5 )


f 3 ( x )=(−1 ) +( 2 ) + (−2 ) +( 3 )
(−2 )(−4 )(−6 ) ( 2 ) (−2 )(−4 ) ( 4 )( 2 ) (−2 ) ( 6 )( 4 ) ( 2 )

(x ¿ ¿ 3−15 x 2+ 71 x−105) ( x 3−13 x2 + 49 x −35 ) ( x 3−11 x 2+ 31 x−21 ) (x ¿ ¿ 3−9 x 2 +


f 3 ( x )=(−1 ) + (2) + (−2 ) +( 3 )
(−48 ) ( 16 ) (−16 ) ( 48 )
Al simplificar la ecuación anterior, ésta se reduce a la siguiente:
f3(x)=0.3333x3-3.785x2+12.6667x-10.125
solución con la hoja de cálculo de Excel para un polinomio de tercer grado:

La coincidencia de los datos, indica que el polinomio interpolante de Lagrange de


tercer grado es aceptable, ya que cumple con los datos de entrada.

Polinomios de Lagrange Página 5 de 5 M en C Jose Luis Parra.

También podría gustarte