- •«Вычислительная математика»
- •Часть 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.
Задачи для самостоятельного решения
1.
Вычислить
,
если
Оценить
погрешность результата.
Ответ:
2.
Даны два приближенных числа:
=
25,1;
=
17,43. Вычислить
.
Оценить погрешности результатов.
Ответы:
=
42,5;
=
0,055;
=
0,001;
=
7,7;
=
0,055;
=
0,007;
=
437;
= 0,002;
=
0,9;
=
1,44;
=
0,002;
=
0,003.
3.
Вычислить
.
Найти число верных знаков в ответе.
Ответ:
.
Нет верных знаков.
4.
Вычислить
с точностью
.
Ответ:
=
6,61.
1.3 Погрешности вычисления значения функции
Если
значение аргумента функции
-приближенное
число
,
то погрешность значения функции оценивают
по формуле:
,
где максимум производной - это ее
наибольшее (по модулю) значение на
промежутке
.
Погрешность
значения функции 2-х аргументов
можно оценить формуле:
,
где максимум частных производных
находят среди всех значений в области
D = (x,y)/ x-x x* x+x ; y-y y* y+y.
Аналогично
определяется погрешность значения
функции n
переменных
:
,
где максимум
частных производных находят среди всех
значений в области
.
Задачи
Объем шара:
=
2,03
0,01.
Найти
.
Решение.
;
. Ответ:
=
35,04
0,52.
.
Найти
.
Решение.
т.е.
=
10,1
0,0101;
=
4,6
0,0046.
Погрешность
функции:
.
. Ответ:
.
Задачи для самостоятельного решения
1.
В прямоугольном треугольнике катет
,
гипотенуза
.
Найти синус угла, лежащего напротив
и оценить погрешность. Ответ:
.
2.
Период колебаний простого маятника
.
Доказать, что
3.
Вычислить удельный вес металла
и оценить погрешность результата, если
металлический цилиндр диаметром 2 см
и высотой 11 см весит 93,4 г. Относительная
погрешность измерения длин равна 0,01, а
относительная погрешность взвешивания
0,001. Ответ:
=2,70
0,08.
Литература:
[1], гл. I, работы №1, №2.
[3], гл X, §1 №3108 – 3124.
[4], т. I, гл. I, §1 – 4.
[5], гл. I, §1.1 – 1.10, упражнения к §1.
[6], гл. I, §1 – 5.
§2 Методы решения нелинейных уравнений
В этом
параграфе рассматриваются методы
решения нелинейных уравнений. В общем
случае требуется найти решения уравнения
(или,
иначе говоря, корень функции
),
где
-
нелинейная функция одной переменной.
Все рассматриваемые в параграфе методы
ориентированы на отыскание одного из
возможных корней функции. Поэтому, если
имеет
более одного корня, то для отыскания
каждого из них потребуется отдельное
применение метода.
Использование
каждого из методов вычислений предполагает,
что известен отрезок
,
внутри которого содержится хотя бы один
корень функции. На практике определение
такого отрезка часто требует отдельной
процедуры и носит название этапа
отделения или изоляции корня. Этот этап
рассмотрен в первой части параграфа.
Полученный на этапе отделения корня отрезок затем последовательно уменьшается и значение корня функции в результате определяется с требуемой точностью. Этот процесс носит название этапа уточнения корня. Именно он реализуется с помощью различных методов, называемых методами решения нелинейных уравнений и рассматриваемых в последующих частях параграфа.
Для простоты в ходе дальнейшего изложения предполагается непрерывность функции , что обеспечивает необходимые свойства большинства рассматриваемых в параграфе методов. Если требуются более жесткие условия на функцию, то они оговариваются особо.
