🦅 Prawdopodobieństwo — lot Lévy'ego i stabilne rozkłady
Większość procesów losowych omawianych na kursach wprowadzających opiera się na rozkładach o dobrze zdefiniowanej średniej i wariancji: rozkład Gaussa, Poissona, dwumianowy. Jednak duża rodzina zjawisk występujących w rzeczywistości podlega rozkładom, w których wariancja — a nawet średnia — jest nieskończona. Nowa symulacja Lot Lévy'ego oddaje to bezpośrednio w Twoje ręce.
Grube ogony i rozkłady α-stabilne
Rozkład nazywamy rozkładem o grubym ogonie, jeśli jego ogon zanika zgodnie z prawem potęgowym: P(X > x) ~ x−α. Dla α ≤ 2 wariancja jest nieskończona; dla α ≤ 1 nieskończona jest także średnia. Rozkłady te nazywa się α-stabilnymi (lub stabilnymi Lévy'ego) i są to jedyne rozkłady zamknięte ze względu na sumowanie — granica sum niezależnych, jednakowo rozłożonych zmiennych o grubych ogonach jest zawsze rozkładem stabilnym (uogólnione centralne twierdzenie graniczne).
- α = 2: rozkład Gaussa. Obowiązuje centralne twierdzenie graniczne. Błądzenie losowe jest ruchem Browna.
- α = 1,5: rozkład Lévy'ego. Wariancja jest nieskończona. Sporadyczne, bardzo długie skoki.
- α = 1,2: bardzo gruby ogon. Skrajnie rzadkie skoki dominują całkowite przemieszczenie.
- α → 1: granica rozkładu Cauchy'ego. Średnia jest nieokreślona; brak charakterystycznej skali.
Algorytm Mantegny
Generowanie próbek z rozkładu stabilnego nie jest trywialne — nie istnieje jawny wzór na odwrotną dystrybuantę. Algorytm Mantegny (1994) dostarcza wydajnego przybliżenia, łącząc dwa standardowe losowania z rozkładu normalnego ze współczynnikiem korekcyjnym wyprowadzonym z indeksu stabilności α. Symulacja implementuje to w czasie rzeczywistym dla maksymalnie 8 jednoczesnych błądzących cząstek.
Zastosowania statystyki Lévy'ego
Lot Lévy'ego i statystyka kroków opisana prawem potęgowym pojawiają się na wielu skalach: wzorce poszukiwania pożywienia przez zwierzęta (albatrosy, rekiny, pszczoły), dyfuzja w turbulentnej plazmie, ruchy cen akcji notowane transakcja po transakcji, mobilność ludzi na podstawie danych z telefonów komórkowych, sekwencje przemieszczeń podczas trzęsień ziemi oraz rozkłady swobodnej drogi fotonów w atmosferach gwiazd.
Powiązane symulacje w kategorii prawdopodobieństwa obejmują Centralne twierdzenie graniczne, Ruch Browna, oraz Błądzenie losowe.
🦅 Otwórz „Lot Lévy'ego” →🌿 Sztuka generatywna — dyfuzyjna agregacja ograniczona
Dyfuzyjna agregacja ograniczona (DLA) to algorytm z 1981 roku autorstwa Wittena i Sandera, który tworzy wizualnie uderzające fraktalne skupiska na podstawie zaskakująco prostych reguł. Umieść zarodek; uwalniaj cząstki pojedynczo, by błądziły losowo, aż się przykleją. Po tysiącach cząstek powstaje rozgałęziona, drzewiasta struktura, która przeczy intuicyjnym oczekiwaniom wynikającym z tak prostej dynamiki.
Od prostych reguł do geometrii fraktalnej
Kluczowa obserwacja polega na tym, że końcówki rosnącego skupiska mają większą szansę przechwycić nadlatujące błądzące cząstki niż wklęsłe zatoki. Ta niestabilność rozszczepiania końcówek wzmacnia każdą początkową asymetrię: gdy gałąź staje się nieco dłuższa, przechwytuje więcej cząstek i rośnie jeszcze bardziej. Rezultatem jest samopodobne rozgałęzienie na każdej skali — definiująca cecha fraktala.
Wymiar fraktalny D_f ≈ 1,71 oznacza, że skupiska DLA znajdują się w 71% drogi między linią (wymiar 1) a wypełnionym dyskiem (wymiar 2). Wartość ta została potwierdzona w wielkoskalowych eksperymentach numerycznych i różni się od dyfuzji ze skończonym prawdopodobieństwem przyklejenia (DLA z redukcją szumu), która zbliża się do D_f = 2, gdy prawdopodobieństwo przyklejenia dąży do zera.
Tryby kolorów
Symulacja oferuje cztery sposoby wizualizacji historii i struktury skupiska:
- Wiek: cząstki, które przykleiły się najwcześniej, są najgorętsze (czerwone), a najnowsze — najchłodniejsze (niebieskie), co ujawnia historię wzrostu warstwa po warstwie.
- Odległość: promień od początku układu jest mapowany na barwę, co uwypukla symetrię kołową i długości gałęzi.
- Ciepło: lokalna gęstość cząstek moduluje kolor, pokazując różnicę w upakowaniu między końcówkami gałęzi a ich trzonami.
- Mono: jednolite, płaskie wypełnienie kolorem, przydatne do druku lub obserwacji czystej geometrii bez rozpraszania kolorem.
Dlaczego DLA ma znaczenie poza sztuką
DLA jest modelem elektroosadzania (tworzenie dendrytów w bateriach), wzorców przebicia dielektrycznego (figury Lichtenberga), lepkiego palcowania w ośrodkach porowatych, tworzenia dendrytów mineralnych w szczelinach skalnych oraz morfologii wzrostu koralowców. Ten sam mechanizm niestabilności rozgałęziania tworzy podobne geometrie fraktalne w tych bardzo różnych układach fizycznych.
Powiązane symulacje w kategorii sztuki generatywnej obejmują Paproć Barnsleya, Zbiór Mandelbrota, oraz Grę w chaos.
🌿 Otwórz fraktal DLA →🌊 Dynamika płynów — równania płytkiej wody
Równania płytkiej wody (SWE) — nazywane również równaniami Saint-Venanta na cześć Adhémara Jeana-Claude'a Barrégo de Saint-Venanta, który wyprowadził je w 1871 roku — są fundamentalnym modelem w geofizycznej dynamice płynów. Opisują one uśredniony po głębokości przepływ poziomy w akwenach, gdzie głębokość jest znacznie mniejsza niż skala pozioma.
Gdzie stosuje się SWE
Równania SWE rządzą: tsunami (których długość fali wynosząca 100–200 km znacznie przewyższa głębokość oceanu); falami sztormowymi i powodziami przybrzeżnymi; propagacją pływów w ujściach rzek; falami powodziowymi po przerwaniu tamy; przepływem w kanałach irygacyjnych i rzekach; oraz dynamiką atmosferyczną w skali synoptycznej. Założenia (uśrednienie po głębokości, ciśnienie hydrostatyczne) zawodzą jedynie wtedy, gdy prędkość pionowa jest znacząca — załamujące się fale, turbulentne progi wodne i przerwania tamy w polu bliskim wymagają pełnych trójwymiarowych równań Naviera-Stokesa.
Schemat Laxa-Friedrichsa
Równania SWE tworzą nieliniowy hiperboliczny układ równań różniczkowych cząstkowych, który może rozwijać fale uderzeniowe (progi wodne, skoki hydrauliczne). Schemat różnic skończonych Laxa-Friedrichsa radzi sobie z nimi w naturalny sposób dzięki wbudowanej lepkości numerycznej, która wygładza nieciągłość na obszarze jednej do dwóch komórek siatki. Kompromisem jest dokładność przestrzenna pierwszego rzędu i nieznaczne nadmierne wygładzanie ostrych frontów — akceptowalne do celów edukacyjnych i znacznie prostsze w implementacji niż schematy wysokiej rozdzielczości, takie jak MUSCL-Hancock czy HLLE.
Cztery wyjaśnione scenariusze
- Przerwanie tamy: problem Riemanna. W chwili t = 0 usuwana jest bariera między dwiema kolumnami wody o różnej wysokości. Dokładne rozwiązanie analityczne składa się z rozchodzącej się w lewo fali rozrzedzenia (ciągłej) oraz propagującego się w prawo progu wodnego (fali uderzeniowej). Symulacja odwzorowuje to jakościowo i zbiega do dokładnego rozwiązania wraz ze wzrostem rozdzielczości siatki.
- Kołysanie: kosinusoidalny początkowy profil głębokości w zamkniętym zbiorniku wywołuje oscylacje fali stojącej. Okres podstawowy wynosi T ≈ 2L / √(g·h̄) i zależy od średniej głębokości. Fala wygasa powoli na skutek tarcia Manninga i lepkości numerycznej.
- Fala pływowa: sinusoidalny dopływ na lewej granicy wywołuje falę postępową w nieruchomej wodzie. Przy podkrytycznych liczbach Froude'a fala rozchodzi się w prawo z prędkością c = √(gh). Częstotliwość i okres są widoczne na panelu prędkości.
- Skok hydrauliczny: wysokoprędkościowy dopływ w lewej jednej trzeciej dziedziny wraz ze schodkowym wzrostem batymetrii testuje sposób, w jaki schemat obsługuje człony źródłowe nachylenia dna oraz przejście między przepływem podkrytycznym a nadkrytycznym.
Liczba Froude'a
Liczba Froude'a Fr = |u| / √(gh) odgrywa w przepływach płytkiej wody taką samą rolę, jak liczba Macha w aerodynamice ściśliwej. Fr < 1 oznacza przepływ podkrytyczny (spokojny); Fr > 1 oznacza przepływ nadkrytyczny (szybki). Przejście między nimi to skok hydrauliczny — gwałtowny wzrost głębokości i spadek prędkości, rozpraszający energię kinetyczną w postaci turbulencji. Panel prędkości oznacza próg Fr = 1 linią przerywaną.
Powiązane symulacje w kategorii dynamiki płynów obejmują Efekt Bernoulliego, Konwekcję Bénarda, oraz Warstwę przyścienną.
🌊 Otwórz Równania płytkiej wody →Odkryj więcej symulacji
Przeglądaj wszystkie symulacje w każdej kategorii wzbogaconej przez Wave 57: