- •Численные методы
- •2012 Введение
- •1. Решение нелинейных уравнений
- •Постановка задачи
- •1.2. Основные этапы отыскания решения
- •1.3. Метод половинного деления
- •1.4. Метод простой итерации
- •Приведение уравнения к видудля обеспечения выполнения неравенства
- •1.5. Метод Ньютона (метод касательных)
- •1.6. Видоизменённый метод Ньютона
- •1.7. Метод хорд
- •1.8. Комбинированный метод
- •Решение систем линейных алгебраических уравнений
- •2.1. Постановка задачи
- •2.2. Метод простой итерации
- •2.3. Метод Зейделя
- •3. Решение систем нелинейных уравнений
- •3.1. Постановка задачи
- •3.2. Метод Ньютона для системы нелинейных уравнений
- •3.3. Метод итерации для нелинейной системы уравнений
- •3.4. Метод скорейшего спуска решения нелинейных систем
- •3.5. Метод скорейшего спуска для случая линейной системы
- •4. Приближение функций
- •4. 1. Метод наименьших квадратов
- •4.2. Построение интерполяционных многочленов
- •Многочлен Лагранжа
- •Многочлен Ньютона с конечными разностями
- •5. Вычисление собственных значений матрицы Методом Данилевского
- •6. Вычисление определённых интегралов. Метод симпсона (метод парабол)
- •7. Численное решение дифференциальных уравнений
- •7.1. Постановка задачи Коши
- •7.2. Метод Эйлера
- •7.3. Модифицированные методы Эйлера
- •7.4. Метод Рунге – Кутта
- •Решение краевой задачи для линейного дифференциального уравнения второго порядка методом прогонки
1.6. Видоизменённый метод Ньютона
Если производная мало изменяется на отрезке, то в расчетной формуле метода можно положить:. Отсюда для корняуравненияполучаем последовательные приближения
.
Геометрически этот способ означает, что касательные заменяются прямыми, параллельными касательной к кривой , в ее фиксированной точке. Этот способ избавляет от необходимости вычислять каждый раз значения производной, поэтому эта формула полезна, еслисложна.
1.7. Метод хорд
Рассмотрим еще одну модификацию метода Ньютона. Пусть известно, что простой корень уравнениянаходится на отрезке, то есть. И предположим, чтопри(если это не так, то будем рассматривать уравнение). Заменим кривуюхордой.
А
х2 х1 b = x0
а х* x
В
Рис. 9
y
B
а = х0 х1 х2 b
х* x
А
Рис. 10
Возможны два случая: 1) (рис. 9); 2)(рис. 10). В первом случае конец неподвижен и последовательные приближения:
(9)
образуют ограниченную монотонно убывающую последовательность, причем .
Во втором случае неподвижен конец , а последовательные приближения:
(10)
образуют ограниченную монотонно убывающую последовательность, причем Итак, в результате получаем следующее
Выбор начального условия:
1. Рассматриваем только случай (иначе).
2. Начальное приближение x0 выбираем из условия
Неподвижен тот конец, для которого знак функции совпадает со знаком ее второй производной.
Критерий окончания. Критерий окончания итераций метода хорд такой же, как и для метода Ньютона. При заданной точности вычисления нужно вести до тех пор, пока не будет выполнено неравенство.
Пример. Найти положительный корень уравнения с точностью . Отделим корень. Так как,, то. Разделим интервал пополам:, тогда.
Найдём производные: ,. Исходя из того, что, тои пользуемся формулой (10):,.
, ,.
Так как , то.
1.8. Комбинированный метод
Пусть , аисохраняют постоянные знаки на отрезке. Соединяя методы хорд и касательных, получаем метод на каждом этапе, которого находим значения по недостатку и значения по избытку точного корняуравнения. Пусть– последовательные приближения метода хорд,– последовательные приближения метода касательных. Пошаговая иллюстрация представлена на рис. 11.
Возможны 4 случая: 1) , 2) ,
3) , 4),
которые можно свести к первому случаю.
y
x0=а х1 х2
x* х
Рис. 11
.
. .
Очевидно, что и.
По окончании процесса за значение корня лучше всего взять среднее арифметическое полученных значений:.
Пример. Вычислить положительный корень уравнения . Так как, то.
, на, поэтому.
.
.
; .
Так как , то
; .
Так как , то.