PythonでのLebenverg-Marquardt
PythonでNRTLの平衡曲線を書くところまではできたので、次はやはり実データとのフィッティングがしたいと思い調べていたら、Scipy.optimizeのleastsq関数が便利そうなので試してみました。
（SciPy.orgに詳細が載っていました。）
NRTLのバイナリデータを合わせこむことを考えると、式形は分かっているのでJacobianを自分で求めてから計算した方が早そうな気がしたものの、まずはマニュアルに従って猿真似をしてみる。
サンプルはy=Asin(2πkx+θ)に対してフィッティングを行う例になっている。（パラメーターはA,k,θ）
ちょっと計算に時間がかかるものの、これは便利です。
コメント