Devlog #43 – Fala 23: Topologia, OTW i kosmologia |…

Statystyki platformy po Fali 23

350

Symulacji

Kategorii

130

Postów na blogu

Fal

Spotlightów

Learning

Devlogów

Posty Fali 23

Opublikowane wraz z tą falą

📊 Spotlight #44 — Ogólna teoria względności i kosmologia (∼17 min) 📖 Learning #33 — Analiza zespolona w fizyce (∼16 min)

Nowe symulacje

📊

Wizualizator PCA / SVD

Analiza głównych składowych (PCA) i rozkład według wartości osobliwych (SVD) na żywym wykresie punktowym 2D. Generuj syntetyczny zbiór danych z regulowaną wariancją, kowariancją i wartościami odstającymi. Obserwuj, jak algorytm oblicza macierz kowariancji, diagonalizuje ją i rysuje osie głównych składowych skalowane przez ich wartości własne. Wykres słupkowy pokazuje wyjaśnioną wariancję na składową; suwak pozwala rzutować chmurę punktów na k głównych składowych i rekonstruować przy różnym k, by zobaczyć jakość kompresji.

Zakładka SVD rozkłada dowolną macierz 3×3 M = U Σ Vᵀ i animuje trzy transformacje geometryczne: obrót V, skalowanie Σ, obrót U. Wartości osobliwe są proporcjonalne do szerokości słupków; rząd macierzy to liczba niezerowych wartości osobliwych.

Canvas 2D PCA / SVD Interaktywne

Notatka inżynierska: przyrostowa kowariancja i dekompozycja własna

Macierz kowariancji C = XᵀX/(n−1) (dane o zerowej średniej) jest obliczana w czasie O(n·d²). Dla d = 2 wartości własne podlegają wzorowi zamkniętemu: λ₁,₂ = (σₓ² + σᵧ²)/2 ± √((σₓ² − σᵧ²)²/4 + σₓᵧ²). Wektory własne są obliczane poprzez rotację Jacobiego 2×2 i aktualizowane na żywo podczas przeciągania punktów — wszystko bez żadnej biblioteki algebry liniowej, czystym JS przeglądarki. Interfejs SVD wykorzystuje iteracyjną bidiagonalizację Goluba-Reinscha dla przypadku 3×3, zbieżną w ≤ 10 iteracjach.

PCA i maksymalizacja wariancji C = (1/n) Xᵀ X (macierz kowariancji danych o zerowej średniej X ∈ ℝⁿˣᵈ)
C vₖ = λₖ vₖ (dekompozycja własna)
Wyjaśniony współczynnik wariancji: EVR(k) = λₖ / ∑ λⱼ
Błąd rekonstrukcji: ‖X − X_rekon‖² = ∑_j>k λⱼ (tw. Eckarta-Younga)

🍩

Torus i rodzaj powierzchni

Parametryczne powierzchnie 3D renderowane w Three.js r160: sfera (rodzaj 0), torus (rodzaj 1), podwójny torus (rodzaj 2) i potrójny torus (rodzaj 3). Obracaj myszą, ściskaj, by powiększyć, przełączaj siatkę/bryłę/oba. Panel informacyjny pokazuje rodzaj g, charakterystykę Eulera χ = 2 − 2g oraz ciąg identyfikacji wielokąta fundamentalnego (aba⁻¹b⁻¹ dla torusa).

Sekcja o kompleksach symplicjalnych pozwala przechodzić między poziomami rodzaju, dodając uchwyty; charakterystyka Eulera V − E + F aktualizuje się w czasie rzeczywistym. Twierdzenie o włochatej kuli jest demonstrowane przez próbę „uczesania" pola wektorowego na sferze (rodzaj 0) w porównaniu z torusem (rodzaj 1): torus dopuszcza pole wektorowe bez zer.

Three.js r160 WebGL Topologia

Notatka inżynierska: torus parametryczny z dołączaniem uchwytów

Powierzchnia rodzaju g jest konstruowana przez sklejenie g torusów uchwytowych z bazową sferą. Dla g = 1 stosuje się standardową parametryzację torusa x = (R + r cos v) cos u, y = (R + r cos v) sin u, z = r sin v. Dla g = 2 dwa uchwyty są dołączane w punktach antypodalnych bazowej sfery, każdy uchwyt generowany przez przesunięty fragment torusa i mieszany z gładką funkcją garbu, by zapewnić ciągłość C¹. BufferGeometry w Three.js indeksuje wierzchołki, by współdzielić normalne wzdłuż szwów, utrzymując wywołanie rysowania WebGL na poziomie jednej indeksowanej siatki na powierzchnię.

Charakterystyka Eulera i twierdzenie o klasyfikacji χ = V − E + F (charakterystyka Eulera, niezmiennik topologiczny)
χ = 2 − 2g dla zamkniętej orientowalnej powierzchni rodzaju g
g = 0 (sfera): χ = 2 | g = 1 (torus): χ = 0
g = 2 (podwójny torus): χ = −2 | g = 3: χ = −4

Twierdzenie Gaussa-Bonneta: ∫∫ K dA = 2π χ (łączy krzywiznę z topologią)

Soczewkowanie grawitacyjne (wyróżniona istniejąca symulacja)

Symulacja soczewkowania grawitacyjnego otrzymała w Fali 23 znaczące ulepszenia: masę soczewki można teraz przemiatać od mikrosoczewki o masie Słońca po gromadę galaktyk o masie 10¹⁵ M_☉, z automatycznym przełączaniem między reżimem mikrosoczewkowania (wzmocnienie źródła punktowego) a reżimem silnego soczewkowania (wielokrotne obrazy przeciwstawne i pierścienie Einsteina). Obok płaszczyzny obrazu rysowana jest powierzchnia potencjału Fermata, pokazująca, że dwa główne obrazy odpowiadają minimom i punktom siodłowym Fermata. Dodano kontrolki eliptyczności źródła do wizualizacji ścinania w słabym soczewkowaniu.

Posty na blogu tej fali

Spotlight #44 — Ogólna teoria względności i kosmologia

17-minutowe dogłębne omówienie rozciągające się od zasady równoważności po problem stałej kosmologicznej. Obejmuje: tensor metryczny i geodezyjne; równania pola Einsteina (G_μν + Λg_μν = 8πG/c⁴ T_μν); metrykę Schwarzschilda, horyzont zdarzeń, sferę fotonową (r = 3r_s/2), ISCO (r = 3r_s), rotację Kerra i proces Penrose'a; soczewkowanie grawitacyjne (kąt ugięcia 4GM/c²b, pierścienie Einsteina, Eddington 1919); grawitacyjne przesunięcie ku czerwieni i korekty GPS (+45,9 − 7,2 = +38,7 μs/dobę); metrykę FLRW, równania Friedmanna, Ω_m ≈ 0,315, Ω_Λ ≈ 0,685; piki akustyczne CMB przy ℓ ≈ 220, tłumienie Silka, polaryzację modu B i inflację; równanie stanu ciemnej energii w ≈ −1 i problem dostrojenia rzędu 10¹²³.

Learning #33 — Analiza zespolona w fizyce

16-minutowy ustrukturyzowany kurs. Tematy: płaszczyzna zespolona, wzór Eulera (e^iπ+1=0), wzór de Moivre'a, pierwiastki z jedności; równania Cauchy'ego-Riemanna i funkcje harmoniczne; twierdzenie Cauchy'ego i wzór całkowy (wartości wewnątrz określone przez granicę); szeregi Laurenta, residua, klasyfikacja biegunów; twierdzenie o residuach i zastosowania do klasycznych całek niewłaściwych; odwzorowania konforemne (Möbiusa, twierdzenie o sile nośnej profilu Żukowskiego L = ρUΓ, Schwarza-Christoffela dla obszarów wielokątnych); zastosowania fizyczne (dwuwymiarowy przepływ potencjalny, zespolona prędkość dF/dz = vₓ − ivᵧ; metoda obrazów w elektrostatyce; fazory impedancji prądu zmiennego; inwersje konturowe Fouriera/Laplace'a i propagatory kwantowe); kontynuacja analityczna i hipoteza Riemanna.

Co dalej: Fala 24

Fala 24 skupi się na dynamice płynów, spalaniu i obliczeniach kwantowych. Planowane symulacje obejmują:

Wizualizator szeregu Taylora — przybliżenia wielomianowe sin/cos/e^x z interaktywnym promieniem zbieżności.
Kinetyka Michaelisa-Menten — kinetyka enzymatyczna Vmax, Km, typy inhibitorów, wykres Lineweavera-Burka.
Ulepszenia wizualizatora łańcuchów Markowa — stany pochłaniające, ergodyczność, iteracja potęgowa PageRank.
Równowaga fazowa termodynamiczna — reguła Le Chateliera, ΔG kontra Q/K, przesunięcia temperatury i ciśnienia.
Świergot fali grawitacyjnej — przebieg fali spiralnego zbliżania h(t), przemiatanie częstotliwości, odpowiedź detektora LIGO.

Posty na blogu: Spotlight #45 o dynamice płynów i turbulencji; Learning #34 o termodynamice i mechanice statystycznej.