CC BY
53
УДК 519.651.5
В.И. Шмойлов, Г.А. Кириченко
ОПРЕДЕЛЕНИЕ ЗНАЧЕНИЙ РАСХОДЯЩИХСЯ НЕПРЕРЫВНЫХ
ДРОБЕЙ И РЯДОВ
Рассматривается иное, нежели традиционное, определение сходимости непрерывных дробей. Новый метод суммирования используется при определении значений расходящихся в классическом смысле непрерывных дробей и рядов. Метод суммирования применим не только к обыкновенным непрерывным дробям, но и к непрерывным дробям иных классов, например, к непрерывным дробям Хессенберга. Предложенный алгоритм суммирования расходящихся непрерывных дробей может быть использован при решении разнообразных задач вычислительной математики, в частности, для построения эффективных итерационных алгоритмов решения СЛАУ. Рассматриваемый в статье метод суммирования позволяет найти не только действительные, но и комплексные решения бесконечных систем линейных алгебраических уравнений, если они имеются, что не обеспечивают известные алгоритмы решения БСЛАУ.
Непрерывные дроби; расходящиеся ряды; сходимость непрерывных дробей; r/q-алгоритм.
V.I. Shmoylov, G.A. Kirichenko
DETERMINATION OF THE VALUES OF DIVERGENT CONTINUOUS FRACTIONS AND SERIES
Is considered other than traditional definition of convergence of the continuing rise of the fractions. A new method of summation is used in determining the values of forking in the classical sense of continued fractions and series. Summation method is applicable not only to ordinary continued fractions, but also to the continued fractions of other classes, for example, to the continued fractions Heisenberg. The offered algorithm of summation of divergent continued fractions can be used in solving various problems of computational mathematics, in particular, to build efficient iterative algorithms for solving SLAE. The article summation method allows you to find not only real, but also complete solutions of infinite systems of linear algebraic equations, if they are available, which do not provide the known algorithms of solution of BSLAE.
Continuous fractions; divergent series; convergence continuous fractions; r/q-algorithm.
Введение. Широкое использования непрерывных дробей в вычислительной математике в значительной степени обуславливается тем, что непрерывные дроби в большинстве случаев дают гораздо более общие представления трансцендентных функций, чем степенные ряды. Непрерывные дроби могут быть с большим эффектом использованы для ускорения сходимости медленно сходящихся рядов. Более того, преобразуя расходящиеся ряды в соответствующие непрерывные дроби, нередко можно просуммировать, расходящиеся ряды. Известно, что непрерывные дроби тесно связаны с аппроксимациями Паде, которые, как отмечается в [1], стали главным вычислительным средством в задачах статической механики и физики твердого тела. Поэтому существенные результаты, полученные в теории непрерывных дробей, в частности, в вопросах сходимости, могут быть использованы и в аппроксимациях Паде.
Бесконечной цепной дробью, или непрерывной дробью, называют выражение
вида
а1
к+-1а
Ъ +-
к2 + ■ , ап
■ + -
Ъ +:
где а; и Ъг , г = 1, 2,... - в общем случае независимые переменные.
Часто непрерывную дробь записывают компактно в форме Гершеля:
Ъ + а1
0 Ъ, + Ъ2 + ... + Ъп + ... .
Непрерывная дробь называется сходящейся, если последовательность ее подходящих дробей имеет конечный предел. Непрерывная дробь называется расходящейся, если последовательность ее подходящих дробей предела не имеет. Имеется большое количество признаков сходимости, при помощи которых можно сказать, что существует предел последовательности подходящих дробей или нет. Наиболее широкое применение, пожалуй, получил достаточный признак Ворпиц-кого. По признаку Ворпицкого непрерывная дробь
а1 а2 ап
1 + 1 +... +1 +...
сходится, если |ап| < 1/4, п = 2,3... .
В статье будет рассмотрено несколько задач из разных разделов вычислительной математики, решенных при помощи так называемого т/ф - алгоритма, -нового метода суммирования расходящихся непрерывных дробей. Этот алгоритм формируется следующим образом:
Непрерывная дробь сходится и имеет своим значением в общем случае комплексное число г = т0 е , если существуют пределы
Иш
1
ПI Р / I = То
г=1 , (1)
яИш = |фо| (2)
где Р1 ^ - значения 1-й подходящей дроби из совокупности, включающей 8 подходящих дробей; к8 - число отрицательных подходящих дробей из 8 подходящих дробей.
Этот алгоритм применим как к обыкновенным непрерывным дробям, так и к непрерывным дробям с графами иных типов, например, к непрерывным дробям Хессенберга и ветвящимся непрерывным дробям.
В случае непрерывных дробей, сходящихся в классическом смысле, аргумент ф0 примет значения 0 или я . Если ф0 = 0, то значение сходящейся непрерывной
дроби будет совпадать со значением модуля т0 :
г 0
г = т е = Т0.
2
Если р0 = 7, то значение сходящейся непрерывной дроби будет отрицательное число:
г = гоe = ~гО .
Проиллюстрируем эффективность предложенного способа суммирования расходящихся в классическом смысле непрерывных дробей решением ряда задач.
Рис. 1. Распределение значений подходящих непрерывной дроби логарифмической
функции
2. Суммирование расходящихся непрерывных дробей и рядов. Известно, что непрерывные дроби целесообразнее использовать для аппроксимации функций, нежели степенные ряды, так как непрерывные дроби зачастую сходятся в более широкой области. Например, ряд Меркатора
X2 х3 х4
1п(1 + х) = X--+---+ ...
2 3 4
представляет логарифмическую функцию в единичном круге, в то время как непрерывная дробь Лагранжа
х х х 2 х 2 х пх пх /о \
1п(1 + х) = — — — - - - --(3)
1 + 2 + 3 + 2 + 5 + ... + 2 + 2 п + 1 + ...
сходится к функции 1п(1 + х) на всей плоскости комплексного переменного, за исключением выреза от -1 до — [2].
При отрицательных значениях аргумента логарифмическая функция имеет комплексное значение и, естественно, что непрерывные дроби, получающиеся из непрерывной дроби Лагранжа при х <—1 будут расходящимися в классическом смысле.
На рис. 1 (а,б,в,г,д,е) показано распределение подходящих значений непрерывной дроби Лагранжа (3) при х = -10, -100, -1000 на начальных участках (п = 1 -100 и п = 1 ■ 500).
В табл. 1 показаны результаты суммирования расходящейся непрерывной дроби
1п(—2) = —3 3 3 6 6 ^ ^ (4)
1 — 2 — 3 — 2 — 5 —... — 2 — 2п +1 —... .
Таблица 1
Определение значения расходящейся непрерывной дроби (4) т0 = 3.217150511-7.. ф0 = 1.3536398454...
№ звена дроби Значение подходящей дроби Модуль комплексного числа, т ' п Погрешность, £т = |Т0 — Тп\ ш1п £т Аргумент комплексного числа, ф Погрешность, £Ф=ф0 —ф\ ш1п
1 2 3 4 5 6 7 8
1 -3,0000000 3,0000000000 0,2171505117 3,1415926535 1,7879528081 ш
2 6,0000000 4,2426406871 1,0254901754 1,5707963267 0,2171564813 ш
4 -3,0000000 3,0000000000 0,2171505117 1,5707963267 0,2171564813
8 -97,5000000 4,9614481602 1,7442976485 1,5707963267 0,2171564813
16 1,4880473 3,5474336503 0,3302831386 1,3744467859 0,0208069405 ш
32 3,1985122 3,6050160485 0,3878655367 1,3744467859 0,0208069405
64 62,8693924 3,3885474566 0,1713969449 ш 1,3744467859 0,0208069405
128 0,9165216 3,1810462758 0,0361042359 ш 1,3499030933 0,0037367521 ш
256 1,7095765 3,2148854739 0,0022650377 ш 1,3621749396 0,0085350941
512 3,9037050 3,2112688498 0,0058816618 1,3499030933 0,0037367521
1024 -15,4772571 3,2219262392 0,0047757275 1,3560390164 0,0023991710 ш
2048 2,6358581 3,2194825453 0,0023320336 1,3529710549 0,0006687905 ш
4096 11,1007665 3,2127253440 0,0044251676 1,3529710549 0,0006687905
8192 -0,6961262 3,2169015620 0,0002489496 ш 1,3533545501 0,0002852953 ш
16384 -1,7591587 3,2167104407 0,0004400709 1,3533545501 0,0002852953
32768 -6,4347291 3,2170964982 0,0000540134 1,3536421715 0,0000023260 ш
65536 5,5879135 3,2171496506 0,0000008610 1,3536421715 0,0000023260
131072 -3,9038315 3,2171884212 0,0000379094 1,3536182030 0,0000216423
262144 16,0431708 3,2171480639 0,0000024477 1,3535942346 0,0000456108
524288 -0,0551483 3,2171287791 0,0000217325 1,3536421715 0,0000023260
1048576 -0,2709104 3,2171427009 0,0000078107 1,3536361793 0,0000036660
2097152 -0,7308612 3,2171496552 0,0000008564 ш 1,3536361793 0,0000036660
4194304 -1,8537413 3,2171502478 0,0000002638 ш 1,3536391754 0,0000006699 ш
8388608 -7,2124648 3,2171495794 0,0000009323 1,3536399244 0,0000000790 ш
При отрицательном аргументе логарифмическая функция имеет комплексное значение
1п(-2) = 3,2171505117 е;13536398454 ,
которое, естественно, не может приближаться непрерывной дробью с вещественными элементами и, тем не менее, суммирование при помощи г/ р - алгоритма позволяет установить значение дроби (4).
В первой колонке таблицы даны номера п подходящих дробей разложения
(4). Номера подходящих дробей составляют степень 2: п = 2', ' = 1 ^ 23 . Значения подходящих дробей с этими номерами приведены в соседней колонке 2. Как и следовало ожидать, значения подходящих дробей \Рп/0.п} с ростом п не стремятся к какому-либо пределу. Для чисел же, расположенных в колонке 3, напротив, стремление к пределу можно без труда обнаружить, - значения приближаются к величине 3.2171505117..., т.е. к модулю комплексного числа 1п(-2). Даже беглого взгляда на колонки 6 и 7 достаточно, чтобы убедиться, что с ростом количества подходящих дробей разложения (4) все более точно устанавливается значение аргумента искомого комплексного числа.
Найдены представления элементарных и некоторых специальных функции в виде непрерывных дробей Хессенберга [3]. Например,
ln| 1 + -
х/1! х/ 2! х/ 3! х/ 4! х/ 5
-1 х/1! х/ 2! х/ 3! х/ 4
0 -1 х/1! х/ 2! х/ 3
0 0 -1 х/1! х/ 2
0 0 0 -1 х/1
х/1! х/ 2! х/ 3! х/ 4
-1 х/1! х/ 2! х/ 3
0 -1 х/1! х/ 2
0 0 -1 х/1
(5)
Непрерывная дробь (5) определяет логарифмическую функцию на всей плоскости комплексного переменного с разрезом от 0 до -1. Комплексное значение логарифмической функции на разрезе определяется из (5) суммированием по формулам (1) и (2). В качестве подходящих дробей (5) будем брать последовательность отношений определителей:
х/1! х/2! х/3!
х/1! х/2! -1 х/1!
Qo
х/1! PL_
1 ' Qi _
х/1!
ß2
-1 0
х/1! х/2! -1 х/1!
х/1! х/2! -1 х/1!
1
В табл. 2 приведены результаты вычисления комплексного числа 1/1п(-2) при помощи функциональной непрерывной дроби Хессенберга (5) и алгоритма суммирования расходящихся непрерывных дробей.
Таблица 2
Определение значения непрерывной дроби (5) при х = -1/3.
х = -1/3, г0 = 0.3108340739..., ф0 =-1.3536398454... .
№ звена дроби Значение подходящей дроби Модуль комплексного числа, гп Погрешность, Ег = |Г0 - Гп\ тп Ег Аргумент Комплексного числа, фп Погрешность, Еф = ф0 -Фп тт еф
1 2 3 4 5 6 7 8
20 0,85127256 0,99198117 0,68114710 т -1,41371669 0,06007684 т
21 0,02044264 0,82455439 0,51372032 т -1,34639685 0,00724299 т
42 -0,02832422 0,50977666 0,19894259 т -1,42119667 0,06755683
84 -0,14445825 0,40513373 0,09429965 т -1,38379676 0,03015691
168 -0,75441728 0,35654051 0,04570644 т -1,36509680 0,01145696
336 0,32678464 0,33285231 0,02201823 т -1,35574682 0,00210698 т
672 2,01078357 0,32176153 0,01092745 т -1,35574682 0,00210698
1344 -0,03019484 0,31569144 0,00485737 т -1,35574682 0,00210698
2688 -0,14954713 0,31342156 0,00258748 т -1,35457808 0,00093823 т
5376 -0,81466750 0,31217034 0,00133626 т -1,35399370 0,00035386 т
10752 0,30407891 0,31149855 0,00066447 т -1,35370152 0,00006167 т
21504 1,28357034 0,31116989 0,00033582 т -1,35370152 0,00006167
43008 -0,09454281 0,31099243 0,00015836 т -1,35370152 0,00006167
86016 -0,38364272 0,31091665 0,00008257 т -1,35366499 0,00002515 т
172032 0,81555456 0,31087600 0,00004193 т -1,35364673 0,00000689 т
344064 -0,22761506 0,31085447 0,00002039 т -1,35364673 0,00000689
688128 -10,07498787 0,31084454 0,00001047 т -1,35364217 0,00000232 т
1376256 0,08515530 0,31083902 0,00000495 т -1,35363988 0,00000004 т
2752512 0,10347124 0,31083657 0,00000249 т -1,35363988 0,00000004
5505024 0,14104327 0,31083533 0,00000126 т -1,35363988 0,00000004
11010048 0,22449739 0,31083471 0,00000064 т -1,35363988 0,00000004
Приведём ещё пример нахождения значения расходящейся в классическом смысле непрерывной дроби при помощи описанного выше алгоритма суммирования.
Известна непрерывная дробь Лагранжа для степенной функции [2]:
(1 + ху = 1 + у (1 -у)х (1 + у)х (т-у)х (т + у)х (6) 1 + 2 + 3 +... + 2 + 2т +1 + ..."
Непрерывная дробь (6) сходится на всей плоскости комплексного переменного, разрезанной по вещественной оси от -1 до . При х < -1 подкоренное выражение функции у = (1 + х)У становится отрицательным, следовательно, значение функции (1 + х)У при х < -1 становится комплексной величиной , которая, очевидно, не может приближаться вещественной последовательностью подходящих дробей разложения (6).
В табл. 3 приведены результаты вычисления значения расходящейся непрерывной дроби для 3- 9 при помощи г/ф -алгоритма.
Таблица 3
Определение значения расходящейся непрерывной дроби
-1 _10 2-10 4-10 5-10 (3n-1)10 (3n +1)10 - 3 _ 2 _ 9 _ 2 _..._ 2 _ 3(2n +1) _..."
r0 - 2.080083..., p0 -1.047197....
Номер Значение Модуль Погрешность, min £r Аргумент Погрешность,
звена дроби подходящей дроби комплексного числа, rn £r - |Г0 _ Гп\ комплексного числа, pn £p - p0 _Pn| min EP
1 2 3 4 5 6 7 8
2 2,4285714 8,095238095238 6,015154272186 m 0,000000000000 1,047197551196 m
4 7,9230769 1,892578548704 0,187505274347 m 0,000000000000 1,047197551196
8 0,0083609 0,946908431432 1,133175391618 0,785398163397 0,261799387799 m
16 2,3793179 1,514212477498 0,565871345553 0,785398163397 0,261799387799
32 3,4334933 1,931910874242 0,148172948809 m 0,687223392972 0,359974158223
64 -6,6957404 2,285828469458 0,205744646406 0,932660319034 0,114537232162 m
128 1,8056646 2,171608471371 0,091524648319 m 0,957204011640 0,089993539555 m
256 2,7046939 2,095947975804 0,015864152753 m 1,043106935762 0,004090615434 m
512 14,0355053 2,077080542151 0,003003280900 m 1,043106935762 0,004090615434
1024 0,4131559 2,063579787259 0,016504035792 1,046174897338 0,001022653858 m
2048 -0,5703198 2,085479528230 0,005395705178 1,044640916550 0,002556634646
4096 -33,5321135 2,082858835017 0,002775011965 m 1,046174897338 0,001022653858
8192 1,1197769 2,079900863342 0,000182959709 m 1,048092373322 0,000894822126 m
16384 1,0914866 2,080396700069 0,000312877017 1,047133635330 0,000063915866 m
32768 1,0349376 2,080307362394 0,000223539343 1,047229509129 0,000031957933 m
65536 0,9216825 2,079865087605 0,000218735446 1,047229509129 0,000031957933
131072 0,6914169 2,080106476321 0,000022653270 m 1,047205540679 0,000007989483 m
262144 0,1869493 2,079930141859 0,000153681192 1,047193556454 0,000003994741 m
524288 -1,4500205 2,080071498760 0,000012324291 m 1,047199548567 0,000001997370 m
1048576 5,5741164 2,080081877069 0,000001945982 m 1,047190560398 0,000006990797
2097152 -0,8763266 2,080075907132 0,000007915919 1,047195054483 0,000002496713
4194304 15,7595707 2,080077041783 0,000006781268 1,047199548567 0,000001997370
8388608 0,4758782 2,080083260804 0,000000562247 m 1,047197301525 0,000000249671 m
Следует отметить, что при использовании г/ф -алгоритма необходимо иметь эффективные алгоритмы для построения так называемых соответствующих непрерывных дробей и для вычисления длинных серий значений подходящих дробей. Использование для счёта классического рекуррентного алгоритма (ЕЯ - алгоритм) приводит к быстрому переполнению разрядной сетки компьютера, а применение естественной процедуры вычисления непрерывной дроби "снизу - вверх" (ВЯ -алгоритм) невозможно из-за недопустимо больших временных затрат при определении серий значений подходящих дробей. В [4] были детально рассмотрены алгоритмы для вычисления длинных серий значений подходящих дробей. В [5] изложены различные методы построения соответствующих непрерывных дробей.
Расходящиеся ряды нередко суммируются через соответствующие непрерывные дроби [6]. Для степенного ряда
2 п
С0 + С1 X + С2 X + ... + Спх + ... (7)
можно построить непрерывную дробь
^ ОДх ОДх тъх ОД2п-1х т2пх (8)
С00 Н--- - - - , (8)
0 1 - 1 + 1 -... + 1 - 1 +...
такую, что разложение п-й подходящей непрерывной дроби будет совпадать с исходным рядом (7) вплоть до члена Спхп включительно:
Р (X)
= С0 + С1X + с2 X +... +Спх + уп+1 X + .
а (х)
Такую непрерывную дробь называют соответствующей ряду или соответствующей непрерывной дробью.
Используя коэффициенты С1 степенного ряда (7), можно построить соответствующую непрерывную дробь (8) по формулам Хейлерманна-Стилтьеса [7]:
с0 = с0, с1 = с1,
Рп =
С2п = <Рп-г¥п+1 Рп¥п С2п+1 = Рп+¥п . Р¥п+1 '
С1 С2 ... Сп С2 Сз . .. Сп
С2 С3 ... Сп+1 , ¥п = С3 С4 . .. Сп+1
Сп Сп+1 ... С2 п—1 Р<> = 1, ¥1 = 1. Сп Сп+1 . .. С2п-2
(9)
Прежде чем рассматривать расходящиеся ряды, приведём пример ускорения медленно сходящихся рядов при помощи соответствующих непрерывных дробей. Предложено много алгоритмов ускорения сходимости рядов и последовательностей. Как показывает практика, пожалуй, самым эффективным и надёжным способом ускорения сходимости медленно сходящихся рядов является их суммирование через соответствующие непрерывные дроби. В табл. 4 приведены результаты вычисления 1п2 с использованием ряда Меркатора.
Таблица 4
Определение значения ряда Меркатора
,„,1111 1п2 = 1--+---+ --....
2 3 4 5
Число членов Значение частичных Погрешность
ряда сумм ряда аппроксимации
2 0,50000000000000 0,19314718055994
10 0,64563492063492 0,04751225992502
100 0,68817217931020 0,00497500124974
1000 0,69264743055982 0,00049975000012
10000 0,69309718305996 0,00004999749998
100000 0,69314218058498 0,00000499997496
1000000 0,69314668056025 0,00000049999969
10000000 0,69314713056010 0,00000004999984
100000000 0,69314717556042 0,00000000499952
В табл. 5 константа 1п2 определена при помощи цепной дроби Лагранжа.
Таким образом, использование непрерывной дроби, содержащей 19 звеньев, обеспечивает точность вычисления 1п2 с 14 десятичными разрядами. Для сравнения: применение для определения 1п2 ряда Меркатора, включающего 100 миллионов членов, позволяет вычислить 1п2 с точностью всего 8 десятичных знаков. Разница в эффективности аппроксимации фантастическая.
Таблица 5
Определение значения непрерывной дроби Лагранжа 1112 2 n n
In 2 -
1 + 2 + 3 + 2 + 5 +... + 2 + 2n +1 +...
Число Значение Погрешность
звеньев подходящих дробей аппроксимации
2 0,66666666666666 0,02648051389328
3 0,70000000000000 0,00685281944005
4 0,69230769230769 0,00083948825225
5 0,69333333333333 0,00018615277339
6 0,69312169312169 0,00002548743825
7 0,69315245478036 0,00000527422042
8 0,69314641744548 0,00000076311446
9 0,69314733235438 0,00000015179444
10 0,69314715785304 0,00000002270690
11 0,69314718496213 0,00000000440219
12 0,69314717988653 0,00000000067341
13 0,69314718068816 0,00000000012822
14 0,69314718054001 0,00000000001993
15 0,69314718056369 0,00000000000375
16 0,69314718055936 0,00000000000058
17 0,69314718056005 0,00000000000012
18 0,69314718055993 0,00000000000001
19 0,69314718055995 0,00000000000001
Области сходимости степенного ряда (7) и соответствующей непрерывной дроби (8) для одной и той же функции могут быть различны. Например, расходящийся ряд Эйлера
1 -1!х + 2!х2 - 3!х3 + 4!х4 -... (10)
имеет соответствующую непрерывную дробь
1 1х 1х 2х 2х пх пх (11)
I+т+т+т+т+...+т+т+....
г е-
Ряд (8) представляет функцию I -йг , которая связана с интегральной
0 1 + хг
показательной функцией Ег(х) [8]. Числовой ряд
С0 + С1 + С2 + ... + Сп + ... (12)
суммируется через значение непрерывной дроби (8) при х=1, если ряд (12) рассматривать как степенной ряд (7) при х = 1.
Расходящийся всюду, кроме х = 0, ряд Эйлера при х = 1 принимает вид
1-1! + 2! - 3! + 4!-... . (13)
Определяя значение расходящегося ряда (13) через соответствующую непрерывную дробь (11) при х = 1, получим:
11."..-..,.. 1112 2 п п
1 -1!+2!—3!+4!—... = - - - — — — — = 1 + 1 + 1 + 1 + 1 +... + 1 + 1 + ...
=0,5963473623231945... (14)
Выше уже отмечалось, что ряд Эйлера (10) связан с интегральной показательной функцией Е1(х) .Действительно известен быстро сходящийся ряд:
-е£г(-1) = - е\ С - 1 + —---— + ... | = 0,5963473623231945... ,
{ 2 • 2! 3 • 3! )
где е - неперово число, равное 2.718281... ,
С - постоянная Эйлера, имеющая значение 0,577215... .
Определив по расходящимся рядам соответствующие непрерывные дроби, можем найти значения других расходящихся рядов. Например, 1 + 1 - 1 + 2 - 5 + 14 - 42 + 132-.=
= 1 +1 1 1 = = 1,618033 ... . (15) 1 + 1 +... + 1 +... 2
1 -1 + 1- 3-1-3-5 + 1- 3 • 5 • 7-1-3-5 • 7 • 9 +... =1 123 п = 0.655679 (16)
1 +1 + 1 + 1 +... + 1 +...
Расходящийся ряд (15) имеет своим значением "отношение золотого сечения", а расходящийся ряд (16) связан с неполной гамма-функцией, так как
112 3 п [ёт/1 1 ^ Лё ,( 1 ^ ---- - =, -Г|-,- I = .—еф\-;=\ = 0.6556795424....
1 +1 + 1 + 1 +... + 1 +... \2 { 2 2 ) \2 {у/2 )
Расходящимся рядам могут соответствовать конечные непрерывные дроби, т.е. эти ряды связаны с рациональными функциями:
1 -1 +1 -1 +1 -... = 1 -1 1 = -1, 1 +1 2
2 2
1 + 2 + 4 + 8 +16 +... = 1 +--= -1,
1 -1
1 ^ О X С ,23121 1 - 2 + 3 - 4 + 5 -... = 1 — - - - = -.
1 - 2 - 3 + 1 4
Нередки случаи, когда соответствующая непрерывная дробь, построенная для расходящегося ряда, не сходится в классическом смысле. Построим для ряда
1 +1! + 2! + 3! + 4!+... (17)
соответствующую непрерывную дробь. По формулам (9) получим:
1 +1! + 2!+ 3! + 4! +... = 1 1 1 2 2 п 1 (18)
1 -1 -1 - 1 - 1 -... - 1 - 1 -... .
Непрерывная дробь (18), в отличие от непрерывной дроби (14), имеет отрицательные частные числители. Непрерывная дробь (18) по признаку Ворпицкого является расходящейся непрерывной дробью. В самом деле, значения подходящих дробей разложения (18) не стремятся ни к какому пределу, т.е., непрерывная дробь (18) является в классическом смысле расходящейся.
Найдём значение расходящейся непрерывной дроби (18) при помощи г /ф-алгоритма. В табл. 6 приведены результаты вычислений непрерывной дроби (18), представляющей Е1(1).
™ (1112233 п п Л £г(1) = е--- — — - - — —
-1 -1 - 1 - 1 - 1 - 1 -... - 1 - 1 -...)
= 3.6689338967 е™017377
Таблица 6
Определение значения расходящейся непрерывной дроби г0 = 3.6689338967?. , ф0 = 1.028001737!..
№ Значение Модуль Погрешность, Аргумент Погрешность,
звена дроби подходящей дроби комплексного числа, rn = ко " Гп\ комплексного числа, <pn £< = К -Фп\
1 2 3 4 5 6
4 5,43656365 3,23260009 0,43633380 0,00000000 1,02800173
6 1,08731273 3,84423102 0,17529713 0,39269908 0,63530265
16 1,26199375 4,12365173 0,44717838 0,78539816 0,24260357
32 -3,22229391 3,21557146 0,45336243 0,88357293 0,14442880
64 -5,44198600 3,40441763 0,26451626 0,98174770 0,04625403
128 3,40902451 3,44911711 0,21981677 1,03063508 0,00283335
256 2,26696092 3,60976107 0,05917282 1,04310693 0,01510519
512 2,46414684 3,65055907 0,01837481 1,04310693 0,01510519
1024 1,67450019 3,69927410 0,03034020 1,02163120 0,00637053
2048 5,89710924 3,70773799 0,03880409 1,02623314 0,00176889
4096 1,76596456 3,69848222 0,02954832 1,03236907 0,00436733
8192 -2,16761507 3,67756787 0,00863398 1,03428654 0,00628480
16384 -5,45891725 3,66492493 0,00400896 1,02968460 0,00168286
32768 -1,02199219 3,66736597 0,00156791 1,02728775 0,00071397
65536 1,98286910 3,67009386 0,00115907 1,02968460 0,00168286
131072 5,57106207 3,66798677 0,00094711 1,02721585 0,00078588
262144 -8,97247936 3,66782858 0,00110530 1,02751545 0,00048627
524288 6,23221177 3,66849563 0,00043826 1,02756938 0,00043234
1048576 5,89027864 3,66844670 0,00048718 1,02773716 0,00026457
2097152 8,12082269 3,66953895 0,00060305 1,02787255 0,00017918
4194304 7,05859556 3,66906517 0,00013127 1,02781431 0,00018742
8388608 2,28589509 3,66933307 0,00039917 1,02805737 0,00005563
Таким образом, удалось просуммировать расходящийся ряд (17), связанный с интегральной показательной функцией:
1 Ei(1) = 1 +1!+ 2!+ 3!+ 4!+ ... =1 1 1 2 2 3 3 П П = е 1 -1 -1 - 1 - 1 - 1 - 1 -... - 1 - 1 -... . (19)
= 1,34970е'1 02800
Аналогичным образом, преобразуя ряды в соответствующие непрерывные дроби с помощью формул (9), а затем находя значения расходящихся в классическом смысле непрерывных дробей при помощи r / < -алгоритма, т.е. формул (1) и (2), можно просуммировать и другие расходящиеся ряды. Например,
111 *
1 +1 +1 + 2 + 5 +14 + 42 +132 + ... = 1--- - = е 3, (20)
1 -1 -... -1 -...
1 +1 +1 • 3 +1 • 3 • 5 +1 • 3 • 5 • 7 +... = I 1 2 3 1 П =
1-1-1-1-1-...-1-...
1,050318 е10 809229. (21)
Непрерывная дробь (20) представляет комплексный корень квадратного уравнения
2 Я
X - 2 С08— X +1 = 0. 3
Непрерывная дробь (21) связана с неполной гамма функцией:
112 3 1 П = ; &(!-!) = , 1^ег/е(-^) =
1-1-1-1-1-...-1-... V 2е 2 2 \2е 42
1,0503178040е'08092294339. (22)
В табл. 7 приведены результаты определения значения расходящейся непрерывной дроби (22).
Таблица 7
Определение значения расходящейся непрерывной дроби
Л! г[1,-11 = Л* ег1с[+\ = 1 1 2 3 4 П .
\2е {2 2) \2е ^ д/2) 1-1-1-1-1 -...- 1 -...
г0 = 1.050317804055... , ф0 = 0.80922943390... .
Номер Значение Модуль Комплексного числа, Гп Погрешность, шш Аргумент комплексного числа, фп ПогрешноСть,
звена дроби подходящей дроби = |Г0 -Гп\ £ф =ф0 -Фп| ЕФ
1 2 3 4 5 6 7 8
8 0,909090 1,018429678599 0,031888125466 ш 0,897597901025 0,088368467122 ш
16 -0,255298 1,017392202149 0,032925601915 0,837758040957 0,028528607054 ш
32 1,232969 1,051369832722 0,001052028657 ш 0,709391889520 0,099837544382
64 5,097151 1,188307933456 0,137990129391 0,847731350968 0,038501917065
128 2,897390 1,137523550255 0,087205746189 0,816319350932 0,007089917029 ш
256 0,485401 1,099518268902 0,049200464837 0,813118098576 0,003888664673 ш
512 0,152295 1,073023887296 0,022706083231 0,786935146104 0,022294287798
1024 0,215155 1,069492823824 0,019175019759 0,798449745780 0,010779688122
2048 -1,789037 1,045731096537 0,004586707527 0,813406988960 0,004177555057
4096 -1,076242 1,043640648366 0,006677155698 0,810906821695 0,001677387792 ш
8192 1,238412 1,046287186312 0,004030617752 0,805821775752 0,003407658150
16384 1,513253 1,050413791934 0,000095987869 ш 0,804813792780 0,004415641123
32768 2,634049 1,053561869047 0,003244064982 0,809583067321 0,000353633417 ш
65536 -18,763236 1,050054380996 0,000263423069 0,811536156750 0,002306722846
131072 0,007903 1,050452714611 0,000134910545 0,807694983945 0,001534449958
262144 0,785107 1,049839458591 0,000478345473 0,810292489546 0,001063055643
524288 -12,041543 1,050061241159 0,000256562905 0,809985345715 0,000755911812
1048576 0,846676 1,049973494523 0,000344309541 0,809607069807 0,000377635904
2097152 0,815962 1,050125044812 0,000192759252 0,809540770488 0,000311336585 ш
4194304 0,975265 1,050106631988 0,000211172077 0,809297147852 0,000067713949 ш
8388608 0,909983 1,050176943924 0,000140860141 0,809317649265 0,000088215362
16777216 1,286952 1,050235875462 0,000081928602 ш 0,809141582712 0,000087851190
Заключение. Метод суммирования расходящихся непрерывных дробей может быть эффективно использован при решении бесконечных нерегулярных систем линейных алгебраических уравнений [9]. Под нерегулярными понимают системы алгебраических уравнений, решения которых не стремятся к пределу в классическом смысле с ростом размерности системы. Известно, что решения трёхдиа-гональных систем линейных алгебраических уравнений можно представить в виде непрерывных дробей. Этот способ записи решений трёхдиагональных систем эк-
вивалентен алгоритму "прогонки". Непрерывные дроби, представляющие решения системы, будут бесконечными, если бесконечна система линейных алгебраических уравнений. Но бесконечные непрерывные дроби могут быть сходящимися и расходящимися. При решении регулярных систем не возникает трудностей, так как вместо бесконечных систем оперируют "усечёнными" системами. Непосредственно метод "усечения", очевидно, не срабатывает в случае нерегулярных бесконечных систем. При некоторых ограничениях решения нерегулярных бесконечных систем находятся при помощи рассмотренного выше метода суммирования расходящихся непрерывных дробей. Можно указать ещё ряд задач в вычислительной математике, эффективное решение которых обеспечивается алгоритмом суммирования расходящихся непрерывных дробей [10-12].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Бейкер Дж., Грейвис-Морис П. Аппроксимация Паде: Пер. с апгл. - М.: Мир, 1986.
- 502 с.
2. Джоунс У., Трон В. Непрерывные дроби. Аналитическая теория и приложения: Пер. с апгл.- М.: Мир, 19S5. - 414 с.
3. Шмойлов В.И. Периодические цепные дроби. - Львов: Академический Экспресс, 199S.
- 219 с.
4. Качмар В.С., Русин Б.П., Шмойлов В.И. Алгоритмы вычисления значений цепных дробей II Вычислительная математика. - 199S. - Т. 3S, № 9. - С. 143б-1451.
5. Шмойлов В.И. Непрерывные дроби. В 3-х т. Т. 2. Расходящиеся непрерывные дроби. НАН Украины, Ип-т прикл. пробл. механики и математики. - Львов: Меркатор, 2004.
- 55S с.
6. Brezlnskl C. History of continued fraction and Pade approximants. - Springer- Verlag, Berlin, 1991. - 547 p.
7. Lorentzen L., Waadeland H. Continued fractions with applications. - Amsterdam - London -New-York - Tokyo, 1992. - б0б p.
S. Эрдейи А. Асимптотические разложения: Пер. с. апгл. - М.: ГИФМЛ, 19б2. - 127 с.
9. Шмойлов В.И. Расходящиеся системы линейных алгебраических уравнений. - Таганрог: Изд-во ТТИ ЮФУ, 2010. - 205 с.
10. Шмойлов В.И. Решение алгебраических уравнений при помощи ^ф-алгоритма. - Таганрог: Изд-во ТТИ ЮФУ, 2012. - 330 с.
11. Шмойлов В.И. Непрерывные дроби и г!ф алгоритм. - Таганрог: Изд-во ТТИ ЮФУ, 2012. - 60s с.
12. Шмойлов В.И., Коваленко В.Б. Некоторое применения алгоритма суммирования расходящиеся непрерывных дробей II Вестпик Южного паучпого центра РАН. - 2012. - Т. S, № 4. - C. 3-13.
Статью рекомендовал к опубликованию д.т.н. И.И. Левин.
Шмойлов Владимир Ильич - Южный научный цептр РАН; e-mail: [email protected]; 347902, г. Таганрог, ул. Свободы, 27-6, кв. 6; тел.: 88634318910, 88634368337, факс: 88634360376; научный сотрудник.
Кириченко Геннадий Анатольевич - Федеральное государственное автономное образова-тельпое учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: [email protected]; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: 88634371428, 89064287987; аспирант
Shmoylov Vladimir Ilyich - Southern Scientific Center of Russian Academy of Sciences (SSC RAS); e-mail: [email protected]; 27-6, Svobody street, kv. 6, Taganrog, 347902, Russia; phones: +78634318910, +78634368337, fax: +78634360376; researcher.
Kirichenko Gennady Anatolievich - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education "Southern Federal University"; e-mail: [email protected]; 44, Nekrasovsky, Taganrog, 347928, Russia; phones: +78634371428, +79064287987; postgraduate student.
