This week we began our study of Fourier Analysis and Fourier Series. Our system of study for this discussion will be a string on which the wave velocity is , and which has fixed endpoints at x=0 and x=L where x is a position along the length of the string.
Recall that the wave equation is
and that we can write any solution to this wave equation as a superposition of normal modes:
where .
We would like to go from a set of initial conditions (position and velocity of every segment of the string) to the values of each and . We were able to do this for systems with a finite number of degrees of freedom, so we should, in principle, be able to do it in the case of systems with a very large number of degrees of freedom.
Let us first consider just a set of positions at time $t =0$. Let’s say we are given these, and they take the form , where is any function with . These are the boundary conditions obeyed by our segment of string.
If we now apply this to the expression for y in terms of the normal modes, we find
Now for the moment, let us combine the coefficients of each into a set of coefficients , in which case we then have
Let us now try to get the coefficients from the function .
To do this, we use a trick involving integrals of sin functions:
multiplying both sides of the equation for f(x) by and integrating over $x$ from $0$ to $L$, we find
This means that if we know the function f(x), and perform some integrals, then we can get the coefficients!
See this weeks Mathematica code to begin to get an idea of how this works:
https://jhubisz.expressions.syr.edu/phy360/wp-content/uploads/sites/5/2016/11/Week-10.nb
We can then go a bit further. Again, the goal is to get both the and . In order to do this, we need not just initial position information, but also initial velocity information. Right now, we have , so we only have fixed a combination of and . Let us say we have that velocity information as a given:
.
In this case, we can write , and thus we can also express the velocity of each point of the string in terms of a sum over sin functions:
We can then use the same trick to get:
So now you have all of the B’s, and all of the C’s (since we presumed we know both initial position, f(x), and initial velocity, g(x)). We then have
From this, we can finally get both and :
.
I will soon be uploading a mathematica notebook that goes through an example of a particular set of initial positions and velocities, and go through how to get the A’s and ‘s for a specific example.