Метод касательных Ньютона для корня уравнения

Метод касательных Ньютона - самый быстрый из классических способов численно найти корень уравнения $f(x) = 0$. Идея геометрическая: в текущей точке $x_n$ проводят касательную к графику функции и берут за следующее приближение ту точку, где эта касательная пересекает ось абсцисс. Если функция гладкая, а старт выбран достаточно близко к корню, метод Ньютона сходится квадратично - число верных знаков примерно удваивается на каждом шаге, и для машинной точности обычно хватает четырёх-шести итераций. За эту скорость приходится платить: нужна производная $f'(x)$ и аккуратное начальное приближение, иначе итерации могут разойтись или зациклиться. Разберём вывод формулы, геометрический смысл касательной, условие сходимости, скорость и подводные камни.

Идея и итерационная формула

Пусть функция $f$ дифференцируема и мы ищем корень уравнения $f(x) = 0$. Разложим $f$ в окрестности приближения $x_n$ в ряд Тейлора и оставим линейную часть:

$f(x) \approx f(x_n) + f'(x_n)\,(x - x_n).$

Это уравнение касательной к графику в точке $x_n$. Приравняв правую часть к нулю и разрешив относительно $x$, получаем точку пересечения касательной с осью - она и становится следующим приближением:

$x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)}.$

Это и есть итерационная формула метода касательных Ньютона. Каждый шаг требует двух вычислений: значения функции $f(x_n)$ и её производной $f'(x_n)$. Условие применимости очевидно из формулы - производная в точке итерации не должна обращаться в ноль: при $f'(x_n) = 0$ касательная горизонтальна и не пересекает ось.

Пошаговый расчёт корня

Прежде чем разбирать сходимость в общем виде, удобно прогнать конкретное уравнение через сам метод и увидеть, как стремительно сужается погрешность. Ниже - небольшая форма: задаёшь уравнение $f(x) = 0$, начальное приближение $x_0$ и точность $\varepsilon$, а в ответ получаешь производную $f'(x)$, итерационную формулу и пошаговую таблицу метода касательных с найденным корнем.

После того как видно, насколько быстро удваивается число верных знаков, разберём, что именно гарантирует такую сходимость и когда метод Ньютона может подвести.

Геометрический смысл касательной

Название «метод касательных» отражает суть наглядно. На графике функции в точке $(x_n, f(x_n))$ строится касательная прямая с угловым коэффициентом $f'(x_n)$. Она пересекает ось абсцисс в точке $x_{n+1}$, которая, как правило, оказывается ближе к корню $x^{*}$, чем исходная $x_n$. Затем процесс повторяется уже из новой точки: строится новая касательная, находится новое пересечение - и так до нужной точности.

Геометрия объясняет и силу, и слабость метода. Когда график круто пересекает ось (производная по модулю велика), касательная «прицеливается» в корень очень точно. Но если функция полога ($f'$ близка к нулю) или имеет точку перегиба между приближением и корнем, касательная может «улететь» далеко в сторону - приближение ухудшится. Именно поэтому метод касательных так чувствителен к выбору стартовой точки.

Условие сходимости

Метод касательных не сходится безусловно. Достаточные условия даёт теорема о сходимости метода Ньютона: если на отрезке, содержащем корень $x^{*}$, функция $f$ дважды непрерывно дифференцируема, производная не обращается в ноль ($f'(x) \ne 0$), а вторая производная $f''$ ограничена, то при начальном приближении $x_0$, достаточно близком к корню, итерации сходятся к $x^{*}$.

Удобный практический критерий - условие Канторовича и условие знака: если на отрезке $[a, b]$ выполнено $f(a)\cdot f(b) < 0$, а $f'$ и $f''$ сохраняют знак (функция монотонна и не меняет выпуклость), то метод Ньютона сходится из того конца отрезка, где знак функции совпадает со знаком второй производной: $f(x_0)\cdot f''(x_0) > 0$. Это правило выбора стартовой точки гарантирует, что касательные будут приближаться к корню монотонно, не перепрыгивая через него.

Метод касательных тесно связан с методом простой итерации: формула Ньютона - это итерация $x_{n+1} = \varphi(x_n)$ с функцией $\varphi(x) = x - \frac{f(x)}{f'(x)}$. В корне $\varphi'(x^{*}) = 0$, и именно обращение производной итерационной функции в ноль даёт не линейную, а квадратичную скорость.

Скорость сходимости

Главное достоинство метода Ньютона - квадратичная сходимость вблизи простого корня. Если обозначить погрешность $e_n = x_n - x^{*}$, то выполняется оценка

$e_{n+1} \approx \frac{f''(x^{*})}{2\,f'(x^{*})}\, e_n^{2}.$

Погрешность на следующем шаге пропорциональна квадрату предыдущей. На практике это означает удвоение числа верных десятичных знаков за итерацию: $10^{-2} \to 10^{-4} \to 10^{-8} \to 10^{-16}$. Поэтому для машинной точности обычно достаточно четырёх-шести шагов - несравнимо быстрее линейной сходимости метода бисекции, которому на тот же результат нужно около пятидесяти итераций.

Квадратичная скорость справедлива лишь для простого корня, где $f'(x^{*}) \ne 0$. Для корня кратности $m \ge 2$ производная в корне тоже обращается в ноль, и метод Ньютона деградирует до линейной сходимости. Скорость восстанавливают модифицированной формулой $x_{n+1} = x_n - m\,\frac{f(x_n)}{f'(x_n)}$, где множитель $m$ - кратность корня.

Критерий остановки

Итерации метода касательных прекращают по одному из условий или их комбинации:

• По приращению: $|x_{n+1} - x_n| < \varepsilon$ - соседние приближения почти совпали. Благодаря квадратичной сходимости это приращение надёжно оценивает погрешность $|x_n - x^{*}|$.
• По значению функции: $|f(x_{n+1})| < \delta$ - приближение близко к нулю по вертикали.
• По числу итераций: жёсткий лимит $n_{\max}$ как страховка от расходимости или зацикливания.

На практике основным берут критерий по приращению $|x_{n+1} - x_n|$, дополняя его лимитом числа шагов. Если за $n_{\max}$ итераций сходимости нет, это сигнал о неудачном старте или о нарушении условий сходимости - нужно сменить начальное приближение или метод.

Метод секущих как разностный аналог

Когда производная $f'(x)$ недоступна аналитически или дорога в вычислении, её заменяют разностным отношением по двум последним точкам:

$x_{n+1} = x_n - f(x_n)\,\frac{x_n - x_{n-1}}{f(x_n) - f(x_{n-1})}.$

Это метод секущих - вместо касательной через точку проводят секущую через две соседние. Производная не нужна, но скорость падает до порядка $\approx 1{,}618$ (золотое сечение) - между линейной бисекцией и квадратичным Ньютоном. Метод секущих удобен, когда $f$ задана таблично или её производную сложно выписать.

Частые ошибки

• Деление на ноль или малую производную. Если $f'(x_n) \approx 0$, касательная почти горизонтальна и шаг $\frac{f(x_n)}{f'(x_n)}$ выбрасывает приближение далеко. Перед делением проверяйте, что производная не мала.
• Неудачное начальное приближение. Далёкий или неправильно выбранный $x_0$ ведёт к расходимости, зацикливанию (например, между двумя точками) или сходимости к чужому корню. Используйте правило знака $f(x_0)\cdot f''(x_0) > 0$.
• Кратный корень без модификации. Для корня кратности $m \ge 2$ обычный метод Ньютона теряет квадратичную скорость и сходится линейно - нужна формула с множителем $m$.
• Игнорирование точек перегиба. Если между $x_0$ и корнем меняется знак $f''$, касательная может перепрыгнуть корень; монотонность сходимости не гарантирована.
• Остановка только по $|f(x_n)|$. Для пологой функции малое значение $f$ не означает близость к корню - надёжнее контролировать приращение $|x_{n+1} - x_n|$.

FAQ

Чем метод касательных отличается от метода бисекции? Бисекция использует только знаки функции, сходится линейно, но безусловно - достаточно непрерывности и смены знака. Метод Ньютона требует производной и хорошего старта, зато сходится квадратично: число верных знаков удваивается за шаг. Часто их комбинируют - бисекция локализует, касательные уточняют.

Что делать, если метод Ньютона расходится? Чаще всего виноват неудачный $x_0$ или малая производная по пути. Сначала грубо локализуйте корень (например, бисекцией), выберите старт по правилу $f(x_0)\cdot f''(x_0) > 0$ и проверьте, что $f'$ не обращается в ноль на отрезке. Если производная недоступна - перейдите к методу секущих.

Нужна ли точная производная или можно численную? Можно и численную: тогда метод касательных превращается в метод секущих с порядком сходимости около $1{,}618$ - медленнее квадратичного, но без аналитической производной. Точная производная даёт максимальную скорость, поэтому её используют, когда выписать $f'(x)$ несложно.

Коротко

Метод касательных Ньютона находит корень уравнения $f(x) = 0$ по итерационной формуле $x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)}$: в точке приближения строится касательная, и за новое приближение берётся её пересечение с осью. Вблизи простого корня метод сходится квадратично - погрешность убывает как $e_{n+1} \approx C\,e_n^{2}$, число верных знаков удваивается за шаг, и обычно хватает четырёх-шести итераций. Платой за скорость служат требования: гладкость функции, ненулевая производная и достаточно близкое начальное приближение, иначе итерации расходятся. Когда производная недоступна, её заменяет разностное отношение - это метод секущих с порядком сходимости около $1{,}618$.