〓 Mathematica 四 〓
■ このスレッドは過去ログ倉庫に格納されています
0645132人目の素数さん
2009/11/08(日) 15:00:02ドキュメントセンターのFindFitの項目に例題がいっぱい載っている。
例えば
data = {{1.0, 12.}, {1.9, 10.}, {2.6, 8.2}, {3.4, 6.9}, {5.0, 5.9}};
model[a_?NumberQ, k_?NumberQ] :=
model[a, k] =
First[x /. NDSolve[{x'[t] + k x[t] == 0, x[0] == a}, x, {t, 0, 6}]]
fit = FindFit[data, model[a, k][t], {a, k}, t]
Show[Plot[model[a, k][t] /. fit, {t, 0, 6}], ListPlot[data]]
のようにすればよい。
■ このスレッドは過去ログ倉庫に格納されています