回答此问题可获得 20 贡献值,回答如果被采纳可获得 50 分。
<p>我正在研究将μ子寿命数据拟合到曲线上,用lmfit函数提取平均寿命。我使用的一般过程是使用直方图函数将13000个数据点分为10个单元,用每个单元中计数的平方根计算不确定性(这是一个指数模型),然后使用lmfit模块确定最佳拟合以及平均值和不确定度。但是,绘制模型.拟合()方法返回此图,其中红线表示拟合(显然不是正确拟合)。<a href="https://i.stack.imgur.com/V2JzB.png" rel="nofollow noreferrer">Fit result output graph</a></p>
<p>我在网上找过,找不到解决办法,我真的很感激能帮我弄清楚到底发生了什么。这是密码。在</p>
<pre><code>import os
import numpy as np
import matplotlib.pyplot as plt
from numpy import sqrt, pi, exp, linspace
from lmfit import Model
class data():
def __init__(self,file_name):
times_dirty = sorted(np.genfromtxt(file_name, delimiter=' ',unpack=False)[:,0])
self.times = []
for i in range(len(times_dirty)):
if times_dirty[i]<40000:
self.times.append(times_dirty[i])
self.counts = []
self.binBounds = []
self.uncertainties = []
self.means = []
def binData(self,k):
self.counts, self.binBounds = np.histogram(self.times, bins=k)
self.binBounds = self.binBounds[:-1]
def calcStats(self):
if len(self.counts)==0:
print('Run binData function first')
else:
self.uncertainties = sqrt(self.counts)
def plotData(self,fit):
plt.errorbar(self.binBounds, self.counts, yerr=self.uncertainties, fmt='bo')
plt.plot(self.binBounds, fit.init_fit, 'k--')
plt.plot(self.binBounds, fit.best_fit, 'r')
plt.show()
def decay(t, N, lamb, B):
return N * lamb * exp(-lamb * t) +B
def main():
muonEvents = data('C:\Users\Colt\Downloads\muon.data')
muonEvents.binData(10)
muonEvents.calcStats()
mod = Model(decay)
result = mod.fit(muonEvents.counts, t=muonEvents.binBounds, N=1, lamb=1, B = 1)
muonEvents.plotData(result)
print(result.fit_report())
print (len(muonEvents.times))
if __name__ == "__main__":
main()
</code></pre>