Округление чисел в ЭВМ

Если предположить, что исходная информация не содержит никаких ошибок и все вычислительные процессы выполняются абсолютно точно, то всегда существует третий тип ошибок — ошибки округления, которые возникают при переводе чисел из одной системы счисления в другую и последующем представлении их в разрядной сетке машины, а также при получении внутри машины чисел, разрядностью большей, чем это допустимо, например, при умножении. В этом случае число А округляют, т. е. заменяют его машинным числом [A] заданной разрядности. Округление (обозначим его знаком ¨) называется оптимальным, если для любого машинного числа [A] справедливо ¨А=[A]. Пусть [A]₁и [A]₂ - два последовательных машинных числа, тогда при оптимальном округлении вещественное число A такое, что [A]₁‹A‹ [A]₂ заменяется либо числом [A]₁, либо числом [A]₂. Если ¨А≤A, то говорят об округлении по недостатку, если ¨А ≥A, то говорят об округлении по избытку. Округление называют симметричным, когда ¨А= -¨(-А). Различают три вида симметричного округления.

1. Округление в направлении к нулю, когда вещественное число округляется до ближайшего к нулю машинного числа.

2. Округление в направлении от нуля, когда округление производится до машинного числа, лежащего дальше от нуля, чем вещественное число А.

3. Округление по дополнению, когда округление производится до ближайшего машинного числа.

В качестве параметров, по которым будут сравниваться способы округления, целесообразно использовать максимальную величину модуля погрешности, т. е. _max , где  = А-[A], и математическое ожидание погрешности округления .

Округление к нулю или усечение.

Для конкретности считаем, что числа в машине представлены в прямом коде с запятой, фиксированной перед старшим разрядом, [A]= а_-1... a _{- n} (правильная дробь). Пусть в результате каких-либо действий над машинными числами внутри машины сформировалось число [А]', имеющее k = n + t разрядов. Очевидно, что самый простой способ округления состоит в отбрасывании хвоста числа [А]', который состоит из лишних разрядов, т.е. разрядов с номерами а_{- n -1} , a _{- n -2} , ... , a _{- n - t} .

Если считать, что появление чисел с абсолютной величиной А, но разных знаков равновероятно и равновероятны все значения хвоста чисел одного знака, то математическое ожидание погрешности в данном случае равно нулю, т.е. = 0.

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

При действиях с числами одного знака погрешность усечения носит систематический характер, что приводит к накоплению погрешности. Это обстоятельство заставляет исследовать другие способы округления, которые рассматриваются пока для прямых кодов.

Округление от нуля.

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

Округление по недостатку.

Реализация данного способа базируется на анализе знака округляемого числа. Если [А]'> 0, то округление заключается в отбрасывании хвоста. Если же [А]' < 0, то хвост также отбрасывается, а к величине оставшейся части добавляется единица в младший разряд, если хвост не равен нулю. Таким образом, реализация данного способа еще более усложнена по сравнению со способом округления от нуля за счет анализа знака числа [А]', хотя величина _max осталась при этом прежней.

Если рассматривать округление чисел только одного знака, то при А' > 0 данный способ совпадает с усечением, а при А' < 0 - с округлением от нуля. Отсюда ясно, что он не может конкурировать с усечением результатов.

Округление по избытку.

Этот способ во всем подобен предыдущему, с тем отличием, что добавление единицы в младший разряд сохраняемой части числа производится, когда оно больше нуля и хвост не равен нулю. При А'<0 хвост просто отбрасывается. Характеристики данного способа точно такие же, как у предыдущего, за исключением знака величины , который меняется на противоположный.

Округление по дополнению.

Данный способ представляет собой объединение способов округления от нуля и к нулю. Его реализация связана с коррекцией сохраняемой части числа А', которая производится по результатам анализа значения старшей цифры отсекаемой части a_-n-1, т. е. цифры дополнительного ( n + 1)-го разряда (ДР). Когда (n + 1 разряд) = 0, происходит округление к нулю, в противном случае - от нуля.

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

По своим характеристикам способ округления по дополнению лучше, чем усечение. Систематические ошибки при округлении чисел одного знака обусловлены в данном случае тем, что округление особенно неточно производится, когда значение отсекаемой части близко к половине единицы младшего сохраняемого разряда. Этот недостаток устраняется в следующем способе округления.

Усовершенствованное округление по дополнению.

В этом случае решение о коррекции сохраняемой части числа А' принимается на основе анализа значения всех разрядов его отсекаемой части, а не только старшего. Таким образом, усовершенствованный способ дает хорошее округление в случае чисел одного знака, однако это достигается за счет усложнения его реализации.

Упрощенное округление по дополнению.

При реализации способов округления по дополнению из-за возникновения переносов при суммировании сохраняемой части числа с единицей округления необходимо выполнить операцию сложения, что требует дополнительного времени, т. е. снижает реальное быстродействие ЭВМ и, кроме того, может повлечь за собой переполнение разрядной сетки. Этот способ округления состоит в том, что младший разряд сохраняемой части числа принудительно устанавливается в единицу, если старший разряд отбрасываемой части равен единице.

При этом если в неокругленном результате разряд а_{- n} равен единице, то он не изменяется при округлении. Если в неокругленном результате операции значение разряда а_{- n} есть нуль, то установка его в единицу вносит погрешность.

Таким образом, при равновероятном появлении нуля и единицы в младшем разряде сохраняемой части для знакопеременных чисел снова получим симметричное распределение погрешностей.

Вероятностное округление.

Для такого округления необходимо иметь датчик случайных величин (0 или 1), единица с выхода которого прибавляется к младшему разряду сохраняемой части числа. Погрешность округления при равновероятном распределении значений отбрасываемой части является случайной величиной с нулевым математическим ожиданием.

Таким образом, самым простым способом округления является усечение, при котором не требуется дополнительных затрат времени и оборудования. Однако на практике важнее всего точность вычислений. Только для трех способов округления по дополнению максимальная ошибка близка к половине единицы младшего разряда машинного числа, т. е. является наименьшей. Наиболее быстродействующим из них является упрощенный способ, а наиболее точным — усовершенствованный. Поэтому предпочтение тому или другому способу округления следует отдавать только после анализа требований, предъявляемых к быстродействию и погрешности вычислений конкретной машины. Будем в дальнейшем пользоваться способом округления по дополнению.