Practica 4.
Transformada de Fourier
Preámbulo
Limpia el workspace
clear
close all;
Declaración de variables simbólicas
syms t
1. Función pulso triangular
Obtenga la transformada de Fourier y luego grafique el espectro de frecuencia y de fase de un
pulso triangular
1.1 Declaración de la función
tao = 1/2;
f_1 = rectangularPulse(-1/2,1/2,t)*(1-abs(t)/tao)
1.2 Grafica de la función
disp(['Gráfica de pulso triangular con Tao =' ' ' num2str(tao)])
function_graph(f_1, 'f', [-1 1])
1.3 Transformada de Fourier
disp('Transformada de Fourier')
ft_1 = fourier(f_1)
1.4 Espectro de amplitud
amp_graph(ft_1, 'F', [-35 35])
1.5 Espectro de fase
angle_graph(ft_1, 'F')
2. Función pulso triangular
Obtenga la transformada de Fourier y luego grafique el espectro de frecuencia y de fase de un
pulso RF, cuando fc esta en la banda de radio frecuencia. Asigne valores a las constantes: A,
τ, y fc
2.1 Declaración de la función
A = 10;
tao = 1;
f_c = 10*(1/tao);
disp('Pulso radio frecuencia z(t)')
z = rectangularPulse(-tao/2,tao/2,t)*cos(2*pi*f_c*t)*A
2.2 Grafica de la función
function_graph(z, 'z', [-tao tao])
2.3 Transformada de Fourier
disp('Transformada de Fourier')
z_ft = fourier(z)
2.4 Espectro de amplitud
amp_graph(z_ft, 'Z')
2.5 Espectro de fase
angle_graph(z_ft, 'Z')
3. Función w(t)
Obtenga la transformada de Fourier y luego grafique el espectro de frecuencia y de fase de la
v(t).
3.1 Declaración de la función
syms t
A = 1;
B = 2;
T = 1;
y_1 = A*triangularPulse(-T+(4*T), 4*T, T+(4*T), t) ;
w_2 = B*triangularPulse(-T-(4*T), -4*T, T-(4*T), t) ;
w = y_1 + w_2
3.2 Grafica de la función
function_graph(w, 'w', [-6*T 6*T])
3.3 Transformada de Fourier
disp('Transformada de Fourier')
w_ft = fourier(w)
3.4 Espectro de amplitud
amp_graph(w_ft, 'W', [-10 10])
3.5 Espectro de fase
angle_graph(w_ft, 'W')
4. Función y(t)
Obtenga la transformada de Fourier y luego grafique el espectro de frecuencia y de fase de la
señal y(t). Asigne valores a las constantes.
4.1 Declaración de la función
syms t
A = 1;
B = 2;
T = 1;
y_1 = A*rectangularPulse(-T+(4*T), T+(4*T), t) ;
y_2 = B*rectangularPulse(-T-(4*T), T-(4*T), t) ;
y = y_1 + y_2;
4.2 Grafica de la función
function_graph(y,'y', [-6*T 6*T])
4.3 Transformada de Fourier
disp('Transformada de Fourier')
y_ft = fourier(y)
4.4 Espectro de amplitud
amp_graph(y_ft, 'Y')
4.5 Espectro de fase
angle_graph(y_ft, 'Y')
5. s(t)
Obtenga la transformada de Fourier y luego grafique el espectro de frecuencia y de fase de la
señal y(t). Asigne valores a las constantes.
5.1 Declaración de la función
syms t
s = cos(20*pi*t)*cos(80*pi*t);
5.2 Grafica de la función
function_graph(s,'s', [-0.2 0.2])
5.3 Transformada de Fourier
disp('Transformada de Fourier')
s_ft = fourier(s)
5.4 Espectro de amplitud
amp_graph(s_ft, 'S')
5.5 Espectro de fase
angle_graph(s_ft, 'S')
function amp_graph(func, name, xlimits)
disp('Espectro de Amplitud')
if ~exist('xlimits','var')
fplot(abs(func), '-r');
else
fplot(abs(func),xlimits, '-r');
end
grid on;
xlabel('\omega');
name_str = strcat('|', name, '(\omega)|');
ylabel(name_str);
title(name_str);
end
function angle_graph(func, name)
disp('Espectro de Fase')
fplot(angle(func));
grid on;
xlabel('\omega');
label = ['\phi' ' ' name ' ' '(\omega)'];
ylabel(label);
title(strcat('Angle \phi', {' '}, name, '(\omega)'));
end
function function_graph(func, name, xlimits)
disp(['Grafica de función ' name '(t)'])
fplot(func, xlimits)
grid on
title([name '(t)'])
xlabel('t')
ylabel([name '(t)'])
end