Andrew Roth's Fourier Series/Laplace Transform Project: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
Created page with '==Problem Statement== Use a real life application of the Fourier Series or the Laplace Transform ==Solution== For my project I will plot a car's wheel running over the rumble st…'
 
No edit summary
 
Line 3: Line 3:


==Solution==
==Solution==
[[File:Squarewave.png|600px|thumb|right]]
For my project I will plot a car's wheel running over the rumble strips on the side of the road using the Fourier Transform. This will take the form of a simple square wave as the tire will fall into the dip and bump back up to the surface of the the rumble strip. The Fourier Transform may be used because the action of the car running over the rumble strips can be considered periodic.
For my project I will plot a car's wheel running over the rumble strips on the side of the road using the Fourier Transform. This will take the form of a simple square wave as the tire will fall into the dip and bump back up to the surface of the the rumble strip. The Fourier Transform may be used because the action of the car running over the rumble strips can be considered periodic.


Line 14: Line 15:




by inspection, <math>a_0 = 1</math>
By inspection, <math>a_0 = 1</math> and <math>T=4</math>
 
 
Also by inspection since the "wave" has even symmetry, we know that <math>b_n = 0</math>
 
 
Now we must find <math>a_n</math>. Recall that
 
 
 
:::<math>a_n = \frac{2}{T} \int_0^T f(t) cos(\frac{2 \pi n t}{T}) dt</math>
 
::::<math>= \frac{2}{T} \int_\frac{-T}{2}^\frac{T}{2} 2 cos(\frac{2 \pi n t}{T}) dt</math>
 
::::<math>=\frac{2}{\pi n} [sin(\frac{\pi n}{2}) - sin(\frac{- \pi n}{2})]</math>
 
::::<math>=\frac{4 sin(\frac{\pi n}{2})}{n \pi}</math>
 
So our Fourier Series is:
 
<math>1 + \Sigma_{n=1}^\infin [\frac{4 sin(\frac{\pi n}{2})}{n \pi} cos(\frac{2 \pi n t}{T})]</math>
 
 
==Octave Script==
clear all;
close all;
t=-10:.01:10;
T=4
M=100 
sum1=1;
for n=1:2:M,
sum1 = sum1+(4*sin(pi*n/2)/(n*pi))*cos(2*pi*n*t/T);
end
plot(t,sum1,'b-')
title('Fourier Series Representation of a car moving over rumble strips on a road.')
xlabel('time (seconds)')
ylabel('position of tire')
grid on;
print("squarewave.png","-dpng")

Latest revision as of 22:41, 1 November 2010

Problem Statement

Use a real life application of the Fourier Series or the Laplace Transform

Solution

Error creating thumbnail: File missing

For my project I will plot a car's wheel running over the rumble strips on the side of the road using the Fourier Transform. This will take the form of a simple square wave as the tire will fall into the dip and bump back up to the surface of the the rumble strip. The Fourier Transform may be used because the action of the car running over the rumble strips can be considered periodic.


For simplicity. I will assume the height and width of each rumble strip is 2 inches with 2 inches gap between. I will also place the origin of the coordinate system in the middle of a rumble strip to give it even symmetry.


Recall that the equation for a Fourier Series is:

a0+Σn=1[ancos(2πntT)+bnsin(2πntT)]


By inspection, a0=1 and T=4


Also by inspection since the "wave" has even symmetry, we know that bn=0


Now we must find an. Recall that


an=2T0Tf(t)cos(2πntT)dt
=2TT2T22cos(2πntT)dt
=2πn[sin(πn2)sin(πn2)]
=4sin(πn2)nπ

So our Fourier Series is:

1+Σn=1[4sin(πn2)nπcos(2πntT)]


Octave Script

clear all;
close all;
t=-10:.01:10;
T=4
M=100   
sum1=1;
for n=1:2:M,
sum1 = sum1+(4*sin(pi*n/2)/(n*pi))*cos(2*pi*n*t/T);
end
plot(t,sum1,'b-')
title('Fourier Series Representation of a car moving over rumble strips on a road.')
xlabel('time (seconds)')
ylabel('position of tire')
grid on;
print("squarewave.png","-dpng")