Как распечатать фит-параметры свернутых функций? - PullRequest
0 голосов
/ 21 января 2020

Я анализирую спектры, и я установил его с помощью функции Фойгта.

Поскольку функция Фойгта является сложной функцией Гаусса и Лоренца, я всегда получу либо параметры Fit Гаусса, либо Fit Фойгта. параметры.

Есть ли способ распечатать фактические амплитуды и центры функции voigt?

Вот функция Fit, которую я использую:

def _4Voigt(x, ampG1, cenG1, sigmaG1, ampL1, cenL1, widL1, ampG2, cenG2, sigmaG2, ampL2, cenL2, widL2, ampG3, cenG3, sigmaG3, ampL3, cenL3, widL3, ampG4, cenG4, sigmaG4, ampL4, cenL4, widL4, off):
return (ampG1*(1/(sigmaG1*(np.sqrt(2*np.pi))))*(np.exp(-((x-cenG1)**2)/((2*sigmaG1)**2)))) +\
          ((ampL1*widL1**2/((x-cenL1)**2+widL1**2)) +\
           (ampG2 * (1 / (sigmaG2 * (np.sqrt(2 * np.pi)))) * (np.exp(-((x - cenG2) ** 2) / ((2 * sigmaG2) ** 2)))) + \
           ((ampL2 * widL2 ** 2 / ((x - cenL2) ** 2 + widL2 ** 2))) +\
            (ampG3 * (1 / (sigmaG3 * (np.sqrt(2 * np.pi)))) * (np.exp(-((x - cenG3) ** 2) / ((2 * sigmaG3) ** 2)))) + \
            ((ampL3 * widL3 ** 2 / ((x - cenL3) ** 2 + widL3 ** 2))) +\
            (ampG4 * (1 / (sigmaG4 * (np.sqrt(2 * np.pi)))) * (np.exp(-((x - cenG4) ** 2) / ((2 * sigmaG4) ** 2)))) + \
            ((ampL4 * widL4 ** 2 / ((x - cenL4) ** 2 + widL4 ** 2)))) + off

PS : Пожалуйста, не жарите, Гринхорн здесь.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...