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

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

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

рд╕реВрддреНрд░ (рдлреЙрд░реНрдореВрд▓рд╛): рд╣реИрд▓реА рд╡рд┐рдзрд┐ рд░реВрдЯ-рдлрд╛рдЗрдВрдбрд┐рдВрдЧ рдХреИрд▓рдХреБрд▓реЗрдЯрд░

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

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

рдЕрдиреБрдорд╛рдирд┐рдд рдореВрд▓ x (f(x) = 0)
0.7390851332151607
рд╣реИрд▓реА рд╡рд┐рдзрд┐ рд╕реЗ рдкреНрд░рд╛рдкреНрдд
рдкреНрд░рдпреБрдХреНрдд рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпрд╛рдБ 3
рдЕрдВрддрд┐рдо рдЕрд╡рд╢реЗрд╖ f(x) 0.0
рд╕реНрдерд┐рддрд┐ Converged

рд╣реИрд▓реА рд╡рд┐рдзрд┐ рдХреНрдпрд╛ рд╣реИ?

рд╣реИрд▓реА рд╡рд┐рдзрд┐ f(x) = 0 рд░реВрдк рд╡рд╛рд▓реЗ рд╕рдореАрдХрд░рдг рдХреЛ рд╣рд▓ рдХрд░рдиреЗ рдХреА рдПрдХ рдкреБрдирд░рд╛рд╡реГрддреНрддреАрдп (iterative) рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рддрдХрдиреАрдХ рд╣реИред рдпрд╣ рдиреНрдпреВрдЯрди-рд░реИрдлрд╕рди рд╡рд┐рдзрд┐ рдХреА рд╣реА рддреАрд╕рд░реЗ рдХреНрд░рдо рд╡рд╛рд▓реА (рдХреНрдпреВрдмрд┐рдХ-рдЕрднрд┐рд╕рд░рдг) рдХрд░реАрдмреА рдорд╛рдиреА рдЬрд╛рддреА рд╣реИ: рдЬрд╣рд╛рдБ рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ рдХреЗрд╡рд▓ рдлрд▓рди рдФрд░ рдЙрд╕рдХреЗ рдкрд╣рд▓реЗ рдЕрд╡рдХрд▓рдЬ рдХрд╛ рдкреНрд░рдпреЛрдЧ рдХрд░рддреА рд╣реИ, рд╡рд╣реАрдВ рд╣реИрд▓реА рд╡рд┐рдзрд┐ рджреВрд╕рд░реЗ рдЕрд╡рдХрд▓рдЬ рдХреЛ рднреА рдЬреЛрдбрд╝рддреА рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдЖрдорддреМрд░ рдкрд░ рдХрд┐рд╕реА рджреА рдЧрдИ рд╕рдЯреАрдХрддрд╛ рддрдХ рдХрдо рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпреЛрдВ рдореЗрдВ рд╣реА рдкрд╣реБрдБрдЪрд╛ рдЬрд╛ рд╕рдХрддрд╛ рд╣реИред рдпрд╣ рдПрдХ рд╕рд╛рд░реНрд╡рднреМрдорд┐рдХ рдЧрдгрд┐рдд / рд╕рдВрдЦреНрдпрд╛рддреНрдордХ-рд╡рд┐рд╢реНрд▓реЗрд╖рдг рдЙрдкрдХрд░рдг рд╣реИ рдФрд░ рд╣рд░ рдЬрдЧрд╣ рд▓рд╛рдЧреВ рд╣реЛрддрд╛ рд╣реИ; рддреНрд░рд┐рдХреЛрдгрдорд┐рддреАрдп рдлрд▓рдиреЛрдВ рдХреЗ рдЕрдВрджрд░ рджрд┐рдП рдЧрдП рдХреЛрдг рд░реЗрдбрд┐рдпрди рдореЗрдВ рдорд╛рдиреЗ рдЬрд╛рддреЗ рд╣реИрдВред

рд╡рдХреНрд░ f(x) рдПрдХ рдореВрд▓ рдкрд░ x-рдЕрдХреНрд╖ рдХреЛ рдХрд╛рдЯрддрд╛ рд╣реБрдЖ, рдЬрд┐рд╕рдХреЗ рдХреНрд░рдорд┐рдХ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдмрд┐рдВрджреБ рдЙрд╕рдХреА рдУрд░ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛрддреЗ рд╣реИрдВ
рд╣реИрд▓реА рд╡рд┐рдзрд┐ рдПрдХ рдЕрдиреБрдорд╛рди рдХреЛ рдмрд╛рд░-рдмрд╛рд░ рдкрд░рд┐рд╖реНрдХреГрдд рдХрд░рддреА рд╣реИ рдЬрдм рддрдХ рд╡рд╣ рдЙрд╕ рдореВрд▓ рдкрд░ рдЕрднрд┐рд╕рд░рд┐рдд рди рд╣реЛ рдЬрд╛рдП рдЬрд╣рд╛рдБ \(f(x)=0\) рд╣реЛрддрд╛ рд╣реИред

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

рдЬрд┐рд╕ рдлрд▓рди рдХрд╛ рдореВрд▓ рдЖрдк рдЦреЛрдЬрдирд╛ рдЪрд╛рд╣рддреЗ рд╣реИрдВ рдЙрд╕реЗ \(f(x)\) рдХреЗ рд░реВрдк рдореЗрдВ рджрд░реНрдЬ рдХрд░реЗрдВ, рдлрд┐рд░ рдЙрд╕рдХрд╛ рдкрд╣рд▓рд╛ рдЕрд╡рдХрд▓рдЬ \(f'(x)\) рдФрд░ рджреВрд╕рд░рд╛ рдЕрд╡рдХрд▓рдЬ \(f''(x)\) x рдореЗрдВ рд╡реНрдпрдВрдЬрдХ рдХреЗ рд░реВрдк рдореЗрдВ рджреЗрдВред рд╕рдорд░реНрдерд┐рдд рд╕рд┐рдВрдЯреИрдХреНрд╕ рдореЗрдВ рдШрд╛рддреЛрдВ рдХреЗ рд▓рд┐рдП +, -, *, /, ^, рдХреЛрд╖реНрдардХ рддрдерд╛ sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, exp, log/ln, log10, sqrt рдФрд░ abs рдЬреИрд╕реЗ рдлрд▓рди, рдФрд░ pi рд╡ e рдЬреИрд╕реЗ рд╕реНрдерд┐рд░рд╛рдВрдХ рд╢рд╛рдорд┐рд▓ рд╣реИрдВред рдЕрдкрдиреЗ рдЗрдЪреНрдЫрд┐рдд рдореВрд▓ рдХреЗ рдкрд╛рд╕ рдПрдХ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдЕрдиреБрдорд╛рди \(x_0\) рдЪреБрдиреЗрдВ, рдЕрдзрд┐рдХрддрдо рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╕рдВрдЦреНрдпрд╛ \(n\) рддрдп рдХрд░реЗрдВ, рдФрд░ рдЪреБрдиреЗрдВ рдХрд┐ рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рдЕрдВрдХ (significant digits) рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд░рдиреЗ рд╣реИрдВред

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

рд╣рд░ рдЪрд░рдг $$x_{n+1} = x_n - \frac{2\,f(x_n)\,f'(x_n)}{2\,[f'(x_n)]^2 - f(x_n)\,f''(x_n)}$$ рдХреА рдЧрдгрдирд╛ рдХрд░рддрд╛ рд╣реИред рдЕрдВрд╢ (numerator) рд╕рд╛рдорд╛рдиреНрдп рдиреНрдпреВрдЯрди рд╕реБрдзрд╛рд░ (рд╕реНрдХреЗрд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ) рд╣реИ, рдЬрдмрдХрд┐ рд╣рд░ (denominator) рдореЗрдВ рдЕрддрд┐рд░рд┐рдХреНрдд рдкрдж \(- f(x_n) f''(x_n)\) рдлрд▓рди f рдХреА рд╡рдХреНрд░рддрд╛ (curvature) рдХреЗ рд▓рд┐рдП рд╕реБрдзрд╛рд░ рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рд╕реЗ рдЕрднрд┐рд╕рд░рдг рддреЗрдЬрд╝ рд╣реЛрддрд╛ рд╣реИред рдЬрдм x рдореЗрдВ рдкрд░рд┐рд╡рд░реНрддрди рдпрд╛ рдЕрд╡рд╢реЗрд╖ (residual) \(f(x)\) рдХрд┐рд╕реА рдЕрддреНрдпрдВрдд рдЫреЛрдЯреЗ рд╕рд╣рди-рд╕реНрддрд░ рд╕реЗ рдиреАрдЪреЗ рдЖ рдЬрд╛рддрд╛ рд╣реИ, рдпрд╛ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдХреА рд╕реАрдорд╛ рдкреВрд░реА рд╣реЛ рдЬрд╛рддреА рд╣реИ, рддрдм рд▓реВрдк рд░реБрдХ рдЬрд╛рддрд╛ рд╣реИред рдпрджрд┐ рд╣рд░ рд╢реВрдиреНрдп рд╣реЛ рдЬрд╛рдП рддреЛ рд╡рд┐рдзрд┐ рд╡рд┐рдлрд▓ рд╣реЛ рдЬрд╛рддреА рд╣реИ рдФрд░ рдЖрдкрдХреЛ рдХреЛрдИ рдЕрд▓рдЧ \(x_0\) рдЖрдЬрд╝рдорд╛рдирд╛ рдЪрд╛рд╣рд┐рдПред

рдПрдХ рдмрд┐рдВрджреБ рдкрд░ рдлрд▓рди рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯрди рдХрд░рддреА рд╕реНрдкрд░реНрд╢рд░реЗрдЦрд╛ рдФрд░ рд╡рдХреНрд░рддрд╛ рд╡рд╛рд▓реА рд╕рдВрд╕реНрдкрд░реНрд╢реА рд╡рдХреНрд░ рджрд┐рдЦрд╛рддрд╛ рдЖрд░реЗрдЦ
рд╣реИрд▓реА рд╡рд┐рдзрд┐ рдврд▓рд╛рди \(f'(x)\) рдФрд░ рд╡рдХреНрд░рддрд╛ \(f''(x)\) рджреЛрдиреЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреА рд╣реИ, рдФрд░ рдПрдХ рдРрд╕рд╛ рд╡рдХреНрд░ рдлрд┐рдЯ рдХрд░рддреА рд╣реИ рдЬреЛ рд╕реАрдзреА рд╕реНрдкрд░реНрд╢рд░реЗрдЦрд╛ рдХреА рддреБрд▓рдирд╛ рдореЗрдВ рдлрд▓рди рд╕реЗ рдмреЗрд╣рддрд░ рд╕рдЯрддрд╛ рд╣реИред

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

\(f(x) = x - \cos(x)\), \(f'(x) = 1 + \sin(x)\), \(f''(x) = \cos(x)\) рдХреЗ рд▓рд┐рдП, \(x_0 = 1\) рд╕реЗ рд╢реБрд░реВ рдХрд░рддреЗ рд╣реБрдП: рдкрд╣рд▓рд╛ рдЪрд░рдг рджреЗрддрд╛ рд╣реИ $$x_1 = 1 - \frac{2 \times 0.4596977 \times 1.8414710}{2 \times 1.8414710^2 - 0.4596977 \times 0.5403023} = 1 - \frac{1.6930504}{6.5336550} = 0.7408769$$ рдпрд╣ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рд╢реАрдШреНрд░ рд╣реА рдбреЙрдЯреА рд╕рдВрдЦреНрдпрд╛ (Dottie number) \(x = 0.7390851332151607\) рдкрд░ рд╕реНрдерд┐рд░ рд╣реЛ рдЬрд╛рддреА рд╣реИ, рдЬреЛ \(x = \cos(x)\) рдХрд╛ рдПрдХрдорд╛рддреНрд░ рд╣рд▓ рд╣реИред

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

рд╣реИрд▓реА рд╡рд┐рдзрд┐ рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ рд╕реЗ рдХреИрд╕реЗ рдЕрд▓рдЧ рд╣реИ? рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ рд╡рдХреНрд░рддрд╛ рдХреЛ рдЕрдирджреЗрдЦрд╛ рдХрд░рддреА рд╣реИ; рд╣реИрд▓реА рд╡рд┐рдзрд┐ рджреВрд╕рд░реЗ-рдЕрд╡рдХрд▓рдЬ рд╡рд╛рд▓рд╛ рдкрдж рдЬреЛрдбрд╝рддреА рд╣реИ, рдЬрд┐рд╕рд╕реЗ рджреНрд╡рд┐рдШрд╛рдд (quadratic) рдХреЗ рдмрдЬрд╛рдп рдХреНрдпреВрдмрд┐рдХ рдЕрднрд┐рд╕рд░рдг рдорд┐рд▓рддрд╛ рд╣реИ, рдФрд░ рдЖрдорддреМрд░ рдкрд░ рд╣рд░ рд╕рд╣реА рдЕрдВрдХ рдХреЗ рд▓рд┐рдП рдХрдо рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпрд╛рдБ рд▓рдЧрддреА рд╣реИрдВред

рдореБрдЭреЗ рдЕрд╡рдХрд▓рдЬ рдХреНрдпреЛрдВ рджрд░реНрдЬ рдХрд░рдиреЗ рдкрдбрд╝рддреЗ рд╣реИрдВ? рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдЖрдкрдХреЗ рджрд┐рдП рдЧрдП рдЕрд╡рдХрд▓рдЬреЛрдВ рдкрд░ рднрд░реЛрд╕рд╛ рдХрд░рддрд╛ рд╣реИред рдпрджрд┐ рд╡реЗ рдЧрд▓рдд рд╣реЛрдВрдЧреЗ, рддреЛ рдЕрднрд┐рд╕рд░рдг рдХрдордЬрд╝реЛрд░ рд░рд╣реЗрдЧрд╛ рдпрд╛ рдЕрд╕рдлрд▓ рд╣реЛ рдЬрд╛рдПрдЧрд╛ред \(f(x)\) рдХрд╛ рд╕рд╛рд╡рдзрд╛рдиреА рд╕реЗ рдЕрд╡рдХрд▓рди рдХрд░рдХреЗ \(f'(x)\) рдФрд░ \(f''(x)\) рдирд┐рдХрд╛рд▓реЗрдВред

рдЕрдЧрд░ рдЕрднрд┐рд╕рд░рдг рди рд╣реЛ рддреЛ рдХреНрдпрд╛ рдХрд░реЗрдВ? рдпрд╣ рд╡рд┐рдзрд┐ \(x_0\) рдХреЗ рд╕рдмрд╕реЗ рдирд┐рдХрдЯ рд╡рд╛рд▓рд╛ рдореВрд▓ рдЦреЛрдЬрддреА рд╣реИред рдЦрд░рд╛рдм рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдмрд┐рдВрджреБ рд╕реЗ рдкрд░рд┐рдгрд╛рдо рдЕрдкрд╕рд░рд┐рдд (diverge) рд╣реЛ рд╕рдХрддрд╛ рд╣реИ рдпрд╛ рдХрд┐рд╕реА рджреВрд╕рд░реЗ рдореВрд▓ рдкрд░ рдЬрд╛ рд╕рдХрддрд╛ рд╣реИ, рдФрд░ рд╣рд░ рдХреЗ рд╢реВрдиреНрдп рд╣реЛрдиреЗ рдкрд░ рдкреНрд░рдХреНрд░рд┐рдпрд╛ рдкреВрд░реА рддрд░рд╣ рд░реБрдХ рдЬрд╛рддреА рд╣реИред рдРрд╕реЗ рдореЗрдВ \(x_0\) рдмрджрд▓реЗрдВ рдпрд╛ рдЕрдкрдиреЗ рдЕрд╡рдХрд▓рдЬ рд╡реНрдпрдВрдЬрдХреЛрдВ рдХреА рдЬрд╛рдБрдЪ рдХрд░реЗрдВред

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