
- •Численные методы
- •Введение
- •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. Метод Рунге – Кутта
- •Решение краевой задачи для линейного дифференциального уравнения второго порядка методом прогонки
- •Приложение
- •Библиографический список
- •Оглавление
- •3.2. Метод Ньютона для системы нелинейных уравнений 26
- •3.3. Метод итерации для нелинейной системы уравнений 30
1. Решение нелинейных уравнений
Постановка задачи
Пусть
дана некоторая функция
и требуется найти все или некоторые
значения
,
для которых
.
Значение
,
при котором
,
называетсякорнем
(или
решением)
уравнения. Относительно функции
часто предполагается, что
дважды непрерывно дифференцируема в
окрестности корня.
Корень
уравнения называетсяпростым,
если первая производная функции
в точке
не равна нулю, т. е.
.
Если же
,
то корень
называетсякратным
корнем.
Геометрически
корень уравнения есть точка пересечения
графика функции
с
осью абсцисс. На рис. 1 изображен график
функции
,
имеющей четыре корня: два простых
и два кратных
.
Рис. 1
Большинство методов решения уравнения ориентировано на отыскание простых корней.
1.2. Основные этапы отыскания решения
В процессе приближенного отыскания корней уравнения обычно выделяют два этапа: локализация (или отделение) корня и уточнение корня.
Локализация
корня заключается в определении отрезка
,
содержащего один и только один корень.
Не существует универсального алгоритма
локализации корня. Иногда удобно бывает
локализовать корень с помощью построения
графика или таблицы значений функции
.
На наличие корня на отрезке
указывает различие знаков функции на
концах отрезка. Основанием для этого
служит следующая теорема.
Теорема.
Если
функция
непрерывна на отрезке
и принимает на его концах значения
разных знаков так что
,
то отрезок
содержит по крайней мере один корень
уравнения.
Однако
корень четной кратности таким образом
локализовать нельзя, так как в окрестности
такого корня функция
имеет постоянный знак. На этапе уточнения
корня вычисляют приближенное значение
корня с заданной точностью
.
Приближенное значение корня уточняют
с помощью различных итерационных
методов. Суть этих методов состоит в
последовательном вычислении значений
,
которые являются приближениями к
корню
.
1.3. Метод половинного деления
Метод
половинного является самым простым и
надежным способом решения нелинейного
уравнения. Пусть из предварительного
анализа известно, что корень уравнения
находится на отрезке
,
т. е.
,
так, что
.
Пусть функция
непрерывна на отрезке
и принимает на концах отрезка значения
разных знаков, т.е.
.
Разделим
отрезок
пополам. Получим точку
.
Вычислим значение функции в этой точке:
.
Если
,
то
– искомый корень, и задача решена. Если
,
то
– число определённого знака:
либо
.
Тогда либо на концах отрезка
,
либо на концах отрезка
значения функции
имеют разные знаки. Обозначим такой
отрезок
.
Очевидно, что
и длина отрезка
в два раза меньше, чем длина отрезка
.
Поступим аналогично с отрезком
.
В результате получим либо корень
,
либо новый отрезок
и т. д. (рис. 2).
Рис. 2
Середина
-го
отрезка
.
Очевидно, что длина отрезка
будет равна
,
а так как
,
то
.
(1)
Критерий
окончания.
Из соотношения (1) следует, что при
заданной точности приближения
вычисления
заканчиваются, когда будет выполнено
неравенство
или
неравенство
.
Таким образом, количество итераций
можно определить заранее. За приближенное
значение корня берется величина
.
Пример.
Найдем
приближенно
с точностью
.
Эта задача эквивалентна решению уравнения
,
или нахождению нуля функции
.
В качестве начального отрезка
возьмем отрезок
.
На концах этого отрезка функция принимает
значения с разными знаками:
.
Найдем число
делений отрезка
,
необходимых для достижения требуемой
точности. Имеем:
.
Следовательно,
не позднее 6-го деления найдем
с требуемой точностью,
.
Результаты вычислений представлены в
таблице 1.
Таблица 1
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
|
1,0000 |
1,0000 |
1,0000 |
1,1250 |
1,1250 |
1,1406 |
1,1406 |
|
2,0000 |
1,5000 |
1,2500 |
1,2500 |
1,1875 |
1,1875 |
1,1562 |
|
1,5000 |
1,2500 |
1,1250 |
1,1875 |
1,1406 |
1,1562 |
1,1484 |
Зн
|
- |
- |
- |
- |
- |
- |
- |
Зн
|
+ |
+ |
+ |
+ |
+ |
+ |
+ |
|
5,5938 |
0,7585 |
-0,2959 |
0,1812 |
-0,0691 |
0,0532 |
-0,0078 |
|
1,0000 |
0,5000 |
0,2500 |
0,1250 |
0,0625 |
0,0312 |
0,0156 |