0

我试图用给定的数据将曲线拟合到下面的等式。方程是Rate=k*Concentration^n。我n在安装时遇到了麻烦-6,这是不可能的,所以我试图设置一个界限min=0。但是,我收到未定义的术语参数错误。任何帮助都会非常感谢。

from IPython import get_ipython 
get_ipython().magic('reset -sf')
import numpy as np
from lmfit import Model
import matplotlib.pyplot as plt



# Homework Problem #2 
Time = np.array ([0, 48, 76, 124, 204, 238, 289])
Concentration =np.array ([19.04, 17.6, 16.9, 15.8, 14.41, 13.94, 13.37])

# Rate Determination 
Rate=Concentration/Time 


# Model Definition 

def rateEq(Concentration, k, n):
    return k*(Concentration)**n

# Model creation
model=Model(rateEq)

# Parameters
params = parameters()
params.add(k=0.001)
params.add(n=0.001)
par.set(min=0)


# Data Fit v
result=model.fit(Rate, params, Concentration=Concentration)
# Print and Plot Results 
print(result.fit_report())
result.plot_fit()
4

1 回答 1

0

parameters是未定义的,因为您没有在任何地方定义它。您将其用作 params = parameters(),可能暗示函数调用,但您没有定义或导入该函数.... 同样,par未定义,因为您没有在任何地方定义它。

你几乎肯定想要

from lmfit import Model, Parameters # explicitly import Parameters

params = Parameters() # note the capitalization
params.add('k', value=0.001)
params.add('n', value=0.001)

不太清楚(因为我猜不出是什么par)是你是否想要

params['k'].min = 0

或者

params['n'].min = 0

另外,作为警告:由于您的time[0]值为 0,因此您Rate[0]将是无限的,这将在运行 fit 时造成麻烦。

于 2018-02-18T23:24:59.910 回答