Introduccin

La solucin de funciones peridicas por medio de series de Fourier se facilita con el uso del concepto de Funciones Ortogonales. Siendo su uso importante para la simplificacin de funciones, como tambin en el uso de sucesiones como las de Cauchy y Hilbert .

Funciones Ortogonales
En matemticas superiores se considera que una funcion es la generalizacion de un vector. En esta seccion veremos como los dos conceptos vectoriales de producto interno (punto) y ortogonalidad se pueden ampliar para abarcar las funciones. Supongamos que u y v son vectores en el espacio tridimensional. El producto interno (u, v) de los vectores, posee las propiedades siguientes: i) (u, v) = (v, u) ii) (ku, v) = k(u, v), donde k es un escalar iii) (u, u)= 0, si u= 0,y (u,u)>0 si u 0 iv) (u + v, w) = (u, w) + (v, w).
2

Producto interno
El producto interno de dos funciones 1 y 2 en un intervalo [a, b] es el numero

Funciones ortogonales Dado que dos vectores u y v son ortogonales cuando su producto interno es cero, definiremos las funciones ortogonales en forma semejante.

Conjuntos ortogonales
Un conjunto de funciones de valor real {F0 (x), F1 (x), F2 (x),...} es ortogonal en un intervalo [a, b] si

Se llama norma cuadrada de n, Si {n(x)} es un conjunto ortogonal de funciones en el intervalo [a, b] y tiene la propiedad que Il n (x)Il = 1 para n = 0, 1, 2, ..., se dice que { n(x)} es un conjunto ortonormal en el intervalo.

Funciones ortogonales
Dos funciones 1 y 2 son ortogonales en un intervalo [a, b] si

A diferencia del anlisis vectorial, en donde la palabra ortogonal es sinnimo de "perpendicular", en el presente contexto el termino ortogonal y la condicin (1) no tienen significado geomtrico.

EJEMPLO 1 Funciones ortogonales Las funciones 1 (x) = x2 y 2 (x)= x3 son ortogonales en el intervalo [-1, 1] porque

Funciones Ortogonales
Se dice que dos funciones reales gm(t) y gn(t) son a t b ortogonales en el intervalo Si la integral del producto gm(t) y gn(t) sobre este intervalo a t b es igual a cero, es decir :
b

gm(t) y gn(t) dt
a

0, m

n
6

Observaciones
Un conjunto de funciones reales g 0, g1, g2,, que satisface la ecuacin 1 para todos lo pares de funciones distintos en le conjunto es un conjunto ortogonal de funciones en ese intervalo.
b

La integral N ( g m (t )) de la funcin gn(t), en el intervalo

2 gm (t )dt se llama Norma a

a t

b
7

Observaciones
Un conjunto ortogonal de funciones, g 0, g1, g2,, en el intervalo a t b tal que cuyas funciones tiene norma 1 satisface las relaciones
b

gm(t) y gn(t) dt
a

0, m

n,n

0, 1, 2, ...

gm(t) y gn(t) dt 1, m
a

n,n

0, 1, 2, ...
8

Aplicaciones
Fourier no solo dedujo la ecuacin de difusin del calor sino que estableci un programa sistemtico de resolucin de ecuaciones en derivadas parciales y el mtodo de separacin de variables utilizado en el desarrollo de circuitos elctricos.

Aplicaciones
El poder extraordinario y la flexibilidad de las series y transformadas de Fourier se ponen de manifest en la asombrosa variedad de las aplicaciones que ellas tienen en diversas ramas de la matemtica y de la fsica matemtica, desde teora de nmeros y geometra hasta mecnica quntica. Para mas informacin sobre las aplicaciones visitar www.emis.de/journals/DM/v5/art6.pdf

10

Introduccin
El anlisis de Fourier fue introducido en 1822 en la Thorie analyitique de la chaleur para tratar la solucin de problemas de valores en la frontera en la conduccin del calor.
Ms de siglo y medio despus las aplicaciones de esta teora son muy bastas: Sistemas Lineales, Comunicaciones, Fsica moderna, Electrnica, ptica y por supuesto, Redes Elctricas entre muchas otras.

Funciones Peridicas
Una Funcin Peridica f(t) cumple la siguiente propiedad para todo valor de t. f(t)=f(t+T) A la constante mnima para la cual se cumple lo anterior se le llama el periodo de la funcin Repitiendo la propiedad se puede obtener: f(t)=f(t+nT), donde n=0, 1, 2, 3,...

Ejemplo: Cul es el perodo de la funcin

Solucin.- Si f(t) es peridica se debe cumplir: Pero como se sabe cos(x+2k )=cos(x) para cualquier entero k, entonces para que se cumpla la igualdad se requiere que T/3=2k1 , T/4=2k2 Es decir, T = 6k1 = 8k2 Donde k1 y k2 son enteros, El valor mnimo de T se obtiene con k1=4, k2=3, es decir,T=24

Grfica de la funcin
3 2 1

f(t)

t t cos( 3 ) cos( 4 )

f(t)=cos(t/3)+cos(t/4)

f(t)

0 -1
-2

24
-3
0 50 100 150 200

Podramos pensar que cualquier suma de funciones seno y coseno produce una funcin peridica.

Esto no es as, por ejemplo, consideremos la funcin f(t) = cos( 1t)+cos( 2t). Para que sea peridica se requiere encontrar dos enteros m, n tales que 1T= 2 m, 2T=2 n De donde Es decir, la relacin racional.
1/
1 2

debe ser un nmero

m n

Ejemplo: la funcin cos(3t)+cos( +3)t no es peridica, 3 no es un nmero racional. ya que 1

2
2

3
f(t)=cos(3t)+cos((3+pi)t)

f(t)
0 -1 -2
0

10

15

20

25

30

Serie Trigonomtrica de Fourier

Algunas funciones peridicas f(t) de periodo T pueden expresarse por la siguiente serie, llamada Serie Trigonomtrica de Fourier f(t) = a0 + a1cos( 0t)+a2cos(2 0t)+... + b1sen( 0t)+b2sen(2 0t)+... Donde 0=2 /T. Es decir,

f (t)

1 a 2 0

n 1

[a n cos(n

0t)

bnsen (n

0 t )]

Es posible escribir de una manera ligeramente diferente la Serie de Fourier, si observamos que el trmino ancos(n 0t)+bnsen(n 0t) se puede escribir como
a2 n b2 n an a2 n b2 n cos(n
0t)

bn a2 n b2 n

sen (n

0t)

Podemos encontrar una manera ms compacta para expresar estos coeficientes pensando en un tringulo rectngulo:

an
b
n

Cn
n

2 n

2 n

a2 n bn a2 n

b2 n b2 n

cos sen

a
n

C n cos

cos(n

t ) sen n sen (n
0

t)

C n cos(n

Si adems definimos C0=a0/2, la serie de Fourier se puede escribir como

f (t) C0
n 1

C n cos(n

As, y

Cn
n

2 n

2 n

bn tan an
1

Ortogonalidad de Senos y Cosenos

Se dice que un conjunto de funciones fk(t) son ortogonales en el intervalo a<t<b si dos funciones cualesquiera fm(t), fn(t) de dicho conjunto cumplen
b

f m (t)f n (t)dt
a

0 rn

para m para m

n n

Ejemplo: las funciones t y t2 son ortogonales en el intervalo 1< t <1, ya que

1 1

tt dt

1 1

t dt

t 4

4 1

0
1

Ejemplo: Las funciones sen t y cos t son ortogonales en el intervalo /2< t < /2, ya que

sentcostdt

sen t 2

Aunque los ejemplos anteriores se limitaron a un par de funciones, el siguiente es un conjunto de una infinidad de funciones ortogonales en el intervalo -T/2<t< T/2. 1,cos 0t, cos2 0t, cos3 0t,...,sen 0t,sen2 0t,sen3 0t,... (para cualquier valor de 0=2 /T).

Para verificar lo anterior podemos probar por pares: 1.- f(t)=1 Vs. cos(m 0t):
T/2

cos(m
T/2

0 t)dt

sen (m 0 t) m 0

T/2 T/2

2sen (m 0T/2) m 0

2sen (m ) m 0

Ya que m es un entero.

2.- f(t)=1 Vs. sen(m

T/2

0t):

sen(m
T/2

0 t)dt

cos(m 0 t) m 0 1 m
0

T/2 T/2

[cos(m 0T/2) - cos(m 0T/2)] 0

0t):

3.- cos(m
T/2

0t)

Vs. cos(n

cos(m
T/2

0 t)cos(n 0 t)dt

0 para m n T / 2 para m n 0

4.- sen(m
T/2

0t)

Vs. sen(n

0t):

sen(m
T/2

0 t)sen(n 0 t)dt

0 para m n T / 2 para m n 0

5.- sen(m
T/2

0t)

Vs. cos(n

0t):

sen(m
T/2

0 t)cos(n 0 t)dt

0 para cualquier m, n

Para calcular las integrales de los casos 3, 4 y 5, son tiles las siguientes identidades trigonomtricas: cos A cos B = [cos(A+B)+cos(A-B)] sen A sen B = [-cos(A+B)+cos(A-B)] sen A cos B = [sen(A+B)+sen(A-B)] Adems: sen2 = (1-cos2 ) cos2 = (1+cos2 )

Clculo de los coeficientes de Fourier

Dada una funcin peridica f(t) cmo se obtiene su serie de Fourier?
f (t)
1 a 2 0 n 1

[a n cos(n

0t)

bnsen (n

0 t )]

Obviamente, el problema se resuelve si sabemos como calcular los coeficientes a0,a1,a2,...,b1,b2,... Esto se puede resolver considerando la ortogonalidad de las funciones seno y coseno comentada anteriormente.

Ejemplo: Encontrar la Serie de Fourier para la siguiente funcin de periodo T:

1
... -T/ 2 T ...

f(t)

T/ 2

-1

Solucin: La expresin para f(t) en T/2<t<T/2 es

f (t) 1 para T t 2 1 para 0 t 0
T 2

Coeficientes an:
2 T 0

an

2 T

T/2

f ( t ) cos(n
T/2

0 t )dt

T/2

cos(n
T/2

0 t )dt 0
0

cos(n

0 t )dt

2 T

1 n
0

sen (n

0t) T/2

1 n
0

T/2

sen (n

0t) 0

0 para n

Coeficiente a0:
2 T

a0
0

2 T

T/2

f ( t )dt
T/2

T/2

dt
T/2
0

dt
0
T/2

2 T

t
T/2

t
0

Coeficientes bn:
2 T 0

bn

2 T

T/2

f ( t )sen (n
T/2

0 t )dt

T/2

sen (n
T/2

0 t )dt 0
0

sen (n

0 t )dt
T/2

2 T

1 n
0

cos(n

0t) T/2

1 n
0

cos(n

0t) 0

1 (1 cos(n )) (cos(n ) 1) n 2 n 1 ( 1) ) para n 0 n

SERIES DE FOURIER DE COSENOS Y DE SENOS

f ( x)

a0 2

n n (an cos x bn sen x) p p 1

ANTECEDENTES
La serie de Fourier de una funcin definida en el intervalo (-p,p) es :
f ( x) a0 2 n n (an cos x bn sen x) p p 1

a0 bn 1 p
p p

1 p

p p

f ( x)dx an 1 p
p p

n f ( x) sen xdx p

f ( x) cos

n xdx p

Series de Fourier de csenos y de senos

Si f es una funcin par en (-p,p), entonces en vista de las propiedades anteriores, los coeficientes de se transforman en :

a0
an 1 p

1 p
p

p p

f ( x)dx

2 p

p 0

n f ( x) cos xdx p

n 2 p f ( x)dx f ( x) cos xdx 0 p p p par p 1 n bn f ( x) sen x dx 0 p p p

im par

En forma parecida, cuando f es impar en el intervalo (-p,p),

an

n=0,1,2,...

bn

2 p

p 0

f ( x) sen

n xdx. p

Resumiendo quedara de la siguiente manera:

a. La serie de Fourier de una funcin par en el intervalo (-p,p) es la serie de csenos

f ( x)

a0 2

n an cos x, p 1
p 0

En que :

a0

2 p 2 p
p 0

f ( x)dx

a0

n f ( x) sen dx p

b) La serie de Fourier de una funcin impar en el intervalo (-p,p) es la serie. de senos

f ( x)
n 1

bn sen

n x, p

En donde :

bn

2 p

p 0

n f ( x) sen xdx. p

PROBLEMAS FOURIER SENOS Y COSENOS

PROBLEMA 1 Considera la funcin f : *0, 2+ R dada por f (x) = 1 x/2. Un clculo elemental muestra que los coeficientes del desarrollo en serie de Fourier de2 senos de f vienen dados por bn = n . Disea una funcin de Matlab que dibuje la suma parcial de los N primeros sumandos de la serie de Fourier de senos de la funcin y muestre simultneamente la grfica de f. Solucin. Basta crear la siguiente funcin en un archivo llamado fousen.m. Obsrvese la estructura de la lnea que define la funcin f. function fousen(N) x=-2:0.05:2; sumparcial=0; b=zeros(1,N); for k=1:N b(k)=2/(k*pi); sumparcial=sumparcial+b(k)*sin(k*x*pi/2); end f=(x<0).*(-1-x/2)+(x>=0).*(1-x/2); plot(x,f,b,x,sumparcial,g),

PROBLEMA 2 Considera la funcin 2-peridica f : R R tal que f (x) = x2 si 1 <x< 1. Disea una funcin de Matlab que aproxime numricamente los N primeros coeficientes de la serie de Fourier de la funcin y que dibuje la suma parcial frente a f. Solucin. Puesto que la funcin es par, tenemos que

De esta forma, la solucin del problema viene dada por la siguiente funcin que debe ser almacenada en un archivo llamado fou.m. function fou(N) x=-1:0.05:1; a=zeros(1,N); sumparcial=1/3; for k=1:N a(k)=quadl(@fun,0,1,1e-9,[],k); sumparcial=sumparcial+a(k)*cos(k*x*pi/2); end f=x.^2; plot(x,f,b,x,sumparcial,g),shg function y=fun(t,n) y=(t.^2).*cos(n*pi*0.5*t);

PROBLEMA 1
Serie de Fourier de cosenos Sea f(x) = x en [0, 2]. Encuentre la serie de Fourier de cosenos para f(x) Se tiene, con L = 2 y f(x) = x
a0
Entonces :

1 2

2 0

xdx 1

an

2 2

2 0

x cos

n x dx 2

4 cos(n ) 1 n2 2
cos

0 par 8 inpar 2 2 n
(2n 1) n 2

f ( x)

x 1
n

8 2 ( 2 n 1 ) 0

PROBLEMA 4
Calcular la serie de Fourier de cosenos de, Solucin. En este caso, L=1 , y aplicando las frmulas de arriba tenemos:

Por lo tanto, la serie de Fourier de cosenos para en el intervalo es:

Obviamente, existe una adaptacin de los criterios de convergencia. Por ser el criterio que nos da ms informacin, enunciamos nicamente el correspondiente al Segundo Teorema de Convergencia.

PROBLEMA 5
Calcular la serie de Fourier de senos de Solucin . En este caso, L=5, y aplicando las frmulas de arriba tenemos:

Por lo tanto, la serie de Fourier de senos de f(X) est dada por:

Esta ltima puede expresarse como:

Problema 6

Determinar una solucin peridica estacionaria de:

..(1) donde F(t) es la funcin con periodo 4, tal que F(t)= 5t, para -2 < t < 2 y la Serie de Fourier
F (t ) 20
n 1

x " 10 x

Ft

( 1) n 1 n t sen n 2

..(2)

Solucin
Sabemos que :
n t xsp (t ) bn sen 2 n 1 Remplazando en (1) tenemos:
bn
n 1

n2 4

n t 10 sen 2

20
n 1

( 1) n 1 n t sen n 2

Al igualar los coeficientes de los trminos semejantes y despejar , obtenemos la solucin peridica estacionaria
xsp (t ) 80
n 1

( 1) n 1 n 40 n 2

n t sen 2

t 2 t 3t xsp (t ) 0.8452sen 24.4111sen 0.1738sen ...... 2 2 2

La magnitud tan grande del segundo termino surge del hecho de que .As el movimiento dominante de un resorte con la ecuacin diferencial (1) ser una oscilacin con frecuencia radianes/segundo, periodo 2s y amplitud cercana a 24, consistente con la grafica de que aparece en la figura sgte.

Grafica de x sp ( t )

Problema 7
Supngase que f(t) = t para 0 < t < L. Determinar la serie de csenos de Fourier y la serie de senos de Fourier para f.

Solucin
Sabemos que :
an 2 L
L

f (t ) cos
0

n t dt L
2 L

Entonces para n = 0 :
a0

2 tdt L0

2 t L 2

L
0

Entonces :
an

2 n t t cos dt L0 L

2L n2 2
4L n2 2

u cos udu
0

si n impar ; par .

2L n usenu cos u 0 n2 2

0 si n

As la serie de csenos de fourier de f es:

t L 2 4L
2

t cos L

1 3 t cos 2 L 3

1 5 t cos 2 L 5

.....

Para 0 < t < L, tenemos que: L 2 n t bn tsen dt L0 L

2L n2 2
n 0

usenudu
0

2L n2 2

u cos u senu

2L ( 1) n n
1 3 t sen 3 L

As la serie de senos de Fourier de f es: 2L t 1 2 t t sen sen L 2 L

.....

Para 0 < t < L, la serie de csenos converge a la extensin par con periodo 2L de f como muestra la fig. 1; y la serie de senos converge al a extensin impar con periodo 2L de f como muestra la fig. 2.

Fig. 1

Fig. 2

Problema 8
Determinar la serie de Fourier de la funcin con periodo 2 definida en un periodo como :

f (t )

0 si t si 2 si

t 0 t t

0; ; .

Si la grafica de f es:

Solucin:
Los valores de f() no son importantes, pues no tienen efecto sobre los valores de las integrales que proporcionan los coeficientes de Fourier. Como f(t)=0 en el intervalo (,0), cada integral desde t =a t =se puede reemplazar por una
integral de t = 0 a t = . Por tanto

a0

1
0

tdt

1 t2 2

an

1
0

t cos ntdt

1 n
2

u cos udu
0

nt , t

u n

1 n2

cos u usenu

n 0

1 n2

En consecuencia, = 0 si n es par y n2

an
A continuacin:

2 n2

si n es impar,

bn

1
0

tsenntdt

1 n2

usenudu
0

1 n senu u cos u 0 2 n

1 cos n n

As,

bn

1 n

n 1

para toda n

1.

Por lo tanto la serie de Fourier de f(t) es:

f (t ) ~

2 4
n

cos nt 2 n im par

( 1) n 1 sennt n n 1