[Semargl]

Czy jest jakiś znawca #python na sali? xD
Chce dofitować nieliniową funkcję do danych które zebrałem (prawo rozpadu naturalnego). ale cos mi kurde nie działa. Używam curve_fit z scipy. Jakby ktoś chciał pomóc to chętnie przyjmę tą pomoc xD
#fizyka #programowanie

[vito199]

@Semargl wstaw kod bo programistom 15 k się nie chce. Leżą brzuchem do góry

[Semargl]

x = np.array([9412.2296835 , 6175.20253865, 4009.2786237 , 2642.9951883 ,
1711.55676409, 1160.18815197])
y = np.array([15, 31, 47, 63, 79, 95])
def modelfunc(t, N0, T, C):
return N_0 * np.exp(-T*t)+C
optparms, parmcov

[Ilythiiri]

@Semargl: @vito199: Ogólnie ciężko się pomaga z kodem którego nie widać. ( ͡° ʖ̯ ͡°)

[Ilythiiri]

@Semargl: Na boga, użyj jakiegoś pastebina albo githuba bo tutaj kod się rozjeżdża. ᕙ(⇀‸↼‶)ᕗ

[Semargl]

@Ilythiiri: wystarczy walnąć jeden tab przy return xD

[Semargl]

@Ilythiiri: a z bibliotek to:
from scipy.optimize import curve_fit
import numpy as np

[Semargl]

@pascal256: tak, wcześniej dopasowywałem krzywe gaussa w taki sposób

[pascal256]

@Semargl: To w takim razie jakie są symptomy? Jakiś błąd leci? Czy zwraca wyniki bez sensu?

[Semargl]

@pascal256: Covariance of the parameters could not be estimated
category=OptimizeWarning)

[Semargl]

@pascal256: walcze z tyjm od wczoraj i nie moge rozwalić tego błędu :<

[pascal256]

@Semargl: A czy x i y nie powinny być czasami na odwrót?

[Semargl]

@Semargl: a przy dopasowywaniu prostej ( sprawdziłem aktywność próbek promieniotwórxczych i zliczenia na konkretnych kanałach, znając energie dla danych izotopów dopasowłąem prostą ax+b która pozwoliła mi wyznaczyć energie na kanałach i śmigało, potem gauss na pomiarach nienznaego izotopu, wyznaczyłem co to za izotop po energii którą interpolowałem funkcją liniową i tezaz zostało mi tylko oliczyć czas połowicznego rozpadu by mieć pewność co to za izotop)

[Semargl]

powinny

@pascal256: nope xD, zmieniłem i porsta jak leżała tak leży xD

[pascal256]

@Semargl: Rzuć jeszcze okiem na:
https://stackoverflow.com/questions/41725377/why-isnt-curve-fit-able-to-estimate-the-covariance-of-the-parameter-if-the-pa#41732757

[mprzemo]

@Semargl: użyj biblioteki lmfit.

[Semargl]

@pascal256: @mprzemo: dzięki za chęć pomocy. ostatecznie pobrałem qtiplota i tam dopasowałem krzywą. czas połowicznego rozpadu izotopu jodu 128 to 25 minut, z mojego dopasowania wyszło, ze jest to 26, także nie jest jakoś bardzo daleko(chociaż to o 1/26 za dużo xD)

[mprzemo]

@Semargl: jakby wyszło 25 to by było podejrzane :D

[Semargl]

@mprzemo: racja, chociaż myślę, ze mógłbym znaleść gdzieś jakieś niepewności pomiarowe i wyliczyć jakaś sigme, może wtedy wynik mógłbym potwierdzić testem 3sigma i nawet miałbym fajny projekcik z fizyki jądrowej na githuba xd Ale i tak wielkie dzięki ;)