- •А.П. Бырдин н.В. Заварзин а.А. Сидоренко л.П. Цуканова основы численных методов Учебное пособие
- •А.П. Бырдин н.В. Заварзин а.А. Сидоренко л.П. Цуканова
- •Введение
- •1. Действия над приближенными числами
- •1.1. Основные источники погрешностей
- •1.2. Приближенные числа. Абсолютная и относительная погрешности
- •1.3. Правила записи приближенных чисел
- •Решение. В нашем случае и . Следовательно,
- •Решение. Имеем .
- •2. Интерполирование функции
- •2.1. Постановка задачи интерполирования
- •2.2. Вычисление значений многочлена по схеме Горнера
- •2.3. Интерполяционный многочлен Лагранжа
- •Интерполяционный многочлен Ньютона
- •Сопоставление интерполяционных формул Лагранжа и Ньютона. Погрешность интерполяции
- •2.6. Интерполирование функции кубическими сплайнами
- •Методы численного решения систем
- •3.1. Метод Гаусса
- •3.2. Метод итерации
- •3.3. Метод Зейделя
- •4. Методы численного решения
- •4.1. Отделение корней
- •Метод половинного деления
- •4.3. Метод хорд
- •4.4. Метод Ньютона
- •Комбинированный метод
- •4.6. Метод итерации
- •4.7. Метод Ньютона для системы двух уравнений
- •Метод итерации для системы двух уравнений
- •5. Численное дифференцирование
- •5.1. Постановка вопроса
- •5.2. Формулы приближенного дифференцирования, основанные на первой интерполяционной формуле Ньютона
- •5.3. Конечно-разностные аппроксимации производных
- •6. Среднеквадратичное приближение функций
- •6.1. Метод наименьших квадратов
- •7. Численное интегрирование
- •Численное решение обыкновенных дифференциальных уравнений
- •8.1. Понятие о численном решении задачи Коши
- •8.2. Метод Эйлера
- •8.3. Методы Рунге-Кутта
- •8.4. Численные решение систем дифференциальных уравнений первого порядка
- •9. Численное решение дифференциальных уравнений в частных производных
- •Решение задачи Дирихле для уравнения Лапласа методом сеток
- •Решение смешанной задачи для уравнения гиперболического типа методом сеток
- •9.3. Решение смешанной задачи для уравнения теплопроводности типа методом
- •Библиографический список
- •Оглавление
- •1. Действия над приближенными числами ..……………...4
- •2. Интерполирование функций ……………………….…... 9
- •3. Методы численного решения систем линейных
- •4. Методы численного решения нелинейных уравнений
- •10. Библиографический список …….…….……..……… 110
- •Составители: Бырдин Аркадий Петрович
- •394026 Воронеж, Московский просп., 14
3.3. Метод Зейделя
Метод Зейделя представляет собой модификацию метода итераций. Основная его идея заключается в том, что при вычислении -го приближения неизвестной учитываются уже вычисленные ранее -е приближения неизвестных .
Запишем систему (3.6) в виде
Выберем произвольно начальные приближения корней . Тогда первое приближение по методу Зейделя вычисляется по формулам
Так же вычисляются следующие приближения
Обычно метод Зейделя дает лучшую сходимость, чем метод простой итерации. Процесс Зейделя может сходиться даже в том случае, если расходиться процесс итерации. Однако, возможны случаи, когда процесс Зейделя сходится медленнее процесса итерации или, когда процесс итерации сходится, а процесс Зейделя расходится.
Пример. Методом Зейделя решить систему уравнений
Решение. Приведем эту систему к виду, удобному для итерации,
В качестве нулевых приближений корней возьмем
Применяя процесс Зейделя, последовательно получим
Результаты вычислений с точностью до четырех знаков помещены в таблице 4.
Таблица 4
|
|
|
|
0 1 2 3 4 5 |
1.2000 1.2000 0.9992 0.9996 1.0000 1.0000 |
0.0000 1.0600 1.0054 1.0001 1.0000 1.0000 |
0.0000 0.9480 0.9991 1.0001 1.0000 1.0000 |
Точные значения корней: .
4. Методы численного решения
НЕЛИНЕЙНЫХ УРАВНЕНИЙ
4.1. Отделение корней
При решении многих алгебраических и трансцендентных уравнений точное значение их корней определить бывает достаточно сложно. Поэтому важное значение приобретают способы приближенного нахождения корней и оценки степени их точности.
Пусть дано уравнение
, (4.1)
где - непрерывная функция от . Всякое значение , обращающее функцию в нуль, т.е. такое, что , называется корнем уравнения (4.1) или нулем функции .
Приближенное нахождение действительных корней уравнения обычно складывается из двух этапов:
отделение корней, т.е. установление промежутков, в которых содержится только один корень уравнения;
уточнение приближенных корней, т.е. доведение их до заданной степени точности.
Известно, что если функция непрерывна и принимает на концах отрезка значения разных знаков, т.е. , то внутри этого промежутка имеется хотя бы один корень уравнения. Отделение корней уравнения для непрерывной в области определения функции можно осуществить различными способами.
Составляют таблицу значений функции на определенном промежутке изменения аргумента , и если окажется, что для соседних значений аргументов значения функции имеют разные знаки, то корень находится между ними.
Уравнение заменяют равносильным . Строят графики функций и ; искомый корень является абсциссой точки пересечения этих графиков.
Строят график функции на промежутке изменения ; тогда абсцисса точки пересечения графика с осью - корень уравнения, т.е. .
Пример. Выяснить, сколько корней имеет уравнение , и найти промежутки, в которых находятся эти корни.
Решение. Рассмотрим три функции:
.
Уравнение эквивалентно уравнению . Отелим его корни двумя способами (таблица 5).
Таблица 5
-
-3.0
-2.0
-1.0
0.0
1.0
-14.05
-4.14
1.63
3.00
-0.72
-14.00
-4.00
2.00
4.00
2.00
0.05
0.14
0.37
1.00
2.72
Из таблицы значений функции на промежутке с шагом изменения , равным 1, видно, что существуют корни на отрезках и , так как значения функции на концах отрезка имеют разные знаки.
Графики функций и пересекаются в двух точках, абсциссы которых являются решениями уравнения , заключенными в указанных промежутках (рис.1).
Рис.1
После отделения корней производится итерационное уточнения каждого корня одним из существующих методов. Рассмотрим простейшие методы уточнения корней.