El ejercicio se basa en la siguiente función, la cual modificaremos bajo el concepto de la sumatoria de fourier para los valores de N=5, N=20 y N=50.
El codigo que se uso en Matlab para graficar las funciones fue el siguiente:
clc
clear all
syms t
syms h
fx=t;
xt=0;
xh=int((fx*exp(-j*pi*h*t)),-1,1)/2;
disp('la función para los coeficientes de fourier es:')
disp(xh);
N=5;
clear all
syms t
syms h
fx=t;
xt=0;
xh=int((fx*exp(-j*pi*h*t)),-1,1)/2;
disp('la función para los coeficientes de fourier es:')
disp(xh);
N=5;
for co=-N:1:N
if co==0
v=0
else
xi=xh*exp(j*pi*h*t);
v=subs(xi,h,co);
end
disp('componente numero:')
disp(co)
disp(v)
xt=xt+v;
end
if co==0
v=0
else
xi=xh*exp(j*pi*h*t);
v=subs(xi,h,co);
end
disp('componente numero:')
disp(co)
disp(v)
xt=xt+v;
end
for ti=-2:0.01:2
xti=subs(xt,t,ti);
plot(ti,xti,'--g.')
grid on
hold on
xlabel('Tiempo')
ylabel('X(t) Aproximada')
title('Serie de Fourier F(X)')
end
xti=subs(xt,t,ti);
plot(ti,xti,'--g.')
grid on
hold on
xlabel('Tiempo')
ylabel('X(t) Aproximada')
title('Serie de Fourier F(X)')
end
No hay comentarios:
Publicar un comentario