- •Понятие о численных методах.
- •Погрешности вычислений.
- •3.Задача теории погрешности.
- •3.1. Абсолютная и относительная погрешность суммы и разности приближенных чисел.
- •Абсолютная и относительная погрешность произведения и частного приближенных чисел.
- •Абсолютная и относительная погрешность при вычислении функции нескольких переменных.
- •Контрольные вопросы.
- •Постановка задачи решения слу.
- •Метод простых итераций решения слу.
- •Метод Зейделя решения слу.
- •Методика решения задач.
- •Контрольные вопросы.
- •Постановка задачи решения уравнений. Отделение корней уравнения.
- •Метод половинного деления.
- •Методика решения задачи.
- •Метод хорд решения алгебраических и трансцендентных уравнений.
- •Методика решения задачи.
- •Метод касательных (метод Ньютона) решения алгебраических и трансцендентных уравнений.
- •Методика решения задач методом Ньютона.
- •Основные понятия аналитического приближения табличных функций.
- •Задача аналитического приближения табличных функций.
- •Интерполирование табличных функций.
- •Оценка погрешности полиномиальной интерполяции.
- •Интерполяционный многочлен Лагранжа.
- •5.Интерполяционный многочлен Ньютона.
- •Конечные разности.
- •Первый интерполяционный многочлен Ньютона.
- •Второй интерполяционный многочлен Ньютона.
- •Линейное интерполирование.
- •Формула линейного интерполирования.
Метод половинного деления.
Пусть
дано уравнение
и отделен
корень х*,
т.е. найден
отрезок
,
что
и на концах
отрезка
функция
имеет
разные знаки,
т.е.
.
Определение: Отрезок называется начальным интервалом неопределенности, потому что известно, что корень ему принадлежит, но его местоположение с требуемой точностью не определено.
Уточнение
корня:
Построение
последовательности вложенных друг в
друга отрезков, каждый из которых
содержит корень уравнения. Для этого
находится середина отрезка текущего
интервала неопределенности:
и в качестве
следующего
интервала
неопределенности
из
двух
возможных
выбирается тот, на концах которого
функция
имеет
разные
знаки.
Процесс завершается, когда длина текущего интервала неопределенности становится меньше заданной величины ɛ, задающей точность нахождения корня. В качестве приближенного значения корня будет середина последнего интервала неопределенности.
Методика решения задачи.
Найти начальный интервал неопределенности
одним из
методов
отделения
корней.
Задать малое положительное число ɛ.
Положить k=0.Найти середину текущего интервала неопределенности: .
Если
,
то положить ak+1=ak,
bk+1=ck.
Если
,
то положить ak+1=ck,
bk+1=bk,
в результате
находится текущий интервал
Если
,
то процесс завершить и положить
.
приближенное значение корня можно
найти по формуле:
.
Если же
,
то положить k=k+1
и перейти к п.2.
Замечание №1: Метод имеет линейную, но обусловленную сходимость, и его погрешность за каждую итерацию уменьшается в два раза:
и тогда число итераций можно найти по формуле:
,
Замечание №2: Метод половинного деления не обобщается на системы нелинейных уравнений и не может использоваться для нахождения корней четной кратности.
Пример:
Найти корень уравнения
методом половинного деления
с точность
ɛ=0,01.
Решение.
Н
айдем
начальный интервал неопределенности
графическим методом. Для этого построим
график функции
.
Имеем, что наше уравнение имеет один корень, который, как видно из графика функции лежит в отрезке [0;0,5].
Положим k=0.
Т.о. имеем, а0=0, b0=0,5. Проверим знаки функции на концах отрезка:
функция
имеет разные знаки на концах отрезка.
Положим k=1
Положим k=2
Положим k=3
Положим k=4
Положим k=5
Положим k=6
процесс
завершим и положим
Метод хорд решения алгебраических и трансцендентных уравнений.
В методе хорд при
тех же предположениях, что и в методе
половинного деления отрезок
делится
в отношении
.
Г
еометрически
метод хорд эквивалентен
замене
кривой
хордой, проходящей через точки
и
.
Уравнение хорды АВ
имеет вид:
,
пологая, что х=х(1)
и у=0,
получим:
Предположим, что
сохраняет постоянный знак, и рассмотрим
два случая:
1 случай: |
2 случай: |
Неподвижен конец а:
|
Неподвижен конец b:
|
Замечание №1:
Если
задача сводится к рассматриваемому,
если уравнение записать в форме
.
Замечание
№2: Для
выявления неподвижного конца используется
условие:
,
где t=а
или t=b.
Если неподвижен конец а, применяется формула:
Если неподвижен конец b, применяется формула:
