% LABORATORIO DE PROCESAMIENTO DIGITAL DE BIOSEÑALES E IMÁGENES
% PRÁCTICA 4
% EXPERIMENTO 2
close all
syms w t
syms F(w)
F(w)=int(exp(-sqrt(-1)*w*t),t,-1,1)
Amp(w)=abs(F(w))
q = subs(Amp,[-6*pi:0.1:6*pi])
figure(1)
plot([-6*pi:0.1:6*pi],q)
title(['Espectro de Amplitud'])
xlabel(['Frecuencia [rad/s]'])
ylabel(['Amplitud'])
grid on
hold on
fase = angle(F(w))
r = subs(fase,[-6*pi:0.1:6*pi])
figure(2)
plot([-6*pi:0.1:6*pi],r,'g')
grid on
title(['Espectro de fase'])
xlabel(['Frecuencia [rad/s]'])
ylabel(['Fase'])