Loading [MathJax]/jax/output/CommonHTML/jax.js

2016年8月19日金曜日

ロドリゲス公式よりも便利な公式

有名な多項式を計算する場合にロドリゲスの公式をつかって計算する人がいるはず。ロドリゲス公式よりもはるかに速く計算できる方法というと漸化式だろうと思う。

漸化式を使わずに、一般ラゲール多項式を計算する公式はある。

L(α)n(z)=nk=0(1)kk!(nk)!zk(α+n)!(α+k)! 便宜上この公式をformula GLとしよう。この公式はロドリゲス公式からライプニッツの公式を使って導出できる。
この公式をつかって L(α)0(z)=1,L(α)1(z)=α+1z,L(α)2(z)=(α+2)(α+1)2(α+2)z+z22,L(α)3(z)=16(α+3)(α+2)(α+1)+(1)2(α+3)(α+2)z+(α+3)2z2+(1)6z3 を計算してみた。 また、4次5次の多項式も L(α)4(z)=124(α+4)(α+3)(α+2)(α+1)+(1)6(α+4)(α+3)(α+2)z+14(α+4)(α+3)z2+(1)6(α+4)z3+124z4L(α)5(z)=1120(α+5)(α+4)(α+3)(α+2)(α+1)+(1)24(α+5)(α+4)(α+3)(α+2)z+112(α+5)(α+4)(α+3)z2+(1)12(α+5)(α+4)z3+124(α+5)z4+(1)120z5 となる。
この公式の利点は漸化式にくらべて各項の係数を把握しやすいというところだと思う。もちろん低次の項に限っての話でしょうが。 こららの他にクンマー級数(クンマーの超幾何関数)をもちいて L(α)n(z)=(α+n)!n!α!M(n,α+1;z) とかくこともできる。この公式をformula Kとすると、もちろんformula GLで計算される結果と一致する。 たとえばα=1.5,N=4α=5.0,N=7の場合を計算しグラフに示すと
といったようになる。

またformula GLを使うことで L(α)n(z)z=L(α+1)n1(z) という公式も証明できる。

単純に、L(α)n(z)の微分を考えてやると、 L(α)n(z)z=nk=1(1)k(k1)!(nk)!zk1(α+n)!(α+k)!=nk=1(1)k1(k1)!((n1)(k1))!zk1(α+1+(n1))!(α+1+(k1))! となる。ここでν=k1 としてやると L(α)n(z)z=n1ν=0(1)νν!(n1ν)!zν(α+1+(n1))!(α+1+ν)!=L(α+1)n1(z) となることがわかる。

0 件のコメント:

コメントを投稿