How to 'shift' Fourier series to match the initial condition of this PDE?

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
4 replies · 2K views
Master1022
Messages
590
Reaction score
116
Homework Statement
How can I match the following Fourier series to match the initial condition for a PDE?
Relevant Equations
Fourier Series
Hi,

Question: If we have an initial condition, valid for [itex]-L \leq x \leq L[/itex]:
[tex]f(x) = \frac{40x}{L}[/tex] how can I utilise a know Fourier series to get to the solution without doing the integration (I know the integral isn't tricky, but still this method might help out in other situations)?

We are given the following Fourier series:
[tex]\frac{2A}{\pi} \left( sin(\phi) + \frac{1}{2} sin(2 \phi) + \frac{1}{3} sin(3\phi) + ... \right)[/tex] where [itex]\phi = \frac{2 \pi x}{T}[/itex] and [itex]A[/itex] is the amplitude. This represents a function that looks like this (this would be sharper with more terms, but I just put a few into Desmos):
Screen Shot 2020-08-24 at 9.08.11 PM.png


My attempt:
I think we want to do the following steps:
1. Reflect the series about the x-axis (multiply by -1)
2. Shift/scale the series have a zero at the origin and have a full period range from [itex]-L[/itex] to [itex]L[/itex]

It is the second step that confuses me and I am not completely sure how to do it, but this is how I approached it:
- we want the point where [itex]\pi[/itex] is to shift to the origin (so [itex]\phi + \pi[/itex])
then we get that [itex]sin(\phi + \pi) = -sin(\phi)[/itex].
- then we let T = 2L to get and A = 40 to get:
[tex]\frac{-2\times 40}{\pi} \left(-sin(\frac{\pi x}{L}) + \frac{1}{2} sin(\frac{2 \pi x}{L}) - \frac{1}{3} sin(\frac{3 \pi x}{L}) + ... +\frac{(-1)^n}{n}sin(\frac{n \pi x}{L}) \right)[/tex]

Does this method seem correct?

However, what if there is a shift and scale required? How does that change the order of the transformation? Should we always do the shift first, for out the trig double angle simplification, and then do the scaling in the x-direction?

Any help is greatly appreciated.
 
  • Like
Likes   Reactions: etotheipi
Physics news on Phys.org
If you need an affine transformation, then you should do that in a single step: [tex] [a,b] \to [0, 2\pi]: x \mapsto \frac{2\pi}{b - a}(x - a).[/tex] Thus here [tex]\phi = \frac{\pi}L (x + L), \qquad x = -L + \frac{L}{\pi} \phi[/tex] and [tex] \frac{40x}{L} = -40 + \frac{40}{\pi} \phi.[/tex]
 
  • Informative
Likes   Reactions: etotheipi and Master1022
pasmith said:
If you need an affine transformation, then you should do that in a single step: [tex] [a,b] \to [0, 2\pi]: x \mapsto \frac{2\pi}{b - a}(x - a).[/tex] Thus here [tex]\phi = \frac{\pi}L (x + L), \qquad x = -L + \frac{L}{\pi} \phi[/tex] and [tex] \frac{40x}{L} = -40 + \frac{40}{\pi} \phi.[/tex]

Thank you for your reply! Okay, this method makes sense. However, why did you substitute in the expression for [itex]x[/itex] into the initial condition? I suppose that might be feasible if the [itex]\phi[/itex] was defined within the problem?

I think for this case, I would just use the expression [itex]\phi = \frac{\pi}L (x + L)[/itex] and substitute into the Fourier series to have a match to my initial condition? I think you know that though, I just wanted to confirm. By substituting in your expression for [itex]\phi[/itex] to the Fourier series, it yields the same trig identity so that is very re-assuring.

Thanks
 
pasmith said:
If you need an affine transformation, then you should do that in a single step: [tex] [a,b] \to [0, 2\pi]: x \mapsto \frac{2\pi}{b - a}(x - a).[/tex] Thus here [tex]\phi = \frac{\pi}L (x + L), \qquad x = -L + \frac{L}{\pi} \phi[/tex] and [tex] \frac{40x}{L} = -40 + \frac{40}{\pi} \phi.[/tex]
Sorry to bring this back up, but what if the original wave described by the Fourier series was shifted by ## \pi ## (or any other phase)? Would I then:
1) use a phase shift ## \phi = \phi + \pi ## to shift the Fourier series wave to where I want
2) expand with double angle formula and simplify
3) use the affine transformation to map from ## \phi ## to ## x ##

Would that process be correct?

Thank you
 
Master1022 said:
Thank you for your reply! Okay, this method makes sense. However, why did you substitute in the expression for [itex]x[/itex] into the initial condition? I suppose that might be feasible if the [itex]\phi[/itex] was defined within the problem?

It's the natural thing to do: To evaluate [tex]\int_a^b f(x) \sin\left(\frac{2n\pi(x - a)}{b - a}\right)\,dx = \frac{b - a}{2\pi}\int_0^{2\pi} f\left(a + \frac{(b-a)}{2\pi}\phi\right) \sin(n\phi)\,d\phi[/tex] one has to. Also, in this special case of a linear boundary condition, it is obvious from [tex] \frac{40x}{L} = -40 + \frac{40}{\pi}\phi[/tex] that the initial condition is -40 plus [itex](40/\pi)[/itex] times the unshifted series for [itex]\phi[/itex]. This is easier to calculate than the shifted series.