FOURIER COEFICIENTES 

clc

clear all

syms t

syms k

fx=t;

xt=0;

xk=int((fx*exp(-j*pi*k*t)),-1,1)/2;

disp('la función para los coeficientes de fourier es:')

disp(xk);

N=5;

for co=-N:1:N

   if co==0

        v=0 

    else

      xi=xk*exp(j*pi*k*t);

      v=subs(xi,k,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,'--b.')

      grid on

      hold on

      xlabel('Tiempo')

      ylabel('x(t) aprox')

      title('Serie de fourier F(x)')

end

N=5

N=20

N=50