- Método de la fase estacionaria
-
En matemáticas, el método de la fase estacionaria o aproximación de fase estacionaria es un principio básico del análisis asintótico, se aplica a las integrales oscilatorias, una clase de integrales de Fourier del tipo.
definidas en el espacio n-dimensional Rn, donde i es la unidad imaginaria. Aquí f y g son funciones continuamente diferenciables que toman valores reales. El papel de la función g es para asegurar la convergencia, es decir, g es una función test. El gran parámetro real k se considera en el límite cuando k → ∞.
Contenido
Conceptos básicos
La idea principal de los métodos de la fase estacionaria se basa en la cancelación de sinusoides cuando la fase varía rápidamente. Si muchas sinusoides tienen la misma fase y son sumadas, se tendrá una suma constructiva. Si, sin embargo, estas mismos sinusoides tienen fases que cambian rápidamente a medida que cambia la frecuencia, se tendrá una suma destructiva.
Un ejemplo
Consideremos la función
El término de fase en esta función, ϕ = kx − ωt es "estacionario" cuando
o de modo equivalente,
Las soluciones de esta ecuación producen frecuencias dominantes ωdom(x,t) for a given x and t. Si expandimos ϕ en una serie de Taylor sobre ωdom y anulamos los términos de orden superior a (ω − ωdom)2,
Cuando x es relativamente grande, incluso una pequeña diferencia ω − ωdom generará oscilaciones rápidas dentro de la integral, lo que lleva a la cancelación. Por lo tanto podemos extender los límites de la integración más allá del límite para un desarrollo de Taylor. Si duplicamos la contribución real de las frecuencias positivas de la transformación para tener en cuenta las frecuencias negativas,
Esto se integra como
Etapas de reducción
La primera afirmación general importante del principio en cuestión establece que el comportamiento asintótico de I(k) depende sólo de los puntos críticos de f. Si por causa de g la integral se localiza en una región del espacio donde f no tiene un punto crítico, el resultado de la integral tiende a 0 cuando la frecuencia de las oscilaciones se toman hasta el infinito. Véase, por ejemplo, el lema de Riemann-Lebesgue.
La segunda afirmación es que, cuando f es una función de Morse, de modo que los puntos singulares de f sean puntos críticos no degenerados y aislados, entonces la cuestión se reduce al caso n = 1. De hecho, a continuación, la elección de g se puede hacer para dividir la integral en casos con un solo punto crítico de P en cada uno de ellos. En ese punto, porque la matriz hessiana en P es por suposición distinta de 0, se aplica el lema de Morse. Mediante un cambio de coordenadas f podrá ser sustituido por
- x12 + ….+ xj2 − xj + 12 − xj + 22 − … − xn2.
El valor de j está dado por la signatura de la matriz hessiana de f en P. En cuanto a g, el caso esencial es que g es un producto de funciones de activación (bump function) de xi. Si se asume ya sin pérdida de generalidad que P es el origen, tomamos una función de activación continuamente diferenciable h con valor 1 en el intervalo [-1,1] y que rápidamente tiende a 0 fuera de ella.
Tomemos
- g(x) = Π h(xi).
Por tanto el teorema de Fubini reduce I(k) a un producto de integrales sobre la línea real, como
con f(x) = x2 o −x2. El caso con el signo menos es el complejo conjugado del caso con el signo más, por lo que hay esencialmente una estimación asintótica necesaria.
De esta manera asintótica se pueden encontrar soluciones para las integrales oscilatorias de las funciones de Morse. El caso degenerado requiere técnicas adicionales. Véase, por ejemplo, la función de Airy.
Caso unidimensional
La definición esencial es esta:
De hecho por integración de contorno se puede demostrar que el término principal en el lado derecho de la ecuación es el valor de la integral en el lado izquierdo, extendida sobre el intervalo [- ∞, ∞]. Por eso, la cuestión es la estimación de la integral sobre, por ejemplo, [1, ∞].[1]
Este es el modelo para todas las integrales de una dimensión I(k) con f teniendo un único punto crítico no degenerado en el que f tiene segunda derivada mayor que 0. De hecho, el caso modelo tiene como segunda derivada 2 en 0. Con el fin de escalar con k, teniendo en cuenta que la sustitución de k por ck donde c es una constante es la misma que escalar x mediante √ c. De ello se deduce que para los valores generales de f "(0)> 0, el factor √ (π / K) se convierte en
Para f "(0) <0 se usa la fórmula conjugada compleja, como se mencionó antes.
Véase también
- Método del descenso más brusco
- Integrales frecuentes en teoría cuántica de campos
Referencias
- ↑ Ver por ejemplo Jean Dieudonné, Infinitesimal Calculus, p. 119.
- Bleistein, N. and Handelsman, R. (1975), Asymptotic Expansions of Integrals, Dover, New York.
- Victor Guillemin and Shlomo Sternberg (1990), Geometric Asymptotics, (see Chapter 1).
- Aki, Keiiti; & Richards, Paul G. (2002). "Quantitative Seismology" (2nd ed.), pp 255–256. University Science Books, ISBN 0-935702-96-2
Enlaces externos
- Hazewinkel, Michiel, ed. (2001), «Stationary phase, method of the» (en inglés), Encyclopaedia of Mathematics, Springer, ISBN 978-1556080104
Categorías:- Análisis matemático
- Teoría de perturbaciones
![f(x, t) = \frac{1}{2\pi} 2 \mbox{Re}\left\{ \exp\left[i\left[k(\omega_{dom})x-\omega_{dom}t\right]\right] \left|F(\omega_{dom})\right| \int_{\mathbb{R}}\exp\left[i\frac{x}{2}\frac{d^2k}{d\omega^2}(\omega-\omega_{dom})^2\right]d\omega\right\}](6/e86e9e87dec31fcba8ce63d29fcac539.png)
![f(x, t) \sim \frac{\left|F(\omega_{dom})\right|}{\pi} \sqrt{ \frac{2\pi}{x\left|\frac{d^2k}{d\omega^2}\right|}} \cos\left[ k(\omega_{dom})x - \omega_{dom}t \pm \frac{\pi}{4}\right]](6/e96a51c4902268e495d035ffcd9dcb18.png)
Wikimedia foundation. 2010.
