■パデ近似(その12)

連分数は関数を近似するときにも有効である。

tanzの連分数から2次近似

  tanz=z(15-z^2)/(15-6z^2)

が得られる。z=π/4のとき、0.9998となるが、

tanz=z+z^3/3+z^5/5ではこれよりも32倍も大きい誤差を生じる。

===================================

誤差積分の連分数近似

∫(0,z)exp(-x^2)dx=(49140z+3570z^3+739z^5)/(49140+19950z^2+2475z^4)

はz=2に対して1.2%の誤差であるが、z^9までのベキ級数展開では正しい値より110%も大きくなる

===================================

exp(-x^2)=1-x^2+x^4/2-x^6/6+・・・

∫(0,z)exp(-x^2)dx=x-x^3/3+x^5/10-x^7/42+・・・=(ax+bx^3+cx^5)/(d+ex^2+fx^4)

(ax+bx^3+cx^5)=(d+ex^2+fx^4)(x-x^3/3+x^5/10-x^7/42+・・・)

=dx+(e-d/3)x^3+(d/10-e/3+f)x^5+・・・

a=49140,b=3570,c=739とおく(あるいはa=49140,b=3570,c=739とおかなくても解けるものなのかもしれない)

a=d=49140

b=e-d/3→e=b+d/3=19950

c=d/10-e/3+f→f=c-d/10+e/3=2475

===================================

exp(-x^2)=1-x^2+x^4/2-x^6/6+・・・

∫(0,z)exp(-x^2)dx=x-x^3/3+x^5/10-x^7/42+・・・=(ax+bx^3+cx^5)/(d+ex^2+fx^4)

(ax+bx^3+cx^5)=(d+ex^2+fx^4)(x-x^3/3+x^5/10-x^7/42+・・・)

=dx+(-d/3+e)x^3+(d/10-e/3+f)x^5+・・・

(あるいはa=49140,b=3570,c=739とおかなくても解けるものなのかもしれない)

a=d

b=-d/3+e

c=d/10-e/3+f→fで表すことはできない→連分数が必要であると思われる 

===================================