lmfit voigtmodel eval()保持对设置不同参数的完全相同值?

发布于 2025-02-10 00:59:23 字数 996 浏览 3 评论 0原文

我有一个有关voigtmodel的eval()函数的快速问题。我有一些不同的参数集,并且正在尝试在一组X值下进行评估。一个简单的演示在下面。不知何故,无论我使用哪种参数集,eval()都会给我相同的值。

from lmfit.models import VoigtModel
import numpy as np
import matplotlib.pyplot as plt

p1={'v0_amplitude': 0.7895507911965995, 'v0_center': 0, 'v0_sigma': 0.3115195639177176, 'v0_gamma': 0.3115195639177176}
p2={'v1_amplitude': 0.2433671058552931, 'v1_center': 0, 'v1_sigma': 0.5154224820790644, 'v1_gamma': -0.10060349484692119}
p3={'v2_amplitude': 0.2560477563760484, 'v2_center': 0, 'v2_sigma': 0.19387399410496087, 'v2_gamma': 0.5953440400428517}

p=[p1, p2, p3]
v1=VoigtModel(prefix='v0_')
v2=VoigtModel(prefix='v1_')
v3=VoigtModel(prefix='v2_')

v = [v1, v2, v3]
x = np.linspace(-10,10, 100)

for i in range(3):
  y = v[i].eval(p[i], x=x)
  plt.plot(x, y)
plt.show()

它给我的正好是三个重叠的配置文件三个重叠线形线的图形。但是,显然,它们都应该不同。我不确定它到底出了什么问题。有什么帮助吗?

I have a quick question regarding the eval() function for the VoigtModel. I have a few different parameters sets and am trying to evaluate under a set of x values. A simple demonstration is down below. Somehow, no matter what parameter set I use, the eval() keeps giving me the same values.

from lmfit.models import VoigtModel
import numpy as np
import matplotlib.pyplot as plt

p1={'v0_amplitude': 0.7895507911965995, 'v0_center': 0, 'v0_sigma': 0.3115195639177176, 'v0_gamma': 0.3115195639177176}
p2={'v1_amplitude': 0.2433671058552931, 'v1_center': 0, 'v1_sigma': 0.5154224820790644, 'v1_gamma': -0.10060349484692119}
p3={'v2_amplitude': 0.2560477563760484, 'v2_center': 0, 'v2_sigma': 0.19387399410496087, 'v2_gamma': 0.5953440400428517}

p=[p1, p2, p3]
v1=VoigtModel(prefix='v0_')
v2=VoigtModel(prefix='v1_')
v3=VoigtModel(prefix='v2_')

v = [v1, v2, v3]
x = np.linspace(-10,10, 100)

for i in range(3):
  y = v[i].eval(p[i], x=x)
  plt.plot(x, y)
plt.show()

what it gives me is exactly three overlapped profiles figure of three overlapping lineshapes. However, obviously, they all should be different. I am not sure where it really goes wrong. Any help maybe?

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。
列表为空,暂无数据
我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
原文