# import dataset
series = pd.read_csv('F1_00000.csv')
data = DataFrame(series)
deltx = data.iat[1, 0] - data.iat[0, 0]
NTOT = sum(data.iloc[:, 1])
prob = []
lenth = len(data)
for i in range(lenth):
prob.append(data.iat[i, 1] / (NTOT * deltx))
# 画图
plt.plot(data.iloc[:, 0], prob)
# 高斯拟合
g_init = models.Gaussian1D()
fit_g = fitting.LevMarLSQFitter()
g = fit_g(g_init, data.iloc[:, 0], prob)
print(g.mean.value, g.stddev.value)
# 拟合画图
plt.plot(data.iloc[:, 0], prob)
plt.plot(data.iloc[:, 0], g(data.iloc[:, 0]), 'r-', linewidth=3, label='astropy.modeling')
plt.xlabel('$x_2$')
plt.ylabel('$y_2$')
plt.legend()
人家那么大的一行 best_fit_gauss = fitter_gauss(model_gauss, x2, y2, weights=1/y2_err**2),你写在哪里了
试下把g改成fit_g看看,应该是这个问题把,不是的话手头没数据我也不知道问题了。
plt.plot(data.iloc[:, 0], fit_g(data.iloc[:, 0]), 'r-', linewidth=3, label='astropy.modeling')
您好,我是有问必答小助手,您的问题已经有小伙伴解答了,您看下是否解决,可以追评进行沟通哦~
如果有您比较满意的答案 / 帮您提供解决思路的答案,可以点击【采纳】按钮,给回答的小伙伴一些鼓励哦~~
ps:问答VIP仅需29元,即可享受5次/月 有问必答服务,了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632
非常感谢您使用有问必答服务,为了后续更快速的帮您解决问题,现诚邀您参与有问必答体验反馈。您的建议将会运用到我们的产品优化中,希望能得到您的支持与协助!
速戳参与调研>>>https://t.csdnimg.cn/Kf0y
你好,知道原因了吗
分子动力学+量子化学+力场开发专题查询:https://mp.weixin.qq.com/s/Kl-DjmaQXa8TvP6SwniABQ