Метод касательных или метод Ньютона

корень функция касательный ньютон Метод касательных или метод Ньютонаэто итерационный численный метод нахождения корня (нуля) заданной функции. Метод был впервые предложен Исааком Ньютоном (1643−1727). Поиск решения осуществляется путем построения последовательных приближений и основан на принципах простой итерации. Метод обладает квадратичной сходимостью. Улучшением метода является метод хорд и касательных. Также метод Ньютона может быть использован для решения задач оптимизации, в которых требуется определить нуль 1й производной либо градиента в случае многомерного пространства.

Описание метода

Чтобы решить уравнение f (x)=0 методом простой итерации, его необходимо привести к следующей форме:

X=ц (x),.

где цсжимающее отображение.

Для наилучшей сходимости метода в точке очередного приближения X*должно выполняться условие ц`(x)=0. Решение данного уравнения ищут в виде ц (x)=x+б (x) f (x), тогда:

В предположении, что точка приближения «достаточно близка» к корню x?, и что заданная функция непрерывна.

При некоторых условиях эта функция в окрестности корня осуществляет сжимающее отображение, и алгоритм нахождения численного р-я уравнения сводится к итерационной процедуре вычисления.

Геометрический смысл метода Ньютона Состоит в том, что на каждом шаге мы строим касательную к графику y=f (x) в точке очередного последовательного приближения Xi, а за следующее приближение xi+1 берет точку пересечения этой касательной с осью OX. Тем самым наклон прямой подстраивается на каждом шаге наилучшим образом.

Затем, что по-другому метод Ньютона мы можем описать так на каждом шаге вместо исходного уравнения f (x)=0 мы решаем приближенное, линеаризованное в точке xi, уравнение.

.

в котором левая часть-это многочлен Тейлера первого порядка для функции f (x) в точке xi., т. е., линейная функция.

Решением линеаризованного уравнения служит следующее приближение xi+1, в то время как решаем исходного точного уравнения f (x)=0 служит искомый корень x*.

Решение нелинейных уравнений. Метод касательных (Ньютона).

Задание: 1) Отделить корни уравнения графически и программно.

2) Уточнить корни уравнения методом касательных с точностью е=0,0001. 3) Нарисовать схему применения метода к каждому корню уравнения.

Интерполяционный полином Ньютона Во многих случаях удобнее пользоваться записью интерполяционного многочлена в другой формев форме Ньютона.

Общее понятие. Разности и их свойства.

Где — обобщенные на область действительных чисел бимолекулярные коэффициенты.

Прямая интерполяционная формула Ньютона или первая интерполяционная формула Ньютона, применяется интерполирования.

Где, , а выражение вида конечный результат.

Обратная интерполяционная формула Где.

Интерполяционные формулы Ньютона-формулы вычислительной математики, применяющие для полиномиальной интерполяции.

Если узлы интерполяции равноотстоящие и упорядочены по величине так, что =const, т. е. то интерполяционный многочлен можно записать в формуле Ньютона.

Интерполирование функции. Полиномы Ньютона.

Задание:

• 1) Найти приближенное значение функции при заданном значении аргумента с помощью соответствующего интерполяционного полинома Ньютона, если функция задана в равноотстоящих узлах;
• 2)Оценить погрешность полученного значения.

Решение:

Численное интегрирование. Общие положения Очень часто в процессе решения конкретных задач, как в научной, так и в инженерной практике возникает необходимость вычисления определенных интегралов вида:

Подынтегральная функция f (x) может быть задана одним из трех способов:

• 1)Задается явная формула для f (x), например, f (x)=cos (x)2
• 2)Функция f (x) явно не задана, но ее значение может быть вычислено при любом x из отрезка a, b. Обычно это значение вычисляется по некоторой подпрограмме.
• 3)Для некоторого фиксированного конечного набора точек xi из отрезка a, b задается таблица значений xi, f (xi) .

Интегралы от функций первого типа иногда удается вычислить аналитически, либо вручную, либо с помощью машинных символьных систем. Интегралы от функций второго и третьего типа (а также первого, если не используются символьные методы) обычно находят численными методами, т. е. методами, позволяющими найти численное значение определенного интеграла приближенно с любой степенью точности.

Все методы приближенного вычисления определенных интегралов основаны на геометрическом смысле интеграла Ньютона-Лейбница. Он заключается в том, что определенный интеграл численно равен площади S криволинейной трапеции, ограниченной графиком функции f (x), осью абсцисс и двумя прямыми x =а и x=b. (аABb), как показано на рисунке 1.

Рис. 1.

Такие методы называют квадратурными формулами.

Процедура численного интегрирования заключается в том, что отрезок [а, b] разбивается на п частичных отрезков, а затем подынтегральная функция f (x) аппроксимируется некоторой другой функцией Ф (x), интеграл от которой вычисляется сравнительно просто. Для аппроксимации f (x) может быть использован любой класс простых функций, таких как полиномы, кусочные полиномы, тригонометрические, экспоненциальные или логарифмические функции. Конкретный выбор класса аппроксимирующих функций может зависеть от некоторых определенных свойств подынтегральной функции, но в наиболее распространенном случае, который здесь и рассматривается, в качестве таких функций используются полиномы.

Заменяя подынтегральную функцию на каждом шаге отрезками линий нулевого, первого и второго порядков, получаем соответственно приближенные формулы для вычисления интеграла:

• — метод прямоугольников;
• — метод трапеций;
• — метод Симпсона.

Метод прямоугольников Простейшим полиномом является константа. В формуле прямоугольников функция f (x) аппроксимируется своим значением в точке a (или в точке b), т. е.

Если значение функции f (x) берется в точке а, то формула (1) носит название формулы левых прямоугольников.

Для подсчета интеграла разделим интервал интегрирования a, b на n равных отрезков длины h = (b-a)/n. На каждом из отрезков функция f (x) заменяется прямоугольником с отрезками как основаниями, равными h и вертикальными боковыми сторонами высотой f (xi). При этом точка Xi выбирается, как середина каждого элементарного отрезка. Метод «средних» прямоугольников (метод средних) является более точным, чем методы «левых» и «правых» прямоугольников, когда в качестве точек xi могут выбираться левые или правые границы элементарных отрезков.

С геометрической точки зрения означает, что площадь криволинейной трапеции (aABb), ограниченной графиком функции f (x), осью абсцисс и двумя прямыми x=a и x=b, принимается приближенно равной площади ступенчатой фигуры, образованной из n прямоугольников с основаниями h = (b-a)/n и высотами f (xi) где хi = а +i *h, i = 1,2… n.

Для интервала a, b и шага интегрирования h полная формула будет записана в виде:

Где n — число разбиений для интервала [a, b], и точка x0 совпадает с a.

II Метод трапеции Следующим простейшим полиномом является линейная функция. Если она выбрана, совпадающей с f (х) в концах отрезка a и b, то получаем трапецию.

Площадь этой трапеции (интеграл от линейной функции), используемая в качестве приближения к значению интеграла от f (х), определяется по формуле:

Эта формула известна как формула трапеции.

Для того чтобы найти приближенное значение площади S, разделим отрезок интегрирования [a, b] на n равных частей длины h = (b-a)/n. В точках разбиения х0 =a, x1=a + h,…xn=b проводим ординаты у0, у1,…уn до пересечения с кривой у = f (х), т. е. уi = f (xi), хi =a + ih, i= 0,1,2…n. Концы ординат соединяем прямолинейными отрезкам, т. е. на каждом отрезке разбиения дугу графика подынтегральной функции у = f (x) заменяем стягивающей ее хордой (линейная интерполяция), и получим трапецию.

Тогда площадь криволинейной трапеции aABb приближенно можно считать равной площади фигуры, ограниченной ломаной линией aACD… Bb. Площадь этой фигуры, которую мы обозначим как S, равна сумме площадей трапеций:

+=.

Таким образом, для интервала a, b и шага интегрирования h полная формула приближенного значения интеграла будет записана в виде:

где n — число разбиений для интервала a, b и точка х0 совпадает с a, а точка хn совпадает с b.