gnuplotでのフィッティングでパラメータを固定したい

Question

とあるサンプルデータをガウス関数でフィッティングしました。（画像）
コードは以下です。

f(x) = a*exp(-(x-b)**2/(2*c**2))+d
a=2000;b=5640;c=120;d=0
fit [0:10000] f(x) "sample.txt" u 1:2 via a,b,c,d
print a,b,c,d
plot "sample.txt" w dots lc "black", f(x) w lines lw 2 lc "red"

フィットすると関数f(x)のa,b,c,dの適切なパラメータを返してくれるのですが、dの値だけ0に完全に固定したいです。

このようなことは可能なのでしょうか？

Answer 1

あまりよく分かっていなくてすみません。

fitの文のところを次のようにして求めるパラメータからdを外せばいいのではないでしょうか?

fit [0:10000] f(x) "sample.txt" u 1:2 via a,b,c

Comments

[mappys]

回答誠にありがとうございます！自分もそれでやってみたのですがダメでした・・・。dを外すとそもそも結果が出力されないんですよね・・・。

[ujimushi_sradjp]

set fit limit 1e-08
とかして収れん演算の終了条件を厳しめにするといいのではないでしょうか
http://takeno.iee.niit.ac.jp/~foo/gp-jman/data/20170823/gnuplot-ja-div/node233.html
とかを参考に

[jbpb0]

> dを外すとそもそも結果が出力されない

適当にデータを用意してやってみたら、できましたよ

「sample.txt」の内容
0 0
1 0
2 3
3 10
4 3
5 0
6 0
7 3
8 10
9 3
10 0
11 0
12 3
13 10
14 3
15 0
16 0

#gnuplotのプログラム
f(x) = a*exp(-(x-b)**2/(2*c**2))+d
a=5;b=14;c=1;d=0
fit [0:16] f(x) "sample.txt" u 1:2 via a,b,c
print a,b,c,d
plot "sample.txt" w points pt 5 lc "black", f(x) w lines lw 2 lc "red"