Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

3 2 3 FormulasdeEuler

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 10

Datos de catalogación bibliográfica

Ecuaciones diferenciales
Isabel Carmona Jover, Ernesto Filio López

Quinta edición

Pearson Educación, México, 2011


ISBN: 978-607-32-0206-0
Área: Matemáticas
Formato: 20 × 25.5 cm. Páginas: 536

Editor: Rubén Fuerte Rivera


e-mail: ruben.fuerte@pearsoned.com
Editor de desarrollo: Felipe Hernández Carrasco
Supervisor de producción: Juan José García Guzmán

QUINTA EDICIÓN, 2011

D.R. © 2011 por Pearson Educación de México, S.A. de C.V.


Atlacomulco 500, 5° piso
Col. Industrial Atoto, CP 53519
Naucalpan de Juárez, Edo. de México

Cámara Nacional de la Industria Editorial Mexicana, Reg. Núm. 1031

Reservados todos los derechos. Ni la totalidad ni parte de esta publicación pueden reproducirse, registrarse o transmitirse,
por un sistema de recuperación de información, en ninguna forma ni por ningún medio, sea electrónico, mecánico, foto-
químico, magnético o electroóptico, por fotocopia, grabación o cualquier otro, sin permiso previo por escrito del editor.

El préstamo, alquiler o cualquier otra forma de cesión de uso de este ejemplar requerirá también la autorización del editor
o de sus representantes.

ISBN VERSIÓN IMPRESA: 978-607-32-0206-0


ISBN VERSIÓN E-BOOK: 978-607-32-0207-7
ISBN E-CHAPTER: 978-607-32-0208-4

Impreso en México. Printed in Mexico.


1 2 3 4 5 6 7 8 9 0 - 13 12 11 10
Addison-Wesley
es una marca de
442 Capítulo 8 Series de Fourier

p
1p
→ ∫ f ( x )sen nxdx = bn p ∴ bn = ∫ f ( x )sen nxdx , n = 1, 2, 3, ...
p −p
−p

Las fórmulas así definidas se llaman fórmulas de Euler:


1 p
a0 = ∫ f ( x )dx
2p − p
1p
an = ∫ f ( x ) cos nxdx , n = 1, 2, 3, ...
p −p
1p
bn = ∫ f ( x )sen nxdx , n = 1, 2, 3, ...
p −p

Definición 8.2
Serie de Fourier. La función:
f ( x ) = a0 + a1 cos x + b1sen x + a2 cos 2 x + ...

= a0 + ∑ (an cos nx + bn sen nx )
n =1
se llama serie de Fourier y los coeficientes obtenidos a partir de a0, an, a0, an
y bn (las fórmulas de Euler para n = 1, 2, 3, ... se llaman coeficientes de
Fourier de f ( x ).

EJEMPLO 1
Hallar la serie de Fourier de la siguiente función periódica con periodo 2p y
trazar la gráfica de las tres primeras sumas parciales.
y

1 2

x
p p p p
2 2

Figura 8-4.

Paso 1. Encontramos los coeficientes de Fourier, mediante las fórmulas de


Euler.
1 p
a0 = ∫ f ( x )dx
2p − p
−p p 
1  2 1 0 2 p
1 
= ∫
2p  − p 2
dx + ∫ 0 dx + ∫ dx + ∫ 
p 2
dx
−p 2 0 2 
Fórmulas de Euler 443

1 x
−p 2
x
p  1
=  +x 02 + = .
p

2p  2 2 p  2
−p 2

1
∴ a0 = .
2
1p
an = ∫ f ( x ) cos nxdx
p −p
−p p 
1 2 1 2 p
1
= ∫
p  − p 2
cos nxdx + ∫ cos nxdx + ∫
p 2
cos nxdx 

0  2

1 1
−p/ 2 p  p
1 2 1
=  sen nx + sen nx + sen nx 
p  2n n 2 p 
−p 0 2

1 1 np 1 np 1 np 
= − sen + sen − sen  = 0
p  2n 2 n 2 2n 2 
∴ an = 0.
1p
bn = ∫ f ( x )sen nxdx
p −p
−p p 
1 2 1 2 p
1
bn = ∫
p  − p 2
sen nxdx + ∫ sen nxdx + ∫ sen nxdx 
p 2 
0  2

1 1 
−p/ 2 p p
1 2 1
= − cos nx − cos nx − cos nx 
p  2n n 2n p 
−p 0 2

1 1 np 1 1 np 1 1 1 np 
= − cos + cos np − cos + − cos np + cos 

p 2n 2 2n n 2 n 2n 2n 2 
1 1 np 1 
= − cos +
p n 2 n 
 1
 np n = 1, 3, 5, 7, 9, 11, ...

1  np   2
= 1 − cos  =  n = 2, 6, 10, ...
np  2   np
0 n = 4, 8, 12, ...

Paso 2. Sustituimos los coeficientes de Fourier en la serie:



f ( x ) = a0 + ∑ bn sen nx
n =1

= a0 + b1sen x + b2sen 2 x + b3sen 3 x + ...


1 1 2 1 1
f ( x ) = + sen x + sen 2 x + sen 3 x + 0 + sen 5 x + ...
2 p 2p 3p 5p
444 Capítulo 8 Series de Fourier

1 1 1 1 
= + sen x + sen 2 x + sen 3 x + sen 5 x + ... .
2 p 3 5 

Paso 3. Graficamos S1, S2 y S3

1 1 1 1 1 1
s1 = , s2 = + sen x , s3 = + sen x + sen 2 x
2 2 p 2 p p

0
p 0.75p 0.5p 0.25p 0 0.25p 0.5p 0.75p p

Figura 8-5.

x −p 3p p p 0 p p 3p p
− − −
4 2 4 4 2 4

s2 1 0.28 0.2 0.3 1 0.72 0.8 0.72 1


2 2 2

s3 1 0.6 0.2 0.02 1 1.04 0.8 0.4 1


2 2 2

Observamos que cada suma parcial se aproxima más a la función original y


en el infinito coincide exactamente. Por ello, la serie (si es convergente)
converge a f ( x ).

EJEMPLO 2
Hallar la serie de Fourier de la siguiente función:

x2
f (x) = , −p < x < p, T = 2p
2
1 p x2
a0 = ∫ dx
2p − p 2
1 p x2
an = ∫ cos nxdx , n = 1, 2, 3, ...
p −p 2
Fórmulas de Euler 445

1 p x2
bn = ∫ sen nxdx , n = 1, 2, 3, ...
p −p 2
p
1 x3 1 p2
a0 =
4p 3
=
12p
( p 3
− (− p) 3
) 6
=
−p

1  2x 
p
p
 x2 2 
an =  2 cos nx +  − 3  sen nx 
2p  n −p n n   −p

1  2p 2p 
=  cos np + 2 cos np (vea página 433)
2p  n 2
n 

 2
1  − n 2 n = 1, 3, 5, 7, ...
= 2 cos np = 
n 2 n = 2, 4, 6, 8, ...
 n 2

1  2x 
p
p
 2 x2 
bn =  2 sen nx +  3 −  cos np
2p  n −p n n  −p

1  2 p2   2 p2  
=  3 −  cos n p −  3 −  cos np = 0
2p  n n  n n  

p2 ∞
→ f (x) = + ∑ an cos nx
6 n =1

p2 2 2
= + (−2) cos x + cos 2 x − cos 3 x + ...
6 4 9


p2 (−1)n
∴ f (x) = + 2∑ 2 cos nx
6 n =1 n

EJEMPLO 3
Hallar la serie de Fourier correspondiente a la función:

 p p
0 −p< x <− y <x<p
2 2
f (x) = 
cos x −
p
<x<
p
 2 2
446 Capítulo 8 Series de Fourier

T = 2p
1 p
2 a0 = ∫ f ( x )dx
2p − p
1p
1
an = ∫ f ( x ) cos nxdx;
p −p
0 n = 1, 2, 3, ...
p 0.75p0.5p 0.25p 0 0.25p 0.5p 0.75p p
1p
Figura 8-6.
bn = ∫ f ( x )sen nxdx;
p −p
n = 1, 2, 3, ...

−p p 
1  2 2 p
a0 = ∫
2p  − p
0 dx + ∫ cos xdx + ∫ 
0 dx
−p −p 2 2

p
1 2 1 1
= sen x = [1 − (−1)] = .
2p −p 2 2p p

p
2
1
an =
p
∫ cos x cos nxdx
−p 2

1
Como cos x cos y = [ cos( x − y) + cos( x + y)]
2

p
2
1
→ an =
2p
∫ [ cos( x − nx ) + cos( x + nx )] dx
−p 2

p
2
1
=
2p
∫ [ cos(1 − n) x + cos(1 + n) x ] dx
−p 2

p
1  1 1  2
=  sen (1 − n) x + sen (1 + n) x 
2 p 1 − n 1+ n  −p
2

1  1 p 1 p 1 p 1 p
=  sen (1 − n) + sen (1 − n) + sen (1 + n) + sen (1 + n) 
2 p 1 − n 2 1− n 2 1− n 2 1+ n 2

1  2 p 2 p
=  sen (1 − n) + sen (1 + n) 
2 p 1 − n 2 1+ n 2
Fórmulas de Euler 447

1 1 p 1 p
=  sen (1 − n) + sen (1 + n)  , n ≠ 1
p 1 − n 2 1+ n 2

1 1  p p p p 1  p p p p 
=  sen cos n − cos sen n  + sen cos n + cos sen n 
p 1 − n  2 2 2 2  1+ n  2 2 2 2 
1 p 1 1 
= cos n  + 
p 2 1 − n 1 + n 
 +2
 p(n 2 − 1) n = 2, 6, ...
−2 2 
= cos n =  0 n = 3, ...
p(n 2 − 1) p 
 −2 n = 4, 8, ...
 p(n 2 − 1)

Para n = 1 vemos:
p
2
1
a1 =
2p
∫ (cos 0 x + cos 2 x )dx
−p 2

p
2
1
=
2p
∫ (1 + cos 2 x )dx
−p 2

p
1  1  2 1 p p 1
=  x + sen 2 x  =  + = .
2p  2  −p 2 p  2 2  2
2
p
2
1
Ahora buscamos: bn =
p
∫ cos x sen nxdx
−p 2

p
1  1 1  2
= − cos( n − 1) x − cos( n + 1) x = 0.
2p  n − 1 n +1 
−p 2

1 1 2 2 2
→ f (x) = + cos x + cos 2 x − cos 4 x + cos 6 x − ...
p 2 3p 15p 35p
1 1 2 1 1 1 
= + cos x +  cos 2 x − cos 4 x + cos 6 x − ... .
p 2 p 3 15 35 

EJEMPLO 4
Hallar la serie de Fourier de:

f (x) = x, −p < x < p

y graficar las tres primeras sumas parciales.


448 Capítulo 8 Series de Fourier

p
1
2p −∫p
a0 = xdx

p
1 x2 1
= ?
2p 2
=
4p
( p2 − (−p)2 ) = 0.
−p

1p
an = ∫ x cos nxdx , n = 1, 2, 3, ...
p −p

1 1 
p
x
=  2 cos nx + sen nx  (vea página 433)

p n n  −p

1 1 1 p p 
=  cos np − 2 cos np + sen np − sen np = 0
p n 2
n n n 

1p
bn = ∫ x sen nxdx , n = 1, 2, 3, ...
p −p

1 1 
p
x
=  sen nx − cos nx 

p n 2
n  −p

1 p p 
= − cos np − cos np
p n n 
 2
 n n = 1, 3, 5, ...
2
= − cos np = 
n − 2 n = 2, 4, 6, ...
 n


∴ f ( x ) = ∑ bn sen nx
n =1

 1 1 1 
= 2 sen x − (sen 2 x ) + (sen 3 x ) − (sen 4 x ) + ... .
 2 3 4 

 1 1 1 
= 2  sen x − sen 2 x + sen 3 x − sen 4 x + ... .
 2 3 4 

Sean: s1 = 2sen x
s2 = 2sen x − sen 2 x

2
s3 = 2sen x − sen 2 x + sen 3 x, entonces,
3
Fórmulas de Euler 449

0
p 0.75p 0.5p 0.25p 0 0.25p 0.5p 0.75p

1

2

3

Figura 8-7.

x −p 3p p p p 0 p p p 3p p
− − − −
4 2 4 6 6 4 2 4
s1 0 −1.4 −2 −1.4 −1 0 1 1.4 2 1.4 0

s2 0 −2.4 −2 −0.4 −0.134 0 0.134 0.4 2 2.4 0


s3 0 −3 −1.33 −0.88 −0.8 0 0.81 1.83 1.33 2.87 0

Esta misma función en el intervalo (1, 1), Mathematica la aproxima y


grafica como:

<<Calculus `FourierTransform`
Clear[f]
f[x_]=x;
fstable=Table[FourierTrigSeries[f[x],{x,-
1,1},n],{n,2,8,2}];
TableForm[fstable]
fstable
g[x_]:=x/;-1“x“1
g[x_]:=g[x-2]/;x>1
plotfourier[i_]:=Plot[{fstable[[i]],g[x]},
{x,-1,3},
PlotStyleq{GrayLevel[0],GrayLevel[.2]},
DisplayFunctionqIdentity];
toshow=Partition[Table[plotfourier[i],
{i,1,4}],2];
Show[GraphicsArray[toshow]]

También podría gustarte