Arquivos para OCTAVE

function[azero,a,b,xsin,xcos, w,theta]=fourier(N,m,time,x,t)

pi = 3.1416;

sumz = 0.0;

w = (2*pi)/0.12;

for i=1:N

  sumz=sumz+x(i);

endfor

azero = (2.0/N)*sumz;

for ii=1:m

  sums = 0.0;

  sumc = 0.0;

  for i=1:N

    theta(ii) = (2.0*pi*t(i)*ii)/time;

    xcos(i)=x(i)*cos(theta(ii));

    xsin(i)=x(i)*sin(theta(ii));

    sums=sums+xsin(i);

    sumc=sumc+xcos(i);

endfor

a(ii)=2.0*sumc/N;

b(ii)=2.0*sums/N;

endfor