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

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

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

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

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

Gauss-Lobatto rule, n = 20
20
points on the interval [-1, 1] ┬╖ sum of weights = 2
i рдиреЛрдб x_i рд╡реЗрдЯ w_i
1 -1 0.005263157894737
2 -0.980743704893914 0.032237123188489
3 -0.935934498812665 0.057181802127567
4 -0.86687797808995 0.08063176399612
5 -0.775368260952056 0.101991499699451
6 -0.663776402290311 0.120709227628675
7 -0.534992864031886 0.136300482358724
8 -0.392353183713909 0.148361554070917
9 -0.239551705922986 0.156580102647475
10 -0.080545937238822 0.160743286387846
11 0.080545937238822 0.160743286387846
12 0.239551705922986 0.156580102647475
13 0.392353183713909 0.148361554070917
14 0.534992864031886 0.136300482358724
15 0.663776402290311 0.120709227628675
16 0.775368260952056 0.101991499699451
17 0.86687797808995 0.08063176399612
18 0.935934498812665 0.057181802127567
19 0.980743704893914 0.032237123188489
20 1 0.005263157894737

рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ

рдпрд╣ рдЯреВрд▓ рд╕рдВрджрд░реНрдн рдЕрдВрддрд░рд╛рд▓ [-1, 1] рдкрд░ рднрд╛рд░ рдлрд▓рди \(w(x) = 1\) рдХреЗ рд╕рд╛рде \(n\)-рдмрд┐рдВрджреБ рдЧреЙрд╕-рд▓реЛрдмрд╛рдЯреЛ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдирд┐рдпрдо рдХреЗ рдиреЛрдбреНрд╕ (рдПрдмреНрд╕рд┐рд╕рд╛) \(x_i\) рдФрд░ рд╡реЗрдЯреНрд╕ \(w_i\) рдХреА рдЧрдгрдирд╛ рдХрд░рддрд╛ рд╣реИред рдорд╛рдирдХ рдЧреЙрд╕-рд▓реЗрдЬреЗрдВрдбреНрд░реЗ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХреЗ рд╡рд┐рдкрд░реАрдд, рдЧреЙрд╕-рд▓реЛрдмрд╛рдЯреЛ рдирд┐рдпрдо рд╣рдореЗрд╢рд╛ рджреЛрдиреЛрдВ рд╕рд┐рд░реЛрдВ \(x = -1\) рдФрд░ \(x = +1\) рдХреЛ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдиреЛрдбреНрд╕ рдХреЗ рд░реВрдк рдореЗрдВ рд░рдЦрддрд╛ рд╣реИред рдпрд╣ рддрдм рдмреЗрд╣рдж рдЙрдкрдпреЛрдЧреА рд╣реЛрддрд╛ рд╣реИ рдЬрдм рд╕реАрдорд╛ (boundary) рдХреЗ рдорд╛рди рдорд╛рдпрдиреЗ рд░рдЦрддреЗ рд╣реЛрдВ тАФ рдЙрджрд╛рд╣рд░рдг рдХреЗ рд▓рд┐рдП рд╕реНрдкреЗрдХреНрдЯреНрд░рд▓ рдПрд▓рд┐рдореЗрдВрдЯ рд╡рд┐рдзрд┐рдпреЛрдВ рдореЗрдВред рдпрд╣ рд╡рд┐рд╢реБрджреНрдз рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╡рд┐рд╢реНрд▓реЗрд╖рдг рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рд╕рдорд╛рди рд░реВрдк рд╕реЗ рд▓рд╛рдЧреВ рд╣реЛрддрд╛ рд╣реИ; рдпрд╣ рдХрд┐рд╕реА рджреЗрд╢ рдпрд╛ рдХреНрд╖реЗрддреНрд░ рд╡рд┐рд╢реЗрд╖ рдкрд░ рдирд┐рд░реНрднрд░ рдирд╣реАрдВ рдХрд░рддрд╛ред

Number line from -1 to 1 with quadrature nodes including both endpoints, each marked by a vertical weight bar
Gauss-Lobatto nodes on [-1, 1] include both endpoints; bar heights suggest the associated weights.

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

рдмрд┐рдВрджреБрдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ \(n\) рдЪреБрдиреЗрдВ (2 рд╕реЗ 100 рдХреЗ рдмреАрдЪ) рдФрд░ рдЪрд╛рд╣реЗрдВ рддреЛ рдкреНрд░рджрд░реНрд╢рди рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рднреА рдирд┐рд░реНрдзрд╛рд░рд┐рдд рдХрд░реЗрдВред рдХреИрд▓рдХреБрд▓реЗрдЯрд░ \(n\) рдкрдВрдХреНрддрд┐рдпреЛрдВ рдХреА рдПрдХ рддрд╛рд▓рд┐рдХрд╛ рд▓реМрдЯрд╛рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рд╣рд░ рдкрдВрдХреНрддрд┐ рдореЗрдВ рдиреЛрдб \(x_i\) рдФрд░ рдЙрд╕рдХрд╛ рд╡реЗрдЯ \(w_i\) рджрд┐рдпрд╛ рд╣реЛрддрд╛ рд╣реИред рдиреЛрдбреНрд╕ 0 рдХреЗ рд╕рд╛рдкреЗрдХреНрд╖ рд╕рдордорд┐рдд рд╣реЛрддреЗ рд╣реИрдВ рдФрд░ рд╡реЗрдЯреНрд╕ рднреА рд╕рдордорд┐рдд рд╣реЛрддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП \(x_i\) рдФрд░ \(-x_i\) рдХрд╛ рд╡реЗрдЯ рдПрдХ рд╣реА рд░рд╣рддрд╛ рд╣реИред рдПрдХ рдЕрдВрддрд░реНрдирд┐рд╣рд┐рдд рдЬрд╛рдБрдЪ рдХреЗ рд░реВрдк рдореЗрдВ, рд╕рднреА рд╡реЗрдЯреНрд╕ рдХрд╛ рдпреЛрдЧ 2 рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рдЕрдВрддрд░рд╛рд▓ рдХреА рд▓рдВрдмрд╛рдИ рдХреЗ рдмрд░рд╛рдмрд░ рд╣реИред

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

рдпрд╣ рдирд┐рдпрдо рд╕рдорд╛рдХрд▓рди (integral) рдХреЛ рдпреЛрдЧ \(w_1 f(x_1) + \ldots + w_n f(x_n)\) рдХреЗ рд░реВрдк рдореЗрдВ рд╕рдиреНрдирд┐рдХрдЯ рдХрд░рддрд╛ рд╣реИ рдФрд░ рдбрд┐рдЧреНрд░реА \(2n-3\) рддрдХ рдХреЗ рдмрд╣реБрдкрджреЛрдВ рдХреЗ рд▓рд┐рдП рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИред

$$\int_{-1}^{1} f(x)\,dx \approx \sum_{i=1}^{n} w_i\, f(x_i)$$

рдЖрдВрддрд░рд┐рдХ рдиреЛрдбреНрд╕ \(x_2, \ldots, x_{n-1}\) рджрд░рдЕрд╕рд▓ \(P_{n-1}'(x)\) рдХреЗ \(n-2\) рд╢реВрдиреНрдп рд╣реИрдВ, рдпрд╛рдиреА рдбрд┐рдЧреНрд░реА \(n-1\) рдХреЗ рд▓реЗрдЬреЗрдВрдбреНрд░реЗ рдмрд╣реБрдкрдж рдХреЗ рдЕрд╡рдХрд▓рдЬ (derivative) рдХреЗ рдореВрд▓ред рд╕рд┐рд░реЛрдВ рдХрд╛ рд╡реЗрдЯ \(\frac{2}{n(n-1)}\) рд╣реЛрддрд╛ рд╣реИ, рдЬрдмрдХрд┐ рдкреНрд░рддреНрдпреЗрдХ рдЖрдВрддрд░рд┐рдХ рдиреЛрдб \(x_i\) рдХрд╛ рд╡реЗрдЯ \(\frac{2}{n(n-1)\left[P_{n-1}(x_i)\right]^{2}}\) рд╣реЛрддрд╛ рд╣реИред

$$\text{where}\quad \left\{ \begin{aligned} x_1 &= -1,\quad x_{n} = 1 \\ x_i &: P_{n-1}^{\prime}(x_i) = 0 \quad(\text{interior}) \\ w_{1} &= w_{n} = \frac{2}{n\,(n-1)} \\ w_i &= \frac{2}{n\,(n-1)\,\left[P_{n-1}(x_i)\right]^{2}} \end{aligned} \right.$$

рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдЖрдВрддрд░рд┐рдХ рдореВрд▓реЛрдВ рдХреЛ рдиреНрдпреВрдЯрди рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ (Newton iteration) рджреНрд╡рд╛рд░рд╛ рдЦреЛрдЬрддрд╛ рд╣реИ, рдЬрд┐рд╕рдХреА рд╢реБрд░реБрдЖрдд рдЪреЗрдмрд┐рд╢реЗрд╡-рдЧреЙрд╕-рд▓реЛрдмрд╛рдЯреЛ рдЕрдиреБрдорд╛рдиреЛрдВ \(\cos\!\left(\frac{\pi j}{n-1}\right)\) рд╕реЗ рд╣реЛрддреА рд╣реИ, рдФрд░ рдЗрд╕рд╕реЗ рдкреВрд░реА рдбрдмрд▓ рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ (рд▓рдЧрднрдЧ 15-16 рд╕рд╛рд░реНрдердХ рдЕрдВрдХ) рдорд┐рд▓рддреА рд╣реИред

Curve f(x) over [-1,1] approximated by weighted samples at Gauss-Lobatto nodes including the endpoints
The integral is approximated by a weighted sum of function values at the nodes, with the two endpoints always included.

рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг (n = 4)

рдЖрдВрддрд░рд┐рдХ рдиреЛрдбреНрд╕ рд╕рдореАрдХрд░рдг \(P_3'(x) = \frac{15x^2 - 3}{2} = 0\) рдХреЛ рд╣рд▓ рдХрд░рддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП \(x = \pm\frac{1}{\sqrt{5}} = \pm 0.4472135955\)ред рд╕рд┐рд░реЗ рдХрд╛ рд╡реЗрдЯ \(\frac{2}{4 \cdot 3} = \frac{1}{6} = 0.1666666667\) рд╣реИред рдЖрдВрддрд░рд┐рдХ рдиреЛрдбреНрд╕ рдХреЗ рд▓рд┐рдП \(P_3\!\left(\frac{1}{\sqrt{5}}\right) = -0.4472135955\) рд╣реИ, рдЬрд┐рд╕рдХрд╛ рд╡рд░реНрдЧ \(0.2\) рд╣реЛрддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рд╡реЗрдЯ \(\frac{2}{4 \cdot 3 \cdot 0.2} = \frac{5}{6} = 0.8333333333\) рдорд┐рд▓рддрд╛ рд╣реИред рдпреЛрдЧ \(\frac{1}{6} + \frac{5}{6} + \frac{5}{6} + \frac{1}{6} = 2\) рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рдирд┐рдпрдо рдХреА рдкреБрд╖реНрдЯрд┐ рдХрд░рддрд╛ рд╣реИред

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

рдпрд╣ рдЧреЙрд╕-рд▓реЗрдЬреЗрдВрдбреНрд░реЗ рд╕реЗ рдХреИрд╕реЗ рдЕрд▓рдЧ рд╣реИ? рдЧреЙрд╕-рд▓реЗрдЬреЗрдВрдбреНрд░реЗ рд╕рднреА рдиреЛрдбреНрд╕ рдХреЛ рдкреВрд░реА рддрд░рд╣ \((-1, 1)\) рдХреЗ рдЕрдВрджрд░ рд░рдЦрддрд╛ рд╣реИ рдФрд░ рдбрд┐рдЧреНрд░реА \(2n-1\) рддрдХ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИред рдЧреЙрд╕-рд▓реЛрдмрд╛рдЯреЛ рджреЛрдиреЛрдВ рд╕рд┐рд░реЛрдВ рдХреЛ рдиреЛрдбреНрд╕ рдХреЗ рд░реВрдк рдореЗрдВ рддрдп рдХрд░рддрд╛ рд╣реИ рдФрд░ рдбрд┐рдЧреНрд░реА \(2n-3\) рддрдХ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИ тАФ рдпрд╛рдиреА рд╕реАрдорд╛ рдХреЛ рд╢рд╛рдорд┐рд▓ рдХрд░рдиреЗ рдХреЗ рдмрджрд▓реЗ рдкрд░рд┐рд╢реБрджреНрдзрддрд╛ рдХреА рджреЛ рдбрд┐рдЧреНрд░реА рдХрд╛ рддреНрдпрд╛рдЧред

рдЗрдиреНрд╣реЗрдВ рдХрд┐рд╕реА рд╕рд╛рдорд╛рдиреНрдп рдЕрдВрддрд░рд╛рд▓ [a, b] рдкрд░ рдХреИрд╕реЗ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд░реВрдБ? рд╣рд░ рдиреЛрдб рдХреЛ \(x \to \frac{b-a}{2}\, x + \frac{a+b}{2}\) рд╕реЗ рд░реВрдкрд╛рдВрддрд░рд┐рдд рдХрд░реЗрдВ рдФрд░ рд╣рд░ рд╡реЗрдЯ рдХреЛ \(\frac{b-a}{2}\) рд╕реЗ рдЧреБрдгрд╛ рдХрд░реЗрдВред рдпрд╣ рдкреЗрдЬ рдХреЗрд╡рд▓ \([-1, 1]\) рдХреЗ рдорд╛рди рджреЗрддрд╛ рд╣реИред

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

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