Трансформация Ван Вейнгаардена

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

Один из алгоритмов вычисления преобразования Эйлера работает следующим образом:

Вычислить строку частичных сумм и сформировать строки средних значений между соседями. Затем первый столбец содержит частичные суммы преобразования Эйлера. с 0 , к = н = 0 к ( 1 ) н а н {\displaystyle s_{0,k}=\sum _{n=0}^{k}(-1)^{n}a_{n}} с дж + 1 , к = с дж , к + с дж , к + 1 2 {\displaystyle s_{j+1,k}={\frac {s_{j,k}+s_{j,k+1}}{2}}} с дж , 0 {\displaystyle s_{j,0}}

Вклад Адриана ван Вейнгаардена состоял в том, что он указал на то, что лучше не доводить эту процедуру до самого конца, а остановиться на двух третях пути. [1] Если доступны, то почти всегда является лучшим приближением к сумме, чем . Во многих случаях диагональные члены не сходятся в одном цикле, поэтому процесс усреднения должен быть повторен с диагональными членами, сводя их в ряд. (Например, это понадобится в геометрической прогрессии с отношением .) Этот процесс последовательного усреднения среднего значения частичной суммы можно заменить, используя формулу для вычисления диагонального члена. а 0 , а 1 , , а 12 {\displaystyle а_{0},а_{1},\ldots ,а_{12}} с 8 , 4 {\displaystyle s_{8,4}} с 12 , 0 {\displaystyle s_{12,0}} 4 {\displaystyle -4}

Для простого, но конкретного примера вспомним формулу Лейбница для числа Пи.

Описанный выше алгоритм создает следующую таблицу:

Вычисление преобразования Эйлера ( 1 ); [2] выделенные значения являются окончательными результатами
1.000000000.666666670.866666670,723809520,834920630,744011540,820934620,754267950,813091480,760459900,808078950,764600690,80460069
0.833333330.766666670,795238100,779365080,789466090,782473080,787601290,783679720,786775690,784269430,786339820,78460069
0.800000000,780952380,787301590,784415580,785969590,785037190,785640500,785227710,785522560,785304630,78547026
0,790476190,784126980,785858590,785192590,785503390,785338840,785434100,785375130,785413590,78538744
0,787301590,784992780,785525590,785347990,785421110,785386470,785404620,785394360,78540052
0,786147190,785259190,785436790,785384550,785403790,785395550,785399490,78539744
0,785703190,785347990,785410670,785394170,785399670,785397520,78539847
0,785525590,785379330,785402420,785396920,785398600,78539799
0,785452460,785390870,785399670,785397760,78539829
0,785421660,785395270,785398710,78539803
0,785408470,785396990,78539837
0,785402730,78539768
0,78540021

Они соответствуют следующим алгоритмическим выводам:

Точность конечных результатов
АлгоритмТермин используемыйЗначение для π / 4 {\displaystyle \пи /4} Относительная погрешность
Наивные частичные суммы с 0 , 12 {\displaystyle s_{0,12}} 0,8046006...+2,4%
преобразование Эйлера с 12 , 0 {\displaystyle s_{12,0}} 0,7854002...+2,6 × 10 −6
результат ван Вейнгаардена с 8 , 4 {\displaystyle s_{8,4}} 0,7853982...+4,7 × 10 −8

Ссылки

  1. ^ А. ван Вейнгаарден , в: Cursus: Wetenschappelijk Rekenen B, Process Analyse, Stichting Mathematich Centrum, (Амстердам, 1965), стр. 51-60
  2. ^ Значения, рассчитанные с помощью выражения J'b11.8'8!:2-:&(}:+}.)^:n+/\(_1^n)*%1+2*n=.i.13

Смотрите также

Взято с "https://en.wikipedia.org/w/index.php?title=Преобразование_Ван_Вейнгаардена&oldid=1187354134"