Enota 4 / 11

Analiza eksperimentalnih podatkov in prilagajanje krivulje

Dobički:

  • Sposobnost prilagajanja eksperimentalnih podatkov s curve_fit z izbiro modela iz fizike in poročanjem parametrov z njihovimi negotovostmi
  • Sposobnost poštenega ocenjevanja kakovosti prileganja s pregledovanjem ostankov in izogibanja pretiranemu opremljanju
  • Sposobnost preverjanja fizične verodostojnosti najdenih parametrov in razlikovanja med "dobrim" prileganjem od "fizično pravilnega" prileganja

Fizika je na koncu odvisna od eksperimentiranja. Ne glede na to, kako elegantna je teorija, je neveljavna, če se ne ujema z izmerjenimi podatki. Prilagajanje krivulje je v središču dela z eksperimentalnimi podatki: prilagajanje modela, ki ga predvideva teorija (črta, eksponenta, nihanje), na hrupne merilne točke, ki jih imate, na najboljši možni način in iskanje parametrov modela (naklon, konstanta, frekvenca). V tej enoti se boste naučili, kako uporabljati umetno inteligenco (AI) pri analizi eksperimentalnih podatkov, kako pravilno nastaviti prileganje krivulje in – kar je najbolj kritično – kako ugotoviti razliko med postavo, ki je »videti dobro« in »biti fizično pravilen«. Standardno orodje je scipy.optimize.curve_fit (funkcija SciPy, ki prilagodi funkcijo modela podatkom).

Logika prilagajanja krivulje in vloga AI

Prilagajanje krivulje naredi to: izberete funkcijo modela (npr. y = a·x + b) in algoritem poišče vrednosti a in b, ki ustrezata podatkom z najmanjšo napako. AI je hiter pri kodiranju pravilne funkcije modela, nastavitvi klica curve_fit in vizualizaciji rezultata. Toda dve ključni odločitvi sta vaši: (1) Kateri model je fizično pravilen? (2) Ali je kakovost prileganja res dobra? Umetna inteligenca bi lahko podatkom prilagodila tudi polinom 10. reda in narisala krivuljo, ki je videti skoraj popolna - vendar bi bilo to fizično nesmiselno pretirano opremljanje.

fizični dogodek

pričakovan model

Najden parameter

gibanje s konstantno hitrostjo

y = a·x + b

naklon = hitrost

radioaktivni razpad

N = N₀·e^(−λt)

λ = konstanta upadanja

pomladna sprostitev

x = A·cos(ωt + φ)

ω = kotna frekvenca

Ohmov zakon

V = R·I

R = upor (naklon)

prosti pad

y = ½·g·t²

g = gravitacijski pospešek

Korak za korakom: prileganje poštene krivulje

1. Izberite model iz fizike, ne iz podatkov. Najprej določite, kateri fizikalni zakon mora upoštevati dogodek. Izbira modela s pregledovanjem podatkov in izbiro "kateri se dobro ujema" je poskus prilagajanja fizike podatkom - nevarno.

2. Vnesite podatke in negotovost. Ob merilnih točkah vnesite tudi merilno napako (sigma parameter) posamezne točke, če obstaja. To je potrebno za natančen izračun prileganja in negotovosti parametrov.

3. Preberite parametre in njihove negotovosti. curve_fit vrne obe vrednosti parametrov in njihove negotovosti, izpeljane iz kovariančne matrike. Poročanje o parametru brez njegove negotovosti je nepopolno: beseda "g = 9,7 ± 0,3 m/s²" je veliko bolj informativna kot le beseda "9,7".

4. Ocenite prileganje z ostanki. Rezidual je odstopanje vsake merilne točke od modela. Naključna porazdelitev ostankov je znak dobrega prileganja; sistematična krivulja (npr. vse na enem mestu nad modelom) kaže, da je model napačen.

5. Preverite fizično razumnost. Ali je vrednost g v vrstnem redu 9,8 ali 50? Ali je upor postal negativen? Fizično nemogoči parametri kažejo na težavo, tudi če prileganje »izgleda dobro«.

Namig: Ne presojajte kakovosti prileganja tako, da se zanašate zgolj na vid. Ostanke vedno narišite na ločen graf. Pri dobrem prileganju se ostanki razpršijo kot hrup okoli ničle; Pri slabem prileganju se pojavi viden vzorec (val, naklon, U-oblika). Ta vzorec je najbolj iskren znak, da je vaš model fizično pomanjkljiv.

trije mini kovčki

1. primer – past za prekomerno opremljanje. En študent je 8 podatkovnim točkam prilagodil polinom 7. stopnje, ki ga je predlagala AI; krivulja je šla skozi vse točke in je bila videti "popolna". Toda med obema točkama je krivulja absurdno nihala. Študent se je obrnil k fiziki: bila je linearna povezava. Ko je namestil preprosto črto, so parametri dobili fizični pomen; polinom višjega reda si je pravkar zapomnil šum.

Primer 2 – Pozabljena negotovost. Raziskovalec je našel in prijavil gravitacijski pospešek g = 9,73 s kodo, ki jo je napisal YZ. Njegov svetovalec je vprašal o negotovosti. Raziskovalec ni uporabil rezultatov kovariance curve_fit; Pri izračunu se je izkazalo, da je negotovost ±0,4, kar pomeni, da je bil rezultat skladen s standardno vrednostjo 9,81. Brez dvoumnosti bi lahko rezultat napačno interpretirali kot "odstopa od 9,81".

Primer 3 – Ostanki so pokazali napako. Učitelj je v eksperiment ohlajanja vgradil ravno črto in zdelo se je razumno. Toda ko je izrisal ostanke, se je pojavil izrazit ukrivljen vzorec. Dejansko se je zgodilo eksponentno ohlajanje. Ko je uporabil pravilen model (eksponentni), so ostanki postali naključni in prileganje je postalo fizično pomembno.

Štiri predloge za kopiranje

1) Prilagoditev fizičnemu modelu:

Napišite kodo curve_fit, ki ustreza fizičnemu modelu = [funkcija modela] na naslednje eksperimentalne podatke ([vrednosti x in y, v enotah]). Natisnite parametre in 1-sigma negotovosti, izpeljane iz kovariance. Narišite podatke, prilagojeno krivuljo in ostanke na LOČENEM grafu. Izbor modela sem naredil iz fizike; samo uporabite.

2) Ocena kakovosti:

Spodaj dodajte kodo, ki ovrednoti kakovost prileganja: narišite ostanke, izračunajte mero, kot je hi-kvadrat ali R², in pisno komentirajte, če je v ostankih sistematičen vzorec. Fit koda: [tukaj]

3) Primerjava modelov:

Napišite kodo, ki ustreza in primerja dva različna fizična modela z istimi podatki: [model A] in [model B]. Pokažite parametre, negotovosti in preostale vzorce obeh. Na meni je, da se odločim, kateri je FIZIČNO bolj primeren; Samo predstavite oba drug ob drugem. Ne pozabite se izogibati prekomernemu opremljanju.

4) Poročanje o negotovosti:

Dodajte kodo, ki natisne vsak parameter naslednjega prileganja z ustreznim številom pomembnih števk v obliki "vrednost ± negotovost (enota)". Prilagodi, če ima parameter fizično nemogočo vrednost (npr. negativna masa). Koda: [tukaj]

Šibek poziv/močan poziv

Slab: "Poiščite krivuljo, ki najbolje ustreza tem podatkom."
Zaključek: AI lahko zagotovi "popolno", a fizično nesmiselno prileganje polinomu visokega reda; negotovost in nič več.
Strong: "Napišite kodo curve_fit, ki ustreza modelu radioaktivnega razpada N = N₀·e^(−λt) tem (t, N) podatkom. Natisnite λ in N₀ ter njune negotovosti. Narišite ostanke v ločen graf in komentirajte, če obstaja sistematičen vzorec. Opozorite, če je λ negativen."
Rezultat: Fizični model, negotovi parametri, preverjanje ostankov in preverjanje verodostojnosti.

Pogoste napake

  • Izbira modela iz podatkov. Izbira "najboljšega prileganja" namesto fizikalnega zakona prisili fiziko, da prilagodi podatke, in vodi do prevelikega prileganja.
  • Ne poročanje o negotovosti. Vrednost parametra brez negotovosti je znanstveno nepopolna; Ali sta dva rezultata združljiva ali ne, je mogoče razumeti le z negotovostjo.
  • Ne praskanje ostankov. Razglasitev skladnosti za "dobro" s pogledom v oči skriva sistematična odstopanja.
  • Precenite pomembno številko. Zapis g = 9,73418 ± 0,4 je napačna natančnost; Negotovost določa, koliko števk je pomembnih.
  • Mimo fizične razumnosti. Ignoriranje nemogočih posledic, kot je negativni upor ali gravitacija 50 m/s².
Pozor: samo zato, ker se krivulja dobro prilega podatkom, še ne dokazuje, da je temeljna fizika pravilna. Z dovolj prostimi parametri se lahko "dobra" krivulja prilega vsem podatkom. Znanstvena integriteta zahteva izbiro modela iz fizike, poročanje o negotovosti in pošten pregled ostankov. Lepa grafika, ustvarjena z umetno inteligenco, ni nadomestilo za te kontrole.

Če povzamem

Pri analizi eksperimentalnih podatkov AI hitro vzpostavi kodo za prileganje krivulje in vizualizacijo; vendar je odgovornost fizika, da izbere fizični model, natančno poroča o negotovosti in pošteno oceni kakovost prileganja. Dobra analiza vzame model iz fizike, poroča o parametrih z njihovo negotovostjo, preuči ostanke in preizkusi fizikalno verodostojnost rezultata. V naslednji enoti bomo podrobneje obravnavali ta koncept negotovosti – merilno napako in širjenje napake.

Aplikacijska naloga

Prilagodite fizični model majhnemu naboru eksperimentalnih podatkov, ki jih imate ali ste jih proizvedli (na primer meritve dolžine-periode nihala ali podatki o tokovni napetosti upora). Natisnite kodo s 1. predlogo v AI, jo zaženite; Preberite parametre z njihovo negotovostjo in narišite ostanke. Preverite, ali obstaja sistematičen vzorec v ostankih, ali so parametri fizikalno sprejemljivi. Zapišite v 5-6 stavkih: ali je bilo prileganje res dobro, kaj je povedala negotovost?

kontrolni seznam

  • [ ] Model sem izbral iz fizike, ne iz podatkov.
  • [ ] Merilno negotovost, če obstaja, sem vključil v prilagajanje.
  • [ ] Poročal sem o parametrih z njihovo negotovostjo in ustrezno pomembno števko.
  • [ ] Ostanke sem narisal na ločen graf in iskal sistematične vzorce.
  • [ ] Preveril sem, ali so parametri fizično razumni.
  • [ ] Izognil sem se prekomernemu opremljanju (nepotrebno veliki svobodi).