Estimating $\pi$ with a Coin

Dit artikel beschrijft een nieuwe Monte Carlo-methode om $\pi$ te schatten door middel van het gooien van een munt, waarbij een nieuwe interpretatie van $\frac{\pi}{4}$ wordt gegeven die voortkomt uit bekende Catalan-getal-identiteiten.

De kern

Stel je voor dat je een muntstuk hebt en je wilt het getal π (pi, ongeveer 3,14) berekenen. Normaal gesproken denk je dan aan cirkels, diameters en meetkunde. Maar in dit leuke artikel beschrijft Jim Propp een manier om π te vinden met alleen maar een muntstuk en een beetje geduld.

Het is alsof je een wiskundig avontuur start waarbij je de munt gooit tot je een specifiek patroon ziet, en dan een simpele rekensom maakt.

Hier is hoe het werkt, stap voor stap:

1. Het Spel: "Wacht tot de koppen winnen"

Je begint met het gooien van een eerlijke munt (50% kans op kop, 50% op munt).

• Je telt hoeveel keer Kop (hoofden) en Munt (staarten) er vallen.
• Je blijft gooien tot het moment waarop het aantal Koppen voor het eerst hoger is dan het aantal Munten.
• Op dat exacte moment stop je en noteer je het totaal aantal worpen.

Voorbeeld:
Stel je gooit: Munt, Kop, Munt, Kop, Kop.

• Na 1 worp: 0 Kop, 1 Munt (Munt wint).
• Na 2 worpen: 1 Kop, 1 Munt (Gelijkspel).
• Na 3 worpen: 1 Kop, 2 Munten (Munt wint nog steeds).
• Na 4 worpen: 2 Kop, 2 Munten (Gelijkspel).
• Na 5 worpen: 3 Kop, 2 Munten (Kop wint!).

Je stopt bij worp 5. Je hebt toen 3 koppen en 2 munten. Je rekent nu de breuk uit: 3 gedeeld door 5 (het aantal koppen gedeeld door het totaal). Dat is 0,6.

2. De Magische Herhaling

Je doet dit spelletje niet één keer, maar heel vaak.

• Je begint telkens opnieuw (alsof je de eerdere worpen vergeet).
• Elke keer stop je op het moment dat Koppen voor het eerst winnen.
• Je noteert elke keer die breuk (bijvoorbeeld 0,6, dan 0,55, dan 0,66, etc.).

Als je dit duizenden keren doet en al die breuken bij elkaar optelt en deelt door het aantal keer dat je hebt gespeeld (het gemiddelde), dan gebeurt er iets wonderlijks:
Het gemiddelde komt steeds dichter bij 0,78539...

En dat getal is precies π gedeeld door 4 ($\frac{\pi}{4}$).

3. Waarom werkt dit? (De "Wiskundige Toer" in het kort)

Het artikel legt uit dat dit te maken heeft met een soort "wiskundige dans" die we een random walk noemen.

• Stel je voor dat Kop een stapje naar boven is en Munt een stapje naar beneden.
• Je begint op de grond (0).
• Je wilt weten hoe lang het duurt voordat je voor het eerst op de eerste verdieping (+1) bent.
• De wiskunde achter deze dans (die te maken heeft met getallenreeksen die we de Catalan-getallen noemen) zorgt ervoor dat de gemiddelde verhouding van "stappen naar boven" tot "totale stappen" precies uitkomt op $\frac{\pi}{4}$.

Het is alsof de natuur zelf een verborgen code heeft in het gooien van munten die de cirkel (π) onthult, zolang je maar wacht tot de "Koppen" de overhand hebben.

4. Is dit een snelle manier om π te vinden?

Nou, niet echt. Het artikel geeft een grappig voorbeeld:

• Als je 10.000 worpen doet, kom je uit op ongeveer 3,22. Dat is niet heel nauwkeurig.
• Om echt dicht bij 3,14 te komen, heb je misschien wel een biljoen worpen nodig.
• Als je één munt per seconde gooit, zou dat 30.000 jaar duren!

Dus, hoewel het een prachtige wiskundige ontdekking is die laat zien hoe diep de verbinding tussen kansrekening en cirkels zit, is het geen praktische manier om je rekenmachine te vervangen.

De Kernboodschap

Het artikel is eigenlijk een mooie ontdekkingstocht. Het laat zien dat als je een simpel spelletje speelt (wacht tot Koppen winnen) en de uitkomsten gemiddeld, je onverwachts het getal π tegenkomt. Het is een bewijs dat wiskunde overal om ons heen zit, zelfs in het simpele geklapper van een muntstuk.

Kort samengevat in één zin:
Als je oneindig vaak een munt gooit en elke keer stopt zodra je net iets meer koppen dan munten hebt, dan is het gemiddelde percentage koppen op dat moment precies π/4.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Estimating π with a Coin" van Jim Propp, geschreven in het Nederlands.

Titel: Schatting van π met een Munt

Auteur: Jim Propp
Datum: November 2024 / Herziening maart 2026

---

1. Probleemstelling

Het artikel introduceert een nieuwe, eenvoudige Monte Carlo-methode om de wiskundige constante $\pi$ te schatten door middel van het gooien van een eerlijke munt. Hoewel er reeds bekende methoden bestaan, zoals de naald van Buffon (waarbij de kans dat een naald een lijn raakt gelijk is aan $2/\pi$), presenteert deze paper een alternatief dat gebaseerd is op het analyseren van de verhouding tussen kop en munt op het moment dat een specifieke stopconditie wordt bereikt. Het doel is om te bewijzen dat de verwachte waarde van deze verhouding exact gelijk is aan$\pi/4$.

2. Methodologie

De methode maakt gebruik van de theorie van willekeurige wandelingen (random walks) en stoptijden (stopping times).

• Het Experiment:

  1. Gooi een eerlijke munt herhaaldelijk en tel het aantal koppen ($H_n$) en munts ($T_n$) na $n$ worpen.
  2. Bepaal het aantal koppen en munts tot het moment dat het aantal koppen voor het eerst strikt groter is dan het aantal munts.
  3. Laat $\tau$ de willekeurige stoptijd zijn: $\tau = \min\{n \ge 1 : H_n > T_n\}$.
  4. Noteer de fractie van koppen op dat moment: $H_\tau / \tau$.
  5. Herhaal dit proces over meerdere runs en bereken het gemiddelde van deze fracties.

• Wiskundige Formulering:
  De situatie wordt gemodelleerd als een simpele symmetrische willekeurige wandeling $S_n = H_n - T_n$ op de gehele getallen $\mathbb{Z}$, beginnend bij $S_0 = 0$. De stapgrootte is $\pm 1$.
  De stoptijd $\tau$ is het eerste moment dat de wandeling $+1$ bereikt. Omdat de wandeling in stappen van 1 beweegt, moet $\tau$ oneven zijn ($\tau = 2k + 1$).

• Katalan-getallen:
  Om op tijd $2k+1$voor het eerst$+1$te bereiken, moet de wandeling de eerste$2k$stappen$\le 0$blijven, op tijd$2k$precies op$0$eindigen, en vervolgens een stap naar$+1$maken. Het aantal van dergelijke paden wordt gegeven door het$k$-de Katalan-getal:
  $$C_k = \frac{1}{k+1}\binom{2k}{k}$$
  De waarschijnlijkheid dat $\tau = 2k + 1$ is:
  $$P(\tau = 2k + 1) = \frac{1}{2} \cdot \frac{1}{4^k} \cdot \frac{1}{k+1} \binom{2k}{k}$$

3. Belangrijkste Bijdragen en Afleiding

De kernbijdrage van het artikel is het bewijs dat de verwachte waarde van de verhouding $H_\tau/\tau$ gelijk is aan $\pi/4$.

• Berekening van de Verwachte Waarde:
  Op het moment van stoppen ($\tau = 2k+1$) geldt dat $H_\tau = k+1$ en $T_\tau = k$. De verhouding is dus $\frac{k+1}{2k+1}$.
  De verwachte waarde wordt berekend als:
  $$E\left[\frac{H_\tau}{\tau}\right] = \sum_{k \ge 0} P(\tau = 2k + 1) \cdot \frac{k+1}{2k+1}$$
  Na invullen van de waarschijnlijkheidsformule en vereenvoudiging:
  $$E\left[\frac{H_\tau}{\tau}\right] = \frac{1}{2} \sum_{k \ge 0} \frac{1}{4^k} \frac{1}{2k+1} \binom{2k}{k}$$

• Verbinding met de Arcsinus:
  De auteur identificeert de bovenstaande reeks als de machtreeksontwikkeling van de arcsinusfunctie voor $x=1$:
  $$\arcsin(x) = \sum_{n=0}^{\infty} \frac{1}{2^{2n}} \binom{2n}{n} \frac{x^{2n+1}}{2n+1}$$
  Door $x=1$ in te vullen, volgt $\arcsin(1) = \pi/2$.
  Hieruit volgt direct:
  $$E\left[\frac{H_\tau}{\tau}\right] = \frac{1}{2} \cdot \frac{\pi}{2} = \frac{\pi}{4}$$

4. Resultaten en Observaties

• Nieuwheid: Hoewel de verwachte waarde van $1/\tau$(gelijk aan$\pi/2 - 1$) al bekend was in de literatuur over Katalan-getallen, is de interpretatie van$\pi/4$als de verwachte stoptijd-waarde van de verhouding$H_\tau/\tau$ nieuw.
• Schatting van $\pi$: Door het gemiddelde van de waargenomen fracties over vele runs te nemen, convergeert deze naar $\pi/4$. Vermenigvuldiging met 4 geeft een schatting van $\pi$.
• Convergentiesnelheid: De nauwkeurigheid van deze methode is laag. De fout neemt af met de snelheid $1/N^{1/4}$, waarbij$N$ het totale aantal muntenworpen is.
  • Voorbeeld: Met 10.000 worpen (zoals gedaan door Matt Parker) werd een waarde van 3,2266 verkregen (fout $\approx 0,1$).
  • Om een nauwkeurige waarde van 3,14 te krijgen, zouden naar schatting een biljoen worpen nodig zijn, wat bij één worp per seconde meer dan 30.000 jaar zou duren.
• Variatie op het thema:
  • Als de stopconditie wordt gewijzigd zodat koppen munts voor het eerst met +2 overtreffen, is de verwachte verhouding gelijk aan $\ln 2$.
  • Voor een algemene overwinning met $+a$ lijkt de verwachte waarde gerelateerd aan $\pi$ als $a$ oneven is, en aan $\ln 2$ als $a$ even is (hoewel dit nog als conjectureel wordt beschouwd).

5. Betekenis en Conclusie

• Probabilistisch Bewijs: Het artikel biedt een elegant probabilistisch bewijs voor de bekende ongelijkheid $3 < \pi < 4$. Aangezien de verhouding$H_\tau/\tau$altijd tussen$1/2$en$1$ligt (en gemiddeld groter is dan$3/4$), impliceert de gelijkheid$E = \pi/4$dat$\pi$ strikt tussen 3 en 4 moet liggen.
• Didactische Waarde: De methode is eenvoudig uit te voeren in een klaslokaal of wiskundeclub, hoewel de statistische onnauwkeurigheid betekent dat het vooral nuttig is als conceptueel voorbeeld en niet als praktische rekenmethode.
• Vergelijking: De auteur merkt op dat andere algoritmen voor het schatten van $\pi$ met munten (zoals die van Flajolet, Pelletier en Soria) waarschijnlijk veel efficiënter zijn.

Samenvattend presenteert dit artikel een elegante verbinding tussen de theorie van willekeurige wandelingen, Katalan-getallen en de constante $\pi$, hoewel de praktische toepassing voor nauwkeurige berekeningen beperkt blijft door de trage convergentie.