- •Методы вычислений
- •Санкт-Петербург
- •Введение
- •Лабораторная работа № 1 Теория приближенных вычислений
- •1.1. Абсолютная и относительная погрешности
- •1.2. Погрешность округленного числа
- •1.3. Погрешности арифметических действий
- •1.4. Погрешности элементарных функций
- •1.5. Способ границ
- •1.6. Обратная задача теории погрешностей
- •Вопросы по теме
- •1.8. Задания к лабораторной работе № 1
- •2. Лабораторная работа №2 Численные методы решения скалярных уравнений
- •2.1. Метод хорд
- •2.2. Метод касательных
- •2.3. Метод простой итерации
- •2.4. Вопросы по теме
- •2.5. Задания к лабораторной работе №2
- •3. Лабораторная работа №3 Численные методы решения систем нелинейных уравнений
- •3.1. Метод Ньютона
- •3.2. Вопросы по теме
- •3.3. Задание к лабораторной работе №3
- •4. Лабораторная работа№4 Численное интегрирование
- •4.1. Метод прямоугольников
- •4.2. Метод Симпсона
- •4.3. Метод трапеций
- •4.4. Метод Монте – Карло
- •4.5. Вопросы по теме
- •4.6. Задание к лабораторной работе №4
- •5. Лабораторная работа №5 Решение обыкновенных дифференциальных уравнений
- •5.1. Метод Пикара
- •5.2. Метод Эйлера и его модификации
- •5.3. Метод Рунге – Кутты
- •5.4. Метод Адамса
- •5.5. Вопросы по теме
- •5.6. Задание к лабораторной работе № 5 Задание 5.1.
- •Литература
- •Содержание
1.4. Погрешности элементарных функций
Пример 1.14.
Исходные числовые значения аргумента заданы цифрами, верными в строгом смысле. Найти абсолютную и относительную погрешности функции. Определить количество верных цифр в строгом смысле по относительной погрешности в следующих элементарных функциях:
а) cos(0,47);
б) у = е-3,1;
в) у =
;
г) у = ln(68,214).
Решение.
а) Находим значение величины х. Оно будет равно 0,891568.
Абсолютная погрешность аргумента е0,47 = 0,005.
Тогда абсолютная и относительная погрешности величины х равны:
еcos(0,47) = │sin(0.47)│ ∙ 0.005 = 0,00226443;
δcos(0,47) = │tan(0.47)│ ∙ 0.005 = 0,00253983 ≤ 0.005 = 0,5 ∙ 10-2.
Это означает, что в числе 0,891568 две цифры после запятой верны в строгом смысле.
Ответ: 0,892.
б) Находим значение величины у. Оно будет равно 0,0450492.
Абсолютная погрешность аргумента еу = 0,05. Тогда абсолютные и относительные погрешности величины у равны:
ех = │е-3,1│ ∙ 0,05 = 0,0225246;
δу = 0,05 = 0,5 ∙ 10-1.
Это означает, что в числе 0,0450492 одна цифра после запятой верна в строгом смысле.
Ответ: 0,04.
в) Находим значение величины у. Оно будет равно 4,6378875.
Абсолютная погрешность аргумента еу = 0,005. Тогда абсолютные и относительные погрешности величины у равны:
еу
=
= 1,078077 ∙ 10-3;
δу
=
= 2,3245002 ∙ 10-4
< 0.5 ∙ 10-3;
Это означает, что в числе 4,6378845 три цифры после запятой верны в строгом смысле.
Ответ: 4,6378.
г) Находим значение величины у. Оно будет равно 4,2226498.
Абсолютная погрешность аргумента еу = 0,0005. Тогда абсолютные и относительные погрешности величины у равны:
еу
=
= 7,3298736 ∙ 10-6;
δу
=
= 1,7358469 ∙ 10-6
< 0,5 ∙ 10-5.
Это означает, что в числе 4,2226498 пять цифр после запятой верны в строгом смысле.
Ответ: 4,222649.
Пример 1.15.
Вычислить значение величины с помощью метода строгого учета границ абсолютных погрешностей после каждой операции.
А =
,
если а = 12,34, b
= 14,3.
Решение.
При пооперационном строгом учете ошибок промежуточные результаты после округления до одной запасной (с учетом вычисленной параллельно величины погрешности) и их погрешности заносят в табл. 1.1.
Таблица 1.1. Расчетная таблица для вычисления
погрешности выражения А =
а b
+
ln(a)
b
+ ln(a)
A
1 2,34 14.3 3.513 3.78 7.30 2.5129 16.81 0.434
ea
eb
e
e
e
e
eb
+ ln(a) eA
0.005 0.05 0.00071 0.0066 0.0073 0.0004 0.05041 0.0017
Значения погрешностей для удобства округлим до двух значащих цифр по избытку и тоже занесем в таблицу.
Цифры даны верными в строгом смысле, значит, еа = 0,005, еb = 0,05.
Найдем
= 3,51283. Абсолютная погрешность равна
(воспользуемся табл. 1.1):
еа
=
= 0,0007117 ≈ 0,00071.
Из полученного значения погрешности видно, что в результате верны две значащие цифры после запятой, т. е.
= 3,51283 ≈ 3,513.
Это число внесем в таблицу.
Найдем абсолютную
погрешность
= 3,781534.
Она будет равна
еb
=
= 0,0066107 ≈ 0,0066.
Значит, в числе b будет одна верная цифра после запятой.
Аналогично находим значения всех остальных действий и функций:
е
= 0,00071 + 0,0066 = 0,00731 ≈ 0,0073;
еln(а)
=
= 0,000405 ≤ 0,0005;
еb + ln(a) = 0.05 + 0.000405 = 0.050405 ≤ 0.5;
еА
=
=
= 0,0017.
Округляя результат А до верной цифры, получаем окончательный ответ.
Ответ: А = 0,434 ± 0,002.
