Pythonの最適化関数

Pocket

PythonでのLebenverg-Marquardt

PythonでNRTLの平衡曲線を書くところまではできたので、次はやはり実データとのフィッティングがしたいと思い調べていたら、Scipy.optimizeのleastsq関数でできそうなことが分かった。

SciPy.orgに詳細が載っていました。)

NRTLのバイナリデータを合わせこむことを考えると、式形は分かっているのでJacobianを自分で求めてから計算した方が早そうな気がしたものの、まずはマニュアルに従って猿真似をしてみる。

サンプルはy=Asin(2πkx+θ)に対してフィッティングを行う例になっている。(パラメーターはA,k,θ)

ちょっと計算に時間がかかるものの、これは便利。

 

Pocket

コメント

タイトルとURLをコピーしました