Mittelineaarne dünaamika (LOFY.04.070)
Kursuse läbiviimine 2022/2023 kevadsemestril
Toimumisaeg ja asukoht
Loengud ja lõpuseminar toimuvad kolmapäeviti kell 12:15-13:45 Physicumis ruumis A101. Semestri viimasel kolmapäeval, 24. mail, toimub loengu asemel üliõpilaste seminar.
Eksami sooritamine ja hinde kujunemine
Kursuse hinne koosneb kolmest osast:
- Koduülesanded: Kokku on neli ülesannet. Iga ülesanne annab kuni 15 punkti. Vastuseid tuleb esitada e-meili teel.
- Lõpuseminari eest antakse kuni 20 punkti.
- Semestri lõpus toimub kirjalik eksam, mis annab kuni 20 punkti.
Teemad
-
Sissejuhatus
- Kaose ajalugu ja põhimõisted.
- Näide: Kaootilised pendlid.
- Näide: Pohl'i ringpendel.
- Näide: Elektrooniline ostsillaator.
- Näited: Lorenzi atraktor ja Rössleri atraktor.
-
Pidevad dünaamilised süsteemid
- Mõisted: dünaamiline süsteem, faasiruum, faasitrajektoor, püsipunkt.
- 1-mõõtmeline süsteem ja püsipunktide stabiilsus.
- Normaalkujul lineaarne süsteem: normaalkujul dünaamilise süsteemi lahendid ja püsipunktide tüübid.
- 2-mõõtmeline üldine lineaarne süsteem: süsteemi viimine normaalkujule ja lahendite leidmine, reaalsete omaväärtuste juht, komplekssete omaväärtuste juht, vastavad püsipunktide tüübid ja faasiportreed.
- N-mõõtmeline üldine lineaarne süsteem: erinevate omaväärtuste juht.
- Mittelineaarsed süsteemid: võrrandite lineariseerimine püsipunktide ümbruses, ohud lineariseerimisel, Hartmani-Grobmani teoreem.
- Perioodilised ja kvaasiperioodilised liikumised. Piirtsüklid. Poincaré lõige.
-
Matemaatiline vahepala maatriksitest
- Maatriksi omaväärtused ja omavektorid.
- Maatriksi diagonaliseerimine, maatriksi Jordani kanooniline kuju.
- Eksponent maatriksist.
-
Diskreetsed dünaamilised süsteemid
- Näide: Logistiline kujutus: püsipunktid, perioodilisus ja itereeritud kujutus.
- Bifurkatsioonid. Feigenbaumi arvud.
- Universaalsus. Teiste süsteemide Feigenbaumi arvud. Üleminek kaosesse.
- Lyapunovi eksponent ja sõltuvus algtingimustest.
-
Kaos
- Lahendite stabiilsus.
- Bifurkatsioonid ja Lyapunovi eksponentid pidevates süsteemides.
- Determineeritud kaos.
- Kaootiliste protsesside uurimismeetodid.
- Kaose klassikalised näited.
- Kaoseteooria rakendused.
- Kaootilise võnkumise uurimine.
- Kaootiline piljard.
-
Fraktalid
- Tuntumad näited: Cantori tolm, Kochi lumehelbeke, Sierpinski vaip.
- Fraktalite mõõtmine: Sarnasusdimensioon, Hausdorff-Besikovitši dimensioon.
- Itereeritud funktsiooni süsteemid (IFS). Enesesarnasus.
- Veidrad atraktorid. Henon atraktor.
- Fraktaalsed basseiniääred.
- Komplekssed kujutused: Julia hulgad, Mandelbroti hulk.
Õppematerjalid
Jupyteri töölehed
Aktiivse sisu Jupyteri formaadis vaatamiseks on vajalik Jupyteri keskkond, mida saab alla laadida Jupyteri kodulehelt. Samuti kasutatakse SciPy pakette SymPy, NumPy ja Matplotlib. Kõige lihtsam on neid kasutada TÜFI Jupyter Hubi keskkonnas, kuhu saab sisse logida TÜ keskse kasutajatunnusega.
Konspektid
Teised materjalid
Videod
Koduülesanded
Pidevad dünaamilised süsteemid
- Valige dünaamilist süsteemi. Süsteem peab olema mittelineaarne ja vähemalt kahemõõtmiline. Lisaks peab tal olema vähemalt kaks püsipunkti. (2 p.)
- Uurige oma süsteemi ja leidke vähemalt kaks püsipunkti. (2 p.)
- Arvutage Jacobi maatriksi nii üldises faasiruumi punktis kui püsipunktides. (3 p.)
- Leidke Jacobi maatriksi omaväärtusi püsipunktides ja iseloomustage püsipunktide stabiilsust. (4 p.)
- Joonistage faasiportree, mis näitab ka püsipunkte, ja seletage, kuidas eelnevalt leitud püsipunktide stabiilsus ilmneb faasiportrees. (4 p.)
Tähtaeg: 7. märts 2023 kell 23:59
Diskreetsed dünaamilised süsteemid
- Valige endale diskreetset dünaamilist süsteemi \(x_{n + 1} = c \cdot f(x_n)\), mis rahuldab tingimusi: (2 p.)
- Funktsioon \(f\) on sile, \(f \in C^{\infty}(\mathbb{R})\).
- \(f(0) = f(1) = 0\).
- Funktsioonil \(f\) on maksimum \(\bar{x} \in (0, 1)\) nii et \(f(\bar{x}) = 1\).
- Funktsioon \(f\) kasvab monotoonselt vahemikus \([0, \bar{x}]\) ja kahaneb monotoonselt vahemikus \([\bar{x}, 1]\).
- Funktsioon \(f\) rahuldab \(\frac{f'''(x)}{f'(x)} - \frac{3}{2}\left(\frac{f''(x)}{f'(x)}\right)^2 < 0\) kogu vahemikus \(x \in [0, 1]\).
- Arvutage Lyapunovi eksponenti \(\lambda\) parameetri vahemikus \(c \in [0, 1]\) ja joonistage \(\lambda(c)\). (3 p.)
- Joonistage Feigenbaumi diagrammi ja selgitage lühidalt seost Lyapunovi eksponendiga. (4 p.)
- Leidke esimest 6 ülistabiilset punkti \(s_1, \ldots, s_6\) ja vastavaid \(\delta_n = \frac{s_{n + 1} - s_n}{s_{n + 2} - s_{n + 1}}\). (4 p.)
- Mis parameetri \(c\) väärtuses toimub üleminek kaosesse? (2 p.)
Tähtaeg: 4. aprill 2023 kell 23:59
Kaos
- Valige endale 3-mõõtmelist dünaamilist süsteemi, mis sõltub vähemalt ühest parameetrist. Näiteks: (3 p.)
- Rössleri süsteem,
- Lorenzi süsteem,
- ise välja mõeldud või teisest allikast leitud süsteem.
- Uurige ja kirjeldage süsteemi käitumist sõltuvalt ühest parameetrist. Kui süsteemil on mitu parameetrit, varieerige ainult ühte, ja hoidke ülejäänud parameetrid konstantsena. (4 p.)
- Joonistage süsteemi trajektoore 3-mõõtmelises ruumis erinevate parameetri väärtuste puhul ja kirjeldage, kas trajektoor on perioodiline, kaootiline, koondub püsipunktile… (5 p.)
- Otsige parameetri väärtusi, kus toimuvad käitumise kvalitatiivsed muutused - näiteks bifurkatsioonid / perioodi muutused, atraktori stabiilsuse muutused, üleminek kaosesse - ja kirjeldage oma tulemusi. (3 p.)
Tähtaeg: 2. mai 2023 kell 23:59
Fraktalid
- Valige fraktali, mis asub kahe või kolmemõõtmelises ruumis. (2 p.)
- Selgitage, kuidas teie valitud fraktal on konstrueeritud. Kas valitud fraktal on diskreetse dünaamika atraktor, pideva dünaamika atraktor, itereeritud funktsioonisüsteem, lõpmatu hulkade jada piirväärtus, atraktiivsuse basseini äär vms? Lisage andmeid (võrrandeid, programmi), millega on võimalik valitud fraktali konstrueerida. (3 p.)
- Joonistage oma fraktali. (5 p.)
- Põhjendage, miks teie valitud fraktal on fraktal. (5 p.)
Tähtaeg: 30. mai 2023 kell 23:59