MCP рдХреЗ рдорд╛рдзреНрдпрдо рд╕реЗ рдХрдиреЗрдХреНрдЯ рдХрд░реЗрдВ тЖТ

рдЧрдгрдирд╛ рджрд░реНрдЬ рдХрд░реЗрдВ

рд╕реВрддреНрд░ (рдлреЙрд░реНрдореВрд▓рд╛)

рд╡рд┐рдЬреНрдЮрд╛рдкрди

рдкрд░рд┐рдгрд╛рдо

20-Point Gauss-Legendre Rule on [-1, 1]
2
рд╡реЗрдЯреНрд╕ рдХрд╛ рдпреЛрдЧ (2 рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП)
i рдиреЛрдб x_i рд╡реЗрдЯ w_i
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-рдмрд┐рдВрджреБ рдЧреЙрд╕-рд▓реАрдЬреЗрдВрдбреНрд░ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдирд┐рдпрдо рдХреЗ рдиреЛрдбреНрд╕ (рдПрдмреНрд╕рд┐рд╕рд╛) рдФрд░ рд╡реЗрдЯреНрд╕ рдХреА рдЧрдгрдирд╛ рдХрд░рддрд╛ рд╣реИред рдЧреЙрд╕-рд▓реАрдЬреЗрдВрдбреНрд░ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдПрдХ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕рдорд╛рдХрд▓рди (numerical integration) рд╡рд┐рдзрд┐ рд╣реИ рдЬреЛ рдХрд┐рд╕реА рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓ рдХреЛ рдлрд▓рди рдорд╛рдиреЛрдВ рдХреЗ рднрд╛рд░рд┐рдд рдпреЛрдЧ рдХреЗ рд░реВрдк рдореЗрдВ рд╕рдиреНрдирд┐рдХрдЯ рдХрд░рддреА рд╣реИ: -1 рд╕реЗ 1 рддрдХ f(x) рдХрд╛ рд╕рдорд╛рдХрд▓ рд▓рдЧрднрдЧ i рдкрд░ w_i рдЧреБрдгрд╛ f(x_i) рдХреЗ рдпреЛрдЧ рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛрддрд╛ рд╣реИред рдХреЗрд╡рд▓ n рдмрд┐рдВрджреБрдУрдВ рдХреЗ рд╕рд╛рде рдпрд╣ 2n-1 рдШрд╛рдд рддрдХ рдХреЗ рдХрд┐рд╕реА рднреА рдмрд╣реБрдкрдж рдХрд╛ рд╕рдЯреАрдХ рд╕рдорд╛рдХрд▓рди рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдпрд╣ рдЯреНрд░реИрдкреАрдЬрд╝реЙрдЗрдбрд▓ рдпрд╛ рд╕рд┐рдореНрдкрд╕рди рдЬреИрд╕реА рд╕рдорд╛рди-рджреВрд░реА рд╡рд╛рд▓реА рд╡рд┐рдзрд┐рдпреЛрдВ рдХреА рддреБрд▓рдирд╛ рдореЗрдВ рдХрд╣реАрдВ рдЕрдзрд┐рдХ рд╕рдЯреАрдХ рд╣реЛ рдЬрд╛рддрд╛ рд╣реИред

рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ

рдСрд░реНрдбрд░ n (рдмрд┐рдВрджреБрдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛, 2 рд╕реЗ 100) рдЪреБрдиреЗрдВ рдФрд░ рдЪрд╛рд╣реЗрдВ рддреЛ рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд┐рдП рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХреЛрдВ (significant digits) рдХреА рд╕рдВрдЦреНрдпрд╛ рднреА рдЪреБрдиреЗрдВред рдХреИрд▓рдХреБрд▓реЗрдЯрд░ n рдкрдВрдХреНрддрд┐рдпреЛрдВ рдХреА рдПрдХ рддрд╛рд▓рд┐рдХрд╛ рд▓реМрдЯрд╛рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рдкреНрд░рддреНрдпреЗрдХ рдкрдВрдХреНрддрд┐ рдореЗрдВ рдПрдХ рдиреЛрдб \(x_i\) рдФрд░ рдЙрд╕рдХрд╛ рд╡реЗрдЯ \(w_i\) рд╣реЛрддрд╛ рд╣реИред рдиреЛрдбреНрд╕ 0 рдХреЗ рдкрд░рд┐рддрдГ рд╕рдордорд┐рдд рд╣реЛрддреЗ рд╣реИрдВ рдФрд░ рд╕рднреА рдкреВрд░реА рддрд░рд╣ (-1, 1) рдХреЗ рднреАрддрд░ рд╕реНрдерд┐рдд рд╣реЛрддреЗ рд╣реИрдВ; рд╡реЗрдЯреНрд╕ рд╕рднреА рдзрдирд╛рддреНрдордХ рд╣реЛрддреЗ рд╣реИрдВ рдФрд░ рдЗрдирдХрд╛ рдпреЛрдЧ рдареАрдХ 2 рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рдЕрдВрддрд░рд╛рд▓ рдХреА рд▓рдВрдмрд╛рдИ рд╣реИред рдХрд┐рд╕реА рдордирдорд╛рдиреЗ рдЕрдВрддрд░рд╛рд▓ [a, b] рдкрд░ рд╕рдорд╛рдХрд▓рди рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП, рдкреНрд░рддреНрдпреЗрдХ рдиреЛрдб рдХреЛ \(t_i = \frac{b-a}{2} x_i + \frac{a+b}{2}\) рд╕реЗ рдореИрдк рдХрд░реЗрдВ рдФрд░ рдкреНрд░рддреНрдпреЗрдХ рд╡реЗрдЯ рдХреЛ \(\frac{b-a}{2}\) рд╕реЗ рдЧреБрдгрд╛ рдХрд░реЗрдВред

рд╕реВрддреНрд░ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛

рдиреЛрдбреНрд╕, рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрдж \(P_n(x)\) рдХреЗ n рдореВрд▓ (roots) рд╣реЛрддреЗ рд╣реИрдВ, рдЬрд┐рдиреНрд╣реЗрдВ рдмреЛрдиреЗрдЯ рдХреА рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ (Bonnet's recurrence) рд╕реЗ рдмрдирд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ: \(P_0=1\), \(P_1=x\), рдФрд░

$$P_k = \frac{(2k-1)\, x\, P_{k-1} - (k-1)\, 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\left(x\, 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\) рдкрд░ рд╡реЗрдЯ \(\frac{8}{9} = 0.888888889\) рд╣реИ, рдФрд░ рдкреНрд░рддреНрдпреЗрдХ \(x = \pm 0.7745966692\) рдкрд░ рд╡реЗрдЯ \(\frac{5}{9} = 0.555555556\) рд╣реИред рд╡реЗрдЯреНрд╕ рдХрд╛ рдпреЛрдЧ

$$\frac{5}{9} + \frac{8}{9} + \frac{5}{9} = 2$$

рд╣реЛрддрд╛ рд╣реИ, рдФрд░ рдпрд╣ 3-рдмрд┐рдВрджреБ рдирд┐рдпрдо 5 рдШрд╛рдд рддрдХ рдХреЗ рдмрд╣реБрдкрджреЛрдВ рдХрд╛ рд╕рдЯреАрдХ рд╕рдорд╛рдХрд▓рди рдХрд░рддрд╛ рд╣реИред

[-1,1] рдкрд░ рддреАрди рд╕рдордорд┐рдд рдЧрд╛рдЙрд╕-рд▓реЗрдЬреЗрдВрдбреНрд░ рдиреЛрдбреНрд╕, рдЬрд┐рдирдореЗрдВ рдХреЗрдВрджреНрд░реАрдп рдиреЛрдб рдХрд╛ рднрд╛рд░ рдЕрдзрд┐рдХ рд╣реИ
3-рдмрд┐рдВрджреБ рдирд┐рдпрдо рдПрдХ рдХреЗрдВрджреНрд░реАрдп рдиреЛрдб рдФрд░ рджреЛ рд╕рдордорд┐рдд рдмрд╛рд╣рд░реА рдиреЛрдбреНрд╕ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рдХреЗрдВрджреНрд░реАрдп рдиреЛрдб рдХреЛ рд╕рдмрд╕реЗ рдмрдбрд╝рд╛ рднрд╛рд░ рджрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред

рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди

рд╡реЗрдЯреНрд╕ рдХрд╛ рдпреЛрдЧ 2 рд╣реА рдХреНрдпреЛрдВ рд╣реЛрддрд╛ рд╣реИ? [-1, 1] рдкрд░ рд╕реНрдерд┐рд░ рдлрд▓рди \(f(x) = 1\) рдХрд╛ рд╕рдорд╛рдХрд▓рди 2 рджреЗрддрд╛ рд╣реИ, рдФрд░ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХреЛ рд╕реНрдерд┐рд░рд╛рдВрдХреЛрдВ рдХрд╛ рд╕рдЯреАрдХ рдкреБрдирд░реБрддреНрдкрд╛рджрди рдХрд░рдирд╛ рд╣реЛрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рд╡реЗрдЯреНрд╕ рдХрд╛ рдХреБрд▓ рдпреЛрдЧ рдЕрдВрддрд░рд╛рд▓ рдХреА рд▓рдВрдмрд╛рдИ рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдПред

рдорд╛рди рдХрд┐рддрдиреЗ рд╕рдЯреАрдХ рд╣реЛрддреЗ рд╣реИрдВ? рдпрд╣ рдПрдХ рдбрдмрд▓-рдкреНрд░рд┐рд╕рд┐рдЬрди рдкреБрдирд░реНрдирд┐рд░реНрдорд╛рдг рд╣реИ рдЬреЛ рд▓рдЧрднрдЧ 15 рд╕рд╛рд░реНрдердХ рдЕрдВрдХ рджреЗрддрд╛ рд╣реИред рдЕрдзрд┐рдХ рдЕрдВрдХ рдорд╛рдБрдЧрдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рджрд░реНрд╢рди рд╡рд┐рдХрд▓реНрдк рдЙрддрдиреЗ рддрдХ рдкреВрд░реНрдгрд╛рдВрдХрд┐рдд рдХрд┐рдП рдЬрд╛рддреЗ рд╣реИрдВ рдЬрд┐рддрдирд╛ рдбрдмрд▓ рдкреНрд░рд┐рд╕рд┐рдЬрди рдирд┐рд░реВрдкрд┐рдд рдХрд░ рд╕рдХрддрд╛ рд╣реИред

рд╕рдЯреАрдХ рд░реВрдк рд╕реЗ рд╕рдорд╛рдХрд▓рд┐рдд рдЕрдзрд┐рдХрддрдо рдШрд╛рдд рдХрд┐рддрдиреА рд╣реИ? рдПрдХ n-рдмрд┐рдВрджреБ рдирд┐рдпрдо \(2n-1\) рдШрд╛рдд рддрдХ рдХреЗ рд╕рднреА рдмрд╣реБрдкрджреЛрдВ рдХреЗ рд▓рд┐рдП рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИред

рдЕрдВрддрд┐рдо рдЕрдкрдбреЗрдЯ: