Universidad Nacional de Ingeniería

Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones

__________________________________________________________________________

INFORME DE LABORATORIO DE
SISTEMA DE COMUNICACIONES I
- EE 430 N -

UNIVERSIDAD NACIONAL DE INGENIERIA

FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

DOCENTE: ING. VIRGINIA ROMERO F.

ALUMNOS:
1. BLAS PERALTA CARLOS SEBASTIÁN (20174538K)
2. CABALLA CASTRO, CHRISTIAN (20172532E)
3. ALEJANDRO AYALA, HENRY (20171280B)
4. AGUILAR ESPINOZA, EMILIO (20171211K)
5. MARÍN SANTARIA, LUIS (20171415E)
ESPECIALIDAD: L2

2021-I
_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

1
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

INFORME FINAL N° 1

SIMULACIÓN DE LA SERIE DE FOURIER MEDIANTE EL

SOFTWARE MATLAB
Docente: Ing. Virginia Romero Fuentes
Grupo N°: 1
Alumnos:
1. Blas Peralta, Carlos Sebastián ; email: [email protected]
2. Caballa Castro, Christian ; email: [email protected]
3. Alejandro Ayala, Henry ; email: [email protected]
4. Aguilar Espinoza, Emilio ; email: [email protected]
5. Marín Santaria, Luis Ángel ; email: [email protected]
Resumen

En el presente Laboratorio se comprobó con ayuda del software Matlab la obtención de la serie de
Fourier como una suma de funciones senoidales y cosenoidales con distintos múltiplos de la
frecuencia fundamental. Asimismo, comprobamos teóricamente la serie de Fourier al igual que las
gráficas de espectro de frecuencia.
Terminada la experiencia coincidimos que el Software de Matlab nos va a ser de mucha ayuda en los
futuros Laboratorios de Sistemas de Comunicaciones 1 por su amplia variedad de aplicaciones que
podemos hacer en ella tanto para graficar como calcular.
I. OBJETIVOS

• Hallar gráficamente la aproximación de una onda por medio de la sumatoria de “n” términos de la
serie de Fourier.
• Graficar el espectro de frecuencias y demostrar analíticamente lo que representa.
• Demostrar que el Software de Matlab es el adecuado para las experiencias de Laboratorio de Sistemas
de Comunicaciones.

II. ANALISIS Y PROCEDIMIENTO

A nuestro grupo se le asignó una función Pulso cuadrado impar con amplitud 20𝑉𝑝𝑝 y periodo de 20
ms.
Con ayuda del software Matlab se realizaron las indicaciones señaladas en la guía de Laboratorio y
comprobaron teóricamente.

a) Visualización de la función en el Programa Matlab.

Pulso impar en Matlab

fs = 10000;
t = 0:1/fs:50;
x1 = 10*square(2*pi*50*t);
_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

2
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

subplot(2,1,1)
plot(t,x1)
axis([0 0.2 -11 11])
xlabel('Time (sec)');
ylabel('Amplitude');
title('Pulso Cuadrado impar');

Figura 1. Gráfica del pulso cuadrado impar

b) Desarrollo Analítico del espectro de frecuencias para la señal asignada.

Para esto podemos usar en Matlab la transformada rápida de Fourier.
clear all, clc, close all
Fs=1000;
T=1/Fs;
L=1500;
t=(0:L-1)*T;
A=10;
A2=0.5;
F1=50;
w0=2*pi*F1;
y=10*square(2*pi*F1*t);
figure (1)
plot(t(1:100),y(1:100))

f=Fs*(0:(L/2))/L;
Y=fft(y);
P2=abs(Y/L);
P1=P2(1:L/2+1)*2;
figure(2)
plot(f,P1);

k=3;
s=0;
for n=1:k
_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

3
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

p=(-20*(cos(n*pi)-1)/(n*pi))*sin(n*w0*t);
s=s+p;

end
figure(3)
plot(t(1:100),s(1:100))

Figura 2. Gráfica del espectro de frecuencia de la función cuadrado impar

Comentario (1): Al tratarse de un pulso impar se esperaba que solo tuviera armónicos cada 100Hz,
con una frecuencia fundamental de 50Hz como también se observa en la gráfica.
También podemos obtener de manera analítica el espectro de frecuencia de la señal cuadrada
usando la serie de Fourier.
40 1 1
𝑥(𝑡) = [sin(𝑛𝑤𝑜 𝑡) + sin(𝑛𝑤𝑜 𝑡) + sin(𝑛𝑤𝑜 𝑡) + ⋯ ]
𝜋 3 5
De la serie de Fourier se observa que:
40
• Para el primer armónico de 50Hz se observa una amplitud de 𝜋
≈ 12.73
• Para el segundo armónico de 100Hz se observa una amplitud de 0 al tratarse de un coeficiente
par.
40
• Para el tercer armónico de 150Hz se observa una amplitud de 3𝜋 ≈ 4.24
• Para el cuarto armónico de 200Hz se observa una amplitud de 0 al tratarse de un coeficiente
par.
40
• Para el quinto armónico de 250Hz se observa una amplitud de 5𝜋 ≈ 12.54
_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

4
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

• Para el sexto armónico de 300Hz se observa una amplitud de 0 al tratarse de un coeficiente par.
Entonces podemos obtener el espectro de frecuencia tanto por la serie de Fourier, así como en el
Software Matlab.

c) Dependencia de “n”

Se desarrolló el programa en Matlab donde se observa la dependencia del valor de los n

armónicos que se desea en la transformada de Fourier.
clear all
clc
%x(t)
T=0.2;
w=(2*pi)/T;
t=(-T)/2:0.001:T/2;
xt=(t>=0).*20+(t>=(-T)/2).*-10;
plot(t,xt,'b','LineWidth',3)
hold on
grid on
%serie de furier
k=5;
s=0;
for n=1:k
p=(-20*(cos(n*pi)-1)/(n*pi))*sin(n*w*t);
s=s+p;

end
plot(t,s)

Figura 2. Gráfica de la serie de Fourier para n=5

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

5
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

Figura 3. Gráfica de la serie de Fourier para n=10

Figura 4. Gráfica de la serie de Fourier para n=15

Figura 5. Gráfica de la serie de Fourier para n=20

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

6
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

Comentario (2): Se observa en el gráfico como a medida que aumentamos “n” o sumandos
a la serie de Fourier se va aproximando a la función cuadrada impar (grafica azul)

d) Especificaciones del programa.

Hallando 𝑎0 , 𝑎𝑛 𝑦 𝑏𝑛 :

Figura 6. Gráfica de la función pulso cuadrado impar

𝑻
−𝟏𝟎 𝒑𝒂𝒓𝒂 −
<𝒕<𝟎 𝟐𝝅
𝒙(𝒕) = { 𝟐 ; 𝒘𝒐 =
𝑻 𝑻
𝟏𝟎 𝒑𝒂𝒓𝒂 𝟎 < 𝒕 <
𝟐

∞
1
𝑥(𝑡) = 𝑎0 + ∑(𝑎𝑛 cos(𝑛𝑤𝑜 𝑡) + 𝑏𝑛 sin(𝑛𝑤𝑜 𝑡))
2
𝑛=1

𝑇 𝑇
0
2 2 2
𝑎0 = ∫ 𝑥(𝑡) = [∫ (−1)𝑑𝑡 + ∫ (1)𝑑𝑡] = 0
𝑇 − 𝑇
−
𝑇
0
2 2

𝒂𝟎 = 𝟎
𝑇 𝑇
2 2 2 0 2
𝑎𝑛 = ∫ 𝑥(𝑡) cos(𝑛𝑤𝑜 𝑡) 𝑑𝑡 = [∫ (−1) cos(𝑛𝑤𝑜 𝑡) 𝑑𝑡 + ∫ (1)cos(𝑛𝑤𝑜 𝑡)𝑑𝑡]
𝑇 − 𝑇 𝑇 − 𝑇
0
2 2

𝟐
𝒂𝒏 = [𝟎] = 𝟎
𝑻
𝑇 𝑇
2 2 2 0 2
𝑏𝑛 = ∫ 𝑥(𝑡) sin(𝑛𝑤𝑜 𝑡) 𝑑𝑡 = [∫ (−1) sin(𝑛𝑤𝑜 𝑡) 𝑑𝑡 + ∫ (1)sin(𝑛𝑤𝑜 𝑡)𝑑𝑡]
𝑇 −𝑇 𝑇 −𝑇 0
2 2

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

7
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

𝑇
2 0 2𝜋 2 2𝜋
𝑏𝑛 = [∫ (−1) sin (𝑛 𝑡) 𝑑𝑡 + ∫ (1)sin(𝑛 𝑡)𝑑𝑡]
𝑇 −𝑇 𝑇 0 𝑇
2

𝟒𝟎
𝟐[𝐜𝐨𝐬(𝒏𝝅) − 𝟏] 𝒑𝒂𝒓𝒂 "𝒏" 𝒊𝒎𝒑𝒂𝒓
𝒃𝒏 = − = { 𝝅 ;𝒏 ≠ 𝟎
𝒏𝝅 𝟎 𝒑𝒂𝒓𝒂 "n" 𝒑𝒂𝒓

4 1 1
𝑥(𝑡) = [sin(𝑤𝑜 𝑡) + sin(3𝑤𝑜 𝑡) + sin(5𝑤𝑜 𝑡) + ⋯ ]
𝜋 3 5
Simulación en Matlab
clear all
clc
%x(t)
T=0.2;
w=(2*pi)/T;
t=(-T)/2:0.001:T/2;
xt=(t>=0).*20+(t>=(-T)/2).*-10;
plot(t,xt,'b','LineWidth',3)
hold on
grid on
%serie de furier
k=50;
a0=0;
An=0;
s=0;

for n=1:k
Bn=-20*(cos(n*pi)-1)/(n*pi);
p=(1/2)*a0+An*cos(n*w*t)+(Bn)*sin(n*w*t);
s=s+p;
end
plot(t,s)

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

8
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

Figura 6. Gráfica de la serie de Fourier para n=50

Figura 7. Gráfica de la serie de Fourier para n=200

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

9
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

Figura 8. Gráfica de la serie de Fourier para n=1000

Comentario 2: Con este segundo método de obtener la serie de Fourier es donde mas es exijido
nuestro computador debido a la cantidad de operaciones que esta haciendo la funcio de Matlab de
integrar, además se aprecia visualmente el fenómeno de Gibbs en los extremos de la onda cuadrada
para “n” menores a 50.
III. RESULTADOS

Valor Teórico Valor Experimental Error % Error

N=3 150Hz 0.86 86%

N=5 250Hz 0.76 76%

N=10 500Hz 0.52 52%

N=15 750Hz 0.29 29%

N=20 1000Hz 0.047 4.76%

Tabla N° 3.1

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

10
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

IV. RESPUESTAS A PREGUNTAS

1. De lo observado en el computador y del desarrollo teórico, aproxime el ancho de banda de

la señal que se le asignó.

El ancho de banda se define como el lugar donde se encuentra la mayor parte de Potencia, en el
espectro. gráficamente el ancho de banda optimo es entre 900Hz y 1100Hz,

En la gráfica de espectro de amplitud se usa el criterio del 95%, ósea se buscará la frecuencia donde
la amplitud sea aproximadamente 5% del pico (12,73), que para este caso es 0,636.

Figura. Criterio del 95% para ancho de banda

Este se genera aproximadamente en el armónico n=21.

Entonces el Ancho de banda=1050Hz

2. ¿Qué diferencias importantes observa entre la función ideal y la obtenida por la serie
de Fourier?
• Se observa que, a pesar de contar una mayor cantidad de armónicos, no se logra eliminar por
completo las ondulaciones que producen las funciones seno y coseno. Sin embargo, es necesario
recalcar que a partir de un número de armónicos se aproxima notablemente a la función ideal.
• Se puede apreciar el fenómeno de Gibbs para los 50 primeros términos en la serie de Fourier.
• El máximo de la onda obtenida mediante la serie de Fourier mantiene una forma ondulada, un
poco diferente a la función ideal donde es un pico.

3. ¿Cuántos armónicos fueron necesarios para una mejor aproximación de la función que se le
asignó?

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

11
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

Figura “a” Figura “b”

Figura “c”

En figura “a” hemos dado n=5 y notamos una lejanía a la función pulso cuadrado impar, en la Figura
“b” para n=50 ya notamos una aproximación más exacta, en la Figura “c” para n=1000 ya se nota
una aproximación exacta.

4. ¿Qué consideraciones adicionales requirió para su simulación?

Se empleó la función escalón en Matlab. Asimismo, la función for para iterar la suma de los
armónicos y el empleo de matrices en Matlab. Se consideró la diferencia entre una gráfica continua
(para la función aproximada) y una gráfica discreta (para el espectro en frecuencias).
5. ¿Tuvo problemas con la nueva versión de Matlab recientemente instalada en su PC?
El principal problema es debido a que Matlab es un software de paga, entonces es necesario instalar
la versión de prueba, la cual tiene una licencia por tiempo limitado. Por el resto de casos, las
funciones utilizadas no presentaron problemas en la versión utilizada.
_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

12
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

6. ¿Logro cumplir los objetivos de la práctica?

Si se logró cumplir con los objetivos, con el gráfico del espectro de frecuencias se aproximó el ancho
de banda. Asimismo, el código empleado permitió aproximar la sumatoria de “n” términos de la serie
de Fourier. El código estaba optimizado para solo modificar el número de armónicos y obtener la
aproximación de Fourier.
V. OBSERVACIONES Y CONCLUSIONES

El grupo presentara una observación y una conclusión por alumno integrante del grupo.
Alejandro Ayala, Henry Iván
Ob.1
Fue necesario tener conocimientos previos de programación para realizar la experiencia.
Con.1
Se logró hallar gráficamente la aproximación de la onda asignada para distintos valores de "n"
en la serie de Fourier.
Caballa Castro, Christian Harold
Ob.1
Para poder alcanzar el objetivo de este proyecto fue necesario trabajar grupalmente, porque
para esta experiencia se necesita conocimientos básicos de Matlab y las series de Fourier.
Con.1
Gracias al trabajo colaborativo, de nuestro grupo, pudimos ejecutar satisfactoriamente las
simulaciones y comparaciones pedidas en este laboratorio.
Con.2
Pudimos apreciar más de cerca el potencial que tiene la herramienta Matlab.

Aguilar Espinoza, Emilio Abraham

Ob.1
MATLAB nos facilitó mucho el tema de las presentaciones de los gráficos en el tiempo y
frecuencia, con herramientas como axis, grid y las opciones de la interfaz gráfica, las cuales
permiten hallar puntos máximos, mínimos y tiempos notables.
Con.1
Se comprobó que a un mayor número de “n” la gráfica se aproxima más a la señal original,
pero ocurre el fenómeno de Gibbs.
Con.2
La herramienta MATLAB es muy útil a la hora de entender y realizar series de Fourier así
como mostrar gráficos detallados de los mismos.

Blas Peralta, Carlos Sebastián

Ob.1
Hallar los coeficientes de la Serie de Fourier mediante Matlab por las integrales definidas es
un método que implica muchos recursos en el computador haciendo que se demore la
visualización de la gráfica.
Con.1
_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

13
 Universidad Nacional de Ingeniería
Facultad de Ingeniería Eléctrica, Electrónica y Telecomunicaciones
__________________________________________________________________________

La función cuadrado impar necesita de un alto ancho de banda en el espectro de frecuencia

para aproximarse a la función real, lo que la hace una función que al ser expresada como serie
de Fourier necesita de muchas funciones senoidales y siempre presentara desviaciones en los
extremos evidenciándose el fenómeno de Gibbs.

Marín Sanitaria, Luis Ángel

Ob.1
Es necesario conocer los coeficientes de la serie de Fourier de la función impulso cuadrada
impar para poder realizar la simulación.
Con.1
He podido notar que mientras el “n” aumenta se acerca cada vez más la función impulso
cuadrada impar.

VI. BIBLIOGRAFÍA USADA

Colocar la bibliografía que han usado para presentar su informe final

[1] Robert L. Boylestad and Louis Nashelsky, Electroni devices and circuit theory, 10th edition.
[2] G. E. TOBEY, "Operational Amplifiers; Design and Applications", Publisher: McGraw-
Hill.
[3] L. STRAUSS, “Wave generation and shapping”, Mc. GRAW HILL 1977.

[4] JACOB MILLMAN, "Microelectronics: Digital and Analog Circuits and S

[5] Fuente del navegador https:// es.wikipedia.org/wiki/Espectro_de_frecuencias
[6] Fuente del navegador https:// es.wikipedia.org/wiki/Joseph_Fourier

[7] José E. Briceño Márquez, “Principios de las Comunicaciones”

VII. FECHA
Lima, 12/05/2021

_____________________________________________________________________________________________________________________________ ________________________________________________________________________________________________________________

Ing. Virginia Romero F. UNI – FIEE

14