させて∑ k = 0 r p k ( n ) y ( n + k ) = f ( n ) {\textstyle \sum _{k=0}^{r}p_{k}(n)\,y(n+k)=f(n)} または同等にL y = f {\textstyle Ly=f} これは多項式係数を持つ漸化式です。この方程式の解を計算するアルゴリズムはいくつか存在します。これらのアルゴリズムは、多項式解、有理式解、超幾何式解、およびダランベール式解を計算できます。同次方程式の解は、線形漸化演算子の核によって与えられます。 カー L = { y ∈ K N : L y = 0 } {\textstyle \ker L=\{y\in \mathbb {K} ^{\mathbb {N} }\,:\,Ly=0\}} 数列空間の部分空間として、このカーネルは基底 を持つ。[ 1 ] とする{ y ( 1 ) 、 y ( 2 ) 、 … 、 y ( m ) } {\textstyle \{y^{(1)},y^{(2)},\dots ,y^{(m)}\}} 基礎となるカー L {\textstyle \ker L} すると形式的な和c 1 y ( 1 ) + ⋯ + c m y ( m ) {\textstyle c_{1}y^{(1)}+\dots +c_{m}y^{(m)}} 任意の定数に対してc 1 、 … 、 c m ∈ K {\textstyle c_{1},\dots ,c_{m}\in \mathbb {K} } 同次問題の一般解と呼ばれる。L y = 0 {\textstyle Ly=0} 。 もしy ~ \textstyle {\tilde {y}}} は、L y = f {\textstyle Ly=f} つまりL y ~ = f {\textstyle L{\tilde {y}}=f} 、 それからc 1 y ( 1 ) + ⋯ + c m y ( m ) + y ~ {\textstyle c_{1}y^{(1)}+\dots +c_{m}y^{(m)}+{\チルダ {y}}} これも非同次問題の解であり、非同次問題の一般解と呼ばれます。
多項式解 1980年代後半、セルゲイ・A・アブラモフは、漸化式の一般多項式解を求めるアルゴリズムを記述した。y ( n ) ∈ K [ n ] {\textstyle y(n)\in \mathbb {K} [n]} 右辺が多項式であるf ( n ) ∈ K [ n ] {\textstyle f(n)\in \mathbb {K} [n]} 彼(そして数年後、マルコ・ペトコフシェク )は多項式解の次数上限を与えた。この方法では、線形方程式系を 考えることで簡単に問題を解くことができる。[ 2 ] [ 3 ] [ 4 ] 1995年、アブラモフ、ブロンシュタイン、ペトコフシェクは、多項式の場合、特定の冪基底(つまり通常の基底ではない)での漸化式の冪級数 解を考えることで、より効率的に解けることを示した。( x n ) n ∈ N {\textstyle (x^{n})_{n\in \mathbb {N} }} [ 5 ]
より一般的な解(例えば、有理数解や超幾何数解)を見つけるための他のアルゴリズムも、多項式解を計算するアルゴリズムに依存している。
超幾何解 シーケンスy ( n ) {\textstyle y(n)} 連続する 2 つの項の比が有理関数である場合、は超幾何関数 と呼ばれます。n {\displaystyle n} つまりy ( n + 1 ) / y ( n ) ∈ K ( n ) {\textstyle y(n+1)/y(n)\in \mathbb {K} (n)} これは、数列が多項式係数を持つ1階漸化式の解である場合に限ります。超幾何数列の集合は加法に関して閉じていないため、数列空間の部分空間ではありません。
1992年、マルコ・ペトコフシェク は、 右辺がf {\displaystyle f} は超幾何数列の和です。このアルゴリズムは、有理関数のゴスパー・ペトコフシェク標準形を利用します。この特定の表現では、変換された方程式の多項式解を再び考慮するだけで十分です。[ 3 ]
マーク・ファン・ホーイによる、より効率的な別のアプローチがある。最初の係数多項式と最後の係数多項式の根を考慮すると、p 0 {\textstyle p_{0}} そしてp r {\textstyle p_{r}} 特異点と呼ばれるものを利用して、超幾何数列はy ( n ) {\textstyle y(n)} 形式を表すy ( n ) = c r ( n ) z n Γ ( n − ξ 1 ) e 1 Γ ( n − ξ 2 ) e 2 ⋯ Γ ( n − ξ s ) e s {\displaystyle y(n)=c\,r(n)\,z^{n}\,\Gamma (n-\xi _{1})^{e_{1}}\Gamma (n-\xi _{2})^{e_{2}}\cdots \Gamma (n-\xi _{s})^{e_{s}}} 一部の人にとってc ∈ K 、 z ∈ K ¯ 、 s ∈ N 、 r ( n ) ∈ K ¯ ( n ) 、 ξ 1 、 … 、 ξ s ∈ K ¯ {\textstyle c\in \mathbb {K} ,z\in {\overline {\mathbb {K} }},s\in \mathbb {N} ,r(n)\in {\overline {\mathbb {K} }}(n),\xi _{1},\dots ,\xi _{s}\in {\overline {\mathbb {K} }}} とξ 私 − ξ j ∉ Z {\textstyle \xi _{i}-\xi _{j}\notin \mathbb {Z} } のために私 ≠ j {\textstyle i\neq j} そしてe 1 、 … 、 e s ∈ Z {\textstyle e_{1},\dots ,e_{s}\in \mathbb {Z} } 。 ここΓ ( n ) \textstyle \Gamma (n) はガンマ関数 を表し、K ¯ {\textstyle {\overline {\mathbb {K} }}} 体の代数的閉包 K {\textstyle \mathbb {K} } それからξ 1 、 … 、 ξ s {\textstyle \xi _{1},\dots ,\xi _{s}} は方程式の特異点(つまり、p 0 {\textstyle p_{0}} またはp r {\textstyle p_{r}} さらに、指数の上限を計算することもできます。e 私 {\textstyle e_{i}} 固定値の場合ξ 1 、 … 、 ξ s 、 e 1 、 … 、 e s {\textstyle \xi _{1},\dots ,\xi _{s},e_{1},\dots ,e_{s}} 候補を与える仮説を立てることが可能である。z {\textstyle z} 特定のz {\textstyle z} 有理関数を得るための仮説を再び立てることができるr ( n ) {\textstyle r(n)} アブラモフのアルゴリズムにより、すべての可能性を考慮すると、漸化式の一般解が得られます。[ 7 ] [ 8 ]
例
内転 内反 の数y ( n ) {\textstyle y(n)} セットのn {\textstyle n} 要素は漸化式で与えられるy ( n ) = ( n − 1 ) y ( n − 2 ) + y ( n − 1 ) 。 {\displaystyle y(n)=(n-1)\,y(n-2)+y(n-1).} 例えばペトコフシェクのアルゴリズム を適用すると、この漸化式には多項式、有理数、超幾何式の解が存在しないことがわかる。[ 4 ]
アプリケーション 関数F ( n 、 k ) {\textstyle F(n,k)} は、以下の条件を満たす場合に超幾何級数と呼ばれます。F ( n 、 k + 1 ) / F ( n 、 k ) 、 F ( n + 1 、 k ) / F ( n 、 k ) ∈ K ( n 、 k ) {\textstyle F(n,k+1)/F(n,k),F(n+1,k)/F(n,k)\in \mathbb {K} (n,k)} どこK ( n 、 k ) {\textstyle \mathbb {K} (n,k)} は有理関数を表す。n {\textstyle n} そしてk {\textstyle k} 超幾何和は、次の形式の有限和である。f ( n ) = ∑ k F ( n 、 k ) {\textstyle f(n)=\sum _{k}F(n,k)} どこF ( n 、 k ) {\textstyle F(n,k)} は超幾何級数です。ツァイルベルガー の独創的なテレスコープアルゴリズムは、このような超幾何級数の和を多項式係数を持つ漸化式に変換できます。この式を解くと、例えば、の閉形式解と呼ばれる超幾何級数の解の線形結合が得られます。f {\textstyle f} [ 4 ]
参考文献 ↑ 数列がほぼすべての項で等しい場合に等しいとみなすならば、この基底は有限である。これについては、Petkovšek、Wilf、Zeilberger共著の書籍『A=B』を参照されたい。 ↑ アブラモフ、セルゲイ A. (1989). 「線形微分方程式および差分方程式の多項式解の 探索に関連するコンピュータ代数の問題」。モスクワ大学計算数学およびサイバネティクス 。3 。 1 2 Petkovšek, Marko (1992). "多項式係数を持つ線形漸化式の超幾何解". Journal of Symbolic Computation . 14 ( 2– 3): 243– 264. doi : 10.1016/0747-7171(92)90038-6 . ISSN 0747-7171 . 1 2 3 4 ペトコフシェク、マルコ。ウィルフ、ハーバート S.ツァイルベルガー、ドロン (1996)。 A=B 。 AKピーターズ。 ISBN 978-1568810638 OCLC 33898705 ↑ Abramov, Sergei A.; Bronstein, Manuel; Petkovšek, Marko (1995). "線形作用素方程式の多項式解について". Proceedings of the 1995 international symposium on Symbolic and algebraic computation - ISSAC '95 . ACM. pp. 290–296 . CiteSeerX 10.1.1.46.9373 . doi : 10.1145/220346.220384 . ISBN 978-0897916998 . S2CID 14963237 . ↑ Abramov, Sergei A. (1989). "多項式係数を持つ線形微分方程式および差分方程式の有理解". USSR Computational Mathematics and Mathematical Physics . 29 (6): 7– 12. doi : 10.1016/s0041-5553(89)80002-3 . ISSN 0041-5553 . ↑ van Hoeij, Mark (1999). "線形漸化式の有限特異点と超幾何解". Journal of Pure and Applied Algebra . 139 ( 1– 3): 109– 131. doi : 10.1016/s0022-4049(99)00008-0 . ISSN 0022-4049 . ↑ Cluzeau, Thomas; van Hoeij, Mark (2006). "線形漸化式の超幾何解の計算". Applicable Algebra in Engineering, Communication and Computing . 17 (2): 83– 115. doi : 10.1007/s00200-005-0192-x . ISSN 0938-1279 . S2CID 7496623 . ↑ Abramov, Sergei A.; Petkovšek, Marko (1994). "線形微分方程式および差分方程式のダランベール解". Proceedings of the international symposium on Symbolic and algebraic computation - ISSAC '94 . ACM. pp. 169–174 . doi : 10.1145/190347.190412 . ISBN 978-0897916387 . S2CID 2802734 . ↑ "A000165 - OEIS" . oeis.org . 2018-07-02 に取得.