Fourier Series Example Script: Difference between revisions
Jump to navigation
Jump to search
No edit summary |
No edit summary |
||
| Line 1: | Line 1: | ||
% This script plots a Fourier Series for a Square Wave | % This script plots a Fourier Series for a Square Wave | ||
clf; | clf; | ||
t=0:.01:10; | t=0:.01:10; | ||
T=2.5 | T=2.5 | ||
M=50 | M=50 | ||
sum1=0; | sum1=0; | ||
for m=1:2:M, | for m=1:2:M, | ||
sum1 = sum1+4/m/pi*sin(m*pi/2)*cos(2*pi*m*t/T); | |||
end | end | ||
plot(t,sum1,'b-',t(1:10:end),sum1(1:10:end),'r*') | plot(t,sum1,'b-',t(1:10:end),sum1(1:10:end),'r*') | ||
title('Fourier Series Representation of a Square Wave') | title('Fourier Series Representation of a Square Wave') | ||
xlabel('time (seconds)') | xlabel('time (seconds)') | ||
ylabel('Function') | ylabel('Function') | ||
grid on; | grid on; | ||
axis([0,10,-2,2]) | axis([0,10,-2,2]) | ||
legend('Five Terms','Five Terms Sampled') | legend('Five Terms','Five Terms Sampled') | ||
print("squarewave.png","-dpng") % Prints the plot to a png file called squarewave.png | print("squarewave.png","-dpng") % Prints the plot to a png file called squarewave.png | ||
Revision as of 14:24, 29 September 2010
% 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