정지 위상 근사

1. 개요

정지 위상 근사는 유클리드 공간에서의 적분을 근사하는 방법으로, 매끄러운 함수와 큰 실수 매개변수 k를 포함하는 적분 I(k)의 값을 k가 매우 클 때 근사하는 데 사용된다. 이 방법은 적분 내에서 빠르게 변동하는 위상을 가진 사인파의 상쇄 원리에 기반하며, 함수 f의 정류점과 헤세 행렬의 고윳값에 의존하는 공식을 제공한다. 정지 위상 근사는 1차원 경우와 다변수 경우로 확장되며, 특히 모스 함수와 같은 특수한 경우에 적용된다. 이 기법은 19세기에 조지 가브리엘 스토크스와 윌리엄 톰슨에 의해 처음 도입되었다.

2. 정의

유클리드 공간 $\backslash mathbb\; R^n$ 위에 다음과 같은 꼴의 적분이 있다고 하자.

:$I(k)=\backslash int\_\{\backslash mathbb\; R^n\}d^n\backslash mathbf\; x\backslash ,g(\backslash mathbf\; x)\backslash exp(ikf(\backslash mathbf\; x))$

여기서 $f$와 $g$는 매끄러운 실함수이며, $g$는 무한대에서 적분이 수렴하도록 충분히 빨리 0으로 간다고 하자. $k$는 임의의 양의 실수 매개변수이다.

$f(\backslash mathbf\; x)$가 $\backslash mathbf\; x=\backslash mathbf\; x\_0$에서 유일한 정류점을 가지고, 이 점에서의 헤세 행렬 $Hf(\backslash mathbf\; x\_0)$이 고윳값 0을 갖지 않는다고 하며, 양의 고윳값들이 $n\_+$개, 음의 고윳값들이 $n\_-$개라고 하자.

:$\backslash nabla\; f(\backslash mathbf\; x\_0)=0$

그렇다면 매우 큰 $k$에 대하여 적분 $I(k)$를 다음과 같이 근사할 수 있다. 이를 정지 위상 근사라고 한다.

:$I(k)\backslash approx\; g(\backslash mathbf\; x\_0)\backslash exp(ikf(\backslash mathbf\; x\_0)+i(n\_+-n\_-)\backslash pi/4)\backslash sqrt\{2\backslash pi/(k|\backslash det\; Hf(\backslash mathbf\; x\_0)|)\}+\backslash mathcal\; O(1/k)$

만약 여러 개의 고립된 정류점들이 존재한다면, 각 점들에서의 값들을 더하면 된다.

함수 $f$의 임계점 집합을 $\backslash Sigma$라고 하고, $g$가 콤팩트하게 지원되거나 지수적으로 감쇠하며, 모든 임계점이 비퇴화($x\_0\; \backslash in\; \backslash Sigma$에 대해 $\backslash det(\backslash mathrm\{Hess\}(f(x\_0)))\backslash neq\; 0$)라고 가정하면, $k\backslash to\; \backslash infty$일 때 다음 점근 공식을 얻는다.

:$\backslash int\_\{\backslash mathbb\{R\}^n\}g(x)e^\{ikf(x)\}\; dx=\backslash sum\_\{x\_0\backslash in\; \backslash Sigma\}\; e^\{ik\; f(x\_0)\}|\backslash det(\{\backslash mathrm\{Hess\}\}(f(x\_0)))|^\{-1/2\}e^\{\backslash frac\{i\backslash pi\}\{4\}\; \backslash mathrm\{sgn\}(\backslash mathrm\{Hess\}(f(x\_0)))\}(2\backslash pi/k)^\{n/2\}g(x\_0)+o(k^\{-n/2\})$

여기서 $\backslash mathrm\{Hess\}(f)$는 $f$의 헤시안을 나타내고, $\backslash mathrm\{sgn\}(\backslash mathrm\{Hess\}(f))$는 헤시안의 부호수, 즉 양의 고유값의 개수에서 음의 고유값의 개수를 뺀 값을 나타낸다.

$n=1$인 경우, 이는 다음과 같이 축소된다.

:$\backslash int\_\backslash mathbb\{R\}g(x)e^\{ikf(x)\}dx=\backslash sum\_\{x\_0\backslash in\; \backslash Sigma\}\; g(x\_0)e^\{ik\; f(x\_0)+\backslash mathrm\{sign\}(f$(x_0))i\pi/4}\left(\frac{2\pi}{k |f(x_0)|}\right)^{1/2}+o(k^{-1/2})

이 경우, $f$에 대한 가정은 모든 임계점이 비퇴화된다는 것으로 축소된다.

이것은 단지 윅 회전된 최급강하법 공식의 버전이다.

3. 기본 원리

정지 위상 기법의 주요 아이디어는 빠르게 변동하는 위상을 가진 사인파의 상쇄에 의존한다. 많은 사인파가 동일한 위상을 가지고 함께 더해지면, 그들은 보강적으로 더해진다. 그러나, 이러한 동일한 사인파들이 주파수가 변함에 따라 빠르게 변하는 위상을 갖는다면, 그들은 비일관적으로 더해져서, 다른 시간에서 보강적 가산과 파괴적 가산 사이에서 변동한다.

핵심적인 내용은 다음과 같다.

:$\backslash int\_\{-1\}^1\; e^\{i\; k\; x^2\}\; \backslash ,\; dx\; =\; \backslash sqrt\{\backslash frac\{\backslash pi\}\{k\}\}\; e^\{i\; \backslash pi\; /\; 4\}\; +\; \backslash mathcal\; O\; \backslash mathopen\{\}\backslash left(\backslash frac\{1\}\{k\}\backslash right)\backslash mathclose\{\}$.

실제로 경로 적분을 통해 방정식의 우변에 있는 주요 항은 적분 구간을 $[-\backslash infty,\; \backslash infty]$로 확장한 좌변의 적분 값과 같다는 것을 보일 수 있다(증명은 프레넬 적분 참조). 예를 들어 $[1,\backslash infty]$에 대한 적분을 추정하는 문제가 남는다.

이것은 $f$가 $f$의 이계 도함수가 $>0$인 비퇴화 임계점 하나를 갖는 모든 1차원 적분 $I(k)$에 대한 모델이다. 사실, 모델 사례는 0에서 이계 도함수 2를 갖는다. $k$를 사용하여 스케일링하기 위해, $k$를 상수 $c$로 $ck$로 대체하는 것은 $x$를 $\backslash sqrt\{c\}$로 스케일링하는 것과 같다는 것을 관찰한다. 따라서 $f$(0)>0의 일반적인 값에 대해, $\backslash sqrt\{\backslash pi/k\}$ 인수는 다음과 같이 된다.

:$\backslash sqrt\{\backslash frac\{2\; \backslash pi\}\{k\; f$(0)}}.

인 경우, 앞서 언급했듯이 복소 켤레 공식을 사용한다.

3.1. 1차원 정지 위상 근사

핵심적인 내용은 다음과 같다.

:$\backslash int\_\{-1\}^1\; e^\{i\; k\; x^2\}\; \backslash ,\; dx\; =\; \backslash sqrt\{\backslash frac\{\backslash pi\}\{k\}\}\; e^\{i\; \backslash pi\; /\; 4\}\; +\; \backslash mathcal\; O\; \backslash mathopen\{\}\backslash left(\backslash frac\{1\}\{k\}\backslash right)\backslash mathclose\{\}$.

실제로 경로 적분을 통해 방정식의 우변에 있는 주요 항은 적분 구간을 $[-\backslash infty,\; \backslash infty]$로 확장한 좌변의 적분 값과 같다는 것을 보일 수 있다(증명은 프레넬 적분 참조). 예를 들어 $[1,\backslash infty]$에 대한 적분을 추정하는 문제가 남는다.

이것은 $f$가 $f$의 이계 도함수가 $>0$인 비퇴화 임계점 하나를 갖는 모든 1차원 적분 $I(k)$에 대한 모델이다. 사실, 모델 사례는 0에서 이계 도함수 2를 갖는다. $k$를 사용하여 스케일링하기 위해, $k$를 상수 $c$로 $ck$로 대체하는 것은 $x$를 $\backslash sqrt\{c\}$로 스케일링하는 것과 같다는 것을 관찰한다. 따라서 $f$(0)>0의 일반적인 값에 대해, $\backslash sqrt\{\backslash pi/k\}$ 인수는 다음과 같이 된다.

:$\backslash sqrt\{\backslash frac\{2\; \backslash pi\}\{k\; f$(0)}}.

인 경우, 앞서 언급했듯이 복소 켤레 공식을 사용한다.

4. 유도

좌표 변환을 통해, 다변수 정지 위상 근사는 1변수의 경우로부터 유도할 수 있다. 1변수의 경우, 정지 위상 근사는 기본적으로 다음과 같은 꼴이다. $|k|$가 매우 크다면, $a$에 상관없이 다음이 성립한다.
:$\backslash int\_\{-a\}^adx\backslash ,\backslash exp(ikx^2/2)=\backslash sqrt\{2\backslash pi/|k|\}\backslash exp((\backslash operatorname\{sgn\}k)i\backslash pi/4)+\backslash mathcal\; O(1/k)$
여기서 $\backslash operatorname\{sgn\}k=k/|k|\backslash in\backslash \{+1,-1\backslash \}$는 $k$의 부호이다. 이 1변수 적분은 쉽게 유도할 수 있다.

I(k)의 점근적 거동이 f의 임계점에만 의존한다. g를 선택하여 적분을 f가 임계점을 갖지 않는 공간 영역으로 국한하면, 진동 주파수가 무한대로 갈 때 결과 적분은 0으로 수렴한다. 예를 들어, 리만-르베그 보조정리를 참조할 수 있다.

f가 모스 함수여서 f의 특이점이 비퇴화이고 고립되어 있을 때, 문제는 n = 1의 경우로 줄일 수 있다. 그러면 g를 선택하여 각 경우에 단 하나의 임계점 P만 갖는 경우로 적분을 나눌 수 있다. 그 지점에서 P에서의 헤세 행렬식이 가정에 의해 0이 아니기 때문에 모스 보조정리가 적용된다. 좌표 변환에 의해 f는 다음과 같이 대체될 수 있다.

:$(x\_1^2\; +\; x\_2^2\; +\; \backslash cdots\; +\; x\_j^2)\; -\; (x\_\{j\; +\; 1\}^2\; +\; x\_\{j\; +\; 2\}^2\; +\; \backslash cdots\; +\; x\_n^2)$.

j의 값은 P에서의 f의 헤세 행렬의 부호수에 의해 주어진다. g의 경우, 본질적인 경우는 g가 x_i의 범프 함수의 곱이라는 것이다. 이제 일반성을 잃지 않고 P가 원점이라고 가정하고, 구간 에서 값 1을 가지고 그 밖에서는 빠르게 0으로 수렴하는 부드러운 범프 함수 h를 취한다.

:$g(x)\; =\; \backslash prod\_i\; h(x\_i)$를 취하면 푸비니 정리에 의해 I(k)는 다음과 같은 실수선에 대한 적분의 곱으로 줄어듭니다.

:$J(k)\; =\; \backslash int\; h(x)\; e^\{i\; k\; f(x)\}\; \backslash ,\; dx$

f(x) = ±x²인 경우이다. 음수 부호가 있는 경우는 양수 부호가 있는 경우의 켤레 복소수이므로 본질적으로 하나의 필요한 점근적 추정이 있다.

이러한 방식으로 모스 함수에 대한 진동 적분에 대한 점근성을 찾을 수 있다. 퇴화된 경우는 추가적인 기술이 필요하다(예: 에어리 함수 참조).

5. 일반화

5.1. 모스 함수

5.2. 퇴화된 임계점

6. 예시

함수를 고려해 보자.

:$f(x,t)\; =\; \backslash frac\{1\}\{2\backslash pi\}\; \backslash int\_\{\backslash mathbb\; R\}\; F(\backslash omega)\; e^\{i\; [k(\backslash omega)\; x\; -\; \backslash omega\; t]\}\; \backslash ,\; d\backslash omega$.

이 함수의 위상 항 $\backslash phi\; =\; k(\backslash omega)\; x\; -\; \backslash omega\; t$는 다음 조건에서 정지한다.

:$\backslash frac\{d\}\{d\backslash omega\}\backslash mathopen\{\}\backslash left(k(\backslash omega)\; x\; -\; \backslash omega\; t\backslash right)\backslash mathclose\{\}\; =\; 0$

또는 다음과 동등하다.

:$\backslash frac\{d\; k(\backslash omega)\}\{d\backslash omega\}\backslash Big|\_\{\backslash omega\; =\; \backslash omega\_0\}\; =\; \backslash frac\{t\}\{x\}$.

이 방정식의 해는 어떤 $x$와 $t$에 대해 지배적인 주파수 $\backslash omega\_0$를 제공한다. $\backslash phi$를 $\backslash omega\_0$에 대한 테일러 급수로 전개하고 $(\backslash omega-\backslash omega\_0)^2$보다 높은 차수의 항을 무시하면 다음을 얻는다.

:$\backslash phi\; =\; \backslash left[k(\backslash omega\_0)\; x\; -\; \backslash omega\_0\; t\backslash right]\; +\; \backslash frac\{1\}\{2\}\; x\; k$(\omega_0) (\omega - \omega_0)^2 + \cdots

여기서 $k$는 $k$의 2차 도함수를 나타낸다. $x$가 비교적 크면 작은 차이 $(\backslash omega-\backslash omega\_0)$조차도 적분 내에서 빠른 진동을 생성하여 상쇄를 일으킨다. 따라서 테일러 전개의 한계를 넘어 적분 한계를 확장할 수 있다. 다음 공식을 사용하면,

:$\backslash int\_\{\backslash mathbb\; R\}\; e^\{\backslash frac\{1\}\{2\}ic\; x^2\}\; d\; x=\backslash sqrt\{\backslash frac\{2i\backslash pi\}\{c\}\}=\backslash sqrt\{\backslash frac\{2\backslash pi\}$