|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| In [[mathematics]], the '''stationary phase approximation''' is a basic principle of [[asymptotic analysis]], applying to '''oscillatory integrals'''
| | Emilia Shryock is my title but you can contact me something you like. Her husband and her reside in Puerto Rico but she will have to move one working day or an additional. For many years he's been operating as a receptionist. Playing baseball is the pastime he will by no means quit doing.<br><br>my homepage [http://Tomjones1190.livejournal.com/ Tomjones1190.livejournal.com] |
| | |
| :<math>I(k) = \int g(x) e^{i k f(x)} \, dx</math>
| |
| | |
| taken over ''n''-dimensional space ℝ<sup>''n''</sup> where ''i'' is the [[imaginary unit]]. Here ''f'' and ''g'' are real-valued [[smooth function]]s. The role of ''g'' is to ensure convergence; that is, ''g'' is a [[test function]]. The large real parameter ''k'' is considered in the limit as {{nowrap|''k'' → ∞}}.
| |
| | |
| This method originates from the 19th century, and is due to [[George Gabriel Stokes]] and [[Lord Kelvin]].<ref>{{Citation
| |
| | last1 = Courant
| |
| | first1 = Richard
| |
| | author1-link = Richard Courant
| |
| | last2 = Hilbert
| |
| | first2 = David
| |
| | author2-link = David Hilbert
| |
| | title = [[Methoden der mathematischen Physik|Methods of mathematical physics]]
| |
| | volume = 1
| |
| | page = 474
| |
| | year = 1953
| |
| | edition = 2nd revised
| |
| | publisher = Interscience Publishers
| |
| | location = New York
| |
| | oclc = 505700
| |
| }}</ref>
| |
| | |
| ==Basics==
| |
| | |
| The main idea of stationary phase methods relies on the cancellation of sinusoids with rapidly varying phase. If many sinusoids have the same phase and they are added together, they will add constructively. If, however, these same sinusoids have phases which change rapidly as the frequency changes, they will add incoherently, varying between constructive and destructive addition at different times.
| |
| | |
| ==An example==
| |
| | |
| Consider a function
| |
| | |
| :<math>f(x,t) = \frac{1}{2\pi} \int_{\mathbb R} F(\omega) e^{i [k(\omega) x - \omega t]} \, d\omega</math>.
| |
| | |
| The phase term in this function, {{nowrap|''ϕ'' {{=}} ''k''(''ω'') ''x'' − ''ω'' ''t''}}, is stationary when
| |
| | |
| :<math>\frac{d}{d\omega}\mathopen{}\left(k(\omega) x - \omega t\right)\mathclose{} = 0</math>
| |
| | |
| or equivalently, | |
| | |
| :<math>\frac{d k}{d\omega} = \frac{t}{x}</math>.
| |
| | |
| Solutions to this equation yield dominant frequencies ''ω''<sub>0</sub> for some ''x'' and ''t''. If we expand ''ϕ'' as a [[Taylor series]] about ''ω''<sub>0</sub> and neglect terms of order higher than {{nowrap|(''ω'' − ''ω''<sub>0</sub>)<sup>2</sup>}},
| |
| | |
| :<math>\phi = \left[k(\omega_0) x - \omega_0 t\right] + \frac{1}{2} x k''(\omega_0) (\omega - \omega_0)^2 + \cdots</math>
| |
| | |
| where ''k''″ denotes the second derivative of ''k''. When ''x'' is relatively large, even a small difference {{nowrap|(''ω'' − ''ω''<sub>0</sub>)}} will generate rapid oscillations within the integral, leading to cancellation. Therefore we can extend the limits of integration beyond the limit for a Taylor expansion. If we double the real contribution from the positive frequencies of the transform to account for the negative frequencies,
| |
| | |
| :<math>f(x, t) \approx \frac{1}{2\pi} \cdot 2 \operatorname{Re} \left\{ e^{i \left[k(\omega_0) x - \omega_0 t\right]} \left|F(\omega_0)\right| \int_{\mathbb R} e^{\frac{1}{2} i x k''(\omega_0) (\omega - \omega_0)^2} \, d\omega \right\}</math>.
| |
| | |
| This integrates to
| |
| | |
| :<math>f(x, t) \approx \frac{\left|F(\omega_0)\right|}{\pi} \sqrt{\frac{2\pi}{x \left|k''(\omega_0)\right|}} \cos\left[k(\omega_0) x - \omega_0 t \pm \frac{\pi}{4}\right]</math>.
| |
| | |
| ==Reduction steps==
| |
| | |
| The first major general statement of the principle involved is that the asymptotic behaviour of ''I''(''k'') depends only on the [[critical point (mathematics)|critical point]]s of ''f''. If by choice of ''g'' the integral is localised to a region of space where ''f'' has no critical point, the resulting integral tends to 0 as the frequency of oscillations is taken to infinity. See for example [[Riemann-Lebesgue lemma]].
| |
| | |
| The second statement is that when ''f'' is a [[Morse function]], so that the singular points of ''f'' are [[non-degenerate critical point|non-degenerate]] and isolated, then the question can be reduced to the case ''n'' = 1. In fact, then, a choice of ''g'' can be made to split the integral into cases with just one critical point ''P'' in each. At that point, because the [[Hessian determinant]] at ''P'' is by assumption not 0, the [[Morse lemma]] applies. By a change of co-ordinates ''f'' may be replaced by
| |
| | |
| :<math>(x_1^2 + x_2^2 + \cdots + x_j^2) - (x_{j + 1}^2 + x_{j + 2}^2 + \cdots + x_n^2)</math>.
| |
| | |
| The value of ''j'' is given by the [[signature of a quadratic form|signature]] of the [[Hessian matrix]] of ''f'' at ''P''. As for ''g'', the essential case is that ''g'' is a product of [[bump function]]s of ''x''<sub>''i''</sub>. Assuming now without loss of generality that ''P'' is the origin, take a smooth bump function ''h'' with value 1 on the interval {{nowrap|[−1, 1]}} and quickly tending to 0 outside it. Take
| |
| | |
| :<math>g(x) = \prod_i h(x_i)</math>,
| |
| | |
| then [[Fubini's theorem]] reduces ''I''(''k'') to a product of integrals over the real line like
| |
| | |
| :<math>J(k) = \int h(x) e^{i k f(x)} \, dx</math>
| |
| | |
| with ''f''(''x'') = ±''x''<sup>2</sup>. The case with the minus sign is the [[complex conjugate]] of the case with the plus sign, so there is essentially one required asymptotic estimate.
| |
| | |
| In this way asymptotics can be found for oscillatory integrals for Morse functions. The degenerate case requires further techniques. See for example [[Airy function]].
| |
| | |
| ==One-dimensional case==
| |
| | |
| The essential statement is this one:
| |
| | |
| :<math>\int_{-1}^1 e^{i k x^2} \, dx = \sqrt{\frac{\pi}{k}} e^{i \pi / 4} + \mathcal O \mathopen{}\left(\frac{1}{k}\right)\mathclose{}</math>.
| |
| | |
| In fact by [[contour integration]] it can be shown that the main term on the right hand side of the equation is the value of the integral on the left hand side, extended over the range {{nowrap|[−∞, ∞]}}. Therefore it is the question of estimating away the integral over, say, {{nowrap|[1, ∞]}}.<ref>See for example [[Jean Dieudonné]], ''Infinitesimal Calculus'', p. 119.</ref>
| |
| | |
| This is the model for all one-dimensional integrals ''I''(''k'') with ''f'' having a single non-degenerate critical point at which ''f'' has [[second derivative]] > 0. In fact the model case has second derivative 2 at 0. In order to scale using ''k'', observe that replacing ''k'' by {{nowrap|''c'' ''k''}}
| |
| where ''c'' is constant is the same as scaling ''x'' by √''c''. It follows that for general values of {{nowrap|''f''″(0) > 0}}, the factor {{nowrap|√(''π'' / ''k'')}} becomes
| |
| | |
| :<math>\sqrt{\frac{2 \pi}{k f''(0)}}</math>.
| |
| | |
| For {{nowrap|''f''″(0) < 0}} one uses the complex conjugate formula, as mentioned before.
| |
| | |
| ==See also==
| |
| * [[Method of steepest descent]]
| |
| * [[Common integrals in quantum field theory]]
| |
| | |
| ==References==
| |
| * Bleistein, N. and Handelsman, R. (1975), ''Asymptotic Expansions of Integrals'', Dover, New York.
| |
| * Victor Guillemin and Shlomo Sternberg (1990), [http://www.ams.org/online_bks/surv14/ ''Geometric Asymptotics''], (see Chapter 1).
| |
| * {{citation|first=L.|last=Hörmander|authorlink=Lars Hörmander|title=Linear Partial Differential Operators, Volume 1|publisher=Springer-Verlag|year=1976|isbn=978-3-540-00662-6}}.
| |
| * Aki, Keiiti; & Richards, Paul G. (2002). "Quantitative Seismology" (2nd ed.), pp 255–256. University Science Books, ISBN 0-935702-96-2
| |
| *Wong, R. (2001), ''Asymptotic Approximations of Integrals'', Classics in Applied Mathematics, Vol. 34. Corrected reprint of the 1989 original. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA. xviii+543 pages, ISBN 0-89871-497-4.
| |
| | |
| ==Notes==
| |
| {{Reflist}}
| |
| | |
| ==External links==
| |
| *{{Springer|id=S/s087270|title=Stationary phase, method of the}}
| |
| | |
| [[Category:Mathematical analysis]]
| |
| [[Category:Perturbation theory]]
| |
Emilia Shryock is my title but you can contact me something you like. Her husband and her reside in Puerto Rico but she will have to move one working day or an additional. For many years he's been operating as a receptionist. Playing baseball is the pastime he will by no means quit doing.
my homepage Tomjones1190.livejournal.com