Fourier Series Example Script

% This script plots a Fourier Series for a Square Wave

clf; t=0:.01:10; T=2.5 M=50 sum1=0; for m=1:2:M,

   sum1 = sum1+4/m/pi*sin(m*pi/2)*cos(2*pi*m*t/T);

end plot(t,sum1,'b-',t(1:10:end),sum1(1:10:end),'r*') title('Fourier Series Representation of a Square Wave') xlabel('time (seconds)') ylabel('Function') grid on; axis([0,10,-2,2]) legend('Five Terms','Five Terms Sampled') print("squarewave.png","-dpng")  % Prints the plot to a png file called squarewave.png