Rekursiv formel för beräkning av pi

Kanske får det vara mitt amatörbidrag till matematiken. Troligare finns formeln redan, eller andra, smartare och elegantare varianter. Var så god och sök på webben 🙂

Formeln bygger på iakttagelsen att en regelbunden månghörning som är inskriven i en cirkel med radien 1, har en omkrets som närmar sig pi när antalet hörn blir större.

Idén är att börja med en kvadrat inskriven i cirkeln. Den har omkretsen $latex 2\sqrt{2}&s=1$ . Dubbla sedan antalet hörn till 8. Dess sida kan relativt enkelt beräknas med hjälp av kvadraten. Nästa steg är en 16-hörning, vars sida kan beräknas med åttahörningen. Etc.

Alltså:
Antag att vi har en cirkel med radie 1.
I cirkeln skriver vi in en regelbunden månghörning med $latex 2^{n}&s=1$ hörn, där $latex n>1&s=1$.
Antag vidare att vi känner till hur lång sidan är på denna månghörning.
Kan vi då beräkna sidan på månghörningen med $latex 2^{n+1}&s=1$ hörn?

Låt oss dela upp den första månghörningen (den med $latex 2^{n}&s=1$ hörn) i likbenta trianglar med toppen i cirkelns mitt. Varje triangel ser då ut något i stil med nedanstående figur. Dela även upp den likbenta triangeln i två rätvinkliga trianglar:

tri6

Vi antar som sagt att för rektangeln med $latex 2^{n}&s=1$ hörn så har vi räknat ut sidan $latex C&s=1$. Omkretsen på denna månghörning kommer då att vara $latex {C_{n}}*2^{n}&s=1$ Omkretsen av cirkeln med radien 1 är enligt skolboken $latex 2 * \pi &s=1$. När antalet hörn ökar kommer alltså $latex {C_{n}}*2^{n – 1}&s=1$ att bli en bättre och bättre approximation av pi.

Låt nästa steg vara en månghörning med dubbelt så många hörn, $latex 2^{n+1}$. Dela även upp denna månghörning i likbenta trianglar med topp i cirkelns mitt, och dela upp den likbenta triangeln i två rätvinkliga trianglar:

tri8

Observera nu att de blåmarkerade trianglarna är kongruenta. Med den iakttagelsen går det att hitta ett samband mellan $latex {C_{n}}&s=1$ och $latex {C_{n+1}}&s=1$. Med lite räknande får man sambandet:

$latex {C_{n+1}} = \sqrt{2 – \sqrt{4-C_{n}^2}}&s=1$

Enkelt uttryckt innebär det att vi har ett samband mellan sidan på en månghörning och sidan på en månghörning med dubbet så många hörn.

Givet att en regelbunden månghörning inskriven i en cirkel med radien 1 har sidan $latex {C_{n}}&s=1$, så vet vi nu att månghörningen med dubbelt så många hörn har omkretsen:

$latex 2^{{n}+1}{C_{n+1}} =2^{{n}+1}\sqrt{2 – \sqrt{4-C_{n}^2}}&s=1$

Efter lite mer räkningar kan vi nu ta fram en rekursiv formel för successiv beräkning av pi:

$latex \pi_{{n}+1}=2^{({n}+1)/2}\sqrt{2^{n}-\sqrt{4^{n}-\pi_{n}^2}}&s=1$

Startvärdet är en kvadrat inskriven i cirkeln, dvs $latex {n}=2&s=1$:

$latex \pi_{1}=2\sqrt{2}&s=1$

Utvecklingen går någorlunda snabbt. Approximation för pi blir:

4-hörning: $latex \pi\approx 2,8&s=1$
8-hörning: $latex \pi\approx 3,06&s=1$
16-hörning: $latex \pi\approx 3,12&s=1$
32-hörning: $latex \pi\approx 3,137&s=1$
64-hörning: $latex \pi\approx 3,140&s=1$
128-hörning: $latex \pi\approx 3,1413&s=1$
256-hörning: $latex \pi\approx 3,14151&s=1$
512-hörning: $latex \pi\approx 3,14157&s=1$
1024-hörning: $latex \pi\approx 3,141588&s=1$
2048-hörning: $latex \pi\approx 3,141591&s=1$
4096-hörning: $latex \pi\approx 3,1415923&s=1$
8192-hörning: $latex \pi\approx 3,1415925&s=1$

/David Armini

 

3 reaktioner till “Rekursiv formel för beräkning av pi”

  1. Som kuriosa borde formeln kunna användas med andra startvärden, tex börja med en triangel, vidare till sexhörning, 12-hörning. Etc.

  2. Som mer kuriosa går det fint att beräkna i Excel. Formeln att använda är =2^((D6)/2)*SQRT(2^(D6-1)-SQRT(4^(D6-1)-E5^2))

    D E
    2 2,828427124746190000000000000000
    3 3,061467458920720000000000000000
    4 3,121445152258050000000000000000
    5 3,136548490545940000000000000000
    6 3,140331156954740000000000000000
    7 3,141277250932760000000000000000
    8 3,141513801144150000000000000000
    9 3,141572940367880000000000000000
    10 3,141587725279960000000000000000
    11 3,141591421504640000000000000000
    12 3,141592345611080000000000000000
    13 3,141592576545000000000000000000
    14 3,141592633463250000000000000000
    15 3,141592654807590000000000000000
    16 3,141592645321220000000000000000
    17 3,141592607375720000000000000000
    18 3,141592910939670000000000000000
    19 3,141591696683680000000000000000
    20 3,141596553704820000000000000000
    Som synes stannar konvergeringen vid 13 – det beror på att Excel ”bara” räknar med 15 signifikanta siffror. Precisionen räcker alltså inte för att mer exakt beräkning av pi. Notera också den lilla Microsoftigheten att Excel räknar med 15 siffror men kan visa upp till 30 decimaler. Why, oh why…?

Lämna ett svar

Din e-postadress kommer inte publiceras. Obligatoriska fält är märkta *