1 minute read

์ˆ˜ํ•™๊ณผ ๋ณต์ˆ˜์ „๊ณต์„ ์œ„ํ•ด ์กธ์—… ๋งˆ์ง€๋ง‰ ํ•™๊ธฐ์— โ€œ์ˆ˜์น˜ํ•ด์„๊ฐœ๋ก โ€ ์ˆ˜์—…์„ ๋“ฃ๊ฒŒ ๋˜์—ˆ์Šต๋‹ˆ๋‹ค. ์ˆ˜ํ•™๊ณผ ์กธ์—…์‹œํ—˜๋„ ๊ฒธ์‚ฌ๊ฒธ์‚ฌ ์ค€๋น„ํ•  ๊ฒธ ํ™”์ดํŒ… ํ•ด๋ด…์‹œ๋‹ค!! ์ „์ฒด ํฌ์ŠคํŠธ๋Š” โ€œNumerical Analysisโ€œ์—์„œ ํ™•์ธํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

โ€œ๋ผ๊ทธ๋ž‘์ฃผ ๋ณด๊ฐ„๋ฒ•โ€œ์˜ ์‘์šฉ์— ๋Œ€ํ•ด ๋‹ค๋ฃจ๋Š” ํฌ์ŠคํŠธ ์ž…๋‹ˆ๋‹ค.

Nodal Polynomial

์ฃผ์–ด์ง„ $n$๊ฐœ์˜ ๋ฐ์ดํ„ฐ ๋…ธ๋“œ์— ๋Œ€ํ•ด ์•„๋ž˜์™€ ๊ฐ™์€ $n$์ฐจ ๋‹คํ•ญ์‹์„ ์ •์˜ํ•ฉ๋‹ˆ๋‹ค.

\[\omega_n(x) = \prod_{j=1}^n (x - x_j)\]

์ด๊ฒƒ์„ โ€œNodal Polynomialโ€์ด๋ผ๊ณ  ํ•ฉ๋‹ˆ๋‹ค. ์ด ๋‹คํ•ญ์‹์€ ๋ชจ๋“  $x_j$๋ฅผ ๊ทผ์œผ๋กœ ๊ฐ–์Šต๋‹ˆ๋‹ค.


๊ทธ๋ฆฌ๊ณ  ์ด ๋‹คํ•ญ์‹์„ ๋ฏธ๋ถ„ํ•˜๋ฉด ์•„๋ž˜์™€ ๊ฐ™์Šต๋‹ˆ๋‹ค.

\[\omega_n'(x) = \sum_{k=1}^{n} \left[ \prod_{j=1, j \ne k} (x - x_j) \right]\]

์‚ฌ์‹ค ๋‹จ์ˆœํžˆ โ€œ๊ณฑ์˜ ๋ฏธ๋ถ„๋ฒ•โ€์„ ์ ์šฉํ•œ ๊ฒฐ๊ณผ์ž…๋‹ˆ๋‹ค. ์ด ๋„ํ•จ์ˆ˜์— $x_j$๋ฅผ ๋Œ€์ž…ํ•˜๋ฉด

\[\omega_n'(x_j) = \prod_{i=1, i \ne j} (x_j - x_i)\]

์ด๋ ‡๊ฒŒ $x_j$๋ฅผ ์ œ์™ธํ•˜๊ณ  ๋งŒ๋“  Nodal Polynomial์„ ์–ป์„ ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

Lagrange Polynomial

์‹ ๊ธฐํ•˜๊ฒŒ๋„ โ€œ๋ผ๊ทธ๋ž‘์ฃผ ๋ณด๊ฐ„๋ฒ•โ€œ๋ฅผ ๊ตฌ์„ฑํ•˜๋Š” ๋ผ๊ทธ๋ž‘์ฃผ ๋‹คํ•ญ์‹ $L_i(x)$์„ Nodal Polynomial๋กœ ํ‘œํ˜„ํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค!

\[L_i(x) = \frac{\omega_n(x)}{(x-x_i) \cdot \omega_n'(x_i)}\]

์ด ํ‘œํ˜„์ด ์ข‹์€ ์ด์œ ๋Š” ๋ชจ๋“  $L_i(x)$์— ๋Œ€ํ•ด์„œ $\omega_n(x)$ ๋ถ€๋ถ„์ด ๊ณตํ†ต์ด๊ธฐ ๋•Œ๋ฌธ์ž…๋‹ˆ๋‹ค! ๊ทธ๋ž˜์„œ ๋ถ„์ž ํŒŒํŠธ์— ๋Œ€ํ•œ ๋ฐ˜๋ณต ๊ณ„์‚ฐ์„ ์ค„์ผ ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

๊ทธ๋ฆฌ๊ณ  $\omega_nโ€™(x_i)$๋Š” ์ƒ์ˆ˜ ์ž…๋‹ˆ๋‹ค! ๊ทธ๋ž˜์„œ ์ฒ˜์Œ ๊ตฌ์„ฑํ•  ๋•Œ ํ•œ๋ฒˆ๋งŒ ๊ณ„์‚ฐํ•˜๋ฉด ๋ฉ๋‹ˆ๋‹ค!

์ •๊ทœํ™” ์ƒ์ˆ˜์˜ ์˜๋ฏธ๋Š” ๋ญ”๊ฐ€?

alt text

alt text

๋‰ดํ„ด ๋ฐฉ๋ฒ•์—์„œ๋„ ์ด ๋ฐฉ์‹์ด ์–ธ๊ธ‰ ๋จ.