Programación y Métodos Numéricos.
Héctor Andrés Granada Díaz
Universidad Nacional de Colombia - Sede Manizales
Abril de 2010
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 1 / 17
¿Existe una relación entre los polinomios de interpolación encontrados con
nodos diferentes en un mismo intervalo y para una misma función y el
error entre ellos y la función original?
La respuesta a este problema fue encontrada por P.L.Chebyshev.
Existe una con…guración de nodos en el intervalo [a, b ] tal que el error
entre el polinomio y la función a interpolar es mínimo entre todas las
con…guraciones de nodos posibles en el intervalo.
Para esto necesitamos conocer algunas propiedades de los polinomios de
Chebyshev.
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 2 / 17
Polinomios de Chebyshev
Inicialmente trabajaremos en el intervalo [ 1, 1] y después realizando un
cambio de variable se generalizara a intervalos [a; b]
1. Estos polinomios se generan de forma recursiva tomando
T0 ( x ) = 1
T1 ( x ) = x
y
Tk (x ) = 2xTk 1 (x ) Tk 2 (x ) para k = 2, 3,
2. El coe…ciente de x N en TN (x ) es 2N 1 para N 1.
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 3 / 17
3. Cuando N = 2M, el polinomio T2M es una función par
T2M ( x ) = T2M (x )
Cuando N = 2M + 1, el polinomio T2M +1 es una función impar
T2M +1 ( x ) = T2M +1 (x )
4. En el intervalo [ 1, 1] poseen representación trigonométrica
1
TN (x ) = cos(N cos (x )) para x 2 [ 1, 1]
5. TN (x ) posee N ceros distintos xk en el intervalo [ 1, 1] dados por
(2k + 1)π
xk = cos para k = 0, 1, 2, N 1
2N
Estos valores se denominan nodos de Chebyshev
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 4 / 17
6. Propiedades de Ortogonalidad
1
N
∑ Ti (xk )Tj (xk ) = 0 cuando i 6= j
k =0
2
N
N +1
∑ T i ( xk ) T j ( xk ) = 2
cuando i = j 6= 0
k =0
3
N
∑ T i ( xk ) T j ( xk ) = N + 1 cuando i = j = 0
k =0
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 5 / 17
Teorema
Supongamos que deseamos aproximar una función en el intervalo [ 1, 1] .
la elección de los nodos que determina el polinomio de interpolación que
mejor se aproxima a la función viene dada por los Nodos de Chebyshev.
Es decir, si deseamos aproximar una función en [ 1, 1] por un polinomio
tomando N muestras de tal manera que el polinomio obtenido al interpolar
sea el que mejor aproxime la función, debemos escoger los Nodos de
Chebyshev.
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 6 / 17
Cambio de Intervalo
Si deseamos aproximar una función en un intervalo [a, b ] utilizando los
nodos de chebyshev, hacemos el cambio de variable
b a a+b
x= t+
2 2
o
x a
t=2 1
b a
tomando a x b y 1 t 1.
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 7 / 17
Los Nodos de Chebyshev en el intervalo [ 1, 1] para N + 1 puntos vienen
dados por
(2k + 1)π
tk = cos para k = 0, 1, ,N
2N + 2
Por lo tanto haciendo el cambio de variable correspondiente, obtenemos
que los Nodos de Chebyshev en el intervalo [a, b ] vienen dados por
b a a+b
xk = tk + para k = 0, 1, ,N
2 2
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 8 / 17
Teorema de aproximación
El polinomio de aproximación de Chebyshev PN (x ) de grado menor o igual
que N para una función f en [ 1, 1] puede escribirse como una suma
ponderada de los polinomios fTj (x )g
N
PN (x ) = ∑ cj T j ( x )
j =1
donde los coe…cientes cj vienen dados por:
N N
1 1
N + 1 k∑ N + 1 k∑
c0 = f ( xk ) T (
0 kx ) = f (xk )
=0 =0
N N
2 2 j π (2k + 1)
N + 1 k∑ N + 1 k∑
cj = f ( xk ) T (
j kx ) = f (xk )
=0 =0 2N + 2
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 9 / 17
Ejemplo
Consideremos la función
1
f (x ) =
(1 + 12x 2 )
Aproximemos la función en el intervalo [ 1, 1] con 10 nodos de Chebyshev
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 10 / 17
Utilizando 20 nodos de Chebyshev
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 11 / 17
Recta de regresión en mínimos cuadrados
Si tenemos un conjunto de puntos fxk , yk gNk =1 donde todas las abscisas
son distintas los cuales son producto de un muestreo experimental (errores
inherentes) entonces la recta de regresión o recta optima (en el sentido de
los mínimos cuadrados) es la recta de ecuación
y = Ax + B
donde los coe…cientes de la recta son la solución del sistema lineal,
conocido como las ecuaciones normales de Gauss
! !
N N N
∑ xk2 A+ ∑ xk B = ∑ xk yk
k =1 k =1 k =1
!
N N
∑ xk A + NB = ∑ yk
k =1 k =1
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 12 / 17
Interpolación Lineal a trozos
Si tenemos un conjunto de puntos fxk , yk gNk =0 La interpolación polinomial
proporciona un buen método para encontrar una función que interpole los
datos, pero la función resultante es de carácter polinomial las cuales
pueden presentar oscilaciones grandes, por eso se busca un método para
aproximar los datos intervalo a intervalo.
Una posible solución es buscar funciones que interpolen solo un par de
datos y después de…nir la función total como una función a trozos.
8 9
>
> f1 (x ) si x 2 [x0 , x1 ] >
>
< f2 (x ) si x 2 [x1 , x2 ] >
> =
f (x ) = ..
>
> . >
>
>
: >
;
fN (x ) si x 2 [xN 1 , xN ]
Además, para cada i 2 f1, 2, , Ng
fi (xi 1) = yi 1
fi (xi ) = yi
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 13 / 17
Con lineas rectas
x 1 3 5 7 9
y 1 10 3 5 1
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 14 / 17
Splines Cúbicos
Supongamos que tenemos N + 1 puntos fxk , yk gN k =0 cuyas abscisas están
ordenadas de forma creciente. Una función s (x ) es una cercha cubica
interpoladora para los datos, si existen N polinomios cúbicos Sk (x ) que
se pueden escribir en términos de unos coe…cientes
sk ,0 , sk ,1 , sk ,2 y sk ,3
Como
Sk (x ) = sk ,0 + sk ,1 (x xk ) + sk ,2 (x xk )2 + sk ,3 (x xk )3
Así
1. Dado x 2 [xk , xk +1 ]
S (x ) = Sk (x ) = sk ,0 + sk ,1 (x xk ) + sk ,2 (x xk )2 + sk ,3 (x xk )3
donde k = 0, 1, N 1 que veri…can las propiedades:
2. La función Interpola los datos
S (xk ) = yk para k = 0, 1, N
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 15 / 17
3. La función es derivable y con derivada continua
Sk (xk +1 ) = Sk +1 (xk +1 ) para k = 0, 1, N 2
4.
0 0
Sk (xk +1 ) = Sk +1 (xk +1 ) para k = 0, 1, N 2
5. La segunda derivada existe
00 00
Sk (xk +1 ) = Sk +1 (xk +1 ) para k = 0, 1, N 2
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 16 / 17
x 1 3 5 7 9
y 1 10 3 5 1
Héctor Granada (Universidad Nacional de Colombia
DEPARTAMENTO
- Sede Manizales)
DE MATEMATICAS Y ESTADISTICA 04/10 17 / 17
