- •16 Челабчи в.Н “Численные методые” Часть 1 Конспект лекций
- •Часть 1 Специальность 6.080400 «Информационные управляющие системы и технологии»
- •Введение в численные методы анализа и линейной алгебры. Http://www.Gasu.Ru/resour/eposobia/metody/posob.Html
- •Введение
- •Точность представления данных при решении прикладных задач
- •Некоторые часто используемые численные методы
- •Интерполяция
- •Линейная интерполяция на отрезке.
- •Квадратичная интерполяция на отрезке.
- •Полиномиальная интерполяция в области определения функции.
- •Аппроксимация
- •Аппроксимация табулированной функции полиномом.
- •Решение систем линейных алгебраических уравнений
- •Прямые методы решения слау.
- •Итерационные методы решения систем алгебраических уравнений.
- •Решение систем нелинейных уравнений
- •Решение алгебраических и трансцендентных уравнений
- •Метод половинного деления
- •М етод хорд.
- •Метод простых итераций.
- •Ковариационный и корреляционный анализ
- •Численное интегрирование и дифференцирование Постоянный шаг аргумента
- •Переменный шаг аргумента
- •Предварительная аппроксимация и сглаживание данных
М етод хорд.
Метод хорд также, как и метод половинного деления, предназначен для уточнения корня на интервале [a,b], на концах которого значения функции f(x) разного знака. Первое приближение берется в точке x1, где пересекается прямая линия, проведенная через точки f(a) и f(b) с осью x. В качестве нового интервала для продолжения итерационного процесса выбирается тот из двух [a,x1] или [x1,b], на концах которого функция f(x) имеет разные знаки. В нашем слчае это интервал [a,x1]. На следующих итерациях после уточнения границ нового интервала расчет проводится по формуле:
x=(f(b)*a - f(a)*b)/(f(b) - f(a));
Итерационный процесс уточнения корня заканчивается, когда выполнится условие |xi-xi-1|< или когда значение функции f(xi) попадет в область шума | f(xi) |<1.
Метод Ньютона (метод касательных).
Рассмотрим графическую иллюстрацию метода на Рис. 1.3.
Предположим, что известно начальное приближение x0 к корню. В точке x0 вычислим левую часть решаемого уравнения f(xQ), а также производную в этой точке f’’(х0). Следующее приближение к корню найдем в точке х1, где касательная к функции f(x), проведенная из точки (x0,f0), пресекает ось абсцисс. Для этого используем формулу (1.19) построенную на основе элементарных геометрических соображений.
(1.19)
Затем считаем точку х1в качестве начальной и продолжаем итерационный процесс (определяем следующее приближение x2) и т.д. Из Рис. 1.3 видно, что таким способом можно быстро приблизиться к корню х*. При этом с каждой итерацией расстояние между очередным xm и предыдущим хm-1 приближениями к корню будет уменьшаться. Расчет уточненного значения корня в общем случае проводим по уравнению (1.20). Процесс уточнения корня закончим, когда выполнится условие (1.21).
(1.20)
xm+1 - хm < ε, (1.21)
где ε – абсолютный итерационный допуск,
m – номер итерации (приближения).
Метод Ньютона обладает высокой скоростью сходимости, его недостатком является необходимость вычисления на каждой итерации не только левой части уравнения, но и ее производной.
Метод простых итераций.
В этом методе решается уравнение (1.22).
x=φ(x) (1.22)
Процесс решения заключается в следующем. Пусть имеется приближенное значение корня xm-1. Для получения уточненного значения используем итерационную формулу:
xm = φ(xm-1) (1.23)
Уточнения (итерации) выполняются до тех пор, пока не выполнится условие (1.19). Процесс уточнений может быть сходящимся, когда от итерации к итерации непрерывно уменьшается значение xm+1 - хm или может быть расходящийся когда значение xm+1-хm непрерывно растет. Условием сходимости итерационного процесс является
φ(x)< 1 (1.24)
Переход от уравнения (1.18) к уравнению (1.22) можно осуществить различными способами в зависимости от вида функции f(x). При таком переходе необходимо построить функцию φ(x) так, чтобы выполнялось условие сходимости (1.24). Чаще всего используется преобразование (1.25)
φ(x)= x+pf(x). (1.25)
Преобразование (1.25) удобно тем, что, изменяя значение параметра p можно управлять скоростью сходимости итерационного процесса.
