SeriesFourierSS

Download Report

Transcript SeriesFourierSS

Series de Fourier
Sensibilización
1
Matemáticas Avanzadas para Ingeniería
Fue en Grenoble ( Sur de Francia ) donde
Jean-Baptiste Joseph Fourier (1768–1830)
condujo sus experimentos sobre la
propagación del calor que le permiten
modelar la evolución de la temperatura a
través de series trigonométricas.
Estos trabajos provocan un adelanto en el
proceso de modelación matemática en
fenómenos físicos y contribuyeron a los
fundamentos de la termodinámica.
Sin embargo, la simplificación excesiva que
proponen estas herramientas fue muy
debatida, principalmente por Pierre-Simon
Laplace y Joseph-Louis Lagrange.
Departamento de Matemáticas
Campus MONTERREY
2
Matemáticas Avanzadas para Ingeniería
Fourier, seguidor de la teoría matemática de la conducción del
calor, estableció la ecuación diferencial parcial que gobierna la
difusión del calor solucionándolo por el uso de series infinitas
de funciones trigonométricas. En esto introduce la
representación de una función como una serie de senos y
cosenos, ahora conocidas como las series de Fourier.
El trabajo de Fourier provee el ímpetu para más tarde trabajar en
series trigonométricas y la teoría de las funciones de variables
reales.
En la obra Théorie analytique de la chaleur (Teoría Analítica del
calor) (1822) de Fourier, los dos primeros capítulos tratan
problemas sobre difusión de calor entre cuerpos disjuntos en
cantidad finita, es decir el problema discreto.
Departamento de Matemáticas
Campus MONTERREY
3
Matemáticas Avanzadas para Ingeniería
En el capítulo III Difusión del calor en un cuerpo
rectangular infinito es donde Fourier introduce su
método original de trabajo con series trigonométricas.
Otro trabajo importante de J. Baptiste J. Fourier fue en el
método de eliminación para la solución de un sistema de
desigualdades, teoría muy usada actualmente para
programación lineal.
Departamento de Matemáticas
Campus MONTERREY
4
Matemáticas Avanzadas para Ingeniería
Conceptos Principales
En 1807, Fourier, establece en los trabajos presentados en el instituto de
Francia que: cualquier señal periódica puede ser representada por una
serie de sumas trigonométricas en senos y cosenos relacionadas
armónicamente.
Los argumentos establecidos por Fourier eran imprecisos y en 1829
Dirichlet proporcionó las condiciones precisas para que una señal
periódica pueda ser representada por una serie de Fourier.
Fourier obtuvo además, una representación para señales no periódicas,
no como suma de senoides relacionadas armónicamente, sino como
integrales de senoides, las cuales no todas están relacionadas
armónicamente.
Al igual que las series de Fourier, la integral de Fourier, llamada
Transformada de Fourier, es una de las herramientas más poderosas
para el análisis de sistemas LTI (Sistema Lineal Invariante en el Tiempo).
Departamento de Matemáticas
Campus MONTERREY
5
Matemáticas Avanzadas para Ingeniería
La mayoría de las señales se distorsionan cuando
pasan a través de un dispositivo lineal e invariante
en el tiempo, y la única señal que no sufre
distorsión es una señal sinusoidal pura.
Un ejemplo de una señal periódica
y sus componentes de frecuencia.
Las primeras componentes de
frecuencia son:
Sumando las primeras 40
Sumando las primeras 3
componentes de
componentes de
frecuencia de la señal periódica.
frecuencia de la señal periódica.
Departamento de Matemáticas
Campus MONTERREY
6
Matemáticas Avanzadas para Ingeniería
Aplicaciones
Contenido: Señales y sistemas .
Señal
Sistema
Respuesta
Señal= f (t ) función real
del tiempo
Señal eléctrica: forma de onda de
voltaje o corriente
Resistencia estándar, para todos los
cálculos de energía y potencia se
asume una resistencia de 1 ohm.
f (t )  A cos( t   )
A : am plitud
 : frecuenciaangular 2f
 : fase, f : frecuencia
Departamento de Matemáticas
Campus MONTERREY
7
Matemáticas Avanzadas para Ingeniería
Justificación: :Señales y sistemas .
Análisis de Fourier descompone una señal en una suma de señales
senoidales y analiza como se distribuye la energía y la potencia en
cada una de esas frecuencias
Las señales se clasifican en:
• Señales de Energía y
• Señales de Potencia
Las Señales de Energía es una señal en forma de pulso que existe
sólo en un intervalo finito de tiempo, o en la que al menos tiene la
mayor parte de la energía concentrada en un intervalo finito de
tiempo
Departamento de Matemáticas
Campus MONTERREY
8
Matemáticas Avanzadas para Ingeniería
Señales y sistemas .
Energía disipada por la señal en el intervalo de tiempo es:
2
t2
E
f (t ) dt
t1

Señal de energía se define como la señal que tiene energía
finita aún cuando el intervalo de tiempo es infinito esto es
cuando
2

E
f (t ) dt ; menorque 


Departamento de Matemáticas
Campus MONTERREY
9
Matemáticas Avanzadas para Ingeniería
Señales y sistemas .
Ejemplos:
Señales de pulso Rectangular
Señales de pulso Senoidal
Señales de pulso Exponencial
Señales de pulso Gaussiano
Un ejemplo de una señal periódica
y sus componentes de frecuencia.
Las primeras componentes de
frecuencia son:
Sumando las primeras 3 componentes de
frecuencia de la señal periódica.
Departamento de Matemáticas
Campus MONTERREY
10
Matemáticas Avanzadas para Ingeniería
Señales y sistemas .
Potencia promedio disipada por la señal
2
t2
1
P
f (t ) dt
(t2  t1) t1

Departamento de Matemáticas
Campus MONTERREY
11
Matemáticas Avanzadas para Ingeniería
Señales y sistemas .
Señal de Potencia, se define como la señal que tiene potencia
promedio finita, pero diferente de cero, aún cuando el
intervalo de tiempo es infinito esto es cuando:
0  lim T  
T /2
T / 2
2
f (t ) dt  
Departamento de Matemáticas
Campus MONTERREY
12
Matemáticas Avanzadas para Ingeniería
Cuando se aplica a un sistema lineal invariante en el
tiempo, una señal sinusoidal no cambia su forma
pero sí cambian:
– Su amplitud.
– Su fase.
• En general, el cambio en la amplitud y en la
fase dependen:
– del sistema.
– de la frecuencia de la señal sinusoidal.
Departamento de Matemáticas
Campus MONTERREY
13
Matemáticas Avanzadas para Ingeniería
Para entender las causas que originan esta distorsión es
necesario analizar el contenido de frecuencias de las
señales utilizadas en ingeniería, el análisis de Fourier
permite conocer el contenido de frecuencias de las señales y
entender las razones para las cuales existe distorsión lineal.
Departamento de Matemáticas
Campus MONTERREY
14
Matemáticas Avanzadas para Ingeniería
Sensibilización:
Otra razón para estudiar a Fourier
Las vibraciones en una membrana o un tambor o las
oscilaciones inducidas en una cuerda de guitarra o violín
son explicadas por una ecuación diferencial parcial
llamada ecuación de onda .
Esta situaciones junto con condiciones iniciales y de
frontera constituyen información para encontrar la
solución única de la ecuación parcial. Pues bien la
solución es una suma infinita de funciones seno, una
forma de expresión de series de Fourier.
Las imágenes obtenidas por los investigadores revelan que las
membranas de los glóbulos rojos pierden flexibilidad, lo cual acaba
conduciendo a la aglomeración de las células, cuando éstas tratan de
navegar por los diminutos vasos sanguíneos. Asimismo, se evidencia la
destrucción de la hemoglobina, la molécula fundamental que los
glóbulos rojos usan para el transporte de oxígeno
Imágenes en 3D de un glóbulo rojo
invadido por el parásito de la
malaria. (Foto: YongKeun Park,
Michael Feld y Subra Suresh)
http://www.falstad.com/m
embrane/
Departamento de Matemáticas
Campus MONTERREY
15
Matemáticas Avanzadas para Ingeniería
Hacia las Series de Fourier ( justifications matemáticas)
La teoría de Series de Fourier trabaja con desarrollos en series trigonométricas.
Primero revisaremos algunas propiedades de las funciones, particularmente importantes
para este estudio: la continuidad por partes, la periodicidad y la simetría par e impar.
Continua por partes
Un función es continua por partes en [a, b] si f es continua en cada punto [a, b], excepto
posiblemente para un número finito de puntos donde f tiene una discontinuidad de salto . Tales
funciones son integrables en cualquier intervalo finito donde sean continuas por partes.
Periodicidad
Una función es periódica con periodo T si f ( x  T )  f ( x) para toda x en el dominio de f . Si se
cumple lo anterior, tambien se cumple f(t)=f(t+2*T)=f(t+3T) etc. El menor valor positivo se llama el
período fundamental.
Las funciones trigonométricas sen x y cos x son ejemplos de funciónes periódicas, con período
fundamental 2π y tan x es períodica con período fundamental π.
Departamento de Matemáticas
Campus MONTERREY
16
Matemáticas Avanzadas para Ingeniería
Función Par
Una función par es aquella que satisface f ( x)  f ( x)
. Tiene una gráfica que es simétrica con respecto al eje y.
para toda x en el dominio de f
2
4
1
,
x
,
x
,....
Ejemplos
cos x
Función Impar
Una función impar es aquélla que satisface f ( x)   f ( x) para toda x en el dominio de
f . Tiene una gráfica que es simétrica con respecto al origen.
Ejemplos
x, x3 , x5
Departamento de Matemáticas
Campus MONTERREY
17
Matemáticas Avanzadas para Ingeniería
Propiedades
El producto de dos funciones pares es una función par.
El producto de dos funciones impares es una función par.
El producto de una función par y una impar es impar
La suma ( resta ) de dos funciones pares es una función par.
La suma ( resta ) de dos funciones impares es una función impar.
a
a
a
0
Si f es una función par, entonces
Si f es una función impar entonces
 f ( x)dx  2  f ( x)dx
a
 f ( x)dx  0
a
Departamento de Matemáticas
Campus MONTERREY
18
Determina si la funciones siguientes son de la forma par o
impar, o ninguna de ellas.
2) f ( x)  x cos x
4) f ( x)  x  4 x
3
6) f ( x)  e  e
x
x
x  5  2  x  0
8) f ( x)  
 x  5 0  x  2
Departamento de Matemáticas
Campus MONTERREY
19
Matemáticas Avanzadas para Ingeniería
Producto interno de Funciones
El producto interno de dos funciones f1, f 2
es el numero obtenido al evaluar la integral
en un intervalo [a,b]
b
f1, f 2 
 f1( x) f 2 ( x)dx
a
Funciones ortogonales
Se dice que dos funciones f1, f 2 son ortogonales en un intervalo [a,b] si el producto
interno entre ellas es cero, es decir si:
b
f1 , f 2   f1 ( x) f 2 ( x) dx  0
a
Departamento de Matemáticas
Campus MONTERREY
20
Matemáticas Avanzadas para Ingeniería
Determina si las funciones dadas son
ortogonales en el intervalo indicado.
2) f1 ( x)  x 3 ,
4) f1 ( x)  cos x,
7) f1 ( x)  senx,
f 2 ( x)  x 2  1
en [1,1]
f 2 ( x)  sen3 x en [ ,  ]

f 2 ( x)  sen3 x en [0, ]
2
nota : el producto entre sen(m * x) con m impar cumple esto. tambien.
x
8) f1 ( x)  1, f 2 ( x)  cos
en [2,2]
2
Departamento de Matemáticas
Campus MONTERREY
21
Matemáticas Avanzadas para Ingeniería
8) f1 ( x)  1,
9) f1 ( x)  1,
f 2 ( x)  cos
x
f 2 ( x)  cos
2
x
en [2,2]
en [3,3]
3
x
10) f1 ( x)  1, f 2 ( x)  sen
en [2,2]
2
x
11) f1 ( x)  1, f 2 ( x)  sen
en [3,3]
3
x
* f1 ( x)  1, f 2 ( x)  sen
en [ p, p]
p
Departamento de Matemáticas
Campus MONTERREY
22
Matemáticas Avanzadas para Ingeniería
EN GENERAL, EL CONJUNTO
FORMA UN CONJUNTO ORTOGONAL, ES DECIR…
LOS PRODUCTOS INTERNOS ENTRE ELLOS SON SIEMPRE
CERO EN EL INTERVALO [-P,P]
Departamento de Matemáticas
Campus MONTERREY
23
Matemáticas Avanzadas para Ingeniería
A continuación algunos lineamientos:
•Norma cuadrada
 f (t )
2
de la función f (t )
Departamento de Matemáticas
Campus MONTERREY
24
Matemáticas Avanzadas para Ingeniería
f (t ) 


f n  n (t )
1
Dos funciones complejas
intervalo [t1, t2] si
* (t )
n (t ) m
t2
t2
t1
t1
* (t )d t 

(
t
)

m
 n
son ortogonales en el
* (t ) (t ) d t  0

n
 m
Dos funciones complejas son mutuamente ortogonales en el
intervalo [t1, t2] si
t2
 n (t ) m (t )d t
*
 0
si
n  m;
t1
 K n si
Departamento de Matemáticas
Campus MONTERREY
n  m
25
Matemáticas Avanzadas para Ingeniería
• Se dice que el conjunto de funciones
normalizado si
Kn 
t2

n (t )
está
n (t ) 2 d t  1, p a ra to d a n
t1
• Si el conjunto de funciones n (t ) es ortogonal y esta
normalizado se llama conjunto ortonormal
• Una función arbitraria f (t ) se puede representar en
una serie de funciones ortogonales como
f (t ) 


n 1
fn
• en donde los coeficientes
como
 n (t )
f n (t )
pueden determinarse
Departamento de Matemáticas
Campus MONTERREY
26
Matemáticas Avanzadas para Ingeniería
• en donde los coeficientes
como sigue: Sea:
f (t ) 
f n (t )


n 1
* 
f (t )m


fn

Por tanto:
fn
 n (t )
*
n (t )m
n 1
t2
t2 
*
f (t )m 
fn
t1
t1 n 1

pueden determinarse
* dt 
n (t )m
1
fn 
Kn
t2


fn 

n 1
t2
fn
 n
* dt  
m
mKm
t1
* (t )dt
f (t )n
t1
t2
O bien:

t1
t2
*d t 
f (t )n
 n m d t 
*
t1
Departamento de Matemáticas
Campus MONTERREY
27
Matemáticas Avanzadas para Ingeniería
• Cálculo del error que se tiene al aproximar f (t ) con
una sumatoria de N términos en lugar de una serie infinita
N
 f n n (t )  N (t )
f (t ) 
n 1
De donde el error cuadrático :
t2 
N
N


 N (t ) dt   f (t ) 
f nn (t )  f (t ) 
f nn (t )



n

1
n

1




t1
t1
t2

2
t2
O bien:

t1


2
 N (t ) dt 
t2

t1
2
f (t ) dt 

N
 fn
n 1
2
Kn
Departamento de Matemáticas
Campus MONTERREY
28
Matemáticas Avanzadas para Ingeniería
• desarrollando se llega a
2
t2

f (t )
2

dt 

f (t )
n 1
t1
1
t

n
K n,
Teorem ade Parseval

1

De donde el error cuadrático para un conjunto ortogonal
completo :
f
(t
)
t 2
d
t
2


f
(t
)
K
n,
Te o
re m
ad
e
P
a
r
se v a
l
2
t2
LimN 
  N (t )
2
dt  0
t1
f (t ) 

 f nn (t )
n 1
Es la representación en una serie
generalizada de Series de Fourier
Departamento de Matemáticas
Campus MONTERREY
29
Matemáticas Avanzadas para Ingeniería
• Cierta función rectangular esta definida
 1
f (t )  
 1,
0  x 1
1 x  2
Se desea aproximar esta función de energía finita empleando un
conjunto de funciones definidas por
n (t )  sen(n t ), n m ayor que 0
Solución: El conjunto sen(n t )es un conjunto ortonormal
en
2
1
sen(n  t ) sen(m t )  
0
0

nm
nm
Departamento de Matemáticas
Campus MONTERREY
30
Matemáticas Avanzadas para Ingeniería
• Por lo tanto
f (t ) 

 f n sen(n  t )
n 1
2
En donde
 f (t ) sen(n  t )dt
fn  0
2
sen2 ( n  t ) dt
0

2
1
2
f n  f (t )sen (n t )dt  sen (n  t )dt  sen (n  t )dt
0
0
1



 cos(n  t ) 1 cos(n  t ) 2 2


 {1  cos (n  )}
n
n
0
1 n
Departamento de Matemáticas
Campus MONTERREY
31
Matemáticas Avanzadas para Ingeniería
• Por lo tanto
1
1


f (t )  4  sent  sen3  t  sen(5  t )  ......
3
5


En donde el error cuadrático integral puede calcularse a partir de
2
t2
 N   N dt 
t1

t2
2
f (t ) dt 
t1

22
4
e  2  ( )  .379
01


2
2 2
4  4
e  2  
0 3
    3

N
 f (t )
n 1
2
Kn
2
2 2
4  4
e  2  
0 5
    3

2

  .119

2
  4
 
  5
2

  .134

2
2
2
2
22
4  4   4   4 
e  2              .101
0 7
    3   5   7 

Departamento de Matemáticas
Campus MONTERREY
32
Matemáticas Avanzadas para Ingeniería
• Por lo tanto
22
4
e  2  ( )  .379/ 2  19%
01


2
2
22
4  4 
e  2        .119 / 2  10%
0 3
    3 

2
2
2
22
4  4   4 
e  2           .134 / 2  6.7%
0 5
    3   5 

2
2
2
2
22
4  4   4   4 
e7  2              .101 / 2  5.1%
0
    3   5   7 

Alrededor del 95% de la energía esta contenida en los primeros
cuatro términos
Departamento de Matemáticas
Campus MONTERREY
33
Matemáticas Avanzadas para Ingeniería
 SERIE DE FOURIER ( Equivalencia en la representación en serie trigonométrica
de f(t) )
Definición 1. Sea f una función continua por partes en el intervalo [-T,T].
La serie de Fourier de f es la serie trigonométrica

a0
nx
nx
f ( x) 
  (an cos
 bn sen
)
2 n 1
p
p
Donde
y
an
an 
están dadas por las fórmulas:
bn
1 p
nx
f
(
x
)
cos
dx


p
p
p
1 p
nx
bn   f ( x)sen
dx

p
p
p
n  0,1,2,....
n  1,2,3,....
Se prueba
integrando ambos
lados de manera
conveniente
Departamento de Matemáticas
Campus MONTERREY
34
Matemáticas Avanzadas para Ingeniería
EJEMPLO 1 Calcular la serie de Fourier de
  x  0
0 x 
0,
f ( x)  
 x,
Solución
En este caso, T=π. Obtenemos los coeficientes con las fórmulas anteriores.
an 
1



f ( x) cos nxdx 

n
1
n 2

1
1


cos
n


1

n 2
n 2
a0 
1


0



f ( x)d x 
1



 x cos nxdx
0
1
n


cos
u

usenu
0
n 2

u cos udu 
1


 1
n
xd x 
0
Departamento de Matemáticas
Campus MONTERREY

1 ,
2
x
2
n  1,2,3,...


0

2
35
Matemáticas Avanzadas para Ingeniería
EJEMPLO 1 (continuación)
Solución
bn 

1



f ( x) sen(nx) dx 

1
n 2

n
0
usen (u ) du 
 cos n
(1) n 1


.
n
n
1


 xsen(nx)dx
0
1
n


sen
(
u
)

u
cos
u
0
n 2
n  1,2,3,..
Por lo tanto,

 1

( 1) n 1
n
f ( x) 
 
(

1
)

1
cos
nx

sen
(
nx
)

2
4
n
n 1  n


2 
1
1



cos5 x  ...
cos x  cos3x 
4
 
9
25




1
1


 senx 
sen2 x  sen3x  ...
2
3


Departamento de Matemáticas
Campus MONTERREY
36
Matemáticas Avanzadas para Ingeniería
EJEMPLO 2 Calcular la serie de Fourier de
 0,
f ( x)  
  x,
  x  0
0 x 
Solución: Usted podrá llegar a esta representación, puede
utilizar el tutorial, si lo hace en la calculadora, puede
acceder al proceso de solución en la pagina del curso.
f ( x) 




1
 1

n
 
1

(

1
)
cos
nx

sen
(
nx
)

2
4
n

n 1  n
Departamento de Matemáticas
Campus MONTERREY
37
Matemáticas Avanzadas para Ingeniería
Gráfica del problema anterior
http://cb.mty.itesm.mx/ma3002/home.htm
Departamento de Matemáticas
Campus MONTERREY
38
Matemáticas Avanzadas para Ingeniería
EJEMPLO 3 Calcular la serie de Fourier de
0,
f ( x)   2
x ,
  x  0
0 x 
Departamento de Matemáticas
Campus MONTERREY
39
Matemáticas Avanzadas para Ingeniería
EJEMPLO 4 Calcular la serie de Fourier de
  x  0
0 x 
 1,
f ( x)  
 1,
De nuevo, T=π. Observe que f es una función impar. Como el producto de una
función impar y una función par es impar, f(x) cos nx también es una función impar.
Así,
1 
an   f ( x) cos nxdx  0,
n  0,1,2,....
 
Solución
Además, f(x) sen nx es el producto de dos funciones impar y por tanto es una función par, de modo que
bn 
1




f ( x) sen(nx)dx 
2



0

2   cos nx 
2  1 (1) n 
sen(nx)dx  
    n  n ,
 
n


0
n  1,2,3,...
 0,

 4
,

n
Así
[1  (1) n ]
4
f ( x) ~ 
sen(nx) 
 n1
n

2

n
par
n
im par
1
1


sen
(
x
)

sen
(
3
x
)

sen
(
5
x
)

...


3
5
Departamento de Matemáticas
Campus MONTERREY
40
Matemáticas Avanzadas para Ingeniería
EJEMPLO 5 Calcular la serie de Fourier de
Solución
f ( x)  x ,
1  x  1
En este caso T=1. Como f es una función par, f(x)sen nπx es una función impar.
Por consiguiente,
1
bn   f ( x) sen(nx)dx  0,
n  1,2,3,....
1
Como f(x) cos nπx es una par, tenemos
1
1
1
1
0
0
a0   f ( x)dx  2 xdx  x 2
1
n
2
u cos udu
0
1
0
 2 n 2 0
2
2
2
n
 2 2 [cos u  usenu ] 0  2 2 (cos n  1)  2 2 [( 1) n  1],
 n
 n
 n
1
1
an   f ( x) cos nxdx  2 x cos nxdx  x 2
1

n  1,2,3,...
Por lo tanto
1  2
1
4
f ( x) ~   2 2 [(1) n  1] cos nx   2
2 n 1  n
2 
1
1


cos5x  ...
cosx  cos3x 
9
25


Departamento de Matemáticas
Campus MONTERREY
41
Matemáticas Avanzadas para Ingeniería
 RECORDEMOS LAS PROPIEDADES DE FUNCIONES PARES E IMPARES
SUPONGAMOS QUE f (x) espar

a0
nx
nx
f ( x) 
  (an cos
 bn sen
)
2
p
p
n 1
ENTONCES
YA QUE EL PRODUCTO DE PARES ES PAR
1 p
nx
an   f ( x) cos
dx
p p
p
Se
prueba
integran
do
ambos
lados de
manera
convenie
nte
2 p
nx
  f ( x) cos
dx
p 0
p
YA QUE EL PRODUCTO
DE PAR POR IMPAR ES
IMPAR
1 p
nx
bn   f ( x) sen
dx  0
p p
p
Departamento de Matemáticas
Campus MONTERREY
42
Matemáticas Avanzadas para Ingeniería
 SERIES DE SENOS Y COSENOS DE FOURIER
Definición 2. Sea f(x) continua por partes en el intervalo [0,T]. La serie de cosenos de
Fourier de f(x) en [0,T] es

a0
nx
  an cos
,
2 n 1
T
Donde
an 
2
T

T
0
f ( x) cos
nx
dx ,
T
n  0,1,...
f ( x) sen
nx
dx ,
T
n  1,2,...
La serie de Fourier de senos de Fourier de f(x) en [0,T] es

 b sen
n 1
n
nx
,
T
Donde
bn 
2
T

T
0
Departamento de Matemáticas
Campus MONTERREY
43
Matemáticas Avanzadas para Ingeniería
EJEMPLO 1 Calcular la serie de Fourier de
Solución
bn 
2


2
2
 x 
Usamos las fórmulas anteriores con T=π, para obtener


0


0 x 

 x,
f ( x)  
  x,


f ( x) sen( nx) dx 
2
n 2

n / 2
0
2

 /2

2
xsen( nx)dx 

0
usenudu  2 

 /2
sen ( nx ) dx 


2
n 2
(  x) sen( nx)dx
/2
n

n/2
usen (u ) du
2
2
n 
2
n / 2
n


sen
(
u
)

u
cos
u

cos

n

cos

[
sen
(
u
)

u
cos(
u
)]
2
0
n / 2
n 2
n
2 

 n
0,

4
n


sen
  4( 1) ( n 1) / 2
2
,
n
2

n 2

Departamento de Matemáticas
Campus MONTERREY
n
n
par
im p a r
44
Matemáticas Avanzadas para Ingeniería
EJEMPLO 1 (continuación)
Solución
Así que al hacer n=2k+1, tenemos que la serie de senos de Fourier para f(x) es
(1) k
4 
1
1


sen
(
2
k

1
)

senx

sen
3
x

sen
5
x

...



 k 0 (2k  1) 2
 
9
25

4

La función f(x) es continua y f ´(x) es continua por partes en (0,π), de modo que el teorema de la
convergencia puntual de las series de Fourier implica que
f ( x) 
4 
1
1

sen3x 
sen5 x  ....
senx 
 
9
25

Para toda x en [0,π]
Departamento de Matemáticas
Campus MONTERREY
45
Matemáticas Avanzadas para Ingeniería
SERIE DE FOURIER COMPLEJA
Sea f una función de variable real, periódica con periodo fundamental p. Supongamos que f es
integrable en [-p/2,p/2]. La serie de Fourier en este intervalo es

1
a0   an cos(no x)  bn sen(no x),
2
n 1
Con 0  2 / p . Se reescriben las ecuaciones como





1
1 in0 x
1


a0   an
e
 e in0 x  bn
e in0 x  e in0 x ,
2
2
2i

n 1 

1
1
1


a0    ( an  ibn )ein0 x 
(an  ibn )e in0 x ,
2
2

n 1  2
En la serie sea
d0 
1
a0
2
y para cada entero positivo
dn 
Departamento de Matemáticas
Campus MONTERREY
1
( a n  ibn )
2
46
Matemáticas Avanzadas para Ingeniería
SERIE DE FOURIER COMPLEJA (continuación)
Entonces la serie llega a ser

d 0   [d n e
in0 x
n 1
 dne
in0 x
Ahora consideramos los coeficientes. Primero,

]  d0   dne

  d n e in0 x
in0 x
n 1
1
1
a0 
2
p
d0 
n 1

p/2
p/2
f (t )dt
Y, para n=1,2,…
dn 
1
1 2
(an  ibn ) 
2
2 p


1
p
p/2
f (t ) cos(n0t )dt 
p/2

p/2
p/2
i 2
2 p

p/2
p/2
f (t ) sen(n0t )dt
f (t )e in0t dt
Y, para n=1,2,…
dn 
1
p

p/2
p/2
f (t )e
in 0 t

1
p

p/2
p/2
Departamento de Matemáticas
Campus MONTERREY
f (t )ein0t dt  d  n
47
Matemáticas Avanzadas para Ingeniería
SERIE DE FOURIER COMPLEJA (continuación 2)
Ponemos estos resultados en la serie para obtener

d0   d ne
in0 x
n 1

 d0   dne
in0 x
n
n   , n  0

  d  n e in0 x
n 1

d
  d n e in0 x
n 1
n 1
 d0 

e
in0 x


d
n  
n
ein0 x
Hemos encontrado esta expresión rearreglando los términos en la serie de Fourier de una función periódica f
.
Departamento de Matemáticas
Campus MONTERREY
48
Matemáticas Avanzadas para Ingeniería
 TEOREMA
Sea f periódica con un periodo fundamental p. Sea f suave a pedazos en [-p/2, p/2].
Entonces, en cada x la serie de Fourier converge a
1
( f ( x )  f ( x ))
2
El espectro de amplitud de la serie de Fourier compleja de una función periódica es la gráfica de
los puntos (n0 , dn ), en donde d n es la magnitud del complejo d n .
Algunas veces este espectro de amplitud es llamado también espectro de frecuencia.
Departamento de Matemáticas
Campus MONTERREY
49
Matemáticas Avanzadas para Ingeniería
EJEMPLO
Sea f ( x)  3x / 4 para
periodo fundamental 8.
0  x  8, y
f ( x  8)  f ( x)
para todo x. Entonces f es periódica con
Aquí p=8 y 0   / 4 . Recordemos que en las fórmulas para los coeficientes se puede realizar la
integración sobre cualquier intervalo de longitud 8. Aquí es conveniente usar [0,8] en lugar de [-4,4]
debido a como está definida f(x). Entonces.
1 0 3
d0 
Si usamos el intervalo [-4,4], entonces podríamos calcular
Ahora,
dn 
1 8 3  ni / 4
3i
te
dt

8 0 4
n
8  4 4
tdt  3
d0 
1 0 3
1 43
(
t

8
)
dt

tdt  3
8  4 4
8 0 4
La serie de Fourier
compleja es
3i
3
n

1 nix / 4
e

n
n   , n  0
Esta serie converge a f(x) para 0< x <8, 8< x <16, 16 < x < 24…. .8< x <0, -16 < x < -8
3
Para trazar el espectro de amplitud, calculamos d 0  3, d n 
. Como n0  n / 4 el espectro de
n
amplitud es un trazo de los puntos 

 n , 3 
 4 n


Departamento de Matemáticas
Campus MONTERREY
50