Naúrar 3 / 11

Kwaikwayo na Tsarin Jiki

Riba:

  • Ikon kwaikwayi tsarin jiki mai ƙarfi tare da masu warwarewa kamar solve_ivp da tabbatar da ingancin simintin ta hanyar saka idanu gabaɗayan makamashi ko kiyayewa mai ƙarfi.
  • Ikon bambance kayan tarihi na lambobi daga kimiyyar lissafi ta ainihi ta zaɓi girman matakin tare da gwajin haɗuwa
  • Ikon gwada sakamakon kwaikwayo ta hanyar kwatanta shi da sanannun yanayi na nazari da kuma iyakokin jiki

Yawancin ilimin kimiyyar lissafi suna hulɗar da tsarin canza lokaci: pendulum mai juyawa, duniyar da ke kewayawa, jiki mai sanyaya, igiyar ruwa mai yaɗawa. Yawancin waɗannan tsarin ana siffanta su ta hanyar ma'auni daban-daban - wato, ƙididdiga waɗanda ke da alaƙa da adadin canjin adadi (matsayi, zafin jiki) zuwa wannan adadin da kansa. Wadannan ma'auni sau da yawa ba su da mafita na nazari (daidaitacce); a maimakon haka muna kwatanta: ciyar da tsarin gaba a cikin ƙananan matakai da kuma lura da halayensa a lambobi. A cikin wannan rukunin, zaku koyi yadda ake amfani da hankali na wucin gadi (AI) don gina lambar simulation da - mafi mahimmanci - yadda ake tantancewa tare da dokokin kiyayewa ko simintin yana nuna daidaitaccen ilimin lissafi.

Asalin ra'ayi na kwaikwayo da kuma rawar AI

Simulation ya ƙunshi wannan sake zagayowar: ɗauki halin yanzu, yi amfani da ma'aunin canji, ci gaba ta ƙaramin Δt, maimaita. AI yana da sauri sosai wajen gina wannan madauki, yana zaɓar madaidaicin mai warwarewa - algorithm wanda ke warware bambance-bambancen mataki-mataki, da gyara lambar. Daidaitaccen kayan aikin physicist shine scipy.integrate.solve_ivp (aikin SciPy wanda ke warware matsalolin ƙimar farko). Amma kafin ku amince da simintin da AI ya ƙirƙira, tambayar da ya kamata ku yi ita ce: Shin wannan simintin yana adana adadin jiki (makamashi, kuzari, kuzarin angular) waɗanda ke buƙatar kiyayewa?

Tsari

Nau'in daidaitawa

Ma'aunin tabbatarwa

pendulum, bazara

Digiri na biyu ODE

Tsayawa na jimlar makamashi

duniyar tazara

Newtonian gravitational ODE

Angular kuzarin + kiyaye makamashi

Sanyi/zafi

Babban darajar ODE

Juyawa zuwa ma'auni zafin jiki

lalatawar rediyo

Bayanin ODE

Ikon rabin rayuwa

Monte Carlo (bazuwar)

kididdiga

Juyawa zuwa sanannen ma'ana/rarrabuwa

Lura: ODE yana nufin "ma'auni na yau da kullun" - ma'auni wanda ya ƙunshi abin da aka samo asali dangane da maɓalli ɗaya (yawanci lokaci).

Mataki-mataki: abin kwaikwaya abin dogaro

1. Bayyana ilimin lissafi da yanayin farko. Menene ma'auni na tsarin? Menene matsayin farawa, gudu, zazzabi? Menene raka'a? Ba da waɗannan a sarari ga AI.

2. Zaɓi mai warwarewa da ya dace da girman mataki. Matsayi mai girma da yawa yana sa simintin ya zama mara tsayayye (sakamakon fashe ko ilimin kimiyyar lissafi ya karya); Ƙananan mataki yana raguwa ba dole ba. AI na iya ba da shawara, amma sarrafa yanke shawara.

3. Bi dokar kiyayewa a lamba. Yi ƙididdigewa da buga jimillar kuzari (ko kuzari) a cikin simintin. Idan ya yi nisa da girman da ya kamata ya kasance dawwama, simulation ɗin ba abin dogaro bane.

4. Kwatanta da sanannen bayani. Misali, lokacin ƙaramin pendulum mai girma ana saninsa ta dabarar T = 2π√(L/g); Kwatanta lokacin da simintin ya bayar da wannan.

5. Canja siga, gwajin hali. Shin motsin motsi yana datsewa lokacin da aka ƙara juzu'i? Shin kewayawa yana canzawa kamar yadda ake tsammani lokacin da yawan taro ya karu? Gwaji tare da fahimtar jiki.

Tukwici: Ƙara “maballin kuzari” zuwa kowane siminti: ƙididdige jimlar makamashi a kowane mataki kuma ƙirƙira shi a matsayin rabo zuwa ƙarfin farko. A cikin tsarin rashin daidaituwa, wannan rabo ya kamata ya ci gaba da kasancewa kusa da 1. Drift na fiye da 1% alamar cewa kana buƙatar rage girman mataki ko canza mai warwarewa.

uku mini lokuta

Case 1 - Fashewar kuzari. Wani ɗalibi ya kwaikwayi kewayawar duniya tare da hanyar Euler mai sauƙi (hanyar matakin farko) wanda AI ta rubuta. A hankali sararin kewayawa ya fadada akan lokaci kuma duniyar "ta tashi". Yayin da ɗalibin yake lura da yawan kuzarin, sai ya ga cewa yana ƙaruwa koyaushe—aibi sananne na hanyar Euler. Ya tambayi AI don mai warwarewa wanda ke adana makamashi mafi kyau (hanyar da ta fi dacewa tare da solve_ivp) kuma sararin samaniya ya zauna cikin kwanciyar hankali.

Case 2 - Kuskuren naúrar, lokacin kuskure. Wani malami ya gano lokacin ya zama 0.2 s a cikin simintin pendulum, amma ƙimar da ake tsammani shine 2 s. Lokacin da ya bincika shi, ya gane cewa AI ta shiga tsayi a cikin santimita kuma yayi amfani da shi azaman mita a cikin dabara. Tare da gyara naúrar, simintin ya yarda da ƙimar nazari na T = 2π√(L/g).

Case 3 - Tabbatar da Monte Carlo. Wani mai bincike ya tambayi AI don lambar da ke ba da samfuran bazuwar (Monte Carlo - simintin ƙididdiga tare da lambobin bazuwar) barbashi suna bugun mai ganowa. Don tabbatar da sakamakon, ya gudanar da shari'ar mai sauƙi da aka sani ta hanyar nazari (misali kimantawa π ta rabon yanki na da'irar) ta hanya guda; Lokacin da lambar tayi annabta daidai π, amincinta ga ainihin simintin ya ƙaru. Ya kuma lura cewa yayin da adadin samfuran ya karu, sakamakon ya ragu kuma ya haɗu zuwa ainihin ƙimar.

Samfura huɗu masu kwafi

1) Simulation mai kula da kiyayewa:

Rubuta Python code wanda ke kwatanta tsarin jiki mai zuwa tare da solve_ivp: [tsarin, daidaito, yanayin farko, raka'a]. A cikin simintin, ƙididdige jimillar makamashi (da ƙarfin aiki, idan akwai) a kowane mataki kuma ƙirƙira rabonsa zuwa ƙimarsa ta farko akan lokaci. Lura cewa a cikin yanayin rashin daidaituwa, wannan rabo ya kamata ya kasance akai ~1.

2) Kwatancen nazari:

Ƙara lambar da ke kwatanta sakamakon simintin mai zuwa zuwa sanannen bayani na nazari na wannan tsarin ([formula]). Sanya su duka akan jadawali ɗaya kuma buga matsakaicin bambanci tsakanin lamba da nazari. Lambar simulation: [nan]

3) Gwajin girman mataki / kwanciyar hankali:

Rubuta lambar da ke gudanar da simintin mai zuwa tare da matakai daban-daban na lokaci (misali dt = 0.1, 0.01, 0.001) kuma yana nuna yadda sakamakon ke canzawa. Ƙayyade bayan wane girman matakin sakamakon ya daidaita (converges). Code: [a nan]

4) Gwajin jiha iyaka iyaka:

Ba da shawarar gwaje-gwajen iyaka na jiki guda 3 don siminti mai zuwa: (misali makamashi ya kamata a kiyaye shi lokacin da gogayya ta kasance sifili, tsayawa da sauri lokacin da juzu'i ya yi girma, yadda lokacin ya kamata ya canza lokacin da taro ya ninka). Rubuta sakamakon da ake tsammanin kowane gwaji a cikin jumla ɗaya. Code: [a nan]

Rauni mai ƙarfi / Ƙarfi mai ƙarfi

Rarrauna: "Rubuta simulation na pendulum."
Sakamako: lamba ba tare da raka'a ba, babu ingantaccen aiki, babu cak ɗin kiyayewa; Babu tabbas ko yana nuna ilimin kimiyyar lissafi ko a'a.
Ƙarfafa: "Rubuta lambar warware_ivp wanda ke kwatanta ƙaramin pendulum na tsawon L = 1 m, fara kusurwar 10 °, tare da g = 9.81 m / s² don sakan 10. Kula da jimillar makamashi a kowane mataki kuma ƙirƙira rabo zuwa farkon. Kwatanta lokacin da aka auna tare da ƙaramin ma'auni mai girma T = 2π√" da (L / bambance-bambancen).
Sakamako: ingantaccen siminti tare da raka'a, kulawar kiyayewa, da kwatance tare da ƙimar nazari.

Kuskuren gama gari

  • Ba bin dokar kiyayewa ba. Idan makamashi ko motsi yana motsawa, simulation ba daidai ba ne; Amincewa da sakamakon ba tare da ganin wannan shine kuskuren da aka fi sani ba.
  • Zaɓin girman matakin ba bisa ka'ida ba. Babban mataki yana haifar da rashin yanke shawara, ƙananan mataki yana haifar da tsadar da ba dole ba; Gwajin haɗuwa yana da mahimmanci.
  • Ba yin kwatancen nazari ba. Ba tare da kwatanta shi da sanannen shari'a ta musamman (ƙaramin girma, yanayin juzu'i ba), simintin ya kasance mara kulawa.
  • Kuskure aibi na lamba don ilimin lissafi. Halayyar kamar fadada orbital sau da yawa aibi ne na hanyar, ba kimiyyar lissafi ta gaske ba.
  • Gyara yawan samfurori a cikin simintin bazuwar. Yana da yaudara don ƙara yawan samfurori a Monte Carlo kuma dogara da sakamako guda ɗaya ba tare da ganin haɗuwa ba.
Tsanaki: Domin kawai simulation yana samar da "jari mai kyau" ba yana nufin yayi daidai ba. Abin gani yana da gamsarwa amma yaudara. Kafin amincewa da simintin, tabbatar da bin dokar kiyayewa kuma kwatanta ta da wani sanannen yanayi. Ƙimar da ba a tabbatar da ita ba shaida ba ce ta zahiri.

A takaice

Kwaikwayo hanya ce mai ƙarfi don fahimtar tsarin jiki mai ƙarfi waɗanda ba su da mafita na nazari kuma da sauri suna gina lambar kwaikwaiyo AI. Amma darajar simintin ya dogara da ko yana nuna daidaitaccen ilimin lissafi. Hanyar tabbatar da hakan ita ce bin dokokin kiyayewa, zaɓi girman matakin da ya dace ta hanyar gwajin haɗuwa, kwatanta sakamakon da sanannun shari'o'in nazari, da gwada iyakokin iyaka na jiki. A cikin raka'a ta gaba za mu matsawa mayar da hankali daga bayanan da aka samar zuwa nazarin ainihin bayanan gwaji.

Aikin aikace-aikace

Zaɓi tsari mai sauƙi mai ƙarfi (pendulum, faɗuwa kyauta + juriyar iska ko sanyaya). Tare da samfuri 1, buga kuma gudanar da lambar siminti wanda ya haɗa da kiyayewa ko duba ma'auni a cikin AI. Bincika ko adadin da aka sa ido (makamashi ko zafin jiki) yana aiki kamar yadda aka zata. Sa'an nan canza siga (gwagwarmaya, taro, girman mataki) kuma duba ko canjin hali ya yarda da hankalin ku na zahiri. Rubuta shi a cikin jimloli 5-6.

jerin abubuwan dubawa

  • [ ] Na fayyace ma'auni, yanayin farko da raka'a na tsarin.
  • [ ] Na bi dokar kiyayewa (makamashi/momentum) a cikin simintin.
  • [ ] Na tabbatar da girman matakin tare da gwajin haɗuwa.
  • [ ] Na kwatanta sakamakon tare da sanannen yanayin nazari.
  • [ ] Na yi aƙalla gwajin jiha iyaka na jiki ɗaya.
  • [ ] Na bambanta lahani na lambobi daga ainihin kimiyyar lissafi.