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

При численном решении математических и прикладных задач почти неизбежно появление на том или ином этапе их решения погрешностей следующих трех типов [1].

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

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

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

Все три описанных типа погрешностей в сумме дают полную по­грешность результата решения задачи. Поскольку первый тип погрешно­стей не находится в пределах компетенции вычислителя, для него он слу­жит лишь ориентиром точности, с которой следует рассчитывать матема­тическую модель. Нет смысла решать задачу существенно точнее, чем это диктуется неопределенностью исходных данных. Таким образом, погреш­ность метода подчиняют погрешности задачи. Наконец, при выводе оценок погрешностей численных методов обычно исходят из предположения, что все операции над числами выполняются точно. Это означает, что погреш­ность округлений не должна существенно отражаться на результатах реа­лизации методов, т.е. должна подчиняться погрешности метода. Влияние погрешностей округлений не следует упускать из виду ни на стадии отбора и алгоритмизации численных методов, ни при выборе вычислительных и программных средств, ни при выполнении отдельных действий и вычисле­нии значений функций [1].

Пусть А и а - два "близких" числа; условимся считать А точным, а — приближенным.

Величина называется абсолютной погрешностью при­ближенного числа а, а его относительной погрешностью.

Числа Δ_a и δ_а такие, что Δ_a > Δа и δ_а = a, называются оценками или границами абсолютной и относительной погрешностей соответствен­но (к Δ_a и δ_а часто применяют также термин " предельные погрешно­сти "). Так как обычно истинные погрешности не известны, то там, где не может возникнуть недоразумений, будем иногда называть Δ_a и δ_а просто абсолютной и относительной погрешностями [1].

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

, а точная абсолютная погрешность результата есть

- модуль полного приращения функции. Главной, т.е. линейной частью этого приращения является, как известно, полный дифференциал du. Та­ким образом, имеем:

, (2.1)

т.е. за границу абсолютной погрешности результата приближенно может быть принята величина

(2.2)

Отсюда легко получается формула приближенной оценки относительной погрешности значения u:

(2.3)

Как частные случаи формул (2.2), (2.3) (точных для функций, линей­ных относительно x_i, или ln x_i соответственно) можно получить известные правила оценивания погрешностей результатов арифметических действий [1].

Действительно, пусть . Тогда и , т.е. при сложении и вычитании приближённых чисел их предельные погрешности складываются.

Пусть теперь , где можно считать все сомножители положительными. Так как и , то, согласно (2.3),

. (2.4)

Если же , то и, значит

. (2.5)

Последнее вместе с (2.4) означает известный результат о сложении пре­дельных относительных погрешностей при умножении и делении прибли­женных чисел.

Возвращаясь к сложению, рассмотрим относительную погрешность суммы n положительных приближенных чисел , имеющих гра­ницы относительных погрешностей соответственно:

,

где . Полученное неравенство говорит о том, что относительная погрешность суммы n положительных приближённых чисел не превосходит максимальной относительной погрешности слагаемых [1].

С вычитанием приближённых чисел дело обстоит хуже: оценка

относительной погрешности разности x₁ - х₂ двух приближенных поло­жительных чисел указывает на возможность сильного возрастания по­грешности при x₁ - х₂ → 0. В этом случае говорят о потере точности при вычитании близких чисел.

2.2.2. Статистический и технический подходы к учёту погрешности действий.

Рассмотренный выше аналитический (или классический) способ учета погрешностей действий, предполагающий точное оценивание по­грешностей, основанное либо на приведенных в предыдущем параграфе правилах подсчета погрешностей арифметических действий, либо на па­раллельной работе с верхними и нижними границами исходных данных, имеет два существенных недостатка. Во-первых, этот способ чрезвычайно громоздок и не может быть рекомендован при массовых вычислениях. Во-вторых, он учитывает крайние, наихудшие случаи взаимодействия погреш­ностей, которые допустимы, но маловероятны. Ясно, что, например, при суммировании нескольких приближенных чисел (полученных в результате измерений, округлений или каким-либо другим путем) среди них почти на­верное будут слагаемые как с избытком, так и с недостатком, т.е. произойдёт частичная компенсация погрешностей. При больших количествах од­нотипных вычислений вступают в силу уже вероятностные или стати­стические законы формирования погрешностей результатов действий. Например, методами теории вероятностей показывается, что математиче­ское ожидание абсолютной погрешности суммы n слагаемых с одинако­вым уровнем абсолютных погрешностей, при достаточно большом n, про­порционально √n. В частности, если n >10 и все слагаемые округлены до m -го десятичного разряда, то для подсчета абсолютной по­грешности суммы 5 применяют правило Чеботарева [1]

(2.6)

Различие в результатах классического и статистического подходов к оцениванию погрешности суммы рассмотрим на примере оценки по­грешности среднего арифметического нескольких приближенных чисел.

Пусть - среднее арифметическое n (>10) приближенных чисел (например, результатов измерений), имеющих одинаковый уровень абсолютных погрешностей . Тогда классическая оценка абсолютной погрешности величины х есть

,

т.е. такая же, как и у исходных данных. В тоже время по формуле (2.6) имеем

.

Как видим, применение правила Чеботарева приводит к естественному вы­воду о том, что арифметическое усреднение результатов измерений или наблюдений увеличивает точность, чего нельзя сказать на основе класси­ческой теории погрешностей [1].