Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
| 수학:안장점_근사 [2026/06/12 14:14] – [예: 스털링 근사] admin | 수학:안장점_근사 [2026/06/15 14:18] (current) – [안장점 방법] admin | ||
|---|---|---|---|
| Line 1: | Line 1: | ||
| - | ======라플라스 방법====== | ||
| - | |||
| - | $$I = \int_a^b dt ~e^{-x f(t)} g(t)$$ | ||
| - | 의 적분을 구하고자 한다. | ||
| - | $x \gg 1$이고 구간 안의 어떤 $t_0$에서 $f$가 최소여서 $f' | ||
| - | $x$가 매우 크므로 $f$의 최소점 부근에서만 주로 적분의 기여가 있을 것으로 기대할 수 있다. | ||
| - | 이 최소점 주위에서 $f(t)$를 이차 함수로 전개하여 | ||
| - | $$f(t) \approx f(t_0) + \frac{1}{2} f'' | ||
| - | 로 적는다. 위 적분은 다음과 같은 근사값을 가진다: | ||
| - | \begin{eqnarray} | ||
| - | I & | ||
| - | & | ||
| - | \end{eqnarray} | ||
| - | 이 때 $t = t_0 + y\sqrt{2/ | ||
| - | $$I \approx e^{-x f(t_0)} g(t_0) \sqrt{\frac{2\pi}{x f'' | ||
| - | 을 얻는다. | ||
| - | |||
| - | =====예: 스털링 근사===== | ||
| - | 감마 함수의 정의에서 시작하자. $t=n u$로 변수변환을 시행할 것이고, $n$은 매우 큰 양수로 가정한다. | ||
| - | \begin{eqnarray*} | ||
| - | \Gamma(n+1) &=& \int_0^\infty t^n e^{-t} dt = \int_0^\infty \exp\left( n \ln t -t \right) dt\\ | ||
| - | &=& \int_0^\infty n^n \exp\left[ n(\ln u - u) \right] du\\ | ||
| - | &=& \int_0^\infty \exp\left[ n(\ln n + \ln u - u) \right] du\\ | ||
| - | &=& \int_0^\infty e^{-n\phi(u)} du. | ||
| - | \end{eqnarray*} | ||
| - | 이때 $\phi(u) = -\ln n - \ln u + u$로서, 최소가 되는 지점을 찾으면 $\phi' | ||
| - | $$\Gamma(n+1) \approx e^{-n\phi(1)} \sqrt{\frac{2\pi}{n\phi'' | ||
| - | |||
| ======안장점 방법====== | ======안장점 방법====== | ||
| 채플링은 다음과 같은 구체적인 예로 설명하고 있다. $x$는 매우 큰 양수이다. | 채플링은 다음과 같은 구체적인 예로 설명하고 있다. $x$는 매우 큰 양수이다. | ||
| $$I(x) = \int_{-\infty}^{\infty} e^{x \phi(t)} dt = \int_{-\infty}^{\infty} e^{ix(t^4-4t)} dt.$$ | $$I(x) = \int_{-\infty}^{\infty} e^{x \phi(t)} dt = \int_{-\infty}^{\infty} e^{ix(t^4-4t)} dt.$$ | ||
| - | 따라서 $\phi(z) \equiv iz^4 - 4iz$이고 $z=p+iq$로 실수부와 허수부를 표현하면, | + | 복소 적분을 활용하여 이 실수 적분의 값을 구하려고 한다. 이를 위해 복소평면상에서 적절한 경로(contour)를 택할 필요가 있다. |
| + | |||
| + | 먼저 | ||
| $$\phi = u + iv = \left( -4p^3 q + 4pq^3 + 4q \right) + i\left(p^4 + q^4 - 6p^2 q^2 - 4p\right)$$ | $$\phi = u + iv = \left( -4p^3 q + 4pq^3 + 4q \right) + i\left(p^4 + q^4 - 6p^2 q^2 - 4p\right)$$ | ||
| 로서 $u(p,q) = -4p^3 q + 4pq^3 + 4q$와 $v(p,q) = p^4 + q^4 - 6p^2 q^2 - 4p$이다. | 로서 $u(p,q) = -4p^3 q + 4pq^3 + 4q$와 $v(p,q) = p^4 + q^4 - 6p^2 q^2 - 4p$이다. | ||
| Line 55: | Line 29: | ||
| {{: | {{: | ||
| - | 경로의 양 끝이 어디를 향해야 하는지 알고 있으므로 그 조건을 만족하는 경로만을 고르면 아래와 같다. 이 경로는 $z_0=1$과 $z_2=e^{-2\pi i/3}$를 지난다. | + | 경로의 양 끝이 어디를 향해야 하는지 알고 있으므로 그 조건을 만족하는 경로만을 고르면 아래와 같다. 이 경로는 $z_0=1$과 $z_2=e^{-2\pi i/3}$를 지난다. 두 곡선이 가까워지는 부분은 피적분함수가 0으로 접근하는 영역이므로 적분 값에 영향을 주지 않고 두 선을 연결할 수 있다. |
| {{: | {{: | ||
| - | $z_0$와 $z_2$의 기여분을 비교해보면 ($\Re$는 실수부 의미) | + | 또 $z_0$와 $z_2$의 기여분을 비교해보면 ($\Re$는 실수부 의미) |
| $$\frac{\exp\left[x \Re{\phi(z_2)}\right]}{\exp\left[x \Re{\phi(z_0)}\right]} = \exp \left( -\frac{3\sqrt{3}}{2}x \right)$$ | $$\frac{\exp\left[x \Re{\phi(z_2)}\right]}{\exp\left[x \Re{\phi(z_0)}\right]} = \exp \left( -\frac{3\sqrt{3}}{2}x \right)$$ | ||
| 로서, $z_2$에서 오는 기여분이 지수함수적으로 작다. 따라서 $z_0=1$ 부근만을 살펴보기로 한다. | 로서, $z_2$에서 오는 기여분이 지수함수적으로 작다. 따라서 $z_0=1$ 부근만을 살펴보기로 한다. | ||
| - | $z_0$를 지나는 경로 $\gamma$를 실수 변수 $t$의 급수해 $F_n(t) = \sum_{k=0}^n c_k t^k$로 표현해보자. 계수 $c_k$는 복소수이다. 이를 따라갈 때에 $\phi = \phi(z_0) - t^2 + O(t^n)$ 형태가 되게끔 계수들을 결정하려고 한다. $c_0=z_0=1$임은 자명하다. 다음 식을 적고 | + | $z_0$를 지나는 경로 $\gamma$를 실수 변수 $t$의 급수 |
| \begin{eqnarray*} | \begin{eqnarray*} | ||
| \phi(1+c_1t+\ldots) &=& -3i + 6i c_1^2 t^2 + 4i(c_1^3+3c_2 c_1) t^3 + i(c_1^4 + 12 c_2 c_1^2 + 12 c_3 c_1 + 6 c_2^2) t^4 + \ldots\\ | \phi(1+c_1t+\ldots) &=& -3i + 6i c_1^2 t^2 + 4i(c_1^3+3c_2 c_1) t^3 + i(c_1^4 + 12 c_2 c_1^2 + 12 c_3 c_1 + 6 c_2^2) t^4 + \ldots\\ | ||
| Line 80: | Line 54: | ||
| &=& e^{-3ix + i\pi/4} \frac{\sqrt{\pi}}{\sqrt{6x}} \left( 1 + \frac{7i}{144x} \right). | &=& e^{-3ix + i\pi/4} \frac{\sqrt{\pi}}{\sqrt{6x}} \left( 1 + \frac{7i}{144x} \right). | ||
| \end{eqnarray*} | \end{eqnarray*} | ||
| - | ======함께 보기====== | ||
| - | [[물리: | ||
| - | [[물리: | + | 만일 큰 고민 없이 곧바로 [[수학: |
| + | $$I(x) \approx e^{-3ix} \frac{\sqrt{\pi}}{\sqrt{6x}}.$$ | ||
| + | ======함께 보기====== | ||
| + | * [[수학: | ||
| + | * [[수학: | ||
| + | * [[물리: | ||
| + | * [[물리: | ||
| ======참고문헌====== | ======참고문헌====== | ||
| - | * Richard Chapling, [[https:// | + | * Richard Chapling, [[https:// |