Pendulum with Dirac Comb excitation l \ddot\theta+b\dot \theta+g\theta=G\,\sum_{-\infty}^\inf

Taliyah Spencer

Taliyah Spencer

Answered question

2022-04-23

Pendulum with Dirac Comb excitation
lθ¨+bθ˙+gθ=G,δ(tnT)
where l, b, g, G are constants and T=2πω.
Show that the resulting motion is given by
θ(t)=GTlω2+2Gcos(ωtπ2)Tbω+ [terms with frequencies 2ω]
and explain why the higher frequency terms are supressed.
My first take was to rearrange to
θ¨+blθ˙+ω2θ=Gl,δ(tnT)
where ω=gl
Then, taking the Laplace transform of both sides I got
Θ(s)=Gl,ω2(b2l)2,ω2(b2l)2(s+b2l)2((b2l)2ω2),11esT
which, as far as I'm concerned transforms to
θ(t)=Gl,ω2(b2l)2n=0,H(tnT),eb2l(tnT)sin(ω2(b2l)2(tnT))
And, assuming that this is a correct form of the solution, I can't see how that is equivalent with the function given in the question. I reckon it has something to do with using Fourier series/transform instead? If so, I'm not sure how to do that. Or, is there a way to convert my solution into the given one?
I've been struggling with this for a good few days now, so any help would be much appreciated.

Answer & Explanation

Ezakhenitne

Ezakhenitne

Beginner2022-04-24Added 10 answers

Step 1
The distinction between apparent solutions via Laplace transform and via Fourier transform depends meaningfully on the function-space that the "target" function is in, and the desired solution is in.
Suppressing some of the constants to reduce clutter, consider the equation θ +aθ+θ=nZδn
Under various external context assumptions, we might suppose that the solution θ is at worst a tempered distribution (the right-hand side is such). This excludes exponentially-growing (either forward or backward) solutions θ. In particular, unless the constants a,b on the left-hand side are such that solutions to the characteristic equation λ2+aλ+b=0, the homogeneous equation has no tempered-distributions solutions (except 0).
Step 2
The Poisson summation formula is the assertion that nδn is its own Fourier transform. Thus, again up to essentially irrelevant constants,
-x2θ^+bixθ^+cθ^ = nδn
which suggest that
θ^ = nδn-x2+bix+c = nδn-n2+bin+c
Transforming back, again suppressing constants, θ(t) = ne2πint-n2+bin+c.

Do you have a similar question?

Recalculate according to your conditions!

Ask your question.
Get an expert answer.

Let our experts help you. Answer in as fast as 15 minutes.

Didn't find what you were looking for?