- •Е.А. Михайлова, в.А. Михайлова Лабораторный практикум по курсу «вычислительная математика»
- •Предисловие
- •Лабораторная работа № 1 нахождение корня нелинейного уравнения
- •1. Постановка задачи
- •2. Методы решения задачи
- •2.1. Метод деления отpезка пополам
- •2.2. Метод простой итерации.
- •2.3. Метод Ньютона
- •Задание
- •Варианты заданий
- •2. Методы решения системы нелинейных уравнений
- •2.1.Метод простой итерации
- •2.2. Метод Ньютона
- •2. Формулы Ньютона – Котесса
- •2.1. Формулы метода прямоугольников
- •2.2. Формула трапеций
- •2.3. Формула Симпсона (формула парабол)
- •Несобственные интегралы 2-го рода
- •3. Пример
- •2. Методы решения
- •Задание
- •Варианты заданий
- •2. Решить уравнения:
- •Лабораторная работа № 6 краевая задача. Метод стрельбы
- •1. Постановка задачи
- •2. Метод стрельбы
- •Задание
- •Варианты заданий
- •Лабораторная работа № 7 метод наименьших квадратов
- •1. Постановка задачи
- •2. Метод наименьших квадратов
- •3. Пример
- •Задание
- •Варианты заданий
- •Лабораторная работа № 8 системы линейных уравнений
- •1. Постановка задачи
- •Метод прогонки
- •Задание
- •Список рекомендуемой литературы
- •Содержание
- •Лабораторный практикум по курсу «вычислительная математика»
- •4 00062, Г. Волгоград, просп. Университетский, 100.
2.2. Метод Ньютона
В переводной литературе можно встретить название метод Ньютона-Рафсона. Этот метод обладает гораздо более быстрой сходимостью, чем метод простой итерации.
Пусть известно некоторое приближение к корню , так что
|
|
|
Тогда исходную систему (2) можно записать следующим образом:
Разлагая уравнение (7) в ряд Тейлора в окрестности точки и ограничиваясь линейными членами по отклонению , получим:
|
, |
|
или в координатной форме:
|
|
(8) |
Систему (8) можно переписать в виде:
|
|
(9) |
Полученная система (9) является системой линейных алгебраических уравнений относительно приращений
|
. |
|
Значение функций F1, F2, …, Fn и их производные в (9) вычисляются при
|
. |
|
Определителем системы (9) является якобиан J:
|
|
(10) |
Для существования единственного решения системы уравнений (9) он должен быть отличен от нуля. Решив систему (9), например, методом Гаусса, найдём новое приближение:
|
. |
|
Проверяем условие (6). Если оно не удовлетворяется, находим и якобиан (10) с новым приближением и опять решаем (9), таким образом, находим 2-е приближение и т.д.
|
|
|
Итерации прекращаются, как только выполнится условие (6).
Задание
Используя метод Ньютона, найдите решения системы нелинейных уравнений с заданной точностью . Исследуйте сходимость итерационного процесса.
Варианты заданий
1 2
3 4
5 6
7 8
9 10
11 12
13 14.
15. 16.
17. 18.
19. 20.
21. 22.
ЛАБОРАТОРНАЯ РАБОТА № 3
ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ
1.Постановка задачи
Пусть необходимо вычислить интеграл
|
|
(1) |
|
Рис.4. Геометрический смысл определенного интеграла (1). |
По определению
|
|
(2) |
где |
|
(3) |
|
|
|
|
|
|
Многие методы приближенного интегрирования основаны на замене интегрируемой функции на всем отрезке или на его частях на более простую функцию , близкую к , легко интегрируемую точно и принимающую в узлах те же значения, что и . Эти формулы полчили название формулы Ньютона-Котесса. Рассмотрим случай равноотстоящих N узлов с шагом
|
|
(4) |