UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE ELECTROTECNIA Y COMPUTACION

Ingeniería Electrónica

5Guía de laboratorio No.: 6

Asignatura: Sistemas y Señales

1. Tema
Series de Fourier

2. Objetivos

3. Contenido

4. Medios a utilizar

Computadora
Simulador MATLAB

5. Introducción

1
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

3.1 Trabajo

A travez de las series de fourier podemos representar una señal como una superposicion pondera-
da de senoides complejas.
Las ecuaciones que nos permiten desarrollar esta representacion son:

∞
x (t )=x ( t+T )= ∑ ak e j w o kt
k=−∞

❑
1
a k= ∫ x (t )e− j w o kt dt
T T

Tomemos la siguiente señal cuadrada para representarla a travez de series de fourier

x(t)
1

0.8

0.6

0.4

0.2

0
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
tiempo t(s)

Calculemos los coeficientes : 0.5 2π

1 −j kt
a k= ∫ x (t )e 1
dt
1 −j
2π
0t 1 −0.5
a 0= ∫ x (t )e 1
dt
T T 0.25

0.25
a k= ∫ − j 2 π kt
e dt
1 −0.25
a 0= ∫ 1 dt
1 −0.25
a k ¿−
1 − j 2 π kt 0.25
jk 2 π
e
−025 |
0.25
a 0 ¿ t|
−0.25 1 2π j
2π
k
a k ¿− (e ¿ ¿− j k −e 4
)¿
jk 2 π 4
a 0=0.5 π
1 π j k
a k ¿− (e ¿ ¿− j k −e 2 )¿
jk 2 π 2

1
a k= ¿)
πk

1 π
a k= sen ( k )
πk 2

2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

Entonces la solucion de la representacion en series de fourier para la onda cuadrada anterior es:
∞
x (t )=a 0+ ∑ a k e j w o kt dt
k=−∞
∞ j 2 π kt
1 π
x (t )=0.5+ ∑ sen( k)e
k=−∞ π k 2
Ahora desarrollaremos un algoritmo en matlab para poder implementar la ecuación anterior:
Fs=10000; %Frecuencia de muestreo
t=-2:1/Fs:2; %Rango de tiempo donde se va a graficar la señal
xt=0.5; %Cargar en xt con el valor de a0
N=3; %Usaremos 3 coeficientes ak
for k = [Link]N %Calculo de la serie de fourier Desde 1 hasta N
akpos=(1/(k*pi))*(sin((k*pi)/2)); %Coeficientes ak positivos
k=-k; %Cambio de k a k negativo para calcular ak negativos
akneg=(1/(k*pi))*(sin((k*pi)/2)); %Coeficientes ak negativos
x=(akpos+akneg)*exp(1i*k*2*pi*t); %carga en x el valor de ak evaluada en N
xt=xt+x; %Almacena en xt el valor x y lo suma con el valor anterior de xt
x=0; %Limpia a x para volver a calcular x en el nuevo N
end
plot(t,xt)
TITLE('Serie de Fourier en TC de x(t)')
xlabel('Tiempo')

Grafica de xt con N=1 Observación

Grafica de xt con N=2 Observación

3
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

Grafica de xt con N=3 Observación

Grafica de xt con N=10 Observación

Grafica de xt con N=20 Observación

4
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

Grafica de xt con N=100 Observación

Ahora cambie el algoritmo anterior para representar la siguiente señal

x(t)
1

0.8

0.6

0.4

0.2

0
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
tiempo t(s)

∞
x (t )=a 0+ ∑ a k e j w o kt dt
k=−∞
∞
x (t )=0.5+ ∑ ¿ ¿
k=−∞

Grafica de xt con N=1 Observación

5
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

Grafica de xt con N=2 Observación

Grafica de xt con N=3 Observación

Grafica de xt con N=10 Observación

Grafica de xt con N=20 Observación

6
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

Grafica de xt con N=100 Observación

Trabajo Independiente:
Determine la representación mediante la serie de Fourier para la onda triangular descrita a continuación:
x(t)
1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2

7
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica

Algoritmo de xt con N=100 Grafica de xt con N=100