- •Методические указания
- •"Численные методы"
- •Лабораторная работа №1. Тема: “Табулирование функции в заданном диапазоне аргумента”.
- •Краткое теоретическое введение. Общие правила вычислений на эвм
- •Вычисление математических функций, отсутствующих в языке программирования.
- •Задание на работу:
- •Контрольные вопросы.
- •Лабораторная работа № 2
- •Краткое теоретическое введение.
- •Задание на работу.
- •Контрольные вопросы
- •Лабораторная работа n 3. Тема: ”Решение алгебраических и трансцендентных уравнений”.
- •Краткое теоретическое введение Алгоритм уточнения корня методом половинного деления
- •Алгоритм уточнения корня методом хорд
- •Алгоритм уточнения корня методом касательных
- •Метод последовательных приближений (итераций)
- •Задание на работу
- •Контрольные вопросы.
- •Лабораторная работа №4 Тема: ”Решение систем алгебраических линейных уравнений”.
- •Краткое теоретическое введение Алгоритм метода Гаусса
- •Правило Крамера
- •Метод итераций и метод Зейделя
- •Задание к работе
- •Контрольные вопросы
- •Лабораторная работа № 5 Тема: «Вычисление определенного интеграла численными методами».
- •Краткое теоретическое введение Формула прямоугольников
- •Формула трапеций
- •Формула Симпсона
- •Задание на работу:
- •Контрольные вопросы
- •Лабораторная работа № 6
- •Краткое теоретическое введение
- •Задание на работу.
- •Контрольные вопросы
- •Лабораторная работа № 7 Тема: “Решение задач интерполяции и экстраполяции”.
- •Краткое теоретическое введение.
- •Задание к работе
- •Контрольные вопросы
- •Лабораторная работа № 8. Тема: “Статистическая обработка результатов эксперимента”.
- •Задание к работе:
- •Контрольные вопросы
- •Лабораторная работа № 9 Тема: “Внутренняя сортировка”.
- •Краткое теоретическое введение Сортировка простыми включениями
- •Сортировка простым выбором
- •Сортировка простым обменом
- •Сортировка включениями с убывающим приращением
- •Сортировка с разделением
- •Контрольные вопросы
- •Лабораторная работа № 10 Тема: “Внешняя сортировка“.
- •Контрольные вопросы.
- •Лабораторная работа № 11. Тема:“Поиск и выбор“.
- •Краткое теоретическое введение Алгоритмы последовательного поиска
- •Поиск в упорядоченной таблице
- •Бинарный поиск
- •Фибоначчиев поиск
- •Контрольные вопросы.
- •Список рекомендованной литературы
Алгоритм уточнения корня методом хорд
Сущность метода. Пусть отрезок изоляции [a; b] корня х уравнения F(х) = 0 найден, причем для определенности пусть F(a) < 0, F(b) > 0 и F'(х) > 0. График функции y = F(х) проходит через точки A(a; F(a)) и B(b; F(b)) (рис. 1). Составим уравнение хорды АВ как прямой, проходящей через точки А и В:
y = kx + l;
F(a) = ka + l;
F(b) = kb + l;
k
=
;
l = F(a) - a
y=
y - F(a) = (x - a).
Рис. 1 – Графическое представление метода хорд
далее находим абсциссу x1 точки пересечения хорды АВ с осью Ох, уравнение которой y = 0:
= x1
- a
x1
= a -
.
Число x1 примем за первое приближение корня х*. Далее, применяя этот же прием к отрезку изоляции [x1; b], на концах которого функция F(x) принимает противоположные знаки, получим второе приближение корня x2:
x2
= x1
-
Этот процесс можно продолжать неограниченно. Описанный процесс называется методом хорд. В результате получим последовательность вложенных отрезков[a; b] [x1; b] [x2; b] ... [xn; b] ...с неподвижным концом b. Последовательные приближения xn (n = 1, 2, ...) к точному значению корня х* вычисляются по формуле
(1)
называемой формулой метода хорд, и образуют монотонно возрастающую последовательность a = x0 < x1 < x2 < ... < xn < xn+1 < ... < b, ограниченную сверху числом b. По теореме Вейерштрасса эта последовательность имеет предел
xn
=
*.
Поскольку
F(х) непрерывна на [а; b], то
F(xn)
= F(
*).
Переходя теперь к пределу в равенстве (1), получаем
*
=
*
-
откуда, так как b -
,
следует, что F(
)
= 0. Но в связи с тем, что уравнение F(x) = 0
на отрезке [а; b]
имеет
единственный корень х*,
то х*
= х*.
Поскольку полученная последовательность (х„) сходится к корню уравнения х*, то любой ее член можно рассматривать в качестве приближенного значения корня. Практически последовательные приближения вычисляют до тех пор, пока не получат приближенное значение корня с требуемой точностью.
Алгоритм уточнения корня методом касательных
Сущность метода. Пусть [а; b] — отрезок изоляции корня х* уравнения F(х) = 0. И пусть для определенности F(a)<0, F(b)>0, F‘(x)>0 и F”(x)>0, x[а; b], то есть производные сохраняют постоянный знак (рис.2). Идея метода касательных, предложенного Ньютоном, сводится к замене небольшой дуги кривой у = F(х) касательной к кривой, проведенной в некоторой точке интервала [а; b]. Выберем, например, x0 = b, так как F(x0) F(x0)>0, и в точке В(x0, F(х0)) проведем касательную к кривой у = F(х). Ее уравнение: y - F(x0) = F(x0)(x - x0).'
Найдем теперь точку пересечения x1 касательной с осью Ох (у = 0):
0 - F(x0)
= F(x0)(x1
- x0)
x1
= x0
-
Рис. 2 – Графическое представление метода касательных
Эту точку x1 принимаем за первое приближение искомого корня х*. Через точку С(x1; F(x1)) снова проведем касательную y - F(x1) = F‘(x1)(x - x1), абсциссу точки пересечения которой с осью Ох примем за второе приближение х2 корня х*.
Имеем: 0
- F(x1)
= F‘(x1)(x2
- x1)
x2
= x1
-
Продолжая этот процесс далее, получим рекуррентную формулу,
называемую формулой метода касательных.
Заметим, что если в рассматриваемом случае (F’(x) > О, F”(x) > О, F(b) > 0, F(а) < 0) касательную провести в точке А, то есть положить x0 = а, то точка пересечения ее с осью абсцисс может оказаться вне отрезка изоляции корня [a; b]. Это значит, что метод касательных неприменим, если в качестве начальной точки x0 выбрать такую, в которой F(x0) F”(x0) < 0.
Как и в методе хорд, можно доказать (предлагаем читателю сделать это самостоятельно), что полученная числовая последовательность
x0>x1>x2>...>xn>xn+1>...>a
сходится к корню уравнения х*.
Для оценки погрешности приближения xn можно воспользоваться, как и в методе хорд, формулой
x* - xn xn - xn-1 .
Анализируя возможные
расположения кривой у = F(х)
на отрезке
изоляции, где последовательные приближения
по методу касательных обозначены
(i
= 0,1,2...), получаем правило для использования
метода касательных:
в качестве начального приближения x0
выбирается тот конец отрезка изоляции
(x0
= а или x0
=
b), в
котором выполняется условие
F(x0) F(x0) > 0
