scipyで最小二乗法
python用の科学計算ライプラリscipyを使った最小二乗法用のテンプレを貼って置きます。
SciPy Cookbook — SciPy Cookbook documentation
を参考にしました。
from scipy.optimize import * fitfunc = lambda p, x: p[0] * x + p[1] #フィットさせる方程式 errfunc = lambda p, x, y: pow((fitfunc(p, x) - y), 2) #方程式とデーター間の距離 p0 = [0, 0] #パラメーターの初期値 p1, success = leastsq(errfunc, p0, args = (x, y)) #最適化するようなパラメーターを受けとる plot(x, y, x, fitfunc(p1, x))