Ошибки округления при использовании метода трапеций

Граф вычислительного процесса для метода трапеций изображен на рис. 3. Предполагается, что все члены суммы, которые необходимо умножить на 2, сначала складываются, а их сумма умножается на 2. При этом не только экономится машинное время, но и уменьшается общая ошибка округления.

ошибки округления при использовании метода трапеций - student2.ru Пусть δi (i = 0, 1,..., п) — относительная ошибка каждой величины yi. Пусть аi (i = 0, 1, ...,n) — относительная ошибка округления при выполнении каждого из n + 1 сложений. На рис. 3 сложениям присвоены порядковые номера, поставленные около кружков со знаком «плюс». Например, а2 есть относительная ошибка округления при выполнении сложения.

Наконец, пусть μi — относительная ошибка округления при выполнении умножения (y4 + у2 + ... + yn-1) на 2, и пусть μ2 — относительная ошибка округления при выполнении умножения (у0 + 2y1 + ... + 2yn-1 + уn) на h/2. Тогда относительная ошибка округления Ih принимает вид

ошибки округления при использовании метода трапеций - student2.ru

ошибки округления при использовании метода трапеций - student2.ru

ошибки округления при использовании метода трапеций - student2.ru

Рис. 3. Граф вычислительного процесса для интегрирования

по правилу трапеций.

Абсолютная ошибка равна

ошибки округления при использовании метода трапеций - student2.ru

Теперь предположим, что ошибки округления при использовании метода трапеций - student2.ru равно среднему из всех значений yi и ошибки округления при использовании метода трапеций - student2.ru , так что можно пренебречь членами, содержащими произведения ошибки округления при использовании метода трапеций - student2.ru . Кроме того, предположим, что

ошибки округления при использовании метода трапеций - student2.ru

где ε - относительная ошибка округления при выполнении арифметической операции (мы уже записывали ее в виде 5•10-t) а φ — некоторая постоянная, определяющая соотношение между ошибками исходной информации и ошибка­ми округления. Тогда

ошибки округления при использовании метода трапеций - student2.ru

Но так как ошибки округления при использовании метода трапеций - student2.ru , то формула для ошибки приобретает следующий вид:

ошибки округления при использовании метода трапеций - student2.ru

Но ведь ошибки округления при использовании метода трапеций - student2.ru , поэтому

ошибки округления при использовании метода трапеций - student2.ru

Очевидно, при малом h первый член из заключенных в скобках значительно превосходит по своей величине оба остальных члена, поэтому в качестве верхней границы возможной ошибки eR можно принять выражение

ошибки округления при использовании метода трапеций - student2.ru (15)

Итак, мы получили интересный результат: верхняя граница для возможной ошибки округления возрастает как 1/h и при некотором h превосходит ошибку ограничения, пропорциональную h2. В действительности сама ошибка округления возрастает медленнее, а именно пропорционально h-p, где 0 < р < 1, но все же при достаточно малом h становится больше ошибки ограничения.

Это очередной поучительный пример расхождения между теорией и практикой вычислений. Теоретически можно сделать Ih сколь угодно близким к I, взяв h достаточно малым. Однако на практике ошибки округления ставят предел достижимой точности результата интегрирования.

В качестве примера рассмотрим интеграл

ошибки округления при использовании метода трапеций - student2.ru

На рис. 4 изображена зависимость общей ошибки EI от числа интервалов n.

ошибки округления при использовании метода трапеций - student2.ru

Рис. 4. Зависимость суммарной ошибки (ограничения и округле­ния) от количества интервалов при интегрировании sin x от 0 до π по правилу трапеций и по формуле Симпсона.

Зависимость эта определена для правила трапеций и правила Симпсона, которое будет рас­смотрено в разд. 6. Заметьте, что с возрастанием nобщая ошибка, равная сумме ошибок ограничения и округления, уменьшается до п = 775. После этого ошибка округления доминирует: увеличение п приводит к росту общей ошибки.

Наконец, заметим, что ошибки округления при использовании метода трапеций - student2.ru не дают вклада в величину члена, зависящего от 1/h в формуле для eI.

Вся ошибка, выражаемая этим членом, возникает во время суммирования

ошибки округления при использовании метода трапеций - student2.ru ,

так что нет смысла повышать точность вычисления самих величин yi

Эту ошибку можно существенно уменьшить, производя сложение с удвоенной точностью, в то время как остальные операции будут выполняться с обычной точностью. Такой метод называется вычислениями с частич­ной удвоенной точностью.

Наши рекомендации