Fourier Transform of f(x) = N*exp(-ax^2/2) and Simplifying Dirac Delta Functions

In summary, the conversation discusses finding the Fourier transform of a function and simplifying expressions involving the Dirac delta function. The solution for the Fourier transform is given, but the method used is questioned. The result for the first Dirac question is clarified, and a sketch is provided for the second and third questions. The conversation ends with a request for help with computing two integrals involving the function and its transform.
  • #1
MathematicalPhysicist
Gold Member
4,699
371
I need to find the fourie transform of f(x)=N*exp(-ax^2/2).
(N and a are constants).
well ofcourse iv'e put into the next integral:
[tex]\int_{-\infty}^{\infty}f(x)exp(-ikx)dx[/tex]
Iv'e changed variables, that i will get instead of exp(-ax^2/2)exp(-ikx),
exp(-z^2)exp(-ikz*sqrt(2/a)), but that didn't helped very much, perhaps I need to use here the definition of dirac delta function to simplify this.

any hints?
there's also a question to simplify the next terms:
1.(x^2+3)d(x+5) where d is dirac's delta function.
2.d(2x-8)
3.d(x^2+x-2)

as for 2, what i did is: d(2x-8)=d(2(x-4))=d(x-4)/2
for 3 i think it equals:
d(x^2+x-2)=d((x-1)(x+2))=d(x-1)+d(x+2)
but i didnt use the integral definition, iv'e just used the fact that at x=1 or x=-2 the argument of the functional is zero and thus by definition it diverges there otherwise it equals zero, is this a sound reasoning for this question?
for 1, not sure, but if we use the integral definition then [tex]\int_{-\infty}^{\infty}f(x)(x^2+3)d(x+5)dx=f(-5)*28=28*\int_{-\infty}^{\infty}f(x)d(x+5)dx[/tex] then (x^2+3)d(x+5)=28*d(x+5), or something like this?

thanks in advance.
 
Physics news on Phys.org
  • #2
For the Fourier transform I got a solution, however I am not certain whether it is the correct one. Please, use this with care. I rewrote the argument of the exponential function as:

[tex]\frac{a}{2}x^2+ikx=\frac{a}{2}\left[x^2+\frac{2}{a}ikx\right]= \frac{a}{2}\left[x^2+2\frac{ik}{a}x+\left(\frac{ik}{a}\right)^2\right]+ \frac{k^2}{2a}= \frac{a}{2}\left(x+\frac{ik}{a}\right)^2+\frac{k^2}{2a}[/tex]

Using then:

[tex]\sqrt{\frac{a}{2}}\left(x+\frac{ik}{a}\right)=t[/tex]

You arrive at:

[tex]F(k)=N\sqrt{\frac{2\pi}{a}}e^{-\frac{k^2}{2a}}[/tex]

Please check this result because the last substitution is a bit awkward and can be regarded as "playing" with real and imaginary numbers. Therefore it might very well be that this is complete and utterly wrong.
 
  • #3
not sure myself if your answer is the correct one, i.e obviously when you turned to complex dumby variable you need to use the residue theorem, but my problem is that i don't think that there's a pole there, so the sum of the residues will be zero, i.e you integrate over a straight line of y=k/a in the complex plane, but there are no poles.
 
  • #4
loop quantum gravity said:
not sure myself if your answer is the correct one, i.e obviously when you turned to complex dumby variable you need to use the residue theorem, but my problem is that i don't think that there's a pole there, so the sum of the residues will be zero, i.e you integrate over a straight line of y=k/a in the complex plane, but there are no poles.

Why the residue theorem? This is for evaluating the value of a contour integral in the complex plane circling around a few poles (if any). The theorem gives a practical way of determining the value of this integral. I don't think that this has something to do with the integral you asked for. I could be mistaken, then I'm missing something.

However I had a second look at the result I presented and it is indeed the correct one. I found this in a table in the book "Mathematical Handbook for Scientists and engineers", by G.A. Korn and T.M. Korn. So, now we can safely go to sleep :smile:

The first one of the dirac questions is (if I remember correctly, it has been a while since I used this) the value of the function that is multiplied with the dirac at the point where the argument of the dirac is zero. This means that it is nothing more than the value of the function in the point x=-5, giving 28. Again, check this please :wink:
 
  • #5
so it should be:
(x^2+3)d(x+5)=28, are you sure?
then why is my way wrong there, i got: (x^2+3)d(x+5)=28d(x+5)
ok i can see what you mean if i integrate this function (without entering f(x), i will get 28), but still i don't see how this means what your'e asserting, i mean if (x^2+3)d(x+5)=28, then if i enter 28 in the the integral, i will get an integral which doesn't converge.
 
  • #6
loop quantum gravity said:
so it should be:
(x^2+3)d(x+5)=28, are you sure?
then why is my way wrong there, i got: (x^2+3)d(x+5)=28d(x+5)
ok i can see what you mean if i integrate this function (without entering f(x), i will get 28), but still i don't see how this means what your'e asserting, i mean if (x^2+3)d(x+5)=28, then if i enter 28 in the the integral, i will get an integral which doesn't converge.

You are correct. The result is indeed zero everywhere except at the point x=-5 where it is 28. That is mathematically 28d(x+5). That is what I meant, but confusion can arise if one considers the function 28 to be the constant function for all x. Shame on me :redface:
 
  • #7
For the third dirac question, I made a sketch and I think it is correct as well.
[tex]\delta(x^2+x-2)=\delta((x-1)(x+2))=\delta(x-1)+\delta(x+2)[/tex]
I have to admit that I can't remember ever to have used this kind of result. The strange thing is that the argument of the function consists of two points where the argument is zero, so it is a kind of "double dirac function". Very nice.

For the second one I assume by considering the definition of the dirac and applying a transformation gets you the result.
 
  • #8
another computation of integrals, after we found g(k) above with the transform, I need to compute the next integrals:
[tex]s=\int_{-\infty}^{\infty}f(x)x^2dx[/tex]
and [tex]t=\int_{-\infty}^{\infty}g(k)k^2dk[/tex]
don't see how to do this really, by parts won't work here.
 
  • #9
Assuming you mean:

[tex]t=\int_{-\infty}^{\infty}g(k)k^2dk= \int_{-\infty}^{\infty} N\sqrt{\frac{2\pi}{a}}e^{-\frac{k^2}{2a}} k^2dk[/tex]

It is solvable by partial integration. Consider the following:

[tex]t=N\sqrt{\frac{2\pi}{a}} \int_{-\infty}^{\infty} e^{-\frac{k^2}{2a}} k^2dk= N\sqrt{\frac{2\pi}{a}} \cdot I[/tex]

With:

[tex]I=\int_{-\infty}^{\infty} e^{-\frac{k^2}{2a}} k^2dk= (-a)\int_{-\infty}^{\infty} k\cdot e^{-\frac{k^2}{2a}}d\left(-\frac{k^2}{2a}\right)[/tex]

Which gives after partial integration:

[tex]I=(-a)\left[k \cdot e^{-\frac{k^2}{2a}}\right]_{-\infty}^{\infty} +a\cdot \int_{-\infty}^{\infty} e^{-\frac{k^2}{2a}}dk= (-a)\cdot (0-0)+a\cdot \int_{-\infty}^{\infty} e^{-\frac{k^2}{2a}}dk[/tex]

Setting now:

[tex]\frac{k}{\sqrt{2a}}=z[/tex]

gives:

[tex]I=a\cdot \int_{-\infty}^{\infty} \sqrt{2a} \cdot e^{-z^2}dz[/tex]

Which is:

[tex]I=\sqrt{2\pi}a^{3/2}[/tex]

The result we are looking for is now:

[tex]t=N2\pi a[/tex]

I do not understand the first one. Can you explain this one a little bit more? Can you write the complete integral?
 
  • #10
Sorry for the flood of questions.
I need to calculate the integral:
[tex]\int_{-\infty}^{\infty}(sin(x)/x)^4dx[/tex]
by using parseval's equation.
now parseval's equation is:
[tex]\int_{-\infty}^{\infty}|f(x)|^2dx=\frac{1}{2\pi}\int_{-\infty}^{\infty}|g(y)|^2dy[/tex] where g(y) is F.T of f(x), obviously here f(x)=(sin(x)/x)^2, so first i need to integrate:
[tex](sin(x)/x)^2e^{-iyx}[/tex]
i tried using residues, but the residue equals zero, and the contour on the big circle goes to zero and so is the small circle around zero (because it equals -pi*Res(z=0), so I don't know what to do here.
mathematica has a nice term for the integral i need to calculate.
 
  • #11
coomast said:
I do not understand the first one. Can you explain this one a little bit more? Can you write the complete integral?
f(x) is the first function i was given in my first post, i.e f(x)=Nexp(-ax^2/2).
anyhow thanks i understand this part.
 
  • #12
OK, so the first integral is thus:

[tex]s=\frac{N\sqrt{2\pi}}{a^{3/2}}[/tex]

Using the same method as the other one. I will have to look a few things up for the the other questions. "I'll be back" :smile:
 
  • #13
I think it's best I will borrow this book by the korns', it could be handy in the open material exam.
 
  • #14
The FT of (sin(x)/x)^2 seems not so easy. I will do my best, however I don't use contour integration that often, so I will need some time. It seems that it is the appropriate way in this case.

The book by Korn and Korn is a thick one giving lot's of mathematical material and formulas, but I would use something else for an open exam. For me the one which gave me the best preparation was the one of the schaum series: "Laplace Transforms" by M. Spiegel. Some people find these books rubbish, I find them an excellent source for problems. Not all of them are good, but the ones written by the late prof. Spiegel are very good. It just takes a bit of time to solve the problems.

[begin edit]
Another book which can be useful is the one "Complex Variables" also from the schaum series and written by Spiegel.
[end edit]

I will further try to find the FT of the function. Just for information, do you have to use Parseval's identity? If I can come up with a different method is that allowed?
 
Last edited:
  • #15
well yes, i need to use parseval's equation.

thanks for the recommendations, I also think that for this theoretical methods for physicists course schaum's will do, I also thought of browsing mary boas' book.
 
  • #16
The integral:

[tex]\int_{-\infty}^{\infty}\left(\frac{sin(x)}{x}\right)^4dx[/tex]

has the solution [tex]\frac{2\pi}{3}[/tex] as you might have obtained from mathematica. I got it from maxima. Anyway it looks as if this is solvable by contour integration, but if you have to use parseval it makes things more complicated. I keep searching for an answer.
 
  • #17
loop quantum gravity, It seems very difficult or I can't see the solution in a practical way. The problem is that I can't find the FT of the sin^2(x)/x^2 function. I tried several different contours but I end up nowhere. Then I tried to look at it as:

[tex]\int_{-\infty}^{\infty}\frac{sin^2(x)}{x^2}e^{-ixy}dx= 2 \cdot \int_{0}^{\infty}\frac{sin^2(x)}{x^2}e^{-sx}dx[/tex]

with the iy replaced by s (another complex variable) and twice the integral of interest because of the even character of it. It looks now as a Laplace transform of the function. So I found this Laplace tranform in a tables book ("table of laplace transforms, by Roberts and Kaufman) to be the following:

[tex]L\left\{\frac{sin^2(at)}{t^2}\right\}=a\cdot atan\left(\frac{2a}{s}\right) -\frac{s}{4}log\left(1+\frac{4a^2}{s^2}\right)[/tex]

Which can be found by considering the property of the Laplace transform that says that if the function is divided by t, you have to integrate the laplace transform. In this case twice. However, if I write this out, replace s by iy, replace the atan of the complex argument, take the absolute value and then the square, I end up with a difficult integral. Maxima gives me a result in the form of a limit.

So, I don't know what else I could use. Perhaps someone else here can help you. I assume that it must be possible using contour integration, but I can't find the contour.

My apologies that I can't be of any more assistance.

[edit]
If you are looking for a good book with formulas, I can recommend "Mathematical Handbook of Formulas and Tables", written by Spiegel, also one out of the schaum series. It does not only cover basic stuff, but goes into advanced functions as well. I use this one very often.
 
Last edited:
  • #18
thanks coomsat, I appreciate your help nontheless.
 

1. What is the Fourier Transform of f(x) = N*exp(-ax^2/2)?

The Fourier Transform of f(x) = N*exp(-ax^2/2) is given by F(k) = N*sqrt(2*pi/a)*exp(-k^2/(4a)). This is a Gaussian function in the frequency domain with a peak at k=0 and a width proportional to the inverse of the square root of the parameter a.

2. How is the Fourier Transform of a Gaussian function derived?

The Fourier Transform of a Gaussian function can be derived using the definition of the Fourier Transform, integration by substitution, and completing the square in the exponent. The result is a Gaussian function in the frequency domain with a peak at k=0 and a width proportional to the inverse of the square root of the parameter a.

3. Can the Fourier Transform of f(x) = N*exp(-ax^2/2) be simplified further?

No, the Fourier Transform of f(x) = N*exp(-ax^2/2) cannot be simplified further. It is already in its simplest form and cannot be expressed in terms of elementary functions.

4. What is the significance of the Dirac Delta function in the Fourier Transform of a Gaussian?

The Fourier Transform of a Gaussian function has a Dirac Delta function at k=0, indicating that the original function f(x) has a constant value N at all points in space. This represents the fact that the Gaussian function is spread out over all frequencies, with no specific frequency dominating.

5. How can the Dirac Delta function be simplified in the context of Fourier Transforms?

The Dirac Delta function can be simplified in the context of Fourier Transforms by using its defining properties, such as shifting, scaling, and sifting. These properties allow us to manipulate the Dirac Delta function and simplify it in certain cases, such as when it is multiplied by a constant or integrated with other functions.

Similar threads

  • Calculus and Beyond Homework Help
Replies
5
Views
355
  • Calculus and Beyond Homework Help
Replies
3
Views
762
  • Calculus and Beyond Homework Help
Replies
31
Views
2K
  • Calculus and Beyond Homework Help
Replies
8
Views
998
  • Calculus and Beyond Homework Help
Replies
6
Views
2K
  • Calculus and Beyond Homework Help
Replies
3
Views
2K
  • Calculus and Beyond Homework Help
Replies
7
Views
1K
  • Calculus and Beyond Homework Help
Replies
5
Views
1K
Replies
5
Views
966
  • Calculus and Beyond Homework Help
Replies
1
Views
2K
Back
Top