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

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

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

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

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

рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрдж P_n(x) рдЯреЗрдмрд▓
P3(x)
101 points computed by Bonnet's recursion
рдШрд╛рдд n 3
рдкрдВрдХреНрддрд┐рдпрд╛рдБ 101
рдкрд╣рд▓рд╛ рдорд╛рди P_n(x_0) -1
рдЕрдВрддрд┐рдо рдорд╛рди P_n(x_last) 1
x P_n(x)
-1 -1
-0.98 -0.88298
-0.96 -0.77184
-0.94 -0.66646
-0.92 -0.56672
-0.9 -0.4725
-0.88 -0.38368
-0.86 -0.30014
-0.84 -0.22176
-0.82 -0.14842
-0.8 -0.08
-0.78 -0.01638
-0.76 0.04256
-0.74 0.09694
-0.72 0.14688
-0.7 0.1925
-0.68 0.23392
-0.66 0.27126
-0.64 0.30464
-0.62 0.33418
-0.6 0.36
-0.58 0.38222
-0.56 0.40096
-0.54 0.41634
-0.52 0.42848
-0.5 0.4375
-0.48 0.44352
-0.46 0.44666
-0.44 0.44704
-0.42 0.44478
-0.4 0.44
-0.38 0.43282
-0.36 0.42336
-0.34 0.41174
-0.32 0.39808
-0.3 0.3825
-0.28 0.36512
-0.26 0.34606
-0.24 0.32544
-0.22 0.30338
-0.2 0.28
-0.18 0.25542
-0.16 0.22976
-0.14 0.20314
-0.12 0.17568
-0.1 0.1475
-0.08 0.11872
-0.06 0.08946
-0.04 0.05984
-0.02 0.02998
0 -0
0.02 -0.02998
0.04 -0.05984
0.06 -0.08946
0.08 -0.11872
0.1 -0.1475
0.12 -0.17568
0.14 -0.20314
0.16 -0.22976
0.18 -0.25542
0.2 -0.28
0.22 -0.30338
0.24 -0.32544
0.26 -0.34606
0.28 -0.36512
0.3 -0.3825
0.32 -0.39808
0.34 -0.41174
0.36 -0.42336
0.38 -0.43282
0.4 -0.44
0.42 -0.44478
0.44 -0.44704
0.46 -0.44666
0.48 -0.44352
0.5 -0.4375
0.52 -0.42848
0.54 -0.41634
0.56 -0.40096
0.58 -0.38222
0.6 -0.36
0.62 -0.33418
0.64 -0.30464
0.66 -0.27126
0.68 -0.23392
0.7 -0.1925
0.72 -0.14688
0.74 -0.09694
0.76 -0.04256
0.78 0.01638
0.8 0.08
0.82 0.14842
0.84 0.22176
0.86 0.30014
0.88 0.38368
0.9 0.4725
0.92 0.56672
0.94 0.66646
0.96 0.77184
0.98 0.88298
1 1

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

рдпрд╣ рдЯреВрд▓ рдЖрдкрдХреЗ рдЪреБрдиреЗ рд╣реБрдП рдШрд╛рдд n рдХреЗ рд▓рд┐рдП рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрдж \(P_n(x)\) рдХреЗ рдорд╛рдиреЛрдВ рдХреА рдПрдХ рдЯреЗрдмрд▓ рддреИрдпрд╛рд░ рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕реЗ x рдХреЗ рдХреНрд░рдорд┐рдХ рдорд╛рдиреЛрдВ рдкрд░ рдирд┐рдХрд╛рд▓рд╛ рдЬрд╛рддрд╛ рд╣реИ, рдФрд░ рд╕рд╛рде рд╣реА рдЙрд╕рд╕реЗ рд╕рдВрдмрдВрдзрд┐рдд рд╡рдХреНрд░ (curve) рднреА рдмрдирд╛рддрд╛ рд╣реИред рдЖрдк рдШрд╛рдд, x рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди, рдмрдврд╝реЛрддрд░реА (step) рдФрд░ рдХрд┐рддрдиреА рдкрдВрдХреНрддрд┐рдпрд╛рдБ рдЪрд╛рд╣рд┐рдП тАФ рдпрд╣ рддрдп рдХрд░рддреЗ рд╣реИрдВ; рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рд╣рд░ рдЬреЛрдбрд╝реА \((x, P_n(x))\) рдХреЗ рд╕рд╛рде рдПрдХ рд░реЗрдЦрд╛ рдЧреНрд░рд╛рдл рд▓реМрдЯрд╛рддрд╛ рд╣реИред рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрдж [-1, 1] рдЕрдВрддрд░рд╛рд▓ рдкрд░ рд▓рд╛рдВрдмрдХреЛрдгреАрдп (orthogonal) рдмрд╣реБрдкрджреЛрдВ рдХрд╛ рдПрдХ рдкреНрд░рд╕рд┐рджреНрдз рдкрд░рд┐рд╡рд╛рд░ рд╣реИрдВ рдФрд░ рднреМрддрд┐рдХреА рддрдерд╛ рдЕрдиреБрдкреНрд░рдпреБрдХреНрдд рдЧрдгрд┐рдд рдореЗрдВ рд╣рд░ рдЬрдЧрд╣ рджрд┐рдЦрддреЗ рд╣реИрдВ тАФ рд▓рд╛рдкреНрд▓рд╛рд╕ рд╕рдореАрдХрд░рдг рдХреЗ рд╣рд▓, рдмрд╣реБрдзреНрд░реБрд╡ рд╡рд┐рд╕реНрддрд╛рд░ (multipole expansion), рдЧреЛрд▓реАрдп рд╣рд╛рд░реНрдореЛрдирд┐рдХреНрд╕ рдФрд░ рдЧреЙрд╕рд┐рдпрди рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдореЗрдВред

x рдХреЗ рдорд╛рдЗрдирд╕ рдПрдХ рд╕реЗ рдПрдХ рддрдХ рдкрд╣рд▓реЗ рдХреБрдЫ рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрджреЛрдВ рдХреЗ рд╡рдХреНрд░
рдЕрдВрддрд░рд╛рд▓ [-1, 1] рдкрд░ рдкрд╣рд▓реЗ рдХреБрдЫ рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрдж \(P_n(x)\)ред

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

n (рдШрд╛рдд) рдХреЛ рдПрдХ рдЕрдЛрдгрд╛рддреНрдордХ рдкреВрд░реНрдгрд╛рдВрдХ рдХреЗ рд░реВрдк рдореЗрдВ рдбрд╛рд▓реЗрдВ (0, 1, 2, тАж)ред x рдХрд╛ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди (рдЕрдХреНрд╕рд░ -1), рдХреНрд░рдорд┐рдХ x рдорд╛рдиреЛрдВ рдХреЗ рдмреАрдЪ рдмрдврд╝реЛрддрд░реА (step) (рдЬреИрд╕реЗ 0.02), рдФрд░ рдЬрд┐рддрдиреА рдкрдВрдХреНрддрд┐рдпрд╛рдБ (rows) рдЖрдк рдЪрд╛рд╣рддреЗ рд╣реИрдВ тАФ рд╕рднреА рд╕реЗрдЯ рдХрд░реЗрдВред i-рд╡реАрдВ рдкрдВрдХреНрддрд┐ рдореЗрдВ \(x = \text{startX} + i \times \text{step}\) рдХрд╛ рдкреНрд░рдпреЛрдЧ рд╣реЛрддрд╛ рд╣реИред рд╣рд╛рд▓рд╛рдБрдХрд┐ рдпреЗ рдмрд╣реБрдкрдж [-1, 1] рдкрд░ рд╕рдмрд╕реЗ рдЕрдзрд┐рдХ рд╕рд╛рд░реНрдердХ рд╣реИрдВ, рдпрд╣ рд╕реВрддреНрд░ рдХрд┐рд╕реА рднреА рд╡рд╛рд╕реНрддрд╡рд┐рдХ x рдХреЗ рд▓рд┐рдП рдХрд╛рдо рдХрд░рддрд╛ рд╣реИ тАФ рдзреНрдпрд╛рди рджреЗрдВ рдХрд┐ рдЗрд╕ рдЕрдВрддрд░рд╛рд▓ рдХреЗ рдмрд╛рд╣рд░ рдорд╛рди рдмрд╣реБрдд рддреЗрдЬрд╝реА рд╕реЗ рдмрдврд╝рддреЗ рд╣реИрдВред

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

рдмрдВрдж-рд░реВрдк (closed form) рд╕реВрддреНрд░реЛрдВ рдХреЛ рд╡рд┐рд╕реНрддреГрдд рдХрд░рдиреЗ рдХреЗ рдмрдЬрд╛рдп, рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕реНрдерд┐рд░рддрд╛ рдХреЗ рд▓рд┐рдП рдмреЙрдиреЗрдЯ рдХреЗ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕реВрддреНрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИ: \(P_0(x) = 1\) рдФрд░ \(P_1(x) = x\) рд╕реЗ рд╢реБрд░реБрдЖрдд рдХрд░реЗрдВ, рдлрд┐рд░ $$P_{k+1}(x) = \frac{(2k+1)\cdot x\cdot P_k(x) - k\cdot P_{k-1}(x)}{k+1}$$ рдХреЛ рджреЛрд╣рд░рд╛рдПрдВред рдкрд╣рд▓реЗ рдХреБрдЫ рдмрдВрдж-рд░реВрдк рд╕реВрддреНрд░ рд╣реИрдВ \(P_2 = \frac{3x^2 - 1}{2}\), \(P_3 = \frac{5x^3 - 3x}{2}\), рдФрд░ \(P_4 = \frac{35x^4 - 30x^2 + 3}{8}\)ред

рдмреЛрдирдЯ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдХреЛ рджрд░реНрд╢рд╛рддрд╛ рдЖрд░реЗрдЦ рдЬреЛ рдкрд┐рдЫрд▓реЗ рджреЛ рдмрд╣реБрдкрджреЛрдВ рдХреЛ рдЕрдЧрд▓реЗ рдореЗрдВ рдЬреЛрдбрд╝рддрд╛ рд╣реИ
рдмреЛрдирдЯ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдкреНрд░рддреНрдпреЗрдХ рдмрд╣реБрдкрдж рдХреЛ рдкрд┐рдЫрд▓реЗ рджреЛ рд╕реЗ рдмрдирд╛рддреА рд╣реИред

рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг

n = 3 рдФрд░ x = 0.5 рдХреЗ рд▓рд┐рдП: \(P_0 = 1\), \(P_1 = 0.5\)ред рдлрд┐рд░ $$P_2 = \frac{3\cdot 0.5\cdot 0.5 - 1}{2} = -0.125,$$ рдФрд░ $$P_3 = \frac{5\cdot 0.5\cdot (-0.125) - 2\cdot 0.5}{3} = \frac{-1.3125}{3} = -0.4375.$$ рдмрдВрдж-рд░реВрдк рд╕реВрддреНрд░ \(\frac{5x^3 - 3x}{2}\) рднреА рд╡рд╣реА рдкрд░рд┐рдгрд╛рдо рджреЗрддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕реВрддреНрд░ рдХреА рдкреБрд╖реНрдЯрд┐ рд╣реЛ рдЬрд╛рддреА рд╣реИред

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

n = 0 рдХрд╛ рдкрд░рд┐рдгрд╛рдо рдХреНрдпрд╛ рд╣реЛрддрд╛ рд╣реИ? рд╣рд░ x рдХреЗ рд▓рд┐рдП рдПрдХ рд╕реНрдерд┐рд░ рдорд╛рди 1, рдЗрд╕рд▓рд┐рдП рдЧреНрд░рд╛рдл рдПрдХ рд╕рдкрд╛рдЯ рдХреНрд╖реИрддрд┐рдЬ рд░реЗрдЦрд╛ рд╣реЛрддрд╛ рд╣реИред рдЕрдВрддрд┐рдо-рдмрд┐рдВрджреБ (endpoint) рдорд╛рди рдХреНрдпрд╛ рд╣реЛрддреЗ рд╣реИрдВ? рд╣рд░ рд▓реАрдЬреЗрдВрдбреНрд░ рдмрд╣реБрдкрдж \(P_n(1) = 1\) рдФрд░ \(P_n(-1) = (-1)^n\) рдХреЛ рд╕рдВрддреБрд╖реНрдЯ рдХрд░рддрд╛ рд╣реИред рд╕реНрдкрд╖реНрдЯ рд╕реВрддреНрд░реЛрдВ рдХреЗ рдмрдЬрд╛рдп рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреНрдпреЛрдВ? рддреАрди-рдкрдж рд╡рд╛рд▓рд╛ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕реВрддреНрд░ рдХрд┐рд╕реА рднреА рдШрд╛рдд рдХреЗ рд▓рд┐рдП рддреЗрдЬрд╝ рдФрд░ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╕реНрдерд┐рд░ рд╣реЛрддрд╛ рд╣реИ, рдФрд░ рдЙрдЪреНрдЪ-рдШрд╛рдд рд╕реНрдкрд╖реНрдЯ рдмрд╣реБрдкрджреЛрдВ рдореЗрдВ рд╣реЛрдиреЗ рд╡рд╛рд▓реА рд░рджреНрджреАрдХрд░рдг (cancellation) рддреНрд░реБрдЯрд┐рдпреЛрдВ рд╕реЗ рдмрдЪрд╛рддрд╛ рд╣реИред

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