GRUPO N°1

Universidad Nacional de Ingenierı́a

Facultad de Ingenierı́a Eléctrica y
Electrónica

CURSO:
SISTEMAS DE COMUNICACIONES I
EE430 - O

INFORME PREVIO N°1

SIMULACIÓN DE LA SERIE DE FOURIER
MEDIANTE EL SOFTWARE MATLAB

Alumno:
• CIPRIANO YURIVILCA, CRISTOPHER JHOSMEL

September 12, 2023

LIMA - PERÚ

1
 SIMULACIÓN DE LA SERIE DE FOURIER
MEDIANTE EL SOFTWARE MATLAB

Objetivos–Encontrar graficamente la aproximación de una Fourier desarrolló esta teorı́a durante su estudio de la
onda por medio de la sumatoria de n términos de la Serie de ecuación del calor y publicó sus resultados en 1807 y 1811, es-
Fourier y estimar su ancho de banda. tableciéndose como el pionero en su investigación sistemática.
Las series de Fourier encuentran aplicación en diversas áreas,
I. F UNDAMENTO T E ÓRICO que incluyen ingenierı́a, por ejemplo en telecomunicaciones,
La serie de Fourier es una representación infinita que se se utilizan para optimizar la concepción de sistemas de trans-
emplea para funciones que son periódicas y continuas o que misión de señales mediante el análisis de los componentes
mantienen continuidad en segmentos especı́ficos. Estas series, espectrales de frecuencia de una señal mediante un analizador
nombradas en honor al matemático francés Jean-Baptiste Jo- de espectro. También existe otra alternativa de representar las
seph Fourier, descomponen una función en una suma infinita funciones periódicas, utilizando una combinación de exponen-
de funciones sinusoidales más simples, como senos y cosenos ciales complejas la cual es la serie compleja de Fourier a
con frecuencias enteras, esta función se expresa en la siguiente diferencia de la serie trigonométrica esta se representa como:
ecuación (1).
+∞
X a0
∝ 
f (t) = Cn ejnωt ∧ C0 = (5)
ao X 2nπ 2nπ

−∞
2
f (t) = + an cos( t) + bn sen( t) (1)
2 n=1
T T Donde el coeficiente Cn se calcula:
Donde f (t) es una función periódica de [−T /2, T /2] y a0 , 1 T
Z
an y bn son respectivamente. Cn = f (t)e−jnωt (6)
T 0
Z T /2 Cada coeficiente esta relacionado con exponencial compleja
2
a0 = f (t)dt (2) la acumulación de varias contribuyen a la reconstrucción de la
T −T /2 función periódica original f(t), gracias a su simplicidad ayuda a
Z T /2  optimizar cálculos como por ejemplo en el estudio del espectro
2 2nπ
an = f (t) cos t dt (3) de frecuencias de señales periódicas.
T −T /2 T

2 T /2
Z  
2nπ
bn = f (t) sen t dt (4)
T −T /2 T

Figura 2. Espectro de señal

II. E QUIPOS Y M ATERIALES

Una computadora
Software de Matlab
Acceso a internet
Figura 1. Función periodica
 III. P ROCEDIMIENTO DE SIMULACI ÓN Entonces el seno es una función periódica de periodo
1. Se utilizo el software MATLAB para la simulación T = 2π
B
previamente se realizaron los cálculos correspondientes f (t + T ) = f (t)
para encontrar la obtención de los coeficientes de su
serie de Fourier. También se le asocia una frecuencia ω = 2π
T = B.
2. Se utiliza un arreglo para definir el tiempo y los saltos Para el caso del coseno sucede lo mismo:
entre que valores (Periodo de muestreo), se definen f (t) = Acos(Bt + C)
valores como periodo, amplitud y frecuencia.
3. Una vez obtenido la serie de Fourier se utiliza una f (t + T ) = Acos(B(t + T ) + C)
sentencia For desde el valor inicial de n y finalizando en
2π
una variable defina k como limite, a su ves guardamos Para T = B
los valores de la suma en la variable S la cual guarda los
2π 2π
valores que toma f para luego graficarlo con el arreglo f (t + ) = Acos(B(t + )+C
creado del tiempo. B B
4. Se asigna diferentes valores a la variable de k con la f (t + T ) = Acos(Bt + 2π + C) = cos(Bt + C)
finalidad de comprobar su aproximación ası́ como el
fenómeno de Gibbs, etc. Entonces el seno es una función periódica de periodo
T = 2π
B
f (t + T ) = f (t)
IV. R ESPUESTAS A PREGUNTAS
1. ¿Cuáles son las caracterı́sticas de una función pe- Gráficamente se logra visualizar que la función
riódica? fundamental sen(t) y cos(t) tienen un periodo de
Las funciones periódicas son funciones continuas que T = 2π.
cumple la siguiente condición:
3. Detallar las propiedades de los coeficientes de Fou-
f (x + T ) = f (x) (7) rier: Identidad de Parseval; y la relación entre los
coeficientes de Fourier y su derivada.
Donde T es el periodo de la función. Mediante la La potencia promedio una señal viene representada
repetición de la ec.3 se puede definir f (x) como: como la suma de cuadrados de los coeficientes de
f (x + 2T ) = f (x + T + T ) = f (x + T ) = f (x) la serie de Fourier, según la identidad de Parseval.

En general: T /2 ∞
a0 2 1 X  2
Z
1
[f (t)] 2 = an + bn 2 (10)

+
f (x) = f (x + nT ), n = 0, ±1, ±2, ... (8) T −T /2 4 2 n=1
Lo cual puede ser representado mediante la serie de La linealidad, si una señal esta compuesta en base a
Fourier. Las funciones exponenciales complejas: la suma de otras señales sus coeficientes presentan
igual comportamiento, es decir, la suma de sus
sk (t) = eikw0 t = ei2πkf0 t , k = 0, ±1; ±2; ... (9) coeficientes de las señales parciales dan origen a
los coeficientes de la señal compuesta.
También son funciones periódicas ya que están relacio- Los coeficientes de Fourier están relacionados con
nadas armónicamente por tener todas un periodo común la derivada de la misma señal, ya que tendrá es
Tp = f10 . mimo periodo de la señal anterior, por ejemplo.
2. ¿La función seno y coseno son funciones periódicas?
La función trigonométrica de seno es periódica por ∝  
satisfacer la ecuación de una función periódica. ao X 2nπ 2nπ
f (t) = + an cos( t) + bn sen( t)
2 n=1
T T

se procede a derivar termino a termino

f (t) = Asen(Bt + C)
f (t + T ) = Asen(B(t + T ) + C) X∝ 
2nπ 2nπ

2π
f´(t) = −an ω0 sen( t) + bn ω0 cos( t)
Para T = B n=1
T T
2π 2π Entonces los nuevos coeficientes son:
f (t + ) = Asen(B(t + )+C
B B
f (t + T ) = Asen(Bt + 2π + C) = Asen(Bt + C) a∗n = bn ω0 ∧ bn ∗ = −an ω0
 Si calculamos lo coeficientes en función a su señal
1
f ´(t) se tiene lo siguiente. fT /2 (x) = (f (x+) + f (x−))
2
2 T /2 ´
Z
f (t)cos(ω0 t) = bn ω0 Donde:
T −T /2 f (x+) = lı́m (f (t))
t→x+
Z T /2
2 f (x−) = lı́m (f (t))
f ´(t)sen(ω0 t) = −an ω0 t→x−
T −T /2

4. Definir y explicar detalladamente la serie de Fourier,

determinar los coeficientes de una función.
La serie de Fourier es una suma infinita que se utiliza
para representar funciones periódicas continuas en el
intervalo del periodo la siguiente expresión es la serie
de Fourier con coeficientes reales.

∝  
ao X 2nπ 2nπ
f (t) = + an cos( t) + bn sen( t)
2 n=1
T T
(11)
Donde f (t) es una función periódica de [−T /2, T /2] de
periodo T > 0 con a0 , an y bn coeficientes de la serie
de Fourier su calculo es el siguiente.
Z T /2
2
a0 = f (t)dt (12)
T −T /2

Z T /2  
2 2nπ
an = f (t) cos t dt (13)
T −T /2 T
Z T /2  
2 2nπ
bn = f (t) sen t dt (14)
T −T /2 T
5. Explicar detalladamente la serie de Fourier, deter-
minar los coeficientes de Dirichilet y el teorema de
convergencia.
Las condiciones de Dirichilet son necesarias para la
representación de una función f (t) mediante la serie de
Fourier:
La función f (t) tiene un número finito de discon-
tinuidades en un periodo.
La función f (t) tiene un numero finito de máximos
y mı́nimos en un periodo.
la integral del valor absoluto de f (t) en un periodo
es finita: Z T /2
|f (t)|dt < ∞ (15)
−T /2

Los coeficientes de Fourier para una función de periodo

T son los siguientes:
Z T /2
2 2nπx
an = f (x) cos dx
T −T /2 T
Z T /2
2 2nπx
bn = f (x) sen dx
T −T /2 T
Por lo que la función convergerá.
6. Explicar el fenómeno de Gibbs
La convergencia de la serie de Fourier de una fun-
ción cerca de una discontinuidad de salto presenta
caracterı́sticas especiales que fueron investigadas por
el matemático J. N. Gibbs alrededor del año 1899.
Si consideramos una función f con perı́odo 2π y la
secuencia Sn de sumas parciales de su serie de Fourier,
Gibbs hizo la siguiente observación: en las proximidades
de una discontinuidad de salto, Sn siempre difiere de
f en una cantidad aproximada al 9 % el tamaño del
salto. Esta propiedad se denomina como el fenómeno
de Gibbs.
Sea la función f de periodo 2π

1
f (x) = (π − x)
2 Figura 5. S9

Figura 3. Funcion periodica

Figura 6. S19
Vemos las aproximaciones en series de Fourier:
n 7. Desarrolle análiticamente el espectro de la frecuen-
X sin(kx)
Sn = (16) cias para la señal-asignada
k
k=1 La señal asigna es de tipo pulso de onda cuadrada impar
con una amplitud de 10 V pp, con periodo de T = 30
(mseg) y la duración del impulso 10 (mseg) lo cual viene
representado gráficamente :

Figura 7. Señal de pulso cuadrado impar

Donde sus coeficientes en la serie de Fourier son los

siguientes:
20
a0 = (17)
3
Figura 4. S5
10 2nπ
an = sin (18)
nπ 3
 10 2nπ
bn = (1 − cos ) (19) R EFERENCIAS
nπ 3
[1] Duoandikoetxea, J. (2003). Lecciones sobre las series y
Para analisar el espectro de frecuencias se utiliza el An transformadas de Fourier. Recuperado de http://www. ugr.
y nω0 donde: es/acanada/docencia/matematicas/analisisdefourier/Duoandiko-
etxeafourier. pdf.
p [2] Bruzual, R., & Dominguez, M. (2003). Series de fourier. Universidad
An = a2n + b2n (20) Central de Venezuela-Año.
[3] Hsu, H. P., Mehra, R. (1982). nálisis de Fourier. México: Sistemas
Entonces el espectro de frecuencias es el siguiente: Técnicos de Edición.
[4] Santos, R. J. (2002). Séries de Fourier. Minas Gerais: Departamento de
Matemática-ICEx.

Figura 8. Espectro de frecuencias

8. Simulación de la serie de Fourier

La serie de Fourier de la señal asignada es de la siguiente
forma:

∞  
10 X 10 2nπ nπ 10 2nπ nπ
f (t) = + sin( ) cos( t) + (1 − cos( ) sin( t))
3 n=1 nπ 3 15 nπ 3 15
(21)
Si tomamos como un valor limite k se tiene:

k  
X 10 10 2nπ nπ 10 2nπ nπ
f (t) = +[ sin( ) cos( t) + (1 − cos( ) sin( t))
n=1
3k nπ 3 15 nπ 3 15
(22)
Por lo cual su simulación es la siguiente:

Figura 9. Serie de Fourier