Fala 108 — 10 nowych symulacji w sztuce generatywnej, prawdopodobieństwie, hydrodynamice i algorytmach

Jest coś eleganckiego w ćwierćkolistym łuku narysowanym na kwadratowym kafelku. Sam w sobie nic nie znaczy. Obróć go losowo i umieść obok kolejnego — i jeszcze jednego, i jeszcze — a pojawią się emergentne wzory na dużą skalę, których nie przewidywał żaden pojedynczy kafelek. To kafelkowanie Truchet: licząca dwa stulecia matematyczna zabawka, która, jak się okazuje, koduje głębokie idee dotyczące losowości, spójności i perkolacji. To także pierwsza z dziesięciu symulacji Fali 108 — wydania, które obejmuje zakres od dziewiętnastowiecznej teorii liczb, przez dwudziestowieczne kody z korekcją błędów, po klasyczną sztukę islamską.

Sztuka generatywna: kafelki, spirale i wielokąty gwiaździste

Kafelki Truchet

Ojciec Sébastien Truchet opisał swój system kafelkowania w 1704 roku: weź kwadratowy kafelek z przekątną linią biegnącą od jednego rogu do przeciwległego, obróć go w jedną z czterech możliwych orientacji i układaj losowo. Powstałe wzory są zaskakująco bogate — obszary połączone łukami tworzą klastry, których rozkład rozmiarów jest powiązany z teorią perkolacji, a granica między połączonymi domenami jest losową krzywą o znanym wymiarze fraktalnym. Współczesny wariant z ćwierćkolistymi łukami (dwa łuki na kafelek zamiast przekątnej) spopularyzowali Cyril Smith i Pauline Boucher w 1987 roku i od tego czasu stał się on standardowym elementem praktyki sztuki generatywnej.

Symulacja pozwala regulować rozdzielczość siatki od 4 do 64 kafelków na bok, przełączać się między wariantem z przekątną a wariantem z łukami oraz ponownie losować układ jednym naciśnięciem klawisza. Nakładka pokazująca spójność na żywo podświetla największy spójny region kontrastującym kolorem, dzięki czemu struktura perkolacji jest natychmiast czytelna.

Cykliczne automaty komórkowe

Cykliczne automaty komórkowe (CCA) należą do najbardziej wizualnie porywających systemów w matematyce dyskretnej. Reguła jest zwodniczo prosta: komórka w stanie s przechodzi do stanu s+1 (mod k), jeśli przynajmniej jeden z jej sąsiadów jest już w stanie s+1. Na losowej siatce początkowej małe jądra uporządkowanych komórek rywalizują, łączą się, aż w końcu organizują całą siatkę w obracające się fale spiralne, które przemierzają siatkę niczym komórkowa zorza polarna.

Liczba stanów k, promień sąsiedztwa i próg (ile sąsiadów musi być w stanie s+1) są regulowalne. Niskie progi dają gęste spirale; wysokie progi dają rzadkie, wolno poruszające się fronty. Przejście między reżimami jest ostre i wykazuje cechy przypominające przejście fazowe w układzie statystyczno-mechanicznym.

Islamskie wzory geometryczne

Islamski ornament geometryczny rozwijał się przez ponad tysiąclecie, stając się jedną z najbardziej wyrafinowanych tradycji sztuki matematycznej w historii ludzkości. Leżąca u podstaw metoda konstrukcji wykorzystuje kafelki girih — zestaw pięciu kształtów wielokątów (dekagon, pentagon, wydłużony sześciokąt, „muszka" i romb), których kąty są wielokrotnościami 36° — do wypełniania płaszczyzny lokalną symetrią dziesięciokrotną i pięciokrotną. Wielokąty gwiaździste postaci {n/k} są wpisywane w kafelki, a ich linie przedłużane, tworząc splecione wzory wstęgowe widoczne w meczetach od Kordoby po Samarkandę.

Symulacja generuje rodziny wielokątów gwiaździstych sparametryzowane przez n (liczbę wierzchołków) i k (odstęp pomijania) i kafelkuje płótno przy użyciu odpowiedniej konstrukcji girih. Możesz wybrać dowolną grupę symetrii spośród 17 grup tapetowych, którymi posługiwali się islamscy rzemieślnicy, oraz wyeksportować wzór jako SVG do wykorzystania w pracach projektowych.

Prawdopodobieństwo: od wykrywania oszustw po finanse stochastyczne

Prawo Benforda

W 1938 roku Frank Benford zauważył, że wiodąca cyfra liczb pochodzących z wielu rzeczywistych zbiorów danych nie jest rozłożona jednostajnie. Cyfra 1 pojawia się jako pierwsza cyfra w mniej więcej 30% przypadków; cyfra 9 pojawia się tylko w około 4,6% przypadków. Rozkład jest logarytmiczny: P(d) = log10(1 + 1/d) dla d ∈ {1, ..., 9}. To nieintuicyjne prawo obowiązuje dla populacji miast, stałych fizycznych, wpisów w rejestrach finansowych, długości rzek i liczb Fibonacciego — zasadniczo dla dowolnego zbioru danych obejmującego kilka rzędów wielkości bez sztucznej granicy górnej.

Prawo Benforda jest dziś rutynowo wykorzystywane do wykrywania oszustw w zeznaniach podatkowych, wynikach wyborów i korporacyjnych zapisach księgowych: sfabrykowane liczby rzadziej zaczynają się od 1 niż dane prawdziwe. Symulacja pozwala wczytać plik CSV lub wybrać jeden z wbudowanych zbiorów danych (ciąg Fibonacciego, potęgi 2, populacje i losowe dane jednostajne) i nakłada obserwowany rozkład pierwszej cyfry na teoretyczną krzywą Benforda, wyświetlając wynik testu chi-kwadrat w czasie rzeczywistym.

Ruina gracza

Problem ruiny gracza to jeden z najstarszych wyników w teorii prawdopodobieństwa. Gracz zaczyna z k dolarami i stawia jednego dolara na rundę z prawdopodobieństwem p wygranej każdego zakładu, grając przeciwko kasynu dysponującemu n - k dolarami. Prawdopodobieństwo, że gracz osiągnie n dolarów, zanim zbankrutuje, wynosi:

P(ruin) = 1 - (1 - (q/p)^k) / (1 - (q/p)^n)   when p ≠ q
P(ruin) = 1 - k/n                               when p = q = 0.5

Symulacja uruchamia tysiące jednoczesnych błądzeń losowych, zestawiając empiryczne prawdopodobieństwo ruiny z krzywą teoretyczną w miarę regulowania k, n i p. Sama zbieżność symulacji z teorią jest pouczająca: nawet przy uczciwej monecie i równych stawkach każdy skończony gracz stający naprzeciw nieskończenie bogatego kasyna jest w dłuższej perspektywie skazany na ruinę.

Proces Wienera i rachunek stochastyczny

Proces Wienera (standardowy ruch Browna) to ciągły w czasie limit prostego błądzenia losowego. Jest on podstawowym elementem nowoczesnego rachunku stochastycznego i leży u podstaw formuły wyceny opcji Blacka-Scholesa, modelu Ornsteina-Uhlenbecka procesów powracających do średniej oraz równania dyfuzji w fizyce. Symulacja renderuje jednocześnie trzy powiązane procesy:

• Standardowy ruch Browna W(t): niezależne przyrosty gaussowskie o wariancji dt.
• Geometryczny ruch Browna S(t) = S(0) exp((μ - σ²/2)t + σW(t)): kanoniczny model cen akcji, ze znosem μ i zmiennością σ.
• Proces Ornsteina-Uhlenbecka dX = θ(μ - X)dt + σdW: powracający do średniej, stosowany w modelach stóp procentowych i fizyce cząstek koloidalnych.

Obwiednie ścieżek próbek, wzrost wariancji empirycznej i rozkład stacjonarny procesu OU są wyświetlane obok ścieżek, dzięki czemu przewidywania teoretyczne można bezpośrednio zweryfikować w tempie interaktywnym.

Regresja do średniej

Francis Galton odkrył regresję do średniej w 1886 roku, badając wzrost rodziców i ich dorosłych dzieci. Wysocy rodzice mają zazwyczaj dzieci niższe od siebie; niscy rodzice zazwyczaj mają wyższe dzieci. Nie dzieje się to z powodu jakiegoś biologicznego mechanizmu korekty, ale dlatego, że wartości ekstremalne w każdym zaszumionym procesie pomiarowym są częściowo dziełem przypadku — a przypadek się nie utrzymuje. Efekt ten jest wszechobecny: tłumaczy, dlaczego „klątwa okładki Sports Illustrated" zdaje się przynosić pecha sportowcom, którzy pojawiają się na okładce, dlaczego uczniowie, którzy uzyskali najwyższe wyniki w teście wstępnym, zazwyczaj poprawiają się mniej niż przeciętnie w teście końcowym, oraz dlaczego wiele interwencji medycznych wydaje się skuteczne, nawet gdy takie nie są.

Symulacja generuje dwuwymiarowe dane normalne o regulowanej korelacji r i pokazuje linię regresji, eliptyczne kontury równej gęstości oraz klasyczną galtonowską strzałkę „regresji ku przeciętności". Towarzyszący panel demonstruje efekt w scenariuszu powtarzanych pomiarów, dzięki czemu mechanizm staje się namacalnie jasny, a nie tylko abstrakcyjnie opisany.

Dynamika płynów: przepływ lepki w rurze

Profile przepływu w rurze

Przepływ lepkiego płynu przez okrągłą rurę pod działaniem stałego gradientu ciśnienia to jedno z niewielu zagadnień mechaniki płynów posiadające dokładne rozwiązanie analityczne. Dla przepływu laminarnego równanie Hagena-Poiseuille'a daje paraboliczny profil prędkości:

u(r) = (ΔP / 4ηL) (R² - r²)

gdzie ΔP to spadek ciśnienia, L to długość rury, η to dynamiczna lepkość, R to promień rury, a r to odległość promieniowa od osi. Maksymalna prędkość w środku jest dokładnie dwukrotnością prędkości średniej — wynik, który pojawia się na każdym kursie mechaniki płynów. Gdy liczba Reynoldsa Re = ρuD/η przekracza w przybliżeniu 2300, przepływ przechodzi w turbulencję, a profil spłaszcza się, przyjmując kształt logarytmicznego prawa ściany.

Symulacja wyświetla obok siebie zarówno laminarną parabolę, jak i profil turbulentny, z suwakiem przesuwającym liczbę Reynoldsa przez obszar przejściowy. Linie prądu pokazują trajektorie poszczególnych cząstek płynu; mapa kolorów przekroju przedstawia lokalną wartość prędkości wewnątrz rury. Ostre przejście między dwoma reżimami — jeden z wielkich nierozwiązanych problemów fizyki klasycznej — jest widoczne jako gwałtowna zmiana kształtu profilu prędkości.

Algorytmy: korekcja błędów i teoria kodowania

Kody Hamminga z korekcją błędów

Richard Hamming wynalazł swoje kody z korekcją błędów w 1950 roku z frustracji wobec komputera w Bell Labs, który stale odrzucał jego weekendowe zadania wsadowe, gdy tylko napotykał błąd parzystości. Kod Hamminga(7,4) koduje 4 bity danych w 7 bitach, dodając 3 bity parzystości umieszczone na pozycjach będących potęgami dwójki. Każdy bit parzystości obejmuje określony podzbiór pozycji danych, tworząc układ równań, który może wskazać dokładnie, który pojedynczy bit został odwrócony — nie tylko wykryć, że wystąpił błąd, ale go skorygować.

Syndrome = P1 P2 P3   (XOR of covered bit positions)
If syndrome = 0: no error
If syndrome = 5: bit 5 is flipped — correct it

Symulacja przeprowadza cię krok po kroku przez kodowanie, transmisję z wprowadzonymi błędami pojedynczego bitu oraz dekodowanie syndromu. Możesz kliknięciem odwrócić poszczególne bity i obserwować, jak obliczenie syndromu aktualizuje się w czasie rzeczywistym. Dostępne są również rozszerzone kody Hamminga z dodatkowym ogólnym bitem parzystości (SECDED — korekcja jednego błędu, wykrycie dwóch) — to właśnie takie kody są używane w pamięci ECC RAM w każdym nowoczesnym serwerze. Zrozumienie kodów Hamminga jest bramą do kodów Reeda-Solomona (używanych w kodach QR i płytach CD), kodów LDPC (używanych w Wi-Fi i 5G) oraz głębokiego związku między teorią kodowania a przepustowością teorioinformacyjną.

Geometria: wielościany w trzech wymiarach

Eksplorator wielościanów

Wzór wielościanu Eulera V - E + F = 2 to jeden z najstarszych i najpiękniejszych wyników w matematyce. Obowiązuje on dla dowolnego wielościanu wypukłego: odejmij krawędzie od wierzchołków, dodaj ściany — a zawsze otrzymasz 2. Pięć brył Platona — czworościan, sześcian, ośmiościan, dwunastościan, dwudziestościan — to jedyne wielościany wypukłe, których wszystkie ściany są przystającymi wielokątami foremnymi. 13 brył Archimedesa rozluźnia ten warunek, dopuszczając wiele typów ścian przy zachowaniu tranzytywności wierzchołków.

Eksplorator wielościanów pozwala obracać i badać wszystkie 18 tych brył w 3D, przełączając się między widokiem szkieletowym, pełnymi ścianami i widokiem z etykietami wierzchołków. Panel statystyk potwierdza wzór Eulera dla każdej bryły w czasie rzeczywistym podczas przełączania między nimi. Podświetlone są relacje ścinania: kubooktaedr to ścięty sześcian, ikosododekaedr to ścięty dwunastościan, a rombokubooktaedr znajduje się w połowie drogi między sześcianem a ośmiościanem w hierarchii grup symetrii. Te zależności łatwiej dostrzec w interaktywnym widoku 3D niż na jakimkolwiek statycznym diagramie.

Wypróbuj sam

Trzy symulacje, od których warto zacząć, jeśli jesteś nowy w tej fali:

Pełny katalog Fali 108 można przeglądać przez indeks kategorii. Wszystkie dziesięć symulacji jest też dostępnych po ukraińsku poprzez prefiks językowy /uk/.

Co dalej

Fala 108 zamyka sprint skupiony na klasycznych modelach matematycznych, które znajdują się w przestrzeni między matematyką czystą a nauką stosowaną — modelach z dokładnymi rozwiązaniami analitycznymi, które symulacja raczej naświetla, niż zastępuje. Obserwowanie, jak paraboliczny profil Poiseuille'a przechodzi w turbulentne prawo logarytmiczne podczas przeciągania suwaka liczby Reynoldsa, przekazuje coś, czego same równania nie przekazują: gwałtowność przejścia, sposób, w jaki płaski wierzchołek profilu turbulentnego wciąż wygina się przy ścianie, wizualną różnicę między laminarnym porządkiem a turbulentnym chaosem.

Fala 109 jest w fazie planowania. Obecne cele obejmują rezonans stochastyczny, wykładniki krytyczne modelu Isinga, triangulację Delaunaya z interaktywnym wstawianiem punktów oraz powrót do układu Lorenza z obliczaniem wykładnika Lapunowa. Jeśli masz symulację, którą chciałbyś zobaczyć zrealizowaną, lub model, który twoim zdaniem jest niedostatecznie reprezentowany w przeglądarkowej edukacji naukowej, strona kontaktowa jest zawsze otwarta.