potass' blog

ポタシウムのことが書いてないブログ。

数値計算

Scipy.odr を使った陰関数フィッティング

実行環境 OS : Windows 8.1 (64 bit) Python : Python 3.2 (Scipy 0.14) やりたいこと 陰関数(implicit function)を fitting したくなった。 陽関数(explicit function)なら最小二乗法(least square)を使えばいいんだがどうしたものか。 そもそも Igor…

Pade 近似の係数を求めるプログラムを C で書いてみた

Pade 近似とは 関数 を有理関数 *1で近似する方法。 次までのマクロリーン展開の係数が両者で一致することが条件。 詳しくは参考 URL まで。 アルゴリズム の 次のマクロリーン展開係数を と書く。 また、 と の多項式の係数をそれぞれ 、 とする。*2 Pade …