MCP๋กœ ์—ฐ๊ฒฐ โ†’

๊ณ„์‚ฐ ์ž…๋ ฅ

n = -1, 0, 1, 2, 3, 4
์ธ์ˆ˜ (์‹ค์ˆ˜, a > 0์—์„œ ์œ ํšจ)

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

Polygamma value ฮจn(a)
0.221322955737
at n = 1, a = 5.0
๋ฏธ๋ถ„ ์ฐจ์ˆ˜ n 1
์ธ์ˆ˜ a 5.0
ํ•จ์ˆ˜ Trigamma ฮจ1(a)

ํด๋ฆฌ๊ฐ๋งˆ ํ•จ์ˆ˜๋ž€?

์ฐจ์ˆ˜ n์˜ ํด๋ฆฌ๊ฐ๋งˆ ํ•จ์ˆ˜ \(\psi^{(n)}(a)\)๋Š” ๊ฐ๋งˆ ํ•จ์ˆ˜์˜ ์ž์—ฐ๋กœ๊ทธ๋ฅผ (n+1)๋ฒˆ ๋ฏธ๋ถ„ํ•œ ๊ฐ’์ž…๋‹ˆ๋‹ค. ์—ฌ๊ธฐ์„œ ์‚ฌ์šฉํ•˜๋Š” ์ธ๋ฑ์‹ฑ ๋ฐฉ์‹์—์„œ๋Š” \(n = -1\)์ด๋ฉด ๋กœ๊ทธ๊ฐ๋งˆ ํ•จ์ˆ˜ \(\ln\Gamma(a)\) ์ž์ฒด๋ฅผ, \(n = 0\)์ด๋ฉด ๋””๊ฐ๋งˆ ํ•จ์ˆ˜ \(\psi(a)\)๋ฅผ, \(n = 1\)์ด๋ฉด ํŠธ๋ฆฌ๊ฐ๋งˆ ํ•จ์ˆ˜๋ฅผ, \(n = 2\)์ด๋ฉด ํ…ŒํŠธ๋ผ๊ฐ๋งˆ ํ•จ์ˆ˜๋ฅผ ๋ฐ˜ํ™˜ํ•˜๋ฉฐ ์ด๋Ÿฐ ์‹์œผ๋กœ ์ด์–ด์ง‘๋‹ˆ๋‹ค. ์ด๋Ÿฌํ•œ ํŠน์ˆ˜ ํ•จ์ˆ˜๋“ค์€ ํ™•๋ฅ , ํ†ต๊ณ„(๊ฐ๋งˆ ๋ถ„ํฌ์™€ ๋ฒ ํƒ€ ๋ถ„ํฌ์˜ ์ตœ๋Œ€์šฐ๋„์ถ”์ •), ์ •์ˆ˜๋ก , ์ ๊ทผ ํ•ด์„ ์ „๋ฐ˜์— ๊ฑธ์ณ ๋“ฑ์žฅํ•ฉ๋‹ˆ๋‹ค. ์ˆœ์ˆ˜ ์ˆ˜ํ•™์— ํ•ด๋‹นํ•˜๋ฏ€๋กœ ์–ด๋””์—์„œ๋‚˜ ๋™์ผํ•˜๊ฒŒ ์ ์šฉ๋ฉ๋‹ˆ๋‹ค.

๋กœ๊ทธ๊ฐ๋งˆ์™€ ๊ทธ ์—ฐ์† ๋„ํ•จ์ˆ˜์ธ ํด๋ฆฌ๊ฐ๋งˆ ํ•จ์ˆ˜๋ฅผ ๋ณด์—ฌ์ฃผ๋Š” ๊ณก์„ 
ํด๋ฆฌ๊ฐ๋งˆ ํ•จ์ˆ˜๋Š” ๋กœ๊ทธ๊ฐ๋งˆ ํ•จ์ˆ˜์˜ ์—ฐ์†์ ์ธ ๋„ํ•จ์ˆ˜๋กœ ๋‚˜ํƒ€๋‚ฉ๋‹ˆ๋‹ค.

๊ณ„์‚ฐ๊ธฐ ์‚ฌ์šฉ ๋ฐฉ๋ฒ•

๋“œ๋กญ๋‹ค์šด์—์„œ ๋ฏธ๋ถ„ ์ฐจ์ˆ˜ n(-1, 0, 1, 2, 3, 4 ์ค‘ ํ•˜๋‚˜)์„ ์„ ํƒํ•˜๊ณ  ์ธ์ˆ˜ a๋ฅผ ์‹ค์ˆ˜๋กœ ์ž…๋ ฅํ•˜์„ธ์š”. ๊ณ„์‚ฐ ๋ฒ„ํŠผ์„ ๋ˆ„๋ฅด๋ฉด \(\psi^{(n)}(a)\)๋ฅผ ์•ฝ 12์ž๋ฆฌ ์œ ํšจ์ˆซ์ž๊นŒ์ง€ ๊ตฌํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค. ์ด ๊ณ„์‚ฐ ๋ฐฉ์‹์€ \(a > 0\)์—์„œ ์™„์ „ํžˆ ์œ ํšจํ•˜๋ฉฐ, a๊ฐ€ 0 ์ดํ•˜์ธ ๊ฒฝ์šฐ์—๋Š” ๋ฐ˜์‚ฌ ๊ณต์‹์ด ํ•„์š”ํ•ฉ๋‹ˆ๋‹ค. ๋˜ํ•œ 0 ์ดํ•˜์˜ ์ •์ˆ˜(๊ฐ๋งˆ ํ•จ์ˆ˜์˜ ๊ทน์ )๋Š” ์ •์˜๋˜์ง€ ์•Š์Œ์œผ๋กœ ํ‘œ์‹œ๋ฉ๋‹ˆ๋‹ค.

๊ณต์‹ ์„ค๋ช…

\(n = -1\)์ผ ๋•Œ๋Š” \(\ln\Gamma\)๋ฅผ ์œ„ํ•œ ๋ž€์ดˆ์Šค(Lanczos) ๊ทผ์‚ฌ๋ฅผ ์‚ฌ์šฉํ•ด ๊ฐ’์„ ๊ณ„์‚ฐํ•ฉ๋‹ˆ๋‹ค. $$\psi^{(-1)}(a) = \ln\Gamma\!\left(a\right)$$ \(n = 0\)์ผ ๋•Œ ๋””๊ฐ๋งˆ๋Š” ์ ํ™”์‹ \(\psi(x) = \psi(x+1) - \frac{1}{x}\)๋กœ ์ธ์ˆ˜๋ฅผ 6๋ณด๋‹ค ํฌ๊ฒŒ ๋งŒ๋“  ๋’ค, $$\psi^{(0)}(a) = \psi\!\left(a\right) = \frac{d}{da}\ln\Gamma\!\left(a\right) = \frac{\Gamma^{\prime}\!\left(a\right)}{\Gamma\!\left(a\right)}$$ ์ ๊ทผ ๊ธ‰์ˆ˜ \(\psi(x) \sim \ln x - \frac{1}{2x} - \frac{1}{12 x^2} + \cdots\) ๋กœ ๊ตฌํ•ฉ๋‹ˆ๋‹ค. \(n \geq 1\)์ผ ๋•Œ๋Š” ์ ํ™”์‹ \(\psi_n(x) = \psi_n(x+1) + (-1)^{n+1} \frac{n!}{x^{n+1}}\)๋กœ ์ธ์ˆ˜๋ฅผ 10 ์ด์ƒ์œผ๋กœ ์ด๋™์‹œํ‚จ ๋‹ค์Œ, ๋ฒ ๋ฅด๋ˆ„์ด ์ˆ˜ ๊ธฐ๋ฐ˜์˜ ์ ๊ทผ ๊ธ‰์ˆ˜๋กœ ๊ณ„์‚ฐ์„ ๋งˆ๋ฌด๋ฆฌํ•ฉ๋‹ˆ๋‹ค. $$\psi^{(n)}(a) = \frac{d^{n+1}}{da^{n+1}} \ln\Gamma\!\left(a\right), \quad n = \text{order}$$

๋กœ๊ทธ๊ฐ๋งˆ๋ฅผ ์—ฐ์†์ ์œผ๋กœ ๋ฏธ๋ถ„ํ•ด ํด๋ฆฌ๊ฐ๋งˆ ์ฐจ์ˆ˜๋ฅผ ๋งŒ๋“œ๋Š” ๋„์‹
๊ฐ ํด๋ฆฌ๊ฐ๋งˆ ์ฐจ์ˆ˜ n์€ ln ฮ“(a)๋ฅผ n+1๋ฒˆ ๋” ๋ฏธ๋ถ„ํ•˜์—ฌ ์–ป์–ด์ง‘๋‹ˆ๋‹ค.

๊ณ„์‚ฐ ์˜ˆ์ œ

\(n = 1\), \(a = 5\)(5์—์„œ์˜ ํŠธ๋ฆฌ๊ฐ๋งˆ)์ธ ๊ฒฝ์šฐ๋ฅผ ์‚ดํŽด๋ด…์‹œ๋‹ค. ํ•ญ๋“ฑ์‹ \(\psi_1(m) = \frac{\pi^2}{6} - \sum_{k=1}^{m-1} \frac{1}{k^2}\)๋ฅผ ์‚ฌ์šฉํ•˜๋ฉด $$1.6449340668 - \left(1 + \frac{1}{4} + \frac{1}{9} + \frac{1}{16}\right) = 1.6449340668 - 1.4236111111 = 0.2213229557$$ ์ด ๋ฉ๋‹ˆ๋‹ค. ๊ณ„์‚ฐ๊ธฐ๋Š” ์•ฝ \(0.221322955737\)์„ ๋ฐ˜ํ™˜ํ•ฉ๋‹ˆ๋‹ค.

์ž์ฃผ ๋ฌป๋Š” ์งˆ๋ฌธ

์™œ a = 0์ด๊ฑฐ๋‚˜ ์Œ์˜ ์ •์ˆ˜๋ฉด ์ •์˜๋˜์ง€ ์•Š๋‚˜์š”? ๊ฐ๋งˆ ํ•จ์ˆ˜๋Š” 0, -1, -2, ...์—์„œ ๊ทน์ ์„ ๊ฐ€์ง€๋ฏ€๋กœ \(\ln\Gamma\)์™€ ๋ชจ๋“  ํด๋ฆฌ๊ฐ๋งˆ ํ•จ์ˆ˜๊ฐ€ ๊ทธ ์ง€์ ์—์„œ ๋ฐœ์‚ฐํ•ฉ๋‹ˆ๋‹ค.

๋””๊ฐ๋งˆ์™€ ํŠธ๋ฆฌ๊ฐ๋งˆ์˜ ์ฐจ์ด๋Š” ๋ฌด์—‡์ธ๊ฐ€์š”? ๋””๊ฐ๋งˆ(\(n = 0\))๋Š” \(\ln\Gamma\)์˜ 1์ฐจ ๋„ํ•จ์ˆ˜์ด๊ณ , ํŠธ๋ฆฌ๊ฐ๋งˆ(\(n = 1\))๋Š” 2์ฐจ ๋„ํ•จ์ˆ˜๋กœ ๋””๊ฐ๋งˆ๋ฅผ ํ•œ ๋ฒˆ ๋” ๋ฏธ๋ถ„ํ•œ ๊ฒƒ๊ณผ ๊ฐ™์Šต๋‹ˆ๋‹ค.

๊ฒฐ๊ณผ๋Š” ์–ผ๋งˆ๋‚˜ ์ •ํ™•ํ•œ๊ฐ€์š”? ์ด๋™-์ ๊ทผ(shift-and-asymptotic) ๋ฐฉ๋ฒ•์„ ๋ฐฐ์ •๋ฐ€๋„(double-precision) ์—ฐ์‚ฐ์œผ๋กœ ์ ์šฉํ•˜๋ฉด \(a > 0\)์—์„œ ๋Œ€๋žต 12~15์ž๋ฆฌ์˜ ์œ ํšจ์ˆซ์ž๋ฅผ ์ •ํ™•ํ•˜๊ฒŒ ์–ป์„ ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

์ตœ์ข… ์—…๋ฐ์ดํŠธ: