UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Guía de laboratorio No.: 6
Asignatura: Sistemas y Señales
1. Tema
Series de Fourier y Filtro RC
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
Si la entrada de un sistema LTI se expresa como una suma ponderada de senoides entonces la salida también es una
suma ponderada de senoides. El k-esimo peso de la suma de la salida está dada por el producto del k-esimo peso en la
suma de entrada y la respuesta en frecuencia del sistema evaluada en la frecuencia de la k-esima senoide. Por tanto si:
∞
x (t )= ∑ a k e j w o kt
k=−∞
Entonces la salida y(t) es:
∞
y ( t ) = ∑ H ( jkwo)a k e j wo kt
k=−∞
Donde H(jw) es la respuesta en frecuencia del sistema
Si consideramos la siguiente señal de entrada:
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)
Y la aplicamos al siguiente circuito:
Convirtiendo a su notación fasorial Reacomodando la ecuación
1 vo 1
ZC = Z R=R =
jwc v i jwRC +1
Haciendo un divisor de voltaje Que es la respuesta en frecuencia
1
zc H ( jw)=
v o= vi jwRC+1
zc+z R
1
vo jwc
=
vi 1
+R
jwc
2
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Ahora usando:
∞
y ( t ) = ∑ H ( jwk )ak e j w k t
k=−∞
1
RC
H ( jkwo)=
1
jk w o+
RC
1 π
a k= sen ( k )
πk 2
Y la salida del sistema seria:
( )
1
( π1k sen ( π2 k )) (e
∞
RC
y ( t ) =a 0+ ∑ 1
j 2 π kt
)
k=−∞
jk w o+
RCen matlab para poder implementar la ecuación anterior:
Ahora desarrollaremos un algoritmo
Fs=10000; %Frecuencia de muestreo
t=-2:1/Fs:2;%Rango de tiempo donde se va a graficar la señal
Fc=1000;%Frecuencia de corte del fitro
RC=1/Fc;%Valor de RC
N=100;%Usaremos 100 coeficientes ak
y=0.5;%Cargar en y con el valor de a0
for k = [Link]N%Calculo de la serie de fourier Desde 1 hasta N
y=(1/(1i*2*pi*k*RC+1))*(1/(k*pi))*(sin((k*pi)/2))*exp(1i*k*2*pi*t)+y;%ak pos
k=-k;%Cambio de k a k negativo para calcular ak negativos
y=(1/(1i*2*pi*k*RC+1))*(1/(k*pi))*(sin((k*pi)/2))*exp(1i*k*2*pi*t)+y;%ak neg
end
plot(t,y)
Grafica de y con Fc=1000 Observación
3
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Grafica de y con Fc=500 Observación
Grafica de y con Fc=100 Observación
Grafica de y con Fc=50 Observación
4
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Grafica de y con Fc=30 Observación
Grafica de y con Fc=10 Observación
Grafica de y con Fc=5 Observación
5
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Grafica de y con Fc=1 Observación
Repita el proceso anterior para una señal diente de sierra.
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 y con Fc=1000 Observación
6
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Grafica de y con Fc=500 Observación
Grafica de y con Fc=100 Observación
Grafica de y con Fc=50 Observación
Grafica de y con Fc=30 Observación
7
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE ELECTROTECNIA Y COMPUTACION
Ingeniería Electrónica
Grafica de y con Fc=10 Observación
Grafica de y con Fc=5 Observación
Grafica de y con Fc=1 Observación
