
- •«Вычислительная математика»
- •Часть 1.
- •§ 1. Действия с приближенными величинами
- •§ 2. Методы решения нелинейных уравнений
- •§ 3. Нормы вектора и матрицы
- •§ 4. Прямые методы решения систем линейных алгебраических уравнений с квадратной матрицей
- •Введение
- •§1. Действия с приближенными величинами
- •1.1 Основные понятия
- •Задачи для самостоятельного решения
- •1.2 Погрешности арифметических операций
- •Задачи для самостоятельного решения
- •1.3 Погрешности вычисления значения функции
- •Задачи для самостоятельного решения
- •Литература:
- •§2 Методы решения нелинейных уравнений
- •2.1 Изоляция корней
- •Задачи для самостоятельного решения
- •2.2 Уточнение корней при помощи методов, не требующих вычисления производных
- •2.2.1 Метод деления отрезка пополам (метод бисекции)
- •2.2.2 Метод простой итерации
- •2.2.3 Метод хорд
- •2.3 Уточнение корней при помощи методов, использующих производные
- •2.3.1 Метод Ньютона – Рафсона (метод касательных)
- •2.3.2 Модификации метода Ньютона
- •2.3.3 Метод хорд и касательных
- •Задачи для самостоятельного решения
- •Литература:
- •§3 Нормы вектора и матрицы
- •Задачи для самостоятельного решения
- •Литература:
- •§4 Прямые методы решения систем линейных алгебраических уравнений с квадратной матрицей
- •4.1 Метод Гаусса
- •4.1.1 Схема единственного деления
- •4.1.2 Схемы с выбором главного элемента
- •4.2 Метод Жордана
- •4.2.1 Классическая схема Жордана
- •4.2.2 Схема оптимального исключения
- •2). Схема оптимального исключения. Записываем в таблицу только рабочую часть расширенной матрицы, т.К. Остальные строки не преобразовываются.
- •4.3 Метод Холецкого.
- •4.4 Прямые методы решения системы линейных уравнений в случае квадратной матрицы специального вида
- •4.4.1 Метод квадратных корней
- •4) Находим решение системы (1):
- •4.4.2 Метод прогонки
- •4.5.1 Использование матриц вращения для решения систем
- •Для решения системы линейных уравнений
- •Литература
- •§5 Итерационные методы решения систем линейных алгебраических уравнений
- •5.1 Метод Якоби
- •5.2 Метод Зейделя
- •5.3 Метод наискорейшего спуска
- •5.4 Метод сопряженных градиентов
- •Литература
- •§6 Вычисление определителя матрицы и построение обратной матрицы
- •6.1 Вычисление определителя
- •6.2 Построение обратной матрицы методом окаймления
- •6.3 Построение обратной матрицы методом пополнения
- •Литература
- •§7 Решение систем линейных уравнений с прямоугольными матрицами
- •7.1 Классификация систем
- •7.2 Нахождение обобщенного решения переопределенной системы с помощью первой трансформации Гаусса
- •Задача 1
- •Задача 2
- •Литература
- •§ 8 Решение проблемы собственных значений
- •8.1 Приведение матрицы к трехдиагональному или почти треугольному виду при помощи метода Гивенса
- •Задача 1
- •8.2 Метод Якоби и схема его реализации для получения собственных значений матрицы
- •1. Классическая схема Якоби. На каждом шаге выбирают максимальный по модулю недиагональный элемент матрицы : и строят матрицу для обнуления этого элемента.
- •6. Копченова н. В., Марон и. А. Вычислительная математика в примерах и задачах - м.: Наука, 1972.
2.3.2 Модификации метода Ньютона
Для решения уравнения в общем случае можно использовать следующие модификации метода Ньютона.
1). Метод Ньютона с регулировкой шага.
(7)
Здесь
- константа, регулирующая длину шага.
Если
,
то выполняется неполный ньютоновский
шаг. Процесс может быть стационарным,
если
не зависит от k.
Регулировка длины шага позволяет
обеспечить сходимость метода для
достаточно произвольно выбранного
значения
.
2). Конечно-разностный метод Ньютона – Рассела.
(8)
Здесь
- шаг конечно-разностной аппроксимации
производной. Эта модификация лишена
одного из недостатков метода Ньютона
– не нужно вычислять производную
.
2.3.3 Метод хорд и касательных
Преимущество
этого метода в том, что при тех же
предположениях относительно
и
,
что и в методе Ньютона, мы получаем
последовательные приближения
и
,
лежащие по разные стороны от
,
поэтому можно следить за достигнутой
точностью в процессе решения. Пусть
и
не меняют знак на
.
Если
,
то расчеты ведут по формулам:
(9)
Если
,
то
(10)
За
начальный промежуток
принимают
.
В результате расчетов получаем
последовательность вложенных промежутков
,
таких, что
.
В качестве критерия остановки расчетов
может быть использовано выполнение
неравенства:
.
Тогда приближенное решение считают
равным
.
Задачи
1.
Уточнить корень уравнения
(
)
методом деления отрезка пополам. Точность
.
Решение.
Так как процесс вычисления этим методом
сходится медленно, уменьшим промежуток
изоляции. Возьмем [-0,5;
-0,2]
[-1;
0]. Проверка:
,
следовательно,
.
Теперь уточняем корень.
Номер
шага
.
.
Затем
.находим
.
Вычисления производим с запасным
десятичным знаком. Результаты расчетов
заносим в таблицу. (Стрелки указывают
перенос значений
при
и
).
|
|
|
|
|
|
-0,5 |
-0,2 |
-0,35 |
-0,331 <0 |
|
-0,5 |
-0,35 |
-0,425 |
0,055 >0 |
|
-0,425 |
-0,35 |
-0,3875 |
-0,139 <0 |
3 |
-0,425 |
-0,3875 |
-0,406 |
-0,044 <0 |
4 |
-0,425 |
-0,406 |
-0,416 |
0,008 >0 |
5 |
-0,416 |
-0,406 |
-0,411 |
|
Процесс
закончен, т.к.
:
Ответ:
.
2.
Составить расчетные формулы для
определения методом простых итераций
корней уравнения
.
Найти один из корней с точностью
.
Корни уравнения лежат в промежутках:
(см. п. 2.1).
Решение.
Для уточнения
и
можно использовать функцию
,
полученную из уравнения. Здесь
;
.
Условие сходимости
выполнено и на промежутке [-1;
-0,5],
и на [1;
1,5], т.к.
,
если
[-1;
-0,5] и
,
если
[1;
1,5]. Значит,
для уточнения этих корней можно
использовать расчетную формулу
,
.
Для
уточнения
эта формула не годится, т.к.
терпит бесконечный разрыв в точке
=
-0,2
.
Можно из исходного уравнения выразить
по-другому:
.
Здесь
.
Условие
выполнено при
[-0,5;
0].
Уточним
корень
,
используя формулу
.
В качестве
можно взять любую точку промежутка
изоляции, например,
.
Дальнейшие расчеты приведены в таблице.
-
0
-0,25
1
-0,20098
2
-0,20033
Расчеты
закончены, т.к.
0,001.
Ответ:
.
3.
Найти один из корней уравнения
методом касательных и методом хорд с
точностью
0,01.
Решение.
Изоляция корня. Уравнение имеет 2 корня,
причем
.
Проверка (для второго корня):
.
Уточнение корня. а) Метод касательных.
Функция
непрерывная,
дважды дифференцируемая,
.
В качестве
нужно взять тот конец промежутка [1;
2], где
.
Т.к.
,
то
.
Результаты расчетов приведены в таблице.
k |
|
|
|
|
0 |
2 |
2 |
4 |
1,5 |
1 |
1,5 |
0,25 |
3 |
1,417 |
2 |
1,417 |
0,008 |
|
|
Вычисления
прекращаем, т.к.
.
Ответ:
.
б)
Метод хорд
.
Производные
сохраняют знак на промежутке
[1;2],
т.к.
и
.
Поскольку выполнено условие
,
считаем,
что
.
Далее
вычисляем по формуле:
.
Результаты расчетов заносим в таблицу.
-
k
0
1
-1
1
1,333
-0,222
2
1,400
-0,04
3
1,412
-0,007
Вычисления
прекращаем, т.к.
.
Ответ:
.