Evaluation of oscillatory integrals

  1. HOw can you evaluate (for big 'u') the exponential integral

    [tex] \int_{0}^{1}dx f(x)e^{iux} [/tex] u-->oo

    if you use Numerical methods, since exp(iux) oscillates rapidly you cannot apply any of them, also since we have inside exponential (iux) you can not apply 'Stationary phase method' unless you make the change.

    [tex] \int f(g(t)exp(iug(t)) [/tex]

    and then apply stationary phase finding [tex] g'(t)=0 [/tex] to evaluate the integral but i am not pretty sure.
  2. jcsd
  3. You should check out Olver's book Asymptotics and Special Functions. I'm pretty sure you'll find something there.
  4. If f is continuous, then taking f out of the integral with the intermediate value theorem, and calculating the integral of the exponential expression, should enable you to prove that this approaches zero when u approches infinity.
  5. And remember that in

    \int\limits_0^1 f(x) e^{iux} dx = f(\xi_u) \int\limits_0^1 e^{iux} dx

    the [itex]\xi_u[/itex] depends on u. I think suprisingly many of my mistakes have been related to first using some intermediate value theorem, and then ignoring on what the xi depends.
  6. Haelfix

    Haelfix 1,723
    Science Advisor

    You can try a physics trick called "wick rotation" by making a change of variables to the real plane (to get rid of the annoying 'i' in the top integral), you just have to check to see if the analytic continuation makes sense there.

    But i'd first go with the method of steepest descent, it seems tailor made for this type of integral after a suitable variable change f(x) --> h(Z) and ix -- > g(Z). I don't see why its not defined in this case. So long as h(Z) and g(Z) are analytic functions in some contour of the complex plane (check this) and u --> infinity, you're set.

    Of course if f(x) is some pathological beast, you're probably out of luck.
    Last edited: Aug 28, 2007
  7. My two previous posts were example of how first making a mistake, and then fixing it, can lead into unnecessarily complicated way. This is what I should have said

    \big| \int\limits_0^1 f(x) e^{iux} dx\big| \leq \big(\underset{x\in [0,1]}{\textrm{sup}} |f(x)|\big) \big| \int\limits_0^1 e^{iux} dx\big| \underset{u\to\infty}{\to} 0

    If the supremum doesn't exist, then you'll have to know something else, and use it somehow.
  8. Please this problem is similar to your problem. But i could'nt solve please help me

    Attached Files:

  9. Holy crap! I've been writing nonsense here in 2007. :surprised

    If you want to use the intermediate value theorem

    \int\limits_0^1 f(x)\phi(x)dx = f(\xi) \int\limits_0^1 \phi(x)dx

    then [itex]\phi[/itex] should be real and non-negative.

    Also, if you set [itex]f(x)=e^{-i\pi x}[/itex] and [itex]u=2\pi[/itex], then

    \int\limits_0^1 f(x) e^{iux}dx = \frac{2i}{\pi}


    \big(\underset{x\in [0,1]}{\textrm{sup}}\; |f(x)|\big) \int\limits_0^1 e^{iux}dx = 0

    So my inequality wasn't fine.
Know someone interested in this topic? Share a link to this question via email, Google+, Twitter, or Facebook