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

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

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

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

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

рдореВрд▓ x, рдЬрд╣рд╛рдБ f(x) = 0
0.73908513321516
Converged
рдЙрдкрдпреЛрдЧ рдХреА рдЧрдИ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпрд╛рдБ 5
рдЕрд╡рд╢реЗрд╖ f(x) 0.000000000000003
рд╡рд┐рдзрд┐ рдмрд┐рдирд╛ рдбреЗрд░рд┐рд╡реЗрдЯрд┐рд╡ рд╡рд╛рд▓реА рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ (рд╕реНрдЯреЗрдлреЗрдиреНрд╕рди рд╢реИрд▓реА)

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

рдпрд╣ рдЯреВрд▓ рдХрд┐рд╕реА рд╕рдореАрдХрд░рдг рдХрд╛ рдореВрд▓ (root) рдЦреЛрдЬрддрд╛ рд╣реИ тАФ рдпрд╛рдиреА x рдХрд╛ рд╡рд╣ рдорд╛рди рдЬрд┐рд╕ рдкрд░ \(f(x) = 0\) рд╣реЛрддрд╛ рд╣реИ тАФ рдФрд░ рдЗрд╕рдХреЗ рд▓рд┐рдП рдпрд╣ рдмрд┐рдирд╛ рдбреЗрд░рд┐рд╡реЗрдЯрд┐рд╡ рд╡рд╛рд▓реА рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕реЗ рд╕реНрдЯреЗрдлреЗрдиреНрд╕рди рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ (Steffensen's iteration) рднреА рдХрд╣рддреЗ рд╣реИрдВред рдкрд╛рд░рдВрдкрд░рд┐рдХ рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ рдореЗрдВ рдкрд╣рд▓реЗ рдЕрд╡рдХрд▓рдЬ (first derivative) \(f'(x)\) рдХреА рдЬрд╝рд░реВрд░рдд рд╣реЛрддреА рд╣реИ; рд▓реЗрдХрд┐рди рдпрд╣ рд╕рдВрд╕реНрдХрд░рдг рдбреЗрд░рд┐рд╡реЗрдЯрд┐рд╡ рдХреА рдЬрдЧрд╣ рд╕рд┐рд░реНрдлрд╝ рдлрд▓рди рдХреЗ рдорд╛рдиреЛрдВ рд╕реЗ рдмрдиреА рдлреЙрд░рд╡рд░реНрдб-рдбрд┐рдлрд░реЗрдВрд╕ рдЕрдиреБрдорд╛рди рдХрд╛ рдЗрд╕реНрддреЗрдорд╛рд▓ рдХрд░рддрд╛ рд╣реИ, рдЗрд╕рд▓рд┐рдП рдЖрдкрдХреЛ рд╣рд╛рде рд╕реЗ рдЕрд╡рдХрд▓рди рдХрд░рдиреЗ рдХреА рдХреЛрдИ рдЬрд╝рд░реВрд░рдд рдирд╣реАрдВ рдкрдбрд╝рддреАред рдпрд╣ рдПрдХ рд╢реБрджреНрдз рдЧрдгрд┐рддреАрдп рдЯреВрд▓ рд╣реИ рдЬреЛ рдХрд┐рд╕реА рднреА рдПрдХрд▓-рдЪрд░ (single-variable) рд╡рд╛рд╕реНрддрд╡рд┐рдХ рдлрд▓рди рдХреЗ рд▓рд┐рдП рдХрд╛рдо рдХрд░рддрд╛ рд╣реИред

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

рдЕрдкрдиреЗ рдлрд▓рди рдХреЛ f(x) рдмреЙрдХреНрд╕ рдореЗрдВ рдЪрд░ x рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП рджрд░реНрдЬ рдХрд░реЗрдВред рдЖрдк + - * / ^, рдХреЛрд╖реНрдардХ (parentheses), рдЕрдЪрд░ pi рдФрд░ e, рддрдерд╛ рд╕рд╛рдорд╛рдиреНрдп рдлрд▓рди рдЙрдкрдпреЛрдЧ рдХрд░ рд╕рдХрддреЗ рд╣реИрдВ: sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, exp, ln, log, log10, sqrt, absред рддреНрд░рд┐рдХреЛрдгрдорд┐рддреАрдп рдлрд▓рди рд░реЗрдбрд┐рдпрди рдореЗрдВ рдХрд╛рдо рдХрд░рддреЗ рд╣реИрдВред рд╢реБрд░реБрдЖрддреА рдЕрдиреБрдорд╛рди x0 рд╕реЗрдЯ рдХрд░реЗрдВ (рдкрд░рд┐рдгрд╛рдо рдЗрд╕реА рдкрд░ рдирд┐рд░реНрднрд░ рдХрд░рддрд╛ рд╣реИ), рдЕрдзрд┐рдХрддрдо рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпреЛрдВ рдХреА рд╕рдВрдЦреНрдпрд╛ n рдЪреБрдиреЗрдВ, рдФрд░ рдЕрднрд┐рд╕рд░рд┐рдд (converged) рдореВрд▓, рдЕрд╡рд╢реЗрд╖ (residual) f(x), рддрдерд╛ рдХрд┐рддрдиреА рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпрд╛рдБ рд▓рдЧреАрдВ тАФ рдпрд╣ рд╕рдм рдкрдврд╝реЗрдВред

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

рдЕрджреНрдпрддрди рдирд┐рдпрдо (update rule) рдпрд╣ рд╣реИ: $$x_{n+1} = x_n - \frac{f(x_n)^2}{f\!\left(x_n + f(x_n)\right) - f(x_n)}$$ рдпрд╣ рдорд╛рдирдХ рдиреНрдпреВрдЯрди рдЪрд░рдг \(x_{n+1} = x_n - f(x_n)/f'(x_n)\) рд╕реЗ рдЖрддрд╛ рд╣реИ, рдЬрд╣рд╛рдБ \(f'(x_n)\) рдХреЛ \(h = f(x_n)\) рд╕реНрдЯреЗрдк рд╕рд╛рдЗрдЬрд╝ рд╡рд╛рд▓реА рдлреЙрд░рд╡рд░реНрдб рдбрд┐рдлрд░реЗрдВрд╕ рд╕реЗ рдЕрдиреБрдорд╛рдирд┐рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред рдЗрд╕ рдЕрдиреБрдорд╛рди рдХреЛ рд░рдЦрдиреЗ рдкрд░ рдКрдкрд░ рд╡рд╛рд▓рд╛ рд╕реВрддреНрд░ рдорд┐рд▓рддрд╛ рд╣реИред рдХрд┐рд╕реА рд╕рд░рд▓ рдореВрд▓ (simple root) рдХреЗ рдкрд╛рд╕ рдпрд╣ рд▓рдЧрднрдЧ рд╡рд░реНрдЧрд╛рдХрд╛рд░ (quadratically) рдЧрддрд┐ рд╕реЗ рдЕрднрд┐рд╕рд░рд┐рдд рд╣реЛрддрд╛ рд╣реИ тАФ рдЕрд╕рд▓реА рдиреНрдпреВрдЯрди рд╡рд┐рдзрд┐ рдХреА рддрд░рд╣ тАФ рдкрд░рдВрддреБ рдХреЗрд╡рд▓ рдлрд▓рди рдХреЗ рдорд╛рдиреЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред

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

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

\(f(x) = x - \cos(x)\) рдФрд░ \(x_0 = 1\) рдХреЗ рд▓рд┐рдП: рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ 1 рдореЗрдВ \(f(1) = 0.45970\), рдкреНрд░реЛрдм \(f(1.45970) = 1.34861\), рд╣рд░ (denominator) \(0.88891\), рдЕрддрдГ $$x_1 = 1 - \frac{0.45970^2}{0.88891} = 0.76224$$ рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рдЬрд▓реНрджреА рд╣реА \(x \approx 0.7390851332\) рдкрд░ рд╕реНрдерд┐рд░ рд╣реЛ рдЬрд╛рддреА рд╣реИ тАФ рдпрд╣реА рд╡рд╣ рдорд╛рди рд╣реИ рдЬрд╣рд╛рдБ \(x = \cos(x)\) рд╣реЛрддрд╛ рд╣реИ (рдбреЙрдЯреА рд╕рдВрдЦреНрдпрд╛ / Dottie number), рдФрд░ \(f(x) \approx 0\) рд░рд╣рддрд╛ рд╣реИред

рдореВрд▓ рдЦреЛрдЬрдиреЗ рд╡рд╛рд▓реЗ рдкреБрдирд░рд╛рд╡реГрддреНрдд рд▓реВрдк рдХрд╛ рдкреНрд░рд╡рд╛рд╣ рдЖрд░реЗрдЦ
рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐ рддрдм рддрдХ рджреЛрд╣рд░рд╛рдИ рдЬрд╛рддреА рд╣реИ рдЬрдм рддрдХ x рдореЗрдВ рдмрджрд▓рд╛рд╡ рдмрд╣реБрдд рдЫреЛрдЯрд╛ рди рд╣реЛ рдпрд╛ рдЕрдзрд┐рдХрддрдо рдкреБрдирд░рд╛рд╡реГрддреНрддрд┐рдпрд╛рдБ рдкреВрд░реА рди рд╣реЛ рдЬрд╛рдПрдБред

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

рдпрд╣ рдЕрднрд┐рд╕рд░рд┐рдд (converge) рдХреНрдпреЛрдВ рдирд╣реАрдВ рд╣реЛ рдкрд╛рддрд╛? рдЦрд░рд╛рдм рд╢реБрд░реБрдЖрддреА рдЕрдиреБрдорд╛рди рд╕реЗ рдпрд╣ рд╡рд┐рд╕рд░рд┐рдд (diverge) рд╣реЛ рд╕рдХрддрд╛ рд╣реИ рдпрд╛ рдХрд┐рд╕реА рджреВрд╕рд░реЗ рдореВрд▓ рдкрд░ рдкрд╣реБрдБрдЪ рд╕рдХрддрд╛ рд╣реИ, рдФрд░ рдпрджрд┐ рд╣рд░ \(f(x+f(x)) - f(x)\) рд╢реВрдиреНрдп рд╣реЛ рдЬрд╛рдП рддреЛ рд╡рд┐рдзрд┐ рд╕реБрд░рдХреНрд╖рд┐рдд рд░реВрдк рд╕реЗ рд░реБрдХ рдЬрд╛рддреА рд╣реИред рдРрд╕реЗ рдореЗрдВ рдХреЛрдИ рджреВрд╕рд░рд╛ x0 рдЖрдЬрд╝рдорд╛рдПрдБред

рдХреЛрдг рдбрд┐рдЧреНрд░реА рдореЗрдВ рд╣реЛрддреЗ рд╣реИрдВ рдпрд╛ рд░реЗрдбрд┐рдпрди рдореЗрдВ? рд░реЗрдбрд┐рдпрди рдореЗрдВ тАФ рдпрд╣реА рдЧрдгрд┐рдд рдХреА рдорд╛рдирдХ рдкрд░рд┐рдкрд╛рдЯреА рд╣реИред рдЬрд╝рд░реВрд░рдд рдкрдбрд╝рдиреЗ рдкрд░ x*pi/180 рд╕реЗ рд░реВрдкрд╛рдВрддрд░рдг рдХрд░реЗрдВред

рдХреНрдпрд╛ рдореБрдЭреЗ рдбреЗрд░рд┐рд╡реЗрдЯрд┐рд╡ рдХреА рдЬрд╝рд░реВрд░рдд рд╣реИ? рдирд╣реАрдВ тАФ рдпрд╣реА рддреЛ рдЗрд╕ рд╡рд┐рдзрд┐ рдХреА рдЦрд╛рд╕ рдмрд╛рдд рд╣реИред рдпрд╣ рдврд▓рд╛рди (slope) рдХрд╛ рдЕрдиреБрдорд╛рди рдХреЗрд╡рд▓ рдлрд▓рди рдХреЗ рдорд╛рдиреЛрдВ рд╕реЗ рд╣реА рд▓рдЧрд╛ рд▓реЗрддреА рд╣реИред

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