Introducción a Series de Fourier
Introducción a Series de Fourier
1.1 Introducción 3
Bibliografía 55
2
1. Series e integrales de Fourier
1.1 Introducción
Al igual que la teoría de Taylor, la teoría de Fourier permite representar funciones en forma
de series. Wikipedia, la enciclopedia libre, menciona:
Una serie de Fourier es una serie infinita que converge puntualmente a una función pe-
riódica y continua a trozos (o por partes). Las series de Fourier constituyen la herramien-
ta matemática básica del análisis de Fourier empleado para analizar funciones periódicas a
través de la descomposición de dicha función en una suma infinita de funciones sinusoidales
mucho más simples (como combinación de senos y cosenos con frecuencias enteras).
Por su parte, como las series de Fourier únicamente se pueden aplicar sobre funciones
periódicas, la integral y transformada de Fourier generalizan la teoría dada Fourier. La trans-
formada de Fourier es aplicada sobre funciones aperiódicas.
La serie de Fourier constituye una teoría matemática sumamente útil para representar
fenómenos vibratorios y su aplicación se extiende a la acústica, óptica, procesamiento de
imágenes y señales, y compresión de datos.
Una función dada f se dice que es periódica si sus imágenes se repiten en intervalos
regulares en su dominio. El periodo de la función, simbolizado por T , viene dada por el inter-
valo entre dos réplicas sucesivas. De esta forma, decimos que una función f es periódica con
periodo T si, para todos los valores t de su dominio y para cualquier entero m se tiene que
f (t + mT ) = f (t)
Así como se muestra en la Figura 1, una función periódica puede dividirse en tiras verticales
que son réplicas una de la otra.
3
Figura 1: Representación gráfica de función periódica
Fuente: [Link]
Una de las propiedades más importantes de las funciones periódicas es que: la combi-
nación lineal de funciones periódicas con el mismo periodo sigue siendo periódicas y conserva
el periodo. Por lo tanto, si f y g en la variable t son dos funciones periódicas, con periodo T ,
entonces la función h es periódica, con periodo T , y está dada por:
.
1.2.2 Teoría de Fourier
La teoría de Fourier consiste en la representación serial trigonométrico de una función
periódica. James (2002) menciona:
La teoría de la serie de Fourier se basa en un teorema que afirma toda función periódica
que satisface ciertas condiciones puede expresarse como la suma de un número de funciones
seno de diferentes amplitudes, fases y periodos. Esto es, si una función periódica f , en la
variable t, con periodo T satisface ciertas condiciones entonces:
4
An sen(nωt + ϕn ) = (An cosϕn ) sen(nωt) + (An senϕn ) cos(nωt) = an cos(nωt) + bn sen(nωt)
a0 ∑ ∑
∞ ∞
f (t) = + an cos(nωt) + bn sen(nωt)
2 n=1 n=1
De las dos maneras en que se puede representar una serie de Fourier, podemos decir que
la representación con la amplitud y fase de la n-ésima armónica puede ser obtenida a partir
de la determinación de los coeficientes de Fourier an y bn . Esto es:
√ (a )
n
An = a2n + b2n y ϕn = arc tg
bn
teniendo cuidado sobre la elección del cuadrante.
A continuación, damos las propiedades de las funciones seno y coseno en las siguientes
2π
integrales, en las cuales el periodo está dado por T = .
ω
∫ d+T 0 si n ̸= 0
1) cos(nωt)dt =
d
T si n=0
5
∫ d+T
2) sen(nωt)dt = 0, ∀n
d
∫
0 si m ̸= n
d+T
3) sen(mωt) sen(nωt)dt =
T
d
si m = n ̸= 0
2
∫
0 si m ̸= n
d+T
4) cos(mωt) cos(nωt)dt =
T
d
si m = n ̸= 0
2
∫ d+T
5) cos(mωt) sen(nωt)dt = 0, ∀m, n
d
∫ d+T ∑∞ ∑∞
a0 T a0 T
f (t)dt = dt + (an · 0) + (bn · 0) =
d 2 n=1 n=1
2
∫
a0 1 d+T
= f (t)dt
2 T d
a0
Por el resultado anterior, se constata que el término constante en la expansión en serie
2
de Fourier representa el valor promedio de la función f (t) sobre un periodo. De esta forma,
para una señal eléctrica representa el nivel de referencia o la componente de corriente directa.
Por lo tanto, la expresión que ayuda a determinar el coeficiente a0 viene dada por:
∫ d+T
2
a0 = f (t)dt
T d
6
Como siguiente paso, obtengamos el coeficiente an . Para ello, supongamos que es posible
integrar término a término los elementos de la serie de Fourier. Esto lo consideraremos válido
debido a las propiedades de convergencia de la serie. Este tema es discutido con más detalles
en textos más avanzados.
∞ ∫
∑ d+T
+ bn sen(nωt)cos(mωt)dt
n=1 d
Considerando las propiedades bajo integración de las funciones seno y coseno dadas por
2, 4 y 5 encontramos que cuando m ̸= 0 la única integral distinta de cero del lado derecho es
la que aparece en la primera sumatoria cuando n = m. Entonces tenemos que:
∫ d+T ∫ d+T
an T
f (t)cos(nωt)dt = an cos(nωt)cos(nωt)dt =
d d 2
7
∫ d+T ∫ d+T ∑∞ ∫ d+T
a0
f (t)sen(mωt)dt = sen(mωt)dt + an cos(nωt)sen(mωt)dt
d d 2 n=1 d
∞ ∫
∑ d+T
+ bn sen(nωt)sen(mωt)dt
n=1 d
Al igual que en el caso an , supongamos que se puede integrar término a término y usando
las propiedades bajo integración de las funciones seno y coseno dadas por 2, 3 y 5 encon-
tramos que la única integral distinta de cero del lado derecho es la que aparece en la segunda
sumatoria cuando n = m. Entonces tenemos que
∫ d+T ∫ d+T
bn T
f (t)sen(nωt)dt = bn sen(nωt)sen(nωt)dt =
d d 2
Las ecuaciones que permiten obtener los coeficientes de Fourier reciben el nombre de las
fórmulas de Euler.
A continuación, damos una observación que resulta importante a tener presente al mo-
mento de considerar el valor de la constante d en la búsqueda de los coeficientes de Fourier.
Observación: En las fórmulas de Euler, los límites de integración pueden ser especificadas
sobre cualquier intervalo que contemple un periodo. Con esto, la elección de d es arbitraria y
puede ser hecha de tal manera que pueda ayudar a facilitar en el cálculo de los coeficientes an
y bn . En la práctica, resulta muy frecuente especificar a f (t) sobre cualquiera de los siguientes
intervalos que contemplan un periodo:
1 1
− T ≤t≤ T o 0≤t≤T
2 2
T
dando los valores d = − y d = 0 respectivamente.
2
8
Paso 2: Evaluar la integral en la siguiente identidad:
∫ d+T
2
an + i b n = f (t)einωt dt
T d
Paso 3: Igualar la parte real e imaginaria de cada lado para obtener los valores de los
coeficientes an y bn .
Este método es particularmente resulta sumamente útil cuando solamente queramos de-
terminar la amplitud, |an + ibn |, de la n-ésima armónica.
a0 ∑ ∑
∞ ∞
f (t) = + an cos(nt) + bn sen(nt)
2 n=0 n=0
En la práctica, en muy raras circunstancias nos encontraremos con una función periódica
que presente una frecuencia angular igual uno. Pero este hecho particular resulta sumamente
útil debido a que reduce la dificultad en el proceso de determinación de los coeficientes de an
y bn para la representación en serie de Fourier de la función con periodo angular uno.
9
Ejemplo: Obtengamos la expansión en series de Fourier de la función periódica f de periodo
2π definida por: ( )
f (t) = t 0 < t < 2π , f (t) = f (t + 2π)
1 ( 2π ) 2
=⇒ bn = − cos(2nπ) = −
π n n
porque sen0 = 0 y ∀n = 1, 2, . . . tenemos sen(2nπ) = 0 y cos(2nπ) = 1.
∑∞
2
f (t) = π − sen(nt)
n=1
n
o en la forma expandida
10
( sen(2t) sen(3t) sen(nt) )
f (t) = π − 2 sent + + + ··· + + ...
2 3 n
11
Por su parte, bn está dado por:
∫ π ∫ π
1 1
bn = f (t)sen(nt)dt = (t2 + t)sen(nt)dt
π −π π −π
π 2 ∑ 4π ∑
∞ ∞
n 2π
f (t) = + 2
(−1) cos(nt) + (−1)n+1 sen(nt)
3 n=1
n n=1
n
π2
f (t) = 3
+ 4π(−cost + cos(2t)
4
+ ··· + cos(nt)
n2
+ . . . ) + 2π(sent + sen(2t)
2
+ ··· + sen(nt)
n
+ ...)
12
1
Dado que = −i y utilizando las funciones de Euler dados por:
i
einπ = cos(nπ) + i sen(nπ) = (−1)n
4π 2π 4π 2π
an + i bn = 2
(−1)n − i (−1)n = 2 (−1)n + i (−1)n+1
n n n n
4π 2π
=⇒ an = 2
(−1)n y bn = (−1)n+1
n n
13
∫ ∫ ∫ 2 ]
1 2
1[ 0 1 [ ]2
a0 = f (t)dt = 0dt + dt = t = 1
2 −2 2 −2 0 2 0
∫ 2 ( nπt ) ∫ ( nπt )
1 1 2 1 [ ( nπt )]2
an = f (t) cos dt = cos dt = sen =0
2 −2 2 2 0 2 nπ 2 0
∫ ( nπt ) ∫ ( nπt )
1 2 1 2 1 [ ( nπt )]2
bn = f (t) sen dt = sen dt = − cos
2 −2 2 2 0 2 nπ 2 0
0 si n es par
1
=⇒ bn = [1 − (−1) ] =
n
nπ
2 si n es impar
nπ
2∑ 1
∞ [ (2n − 1)πt ]
f (t) = 1 + sen
π n=1 2n − 1 2
14
Solución: En la Figura 5 se muestra la gráfica de la función f dada en el intervalo −4 < t < 4.
Tomando en consideración el periodo de la función, los coeficientes a0 , an y bn para todo
n = 1, 2, 3, . . . vienen dados por:
∫ 2 ∫ 1 ∫ 2 [ 3t2 ]1
[ ]2 9
2
a0 = f (t)dt = 3tdt + 3dt = + 3t =
2 0 0 1 2 0 1 2
∫ ∫ 1 ∫ 2
2 2
an = f (t) cos(nπt)dt = 3t cos(nπt)dt + 3cos(nπt)dt
2 0 0 1
[ 3tsen(nπt) 3cos(nπt) ]1 [ 3sen(nπt) ]2 3
= + + = [cos(nπ) − 1]
nπ (nπ)2 0 nπ 1 (nπ)2
0 si n es par
=⇒ an =
6
− si n es impar
(nπ)2
∫ ∫ 1 ∫ 2
2 2
bn = f (t) sen(nπt)dt = 3tsen(nπt)dt + 3sen(nπt)dt
2 0 0 1
[ 3tcos(nπt) 3sen(nπt) ]1 [ 3cos(nπt) ]2 3 3
= − + 2
+ − = − cos(2nπ) = −
nπ (nπ) 0 nπ 1 nπ nπ
Finalmente, la serie viene dada por
A fin de comprender a las funciones periódicas continuas a tramos sobre un periodo, ana-
licemos la siguiente la situación y con su ejemplo respectivo, que fueron tomados de James
(2002).
Sea la función f de periodo 2π , cuya gráfica en el intervalo (−π, π) está dada en la Figura
6, y que viene dada por la expresión:
15
f1 (t) si −π < t < −p
f (t) = f2 (t) si −p < t < q
f3 (t) si q<t<π
Las fórmulas de Euler para los coeficientes de Fourier, según la forma como está definida
la función periódica f , se convierten en:
∫ ∫ ∫
1[ −p q π ]
an = f1 (t)cos(nt)dt + f2 (t)cos(nt)dt + f3 (t)cos(nt)dt
π −π −p q
∫ ∫ ∫
1[ −p q π ]
bn = f1 (t)sen(nt)dt + f2 (t)sen(nt)dt + f3 (t)sen(nt)dt
π −π −p q
16
Solución: La Figura 7 muestra la representación gráfica de la función dada f en el intervalo
[−2π, 3π].
Además, ∀ n = 1, 2, . . . se tiene:
∫ 2π
1
an = f (t)cos(nt)dt
π 0
∫ 1 ∫ π ∫ 2π (
1 [ 2π π 1 ) ]
= tcos(nt)dt + cos(nt)dt + π − t cos(nt)dt
π 0 1
π 2 π 2
( 2
)
1 [ tsen(nt) cos(nt) ] π2 [ πsen(nt) ]π [ (2π − t)sen(nt) cos(nt) ]2π
= + + + −
π n n2 0 2n π
2
2n 2n2 π
( )
1 π ( nπ ) 1 ( nπ ) 1 π ( nπ ) cos(2nπ) cos(nπ)
= sen + 2 cos − 2 − sen − +
π 2n 2 n 2 n 2n 2 2n2 2n2
( )
1 ( nπ )
= 2cos − 3 + cos(nπ)
2πn2 2
Por lo tanto,
1 [ n
]
(−1) 2 − 1 si n es par
πn2
an =
−
2
si n es impar
πn2
17
Encontremos ahora el coeficiente bn ∀ n = 1, 2, . . .
∫ 2π
1
bn = f (t)sen(nt)dt
π 0
∫ 1 ∫ π ∫ 2π (
1 [ 2π π 1 ) ]
= tsen(nt)dt + sen(nt)dt + π − t sen(nt)dt
π 0 1
π 2 π 2
( 2
)
1 [ tcos(nt) sen(nt) ] π2 [ πcos(nt) ]π [ (t − 2π)cos(nt) sen(nt) ]2π
= − + + − + −
π n n2 0 2n π
2
2n 2n2 π
( )
1 π ( nπ ) 1 ( nπ ) π π ( nπ ) π
= − cos + 2 sen − cos(nπ) + cos + cos(nπ)
π 2n 2 n 2 2n 2n 2 2n
1 ( nπ )
= sen
πn2 2
Propiedad: Sea f una función dada por f (t) = l g(t) + m h(t), donde g y h son funciones
periódicas de periodo T , l y m son constantes arbitrarias, entonces la función f tiene repre-
sentación en serie de Fourier en donde los coeficientes son la suma de los coeficientes de las
representaciones en series de Fourier de las funciones g y h multiplicados por las constantes
l y m respectivamente.
18
Demostración Como la combinación lineal de funciones periódicas bajo el mismo periodo
siguen siendo periódicas y conservan dicho periodo, la función f es periódica con periodo T . Si
las representaciones en series de Fourier de las funciones g y h están dadas respectivamente
por las expresiones:
a0 ∑ ∑
∞ ∞
g(t) = + an cos(nωt) + bn sen(nωt)
2 n=1 n=1
α0 ∑ ∑
∞ ∞
h(t) = + αn cos(nωt) + βn sen(nωt)
2 n=1 n=1
entonces las evaluaciones de los coeficientes a∗0 , a∗n y b∗n , para n = 1, 2, . . . , en la repre-
sentación en serie de Fourier de f son:
∫ ∫
2 d+T 2 d+T
a∗0 = f (t)dt = [l g(t) + m h(t)]dt
T d T d
∫ ∫
2l d+T 2m d+T
= g(t)dt + h(t)dt = la0 + mα0
T d T d
∫ ∫
2 d+T 2 d+T
a∗n = f (t)cos(nωt)dt = [l g(t) + m h(t)]cos(nωt)dt
T d T d
∫ ∫
2l d+T 2m d+T
= g(t)cos(nωt)dt + h(t)cos(nωt)dt = l an + m αn
T d T d
∫ ∫
2 d+T 2 d+T
b∗n = f (t)sen(nωt)dt = [l g(t) + m h(t)]sen(nωt)dt
T d T d
∫ ∫
2l d+T 2m d+T
= g(t)sen(nωt)dt + h(t)sen(nωt)dt = l bn + m βn
T d T d
l a0 + m α 0 ∑ ∑
∞ ∞
f (t) = + (l an + m αn ) cos(nωt) + (l bn + m βn ) sen(nωt)
2 n=1 n=1
Veamos el siguiente ejemplo, tomado de James (2002), para fijar la proposición anterior.
Ejemplo: Sean las funciones periódicas g y h, con periodo 2π , definidas dentro del intervalo
[−π, π] por:
g(t) = t2 y h(t) = t
f (t) = 2t2 + 3t
19
Solución: Como la función dada por g(t) = t2 es una función periódica par entonces su
representación en serie de Fourier constará de terminos cosenos solamente. Es decir, dado
T = 2π se tiene:
a0 ∑
∞
g(t) = + an cos(nt)
2 n=1
π2 ∑∞
(−1)n
g(t) = +4 2
cos(nt)
3 n=1
n
∑
∞
h(t) = βn sen(nt)
n=1
∑
∞
(−1)n
h(t) = −2 sen(nt)
n=1
n
1 ( 2 2 ) ∑ [ 4(−1)n ]
∞ [
∑ 2(−1)n ]
∞
f (t) = 2· π + 2· cos(nt) − 3· sen(nt)
2 3 n=1
n2 n=1
n
2 ∑∞
(−1)n ∑∞
(−1)n
=⇒ f (t) = π 2 + 8 2
cos(nt) − 6 sen(nt)
3 n=1
n n=1
n
20
1.2.10 Serie de Fourier en su forma compleja
Tomando en cuenta las funciones de Euler, podemos presentar la serie de Fourier en forma
compleja. Esto es, de las funciones complejas:
a0 ∑ ∑
∞ ∞
f (t) = + an cos(nt) + bn sen(nt)
2 n=1 n=1
a0 ∑∞ [ e inωt
+e −inωt ] ∑∞ [ einωt − e−inωt ]
= + an + bn
2 n=1
2 n=1
2i
a0 1 ∑ [ ]
∞
−inωt −inωt
= + an (e inωt
+e ) + bn (ie − ie )
inωt
2 2 n=1
a0 ∑ [ 1 ]
∞
1
= + (an − ibn )einωt + (an + ibn )e−inωt
2 n=1
2 2
Vemos que an − ibn es el conjugado de a + ibn . Por lo que, si designamos los nuevos
coeficientes por las expresiones:
a0 1 1
c0 = , cn = (an − ibn ) y c−n = cn = (an + ibn )
2 2 2
podremos reescribir la serie de Fourier como:
∑
∞ ∑
∞ ∑
∞ ∑
−∞ ∑
∞
inωt −inωt inωt inωt
f (t) = c0 + cn e + c−n e = c0 + cn e + cn e = cn einωt
n=1 n=1 n=1 n=−1 n=−∞
Así tendremos que la serie de Fourier queda representada en forma compleja, también
llamada forma exponencial, mediante la expresión:
∑
∞
f (t) = cn einωt
n=−∞
21
Si tomamos en consideración la forma alternativa de obtención de los coeficientes an y bn ,
de además del hecho de que cn es el conjugado de c−n tenemos:
∫ d+T
a0 1
c0 = = f (t)dt
2 T d
∫
1 d+T
cn = a − ibn = f (t)e−inωt dt
T d
∫ d+T
1
c−n = f (t)einωt dt
T d
con lo cual, se puede apreciar que podemos determinar los coeficientes en la forma compleja
de Fourier mediante
∫ d+T
1
cn = f (t)e−inωt dt, ∀ n = 0, ±1, ±2, ±3, . . .
T d
cn = |cn |eϕn
1√ 2
donde |cn | = an + b2n representa la amplitu de n-ésima armónica y el ángulo ϕn de cn está
2
relacionado con la fase de la n-ésima armónica.
22
Entonces, la serie de Fourier compleja de la función f es:
2a sen(aπ) 2sen(aπ)
a0 = 2
=
aπ aπ
[ 2(−1)n a sen(aπ) ]
an =
(a2 − n2 )π
bn = 0
Si dos funciones dadas por f (t) y g(t) son periódicas y tienen el mismo periodo T entonces
∫ c+T ∑
∞
1
f (t)g(t)dt = c n dn
T c n=−∞
donde cn y dn son los coeficientes en las expansiones en serie de Fourier complejas de las
funciones f y g respectivamente.
23
a0 ∑ ∑
∞ ∞
f (t) = + an sen(nωt) + bn cos(nωt)
2 n=1 n=1
α0 ∑ ∑
∞ ∞
g(t) = + αn sen(nωt) + βn cos(nωt)
2 n=1 n=1
1 1 1 1
cn = (an − ibn ), c−n = (an + ibn ), dn = (αn − iβn ) y d−n = (αn + iβn )
2 2 2 2
1 1 ∑
∞ ∑
∞
= a0 α0 + c−n d−n + cn dn
4 4 n=0 n=0
1 1 ∑∞
1∑
∞
= a0 α0 + (an − ibn )(α + iβn ) + (an + ibn )(α − iβn )
4 4 n=0
4 n=0
1∑
∞
1
= a0 α0 + (2an αn + 2bn βn )
4 4 n=0
Finalmente se tiene, la expresión:
∫
1∑
c+T ∞
1 1
f (t)g(t)dt = a0 α0 + (an αn + bn βn )
T c 4 2 n=0
Si una función dada por f (t) es una función periódica con periodo T , entonces:
∫ c+T ∑
∞ ∑
∞
1 2
[f (t)] dt = cn c∗n = |cn |2
T c n=−∞ n=−∞
24
Ejemplo: Apliquemos el teorema de Parceval a la función
y demostremos que
∑
∞
1
π −8=
2
n=1
(n2 − 41 )2
∑
∞
1 π [ sen(2aπ) ]
= 2π +
−∞
(n2 − a2 )2 4a2 sen2 (aπ) a
1 ∑∞
1 π [ sen(2aπ) ]
+2 = 2 2π +
a4 n=1
(n2 − a2 )2 4a sen2 (aπ) a
∑
∞
1 π [ sen(2aπ) ] 1
= 2π + − 4
n=1
(n − a )
2 2 2 2 2
8a sen (aπ) a 2a
1
Para finalizar, tomemos a la constante a =
2
∑
∞
1 π ∑ 1
∞
= [2π + 0] − 8 =⇒ 1 2 = π −8
2
n=1
(n − 4 )
2 1 2
2 n=1
(n 2 − 4
)
25
1.2.12 Análisis de la convergencia de las series de Fourier
En el análisis que hemos realizado para la representación en serie de Fourier de funciones
periódicas con periodo 2π , hemos supuesto que la serie construida es convergente. En este
apartado veremos ciertas condiciones que se deben dar para que una función pueda ser re-
presentada mediante una serie de Fourier. Kreyszing (2003) enuncia y demuestra la siguiente
proposición que describe las condiciones que deben darse para la representación en serie de
Fourier de una función periódica de periodo 2π .
Los límites por la izquierda y por la derecha, se evaluan en los puntos en los cuales la
función sufre un corte, tal como lo muestra la Figura 8. Por lo tanto, si x = a es un punto
de discontinuidad de la función dada por y = f (x), los límite por la izquierda y la derecha
respectivamente vienen dados por:
f (a − h) − f (a − 0) f (a + h) − f (a + 0)
lı́m y lı́m
h−→0 −h h−→0 h
mediante valores de x positivos.
Figura 8: Representación gráfica de los límites por la izquierda y por la derecha de una función.
Fuente: [Link]
26
[Link] La condición débil de Dirichlet para las series
La condición débil de Dirichlet: garantiza la existencia de la serie de Fourier si los coefi-
cientes de Fourier; a0 , an y bn para n = 1, 2, 3, . . . , son finitos. Por lo tanto, la integral del valor
absoluto de la forma funcional de la señal periódica en un periodo debe ser finita. Es decir,
sea una señal periódica, con periodo T , dada por f (t) entonces
∫ d+T
|f (t)| < ∞
d
Las condiciones fuertes de Dirichlet, junto con la condición débil, establecen que una
función priódica y con periodo T dada por f (t) podrá ser representada en serie Fourier si se
satisfacen las siguientes dos condiciones:
Ejemplo: Analicemos si las siguientes funciones satisfacen las condiciones débil y fuertes de
Dirichlet en el intervalo [−π, π].
1
1. La función dada por f (t) =
t2
2. La función dada por f (t) = sen(ln t)
1
1) La función dada por f (t) = presenta una discontinuidad infinita en t = 0, que esta
t2
dentro del intervalo [−π, π], y por lo tanto no satisface la condición fuerte de Dirichlet
que que establece que la función periódica solo debe tener un número finito de discon-
tinuidades finitas dentro del periodo considerado.
b) La función dada por f (t) = sen(ln t) tiene un número infinito de máximos y mínimos en
una vecindad de t = 0, que está dentro del intervalo [−π, π], y por lo tanto no satisface
la condición fuerte de Dirichlet que establece que la función debe tener sólo un número
finito de máximos y mínimos aislados dentro de un periodo.
Con lo expuesto se puede decir que las condiciones de Dirichlet son suficientes para ase-
gurar la representación en serie de Fourier de una función periódica, pero no se puede concluir
acerca de si estas condiciones sean necesarias para la convergencia. En la práctica, debido
a que resulta muy dificil crear señales como las del ejemplo anterior en un laboratorio, las
condiciones de Dirichlet se satisfacen para casi todas las señales y por lo tanto admiten re-
presentación en serie de Fourier.
27
1.2.13 Derivación e integración de series de Fourier
En el análisis matemático, la derivada e integración constituyen dos condiciones desables
y en ocasiones sumanente útiles en las aplicaciones de las series de Fourier sobre señales.
Según James (2002):
Propiedad: Una expansión en serie de Fourier de una función periódica dada por f (t) que
satisface las condiciones de Dirichlet puede integrarse término a término, y la serie integrada
converge a la integral de dicha función.
Según la proposición anterior, si la función periódica, de periodo T , dada por f (t) satisface
las condiciones de Dirichlet en un intervalo [t1 , t] ∈ [d, d + T ], con t1 < t, y tiene una expansión
en serie de Fourier dada por:
a0 ∑ ∑
∞ ∞
f (t) = + an cos(nt) + bn sen(nt)
2 n=1 n=1
La expresión anterior, estando de esa menera, no constituye una seire de Fourier. Esto
a0 t
hecho se da debido al término en el lado derecho de la identidad. Al reordenar la identidad
2
podemos convertirla nuevamente en una representación en serie de Fourier de la función dela
siguiente manera:
∫
a0 t1 ∑ an [sen(nt) − sen(nt1 )] ∑ bn [cos(nt1 ) − cos(nt)]
t ∞ ∞
a0 t
f (t)dt − =− + +
t1 2 2 n=1
n n=1
n
Para fijar la proposición anterior, veamos el siguiente ejemplo tomado de James (2002).
28
f (t + 2π) = f (t)
y su representación en serie de Fourier:
π2 ∑∞
(−1)n cos(nt)
f (t) = t2 = +4 para −π ≤t≤π
3 n=1
n2
[ u3 ]t π 2 [ ]t ∑∞
(−1)n [ ]t
= u +4 sen(nu)
3 −π 3 −π
n=1
n2 −π
t3 π2t ∑∞
(−1)n sen(nt)
= +4
3 3 n=1
n3
∑
∞
(−1)n sen(nt)
t − π t = 12
3 2
n=1
n3
Por lo tanto, el lado derecho puede tomarse como la expansión en serie de Fourier de la
función dada por:
g(t) = t3 − π 2 t para − π ≤ t ≤ π
g(t + 2π) = g(t)
Analicemos ahora la derivación de una serie de Fourier. Para ello, veamos la siguiente
porposición mencionada por James (2002).
Propiedad: Si una función periódica dada por f (t) satisface las condiciones de Dirichlet,
entonces su derivada, f ′ (t) siempre que exista, puede encontrarse por derivación término a
término de su representación en serie de Fourier si y solo si dicha función es continua en todas
partes y la función derivada, f ′ (t), tiene expansión en serie de Fourier, es decir que también
satisface las condiciones de Dirichlet.
De la proposición anterior, si la función priodica dada por f (t), de periodo T , satisface las
condiciones de Dirichlet en en el intervalo [d, d + T ] y su representación en serie de Fourier es
a0 ∑ ∑
∞ ∞
f (t) = + an cos(nωt) + bn sen(nωt)
2 n=1 n=1
29
entonces la serie anterior puede derivarse miembro a miembro si la derivada de la función f
satisface las condiciones de Dirichlet de la siguiente manera:
(∑
∞ ∑
∞ )′
′
f (t) = an cos(nωt) + bn sen(nωt)
n=1 n=1
[∑
∞ ∑
∞ ]
= ω n bn cos(nωt) − n an sen(nωt)
n=1 n=1
Se puede apreciar que los coeficientes nbn y nan son los coeficientes de Fourier en la
representación en serie de la derivada de la función f , así, en contraste con la serie integrada,
la serie derivada converge más lentamente que la expansión en serie original de f (t).
π2 ∑∞
(−1)n cos(nt)
2
t = +4
3 n=1
n2
Solución: Vemos que la función periódica f es continua dentro y en los extremos del intervalo
−π ≤ t ≤ π y satisface en dicho intervalo las condiciones de Dirichlet, por lo tanto podemos
aplicar la proposición anterior entonces:
∑
∞
(−1)n+1 sen(nt)
t=2 para −π ≤t≤π
n=1
n
30
1.3.1 Funciones pares e impares
Las funciones de acuerdo a su paridad pueden clasificarse en funciones pares, impares
o ninguna de las dos. Comencemos el análisis de las formas funcionales periódicas, viendo
como vienen dadas las funciones pares e impares.
A partir de este momento será importante tener presente que si f es una función par
entonces paro todo valor de la variable t se tiene que f (t) = f (−t) y en su representación
gráfica, Figura 9, se observará que la función presenta simétrica con respecto al eje vertical.
Uno de los casos más relevalentes de funciones pares resulta ser f (t) = cost. Una propiedad
importante de las funciones pares se da en la integración. Por ejemplo, si la función f es
simétrica con respecto al eje vertical, de la definición de la integral definida se sigue que para
toda b > 0: ∫ ∫
b b
f (t)dt = 2 f (t)dt
−b 0
Por su parte, si f es una función impar entonces para todo valor de la variable t se tiene
que f (t) = −f (−t) y la representación gráfica de la función, Figura 10, es simétrica con
respecto a un punto sobre el eje horizontal, esto es, se observará una simetría de cuadrante
opuesto. De la misma manera que el caso de funciones pares, mencionemos que la función
f (t) = sent es una de las funciones impares más relevantes y que una propiedad importante
de las funciones impares se da en la integración. Por ejemplo, si f una función impar simétrica
con respecto al origen, considerando la definición de integral definida se tiene que para todo
b > 0: ∫ b
f (t)dt = 0
−b
31
Figura 10: Representación gráfica de una función f impar en torno al origen de coordenadas.
Fuente: Elaboración propia
Propiedad: Cualquier función f se puede expresar como la suma de dos funciones, de las
cuales una de ellas es una función par y la otra es impar.
Si f es una función en la variable t, entonces para todo t del dominio de definición, se tiene
que:
1 1 1 1
f (t) = f (t) + f (t) + f (−t) − f (−t)
2 2 2 2
1 1
= [f (t) + f (−t)] + [f (t) − f (−t)]
2 2
= fp (t) + fi (p)
1 1
fp (−t) = [f (−t) + f (− − t)] = [f (t) + f (−t)] = fp (t)
2 2
1 1
fi (−t) = [f (−t) − f (− − t)] = − [f (t) − f (−t)] = −fi (t)
2 2
con lo cual fp y fi son funciones pares e impares respectivamente.
32
Las funciones pares e impares presentan la siguientes propiedades:
1) La única función que es par e impar a la vez, es la función idénticamente nula, f (t) = 0
para todo t.
2) La suma de una función par y una impar no será ni par ni impar, a menos de que una de
las ellas sea la función nula.
10) El cociente de una función par y una función impar es una función impar
Por lo tanto, la expansión en serie de Fourier de una función periódica par con periodo
T contiene solo términos cosenos, es decir que para todo n = 0, 1, 2, . . . se tiene:
∫ T
4 2
an = f (t)cos(nωt)dt
T 0
a0 ∑
∞
f (t) = + an cos(nωt)
2 n=1
33
Si la función periódica f es impar y tiene periodo T , entonces usando la propiedades 6
y 7, tenemos los coeficientes a0 , an y bn para todo n = 1, 2, . . . vienen dada por:
∫ T
2 2
an = f (t)cos(nωt)dt = 0
T − T2
∫ T ∫ T
2 2 4 2
bn = f (t)sen(nωt)dt = f (t)sen(nωt)dt
T − T2 T 0
De esta manera, la expansión en serie de Fourier de una función periódica impar con
periodo T contiene términos senos solamente, es decir:
∑
∞
f (t) = bn sen(nωt)
n=1
Ejemplo: Una función periódica f con periodo 2π está definida para t ∈ [−π, π] por:
−1 si −π ≤ t < 0
f (t) =
1 si 0≤t≤π
∫ ∫ [ ]π
2 ( )
T
π
4 2 4 2 cos(nt)
bn = f (t)sen(nt)dt = sen(nt)dt = − =− con(nπ) − 1
T 0 2π 0 π n πn
0
34
Figura 11: Representación gráfica de una onda cuadrada impar en el intervalo [−4π, 4π]
Fuente: Elaboración propia
4( ) 4∑ ∞
cos(3t) cos(5t) sen[(2k − 1)t]
f (t) = cost + + + ... =
π 3 5 π k=1 (2k − 1)
35
Figura 12: Representación gráfica de una forma de onda par en el intervalo [−3π, 3π]
Fuente: Elaboración propia
π2 ∑∞
(−1)n cos(nt)
f (t) = +4
3 n=1
n2
De esta manera, una función f , en la variable t, periódica y con periodo T se dice que
tiene simetría de media onda si cumple con la condición dada por:
T
f (t) = f (t ± )
2
36
Figura 13: Representación gráfica de una función con simetría de media onda
Fuente: Elaboración propia
Figura 14: Representaciones gráficas de funciones con simetría de cuarto de onda par e impar.
Fuente: [Link]
37
La simetría escondida en una cierta señal se puede evidenciar mediante la realización de
desplazamientos sobre uno de los ejes, según el caso lo requiera, de manera a conseguir una
transformación de la señal periódica dada en otra señal que tenga una de las simetrías, par o
impar. Con este procedimiento, podemos evaluar de manera más sencilla los coeficientes de
Fourier para luego volver a la señal original y así obtener su serie de Fourier correspondiente.
El siguiente ejemplo servirá para ilustrar este situación.
Solución: Observando la Figura 15, constatamos que la función de onda f , con periodo T ,
es una función impar y se define de la siguiente manera:
1 t
f (t) = − con 0 < t < T
2 T
38
Por otra parte, su expansión en serie de Fourier será de la forma:
∑
∞
f (t) = bn sen(nωt)
n=1
1 ∑ sen(nωt)
∞
f (t) =
π n=1 n
Observando la Figura 16, podemos notar que la función g , con periodo T , viene dada por:
t 1
g(t) = 1 − = + f (t), con 0 < t < T
T 2
y su representación en serie de Fourier será:
1 1 ∑ sen(nωt)
∞
g(t) = +
2 π n=1 n
39
Si una función f es periódica, con periodo T , tal que:
(T)
f t+ = −f (t)
2
entonces solo las armónicas impares están presentes en su expansión en serie de Fou-
rier. Por lo tanto, para n = 1, 3, 5, . . . tenemos:
∫ T ∫ T
4 2 4 2
an = f (t)cos(nωt)dt y bn = f (t)sen(nωt)dt
T 0 T 0
∫ T
8 4
a2n−1 = f (t)cos[(2n − 1)ωt]dt
T 0
∫ T
8 4
b2n−1 = f (t)sen[(2n − 1)ωt]dt
T 0
Solución: Analizando la gráfica de la función de onda cuadrada, podemos notar que presenta
simetría de cuarto de onda par y que se define de la siguiente manera:
40
Figura 17: Representación gráfica de una función de onda cuadrada con simetría de orden par.
Fuente: Elaboración propia
T T
−1 si − ≤t<−
2 4
T T
f (t) = 1 si − ≤t<
4 4
T T
−1 si ≤t<
4 2
2π
Con lo cual, la frecuencia angular está dada por ω = y
T
∫ T ∫ T
8 4 8 4
a2n−1 = f (t)cos[(2n − 1)ωt]dt = 1 cos[(2n − 1)ωt]
T 0 T 0
( )
[ 8 sen[(2n − 1)ωt] ] T4 4 [ π]
= = sen (2n − 1) − 0
(2n − 1)ωT 0 (2n − 1)π 2
Entonces: 4
si n es impar
(2n − 1)π
a2n−1 =
−
4
si n es par
(2n − 1)π
41
Ejemplo: Encontremos la expansión en serie de Fourier de la función de onda cuadrada de
T T
periodo T definida el rango − < t < y que se muestra en la Figura 18.
2 2
Figura 18: Representación gráfica de una función de onda cuadrada con simetría de cuarto de onda impar.
Fuente: Elaboración propia
Solución: Analizando la gráfica de la función de onda cuadrada, podemos notar que presenta
simetría de cuarto de onda impar y que se define de la siguiente manera:
T
−1
si −
2
≤t<0
f (t) =
1 si 0≤t<
T
2
2π
Por lo que, la frecuencia angular está dada por ω = y
T
∫ T ∫ T
8 4 8 4
b2n−1 = f (t)sen[(2n − 1)ωt]dt = 1 sen[(2n − 1)ωt]
T 0 T 0
( )
[ 8 cos[(2n − 1)ωt] ] T4 4 [ π]
= − = cos (2n − 1) − 1
(2n − 1)ωT 0 (2n − 1)π 2
4
= −
(2n − 1)π
4 ∑ sen[(2n − 1)ωt]
∞
f (t) = −
π n=1 2n − 1
42
1.5 Aproximación mediante una serie finita de Fourier
En la práctica, cuando estamos frente a señales periódicas se trata de buscar su repre-
sentación en serie de Fourier de tal forma a tomar después solamente algunos términos de la
serie de tal forma que dicha porción serial genere lo más razonablemente el comportamiento
de la señal. Por esta razón, uno de los aspectos importantes del análisis de la convergengia
de las series de Fourier en las aplicaciones prácticas es la llamada razón de convergen-
cia. Esta taza, proporciona una idea de cuantos términos resulta adecuado tomar de la serie
trigonométrica infinita para representar de manera adecuada la señal tomada.
Como una serie de Fourier queda especificada por sus coeficientes an y bn , la cantidad
de téminos a considerar de dicha serie para aproximar el comportamiento de la señal consid-
erada depende exclusivamente del comportamiento de ellos. En este sentido, James (2002)
menciona:
Si la función f y todas sus derivadas hasta de r-ésimo orden son continuas pero la (r+1)-
ésima derivada es discontinua entonces los coeficientes en su representación en serie
1
de Fourier decrecen conforme y por lo tanto, la serie converge más rápidamente y
nr+2
se podrá tomarse unos pocos términos de la serie para representar a la función.
donde ϵ es el error cometido al aproximar f por medio de fN . Una de las maneras de estimar
el valor del error cometido es mediante el error cuadrático total que viene dado por:
43
∫ d+T
ET = [f (t) − fN (t)]2 dt
d
y que origen al error cuadrático medio:
∫ d+T
1
EM = [f (t) − fN (t)]2 dt
T d
Se puede demostrar que si se aproxima una función f por una serie finita de Fourier fN ,
entonces esta aproximación tiene la propiedad de generar el mínimo error cuadrático medio y
este error viene dada por la expresión:
∫
a2 1 ∑ 2
d+T N
∗ 1
EM (t) = [f (t)] dt − 0 −
2
(a + b2n )
T d 4 2 n=1 n
∑
∞
sen(nt)
f (t) = π − 2
n=1
n
2
con a0 = 2π , an = 0 y bn = − para n = 1, 2, . . . . Con lo cual, el error cuadrático medio
n
viene dado por:
∫
4π 2 1 ∑ 4
2π 4
∗ 1
EM = t2 dt − −
2π 0 4 2 n=1 n2
4π 2 ( 1 1 1 )]
= −π −2 1+ + +
2
3 4 9 16
π 2 205
= − = 0, 443
3 72
44
g(t) = f (t) en d ≤ t ≤ d + T, g(t + T ) = g(t)
que tendrá una representación en serie de Fourier convergente.
Podemos apreciar que dentro del intervalo [d, d + T ] las funciones g y f tienen compor-
tamientos idénticos, por lo que la representación en serie de la función g en dicho intervalo
será también válido para f .
f (t) = t
Figura 19: Representación gráfica de la función f en el intervalo [0, 4] y una posible extensión periódica.
Fuente: Elaboración propia
45
∫
1 4
1 [ t2 ]4
a0 = tdt = =4
2 0 2 2 0
( nπt ) ( nπt )
∫ 4 ( ) [ 2tsen 4cos ]4
1 nπt 2 2
an = t cos dt = + = 0, ∀ n = 1, 2, . . .
2 0 2 nπ (nπ)2 0
( nπt ) ( nπt )
∫ 4 ( ) [ 2tcos 4sen ]4
1 nπt 2 2 4
bn = t sen dt = − + 2
= − , ∀ n = 1, 2, . . .
2 0 2 nπ (nπ) 0 nπ
Además, como g(t) = f (t) en el intervalo 0 < t < 4 podemos notar que la serie de Fourier
obtenida converge a la función f solamente este intervalo, por lo que:
( )
4 ∑∞ sen nπt
2
f (t) = 2 − si 0<t<4
π n=1 n
Si la función dada por f (t) está definida sobre un intervalo finito [0, λ], para poder con-
seguirle una representación en serie de Fourier que contenga términos cosenos solamente,
primero f debe cumplir las condiciones de Dirichlet en el intervalo dado y segundo tendremos
que extenderla mediante la función periódica par g dada por:
f (t) si 0<t<λ
g(t) = g(t + 2λ) = g(t)
f (−t) si −λ < t < 0
Como la función g es periódica par y de periodo 2λ, entonces tiene una representación
convergente en serie de Fourier que consista solamente de términos cosenos y viene dada
por:
a0 ∑
∞ ( nπt )
g(t) = + an cos
2 n=1
λ
46
Los coeficientes an se determinan mediante la expresión:
∫ τ ( nπt )
2
an = f (t)cos dt
λ 0 λ
Debido a que, dentro del intervalo 0 < t < λ, la función g es idéntica a f , se sigue que la
serie de términos cosenos converge a f dentro de dicho intervalo.
De manera análoga al caso anterior, si la función dada por f (t) está definida sobre un inter-
valo finito [0, λ], para poder conseguirle una representación en serie de Fourier con términos
senos solamente, primero f debe cumplir las condiciones de Dirichlet en el intervalo dado y
segundo tendremos que extenderla mediante la función periódica impar G dada por:
f (t) si 0<t<λ
G(t) = G(t + 2λ) = G(t)
−f (−t) si −λ < t < 0
∑
∞ ( nπt )
G(t) = bn sen
n=1
λ
Debido a que dentro del intervalo [0, λ], G(t) es idéntica a f (t), se tiene que la serie obteni-
da converge a f dentro de dicho intervalo.
Como las series de Fourier de términos cosenos y senos se han construido en base a la
función no periódica f son llamadas expansiones en series de Fourier de medio recorrido de la
función f . Así, la expansión par de medio recorrido g es llamada expansión en serie de medio
recorrido en cosenos de la función f y la expansión impar de medio recorrido G es llamada
expansión en serie de medio recorrido en senos de la función f .
47
Dibujemos además las gráficas de la función f y de las funciones periódicas representadas
por las dos series obtenidas en el intervalo −8 < t < 12.
Solución: a) Para obtener la serie de medio recorrido en cosenos, definamos la función pe-
riódica g dada por:
f (t) = t si 0<t<4
g(t) = F (t + 8) = F (t)
f (−t) = −t si −4 < t < 0
Como g es una función periódica par con periodo 8 y satisface las condiciones de Dirichlet,
tiene expansión en serie de Fourier convergente con a0 y an para n = 1, 2, . . . dadas por:
∫
2 4
a0 = tdt = 4
4 0
( ) ( )
∫ 4 ( ) [ 4tsen nπt
16cos nπt ]
2 nπt 1 4 4 4 8
an = tcos dt = + = [cos(nπ) − 1]
4 0 4 2 nπ (nπ)2 0 (nπ)2
16
− si n es impar
(nπ)2
=⇒ an =
0 si n es par
Como g(t) = f (t) en el intervalo [0, 4], tenemos que la serie de Fourier converge a f
solamente dentro de este intervalo y por lo cual:
[ ]
(2n−1)πt
16 ∑ cos
∞
4
f (t) = t = 2 −
π2 n=1
(2n − 1)2
48
( ) ( )
∫ 4 ( nπt ) [ nπt
4tcos 4 16sen nπt ]4
2 1 4 8cos(nπ)
bn = tsen dt = − + 2
=−
4 0 4 2 nπ (nπ) 0 (nπ)2
8
=⇒ bn = (−1)n+1
nπ
Como G(t) = f (t) en el intervalo [0, 4], tenemos que la serie de Fourier es convergente a
f dentro de este intervalo. Así, la expansión en serie de medio recorrido en senos de f es:
( )
8 ∑ (−1) sen 4
∞ n+1 nπt
49
1.6 Integral de Fourier y espectros continuos
Hemos visto que en la práctica, la representación en serie de Fourier de una señal pe-
riódica nos permite realizar importantes análisis de comportamiento de los mismos. Pero en la
realidad se pueden presentar también señales aperiódicas. De este hecho, surge la necesidad
de extender la teoría de Fourier mediante la llamada integral de Fourier.
Antes de entrar a detallar todo lo que concierne a la integral de Fourier resulta importante
comprender cuando una función es absolutamente integrable y además, en que consiste una
transformación integral.
es absolutamente integrable
Como la función dada es par, ya que f (−t) = f (t), entonces al integrar el valor absoluto
de la función tenemos:
∫ ∞ ∫ ∞ ∞
|f (t)|dt = 2 e−t dt = −2e−t =2<∞
−∞ 0 0
Una transformación integral consiste en la operación matemática que asocia a cada función
f , en la variable t en un espacio real o corrido, otra función F , en la variable τ en el espacio
recíproco mediante la siguiente relación:
∫ b
F (τ ) = k(τ, t)f (t)dt
a
donde la función k(τ, t) recibe el nombre de kernel de la transformación, y los límites de inte-
gración a y b están dados por la transformada correspondiente.
50
Las transformadas integrales constituyen herramientas matemáticas sumamente útiles en
situaciones en las que un determinado problema resulta ser muy difícil de resolver en su
espacio real o directo; pero al transformarlo al espacio recíproco ya es más sencillo resolverlo.
Después de ello, la transformada inversa nos devuelve la solución en el espacio original.
Si una función, absolutamente integrable, dada por f (t) es continua por secciones en
todo intervalo finito y tiene una derivada por la izquierda y una por la derecha en todo punto
entonces puede representarse por una integral de Fourier. En cualquier punto donde la función
f sea discontinua, el valor de la integral es igual al promedio de los límites por la izquierda y
por la derecha de f en ese punto.
51
y su representación en integral coseno de Fourier es:
∫ ∞
f (t) = A(ω)cos(ωt) dω
0
por lo tanto, admite representación en integral de Fourier y el coeficiente A(ω) viene dado por:
∫
2 10
2 [( t2 2) 2t ]10
A(ω) = 2
t cos(ωt)dt = − 3 sen(ωt) + 2 cos(ωt)
π 0 π ω ω ω 0
2 [( 100 2) 20 ]
=⇒ A(ω) = − 3 sen(10ω) + 2 cos(10ω)
π ω ω ω
Con el resultado anterior, la representación en integral de Fourier es:
∫ [( 100
2 ∞
2) 20 ]
f (t) = − 3 sen(10ω) + 2 cos(10ω) cos(ω t)dω
π 0 ω ω ω
52
Ejemplo: Encontremos la integral de Fourier sobre el pulso encendido-apagado, con A y T
finitos y mayores a cero, dada por:
A si − T ≤ t ≤ 0
f (t) = −A si 0 ≤ t ≤ T
0 en otro caso
por lo tanto, admite representación en integral de Fourier y el coeficiente B(ω) viene dado por:
∫ T ]T
2A 2A
B(ω) = − sen(ωt)dt = cos(ωt)
π 0 πω 0
2A
=⇒ B(ω) = [cos(ωT ) − 1]
πω
53
∫ ∞
cos(ω) − 1 π
sen(ω)dω = −
0 ω 4
∫ ∞ ∫ ∞
cos(ω)sen(ω) sen(ω) π
dω − dω = −
0 ω 0 ω 4
∫ ∞ ∫ ∞
sen(2ω) sen(ω) π
dω − dω = −
0 2ω 0 ω 4
∫ ∞ ∫ ∞
1 sen(τ ) sen(ω) π
dτ − dω = −
2 0 τ 0 ω 4
∫ ∞
1 sen(ω) π
− dω = −
2 0 ω 4
∫ ∞
sen(ω) π
dω =
0 ω 2
54
Bibliografía
1 James, G. 2002. Matemáticas avanzadas para ingeniería. 2da edición. México: Pearson
Educación.
2 Spiegel, M., Lipschutz, S. & Schiller, J. (2011). Variable compleja (Ebook). (2◦ Ed.). Mé-
xico: McGraw-Hill Interamericana.
3 Brown, J.W. & Churchill, R. V. (2005).Variable compleja y aplicaciones. (7◦ Ed.). Madrid:
McGraw-Hill.
5 Kreyszing, E. (2003). Matemática avanzada para ingeniería. (Vol. I). México: Limusa S.A.
55