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

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

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

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

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

рдЕрдиреБрдорд╛рдирд┐рдд рд╕рдорд╛рдХрд▓рди рдорд╛рди
0.9726059074
Gauss-Laguerre quadrature, n = 10, alpha = 0
рд╡рд┐рдзрд┐ рд╕рд╛рдорд╛рдиреНрдпреАрдХреГрдд рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░реЗ (рдЧреЛрд▓рдм-рд╡реЗрд▓реНрд╢)
рдиреЛрдбреНрд╕ рдХреА рд╕рдВрдЦреНрдпрд╛ n 10
рдкреИрд░рд╛рдореАрдЯрд░ alpha 0

рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░реЗ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХреНрдпрд╛ рд╣реИ?

рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░реЗ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдПрдХ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╡рд┐рдзрд┐ рд╣реИ рдЬрд┐рд╕рд╕реЗ (0, рдЕрдирдВрдд) рдЬреИрд╕реЗ рдЕрд░реНрдз-рдЕрдирдВрдд рдЕрдВрддрд░рд╛рд▓ рдкрд░ рдРрд╕реЗ рдЕрдиреБрдЪрд┐рдд рд╕рдорд╛рдХрд▓рдиреЛрдВ (improper integrals) рдХрд╛ рдЕрдиреБрдорд╛рди рд▓рдЧрд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ рдЬрд┐рдирдХрд╛ рд╕рдорд╛рдХрд▓реНрдп (integrand) рдЪрд░рдШрд╛рддрд╛рдВрдХреА рд░реВрдк рд╕реЗ (exponential рдХреА рддрд░рд╣) рдШрдЯрддрд╛ рд╣реИред рдпрд╣ рд╡рд┐рдзрд┐ рд╕рдорд╛рдХрд▓рди рдХреЛ рдХреБрдЫ рдЦрд╛рд╕ рддрд░рд╣ рд╕реЗ рдЪреБрдиреЗ рдЧрдП рдмрд┐рдВрджреБрдУрдВ тАФ рдЬрд┐рдиреНрд╣реЗрдВ рдиреЛрдб (nodes) рдХрд╣рддреЗ рд╣реИрдВ тАФ рдкрд░ рдореВрд▓реНрдпрд╛рдВрдХрд┐рдд рдХрд┐рдП рдЧрдП рднрд╛рд░рд┐рдд рдпреЛрдЧ (weighted sum) рд╕реЗ рдмрджрд▓ рджреЗрддреА рд╣реИред рдЪреБрдиреА рд╣реБрдИ рдХреЛрдЯрд┐ n рдХреЗ рд▓рд┐рдП рдпрд╣ рдирд┐рдпрдо \(2n-1\) рддрдХ рдХреА рдШрд╛рдд рд╡рд╛рд▓реЗ рдХрд┐рд╕реА рднреА рдмрд╣реБрдкрдж рдХреЗ рд▓рд┐рдП (рднрд╛рд░ \(x^{\alpha} e^{-x}\) рдХреЗ рд╕рд╛рдкреЗрдХреНрд╖) рдкреВрд░реНрдгрддрдГ рд╕рдЯреАрдХ рд╣реЛрддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдпрд╣ рд╕рд╣рдЬ (smooth) рд╕рдорд╛рдХрд▓реНрдпреЛрдВ рдХреЗ рд▓рд┐рдП рдорд╛рддреНрд░ рдХреБрдЫ рд╣реА рдореВрд▓реНрдпрд╛рдВрдХрдиреЛрдВ рдореЗрдВ рдЙрд▓реНрд▓реЗрдЦрдиреАрдп рд░реВрдк рд╕реЗ рд╕рдЯреАрдХ рдкрд░рд┐рдгрд╛рдо рджреЗрддрд╛ рд╣реИред

рднрд╛рд░ рдлрд▓рди x^alpha рдЧреБрдгрд╛ e^рдЛрдг x рдЬреЛ рдзрдирд╛рддреНрдордХ x-рдЕрдХреНрд╖ рдкрд░ рдШрдЯрддрд╛ рд╡рдХреНрд░ рд╣реИ, рдЪрд┐рд╣реНрдирд┐рдд рдирдореВрдирд╛ рдиреЛрдбреНрд╕ рдХреЗ рд╕рд╛рде
рдЧреЙрд╕-рд▓рд╛рдЧреНрд╡реЗрд░ рд╕рдорд╛рдХрд▓рди (0, рдЕрдирдВрдд) рдкрд░ рднрд╛рд░ \(x^{\alpha} e^{-x}\) рдХреЗ рд╕рд╛рде рдХрд░рддрд╛ рд╣реИ, рдЬреЛ рд╢реВрдиреНрдп рдХреА рдУрд░ рдШрдЯрддрд╛ рд╣реИред

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

рд╕рдмрд╕реЗ рдкрд╣рд▓реЗ рдЗрдирдкреБрдЯ рдореЛрдб рдЪреБрдиреЗрдВред рдпрджрд┐ рдЖрдкрдХрд╛ рд╕рдорд╛рдХрд▓рди рдкрд╣рд▓реЗ рд╕реЗ рд╣реА \(\int x^{\alpha} e^{-x} f(x)\, dx\) рдХреЗ рд░реВрдк рдореЗрдВ рд╣реИ рдФрд░ рдЖрдк рдХреЗрд╡рд▓ рдЧреБрдгрдХ f рдЯрд╛рдЗрдк рдХрд░рдирд╛ рдЪрд╛рд╣рддреЗ рд╣реИрдВ, рддреЛ f(x) рдЪреБрдиреЗрдВред рдпрджрд┐ рдЖрдкрдХреЗ рдкрд╛рд╕ (0, рдЕрдирдВрдд) рдкрд░ рдкреВрд░рд╛ рд╕рдорд╛рдХрд▓реНрдп g(x) рд╣реИ, рддреЛ g(x) рдЪреБрдиреЗрдВ; рдРрд╕реА рд╕реНрдерд┐рддрд┐ рдореЗрдВ рдЯреВрд▓ рдЕрдВрддрд░реНрдирд┐рд╣рд┐рдд рднрд╛рд░ (built-in weight) рдХреЛ рд╕реНрд╡рддрдГ рд╣рдЯрд╛ рджреЗрддрд╛ рд╣реИред рдлрд▓рди рдХреЛ рдЪрд░ x рдХреЗ рд░реВрдк рдореЗрдВ рдорд╛рдирдХ рд╕рдВрдХреЗрддрди (+, -, *, /, ^, sqrt, exp, ln, sin, cos, tan, рдЖрджрд┐) рдореЗрдВ рджрд░реНрдЬ рдХрд░реЗрдВ, рдиреЛрдбреНрд╕ рдХреА рд╕рдВрдЦреНрдпрд╛ n рд╕реЗрдЯ рдХрд░реЗрдВ, рдФрд░ рднрд╛рд░ рдкреИрд░рд╛рдореАрдЯрд░ alpha рдирд┐рд░реНрдзрд╛рд░рд┐рдд рдХрд░реЗрдВ (рд╕рд╛рдорд╛рдиреНрдп рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░реЗ рдХреЗ рд▓рд┐рдП 0 рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░реЗрдВ)ред n рдмрдврд╝рд╛рдиреЗ рд╕реЗ рд╕рд╣рдЬ рдлрд▓рдиреЛрдВ рдХреЗ рд▓рд┐рдП рд╕рдЯреАрдХрддрд╛ рдмреЗрд╣рддрд░ рд╣реЛрддреА рд╣реИред

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

рдиреЛрдбреНрд╕ \(x_i\), рд╕рд╛рдорд╛рдиреНрдпреАрдХреГрдд рд▓рд╛рдЧреБрдПрд░реЗ рдмрд╣реБрдкрдж \(L_n^{(\alpha)}(x)\) рдХреЗ рдореВрд▓ (roots) рд╣реИрдВ, рдФрд░ рднрд╛рд░ \(w_i\) рдХреЛ рдЧреЛрд▓рдм-рд╡реЗрд▓реНрд╢ (Golub-Welsch) рдПрд▓реНрдЧреЛрд░рд┐рдердо рд╕реЗ рдкреНрд░рд╛рдкреНрдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ: рдПрдХ рд╕рдордорд┐рдд рддреНрд░рд┐рд╡рд┐рдХрд░реНрдгреАрдп (symmetric tridiagonal) рдЬреИрдХреЛрдмреА рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЗ рдЖрдЗрдЧрди-рдорд╛рди (eigenvalues) рдиреЛрдбреНрд╕ рджреЗрддреЗ рд╣реИрдВ, рдЬрдмрдХрд┐ \(w_i = \Gamma(\alpha+1)\) рдЧреБрдгрд╛ рдкреНрд░рддреНрдпреЗрдХ рд╕рд╛рдорд╛рдиреНрдпреАрдХреГрдд рдЖрдЗрдЧрди-рд╕рджрд┐рд╢ (normalized eigenvector) рдХреЗ рдкрд╣рд▓реЗ рдШрдЯрдХ рдХрд╛ рд╡рд░реНрдЧ рд╣реЛрддрд╛ рд╣реИред рдлрд┐рд░ рд╕рдорд╛рдХрд▓рди рдХрд╛ рдЕрдиреБрдорд╛рди рдКрдкрд░ рджрд░реНрд╢рд╛рдП рдЧрдП рднрд╛рд░рд┐рдд рдпреЛрдЧ рд╕реЗ рд▓рдЧрд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ:

$$\int_{0}^{\infty} x^{\alpha} e^{-x}\, g(x)\, dx \approx \sum_{i=1}^{n} w_i\, g(x_i)$$
рдЕрдиреБрдЪрд┐рдд рд╕рдорд╛рдХрд▓ рдХреЛ рдиреЛрдбреНрд╕ рдкрд░ рдлрд▓рди рдорд╛рдиреЛрдВ рдХреЗ рднрд╛рд░рд┐рдд рдпреЛрдЧ рджреНрд╡рд╛рд░рд╛ рд╕рдиреНрдирд┐рдХрдЯрд┐рдд рдХрд┐рдпрд╛ рдЧрдпрд╛
рд╕рдорд╛рдХрд▓ рдХреЛ рдиреЛрдбреНрд╕ \(x_i\) рдкрд░ f рдХреЗ рдорд╛рдиреЛрдВ рдХреЗ рднрд╛рд░рд┐рдд рдкрд░рд┐рдорд┐рдд рдпреЛрдЧ рд╕реЗ рдмрджрд▓ рджрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рднрд╛рд░ \(w_i\) рд╣реИрдВред

рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг

рдорд╛рди рд▓реАрдЬрд┐рдП \(\alpha = 0\), \(n = 2\), рдореЛрдб f(x), \(f(x) = x^2\) (рдЕрд░реНрдерд╛рдд рд╣рдо \(\int x^2 e^{-x}\, dx\) рдХрд╛ рдЕрдиреБрдорд╛рди рд▓рдЧрд╛ рд░рд╣реЗ рд╣реИрдВ, рдЬрд┐рд╕рдХрд╛ рд╕рдЯреАрдХ рдорд╛рди \(\Gamma(3) = 2\) рд╣реИ)ред рджреЛрдиреЛрдВ рдиреЛрдбреНрд╕ рд╣реИрдВ \(x_1 = 2 - \sqrt{2} = 0.585786\) рдЬрд┐рд╕рдХрд╛ \(w_1 = 0.853553\), рдФрд░ \(x_2 = 2 + \sqrt{2} = 3.414214\) рдЬрд┐рд╕рдХрд╛ \(w_2 = 0.146447\)ред рдпреЛрдЧ рд╣реИ $$0.853553 \times 0.343146 + 0.146447 \times 11.656854 = 0.292893 + 1.707107 = 2.000000,$$ рдЬреЛ рд╕рдЯреАрдХ рдЙрддреНрддрд░ рд╕реЗ рдмрд┐рд▓реНрдХреБрд▓ рдореЗрд▓ рдЦрд╛рддрд╛ рд╣реИ рдХреНрдпреЛрдВрдХрд┐ \(x^2\) рдШрд╛рдд 2 рдХрд╛ рдмрд╣реБрдкрдж рд╣реИ рдФрд░ \(2 \le 2n-1 = 3\) рд╣реИред

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

рдкреИрд░рд╛рдореАрдЯрд░ alpha рдХрд╛ рдХреНрдпрд╛ рдХрд╛рдо рд╣реИ? рдпрд╣ рднрд╛рд░ \(x^{\alpha} e^{-x}\) рдореЗрдВ рдШрд╛рддрд╛рдВрдХ рддрдп рдХрд░рддрд╛ рд╣реИред рдорд╛рдирдХ рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░реЗ рдХреЗ рд▓рд┐рдП \(\alpha = 0\) рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░реЗрдВред рдорд╛рди \(\alpha > -1\) рд╣реЛрдирд╛ рдЬрд╝рд░реВрд░реА рд╣реИ рддрд╛рдХрд┐ рднрд╛рд░ рд╕рдорд╛рдХрд▓рдиреАрдп (integrable) рдмрдирд╛ рд░рд╣реЗред

рдореЗрд░рд╛ рдкрд░рд┐рдгрд╛рдо рд╕рдЯреАрдХ рдХреНрдпреЛрдВ рдирд╣реАрдВ рд╣реИ? рдпрд╛ рддреЛ рд╕рдорд╛рдХрд▓реНрдп рд╕рд╣рдЬ рдирд╣реАрдВ рд╣реИ, рдпрд╛ рд╡рд╣ (0, рдЕрдирдВрдд) рдкрд░ рдкрд░реНрдпрд╛рдкреНрдд рддреЗрдЬрд╝реА рд╕реЗ рдирд╣реАрдВ рдШрдЯрддрд╛ред рдпрд╣ рдирд┐рдпрдо рд╣рдореЗрд╢рд╛ рдПрдХ рдкрд░рд┐рдорд┐рдд рд╕рдВрдЦреНрдпрд╛ рд▓реМрдЯрд╛рддрд╛ рд╣реИ, рд▓реЗрдХрд┐рди рд╡рд╣ рддрднреА рд╕рд╛рд░реНрдердХ рд╣реЛрддреА рд╣реИ рдЬрдм рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдорд╛рдХрд▓рди рдЕрднрд┐рд╕рд░рд┐рдд (converge) рд╣реЛ рдФрд░ рд╕рдорд╛рдХрд▓реНрдп рдХреЛ рднрд╛рд░ рдЧреБрдгрд╛ рдмрд╣реБрдкрджреЛрдВ рджреНрд╡рд╛рд░рд╛ рдЕрдЪреНрдЫреА рддрд░рд╣ рд╕реЗ рдЕрдиреБрдорд╛рдирд┐рдд рдХрд┐рдпрд╛ рдЬрд╛ рд╕рдХреЗред рдЕрднрд┐рд╕рд░рдг рдЬрд╛рдБрдЪрдиреЗ рдХреЗ рд▓рд┐рдП n рдмрдврд╝рд╛рдПрдБред

f рдФрд░ g рдореЛрдб рдореЗрдВ рдХреНрдпрд╛ рдЕрдВрддрд░ рд╣реИ? f рдореЛрдб рдореЗрдВ рдЖрдк рдХреЗрд╡рд▓ рд╡рд╣ рдЧреБрдгрдХ рджреЗрддреЗ рд╣реИрдВ рдЬреЛ рдЕрдВрддрд░реНрдирд┐рд╣рд┐рдд рднрд╛рд░ рд╕реЗ рдЧреБрдгрд╛ рд╣реЛрддрд╛ рд╣реИ; g рдореЛрдб рдореЗрдВ рдЖрдк рдкреВрд░рд╛ рд╕рдорд╛рдХрд▓реНрдп рджреЗрддреЗ рд╣реИрдВ рдФрд░ рдпреЛрдЧ рдХреЗ рднреАрддрд░ рднрд╛рд░ рд╣рдЯрд╛ рджрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред рджреЛрдиреЛрдВ рдХреЛ рд╕рд╣реА рдврдВрдЧ рд╕реЗ рд╕реЗрдЯ рдХрд░рдиреЗ рдкрд░ рдПрдХ рд╣реА рдЙрддреНрддрд░ рдорд┐рд▓рддрд╛ рд╣реИред

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