Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Fourier Transform: Why is this so?

  1. Aug 18, 2008 #1
    Why is this so?

    [tex]\displaystyle F\left[ \frac{1}{1-e^{-\pi x}} \right] = i \frac{1+e^{-2k}}{1-e^{-2k}}[/tex]

    Here, [tex]-\infty < x < \infty[/tex]. It has to be done by contour integration, by the way. Unfortunately, I'm having difficulty with the whole thing.
  2. jcsd
  3. Aug 18, 2008 #2
    Are you sure that is so? To me it seems that the Fourier transform isn't defined for that function. The integral diverges at x=0.
  4. Aug 18, 2008 #3
    Yes, I'm (99%) sure.

    The function diverges (not only at x = 0, but also at x = infinity), but not faster than [tex]e^{\alpha x}[/tex]. That should be enough to assure existence of a Fourier transform.
  5. Aug 18, 2008 #4
    If I try to keep things simple, then the situation is that the integral

    \int\limits_{-\infty}^{\infty} \frac{e^{ikx}}{1-e^{-\pi x}} dx

    does not converge, and the Fourier transform doesn't get defined. I don't know what rate of divergence could mean in this context.

    Actually there is two divergence problems. One is at [tex]x=0[/tex], and other one at [tex]x=\infty[/tex].

    It could be there is something going on that I'm not aware of. I'm probably not wrong to guess that this has something to do with some regularizations?
  6. Aug 18, 2008 #5
    To be perfectly honest, I'm not THAT sure of anything.

    The standard idea I've known is to instead examine the fourier transform of [tex]e^{\alpha x}f(x)[/tex]. When one does the fourier inversion, it amounts to integrating not on the real k axis, but along the line [tex][-\infty + i\alpha,-\infty + i\alpha][/tex] (there may be a sign off there). This is why we don't need functions that have convergent integrals -- so long as they don't shoot off faster than an exponential.

    The problem is that this gives you the fourier transform [tex]F(k+\alpha)[/tex], not F(k). It of course isn't a problem once you invert (you get f(x) in terms of an integral -- just not the one on the real k-axis). But the actual F(k) expression...?

    Indeed, when I went down this line of thought, my impression was F(k) doesn't exist. But this is a published result. Moreover, the rest of the results makes sense, which is why I don't doubt this claim. The author writes that it follows from contour integration. No other hint is given.
  7. Aug 18, 2008 #6


    User Avatar
    Science Advisor
    Homework Helper

    What is your definition of "Fourier transform"?
  8. Aug 18, 2008 #7
    Standard definition. Possibly extended to the complex plane.
  9. Aug 18, 2008 #8
    Nevermind. It's been solved at SOS.

    Idiot mistake I made not remembering [tex]1-e^{\pi z}[/tex] has more than the one root at z = 0.
  10. Aug 18, 2008 #9
    What does the notation P.V. mean?
  11. Aug 18, 2008 #10

    Basically, the integral without the singularity. The poster misused some of the notation. The fourier integral is not defined into a principal value, but as "the integral over the washer". It's by virtue of the fact the outer radial integral tends to zero that the principal value is equal to the residues and thus equal to the fourier integral (modulo signs and stuff).

    But of course, I'm not complaining. I got my answer.
  12. Aug 18, 2008 #11
    Ok. It should be stated in the original problem, that the divergence at the origo shall be dealt with the Cauchy principle value. I would be curious to know how reliable the inverse Fourier transform is, once the Fourier transform has been calculated with Cauchy principle value choice. You know, you can deal with the divergence in many different ways, and the inverse formula is not going to hold for them all, at least.
  13. Aug 18, 2008 #12
    I'm still skeptical over this entire result. The expression

    \frac{e^{-ikz}}{1 - e^{-\pi z}}

    does not approach zero on the limit [tex]z\to (+\infty, 0)[/tex], and does not oscillate rapidly either to make the integral go zero. So the vanishing of the integral on the outer semi-circle looks suspicious too.
  14. Aug 18, 2008 #13
    That's standard...

    The Fourier transform was never calculated with just the Cauchy principal value (as I said in my post above, correcting the solution).

    It was calculated by applying the definition: an integral from x = -infinity to +infinity. The contour started at -infinity, ran along the x-axis (this is the principal value), plus the indentations above the singularities, and over the outside of the "washer". It just so happens that the function decays in the outside washer.

    There's no ambiguity here. It's the standard way of dealing with fourier transforms in the complex plane.
  15. Aug 18, 2008 #14
    As z tends to infinity, the integrand behaves like exp(-ik + pi)z. Thus for Im(k) < -pi, the integral tends to zero. I think that's enough, but I'm not quite sure about the portion near the real-axis. Gotta go, but I'll think about it.
  16. Aug 18, 2008 #15
    \lim_{z\to (+\infty, 0)}\frac{1}{1-e^{-\pi z}} \to 1

    so the integrand behaves like

    \frac{e^{-ikz}}{1-e^{-\pi z}} \to e^{-ikz}.

    There is no chance this integral is converging, and also no chance that the residue trickery could have been valid.

    You could interpret the integral as a distribution though, like first fixing the integration domain as finite, and then taking the limit when the domain extends to all real axis, and solving how the integral function behaves as a distribution on this limit.
  17. Aug 18, 2008 #16

    You're right. I was too hasty writing what I wrote when I wrote it. Here is a reply, but it may have an error. I'm in a hurry though.


    [tex]f_{+}(x) = \frac{1}{1-e^{-\pi x}}[/tex] for [tex]x > 0[/tex] and,

    [tex]f_+ = 0[/tex] otherwise.

    Let [tex]F_+(x) = e^{\alpha x}f_+(x)[/tex] where [tex]\alpha < 0[/tex] and real.

    Now calculate,

    [tex]\hat{F_+}(k) = \textbf{F}[F_+][/tex]

    Notice now the integrand tends to zero for,

    [tex]\text{Im}(k) = 0[/tex] and [tex]z \to \infty[/tex]

    By the same method as above,

    [tex]\hat{F_+}(k) = i\left(\frac{1+e^{-2(k+\alpha i)}}{1-e^{-2(k+\alpha i)}}\right)[/tex]

    It is consistent (by extension of the Inverse Fourier Transform to the complex domain) to define,

    [tex]\hat{F_+}(k) = \hat{f_+}(k+i\alpha)[/tex]


    [tex]\hat{f_+} =i\left(\frac{1+e^{-2x}}{1-e^{-2x}}\right)[/tex]

    The same follows for [tex]f_-(x)[/tex] defined similarly and the result follows from superposition.
  18. Aug 19, 2008 #17
    Where I come from, what you're talking about is called the Laplace transform. The idea is that you replace the [itex]i\omega[/itex] in the Fourier transform with a complex variable [itex]s=\sigma+i\omega[/itex]. Then, the transform integral is exactly the Fourier transform of [itex]e^{\sigma x}f(x)[/itex]. In this setting, [itex]\sigma[/itex] is also considered a variable, and so you get a function of [itex]s[/itex], along with an associated region of convergence (typically a range of [itex]\sigma[/itex] values that the integral converges for). The [itex]\sigma[/itex] component is sort of "extra" information in the sense that the time-domain function can be recovered by a line integral along [itex]\omega[/itex] for any value of [itex]\sigma[/itex] in the region of convergence, as you are discussing here.

    It bears emphasizing that when the [itex]i\omega[/itex] axis is not included in the region of convergence, the Fourier transform does not, strictly speaking, exist. The Laplace transform will exist, though, provided f(x) can be made convergent through exponential damping, as you say. Note that, for example, an f(x) that blows up for both large positive and large negative x cannot be made to converge, as the end will still blow up with [itex]\sigma[/itex] positive, and vice-versa with [itex]\sigma[/itex] negative. You can't get around this by chopping it into two pieces in the time domain, as the corresponding regions of convergence of the transforms will not overlap, and so there is no way to put the two transforms back together.
  19. Aug 19, 2008 #18
    Well, not exactly the Laplace Transform. The Laplace Transform is the Fourier Transform for the [itex]f_+[/itex] I defined, with [itex]k = ip[/itex] and thus the inversion contour inversion is taken along [itex]\alpha -i\infty[/itex] to [itex]\alpha +i\infty[/itex] -- a vertical line instead of a horizontal one. Thanks for pointing this connection out

    Hmm...Interesting. I never thought about that.

    But what if I did tried this for the fourier transform of [itex]f(x) = e^{|x|}[/itex], clearly divergent at both ends.

    Examine [itex]F(x) = e^{-|2x|}f(x)[/itex] and define as before [itex]F_{\pm}(x)[/itex] as the chopped-in-two-pieces version with the other side set to zero.

    Then for [itex]\hat{F}(k_1 + ik_2)[/tex] is defined (I believe(?)) for [itex]x>0[/itex] when [itex]k_2 > -1[/itex] and for [itex]x<0[/itex] for [tex]k_2 < 3[/tex] -- an overlapping strip.

    Why would this not work?
  20. Aug 19, 2008 #19
    I'm not sure I see the difference between that and a Laplace transform...?

    Well in this case, you wouldn't even need to chop the function apart, right? You'd be taking the Fourier transform of [itex]e^{-|x|}[/itex], which is a convergent function. The trouble with this approach is that the transform of [itex]e^{-|\sigma x|}f(x)[/itex] is not related to the transform of [itex]f(x)[/itex] in a nice way, as is the case for the transform of [itex]e^{-\sigma x}f(x)[/itex]. So it's not clear to me that this approach would get you very far, outside of very specific functions like [itex]f(x)=e^{-|ax|}[/itex]...
  21. Aug 19, 2008 #20
    Yes, but your original statement was that we can't take fourier transforms of functions which have divergence at both tails.

    I was wondering whether the statement was a pragmatic one or a general one.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook