
- •Итерационные циклические программы
- •Курск 2000
- •1. Цель работы
- •2. Методические рекомендации по теоретической части
- •2.1. Вычисление корней уравнения с одним неизвестным
- •2.1.1. Метод дихотомии
- •2.1.2. Метод простых итераций
- •2.1.3. Метод касательных
- •2.1.4. Метод хорд
- •2.1.5. Комбинированный метод
- •2.2. Приближенное вычисление определенного интеграла
- •2.2.1. Формула прямоугольников
- •2.2.2. Формула трапеций
- •2.2.3. Формула Симпсона
- •2.3. Численные методы решения дифференциальных уравнений первого порядка
- •2.3.1. Метод Эйлера
- •2.3.2. Метод Рунге-Кутта
- •3. Пример решения задачи
- •4. Задачи, выносимые на практические занятия
- •5. Задачи для самостоятельного решения
- •6. Содержание отчета
2.1.3. Метод касательных
Пусть х0 - нулевое приближение корня. Разложим f(x) по степеням (x-x0) в ряд Тейлора:
.
Линеаризуем уравнение:
Тогда
или в общем виде
(2)
Метод касательных имеет следующий геометрический смысл (рис.2).
Рис.2. Метод касательных
Построим график f(x).
Рассмотрим на графике точку M0(x0,f(x0)) и проведем касательную к М0. Уравнение касательной
Найдем точку
х1. Тогда (так как y=0)
и отсюда
Метод касательных всегда приводит к цели, если х0 не лежит слишком далеко от . Возможно применение такого варианта метода:
где все наклонные прямые проводят параллельно касательной в точке М0.
Метод касательных - частный случай
метода простых итераций, где
Поэтому метод сходится, когда
т.е.
Можно доказать, что итерации сходятся
к корню с той стороны, с которой
(3)
Таким образом, алгоритм данного метода имеет следующий вид. Выделяется корень в [a,b]. Проверяется (3). Рассчитывается (2). Если |xn+1 - xn| < , то вычисления прекращаются.
2.1.4. Метод хорд
Если заменить производную в (2) разностью, получим итерационную формулу метода хорд:
(4)
Геометрический смысл метода хорд представлен на рис.3.
Метод хорд является двухшаговым, так как требует запоминания хn и xn-1. Можно доказать, что итерации сходятся к корню с той стороны, с которой
(5)
Таким образом алгоритм данного метода имеет следующий вид. Выделяется корень в [a,b]. Задается x0=a, x1=b и вычисляется x2 по (4). Проверяется x2 на (5). Если (5) выполняется, то берется новый интервал [x2,b], иначе [a,x2] и x1=a. Если |xn+1 - xn| < , то вычисления заканчиваются.
Рис.3. Метод хорд
2.1.5. Комбинированный метод
Если рассматриваемый участок графика не имеет изломов и точек перегиба, то метод хорд и метод касательных дают точки, расположенные по разные стороны от искомого корня.
Рис.6 поясняет комбинированный метод, где а1 строится по методу хорд, а b1 - по методу касательных.
Таким образом, получается двустороннее приближение к корню и новый интервал [a1,b1]. Сходимость метода очень высока. Критерий прекращения вычислений:|bn - an| < 2.
Рис.4. Комбинированный метод
2.2. Приближенное вычисление определенного интеграла
2.2.1. Формула прямоугольников
Пусть требуется вычислить
.
Разобьем отрезок [a,b] на
n частей. Если в качестве
исходных точек выбрать левые концы этих
частей (рис.5), то приближенное значение
Jn определенного
интеграла:
где yi = f(xi) и xi+1 - xi = x для любых i.
Рис.5. Метод прямоугольников
Если в качестве исходных точек выбрать правые концы частей разбиения, то
На практике
обычно принимают
Если соответствующую ординату
обозначить как
,
то
,
где Rn - остаточный член, определяющий погрешность и вычисляемый по формуле
,
где a
z
b и
- наибольшее значение второй производной
подынтегральной функции.
2.2.2. Формула трапеций
Разобьем отрезок [a,b] на n равных частей. Подсчитаем значение f(x) в точках xk: f(xk)=yk. Тогда по формуле площади трапеции (рис.6)
Рис.6. Метод трапеций
Суммируя, получаем
где
,
a
z
b и
имеет тот же смысл, что и в формуле
прямоугольников.