使用 lmfit 创建更精细的复合模型答案

【问题标题】：Creating more elaborate composite models with lmfit使用 lmfit 创建更精细的复合模型
【发布时间】：2020-10-11 08:24:22
【问题描述】：

我想用 lmfit 从简单模型创建更精细/复杂的模型。

我有两个功能，例如一个高斯（规范。峰值统一）和洛伦兹（规范。峰值统一）并希望适合例如它们的线性组合。仍然增加了团结。所以我可以写一个像

这样的新函数

def voigt(*all the parameters*, alpha)
   return alpha*gaussian(...) + (1-alpha)*lorentzian(...)

但这不是很适应。所以我现在改为：

mod = ConstantModel(prefix = 'a1_') + Model(gauss) + ConstantModel(prefix='a2_') * Model(lorentz)
pars = Parameters()
pars.add('a1_c', value = 1, min = 0, max = 1)
pars.add('a2_c', expr = '1-a1_c')

还是觉得有点笨拙。有没有更优雅的方式？

【问题讨论】：

标签： python model-fitting lmfit

【解决方案1】：

我认为这主要取决于您正在寻找什么样的复杂性和适应性。如果我理解了这个问题（可能是您真正目标的示例或简化），您希望将两个峰值函数限制为具有总和为固定值（例如 1）的最大高度。为什么不定义随心所欲扩展的基函数呢，也许

 import numpy as np
 def gauss(x, amp, cen, wid):
      return amp * np.exp(-(x-cen)**2/wid)  # deliberately Gaussian-like

 def loren(x, amp, cen, wid): 
      return amp * wid / ( (x-cen)**2 + wid) # deliberately Lorentzian-like

然后将它们相加并将两个amp 参数的和限制为1：

from lmfit import Model
mod = Model(gauss, prefix='g_') + Model(loren, prefix='l_')

pars = mod.make_params(g_cen=1, g_wid=1, l_cen=1, l_wid=1)
pars.add('g_amp', value=0.5, min=0, max=1)
pars.add('l_amp', expr='1 - g_amp', min=0, max=1)

【讨论】：