Основные задачи вычислительной математики. Теория погрешностей. Приближённое вычисление значений функций заданных аналитически. Оценка погрешности вычислений.
Работа современного инженера, физика и любого другого исследователя связана с моделированием сложных процессов, происходящих в разных областях знаний и деятельности человека. Зачастую, моделирование является средним звеном в разработке проекта и его внедрения в производство. Процесс проектирования можно представить схематически: (рис 1).
рис 1.
Для исследования свойств построенной математической модели, в большинстве случаев, не удаётся аналитически решить задачу. Поэтому, вступают в силу методы вычислительной математики, которые позволяют решение каждой задачи довести до числового результата и оценить точность производимых вычислений.
При работе с приближёнными величинами приходится решать следующие задачи:
а) давать математические характеристики точности приближённых величин;
б) оценивать точность результата, когда известна точность исходных данных;
в) находить точность исходных данных, обеспечивающую заданную точность результата;
г) согласовывать точность исходных данных с тем, чтобы не затрачивать излишней работы при отыскании или вычислении одних данных, если другие данные слишком грубы;
а) Определение: абсолютная погрешность - это абсолютная величина разности между точным значением величины и её приближённым значением :
(1.1)
Здесь следует различать два случая:
- точное значение числа нам известно, что на практике очень редко, тогда пользуемся формулой (1.1).
Пример 1: а=5.129 а*=5.128, тогда ;
- точное значение числа неизвестно, тогда вводят понятие предельной абсолютной погрешности.
Определение: предельной абсолютной погрешностью приближённого числа называют всякое число, не меньшее абсолютной погрешности этого числа.
Таким образом, если - предельная абсолютная погрешность приближённого числа , то
(1.2)
отсюда следует, что
(1.3)
Значение предельной абсолютной погрешности, обычно, подбирается интуитивно по смыслу задачи.
Пример 2: Определить предельную абсолютную погрешность числа , заменяющего число , точное значение которого нам неизвестно.
Так как мы знаем, что , то можем утверждать:
(1.4)
и, следовательно, , т.е. можем сказать, что
(1.5)
Понятия абсолютной погрешности и предельной абсолютной погрешности, хотя и дают представление о точности вычислений, однако не всегда достаточны.
Например: если при измерении длины стержней получены результаты: <l1 и l2>, то, несмотря на совпадение предельных абсолютных погрешностей, качество первого измерения выше второго, т.к. если погрешность близка по величине от самого приближённого числа, то точность этого измерения недостаточна. Изданного примера понятно, что для оценки качества измерения, нам нужна абсолютная погрешность, приходящаяся на единицу длины. Такая погрешность носит название относительной погрешности.
Определение: относительной погрешностью приближённого числа называется отношение абсолютной погрешности этого числа к модулю соответствующего точного числа :
(1.6)
Поскольку точное значение величины нам часто не известно, то рассмотрим понятие предельной относительной погрешности .
Определение: предельной относительной погрешностью данного приближённого числа называется всякое число, не меньшее относительной погрешности этого числа:
(1.7)
Отсюда следует, что
(1.8)
т.е.
(1.9)
но, как известно:
(1.10)
Сопоставление формул (1.9) и (1.10) даёт соотношение между предельной абсолютной погрешностью и предельной относительной погрешностью :
(1.11)
Из этой формулы иногда выражают и пишут:
(1.12)
Рассмотрим примеры:
Пример 3: Вес 1 дм3 воды при равен г. Определить предельную относительную погрешность результата взвешивания.
Решение: очевидно, что предельная абсолютная погрешность г. и , следовательно:
(1.13)
Пример 4: При определении газовой постоянной для воздуха, получили . Зная, что относительная погрешность этого значения , найти пределы, в которых заключается R.
Решение: имеем: , тогда , т.е.
(1.14)
Теперь займёмся изучением распространения погрешностей из-за арифметических действий.
б) Рассмотрим функцию , пусть значения переменных , вычислены приближённо, где соответствующие абсолютные погрешности.
Нас интересует абсолютная и относительная погрешности вычисленных значений функции .
По определению видно, что абсолютная погрешность функции имеет вид:
обычно , поэтому, раскладывая в ряд Тейлора, можно ограничиться лишь линейными членами по . Получаем:
(1.15)
Отсюда получаем оценку:
(1.16)
Тогда для предельных абсолютных погрешностей имеем:
(1.17)
Разделив обе части (1.17) на , получаем предельную относительную погрешность при вычислении функции , в точке :
(1.18)
Или записывая более компактно:
(1.19)
Эту формулу можно переписать в виде:
(1.20)
в) Рассмотрим частные случаи:
1. Пусть . Изучим абсолютные и относительные погрешности суммы.
Решение: т.к.
(1.21)
то из (1.17) получаем
(1.22)
Также из (1.18) получаем:
(1.23)
2. Пусть, . Изучим абсолютные и относительные погрешности разности
Решение: ; , поэтому из (1.17) имеем
(1.24)
А из (1.18) получаем:
(1.25)
Ясно, что если и близкие друг к другу числа, то очень малое число, т.е. абсолютная погрешность разности будет очень большим числом. Поэтому при вычислениях, где это возможно, нужно избегать вычитания близких друг к другу чисел.
Например, если нам нужно вести вычисления по формуле: - объём между двумя сферами, где - очень малое число. Здесь лучше избавиться от вычитания и пользоваться аналогичной формулой , тем самым, обходя вычитание близких чисел, которое может быть больше относительной погрешности вычислений.
... (5.16) Непосредственное использование оценок погрешности (5.4), (5.8) и (5.12) неудобно, так как при этом требуется вычисление производных функции f(x). В вычислительной практике используются другие оценки. Вычтем из равенства (5.15) равенство (5.16): Ih/2 – Ih » Chk(2k – 1). (5.17) Учитывая приближенное равенство (5.16), получим следующее приближенное ...
... . Рассмотрение метода ветвей и границ для решения задачи о коммивояжере удобнее всего проводить на фоне конкретного примера. Пользуясь введенными здесь обозначениями, мы проводим это описание в следующей лекции. Введем некоторые термины. Пусть имеется некоторая чис- ловая матрица. Привести строку этой матрицы означает выде-лить в строке минимальный элемент (его называют константой приведения) ...
... на языке Turbo Pascal 7.0 для решении систем линейных алгебраических уравнений, используя метод простой итерации. 1.2 Математическая формулировка задачи Пусть А – невырожденная матрица и нужно решить систему где диагональные элементы матрицы А ненулевые. 1.3 Обзор существующих численных методов решения задачи Метод Гаусса В методе Гаусса матрица СЛАУ с помощью равносильных ...
... Вычисли: + + 1 – ( + ) 9. Сложите пять чисел, первое из которых равно 5, а каждое следующее на больше предыдущего. Самостоятельная деятельность учащихся на уроке В настоящее время в помощь учителю математики выпущено много разнообразных таблиц, материалов для организации самостоятельной работы. Всё это призвано обеспечить самостоятельную работу каждого ученика, лучше ...
0 комментариев