Asymptotic expansion of an integral

Click For Summary
SUMMARY

The forum discussion focuses on finding the leading order term in the asymptotic expansion of the integral ∫_0^1 sin(x(t - sin t)) sinh^a t dt. The method of stationary phase is employed, identifying t = 0 as the stationary point. Complications arise due to sinh(0) = 0, leading to a zero contribution from f(a). The discussion suggests integrating by parts to derive a new integrand that may yield a non-zero value at the stationary point.

PREREQUISITES
  • Understanding of Fourier integrals and the method of stationary phase
  • Familiarity with Taylor series expansion and derivatives
  • Knowledge of integration techniques, particularly integration by parts
  • Basic concepts of asymptotic analysis in mathematical physics
NEXT STEPS
  • Study the method of stationary phase in greater detail
  • Learn about asymptotic expansions in integrals
  • Explore advanced integration techniques, including integration by parts
  • Investigate the implications of stationary points in Fourier analysis
USEFUL FOR

Mathematics students, researchers in applied mathematics, and physicists dealing with asymptotic analysis and Fourier integrals.

Elec Man
Messages
1
Reaction score
0

Homework Statement



Find the leading order term in the expansion of

[tex]\int_0^1\ \sin(x(t-\sin t))\ \sinh^a t\ dt[/tex]


Homework Equations



For a general Fourier integral, the method of stationary phase gives

[tex]\int_a^b\ f(t) e^{i x \psi(t)}\ dt ~ f(a)e^{ix\psi(a) \pm i\pi/2p}(p!/(x|\psi^{(p)}(a)))^(1/p) \Gamma(1/p)/p[/tex]

where p is the order of the first nonzero derivative of [itex]\psi[/itex]

The Attempt at a Solution



The first thing to do is to put it into an easier form:

[tex]\int_0^1\ e^{ix(t-\sin t)} \sinh^a t\ dt[/tex]

Now it seems like it'd just be an ordinary Fourier integral. To evaluate a Fourier integral, you search for the point of stationary phase (where the derivative of the imaginary part is zero). That point is clearly t = 0. From here, you'd normally expand around t=0 and integrate. Unfortunately, it has two complications that I don't know how to deal with. The first is that [itex]\sinh 0 = 0[/itex], so in the above formula, f(a) = 0. The other problem is that you have to go out to the third derivative of (t-\sin t) to find a nonzero value, which makes brute forcing via Taylor expansion harder.

If you expand around t = 0, the integral becomes

[tex]\int_0^1\ t^a e^{ix(t^3/6)}[/tex]

which doesn't seem any easier to integrate. If you do integration by parts, you get another integral which increases with x, which is the opposite of what we want. Also, in this method, you usually expand the region of integration to 0 to ∞ on the basis that the largest contribution to the integral comes from the area near the stationary point (everything else oscillates rapidly and cancels). I'm not even sure that works here; there's a warning about that if f(a) = 0.

Does anyone have any experience with this that can help me out? Thanks!
 
Physics news on Phys.org
Since your f(0)=0 at the stationary point, what you can do is integrate by part
∫ f(t)exp(ix g(t))dt = f(t)/(ix g'(t))-1/(ix)*∫ (f(t)/g'(t))' exp(ix*g(t)) dt
and hopefully in the new integrand, (f(t)/g'(t))' is non-zero at the stationary point. Also note your stationary point contribution is only half, since your range of integral starts at 0. Also account for the endpoint contribution from 1.
 

Similar threads

Replies
2
Views
2K
Replies
2
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 23 ·
Replies
23
Views
2K
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 105 ·
4
Replies
105
Views
11K
Replies
12
Views
2K
  • · Replies 2 ·
Replies
2
Views
4K
Replies
4
Views
2K