Expected perimeter of the convex hull of planar Brownian motion stopped upon exiting the unit disk

Deze paper berekent de exacte verwachte omtrek van de convexe omhulling van een planaire Brownse beweging die stopt bij het verlaten van de eenheidsschijf, door het probleem te reduceren tot het bepalen van de verwachte maximale horizontale verplaatsing en deze uit te drukken in termen van harmonische maat, terwijl ook niet-triviale grenzen worden afgeleid voor de verwachte oppervlakte.

De kern

Stel je voor dat je een dronken wandelaar bent die in het midden van een perfect rond plein (een eenheidsschijf) begint te lopen. Je loopt willekeurig, links, rechts, vooruit, achteruit, zonder plan. Je blijft lopen tot je per ongeluk de rand van het plein raakt. Op dat moment stopt de wandeling.

Nu komt de vraag: als je al die willekeurige stappen op een kaart tekent en er een elastiekje omheen spant dat precies om al je stappen past, hoe groot is de omtrek van die vorm? En hoe groot is het oppervlak erin?

Dit is precies wat de auteurs van dit artikel, Hugo Panzo en Stjepan Šebek, hebben onderzocht. Ze kijken naar de wiskundige "dronken wandeling" (Browse-beweging) in een cirkel en proberen de gemiddelde grootte van de vorm te vinden die deze wandeling omsluit.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen:

1. Het Probleem: De "Dronken" Omtrek

In de wiskunde noemen ze de vorm die ontstaat door een elastiekje om een verzameling punten te spannen, de convexe hull (het omhulsel).

• De uitdaging: Meestal kijken wiskundigen naar wandelingen die een vast aantal stappen doen. Maar hier kijken ze naar een wandeling die stopt op een willekeurig moment: het moment dat de wandelaar de muur raakt.
• Het doel: Ze willen weten wat de gemiddelde omtrek is van die vorm.

2. De Oplossing: De Omtrek is Eenvoudiger dan Je denkt

Het verrassende nieuws is dat ze een exact antwoord hebben gevonden voor de omtrek!

De Metafoor van de Schaduw:
Stel je voor dat je de wandeling bekijkt vanuit verschillende hoeken. Als je een lichtbron heel ver weg zet en de wandeling projecteert op een muur, zie je een schaduw. De lengte van die schaduw is de "steunfunctie".
De auteurs gebruiken een slimme wiskundige truc (Cauchy's formule) die zegt: De totale omtrek van de vorm is gewoon de som van al die schaduwen die je krijgt als je om de vorm heen loopt.

Omdat de wandeling en het plein perfect rond en symmetrisch zijn, is het gemiddelde van die schaduwen in elke richting hetzelfde. Het probleem reduceert zich dus tot één simpele vraag:
"Hoe ver komt de wandelaar gemiddeld naar rechts (de horizontale as) voordat hij de muur raakt?"

Het Resultaat:
Ze hebben een formule gevonden die precies beschrijft hoe de kans is dat de wandelaar een bepaalde afstand naar rechts heeft afgelegd.

• De gemiddelde afstand naar rechts is ongeveer 0,51.
• Als je dit vermenigvuldigt met de cirkel (2π), krijg je de gemiddelde omtrek van de vorm.
• Het antwoord: De gemiddelde omtrek is ongeveer 3,21. (Ter vergelijking: de omtrek van het plein zelf is ongeveer 6,28. De vorm is dus ongeveer de helft zo groot als de cirkel, wat logisch is omdat de wandelaar niet de hele cirkel vult).

3. De Moeilijke Deel: Het Oppervlak

Nu wordt het lastig. Ze hebben ook gekeken naar het oppervlak (de ruimte binnenin het elastiekje).

• De vergelijking: Als de omtrek berekenen als het "rekenen van de randen van een hek" is, dan is het oppervlak berekenen als het "meten van het gras erin".
• Het probleem: Om het oppervlak exact te berekenen, moet je niet alleen weten hoe ver de wandelaar is gekomen, maar ook wanneer hij daar was en hoe hij daar is gekomen. Het tijdstip waarop de wandelaar zijn verste punt bereikt, is een heel lastig, willekeurig moment dat koppelt aan zijn positie.
• De conclusie: Ze hebben geen exacte formule gevonden voor het oppervlak. Het is te ingewikkeld, als een puzzel waarbij de stukjes niet goed in elkaar passen.

Wat ze wel deden:
Ze hebben een "ondergrens" gevonden. Ze bedachten een kleinere vorm, de "ster-hull" (stel je voor dat je een ster tekent door lijnen te trekken van het midden naar alle wandelpunten). Deze ster zit altijd binnen het elastiekje.

• Ze konden de gemiddelde grootte van die ster wel precies uitrekenen (ongeveer 0,47).
• Dus weten ze zeker: het oppervlak van het elastiekje is minstens 0,47.
• Via computersimulaties (waar ze 100.000 keer een virtuele wandelaar lieten lopen) schatten ze dat het echte oppervlak waarschijnlijk rond de 0,66 ligt.

Samenvatting in het Kort

1. Omtrek: Ze hebben een perfecte formule gevonden. De gemiddelde omtrek van de vorm die een willekeurige wandeling in een cirkel omsluit, is 3,21.
2. Oppervlak: Dit is veel lastiger. Ze hebben geen exacte formule, maar weten wel dat het ergens tussen de 0,47 en 1,14 ligt (waarschijnlijk rond de 0,66).
3. Waarom is dit cool? Het laat zien dat zelfs bij volledig willekeurige bewegingen, er soms mooie, vaste patronen en getallen te vinden zijn, zolang je maar de juiste wiskundige brillen opzet (zoals het gebruik van "harmonische maten" en "conforme afbeeldingen", wat in het kort betekent: het vervormen van de ruimte om het probleem makkelijker te maken).

Kortom: De omtrek is een geslaagde wiskundige ontdekking, maar het oppervlak blijft een mysterie dat nog even op zijn eigen manier moet worden ontrafeld.

Technische samenvatting

Titel en Context

Het artikel onderzoekt de convexe omhulling (convex hull) van een standaard planaire Brownse beweging $W_t$ die start in de oorsprong en stopt op het moment dat het de eenheidsschijf $D$ verlaat (de eerste uitstoptijd $\tau_D$). De hoofdfocus ligt op het berekenen van de verwachte omtrek van deze convexe omhulling. Dit vult recent onderzoek aan over de statistiek van Brownse bewegingen in beklemde geometrieën, specifiek in contrast met reflecterende Brownse bewegingen (Neumann-randvoorwaarden) die eerder zijn bestudeerd. Hier wordt gewerkt met een "gedood" proces (Dirichlet-randvoorwaarden) waarbij de beweging stopt bij het raken van de rand.

Het Probleem

Laat $H_{\tau_D}$ de convexe omhulling zijn van de trajectorie van de Brownse beweging $W$ op het interval $[0, \tau_D]$. De auteurs willen de verwachte waarde van de omtrek $P_{\tau_D} = \text{Per}(H_{\tau_D})$ berekenen.
Hoewel de verwachte omtrek voor een Brownse beweging op een vast tijdstip $t=1$ bekend is ($\mathbb{E}[P_1] = \sqrt{8\pi}$), was er geen gesloten vorm bekend voor het geval dat de beweging stopt bij een willekeurige uitstoptijd uit een domein.

Methodologie

De aanpak combineert drie krachtige wiskundige concepten:

1. Cauchy's oppervlakteformule: Deze stelt een link tussen de omtrek van een convex lichaam en de integraal van zijn ondersteuningsfunctie (support function).
2. Conforme invariantie: De eigenschap dat de wet van de Brownse beweging behouden blijft onder conforme afbeeldingen.
3. Harmonische maat en de Poisson-kern: Het gebruik van de Poisson-kern in het bovenste halfvlak om uitstoppingskansen te berekenen.

Stap 1: Reductie tot maximale horizontale verplaatsing
Met behulp van Cauchy's formule wordt de omtrek $P_{\tau_D}$ uitgedrukt als de integraal van de ondersteuningsfunctie $h(\theta)$ over alle richtingen. Door rotatie-invariantie van de Brownse beweging en de schijf, is de verwachte waarde van de ondersteuningsfunctie in elke richting gelijk aan de verwachte waarde van de maximale horizontale verplaatsing $M = \sup_{0 \le t \le \tau_D} X_t$.
Hieruit volgt direct:
$$\mathbb{E}[P_{\tau_D}] = 2\pi \mathbb{E}[M]$$

Stap 2: Berekening van de verdeling van $M$
Om $\mathbb{E}[M]$ te vinden, analyseren de auteurs de cumulatieve verdelingsfunctie (CDF) van $M$.

• Ze definiëren een "afgeknotte schijf" $D_a = D \cap \{z : \text{Re}(z) < a\}$.
• De gebeurtenis $\{M \ge a\}$ komt overeen met het feit dat de Brownse beweging $D_a$ verlaat via de verticale koorde $V_a$ in plaats van via de cirkelboog $C_a$.
• Ze construeren een conforme afbeelding $f_a$ die het domein $D_a$ afbeeldt op het bovenste halfvlak $\mathbb{H}$. Dit gebeurt in twee stappen:
  1. Een lineaire fractionele transformatie die $D_a$ afbeeldt op een wig (wedge) met een specifieke opening.
  2. Een machtsfunctie die de wig opent tot het bovenste halfvlak.
• Door conforme invariantie wordt het probleem gereduceerd tot het berekenen van de harmonische maat van een straal in het bovenste halfvlak gezien vanuit het beeld van de oorsprong. Dit leidt tot een exacte formule voor de CDF van $M$.

Stap 3: Berekening van de verwachte waarde
De verwachte waarde $\mathbb{E}[M]$ wordt berekend door de integraal van de overlevingsfunctie $P(M \ge a)$ over $a \in [0,1]$. Dit resulteert in een integraal die numeriek kan worden opgelost.

Belangrijkste Resultaten

1. Exacte uitdrukking voor de omtrek
De auteurs leiden een exacte formule af voor de verwachte omtrek:
$$\mathbb{E}[P_{\tau_D}] = 2\pi \int_0^{\pi/2} \left( 1 - \frac{2t}{\frac{\pi}{2} + t} \right) \cos t \, dt \approx 3.214826$$
De cumulatieve verdelingsfunctie van de maximale horizontale verplaatsing $M$ is:
$$P(M < a) = \frac{2 \arcsin a}{\pi - \arccos a}, \quad 0 < a < 1$$
De verwachte waarde van $M$ zelf is $\mathbb{E}[M] \approx 0.511655$.

2. Grenzen voor de verwachte oppervlakte
Voor de oppervlakte $A_{\tau_D}$ van de convexe omhulling vinden de auteurs geen gesloten vorm. De complexiteit ontstaat doordat de oppervlakte afhankelijk is van zowel de maximale verplaatsing als het tijdstip waarop dit maximum wordt bereikt, wat leidt tot een complexe koppeling tussen de coördinaten.
Ze bieden echter strenge grenzen:

• Bovenste grens: $\mathbb{E}[A_{\tau_D}] \le \pi \mathbb{E}[M^2] \approx 1.1397$.
• Onderste grens: Ze gebruiken de "ster-omhulling" (star hull) van het pad. Omdat de ster-omhulling altijd binnen de convexe omhulling ligt, geldt $\mathbb{E}[A_{\tau_D}] \ge \mathbb{E}[\text{Area}(\text{star hull})]$.
  De verwachte oppervlakte van de ster-omhulling wordt exact berekend als:
  $$\mathbb{E}[\text{Area}(\text{star hull})] = \pi - \frac{8}{3} \approx 0.474925$$
  Dus: $0.474925 \le \mathbb{E}[A_{\tau_D}] \le 1.1397$.

3. Simulaties
Monte Carlo-simulaties met $10^5$ trajectorieën bevestigen de theoretische resultaten voor de omtrek (geschatte waarde $3.2136$ vs. theoretisch $3.2148$). Voor de oppervlakte geven de simulaties een schatting van ongeveer $0.6612$, wat binnen de afgeleide grenzen valt.

Significantie en Bijdrage

• Nieuw inzicht: Dit is een van de eerste werken dat een exacte oplossing biedt voor geometrische functionalen van een Brownse beweging die stopt bij een willekeurige uitstoptijd (Dirichlet-voorwaarde), in tegenstelling tot vaste tijdstippen of reflecterende randen.
• Methodologische innovatie: De combinatie van conforme afbeeldingen met specifieke domeinen (afgeknotte schijven) en harmonische maat biedt een krachtig kader voor het analyseren van extremen van Brownse bewegingen in begrensde gebieden.
• Verband met bestaande theorie: Het werk bouwt voort op klassieke resultaten van Letac en Takács (voor vaste tijd) en recent werk over reflecterende bewegingen, maar toont aan dat de "Dirichlet-omhulling" een unieke statistische structuur heeft.
• Open probleem: Het artikel benadrukt dat het vinden van een exacte formule voor de oppervlakte aanzienlijk moeilijker is dan voor de omtrek, voornamelijk door de afhankelijkheid van de verdeling van het tijdstip van het maximum, wat de symmetrie-eigenschappen van de klassieke arcsine-wet doorbreekt.

Samenvattend levert dit artikel een elegante, exacte oplossing voor de verwachte omtrek van een Brownse convexe omhulling in een eenheidsschijf en biedt het robuuste grenzen voor de oppervlakte, waarbij het een brug slaat tussen stochastische processen en meetkunde.