CLASE 16

SERIE COMPLEJA DE FOURIER Y APLICACIONES

Problematización inicial.
¿Qué es una señal exponencial compleja?
¿En qué consiste la serie compleja de Fourier?
¿Cómo calcular los coeficientes de la serie compleja?
¿En qué consisten los espectros discretos de magnitud y fase?
¿Cómo dibujar los espectros de magnitud y fase?
¿En qué consisten las funciones circuitales en el dominio de la frecuencia compleja?
¿Cómo determinar la respuesta de estado estacionario de un sistema lineal invariante?

Contenido.
En este módulo se presentan los fundamentos de las series complejas de Fourier y sus
aplicaciones en el cálculo de la respuesta de estado estacionario de un sistema lineal invariante.
Se hace énfasis en las representaciones gráficas de los espectros discretos de magnitud y fase.
Los temas a desarrollar son los siguientes:

1. La función exponencial compleja.

2. Serie compleja de Fourier.
3. Espectros discretos de magnitud y fase.

1. La función exponencial compleja.

La función exponencial compleja es de uso frecuente en el análisis de los sistemas lineales
invariantes en el tiempo. Con base en la identidad de Euler, dicha función se puede expresar en la
forma:
jωt
e =cos(¿ ωt)+ jsen( ωt) ¿

La conjugada asociada a la función exponencial viene dada por:

− jωt
e =cos (¿ ωt)− jsen(ωt )¿

A partir de las dos funciones descritas previamente se pueden expresar las señales elementales
seno y coseno, así:

1 jωt − jωt
cos (¿ ωt )= [ e +e ] ; sen (ωt )= 1 [ e jωt −e− jωt ] ¿
2 j2

En general, si ω=n ω 0 es un entero, se verifica que:

1 jn ω t − jn ω t
cos (¿ n ω 0 t)=[e +e 0
]¿0

2
1 jn ω t − jn ω t
sen(n ω 0 t)= [ e −e 0
] 0

j2

2. Serie compleja de Fourier.

Con las expresiones anteriores es posible determinar la serie de Fourier compleja de una función
periódica, así:
Dada una función periódica f (t) que satisface las condiciones de Dirichlet, se estudió en el
capítulo anterior que su serie trigonométrica está dada por:
 a0 ∞ ∞
f (t)= + ∑ an cos(¿ n ω0 t)+ ∑ b n sen(n ω0 t)¿
2 n=1 n=1

Los coeficientes de la serie se calculan como:

t 0+T
2
a 0=
T ∫ f (t) dt
t0

t 0+T
2
a n=
T
∫ f (t) cos(¿ n ω0 t)dt ; n ≥ 1¿
t0
t 0+T
2
b n=
T ∫ f (t) sen (n ω0 t)dt ; n ≥ 1
t0

Teniendo en cuenta las funciones complejas: seno y coseno, descritas anteriormente, se puede
escribir la serie de la siguiente manera:

a0 ∞ 1
∞
1
f (t)= + ∑ an ⋅ [ e jnω t + e− jn ω t ] + ∑ b n ⋅ [ e jn ω t−e− jn ω t ]
0 0 0 0

2 n=1 2 n=1 j 2

Teniendo en cuenta el hecho: 1/ j=− j y reuniendo términos semejantes, resulta:

a0 1 ∞ jnω t 1
∞
f (t)= + ∑ ( an− jb n ) e + ∑ ( an + jb n ) e− jn ω t
0 0

2 2 n=1 2 n=1

Si denotamos los coeficientes, por:

a0 an− j bn an + j bn
c 0= , c n= , c −n =
2 2 2
Resulta la forma compleja de la serie de Fourier, así:

∞
f (t)= ∑ c n e jn ω t 0

n=−∞

∞
f ( t )=c0 + ∑ c n e jn ω t
0

n =−∞
n≠ 0

2π
ω 0=
T

1
c n= ¿
2

1
c n= ¿
T
 [∫ ]
t 0 +T
1 − jnω0 t
c n= f (t )e dt
T t0

[∫ ]
t 0+T
1
c 0= f ( t ) dt
T t0

El coeficiente c n de la serie compleja es un número complejo, es decir, se puede caracterizar

mediante dos maneras diferentes:

Forma estándar:
c n= [ ℜ ( n ) ] + j [ ℑ(n) ]

Forma polar:
c n=|c n| e ; n ∈ Z
j θn

|cn|=√ [ ℜ(n)] + [ ℑ(n)]

2 2
; θn =arctan ( ℜ(n)
ℑ(n)
)
Programa en Matlab para calcular los coeficientes de la serie compleja.
clear all
%seriedeFouriercompleja
%entrada de información
syms t n
T=input('periodo=');
ti=input('tinicial=');
f=input('entre la función f(t)=');
w=2*pi/T;
c0=(1/T)*int(f,ti,T+ti);
cn=(1/T)*int(f*exp(-i*n*w*t),ti,T+ti);
c0
cn
pretty(c0)
pretty(cn)

Ejemplo 1.
Una señal periódica está definida en un periodo como:

f ( t )=t [ u ( t )−u (t−2) ]

2

a) Dibuje dos ciclos de la señal.

b) Escriba la serie sin calcular los coeficientes.
c) Calcule los coeficientes usando el programa anexo y usar el editor de ecuaciones.
d) Escriba los coeficientes en forma simplificada (parte real y parte imaginaria).
e) Escriba magnitud y fase.

Solución.
a) La gráfica la hacemos con Matlab, así:
t=0:0.01:2;
y=t.^2.*(heaviside(t)-heaviside(t-2));
plot(t,y,'k')
hold on
t=-2:0.01:0;
y=(t+2).^2.*(heaviside(t+2)-heaviside(t));
plot(t,y,'k')
grid on
La figura muestra la gráfica pedida.

b) La serie, sin calcular los coeficientes, es:

∞
f ( t )=c0 + ∑ c n e jn π t
n =−∞
n≠ 0
c) Ejecutando el programa, se tiene:
>> fouriercompleja2024
periodo=2
tinicial=0
entre la función f(t)=t^2*(heaviside(t)-heaviside(t-2))
c0 =
4/3
cn =
1i/(n^3*pi^3) + (exp(-pi*n*2i)*(pi^2*n^2*4i + 4*pi*n - 2i))/(2*n^3*pi^3)

2 2
1i exp(-pi n 2i) (pi n 4i + pi n 4 - 2i)
------ + ---------------------------------------
3 3 3 3
n pi 2 n pi

Teniendo en cuenta que: exp (− j2 nπ ) =1, se tiene que:

c 0=4 /3
2+ j2 n π
c n= 2 2
,n≠0
n π
d) Puede verse que:
2
ℜ ( n )= 2 2
n π
2nπ
ℑ ( n )= 2 2
n π
e) Magnitud y fase.
 mag ( n )= √
4 +4 n2 π 2
n2 π2
fase ( n )=atan ⁡(n π )

3. Espectros discretos de magnitud y fase.

Las series trigonométricas de Fourier expresan a las señales periódicas como la suma de señales
armónicas. La serie compleja es una suma de fasores, armónicamente relacionados, donde es
importante analizar la amplitud y la fase de cada armónica. Estas armónicas existen para cada una
de las frecuencias armónicas:

ω=0 , ± ω0 , ±2 ω 0 ,± 3 ω 0 , ± 4 ω0 , ⋯ ⋯ ⋯ ,± n ω 0

En la sección anterior se encontró que una función periódica, que satisface las condiciones de
Dirichlet, se puede expresar mediante la serie compleja de Fourier:

∞
f (t)= ∑ c n e jn ω t0

n=−∞

Los coeficientes de la serie son números complejos y viene dados por:

t 0 +T
1
c n=
T
∫ t 2 e− jn ω t dt 0

t0

Dado que los coeficientes son números complejos, se podrán expresar mediante su magnitud y su
fase, así:
an bn
− j =|c n|e ; n ∈ Z
jθ
c n= n

2 2

Tanto |cn| como θn son funciones de la variable discreta n cuyo dominio es el conjunto de los
números enteros.

Espectro discreto de magnitud.

El espectro discreto de magnitud es una gráfica de la magnitud en función de: n ; n ∈ Z , y se
caracteriza por ser una función par.

Espectro discreto de fase.

El espectro discreto de fase es una gráfica de la fase en función de: n ; n ∈ Z y se caracteriza por
ser una función impar.
Las gráficas se obtienen con Matlab de acuerdo con el siguiente código:

Programa en Matlab para dibujar los espectros discretos de magnitud y fase.

%espectrosdiscretos
clear all
c0=input('c0=');
n=0;
stem(n,c0)
hold on
N=input('número de armónicos=');
n=1:N;
R=input('real=');
X=input('imaginaria=');
mag=sqrt(R.^2+X.^2);
fase=atan(X./R);
stem(n,mag)
hold on
stem(-n,mag)
grid on
figure
stem(n,fase)
hold on
stem(-n,-fase)
grid on

Ejemplo.
Para la función del ejemplo anterior, se tiene:
>> espectrosdiscretos2024
c0=4/3
número de armónicos=10
real=2./(pi^2*n.^2);
imaginaria=2./(pi*n);
Las figuras se muestran a continuación.

TAREA
Una señal periódica está definida en un periodo como:

f ( t )=t [ u (t )−u (t−1) ] +u ( t−1 )−u (t−2)

a) Dibuje dos ciclos de la señal.

b) Escriba la serie sin calcular los coeficientes.
c) Calcule los coeficientes usando el programa anexo y usar el editor de ecuaciones.
d) Escriba los coeficientes en forma simplificada (parte real y parte imaginaria).
e) Escriba magnitud y fase.