Bölüm 3 / 11

Fiziki ulgamlaryň simulýasiýasy

Girdeji:

  • Solve_ivp ýaly çözüjiler bilen dinamiki fiziki ulgamy simulirlemek we umumy energiýa ýa-da tizligi tygşytlamak arkaly simulýasiýanyň dogrulygyny barlamak ukyby
  • San artefaktlaryny hakyky fizikadan tapawutlandyrmak ukyby, konwergensiýa synagy bilen ädim ululygyny saýlamak
  • Simulýasiýa netijesini belli analitiki ýagdaýlar we fiziki çäk ýagdaýlary bilen deňeşdirip synap bilmek

Fizikanyň köpüsi wagt üýtgeýän ulgamlar bilen iş salyşýar: yrgyldama maýatnik, orbita planetasy, sowadyjy beden, köpelýän tolkun. Bu ulgamlaryň köpüsi diferensial deňlemeler bilen düşündirilýär, ýagny mukdaryň üýtgemeginiň tizligini (pozisiýa, temperatura) şol mukdaryň özi bilen baglanyşdyrýan deňlemeler. Bu deňlemelerde köplenç analitiki (formula-takyk) çözgüt ýok; munuň ýerine simulasiýa edýäris: ulgamy kiçi wagt ädimlerinde öňe sürmek we özüni alyp barşyna san taýdan gözegçilik etmek. Bu bölümde, simulýasiýa koduny gurmak üçin emeli intellekt (AI) ulanmagy we iň esasysy - simulýasiýanyň fizikany takyk görkezýändigini ýa-da gorap saklamak kanunlary bilen nädip barlamalydygyny öwrenersiňiz.

Simulýasiýa we AI-iň roly baradaky esasy pikir

Simulýasiýa bu siklden ybarat: häzirki ýagdaýy alyň, üýtgeşme deňlemesini ulanyň, kiçijik Δt bilen öňe hereket ediň, gaýtalaň. AI bu aýlawy gurmakda, dogry çözüjini - diferensial deňlemäni ädimme-ädim çözýän algoritmi we kody redaktirlemekde gaty çalt. Fizikiň adaty guraly scipy.integrate.solve_ivp (başlangyç baha meselelerini çözýän SciPy funksiýasy). Emma AI tarapyndan döredilen simulýasiýa ynanmazdan ozal bermeli soragyňyz şu: bu simulýasiýa saklanylmaly fiziki mukdarlary (energiýa, tizlik, burç tizligi) saklaýarmy?

Ulgam

Deňlik görnüşi

Barlamak ölçegleri

maýatnik, bahar

Ikinji dereje ODE

Jemi energiýanyň yzygiderliligi

planetanyň orbitasy

Nýuton agyrlyk güýji ODE

Burç tizligi + energiýany tygşytlamak

Sowuk / ýylylyk

Birinji dereje ODE

Deňagramlylygyň temperaturasyna ýakynlaşma

radioaktiw çüýremek

Ekspensial ODE

Halfarym ömri dolandyrmak

Monte Karlo (tötänleýin)

statistika

Belli ortaça / paýlanyşa ýakynlaşma

Bellik: ODE "adaty diferensial deňleme" - ýeke üýtgeýjä (adatça wagt) degişlilikde emele gelen deňlemäni aňladýar.

Stepdimme-ädim: ygtybarly simulýasiýa

1. Fizikany we başlangyç şertleri düşündiriň. Ulgamyň deňlemesi näme? Başlangyç ýagdaýy, tizligi, temperaturasy näme? Bölümler näme? Bulary AI-e aç-açan beriň.

2. Degişli çözüjini we ädim ululygyny saýlaň. Wagtyň gaty uly bolmagy simulýasiýany durnuksyz edýär (netije ýarylýar ýa-da fizika döwülýär); Örän kiçijik bir ädim gereksiz haýallaýar. AI teklip edip biler, ýöne karara gözegçilik edip biler.

3. Kodda gorap saklamak kanunyny ýerine ýetiriň. Simulýasiýa boýunça umumy energiýany (ýa-da tizligi) hasaplaň we çap ediň. Üznüksiz galmaly ululygy süýşürse, simulýasiýa ygtybarly däl.

4. Belli bir çözgüt bilen deňeşdiriň. Mysal üçin, kiçijik amplituda maýatnik döwri T = 2π√ (L / g) formulasy bilen belli; Simulýasiýa bilen berlen döwri munuň bilen deňeşdiriň.

5. Parametri üýtgediň, özüni alyp barşyňyzy synap görüň. Sürtülme goşulanda yrgyldamalar azalýarmy? Massa köpelende orbitanyň garaşylyşy ýaly üýtgeýärmi? Fiziki duýgurlyk bilen synag.

Maslahat: Her simulýasiýa "energiýa yzarlaýjy" goşuň: her ädimdäki umumy energiýany hasaplaň we başlangyç energiýa gatnaşygy hökmünde meýilleşdiriň. Sürtülmeýän ulgamda bu gatnaşygy hemişe 1-e ýakyn saklamaly. 1% -den gowrak süýşmek, ädimiň ululygyny azaltmaly ýa-da çözüjini üýtgetmeli.

üç sany kiçi halat

1-nji waka - energiýanyň ýarylmagy. Talyp, AI tarapyndan ýazylan ýönekeý Euler usuly (iň başlangyç ädim usuly) bilen planeta orbitasyny simulasiýa etdi. Orbit wagtyň geçmegi bilen kem-kemden giňeldi we planeta "uçdy". Talyp umumy energiýa gözegçilik edende, onuň yzygiderli artýandygyny, Euleriň usulynyň belli bir kemçiligini gördi. AI-den energiýany has gowy tygşytlaýan çözgüdi (solve_ivp bilen has amatly usul) we orbitanyň durnukly ellipsde ýerleşmegini isledi.

2-nji waka - Nädogry birlik, nädogry döwür. Maýatnik simulýasiýasynda mugallym 0,2 s diýip tapdy, ýöne garaşylýan bahasy 2 s. Gözden geçirende, AI-iň santimetre uzynlygyna girendigini we formulada metr hökmünde ulanandygyny gördi. Bölümi düzetmek bilen, simulýasiýa T = 2π√ (L / g) analitiki bahasy bilen ylalaşdy.

3-nji waka - Monte Karlo barlagy. Gözlegçi AI-den detektor urýan tötänleýin nusgalary (Monte Karlo - tötänleýin sanlar bilen statistik simulýasiýa) kod sorady. Netijäni barlamak üçin, şol bir usul bilen analitik taýdan bilinýän ýönekeý bir işi (meselem, tegelegiň meýdanynyň gatnaşygy boýunça baha bermek) işledi; Kod π dogry çaklanda, hakyky simulýasiýa bolan ynam artdy. Şeýle hem, nusgalaryň sanynyň köpelmegi bilen netijäniň daralandygyny we hakyky bahasyna öwrülendigini gördi.

Göçürip boljak dört şablon

1) Gorap saklaýyş gözegçiligi:

Solution_ivp bilen aşakdaky fiziki ulgamy simulýasiýa edýän Python koduny ýazyň: [ulgam, deňleme, başlangyç şertler, birlikler]. Simulýasiýanyň dowamynda her ädimde umumy energiýany (we bar bolsa) hasaplaň we wagtyň geçmegi bilen başlangyç bahasyna bolan gatnaşygy meýilleşdiriň. Sürtülmedik ýagdaýda bu gatnaşygy hemişelik ~ 1 saklamalydygyny unutmaň.

2) Analitik deňeşdirme:

Aşakdaky simulýasiýanyň netijesini bu ulgamyň belli analitiki çözgüdi bilen deňeşdirýän kod goşuň ([formula]). Ikisini-de şol bir grafada ýerleşdiriň we san we analitikanyň arasyndaky iň ýokary tapawudy çap ediň. Simulýasiýa kody: [şu ýerde]

3) Stepdimiň ululygy / durnuklylyk synagy:

Dürli wagt ädimleri bilen (mysal üçin dt = 0.1, 0.01, 0.001) aşakdaky simulýasiýany işleýän we netijeleriň üýtgemegini görkezýän kod ýazyň. Netijäniň haýsy ädim uludygyny kesgitläň (öwrüler). Kod: [şu ýerde]

4) Fiziki çäk döwlet synagy:

Aşakdaky simulýasiýa üçin 3 sany fiziki çäkli döwlet synagyny teklip ediň: (meselem, sürtülme nol bolanda energiýa tygşytlanmaly, sürtülme uly bolanda çalt durmaly, köpçülik iki esse artanda döwür nähili üýtgemeli). Her synagyň garaşylýan netijesini bir sözlemde ýazyň. Kod: [şu ýerde]

Gowşak haýyş / Güýçli haýyş

Gowşak: "Maýatnik simulýasiýasyny ýazyň."
Netije: birlik ýok kod, tassyklama ýok, gorap saklamak barlagy ýok; Fizikany görkezýändigi ýa-da ýokdugy belli däl.
Güýçli: "L = 1 m uzynlykdaky maýatnik, 10 ° burçdan başlap, g = 9.81 m / s² bilen 10 sekundyň dowamynda ýönekeý maýatnik şekillendirýän çözgüt_ivp kod ýazyň. Her ädimde umumy energiýany gözegçilikde saklaň we başlangyç gatnaşygyny meýilleşdiriň. Ölçelenen döwri kiçi amplituda formulasy T = 2π√ (L / g) bilen deňeşdiriň we tapawudy çap ediň."
Netije: Bölümler bilen ygtybarly simulýasiýa, gorap saklamak gözegçiligi we analitiki baha bilen deňeşdirmek.

Adaty ýalňyşlyklar

  • Gorap saklamak kanunyna boýun egmezlik. Energiýa ýa-da tizlik azalýan bolsa, simulýasiýa nädogry; Netijäni görmän ynanmak iň köp ýaýran ýalňyşlyk.
  • Basgançagyň ululygyny saýlamak. Örän uly ädim karar bermeýär, gaty kiçi ädim gereksiz çykdajylary döredýär; Konwergensiýa synagy zerurdyr.
  • Analitik deňeşdirmeler etmezlik. Belli bir ýörite ýagdaý (kiçijik amplituda, sürtülmesiz ýagdaý) bilen deňeşdirmän, simulýasiýa gözegçilik edilmän galýar.
  • Fizika üçin san kemçiligini ýalňyşdyrmak. Orbital giňelmek ýaly özüni alyp barş, köplenç hakyky fizika däl-de, usulyň kemçiligi.
  • Tötänleýin simulýasiýada nusgalaryň sanyny düzmek. Monte Karlodaky nusgalaryň sanyny köpeltmek we ýakynlygy görmän ýekeje netijä bil baglamak ýalňyşdyr.
Üns beriň: Simulýasiýanyň "gowy grafik" öndürýändigi sebäpli, munuň dogrydygyny aňlatmaýar. Görüş ynandyryjy, ýöne aldawçy. Simulýasiýa ynanmazdan ozal, tebigaty goramak kanunyna eýeriň we ony belli bir ýagdaý bilen deňeşdiriň. Barlanmadyk simulýasiýa fiziki subutnama däl.

Gysgaça

Simulýasiýa, analitiki çözgütleri bolmadyk we AI simulýasiýa koduny çalt döredýän dinamiki fiziki ulgamlara düşünmegiň güýçli usulydyr. Emma simulýasiýanyň gymmaty, fizikany takyk görkezýändigine baglydyr. Muny üpjün etmegiň usuly tebigaty goramak kanunlaryna eýermek, konwergensiýa synagy arkaly degişli ädim ululygyny saýlamak, netijäni belli analitiki ýagdaýlar bilen deňeşdirmek we fiziki çäk ýagdaýlaryny barlamakdyr. Indiki bölümde ünsi döredilen maglumatlardan hakyky tejribe maglumatlarynyň derňewine geçireris.

Programma meselesi

Simpleönekeý dinamiki ulgamy saýlaň (maýatnik, erkin ýykylmak + howa garşylygy ýa-da sowadyş). 1-nji şablon bilen, AI-de tygşytlylygy ýa-da deňagramlylygy barlaýan simulýasiýa koduny çap ediň we işlediň. Gözegçilik edilýän mukdaryň (energiýa ýa-da temperatura) garaşylşy ýaly hereket edýändigini barlaň. Soňra bir parametr üýtgediň (sürtülme, massa, ädim ululygy) we özüni alyp barşyň üýtgemeginiň fiziki duýgurlygyňyza laýyk gelýändigini ýa-da ýokdugyny görüň. 5-6 sözlemde ýazyň.

gözegçilik sanawy

  • [] Ulgamyň deňlemesini, başlangyç şertlerini we birliklerini düşündirdim.
  • [] Simulýasiýada gorap saklamak kanunyna (energiýa / tizlik) eýerdim.
  • [] Basgançagyň ululygyny konwergensiýa synagy bilen barladym.
  • [] Netijäni belli analitiki ýagdaý bilen deňeşdirdim.
  • [] Iň azyndan bir fiziki çäkli döwlet synagyny geçirdim.
  • [] San kemçiliklerini hakyky fizikadan tapawutlandyrdym.