6. Численный эксперимент

В табл. 2.1 представлены результаты расчетов сумм (2.2) при =1.1 дляn=2^k. Поскольку сумма этого ряда может быть вычислена независимо через интегралы [6]

,

этот пример взят в качестве тестового. (При =1.1 сумма ряда (2.1)z_точн10.5844484649508098).

Для ряда (2.1) справедливо разложение

, (2.25)

где B_k– числа Бернулли.

Отметим, что коэффициенты c_iприне убывают приi, что, тем не менее, не мешает проводить весьма точные экстраполяции поnпри фиксированныхi.

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

Как видно из таблицы 2.1, погрешность z_nz_точнвычисленных значений весьма велика (~3) даже для достаточно большихn. Оценка погрешностиz_nz_n^*по формуле (2.6) при всех приведенных в таблицеnвполне приемлема, а дляn>500 совпадает с самой погрешностью с точностью не ниже 3-х десятичных знаков. Погрешностьz_n^*z_точнуточненных значений при увеличенииnбыстро уменьшается.

Отношение погрешностей , характеризующее размытость оценок, не превышает 0.3 и уменьшается в два раза при удвоенииn. Погрешностьz_n^**z_точндважды экстраполированных значений существенно меньше, чем погрешности вычисленных чисел и результатов первой экстраполяции.

В результате первая экстраполяция при N10⁵уменьшает погрешность результатов до 10^-5, а вторая - до 10^-10. Нетрудно подсчитать, что для получения такой точности путем прямого вычисления сумм потребуется более 10⁵⁰- 10¹⁰⁰слагаемых. На современном компьютере это заняло бы более 10³⁰лет.

Табл. 2.1

n	znzточн	znzn*	zn*zточн		zn**zточн
2	-9.11	-6.49	-2.61100	2.8710-1	-	-
4	-8.60	-7.19	-1.40100	1.6310-1	-3.5010-1	2.4810-1
8	-8.07	-7.36	-7.0610-1	8.7510-2	-9.2510-2	1.3110-1
16	-7.55	-7.21	-3.4110-1	4.5210-2	-2.2410-2	6.5710-2
32	-7.06	-6.89	-1.6210-1	2.2910-2	-5.2910-3	3.2610-2
64	-6.59	-6.51	-7.6310-2	1.1510-2	-1.2310-3	1.6210-2
128	-6.15	-6.11	-3.5710-2	5.8110-3	-2.8910-4	8.0810-3
256	-5.74	-5.72	-1.6710-2	2.9110-3	-6.7410-5	4.0310-3
512	-5.35	-5.35	-7.8010-3	1.4510-3	-1.5710-5	2.0110-3
1024	-4.99	-4.99	-3.6410-3	7.2810-4	-3.6710-6	1.0010-3
2048	-4.66	-4.66	-1.7010-3	3.6410-4	-8.5610-7	5.0310-4
4096	-4.35	-4.35	-7.9310-4	1.8210-4	-1.9910-7	2.5110-4
8192	-4.06	-4.06	-3.7010-4	9.1110-5	-4.6510-8	1.2510-4
16384	-3.78	-3.78	-1.7210-4	4.5510-5	-1.0810-8	6.2910-5
32768	-3.53	-3.53	-8.0510-5	2.2710-5	-2.5310-9	3.1410-5
65536	-3.29	-3.29	-3.7510-5	1.1310-5	-5.9110-10	1.5710-5
131072	-3.07	-3.07	-1.7510-5	5.6910-6	-1.3710-10	7.8610-6

Результаты расчетов и оценок удобно представить в виде графика (рис. 2.1), где по оси абсцисс отложены десятичные логарифмы n, а по оси ординат – десятичные логарифмы абсолютных величин погрешностей (десятичные логарифмы более удобны с точки зрения наглядности, так как легко определяются порядки погрешностей и чиселn). В таком представлении зависимости близки к линейным.

На рис. 2.1 цифрой 0обозначена зависимость погрешностейz_nzрассчитанных результатов, цифрами (1)-(5) обозначены результаты последовательного вычитания изz_nzодного, двух и трех и т.д. слагаемыхсогласно (2.25). Тем самым, рис. 2.1 есть иллюстрация того, как должен работать идеальный метод экстраполяции.

Ограничение графиков на уровне, соответствующем величине погрешности порядка нескольких единиц 17-го разряда объясняется ограниченной разрядностью чисел в машинном преставлении.

а)

б)

Рис. 2.1. Результаты «идеальной экстраполяции» а) – при удвоении числа слагаемых; б) – при увеличении числа слагаемых на единицу.

а)

б)

Рис. 2.2. Результаты экстраполяции данных, полученных при удвоении числа слагаемых, с помощью метода Ромберга: а) – сравнение с эталоном; б) – оценки по правилу Рунге.

На рис. 2.2а приведены результаты применения повторной экстраполяции по методу Ромберга (2.11) при Q=2. Цифрой1обозначена зависимость погрешностей результатов, экстраполирован­ных один раз, цифрой2– два раза и т. д. Величинаk₁при оценке равнялось 0.1 и далее при возрастании номераjзначениеk_jувеличивалось на 1. Поскольку согласно (2.25) слагаемое сk_j=3.1 отсутствует, прямые3и4совпадают. Из рисунка видно, что путем повторной экстраполяции получаются хотя и не идеальные, но весьма точные результаты. При этом наклон прямых на рис. 1а и 2а совпадает, однако прямые, полученные экстраполяцией, смещены вниз относительно идеальных. Это объясняется погрешностью экстраполяции, оценка которой дается формулой (2.22).

Примечание.Если точное значение искомой величины неизвестно (а это при практических расчетах так и есть), то в таблицах и графиках для сравнения вместо точного можно использовать наиболее точное, полученное на последнем этапе экстраполированное значение. Более того, графики зависимостей логарифма оценок погрешности по Рунге lgz_n^**z_n^* (рис. 2.2б) практически не отличаются от приведенных на рис. 2.2а. Тем не менее, имеет место разница (на один - два десятичных порядка) в величине предельного уровня погрешности округления, связанной с машинным представлением чисел.

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

Если величина погрешности округления сравнима с погрешностью численного метода, то предположение о малости (n) в (2.4) теряет силу. Относительная размытостьоценки погрешности в этом случае имеет значение около единицы и условие (2.16) не выполняется. Однако, погрешности аппроксимации, превышающие погрешность округления на один порядок или больше оцениваются весьма точно.

Для изучения возможности получения наибольшего уточнения результатов при малых nпредставляет интерес рассмотреть последовательность сумм, полученную при увеличении числа слагаемых на единицу. Для экстраполяции при этом необходимо численное решение системы линейных алгебраических уравнений типа (2.13). На рис. 2.3а представлены результаты такой экстраполяции расчетов сумм (2.2) при=1.1. Как видно из рисунка, таким способом возможно получение уточнения до 6 знаков приn=10 и до 9 знаков приn=20-30. При удвоении числа слагаемых при такихnполучается только 2-4 верных знака. Однако дальнейшее уточнение при увеличе­нииnна единицу невозможно из-за ошибки округления. В то же время при удвоенииnошибка округления не растет так быстро и позволяет при той же длине мантиссы (порядка 17 десятичных разрядов) уменьшить погрешность до 10^-14- 10^-15приn=1000-10000.

Как промежуточный вариант, возможно строить последователь­ность при увеличении числа слагаемых на некоторое постоянное число большее единицы. В качестве примера на рис. 2.3б приводятся результаты, полученные при увеличении nна 100. Видно, что таким способом можно получить 12 точных знаков приn1000.

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

а)

б)

Рис. 2.3. Результаты экстраполяции данных с помощью метода Ромберга при увеличении числа слагаемых на постоянную величину:

а) – на единицу; б) – на 100.

Таблица Нэвилла (2.12), применяется для тех случаев, когда показатели разложения (2.13) k_jувеличиваются пропорциональноj, т.е.. В этом случае обозначивt=n^-a, приходим к алгебраическому многочлену или степенному ряду. Номераnчленов последовательности должны увеличиваться на константу. Поэтому рассмотрим следующий пример

(2.26)

Результаты применения экстраполяции с помощью таблицы Нэвилла к этой последовательности показаны на рис. 2.4а. Поскольку последовательность (2.26) представляет алгебраический многочлен относительно n^-1, то результаты второй экстраполяции дают теоретически точное значение. Остаточная погрешность вызвана ошибкой округления.

Следует отметить следующую особенность методов экстраполяции, предназначенных для использования членов последовательностей, номера которых увеличиваются на константу. Если не рассматривать один или несколько первых членов последовательности и перенумеровать остальные, то результаты экстраполяции могут сильно измениться (рис. 2.4б).

а)

б)

Рис. 2.4. Результаты экстраполяции последовательности (2.26) с помощью таблицы Нэвилла, начинающейся: а) –с 1-го номера; б) – со 2-го номера.

В результате рассмотрения примеров можно сформулировать следующие выводы.

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

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

Процесс экстраполяции должен контролироваться критерием, выполнение которого позволяет в некотором смысле гарантировать справедливость результатов.