- •Основы вычислительной математики Учебно-методическое пособие
- •Оглавление
- •Предисловие
- •I. Введение
- •Контрольные вопросы
- •II. Системы линейных алгебраических уравнений
- •II.1. Основные понятия и определения
- •II.2. Прямые методы решения систем линейных уравнений
- •II.3. Метод Гаусса решения слау
- •II.4. Решение слау с помощью -разложения
- •II.5. Метод итераций решения слау
- •II.6. Метод Зейделя решения слау
- •II.7. Сходимость итерационных процессов для систем линейных уравнений
- •III. Методы решения нелинейных скалярных уравнений
- •III.1. Локализация корней
- •III.2. Метод простых итераций
- •III.3. Приведение нелинейного уравнения к виду , допускающему сходящиеся итерации
- •III.4. Метод Ньютона (Метод касательных)
- •III.5. Контрольные вопросы
- •III.6. Задачи для самостоятельного решения
- •IV. Приближенное решение систем нелинейных уравнений
- •IV.1. Метод Ньютона
- •IV.2. Метод итераций
- •IV.3. Контрольные вопросы
- •IV.4. Задачи для самостоятельного решения
- •V. Интерполяция и аппроксимация
- •V.1. Интерполяция
- •V.2. Интерполяционный многочлен в форме Лагранжа
- •V.3. Интерполяционный многочлен в форме Ньютона
- •V.4. Приближение методом наименьших квадратов
- •V.5. Контрольные вопросы
- •V.6. Задания для самостоятельного решения
- •VI. Равномерное приближение функций
- •V.1. Полиномы Чебышева
- •VI.2. Понятие о равномерном приближении функций
- •VI.3. Контрольные вопросы
- •VII. Приближенное дифференцирование
- •VII.1. Постановка вопроса
- •VII.2. Формулы приближенного дифференцирования, основанные на первой интерполяционной формуле Ньютона
- •VII.3. Формулы приближенного дифференцирования, основные на формуле Стирлинга
- •VII.4. Контрольные вопросы
- •VII.5. Задачи для самостоятельного решения
- •VIII. Численное интегрирование функций
- •VIII.1. Общие замечания
- •VIII.2.Квадратурные формулы Ньютона-Котеса
- •VIII.3. Формула трапеций и ее остаточный член
- •Общая формула трапеций (правило трапеций).
- •VIII.4. Формула Симпсона и ее остаточный член
- •Общая формула Симпсона (параболическая формула).
- •VIII.5. Формулы Ньютона-Котеса высших порядков
- •VIII.6. Контрольные вопросы
- •VIII.7. Задачи для самостоятельного решения
- •IX. Приближенное решение обыкновенных дифференциальных уравнений
- •IX.1. Интегрирование дифференциальных уравнений с помощью степенных рядов
- •IX.2. Метод последовательных приближений
- •IX.3. Метод Эйлера
- •IX.4. Метод Эйлера с итерациями
- •IX.5. Метод Рунге-Кутта
- •IX.6. Контрольные вопросы
- •IX.7. Задачи для самостоятельного решения
- •X. Преобразование фурье
- •X.1. Разложение периодических функций в ряд Фурье
- •X.2. Эффект Гиббса
- •X.3. Контрольные вопросы
- •IX.4. Задачи для самостоятельного решения
- •Список литературы
III.5. Контрольные вопросы
Определение абсолютной и относительной погрешности.
Теорема Больцано-Коши.
Достаточное условие существования единственного корня уравнения на .
Метод дихотомии (уточнение корней уравнение на ).
Метод простых итераций (итерационная формула).
Теорема о сходимости процесса итерации.
Приведение нелинейного уравнения к виду .
Метод Ньютона решения уравнения .
Достаточное условие сходимости метода Ньютона.
Рабочая формула модификации метода Ньютона.
III.6. Задачи для самостоятельного решения
1.
Доказать графическим и аналитическим
методами решения нелинейных уравнений,
что на отрезке
уравнение
имеет единственный корень; построить
рабочую формулу метода простых итераций
для решения нелинейного уравнения.
№ |
|
|
№ |
|
|
1. |
|
|
11 |
|
|
2. |
|
|
12 |
|
|
3. |
|
|
13 |
|
|
4. |
|
|
14 |
|
|
5. |
|
|
15 |
|
|
6. |
|
|
16 |
|
|
7 |
|
|
17. |
|
|
8 |
|
|
18. |
|
|
9 |
|
|
19. |
|
|
10 |
|
|
20. |
|
|
2. Найти по методу Ньютона наименьший
положительный корень уравнения
с точностью до 0,0001.
3.
Вычислить с точностью до 0,0005 единственный
положительный корень уравнения
4. Найти
наибольший положительный корень
уравнения
с точностью до
.
5. Найти
действительные корни уравнения
с точностью до трёх значащих цифр.
IV. Приближенное решение систем нелинейных уравнений
IV.1. Метод Ньютона
Рассмотрим нелинейную систему уравнений
(4.1.)
с действительными левыми частями.
Запишем (4.1) в векторном виде
,
(4.2)
где
-мерный вектор,
-мерный вектор-функция.
Для решения системы (4.1) воспользуемся методом последовательных приближений.
Предположим,
что найдено
-е
приближение
одного из изолированных корней
векторного уравнения (4.2). Тогда точный
корень (4.2) можно представить в виде
(4.3)
где
правка (погрешность корня). Подставляя (4.3) в (4.2), получим
(4.4)
Предполагая,
что функция
непрерывно дифференцируема в некоторой
области, содержащей
и
,
разложим левую часть (4.4) по степеням
малого вектора
,
ограничиваясь линейными членами:
(4.5)
где
– матрица Якоби системы функций
относительно переменных
,
т.е.
или в краткой записи
Формулу (4.5) можно записать в виде
Полагая,
что матрица
– не особенная, получим
Следовательно,
(4.6)
(Метод Ньютона).
За нулевое приближение можно принять грубое значение искомого корня.
Пример 4.1. Приближенно найти положительные решение системы уравнений
(4.7)
Решение.
Построим графики функций
и
.
Рис. 4.1
Кривые,
определяемые системой (4.7) пересекаются
приблизительно в точках
и
.
Исходя из начального приближения
,
вычислим первое приближение корней,
производя вычисления с точностью до
4-х десятичных знаков. Полагая
,
имеем
Составим матрицу Якоби
где
.
Отсюда
причём
т.е. матрица
– неособенная. Составим обратную ей
матрицу
где
.
Используя формулу (4.6), получаем первое приближение
Аналогично находят дальнейшие приближения. Результаты вычислений приведены в таблице.
|
|
|
|
|
0 |
3,4 |
0,0899 |
2,2 |
0,0633 |
1 |
3,4899 |
–0,0008 |
2,2633 |
–0,0012 |
2 |
3,4891 |
–0,0016 |
2,2621 |
–0,0005 |
3 |
3,4875 |
– |
2,2616 |
– |
Останавливаясь на приближении
,
имеем
;
,
причем
Пример 4.2. Методом Ньютона приближенно найти положительное решение системы уравнений
исходя из начального
приближения
.
Решение. Имеем:
Отсюда
Составим матрицу Якоби
Имеем
и
Находим обратную матрицу
По формуле (4.6) получаем первое приближение
Далее
вычислим второе приближение
имеем
и
Отсюда
Используя (4.6), получим
Аналогично находим дальнейшие приближения:
и т.д.
Ограничиваясь третьим приближением, получим:
.
