Cálculo Numérico.

Mínimos Cuadrados Continuo.

Polinomios Trigonométricos
Prof. Jorge Omar

Junio 2020

1 Mínimos Cuadrados Continuo. Polinomios Trigonométri-

cos.

Supongamos que el conjunto de funciones f 0 ; 1 ; :::; ng es linealmente en [a; b] : Dada la fun-

ción f; buscamos una combinacion lineal
X
n
Pn (x) = ak k
k=1

tal que !2
Z b X
n
E (a0 ; a1 ; :::; an ) = f (x) ak k dx
a k=1

sea mínimo. Los coe…cientes se obtienen resolviendo el sistema de ecuaciones normales

@E (a0 ; a1 ; :::; an )
= 0; (j = 0; 1; 2; :::; n)
@aj
Z !
b X
n
2 f (x) ak k j dx = 0; (j = 0; 1; 2; :::; n)
a k=1

1
 X
n Z b Z b
ak k (x) j (x) dx = j (x) f (xi ) dx; (j = 0; 1; 2; :::; n)
k=1 a a

Si f 0 ; 1 ; :::; ng es una base ortogonal en [a; b] ; es decir,

8
Z b < 0; k 6= j
k ; j = k (x) j (x) dx =
a : > 0; k = j j

P
n
Entonces los coe…cientes del polinomio Pn (x) = ak k; que mejor aproxima a la función
k=1
f (x) en [a; b] ; en el sentido del método de mínimos cuadrados, se obtienen con la fórmula
Rb
hf; k i f (x) k (x) dx
ak = = a Rb 2 ; k = 0; 1:::; n
h k; ki k (x) dx a

Las funciones

0 (x) = 1; 1 (x) = cos (x) ; 2 (x) = sin (x) ; :::

3 (x) = cos (2x) ; 4 (x) = sin (2x) ; :::

R
son ortogonales en [ ; ] ; es decir k; j = k (x) j (x) dx = 0; para k 6= j:
El conjunto f 0 ; 1 ; :::; ng forman una base ortogonal del subespacio generado
Dada la función f (x) de…nida en [ ; ] ; vamos a obtener la aproximación de mínimos
cuadrados continua de la forma

a0 X
n
Sn (x) = + (ak cos (kx) + bk sin (kx))
2 k=1

A Sn (x) se le llama polinomio trigonométrico de grado menor o igual que n

Los coe…cientes se calculan de manera que se satisfaga la condición
Z
2
E = kf (x) Sn (x)k = [f (x) Sn (x)]2 dx ! min

@E @E
es decir, = 0; k = 0; 1; :::; n y = 0; k = 1; :::; n 1
@ak @bk
Coe…cientes de Fourier:

hf (x) ; cos (kx)i hf (x) ; sin (kx)i

ak = ; bk =
hcos (x) ; cos (kx)i hsin (x) ; sin (kx)i

2
Ejemplo 1 Determine el polinomio trigonométrico Sn (x) para la función f (x) = jxj ; x 2
[ ; ]

Solución:
1R 1 R0 1R 2R
a0 = jxj dx = xdx + 0
xdx = 0
xdx =

1R 2R 2 h i
ak = jxj cos (kx) dx = 0
x cos (kx) dx = ( 1)k 1
k2
1R
bk = jxj sin (kx) dx = 0

2 Pn ( 1)k 1
Sn (x) = + cos (kx)
2 k=1 k2

5
5 y
y
4
4

3
3

2 2

1 1

-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
x x

n=1 n=5

Ejercicios:
Determinar el polinomio trigonométrico Sn (x) para la función f (x) ; con x 2 [ ; ]:
Gra…que
(a) f (x) = 8
x
< x; si x<0
2
(b) f (x) =
: + x; si 0 x <
8 2
< 0; si x<0
(c) f (x) =
: x; si 0 x <

3
 Caso Discreto.
Supongamos que tenemos un conjunto de N+1 puntos equiespaciados en el interalo [ ; ]

2
xi = + kh; h= ; i + 0; 1; :::; N
N

Si 2n < N; entonces existe

a0 X
n
Sn (x) = + (ak cos (kx) + bk sin (kx))
2 k=1

tal que

2 X
N
aj = f (xk ) cos (kx) ; j = 0; 1; 2; :::; n
N k=1

2 X
N
bj = f (xk ) sin (kx) ; j = 1; 2; :::; n
N k=1

Ejemplo 2 Dada la función f (x) = x2 ; x 2 ( ; ] y periódica de período 2 : Encontrar el

polinomio trigonométrico discreto de grado 5 usando una malla de 12 puntos

y 1.4
1.2
1.0
0.8
0.6
0.4
0.2

-5 -4 -3 -2 -1 1 2 3 4 5
-0.2 x
-0.4
-0.6
-0.8
-1.0
-1.2
-1.4

2 k
Solución: xk = + 12
; k = 1; :::; 12

4
 2 k + 212k
(xk ; f (xk )) = + 12
; 2
; k = 1; :::; 11; f( )=0
f (x) es impar ) aj = 0; j = 0; 1; 2; :::; 12
1
P11
b1 = 12 xk sin (xk ) = 0:977 049
k=1
1
P11
b2 = 12
xk sin (2xk ) = 0:453 450
k=1
1
P11
b3 = 12
xk sin (3xk ) = 0:261 799
k=1
1
P11
b4 = 12
xk sin (4xk ) = 0:151 150
k=1
1
P11
b5 = 12
xk sin (5xk ) = 0:07014 89
k=1
S5 (x) = 0:977 049 sin (x) 0:453 450 sin (2x)+0:261 799 sin (3x) 0:151 150 sin (4x)+0:07014 89 sin (5x)

y
1.2
1.0
0.8
0.6
0.4
0.2

-5 -4 -3 -2 -1 1 2 3 4 5
-0.2 x
-0.4
-0.6
-0.8
-1.0
-1.2