Laboratorio N° 1: Trabajo en Simulink - Serie de Fourier 30/08/2018 1

UNIVERSIDAD SURCOLOMBIANA
INGENIERIA ELECTRÓNICA
ELECTRÓNICA INDUSTRIAL
INFORME N° 1
TRABAJO EN SIMULINK - SERIE DE FOURIER
JUAN SEBASTIAN CASTILLO TOVAR
CODIGO 20152140522
ALVARO ANDRES QUINTERO DIAZ
CODIGO 20151134150
BRAYAN BRAVO FIERRO
CODIGO 20151135393
GRUPO 02

 red corrientes periódicas no senoidales. Estas corrientes están

Resumen—En esta práctica se analizarán cinco señales de formadas por un componente fundamental de frecuencia 50 ó
diferentes formas de onda, a las cuales se les encontrarán los 60 Hz, más una serie de corrientes superpuestas de frecuencias,
diferentes coeficientes de Fourier, dependiendo si la función es par múltiplos de la fundamental, que denominamos Armónicos
o impar, con el fin de encontrar los diferentes armónicos de dichas eléctricos, que generan costes técnicos y económicos
señales y poder graficarlas con ayuda del software Simulink.
importantes.
Simulink: Es una toolbox especial de Matlab que sirve para
Palabras claves—Serie de Fourier, Armónico, Función par,
Función impar, Simulink. simular el comportamiento de los sistemas dinámicos. Puede
simular sistemas lineales y no lineales, modelos en tiempo
continuo y tiempo discreto y sistemas híbridos de todos los
I. OBJETIVOS anteriores. Es un entorno gráfico en el cual el modelo a simular
se construye clicando y arrastrando los diferentes bloques que
 Encontrar los diferentes coeficientes de Fourier de cada lo constituyen.
señal, para poder simularla en simulink. Funcion par: Una función es par si, para cada x en
 Encontrar la serie de Fourier de las diferentes señales el dominio de f , f (– x ) = f (x). Las funciones pares tienen
dadas. simetría reflectiva a través del eje de las y.
 Ganar destreza encontrando cuantos y que armónicos tiene Funcion impar: Un función es impar si, para cada x en el
una señal. dominio de f, f (– x) = – f (x). Las funciones impares tienen
 Identificar cuando una señal es par o impar. simetría rotacional de 180º con respecto del origen.
 Adquirir conocimiento y experiencia en el uso de simulink
como software de simulación. III. ELEMENTOS MATERIALES Y EQUIPOS

 Simulink
II. MARCO TEORICO
Series de Fourier: Es una serie infinita que converge PROCEDIMIENTO
puntualmente a una función periódica y continua a trozos (o
por partes). Las series de Fourier constituyen la herramienta Para la señal 1:
matemática básica del análisis de Fourier empleado para Se realizarán los cálculos de los coeficientes de la señal 1:
analizar funciones periódicas a través de la descomposición de
dicha función en una suma infinita de funciones sinusoidales 𝑎0 = 0
mucho más simples (como combinación de senos y cosenos con Para coeficientes a=
frecuencias enteras).
Armónicos: Las cargas no lineales tales como: rectificadores, −4 2𝜋 ∗ 𝑛
𝑎𝑛 = ∗ (𝑠𝑒𝑛 ( ))
inversores, variadores de velocidad, hornos, etc; absorben de la 𝜋∗𝑛 3
Laboratorio N° 1: Trabajo en Simulink - Serie de Fourier 30/08/2018 2

Reemplazando:
n=1
−4 2𝜋 ∗ 1
𝑎1 = ∗ (𝑠𝑒𝑛 ( )) = −1.102657791
𝜋∗1 3

n=3
−4 2𝜋 ∗ 3
𝑎3 = ∗ (𝑠𝑒𝑛 ( )) = 0
𝜋∗3 3

n=5
−4 2𝜋 ∗ 5
𝑎5 = ∗ (𝑠𝑒𝑛 ( )) = 0.2205315582
𝜋∗5 3

Para coeficientes b = Figura 2. Grafica de la señal 1 simulada en simulink

4 2𝜋 ∗ 𝑛
𝑏𝑛 = ∗ (1 + cos⁡( )
𝑛∗𝜋 3
n=1
4 2𝜋 ∗ 1
𝑏1 = ∗ (1 + cos ( ) = 0.6366197724
1∗𝜋 3

n=3

4 2𝜋 ∗ 3
𝑏3 = ∗ (1 + cos ( ) = 0.8488263632
3∗𝜋 3
n=5

4 2𝜋 ∗ 5 Para la señal 2:
𝑏5 = ∗ (1 + cos ( ) = 0.1273239545 𝜋
5∗𝜋 3 1 2 1
𝑎0 = ∫ cos(𝜃) 𝑑𝜃 =
𝑇 0 3𝜋
La función quedaría:
∞ ∞
𝜋
𝑓(𝑥) = 0 + ∑ 𝑎𝑛(cos⁡(𝜃))𝑑𝜃 + ∑ 𝑏𝑛(𝑠𝑒𝑛(𝜃))𝑑𝜃 2 2
𝑎𝑛 = ∫ cos(𝜃) ∗ cos(𝑛𝜃)
𝑛=1 𝑛=1 3𝜋 0

2𝜋 Reemplazando: 𝜃 = 𝜋/2
𝑓(𝑥) = 0 + ((−1.102657791 ∗ cos⁡( ) + (0 ∗ cos(𝜋))
3 n=1
10 ∗ 𝜋 2 1 𝑠𝑒𝑛(2𝜃)
+ (0.2205315582 ∗ cos ( )) 𝑎1 = ∗ ( (𝜃 + )) = 0.3333333
3 3𝜋 2 2
2𝜋
+ ((0.6366197724 + 𝑠𝑒𝑛 ( ))
3 n=2
+ (0.8488263632 ∗ 𝑠𝑒𝑛(2𝜋))
10𝜋 1 1 𝑠𝑒𝑛(3𝜃)
+ (0.1273239545 ∗ 𝑠𝑒𝑛 ( )) 𝑎2 = ∗ (( ( )) + 𝑠𝑒𝑛(𝜃)) = 0.1414710
3 3𝜋 3 2
Luego de calcular los coeficientes de la señal procedemos a
simular en Simulink la señal con sus diferentes armónicos. n=3
1
Figura 1. Circuito de la señal 1 simulado en simulink 𝑎3 = ∗ (0 + 0) = 0
3𝜋

n=4
2 𝑠𝑒𝑛(5𝜃) 𝑠𝑒𝑛(3𝜃)
𝑎4 = ∗ (( )+( )) = −0.0282942
3𝜋 5 3
n=5
1
𝑎5 = ∗ (0 + 0) = 0
3𝜋

La función quedaría:
Laboratorio N° 1: Trabajo en Simulink - Serie de Fourier 30/08/2018 3

∞ ∞
𝑛=1
𝑓(𝑥) = 0 + ∑ 𝑎𝑛(cos⁡(𝜃))𝑑𝜃 + ∑ 𝑏𝑛(𝑠𝑒𝑛(𝜃))𝑑𝜃 −2 cos(1𝜋)
𝑛=1 𝑛=1
𝑏1 = =2
1
1 𝑛=2
𝑓(𝑥) = + 0.3333333 ∗ cos⁡(𝜋) + 0.1414710 ∗ cos(𝜋)
3𝜋 −2 cos(2𝜋)
3𝜋 𝑏2 = = −1
+ 0 − 0.0282942cos⁡( ) + 0 2
2
Luego de calcular los coeficientes de la señal procedemos a 𝑛=3
simular en Simulink la señal con sus diferentes armónicos. −2 cos(3𝜋) 2
𝑏3 = =
3 3

𝑛=4
−2 cos(4𝜋) 1
𝑏4 = =−
4 2

Figura 3. Circuito de la señal 2 simulado en simulink 𝑛=5

−2 cos(5𝜋) 2
𝑏5 = =
5 5
∞ ∞

𝑓(𝑥) = 𝑎0 + ∑ 𝑎𝑛(cos⁡(𝜃))𝑑𝜃 + ∑ 𝑏𝑛(𝑠𝑒𝑛(𝜃))𝑑𝜃

𝑛=1 𝑛=1

1 2 2 3
𝑓(𝑥) = 0 + 2 (sen⁡( )) − 1 (sen ( )) + (sen ( ))
2𝜋 2𝜋 3 2𝜋
1 4 2 5
− (sen ( )) + (sen ( ))
2 2𝜋 5 2𝜋

Luego de calcular los coeficientes de la señal procedemos a

simular en Simulink la señal con sus diferentes armónicos.
Figura 4. Grafica de la señal 2 simulada en simulink
Figura 5. Circuito de la señal 3 simulado en simulink

Para la señal 3:
Se realizarán los cálculos de los coeficientes de la señal 3:
Como el area de la porcion positiva es igual al area de la porcion
negative el componente DC sera igual a 0.
En este caso tomamos un periodo 2π
Como es una función impar contiene solo coeficientes de tipo
sen (𝜃) 𝑎𝑛 por lo general es 0; para todo n
𝑎0 = 0
𝑎𝑛 = 0
2 𝑇 2 𝜋
𝑏𝑛 = ∫ 𝑥(𝑡) sin(𝑛𝑡)𝑑𝑡 = ∫ 𝑥(𝑡) sin(𝑛𝑡)𝑑𝑡 × 2
𝑇 0 2𝜋 0
Figura 6. Grafica de la señal 3 simulada en simulink
−2 cos(𝑛𝜋)
𝑏𝑛 =
𝑛

Para sacar los primeros cinco armónicos:

Laboratorio N° 1: Trabajo en Simulink - Serie de Fourier 30/08/2018 4

Para la señal 4:
Se realizarán los cálculos de los coeficientes de la señal 4:
En este caso tomamos un periodo 2π.
Como es una función par el coeficiente de Fourier bn por lo
general es 0
𝑇
1 2
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥 = 3 × 2 = 6
𝑇 0
2 𝑇 12 𝜋
𝑎𝑛 = ∫ 𝑓(𝑥) cos(𝑛𝜃) = sin(𝑛 )
𝑇 0 𝑛𝜋 2
𝑏𝑛 = 0
Para sacar los primeros cinco armónicos
𝑛=1
12
𝑎1 = = 3.8197
𝜋
𝑛=2 Figura 8. Grafica de la señal 4 simulada en simulink
6
𝑎2 = sin(𝜋) = 0
𝜋
𝑛=3
4 3𝜋 4
𝑎3 = sin ( ) = − = −1.2732
𝜋 2 𝜋
𝑛=4
3
𝑎4 = sin(2𝜋) = 0
𝜋
𝑛=5
12 5𝜋 12
𝑎5 = sin ( ) = = 0.7639
5𝜋 2 5𝜋
∞ ∞

𝑓(𝑥) = 𝑎0 + ∑ 𝑎𝑛(cos⁡(𝜃))𝑑𝜃 + ∑ 𝑏𝑛(𝑠𝑒𝑛(𝜃))𝑑𝜃

𝑛=1 𝑛=1
1 3 Para la señal 5:
𝑓(𝑥) = 6 + 3.8197 (cos⁡( )) − 1.2732 (cos ( ))
2𝜋 2𝜋
5 Para este caso tenemos una señal impar
+ 0.7639 (cos⁡( )) 𝑇
2𝜋 1 2 1
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥 =
𝑇 0 𝜋
2 𝑇 1 1 + (−1𝑛 )
𝑎𝑛 = ∫ 𝑓(𝑥) cos(𝑛𝜃)𝑑𝜃 = ( 2 )
𝑇 0 𝜋 𝑛 −1
Luego de calcular los coeficientes de la señal procedemos a
simular en Simulink la señal con sus diferentes armónicos.
2 𝑇
𝑏𝑛 = ∫ 𝑓(𝑥) sin(𝑛𝜃) 𝑑𝜃 = 0⁡𝑝𝑎𝑟𝑎⁡𝑡𝑜𝑑𝑜𝑠⁡𝑚𝑒𝑛𝑜𝑠⁡𝑝𝑎𝑟𝑎⁡𝑏1
𝑇 0
Figura 7. Circuito de la señal 4 simulado en simulink
Al evaluar la integral para hallar el bn nos encontramos con un
caso en el cual se puede lograr una indeterminación y es el
caso cuando n=1 entonces se puede decir que el b1 va a ser
diferente de 0, por lo tanto evaluamos de nuevo la integra
teniendo en cuenta n=1. Y obtenemos:
Laboratorio N° 1: Trabajo en Simulink - Serie de Fourier 30/08/2018 5

1
𝑏1 =
2

Sacamos los armónicos de con el an

𝑛=1

1 1 + (−1𝑛 )
𝑎1 = ( 2 )=0
𝜋 𝑛 −1

𝑛=2

1 1 + (−1𝑛 ) 2
𝑎2 = ( 2 )= = 0.212206
𝜋 𝑛 −1 3𝜋
IV. ANALISIS DE RESULTADOS
𝑛=3 Como se pudo observar al representar cada función a traves del
matlab se obtenian con cierta curvatura en sus extremos, esto se
1 1 + (−1𝑛 ) debe a que el simulink solo representa funciones seno y coseno
𝑎3 = ( 2 )=0 dependiendo la fase de cada una de las graficas. Tambien se
𝜋 𝑛 −1
puede apreciar como diferentes propiedades de las graficas,
como determinar si son pares o no; ayudaban en la
𝑛=4
simplificacion a la hora de encontrar armonicos
1 1 + (−1𝑛 ) 2
𝑎4 = ( 2 )= = 0.0424
𝜋 𝑛 −1 15𝜋 V. CONCLUSIONES
 Simulink permite un analisis detallado de la funcion para
𝑛=5 determinar que tan senoidal es y que tanta distorsión hay
en ella ósea que cantidad de armónicos posee.
1 1 + (−1𝑛 )  La señal 1 presenta una mayor taza de distorcion armonica
𝑎5 = ( 2 )=0
𝜋 𝑛 −1 (THD) con respecto a la señal 2 debido a que contiene
∞ ∞
mayor presencia de armónicos.
𝑓(𝑥) = 𝑎0 + ∑ 𝑎𝑛(cos⁡(𝜃))𝑑𝜃 + ∑ 𝑏𝑛(𝑠𝑒𝑛(𝜃))𝑑𝜃 
𝑛=1 𝑛=1
1 1 1 1 VI. REFERENCIAS
𝑓(𝑥) = + (sin⁡( )) + 0.212206 (cos ( ))
𝜋 2 2𝜋 𝜋 [3] disponible en la web: http://circutor.es/es/formacion/armonicos-
2
+ 0.0424 (cos⁡( )) electricos; consultada el 30/08/2018
𝜋
disponible en la web: https://es.wikipedia.org/wiki/Serie_de_Fourier;
Luego de calcular los coeficientes de la señal procedemos a consultada el 30/08/2018
simular en Simulink la señal con sus diferentes armónicos.
disponible en la web:
https://ocw.upc.edu/sites/all/modules/ocw/estadistiques/download.ph
Figura 9. Circuito de la señal 5 simulado en simulink p?file=51427/2011/1/54513/tema_5_simulink-5156.pdf; consultada
el 30/08/2018

disponible en la web
https://www.varsitytutors.com/hotmath/hotmath_help/spanish/topics/
even-odd-functions; consultada el 30/08/2018

Figura 10. Grafica de la señal 5 simulada en simulink