sander/rowsandall
Private
Public Access
1
0
Fork 0
Code Activity

constrain fit params to positive vals in CP

Browse Source
This commit is contained in:
Sander Roosendaal
2017-06-16 08:17:30 +02:00
parent b7487f1f77
commit 4b8d9d02f1
1 changed files with 2 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
rowers/interactiveplots.py
View File

@@ -619,7 +619,7 @@ def interactive_otwcpchart(powerdf,promember=0):
# there is no Paul's law for OTW # there is no Paul's law for OTW
# Fit the data to thee parameter CP model # Fit the data to thee parameter CP model
fitfunc = lambda pars,x: pars[0]/(1+(x/pars[2])) + pars[1]/(1+(x/pars[3])) fitfunc = lambda pars,x: pars[0]/(1+(x/pars[2])) + pars[1]/(1+(x/abs(pars[3])))
errfunc = lambda pars,x,y: fitfunc(pars,x)-y errfunc = lambda pars,x,y: fitfunc(pars,x)-y
p0 = [500,350,10,8000] p0 = [500,350,10,8000]
@@ -636,6 +636,7 @@ def interactive_otwcpchart(powerdf,promember=0):
p1 = [p0[0]/factor,p0[1]/factor,p0[2],p0[3]] p1 = [p0[0]/factor,p0[1]/factor,p0[2],p0[3]]
p1 = [abs(p) for p in p1]
fitt = pd.Series(10**(4*np.arange(100)/100.)) fitt = pd.Series(10**(4*np.arange(100)/100.))
fitpower = fitfunc(p1,fitt) fitpower = fitfunc(p1,fitt)