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

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

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

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

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

рд╕рдорд╛рдХрд▓рди рдХрд╛ рдорд╛рди
0.7468241328
[a, b] рдкрд░ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓рди рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯ рдорд╛рди
рд╡рд┐рдзрд┐ Tanh-Sinh (рдбрдмрд▓ рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓)
рдЙрдкрд╡рд┐рднрд╛рдЬрди n 100
рдЙрдкрдпреЛрдЧ рдХрд┐рдП рдЧрдП рдлрд▓рди рдореВрд▓реНрдпрд╛рдВрдХрди 201

Tanh-Sinh рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХреНрдпрд╛ рд╣реИ?

Tanh-Sinh рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░, рдЬрд┐рд╕реЗ рдбрдмрд▓ рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ (DE) рдирд┐рдпрдо рднреА рдХрд╣рд╛ рдЬрд╛рддрд╛ рд╣реИ, рдПрдХ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕рдорд╛рдХрд▓рди (numerical integration) рд╡рд┐рдзрд┐ рд╣реИ рдЬреЛ рдХрд┐рд╕реА рдкрд░рд┐рдорд┐рдд рдЕрдВрддрд░рд╛рд▓ [a, b] рдкрд░ рдирд┐рд╢реНрдЪрд┐рдд рд╕рдорд╛рдХрд▓рди рдирд┐рдХрд╛рд▓рдиреЗ рдореЗрдВ рдмреЗрд╣рдж рдХрд╛рд░рдЧрд░ рд╣реИтАФрдЦрд╛рд╕рдХрд░ рддрдм рдЬрдм рдЗрдВрдЯреАрдЧреНрд░реИрдВрдб рдХреЗ рд╕рд┐рд░реЛрдВ (endpoints) рдкрд░ рд╕рд┐рдВрдЧреБрд▓реИрд░рд┐рдЯреА рд╣реЛред рдпрд╣ рд╡рд┐рдзрд┐ рдЪрд░ рдкрд░рд┐рд╡рд░реНрддрди \( u = \tanh\!\left(\tfrac{\pi}{2}\sinh(t)\right) \) рд▓рд╛рдЧреВ рдХрд░рддреА рд╣реИ, рдЬреЛ рджреЛрдиреЛрдВ рд╕рд┐рд░реЛрдВ рдХреЛ \( t = \pm\infty \) рдкрд░ рднреЗрдЬ рджреЗрддрд╛ рд╣реИред рдЗрди рд╕рд┐рд░реЛрдВ рдХреЗ рдкрд╛рд╕ рдЗрдВрдЯреАрдЧреНрд░реИрдВрдб рдХрд╛ рдпреЛрдЧрджрд╛рди рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ рджрд░ рд╕реЗ рдШрдЯрддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рдРрд╕реЗ рдлрд▓рди рднреА рдЬреЛ a рдпрд╛ b рдкрд░ рдЕрдирдВрдд рд╣реЛ рдЬрд╛рддреЗ рд╣реИрдВ (рдЬреИрд╕реЗ \( 1/\sqrt{1-x^2} \)) рд╕рдЯреАрдХрддрд╛ рд╕реЗ рд╕рдорд╛рдХрд▓рд┐рдд рдХрд┐рдП рдЬрд╛ рд╕рдХрддреЗ рд╣реИрдВред рдпрд╣ рдЯреВрд▓ рд╕рд╛рд░реНрд╡рднреМрдорд┐рдХ рдЧрдгрд┐рдд рдкрд░ рдЖрдзрд╛рд░рд┐рдд рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рд▓рд╛рдЧреВ рд╣реЛрддрд╛ рд╣реИред

S-рдЖрдХрд╛рд░ рдХрд╛ tanh-sinh рд░реВрдкрд╛рдВрддрд░рдг рд╡рдХреНрд░ рдЬреЛ рдЛрдг рдПрдХ рдФрд░ рдзрди рдПрдХ рдкрд░ рдЕрдирдВрддрд╕реНрдкрд░реНрд╢рд┐рдпреЛрдВ рдХреА рдУрд░ рд╕рдкрд╛рдЯ рд╣реЛрддрд╛ рд╣реИ
рджреНрд╡рд┐-рдШрд╛рддрд╛рдВрдХреА рдЪрд░ рдкрд░рд┐рд╡рд░реНрддрди рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд░реЗрдЦрд╛ рдХреЛ (-1, 1) рдкрд░ рдореИрдк рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдиреЛрдбреНрд╕ рд╕рд┐рд░реЛрдВ рдХреА рдУрд░ рдПрдХрддреНрд░рд┐рдд рд╣реЛрддреЗ рд╣реИрдВред

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

рдЕрдкрдирд╛ рдлрд▓рди \( f(x) \) рдорд╛рдирдХ рд╕рдВрдХреЗрддрди рдореЗрдВ рджрд░реНрдЬ рдХрд░реЗрдВ (рдСрдкрд░реЗрдЯрд░ + - * / ^, рдХреЛрд╖реНрдардХ, рдФрд░ sin, cos, exp, log, sqrt, abs рдЬреИрд╕реЗ рдлрд▓рди, рд╕рд╛рде рд╣реА рд╕реНрдерд┐рд░рд╛рдВрдХ pi рдФрд░ e)ред рдирд┐рдЪрд▓реА рд╕реАрдорд╛ a, рдКрдкрд░реА рд╕реАрдорд╛ b, рдФрд░ рдЙрдкрд╡рд┐рднрд╛рдЬрдиреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ n рджреЗрдВ, рдЬреЛ рдиреЛрдб рдШрдирддреНрд╡ рдХреЛ рдирд┐рдпрдВрддреНрд░рд┐рдд рдХрд░рддреА рд╣реИред рдмрдбрд╝рд╛ n рд╕рдЯреАрдХрддрд╛ рдмрдврд╝рд╛рддрд╛ рд╣реИ рдкрд░ рдЧрдгрдирд╛ рдХреА рд▓рд╛рдЧрдд рднреА рдмрдврд╝рд╛рддрд╛ рд╣реИ; рд╡реНрдпрд╡рд╣рд╛рд░ рдореЗрдВ 50тАУ400 рдХрд╛ рджрд╛рдпрд░рд╛ рд╕рд╛рдорд╛рдиреНрдп рд╣реИред рдЗрдВрдЯреАрдЧреНрд░реИрдВрдб рдЦреБрд▓реЗ рдЕрдВрддрд░рд╛рд▓ рдкрд░ рд╡рд┐рд╢реНрд▓реЗрд╖реА (analytic) рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП (рд╕рд┐рд░реЛрдВ рдкрд░ рд╕рд┐рдВрдЧреБрд▓реИрд░рд┐рдЯреА рдареАрдХ рд╣реИ) рдФрд░ рдЖрд╡рд░реНрддреА (periodic) рдирд╣реАрдВ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдПред

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

рд╕рдмрд╕реЗ рдкрд╣рд▓реЗ рдЕрдВрддрд░рд╛рд▓ рдХреЛ \( x = \tfrac{b-a}{2}u + \tfrac{a+b}{2} \) рдФрд░ \( dx = \tfrac{b-a}{2}\,du \) рдХреЗ рджреНрд╡рд╛рд░рд╛ [-1, 1] рдкрд░ рд╕рд╛рдорд╛рдиреНрдпреАрдХреГрдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред рдлрд┐рд░ DE рдирд┐рдпрдо рдиреЛрдб \( t_k = k h \), рднреБрдЬ (abscissa) \( u_k = \tanh\!\left(\tfrac{\pi}{2}\sinh(t_k)\right) \) рдФрд░ рднрд╛рд░ (weights) \( w_k = \dfrac{\tfrac{\pi}{2}\cosh(t_k)}{\cosh^{2}\!\left(\tfrac{\pi}{2}\sinh(t_k)\right)} \) рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред рд╕рдорд╛рдХрд▓рди рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯрди рдирд┐рдореНрди рд╕реВрддреНрд░ рдХреЗ рдпреЛрдЧ рд╕реЗ рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ:

$$\int_{a}^{b} f(x)\,dx \approx \frac{b-a}{2}\, h \sum_{k=-N}^{N} w_k\, f\!\left(\frac{a+b}{2} + \frac{b-a}{2}\,x_k\right)$$

рдЬрд┐рди рдиреЛрдб рдХрд╛ рднрд╛рд░ рд╢реВрдиреНрдп рддрдХ рдШрдЯ (underflow) рдЬрд╛рддрд╛ рд╣реИ (рд╕рдВрддреГрдкреНрдд рд╕рд┐рд░реЗ) рдЙрдиреНрд╣реЗрдВ рдЫреЛрдбрд╝ рджрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рд╕рдВрднрд╛рд╡рд┐рдд рд╕рд┐рдВрдЧреБрд▓рд░ рд╕реАрдорд╛рдУрдВ рдкрд░ рдореВрд▓реНрдпрд╛рдВрдХрди рд╕реЗ рдмрдЪрд╛ рдЬрд╛ рд╕рдХреЗред

a рд╕реЗ b рддрдХ рдХрд╛ рдЕрдВрддрд░рд╛рд▓, рдЬрд┐рд╕рдореЗрдВ рдЫрд╛рдпрд╛рдВрдХрд┐рдд рдХреНрд╖реЗрддреНрд░ рдХреЗ рдиреАрдЪреЗ рджреЛрдиреЛрдВ рд╕рд┐рд░реЛрдВ рдХреЗ рдкрд╛рд╕ рд╕рдШрди рд░реВрдк рд╕реЗ рдПрдХрддреНрд░рд┐рдд рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдиреЛрдбреНрд╕ рд╣реИрдВ
рдиреЛрдбреНрд╕ рд╕рд┐рд░реЛрдВ a рдФрд░ b рдХреЗ рдкрд╛рд╕ рдПрдХрддреНрд░рд┐рдд рд╣реЛрддреЗ рд╣реИрдВ, рдЗрд╕рд▓рд┐рдП рд╕рд┐рд░реЛрдВ рдХреА рд╡рд┐рдЪрд┐рддреНрд░рддрд╛рдПрдБ рд╕рдЯреАрдХрддрд╛ рд╕реЗ рд╕рдВрднрд╛рд▓реА рдЬрд╛рддреА рд╣реИрдВред

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

\( f(x) = e^{-x^2} \) рдХреЛ [0, 1] рдкрд░ рд╕рдорд╛рдХрд▓рд┐рдд рдХрд░реЗрдВред рд╕рдЯреАрдХ рдорд╛рди $$\frac{\sqrt{\pi}}{2}\cdot\operatorname{erf}(1) \approx 0.7468241$$ рд╣реИред рдореЛрдЯреЗ рдЪрд░рдг (\( h = 0.5 \), \( N = 4 \)) рдкрд░ рднреА рдпрд╣ рдирд┐рдпрдо рд▓рдЧрднрдЧ 0.7467 рджреЗ рджреЗрддрд╛ рд╣реИ; рдФрд░ рдбрд┐рдлрд╝реЙрд▓реНрдЯ \( n = 100 \) рдХреЗ рд╕рд╛рде рдпрд╣ рдХрд░реАрдм рдмрд╛рд░рд╣ рдЕрдВрдХреЛрдВ рддрдХ рд╕рд╣реА рдмреИрдарддрд╛ рд╣реИред

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

рдХреНрдпрд╛ рдореИрдВ рд╕рд┐рд░реЛрдВ рдкрд░ рд╕рд┐рдВрдЧреБрд▓реИрд░рд┐рдЯреА рд╡рд╛рд▓реЗ рдлрд▓рди рд╕рдорд╛рдХрд▓рд┐рдд рдХрд░ рд╕рдХрддрд╛ рд╣реВрдБ? рд╣рд╛рдБтАФрдпрд╣реА рдЗрд╕ рд╡рд┐рдзрд┐ рдХреА рд╕рдмрд╕реЗ рдмрдбрд╝реА рдЦреВрдмреА рд╣реИред a рдпрд╛ b рдкрд░ рд╕рдорд╛рдХрд▓рдиреАрдп рд╕рд┐рдВрдЧреБрд▓реИрд░рд┐рдЯреА рдХреЛ рдпрд╣ рд╕рд╣рдЬрддрд╛ рд╕реЗ рд╕рдВрднрд╛рд▓ рд▓реЗрддреА рд╣реИред

рдЖрд╡рд░реНрддреАрдпрддрд╛ (periodicity) рдХреНрдпреЛрдВ рдорд╛рдпрдиреЗ рд░рдЦрддреА рд╣реИ? рдбрдмрд▓-рдПрдХреНрд╕рдкреЛрдиреЗрдВрд╢рд┐рдпрд▓ рдирд┐рдпрдо рдЧреИрд░-рдЖрд╡рд░реНрддреА рдЗрдВрдЯреАрдЧреНрд░реИрдВрдб рдХреЗ рд▓рд┐рдП рдмрдирд╛ рд╣реИ; рдЖрд╡рд░реНрддреА рдлрд▓рдиреЛрдВ рдХреЗ рд▓рд┐рдП рдЯреНрд░реИрдкреЗрдЬрд╝реЙрдЗрдбрд▓ рдирд┐рдпрдо рддреЗрдЬрд╝реА рд╕реЗ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛрддрд╛ рд╣реИ рдФрд░ DE рдЧрд▓рдд рд╣реЛ рд╕рдХрддрд╛ рд╣реИред

рдЕрдЧрд░ a, b рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛ рддреЛ? рд╕рдорд╛рдХрд▓рди рд╢реВрдиреНрдп рд╣реЛрдЧрд╛ред рдЕрдЧрд░ a > b рд╣реЛ рддреЛ рдкрд░рд┐рдгрд╛рдо рд╕рд╣реА рдврдВрдЧ рд╕реЗ рдЛрдгрд╛рддреНрдордХ рд╣реЛ рдЬрд╛рддрд╛ рд╣реИ, рдХреНрдпреЛрдВрдХрд┐ \( \tfrac{b-a}{2} \) рдЧреБрдгрд╛рдВрдХ рдЪрд┐рд╣реНрди рдХреЛ рдЕрдкрдиреЗ рд╕рд╛рде рд▓реЗ рдЬрд╛рддрд╛ рд╣реИред

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