Matemáticas Aplicadas a la

Electrónica
Laboratorio N°14
Serie de Fourier
Alumno:
Navarro Quispe, Alexander
Ramos Huamán, Rogger Smit
Parco Guerra, Alonzo
Ciclo: III Sección: 3C5

Profesor: Godinez De La Cruz, Ernesto.

Fecha de realización: 19/06/23
Fecha de presentación: 25/06/23

2023-1

1
Objetivo

El objetivo del presente laboratorio consiste en obtener la serie de Fourier de una

señal periódica

Fundamentos Teóricos

Si es una señal periódica, de periodo T, frecuencia f, amplitud A y frecuencia

angular ω. La serie de Fourier es una función expresada en términos seno y

coseno, la cual es aproximadamente igual a la función

Los coeficientes de la serie de Fourier son

Procedimiento

Ejercicio 1:

Se tiene una señal periódica , cuya regla de correspondencia para un periodo

es

a) Encontrar los coeficientes de la serie de Fourier para n = 1, 2, …,10

2
 b) Graficar la función y la serie de Fourier con 5, 8 y 10 armónicos
5 armónicos:

3
8 armónicos:

10 armónicos:

4
 c) Expresar la serie de Fourier en términos de t con 5 armónicos

Código:

//Laboratorio 14
//Transformada de Laplace
//Parco Guerra, Alonzo Fernando
//Pregunta 1
clc
clear
close
//***************************
ti=0
tf=6
T=tf-ti //periodo
w=2*%pi/T //frecuencia angular
h=0.00001
t=ti:h:tf
f1=5
f2=-3
t1=(t>=0)&(t<=4)
t2=(t>4)&(t<=6)
f=f1.*t1+f2.*t2 //señal
//Cálculo de a0
a0=(1/T)*inttrap(t,f)
disp('Coeficientes de la serie de Fourier')
mprintf('\n a0 = %1.4f',a0)
//Cálculo de an y bn
na=5 //numero de armónicos
s=0
for n=1:na
a(n)=(2/T)*inttrap(t,f.*cos(n*w*t))
b(n)=(2/T)*inttrap(t,f.*sin(n*w*t))
s=s+a(n)*cos(n*w*t)+b(n)*sin(n*w*t)
end
5
 ff=a0+s
for n=1:na
mprintf('\n a(%d)=%1.4f',n,a(n))
end
for n=1:na
mprintf('\n b(%d)=%1.4f',n,b(n))
end
//Gráfica de f(t) y la serie de Fourier
plot(t,f,t,ff)
xlabel('t')
title('Gráficas de f y ff')
xgrid
legend('señal', 'serie de Fourier')
//Expresión de la serie de Fourier
na=5
disp('Serie de Fourier')
mprintf('\n ff(t) = %1.4f',a0)
for n=1:na
mprintf(' %+1.4fcos(%dwt)+%1.4fsen(%dwt)',n,a(n),b(n),n)
end

Ejercicio 2:

Se tiene una señal periódica , cuya regla de correspondencia para un periodo

es

6
 a) Encontrar los coeficientes de la serie de Fourier para n = 1, 2, …,10

b) Graficar la función y la serie de Fourier con 5, 8 y 10 armónicos

5 armónicos:

7
8 armónicos:

10 armónicos:

c) Expresar la serie de Fourier en términos de t con 5 armónicos

8
Código:
//Laboratorio 14
//Transformada de Laplace
//Ramos Huamán,Rogger Smit
//Pregunta 2
clc
clear
close
//***************************
ti=0
tf=6
T=tf-ti //periodo
w=2*%pi/T //frecuencia angular
h=0.00001
t=ti:h:tf
f1=2
f2=5
f3=3
t1=(t>=0)&(t<=2)
t2=(t>2)&(t<=4)
t3=(t>4)&(t<=6)
f=f1.*t1+f2.*t2+f3*t3 //señal
//Cálculo de a0
a0=(1/T)*inttrap(t,f)
disp('Coeficientes de la serie de Fourier')
mprintf('\n a0 = %1.4f',a0)

9
//Cálculo de an y bn
na=5 //numero de armónicos
s=0
for n=1:na
a(n)=(2/T)*inttrap(t,f.*cos(n*w*t))
b(n)=(2/T)*inttrap(t,f.*sin(n*w*t))
s=s+a(n)*cos(n*w*t)+b(n)*sin(n*w*t)
end
ff=a0+s
for n=1:na
mprintf('\n a(%d)=%1.4f',n,a(n))
end
for n=1:na
mprintf('\n b(%d)=%1.4f',n,b(n))
end
//Gráfica de f(t) y la serie de Fourier
plot(t,f,t,ff)
xlabel('t')
title('Gráficas de f y ff')
xgrid
legend('señal', 'serie de Fourier')
//Expresión de la serie de Fourier
na=5
disp('Serie de Fourier')
mprintf('\n ff(t) = %1.4f',a0)
for n=1:na
mprintf('\n %+1.4fcos(%dwt)+%1.4fsen(%dwt)',n,a(n),b(n),n)

10
end

Ejercicio 3:

Se tiene una señal periódica , cuya regla de correspondencia para un periodo

es

a) Encontrar los coeficientes de la serie de Fourier para n = 1, 2, …,10

11
 b) Graficar la función y la serie de Fourier con 5, 8 y 10 armónicos
5 armónicos:

8 armónicos:

10 armónicos:

12
 c) Expresar la serie de Fourier en términos de t con 5 armónicos

Código:

//Laboratorio 14
//Transformada de Laplace
//Navarro Quispe, Alexander
//Pregunta 3
clc
clear
close
//***************************
ti=-2
tf=2
T=tf-ti //periodo
w=2*%pi/T //frecuencia angular
h=0.00001
t=ti:h:tf
f1=exp(t)
f2=exp(-t)
t1=(t>=-2)&(t<=0)
t2=(t>0)&(t<=2)

13
f=f1.*t1+f2.*t2 //señal
//Coeficientes de Fourier
a0=(1/T)*inttrap(t,f)
na=5 //Número de armónicos
s=0
for n=1:na
a(n)=(2/T)*inttrap(t,f.*cos(n*w*t))
b(n)=(2/T)*inttrap(t,f.*sin(n*w*t))
s=s+a(n)*cos(n*w*t)+b(n)*sin(n*w*t)
end
ff=a0+s
disp('')
disp('Valores de an')
for n=1:na
mprintf('\n a(%d)=%1.5f',n,a(n))
end
mprintf('\n')
disp('Valores de bn')
for n=1:na
mprintf('\n b(%d)=%1.5f',n,b(n))
end
//Gráficas de f y ff
plot(t,f,t,ff,'linewidth',3)
xlabel('t','fontsize',4)
title('Gráficas de f y ff','fontsize',4)
legend('Señal','Fourier','fontsize',4)
xgrid

14
//Serie d Fourier
mprintf('\n')
disp('Serie de Fourier')
mprintf('\n ff(t) = %1.5f',a0)
for n=1:na
mprintf('\n %+1.5fcos(%dwt)+%1.5fsen(%dwt)',a(n),n,b(n),n)
end

Aplicaciones

1) Se tiene una señal periódica , cuya regla de correspondencia para un

periodo es

a) Encontrar los coeficientes de la serie de Fourier para n = 1, 2, ...,10

15
 b) Graficar la función y la serie de Fourier con 5, 8 y 10 armónicos

5 armónicos:

8 armónicos:

10 armónicos:

16
 c) Expresar la serie de Fourier en términos de t con 5 armónicos

Código:

//Laboratorio 14
//Transformada de Laplace
//Parco Guerra, Alonzo Fernando
//Aplicación 1
clc
clear
close
//***************************
e=%e
ti=0
tf=2
T=tf-ti //periodo
w=2*%pi/T //frecuencia angular
h=0.00001
t=ti:h:tf
f1=t^2
f2=(t-2)^2
t1=(t>=0)&(t<=1)
t2=(t>1)&(t<=2)
17
f=f1.*t1+f2.*t2 //señal
//Coeficientes de Fourier
a0=(1/T)*inttrap(t,f)
na=5 //n° de armónicos
s=0
for n=1:na
a(n)=(2/T)*inttrap(t,f.*cos(n*w*t))
b(n)=(2/T)*inttrap(t,f.*sin(n*w*t))
s=s+a(n)*cos(n*w*t)+b(n)*sin(n*w*t)
end
ff=a0+s

disp('a0')
disp(a0)
disp('Valores de an')
for n=1:na
mprintf('\n a(%d)= %1.5f',n,a(n))
end
disp('\n')
disp('Valores de bn')
for n=1:na
mprintf('\n b(%d)= %1.5f',n,b(n))
end
// Gráficas de f y ff

plot(t,f,t,ff,'linewidth',3)
xlabel('t','fontsize',4)

18
title('Gráficas de f y ff','fontsize',4)
legend('Señal','Fourier','fontsize',4)
xgrid

//Serie de Fourier
mprintf('\n')
disp('Serie de Fourier')
mprintf('\n')
mprintf('\n ff(t) = %1.5f',a0)
for n=1:na
mprintf('\n %1.5fcos(%dwt)%+1.5fsen(%dwt)',a(n),n,b(n),n)
end

2) Se tiene una señal periódica , cuya regla de correspondencia para un

periodo es

a) Encontrar los coeficientes de la serie de Fourier para n = 1, 2, ...,10

19
 b) Graficar la función y la serie de Fourier con 5, 8 y 10 armónicos

5 armónicos:

8 armónicos:

10 armónicos:

20
 c) Expresar la serie de Fourier en términos de t con 5 armónicos

Código:

//Laboratorio 14
//Transformada de Laplace
//Navarro Quispe,Alexander
//Aplicación 2
clc
clear
close
//***************************
ti=0
tf=6
T=tf-ti //periodo de la señal
w=2*%pi/T //Frecuencia angular
h=0.00001
t=ti:h:tf
f1=t
f2=2
f3=-2*t+12
t1=(t>=0)&(t<=2)
t2=(t>2)&(t<=5)
t3=(t>5)&(t<=6)

21
f=f1.*t1+f2.*t2+f3.*t3 //señal
//Coeficientes de Fourier
a0=(1/T)*inttrap(t,f)
na=5 //n° de armónicos
s=0
for n=1:na
a(n)=(2/T)*inttrap(t,f.*cos(n*w*t))
b(n)=(2/T)*inttrap(t,f.*sin(n*w*t))
s=s+a(n)*cos(n*w*t)+b(n)*sin(n*w*t)
end
ff=a0+s

disp('a0')
disp(a0)
disp('Valores de an')
for n=1:na
mprintf('\n a(%d)= %1.5f',n,a(n))
end
disp('\n')
disp('Valores de bn')
for n=1:na
mprintf('\n b(%d)= %1.5f',n,b(n))
end
// Gráficas de f y ff

plot(t,f,t,ff,'linewidth',3)
xlabel('t','fontsize',4)

22
title('Gráficas de f y ff','fontsize',4)
legend('Señal','Fourier','fontsize',4)
xgrid

//Serie de Fourier
mprintf('\n')
disp('Serie de Fourier')
mprintf('\n')
mprintf('\n ff(t) = %1.5f',a0)
for n=1:na
mprintf('\n %1.5fcos(%dwt)%+1.5fsen(%dwt)',a(n),n,b(n),n)
end

3) Se tiene una señal periódica , cuya regla de correspondencia para un

periodo es

a) Encontrar los coeficientes de la serie de Fourier para n = 1, 2, ...,10

23
 b) Graficar la función y la serie de Fourier con 5, 8 y 10 armónicos

5 armónicos:

8 armónicos:

10 armónicos:

24
 c) Expresar la serie de Fourier en términos de t con 5 armónicos

Código:

//Laboratorio 14
//Transformada de Laplace
//Parco Guerra, Alonzo Fernando
//Aplicación 3
clc
clear
close
//***************************
ti=0
tf=8
T=tf-ti //periodo de la señal
w=2*%pi/T //Frecuencia angular
h=0.00001
t=ti:h:tf
f1=1.5*t
f2=-1.5*t+6
f3=1.5*t-12
t1=(t>=0)&(t<2)
t2=(t>=2)&(t<6)
25
t3=(t>=6)&(t<=8)
f=f1.*t1+f2.*t2+f3.*t3 //señal
//Coeficientes de Fourier
a0=(1/T)*inttrap(t,f)
na=5 //n° de armónicos
s=0
for n=1:na
a(n)=(2/T)*inttrap(t,f.*cos(n*w*t))
b(n)=(2/T)*inttrap(t,f.*sin(n*w*t))
s=s+a(n)*cos(n*w*t)+b(n)*sin(n*w*t)
end
ff=a0+s

disp('a0')
disp(a0)
disp('Valores de an')
for n=1:na
mprintf('\n a(%d)= %1.5f',n,a(n))
end
disp('\n')
disp('Valores de bn')
for n=1:na
mprintf('\n b(%d)= %1.5f',n,b(n))
end
// Gráficas de f y ff

plot(t,f,t,ff,'linewidth',3)

26
xlabel('t','fontsize',4)
title('Gráficas de f y ff','fontsize',4)
legend('Señal','Fourier','fontsize',4)
xgrid

//Serie de Fourier
mprintf('\n')
disp('Serie de Fourier')
mprintf('\n')
mprintf('\n ff(t) = %1.5f',a0)
for n=1:na
mprintf('\n %1.5fcos(%dwt)%+1.5fsen(%dwt)',a(n),n,b(n),n)
end

Observaciones

● Lo aprendido en clase puede aplicarse en distintos rubros que requieren de

cálculos como los trabajados.
● Se requiere de lo aprendido con anterioridad para poder obtener
respuestas adecuadas a los ejercicios propuestos.
● La función de transferencia ha de realizarse con mucho cuidado, para así
poder tener un resultado correcto.
● Se realizó el ajuste de las gráficas luego de ejecutar el programa, para así
poder analizar mejor la imagen.

Conclusiones

● Se analizó la función de transferencia, de esta manera se realizó las

gráficas para un mejor entendimiento.
● Se comprendió el proceso de solución de las transformadas de Laplace.

27
 ● Se utilizaron los comandos de Scilab para obtener y calcular los resultados
de los ejercicios propuestos.
● Se encontró la transformada de Laplace de una ecuación diferencial en el
programa Scilab con los comandos aplicados.

Recomendaciones

● Identificar las funciones del Scilab favorece el desarrollo del laboratorio

realizado.
● Determinar la ecuación diferencial favorece la resolución de los ejercicios
propuestos.
● Revisar los temas anteriores, así como sus funciones propuestas, para
aplicar lo aprendido en el laboratorio actual junto a lo anterior. Así poder
manejar el Scilab de la mejor manera posible.
● Presentar una adecuada explicación de cada variable o función aplicada,
facilita la comprensión y análisis del programa realizado. Así como la
búsqueda y reconocimiento de errores.

28