Fourier transform to solve diff equation

dankaroll
Messages
13
Reaction score
0

Homework Statement



Use Fourier transform to find the solution of the following differential equation:

\frac{\mathrm{d^3}y }{\mathrm{d} x^3}+ \lambda \frac{\mathrm{dy} }{\mathrm{d} x} - xy = 0, \lim_{x \to \infty } y(x)=0

Find the asymptotic of the solution for lambda>> 1. Normalize the solution so y(0) =1.

Homework Equations



Using differentiation properties of Fourier transform,

\frac{\mathrm{d^n}y }{\mathrm{d} x^n} = {(ik)^n }Y[f]

The Attempt at a Solution



using the property

(ik)^3Y[f]+\lambda (ik)Y[f]-xY[f] =0

Y[f]((ik)^3+\lambda ik-x) = 0

so I'm slightly stuck here, since its zero on the right hand side...
 
Physics news on Phys.org
dankaroll said:
Using differentiation properties of Fourier transform,

\frac{\mathrm{d^n}y }{\mathrm{d} x^n} = {(ik)^n }Y[f]
This equation is wrong, look it up.
 
Note that the Fourier transform of xy isn't xY(k). You can't treat x like a constant.

Also, how is the Fourier transform defined in your class? There are several common conventions in use. It would help to know which one you're using to avoid confusion.
 
Klockan3 said:
This equation is wrong, look it up.

<br /> F\left[ \frac{\partial f}{\partial x} \right] = ikF[f] <br />

http://www.thefouriertransform.com/transform/properties.php#derivative

where 2*f*pi = k

vela said:
Note that the Fourier transform of xy isn't xY(k). You can't treat x like a constant.

Looking over some properties, I can't seem to find what this would transform to. I think we can agree that x needs to go on the RHS and y needs to go away somehow...
 
How are you defining the Fourier transform?
 
Let me try this problem.
Let's say I use the definition of the Fourier transform f(x)=\frac{1}{\sqrt {2\pi }}\int _{-\infty}^{\infty } g(y)e^{ixy} dy, g(y)=\frac{1}{\sqrt{2\pi}}\int _{-\infty}^{\infty } f(x)e^{-ixy}dx.
I want to take the Fourier transform of the DE with respect to x.
So I get that it's worth -ik^3 \mathbb{F}(y)+\lambda \mathbb{F}(y)-\mathbb{F}(xy)=0.
\mathbb{F}(xy)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}xye^{-ixy}dx. I'm stuck here. Not sure how to solve this integral. Any idea?
 
Hint: What's F'(y) equal to?
 
vela said:
Hint: What's F'(y) equal to?

Hmm I see I made a mistake, there shouldn't be any k's.
I'm a bit confused because y depends on x.
\mathbb{F}(y&#039;)=\frac{1}{\sqrt{2\pi }} \int _{\infty}^{\infty} y&#039;e^{-ixy}dx.
I was thinking about integration by parts but the fact that y depends on x troubles me. I know I'm almost sure I should express this in terms of \mathbb{F}(y)=\frac{1}{\sqrt{2\pi }} \int _{\infty}^{\infty} ye^{-ixy}dx.
 
I didn't even notice the k's. Since the original problem uses y(x), let's use k to be the variable conjugate to x:
\begin{align*}
f(x) &= \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} F(k)e^{ikx}\, dk \\
F(k) &= \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} f(x)e^{-ikx}\,dx
\end{align*}
So consider F'(k).
 
  • #10
vela said:
I didn't even notice the k's. Since the original problem uses y(x), let's use k to be the variable conjugate to x:
\begin{align*}
f(x) &= \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} F(k)e^{ikx}\, dk \\
F(k) &= \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} f(x)e^{-ikx}\,dx
\end{align*}
So consider F'(k).

Ok!
F&#039;(k)=-\frac{i}{\sqrt {2\pi}} \int _{-\infty}^{\infty}xf(x)e^{-ikx}dx=-i \mathbb{F}(xf(x)). Does this mean that if f(x)=y, I get that F&#039;(k)=-i \mathbb{F}(xy)?
I must say the notation of Mathews and Walker's book somehow confuses me, because the variable conjugate of x is y there.
 
  • #11
Yes, that's right.
 
  • #12
Ok thank you.
I get: \mathbb{F}[xf(x)]=-\frac{d}{idk}\mathbb{F}[f(x)]. I take f(x)=y(x)=y (for notation).
Thus \mathbb{F}(xy)=-\frac{d}{idk}\mathbb{F}(y). Using this, I tansform the original DE into -ik^3\mathbb{F}(y)+\lambda ik \mathbb{F}(y)+\frac{1}{i}\frac{d\mathbb{F}(y)}{dk}=0.
I multiply by i and factorize some terms to get F&#039;(k)+(k^3-\lambda k)F(k)=0.
I used separation of variables to get F(k)=Ae^{\frac{k^2}{2}\left ( \lambda - \frac{k^2}{2} \right ) }.
I doubt this is right because the question says to get the solution for large lambda.
If that is right, now the next step is to invert that expression to get y(x), right?
 
  • #13
I got the opposite sign in the exponent, but there's still the problem when k goes to infinity.
 
  • #14
vela said:
I got the opposite sign in the exponent, but there's still the problem when k goes to infinity.
You mean you have F(k)=Ae^{\frac{k^2}{2}\left ( \frac{k^2}{2}-\lambda \right ) }?
Hmm and what does this problem mean? That the Fourier transform isn't appropriate or...
 
  • #15
Found a sign error. Now my answer matches yours.
 
Back
Top