Articulo de referencia

El algoritmo de Borwein

El algoritmo de Borwein fue ideado por Jonathan y Peter Borwein para calcular el valor de 1 / π {\displaystyle 1/\pi } Este y otros algoritmos se pueden encontrar en el libro Pi...

El algoritmo de Borwein fue ideado por Jonathan y Peter Borwein para calcular el valor de1/π{\displaystyle 1/\pi }Este y otros algoritmos se pueden encontrar en el libro Pi y la AGM: un estudio de teoría analítica de números y complejidad computacional . [ 1 ]

Serie Ramanujan-Sato

Estos dos son ejemplos de una serie de Ramanujan-Sato . El algoritmo de Chudnovsky relacionado utiliza un discriminante con número de clase 1.

Clase número 2 (1989)

Comience configurando [ 2 ]

A=21217571091261+1657145277365B=1377398089267261+107578229802750do=(5280(236674+3030361))3{\displaystyle {\begin{aligned}A&=212175710912{\sqrt {61}}+1657145277365\\B&=13773980892672{\sqrt {61}}+107578229802750\\C&=\left(5280\left(236674+30303{\sqrt {61}}\right)\right)^{3}\end{aligned}}}

Entonces

1π=12norte=0(1)norte(6norte)¡(A+norteB)(norte¡)3(3norte)¡donorte+12{\displaystyle {\frac {1}{\pi }}=12\sum _{n=0}^{\infty }{\frac {(-1)^{n}(6n)!\,(A+nB)}{(n!)^{3}(3n)!\,C^{n+{\frac {1}{2}}}}}}

Cada término adicional de la suma parcial produce aproximadamente 25 dígitos.

Clase número 4 (1993)

Comience configurando [ 3 ]

A=63365028312971999585426220+283377021408008420468256005+3845(10891728551171178200467436212395209160385656017+48709290865788102250773385345416887213512550405)12B=7849910453496627210289749000+35105866782609320289656064005+25159683110(6260208323789001636993322654444020882161+27996502730604442965772068907188251902355)12do=21477299506351224096049403338648032512965(10985234579463550323713318473+49127462536923627546073959125)12{\displaystyle {\begin{aligned}A={}&63365028312971999585426220\\&{}+28337702140800842046825600{\sqrt {5}}\\&{}+384{\sqrt {5}}{\big (}10891728551171178200467436212395209160385656017\\&{}+\left.4870929086578810225077338534541688721351255040{\sqrt {5}}\right)^{\frac {1}{2}}\\B={}&7849910453496627210289749000\\&{}+3510586678260932028965606400{\sqrt {5}}\\&{}+2515968{\sqrt {3110}}{\big (}6260208323789001636993322654444020882161\\&{}+\left.2799650273060444296577206890718825190235{\sqrt {5}}\right)^{\frac {1}{2}}\\C={}&-214772995063512240\\&{}-96049403338648032{\sqrt {5}}\\&{}-1296{\sqrt {5}}{\big (}10985234579463550323713318473\\&{}+\left.4912746253692362754607395912{\sqrt {5}}\right)^{\frac {1}{2}}\end{aligned}}}

Entonces

do3π=norte=0(6norte)¡(3norte)¡(norte¡)3A+norteBdo3norte{\displaystyle {\frac {\sqrt {-C^{3}}}{\pi }}=\sum _{n=0}^{\infty }{{\frac {(6n)!}{(3n)!(n!)^{3}}}{\frac {A+nB}{C^{3n}}}}}

Cada término adicional de la serie produce aproximadamente 50 dígitos.

Algoritmos iterativos

Convergencia cuadrática (1984)

Comience configurando [ 4 ]

a0=2b0=0pag0=2+2{\displaystyle {\begin{aligned}a_{0}&={\sqrt {2}}\\b_{0}&=0\\p_{0}&=2+{\sqrt {2}}\end{aligned}}}

Luego itera

anorte+1=anorte+1anorte2bnorte+1=(1+bnorte)anorteanorte+bnortepagnorte+1=(1+anorte+1)pagnortebnorte+11+bnorte+1{\displaystyle {\begin{aligned}a_{n+1}&={\frac {{\sqrt {a_{n}}}+{\frac {1}{\sqrt {a_{n}}}}}{2}}\\b_{n+1}&={\frac {(1+b_{n}){\sqrt {a_{n}}}}{a_{n}+b_{n}}}\\p_{n+1}&={\frac {(1+a_{n+1})\,p_{n}b_{n+1}}{1+b_{n+1}}}\end{aligned}}}

Entonces, p k converge cuadráticamente a π ; es decir, cada iteración duplica aproximadamente el número de dígitos correctos. El algoritmo no es autocorrectivo; cada iteración debe realizarse con el número deseado de dígitos correctos para el resultado final de π .

Convergencia cúbica (1991)

Comience por configurar

a0=13s0=312{\displaystyle {\begin{aligned}a_{0}&={\frac {1}{3}}\\s_{0}&={\frac {{\sqrt {3}}-1}{2}}\end{aligned}}}

Luego itera

rk+1=31+2(1sk3)13sk+1=rk+112ak+1=rk+12ak3k(rk+121){\displaystyle {\begin{aligned}r_{k+1}&={\frac {3}{1+2\left(1-s_{k}^{3}\right)^{\frac {1}{3}}}}\\s_{k+1}&={\frac {r_{k+1}-1}{2}}\\a_{k+1}&=r_{k+1}^{2}a_{k}-3^{k}\left(r_{k+1}^{2}-1\right)\end{aligned}}}

Entonces , k converge cúbicamente a 1 / π ; es decir, cada iteración triplica aproximadamente el número de dígitos correctos .

Convergencia cuártica (1985)

Comience configurando [ 5 ]

a0=2(21)2y0=21{\displaystyle {\begin{aligned}a_{0}&=2\left({\sqrt {2}}-1\right)^{2}\\y_{0}&={\sqrt {2}}-1\end{aligned}}}

Luego itera

yk+1=1(1yk4)141+(1yk4)14ak+1=ak(1+yk+1)422k+3yk+1(1+yk+1+yk+12){\displaystyle {\begin{aligned}y_{k+1}&={\frac {1-\left(1-y_{k}^{4}\right)^{\frac {1}{4}}}{1+\left(1-y_{k}^{4}\right)^{\frac {1}{4}}}}\\a_{k+1}&=a_{k}\left(1+y_{k+1}\right)^{4}-2^{2k+3}y_{k+1}\left(1+y_{k+1}+y_{k+1}^{2}\right)\end{aligned}}}

Entonces , k converge cuárticamente con respecto a 1 / π ; es decir, cada iteración cuadruplica aproximadamente el número de dígitos correctos. El algoritmo no es autocorrectivo; cada iteración debe realizarse con el número deseado de dígitos correctos para el resultado final de π .

Una iteración de este algoritmo es equivalente a dos iteraciones del algoritmo de Gauss-Legendre . Una demostración de estos algoritmos se puede encontrar aquí: [ 6 ].

Convergencia quíntica

Comience por configurar

a0=12s0=5(52)=5ϕ3{\displaystyle {\begin{aligned}a_{0}&={\frac {1}{2}}\\s_{0}&=5\left({\sqrt {5}}-2\right)={\frac {5}{\phi ^{3}}}\end{aligned}}}

dóndeϕ=1+52{\displaystyle \phi ={\tfrac {1+{\sqrt {5}}}{2}}}es la proporción áurea . Luego itera

incógnitanorte+1=5snorte1ynorte+1=(incógnitanorte+11)2+7znorte+1=(12incógnitanorte+1(ynorte+1+ynorte+124incógnitanorte+13))15anorte+1=snorte2anorte5norte(snorte252+snorte(snorte22snorte+5))snorte+1=25(znorte+1+incógnitanorte+1znorte+1+1)2snorte{\displaystyle {\begin{aligned}x_{n+1}&={\frac {5}{s_{n}}}-1\\y_{n+1}&=\left(x_{n+1}-1\right)^{2}+7\\z_{n+1}&=\left({\frac {1}{2}}x_{n+1}\left(y_{n+1}+{\sqrt {y_{n+1}^{2}-4x_{n+1}^{3}}}\right)\right)^{\frac {1}{5}}\\a_{n+1}&=s_{n}^{2}a_{n}-5^{n}\left({\frac {s_{n}^{2}-5}{2}}+{\sqrt {s_{n}\left(s_{n}^{2}-2s_{n}+5\right)}}\right)\\s_{n+1}&={\frac {25}{\left(z_{n+1}+{\frac {x_{n+1}}{z_{n+1}}}+1\right)^{2}s_{n}}}\end{aligned}}}

Entonces, k converge quínticamente a 1 / π ( es decir , cada iteración quintuplica aproximadamente el número de dígitos correctos), y se cumple la siguiente condición:

0<anorte1π<165nortemi5norteπ{\displaystyle 0<a_{n}-{\frac {1}{\pi }}<16\cdot 5^{n}\cdot e^{-5^{n}}\pi \,\!}

Convergencia nónica

Comience por configurar

a0=13r0=312s0=(1r03)13{\displaystyle {\begin{aligned}a_{0}&={\frac {1}{3}}\\r_{0}&={\frac {{\sqrt {3}}-1}{2}}\\s_{0}&=\left(1-r_{0}^{3}\right)^{\frac {1}{3}}\end{aligned}}}

Luego itera

tnorte+1=1+2rnortenorte+1=(9rnorte(1+rnorte+rnorte2))13vnorte+1=tnorte+12+tnorte+1norte+1+norte+12wnorte+1=27(1+snorte+snorte2)vnorte+1anorte+1=wnorte+1anorte+32norte1(1wnorte+1)snorte+1=(1rnorte)3(tnorte+1+2norte+1)vnorte+1rnorte+1=(1snorte+13)13{\displaystyle {\begin{aligned}t_{n+1}&=1+2r_{n}\\u_{n+1}&=\left(9r_{n}\left(1+r_{n}+r_{n}^{2}\right)\right)^{\frac {1}{3}}\\v_{n+1}&=t_{n+1}^{2}+t_{n+1}u_{n+1}+u_{n+1}^{2}\\w_{n+1}&={\frac {27\left(1+s_{n}+s_{n}^{2}\right)}{v_{n+1}}}\\a_{n+1}&=w_{n+1}a_{n}+3^{2n-1}\left(1-w_{n+1}\right)\\s_{n+1}&={\frac {\left(1-r_{n}\right)^{3}}{\left(t_{n+1}+2u_{n+1}\right)v_{n+1}}}\\r_{n+1}&=\left(1-s_{n+1}^{3}\right)^{\frac {1}{3}}\end{aligned}}}

Entonces , k converge nonológicamente a 1 / π ; es decir, cada iteración multiplica aproximadamente por nueve el número de dígitos correctos . [ 7 ]

Véase también

Referencias

  1. Jonathan M. Borwein, Peter B. Borwein, Pi and the AGM – A Study in Analytic Number Theory and Computational Complexity , Wiley, Nueva York, 1987. Muchos de sus resultados están disponibles en: Jorg Arndt, Christoph Haenel, Pi Unleashed, Springer, Berlín, 2001, ISBN 3-540-66572-2
  2. Bailey, David H (1 de abril de 2023). "Peter Borwein: Un matemático visionario". Notices of the American Mathematical Society . 70 (4): 610– 613. doi : 10.1090/noti2675 . ISSN 0002-9920 . 
  3. Borwein, JM; Borwein, PB (1993). "Serie tipo Ramanujan de clase tres para 1/π" . Journal of Computational and Applied Mathematics . 46 ( 1–2 ): 281–290 . doi : 10.1016/0377-0427(93)90302-R .
  4. Arndt, Jörg; Haenel, Christoph (1998). π desatado . Springer-Verlag. pag. 236.ISBN  3-540-66572-2.
  5. Mak, Ronald (2003). The Java Programmers Guide to Numerical Computation . Pearson Educational. p. 353. ISBN  0-13-046041-9.
  6. Milla, Lorenz (2019), Prueba sencilla de tres algoritmos π recursivos , arXiv : 1907.04110
  7. Henrik Vestermark (4 de noviembre de 2016). "Implementación práctica de Algoritmos π" (PDF) . Consultado el 29 de noviembre de 2020 .
  • Fórmulas de Pi de Wolfram MathWorld
Obtenido de " https://en.wikipedia.org/w/index.php?title=Borwein%27s_algorithm&oldid=1321323705 "