
- •Е.А. Михайлова, в.А. Михайлова Лабораторный практикум по курсу «вычислительная математика»
- •Предисловие
- •Лабораторная работа № 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) |