Fourier transform with mixed derivatives/ 2nd order ODE

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
3 replies · 5K views
jianxu
Messages
91
Reaction score
0

Homework Statement


Hi, So I'm suppose to solve the following problem:

[tex]\left.\frac{d^{2}u}{dt^{2}}-4\frac{d^{3}u}{dt dx^{2}}+3\frac{d^{4}u}{dx^{4}}=0[/tex]

[tex]\left.u(x,0) = f(x)[/tex]
[tex]\left.\frac{du}{dt}(x,0) = g(x)[/tex]

Homework Equations



The Attempt at a Solution


First I use Fourier transform on the given expression so that I get the following:

Fourier transform of [tex]\left.\frac{d^{2}u}{dt^{2}}(x,t) = \frac{d^{2}\widehat{u}}{dt^{2}}(\omega ,t)[/tex]

Fourier transform of [tex]\left.\frac{du}{dt}(x,t) = \frac{d\widehat{u}}{dt}(\omega ,t)[/tex]

Fourier transform of [tex]\left.\frac{d^{2}u}{dx^{2}}(x,t) = \left(i\omega\right)^{2}\widehat{u}(\omega ,t) = -\left(\omega\right)^{2}\widehat{u}(\omega ,t)[/tex]

Fourier transform of [tex]\left.\frac{d^{4}u}{dx^{2}}(x,t) = \left(i\omega\right)^{4}\widehat{u}(\omega ,t) = \left(\omega\right)^{4}\widehat{u}(\omega ,t)[/tex]

Which means me overall expression after transform is:
[tex]\left.\frac{d^{2}\widehat{u}}{dt^{2}}(\omega ,t)+4\left(\omega\right)^{2}\frac{d\widehat{u}}{dt}(\omega ,t)+3\left(\omega\right)^{4}\widehat{u}(\omega ,t)=0[/tex]

Now assuming I did that correctly, the next step I think I should proceed with is to solve for [tex]\left.\widehat{u}(\omega ,t)[/tex]. I don't remember how to solve this type of ODE, I was reading a couple of sites and it says I should use a characteristic equation which would I assume then be, [tex]\left.\lambda^{2}+4\omega^{2} \lambda +3\omega[/tex] where [tex]\lambda[/tex] is just an arbitrary symbol to denote a quadratic equation. I looked for the roots and used it along with the general expression of the 2nd order ODE to get
[tex]\left.\widehat{u}(\omega ,t)=c_{1}+c_{2}e^{-4\omega^{2}t}[/tex]
But it seems to be incorrect since I took the derivative and plugged it back into my Fourier transform expression and did not get a 0 for my answer so...Any guidance would be much appreciated! Thanks!
 
Physics news on Phys.org
HI jianxu! :smile:

(have an omega: ω :wink:)
jianxu said:
Which means me overall expression after transform is:
[tex]\left.\frac{d^{2}\widehat{u}}{dt^{2}}(\omega ,t)+4\left(\omega\right)^{2}\frac{d\widehat{u}}{dt}(\omega ,t)+3\left(\omega\right)^{4}\widehat{u}(\omega ,t)=0[/tex]

… I should use a characteristic equation which would I assume then be, [tex]\left.\lambda^{2}+4\omega^{2} \lambda +3\omega^4[/tex] where [tex]\lambda[/tex] is just an arbitrary symbol to denote a quadratic equation. I looked for the roots and used it along with the general expression of the 2nd order ODE to get
[tex]\left.\widehat{u}(\omega ,t)=c_{1}+c_{2}e^{-4\omega^{2}t}[/tex]

No, the roots are ω2 = -1 and -3, so the general solution is c1e2t + c2e-3ω2t :smile:
 
Hi TinyTim, thanks for the reply! I just realized where I made my mistake! Thanks very much for the help!
 
Hello!

I've been working on this problem and was wondering if someone could check if I've done the rest of this problem correctly!

So after finding the roots, I apply the initial conditions where:
[tex]\left.\widehat{u}\left(\omega,0\right) = \widehat{f}\left(\omega\right)[/tex]

since t = 0, I have:
[tex]\left.\widehat{u}\left(\omega,0\right) = \widehat{f}\left(\omega\right) = C_{1}+C_{2}[/tex]

For [tex]\left.\frac{d\widehat{u}}{dt}[/tex]:
[tex]\left.\frac{d\widehat{u}}{dt}= -\omega C_{1}e^{-\omega^{2}t}-3\omega C_{2}e^{-3\omega^{2}t}[/tex]
applying initial conditions:
[tex]\left.\frac{d\widehat{u}}{dt}\left(\omega,0\right)= -\omega C_{1}-3\omega C_{2} = \widehat{g}\left(\omega\right)[/tex]

Now I solved for the constants using elimination and got:
[tex]\left.C_{1}= \frac{3\widehat{f}(\omega)}{2} + \frac{\widehat{g}(\omega)}{2\omega^{2}}[/tex]
and:
[tex]\left.C_{2}= -\frac{\widehat{f}(\omega)}{2} - \frac{\widehat{g}(\omega)}{2\omega^{2}}[/tex]

Therefore our [tex]\left.\widehat{u}(\omega ,t) = \left(\frac{3\widehat{f}(\omega)}{2} + \frac{\widehat{g}(\omega)}{2\omega^{2}}\right)e^{-\omega^{2}t}+ \left(-\frac{\widehat{f}(\omega)}{2} - \frac{\widehat{g}(\omega)}{2\omega^{2}}\right)e^{-3\omega^{2}t}[/tex]

Now I apply inverse Fourier transformation so that:
[tex]\left.u(x,t) = \widehat{f^{-1}}(\omega)[/tex]

Which means:
[tex]\left.u(x,t) = \frac{1}{\sqrt{2\pi}}\int^{\infty}_{-\infty}(\widehat{u}(\omega ,t))e^{i\omega x}d\omega[/tex]

that would be the solution seeing that we don't know what [tex]\left.\widehat{f}(\omega )[/tex] or [tex]\left.\widehat{g}(\omega )[/tex] are.

Thanks!