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

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

20-Point Gauss-Legendre Rule on [-1, 1]
2
๊ฐ€์ค‘์น˜ ํ•ฉ (2๊ฐ€ ๋˜์–ด์•ผ ํ•จ)
i ๋…ธ๋“œ xแตข ๊ฐ€์ค‘์น˜ wแตข
1 -0.993128599185095 0.017614007139152
2 -0.963971927277914 0.040601429800387
3 -0.912234428251326 0.062672048334109
4 -0.839116971822219 0.083276741576705
5 -0.746331906460151 0.10193011981724
6 -0.636053680726515 0.118194531961518
7 -0.510867001950827 0.131688638449176
8 -0.37370608871542 0.142096109318382
9 -0.227785851141645 0.149172986472604
10 -0.076526521133497 0.152753387130726
11 0.076526521133497 0.152753387130726
12 0.227785851141645 0.149172986472604
13 0.37370608871542 0.142096109318382
14 0.510867001950827 0.131688638449176
15 0.636053680726515 0.118194531961518
16 0.746331906460151 0.10193011981724
17 0.839116971822219 0.083276741576705
18 0.912234428251326 0.062672048334109
19 0.963971927277914 0.040601429800387
20 0.993128599185095 0.017614007139152

๊ฐ€์šฐ์Šค-๋ฅด์žฅ๋“œ๋ฅด ๊ตฌ์ ๋ฒ• ๊ณ„์‚ฐ๊ธฐ๋ž€?

์ด ๋„๊ตฌ๋Š” ๊ธฐ์ค€ ๊ตฌ๊ฐ„ [-1, 1]์—์„œ n์  ๊ฐ€์šฐ์Šค-๋ฅด์žฅ๋“œ๋ฅด ๊ตฌ์ ๋ฒ•์˜ ๋…ธ๋“œ(๋ถ„์ )์™€ ๊ฐ€์ค‘์น˜๋ฅผ ๊ณ„์‚ฐํ•ฉ๋‹ˆ๋‹ค. ๊ฐ€์šฐ์Šค-๋ฅด์žฅ๋“œ๋ฅด ๊ตฌ์ ๋ฒ•์€ ์ •์ ๋ถ„์„ ํ•จ์ˆ˜๊ฐ’์˜ ๊ฐ€์ค‘ํ•ฉ์œผ๋กœ ๊ทผ์‚ฌํ•˜๋Š” ์ˆ˜์น˜ ์ ๋ถ„ ๊ธฐ๋ฒ•์œผ๋กœ, -1๋ถ€ํ„ฐ 1๊นŒ์ง€์˜ \(\int f(x)\,dx\) ๊ฐ’์„ \(\sum w_i\, f(x_i)\) ํ˜•ํƒœ๋กœ ๋‚˜ํƒ€๋ƒ…๋‹ˆ๋‹ค. ๋‹จ n๊ฐœ์˜ ์ ๋งŒ์œผ๋กœ ์ฐจ์ˆ˜๊ฐ€ 2n-1 ์ดํ•˜์ธ ๋ชจ๋“  ๋‹คํ•ญ์‹์„ ์ •ํ™•ํ•˜๊ฒŒ ์ ๋ถ„ํ•˜๋ฏ€๋กœ, ์‚ฌ๋‹ค๋ฆฌ๊ผด ๊ณต์‹์ด๋‚˜ ์‹ฌํ”„์Šจ ๊ณต์‹์ฒ˜๋Ÿผ ๋“ฑ๊ฐ„๊ฒฉ ์ ์„ ์“ฐ๋Š” ๋ฐฉ๋ฒ•๋ณด๋‹ค ํ›จ์”ฌ ์ •๋ฐ€ํ•ฉ๋‹ˆ๋‹ค.

์‚ฌ์šฉ ๋ฐฉ๋ฒ•

๋จผ์ € ์ฐจ์ˆ˜ n(์ ์˜ ๊ฐœ์ˆ˜, 2~100)์„ ์„ ํƒํ•˜๊ณ , ํ•„์š”ํ•˜๋ฉด ํ‘œ์‹œํ•  ์œ ํšจ์ˆซ์ž ์ž๋ฆฟ์ˆ˜๋„ ์ง€์ •ํ•ฉ๋‹ˆ๋‹ค. ๊ณ„์‚ฐ๊ธฐ๋Š” n๊ฐœ์˜ ํ–‰์œผ๋กœ ์ด๋ฃจ์–ด์ง„ ํ‘œ๋ฅผ ์ถœ๋ ฅํ•˜๋ฉฐ, ๊ฐ ํ–‰์—๋Š” ๋…ธ๋“œ \(x_i\)์™€ ๊ทธ์— ๋Œ€์‘ํ•˜๋Š” ๊ฐ€์ค‘์น˜ \(w_i\)๊ฐ€ ๋“ค์–ด ์žˆ์Šต๋‹ˆ๋‹ค. ๋…ธ๋“œ๋Š” 0์„ ๊ธฐ์ค€์œผ๋กœ ๋Œ€์นญ์ด๋ฉฐ ๋ชจ๋‘ ๊ฐœ๊ตฌ๊ฐ„ (-1, 1) ๋‚ด๋ถ€์— ์œ„์น˜ํ•ฉ๋‹ˆ๋‹ค. ๊ฐ€์ค‘์น˜๋Š” ๋ชจ๋‘ ์–‘์ˆ˜์ด๊ณ  ๊ทธ ํ•ฉ์€ ๊ตฌ๊ฐ„ ๊ธธ์ด์™€ ๊ฐ™์€ ์ •ํ™•ํžˆ 2๊ฐ€ ๋ฉ๋‹ˆ๋‹ค. ์ž„์˜์˜ ๊ตฌ๊ฐ„ [a, b]์—์„œ ์ ๋ถ„ํ•˜๋ ค๋ฉด ๊ฐ ๋…ธ๋“œ๋ฅผ \(t_i = \frac{b-a}{2}\cdot x_i + \frac{a+b}{2}\) ๋กœ ๋ณ€ํ™˜ํ•˜๊ณ , ๊ฐ ๊ฐ€์ค‘์น˜์— \(\frac{b-a}{2}\) ๋ฅผ ๊ณฑํ•˜๋ฉด ๋ฉ๋‹ˆ๋‹ค.

๊ณต์‹ ํ’€์ด

๋…ธ๋“œ๋Š” ๋ฅด์žฅ๋“œ๋ฅด ๋‹คํ•ญ์‹ \(P_n(x)\)์˜ n๊ฐœ ๊ทผ์œผ๋กœ, ๋ณด๋„ค(Bonnet) ์ ํ™”์‹์œผ๋กœ ๊ตฌ์„ฑ๋ฉ๋‹ˆ๋‹ค: \(P_0 = 1\), \(P_1 = x\), ๊ทธ๋ฆฌ๊ณ 

$$P_k = \frac{(2k-1)\cdot x\cdot P_{k-1} - (k-1)\cdot P_{k-2}}{k}$$

์ž…๋‹ˆ๋‹ค. ๊ฐ ๊ฐ€์ค‘์น˜๋Š”

$$w_i = \frac{2}{\left(1 - x_i^{2}\right)\left[P_n^{\prime}(x_i)\right]^{2}}$$

๋กœ ์ฃผ์–ด์ง€๋ฉฐ, ์—ฌ๊ธฐ์„œ ๋„ํ•จ์ˆ˜๋Š”

$$P_n^{\prime}(x) = \frac{n\cdot\left(x\cdot P_n(x) - P_{n-1}(x)\right)}{x^{2} - 1}$$

์ž…๋‹ˆ๋‹ค. ๊ทผ์€ \(x = \cos\!\left(\frac{\pi(i - 0.25)}{n + 0.5}\right)\) ๋ฅผ ์ดˆ๊ธฐ๊ฐ’์œผ๋กœ ํ•œ ๋‰ดํ„ด ๋ฒ•์œผ๋กœ ์ฐพ์œผ๋ฉฐ, ๋ช‡ ๋ฒˆ์˜ ๋ฐ˜๋ณต๋งŒ์— ์ˆ˜๋ ดํ•ฉ๋‹ˆ๋‹ค.

[-1,1] ๊ตฌ๊ฐ„์˜ ๋งค๋„๋Ÿฌ์šด ํ•จ์ˆ˜๋ฅผ ๊ฐ€์ค‘ ๊ธฐ์—ฌ๋ฅผ ๊ฐ–๋Š” ์—ฌ๋Ÿฌ ๋ถˆ๊ท ์ผํ•œ ๋Œ€์นญ ์ ˆ์ ์—์„œ ํ‘œ๋ณธํ™”
๊ฐ€์šฐ์Šค-๋ฅด์žฅ๋“œ๋ฅด ๊ตฌ์ ๋ฒ•์€ ํŠน๋ณ„ํžˆ ์„ ํƒ๋œ ์ ˆ์ ์—์„œ์˜ ํ•จ์ˆ˜๊ฐ’์„ ๊ฐ€์ค‘ํ•ฉ์œผ๋กœ ๋”ํ•ด ์ ๋ถ„์„ ๊ทผ์‚ฌํ•ฉ๋‹ˆ๋‹ค.

๊ณ„์‚ฐ ์˜ˆ์‹œ (n = 3)

\(P_3\)์˜ ๊ทผ์€ \(x = 0\) ๊ณผ \(x = \pm\sqrt{3/5} = \pm 0.7745966692\) ์ž…๋‹ˆ๋‹ค. \(x = 0\) ์—์„œ์˜ ๊ฐ€์ค‘์น˜๋Š” \(8/9 = 0.888888889\) ์ด๊ณ , \(x = \pm 0.7745966692\) ๊ฐ๊ฐ์˜ ๊ฐ€์ค‘์น˜๋Š” \(5/9 = 0.555555556\) ์ž…๋‹ˆ๋‹ค. ๊ฐ€์ค‘์น˜ ํ•ฉ์€ $$\frac{5}{9} + \frac{8}{9} + \frac{5}{9} = 2$$ ์ด๋ฉฐ, ์ด 3์  ๊ณต์‹์€ ์ฐจ์ˆ˜ 5 ์ดํ•˜์˜ ๋‹คํ•ญ์‹์„ ์ •ํ™•ํ•˜๊ฒŒ ์ ๋ถ„ํ•ฉ๋‹ˆ๋‹ค.

[-1,1] ๊ตฌ๊ฐ„์˜ ๋Œ€์นญ์ธ ์„ธ ๊ฐ€์šฐ์Šค-๋ฅด์žฅ๋“œ๋ฅด ์ ˆ์ , ์ค‘์•™ ์ ˆ์ ์ด ๋” ํฐ ๊ฐ€์ค‘์น˜๋ฅผ ๊ฐ€์ง
3์  ๊ณต์‹์€ ํ•˜๋‚˜์˜ ์ค‘์•™ ์ ˆ์ ๊ณผ ๋Œ€์นญ์ธ ๋‘ ๊ฐœ์˜ ์™ธ์ธก ์ ˆ์ ์„ ์‚ฌ์šฉํ•˜๋ฉฐ, ์ค‘์•™ ์ ˆ์ ์— ๊ฐ€์žฅ ํฐ ๊ฐ€์ค‘์น˜๋ฅผ ๋ถ€์—ฌํ•ฉ๋‹ˆ๋‹ค.

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

๊ฐ€์ค‘์น˜์˜ ํ•ฉ์ด ์™œ 2์ธ๊ฐ€์š”? ์ƒ์ˆ˜ ํ•จ์ˆ˜ \(f(x) = 1\) ์„ [-1, 1]์—์„œ ์ ๋ถ„ํ•˜๋ฉด 2๊ฐ€ ๋˜๊ณ , ๊ตฌ์ ๋ฒ•์€ ์ƒ์ˆ˜๋ฅผ ์ •ํ™•ํžˆ ์žฌํ˜„ํ•ด์•ผ ํ•˜๋ฏ€๋กœ ๊ฐ€์ค‘์น˜์˜ ์ดํ•ฉ์€ ๊ตฌ๊ฐ„ ๊ธธ์ด์™€ ๊ฐ™์•„์•ผ ํ•ฉ๋‹ˆ๋‹ค.

๊ฐ’์€ ์–ผ๋งˆ๋‚˜ ์ •ํ™•ํ•œ๊ฐ€์š”? ์ด ๊ณ„์‚ฐ๊ธฐ๋Š” ๋ฐฐ์ •๋ฐ€๋„(double precision)๋กœ ๊ตฌํ˜„๋˜์–ด ์•ฝ 15์ž๋ฆฌ์˜ ์œ ํšจ์ˆซ์ž๋ฅผ ์ œ๊ณตํ•ฉ๋‹ˆ๋‹ค. ๊ทธ๋ณด๋‹ค ๋งŽ์€ ์ž๋ฆฟ์ˆ˜๋ฅผ ์š”์ฒญํ•˜๋ฉด ๋ฐฐ์ •๋ฐ€๋„๊ฐ€ ํ‘œํ˜„ํ•  ์ˆ˜ ์žˆ๋Š” ๋ฒ”์œ„๋กœ ๋ฐ˜์˜ฌ๋ฆผ๋ฉ๋‹ˆ๋‹ค.

์ •ํ™•ํ•˜๊ฒŒ ์ ๋ถ„๋˜๋Š” ์ตœ๋Œ€ ์ฐจ์ˆ˜๋Š” ์–ผ๋งˆ์ธ๊ฐ€์š”? n์  ๊ณต์‹์€ ์ฐจ์ˆ˜๊ฐ€ 2n-1 ์ดํ•˜์ธ ๋ชจ๋“  ๋‹คํ•ญ์‹์— ๋Œ€ํ•ด ์ •ํ™•ํ•ฉ๋‹ˆ๋‹ค.

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