Мне нужно подогнать некоторые пики с помощью функции Гаусса (bg - это фон)
bg45(x) = a0 + a1*x
k4(x) = (x-x4)/w4
Gauss_r4(x) = A4*exp(-log(2)*k4(x)**2)
Cauchy_r4(x) = A4/(1 + k4(x)**2)
pV4(x) = alpha*Gauss_r4(x) + (1-alpha)*Cauchy_r4(x)
Я определяю исходные параметры подгонки как:
alpha = 0
beta = 0
A4 = 0.71
x4 = 4.28
w4 = 0.49
a0 = -5.67
a1 = 0.51
Тогда я подхожу так:
fit [3.68:5.45] f45(x) filename u 1:2 via x4,A4,w4,a0,a1
Я получаю правильный график с идеальной функцией подгонки, но когда я открываю файл .dat с найденными параметрами, они совершенно неверны.
#id r4 r4e A4 A4e w4 w4e
0 4.186 0.002 4.1 0.0 0.615 0.005
Амплитуда (высота пика А4) должна составлять около 1, а не 4. У вас есть идеи, почему это происходит?