漸化式ソルバー

漸化式ソルバー

漸化式ソルバーは、特性方程式を解くことで、任意の定数係数線形同次漸化式の一般項（閉じた形の解）を計算し、複素平面上に根をプロットし、数列の最初のN項を生成します。漸化式は、順序付けられた係数リストまたは a(n) = 3·a(n−1) − 2·a(n−2) のような自然な数式のいずれかで入力でき、異なる実数解、重解、共役複素数解を自動的に処理します。

線形漸化式とは？

次数 k の定数係数線形同次漸化式は、以下の形式を持ちます。

ここで、c₁, c₂, …, c_k は固定された実数であり、k は次数です。k 個の初期値 a(0), a(1), …, a(k−1) と合わせることで、漸化式はそれ以降のすべての項を一意に定義します。代表的な例は以下の通りです。

• フィボナッチ数列: a(n) = a(n−1) + a(n−2), 初期値 0, 1 → 0, 1, 1, 2, 3, 5, 8, 13, …
• リュカ数列: a(n) = a(n−1) + a(n−2), 初期値 2, 1 → 2, 1, 3, 4, 7, 11, 18, 29, …
• ペル数: a(n) = 2·a(n−1) + a(n−2), 初期値 0, 1 → 0, 1, 2, 5, 12, 29, 70, …
• トリボナッチ数列: a(n) = a(n−1) + a(n−2) + a(n−3), 初期値 0, 0, 1 → 0, 0, 1, 1, 2, 4, 7, 13, 24, …

特性方程式法

a(n) の一般項を求めるために、a(n) = rⁿ の形の解を探します。これを漸化式に代入し、r^n−k で割ると以下のようになります。

これは特性方程式と呼ばれ、r に関する k 次の多項式です。代数学の基本定理により、この方程式は（重複を含めて）ちょうど k 個の複素数解を持ちます。漸化式の一般解は、これらの解の構造によって決まります。

ケース1：異なる実数解 r₁, …, r_k

定数 A₁, …, A_k は、n = 0, 1, …, k−1 を代入し、初期値に対する線形連立方程式を解くことで決定されます。

ケース2：重複度 m の解 r

各重解は、m 個の線形独立な基底数列 rⁿ, n·rⁿ, n²·rⁿ, …, n^m−1·rⁿ を提供します。

ケース3：共役複素数解 r = ρ·e^iθ, r̄ = ρ·e^−iθ

漸化式が実数係数を持つ場合、複素数解は常に共役なペアで現れます。各ペアは、幾何学的なエンベロープ ρⁿ と周波数 θ を持つ実数の振動項に統合されます。

最大根による増加の分類

ρ = max|r_i| を解の最大の絶対値（スペクトル半径）とします。a(n) の長期的な挙動は以下によって支配されます。

フィボナッチ数列の計算例

フィボナッチ漸化式 a(n) = a(n−1) + a(n−2)、初期値 a(0) = 0, a(1) = 1 を考えます。

1. 特性方程式: r² − r − 1 = 0
2. 解（解の公式）: r = (1 ± √5) / 2、すなわち φ ≈ 1.6180 および ψ ≈ −0.6180
3. 一般形: a(n) = A·φⁿ + B·ψⁿ
4. 初期条件の適用: A + B = 0 および A·φ + B·ψ = 1 より、A = 1/√5, B = −1/√5
5. ビネの公式: a(n) = (φⁿ − ψⁿ) / √5

|ψ| < 1 であるため、n → ∞ において第2項は消失します。したがって、a(n) は近似的に φⁿ / √5 となり、これがフィボナッチ数が各ステップで約 φ 倍ずつ成長する理由です。

この電卓の使い方

1. 入力モードを選択: 「ガイド付き」では次数を選択し、カンマ区切りの係数を入力できます。「自由形式の数式」では a(n) = a(n-1) + 6*a(n-2) - 8*a(n-3) のような完全な漸化式を入力できます。
2. 係数または数式を入力: 小数 (0.5) と分数 (1/2) の両方が受け付けられます。
3. 初期値を入力: 漸化式の次数に合わせ、正確に k 個の値（a(0), a(1), …, a(k−1)）を入力してください。
4. 表示する項数を選択: 最大 60 項まで指定できます。
5. 「解く」をクリック: 結果ページに、特性方程式、複素平面上の解の位置、一般項の公式、数列のアニメーション棒グラフが表示されます。

対応しているケースと制限事項

• 次数: 1から6まで（3次以上の特性多項式は Durand–Kerner 法により数値的に解かれます）。
• 実定数係数: 複素数係数には対応していません。c_i は実数である必要があります。
• 同次式のみ: このツールは同次漸化式を解きます（+ n や + 2ⁿ のような強制項は含まれません）。非同次漸化式の場合は、ここで同次部分を解き、別途特解を加えてください。
• 数値精度: 計算は IEEE-754 倍精度で行われます。非常に条件の悪い漸化式（解の大きさの差が非常に大きい場合）では、検証バナーが一般項と反復値の乖離を警告することがあります。

応用分野

• アルゴリズム解析: 分割統治アルゴリズムの実行時間はしばしば線形漸化式に従います（マスター定理）。
• 組合せ数学: カタラン数、攪乱順列、タイル張りなどの数え上げ数列は、頻繁に漸化式で与えられます。
• 信号処理: フィードバックを持つ離散時間 LTI システムは線形漸化式です。その安定性は解の位置（単位円内 ⇔ 安定）によって決定されます。
• 個体群動態と金融: 複利計算、年齢構造化人口モデル、自己回帰 AR(p) 時系列。
• 物理学: 1次元格子モデル、密結合ハミルトニアン、トランスファー行列法。

よくある質問

定数係数線形漸化式とは何ですか？

定数係数線形漸化式とは、a(n) = c₁·a(n−1) + c₂·a(n−2) + … + c_k·a(n−k) の形式の方程式です。ここで c₁, c₂, …, c_k は固定された実数であり、k は次数です。数列の各項は、直前の k 個の項の線形結合です。一般的な例には、フィボナッチ漸化式や、異なる初期値を持つリュカ漸化式があります。

漸化式の特性方程式とは何ですか？

漸化式 a(n) = c₁·a(n−1) + c₂·a(n−2) + … + c_k·a(n−k) に対し、その特性方程式は r^k − c₁·r^k−1 − c₂·r^k−2 − … − c_k = 0 です。この多項式方程式は（重複を含めて）ちょうど k 個の複素数解を持ち、漸化式のすべての解は、r を解、j をその重複度マイナス 1 までの指数としたときの n^j·rⁿ の形式の数列の線形結合となります。

a(n) の一般項を求めるにはどうすればよいですか？

特性方程式を解いて、その根 r₁, r₂, …, r_k を求めます。すべての根が異なる場合、一般項は a(n) = A₁·r₁ⁿ + A₂·r₂ⁿ + … + A_k·r_kⁿ となります。定数 A_i は初期値を代入して連立方程式を解くことで求められます。根 r が重複度 m を持つ場合、rⁿ, n·rⁿ, n²·rⁿ, …, n^m−1·rⁿ の m 個の基底項が寄与します。この電卓は、この全手順を自動的に行います。

複素数解は数列にとって何を意味しますか？

実数係数の漸化式では、複素数解は常に共役なペア r = ρ·e^iθ および r̄ = ρ·e^−iθ として現れます。このペアは振動挙動を引き起こし、一般項には 2·ρⁿ·[α·cos(nθ) − β·sin(nθ)] の形の項が含まれます。ρ = 1 のとき数列は一定振幅で振動し、ρ < 1 のとき振動は減衰し、ρ > 1 のとき振幅は幾何学的に増加します。

なぜ最大根が数列の増加を決定するのですか？

n が大きくなると、|r| が最大の項が他の項よりも速く成長するため、その項が支配的になります。ρ = max|r_i| とすると、|a(n)| は漸近的に ρⁿ に比例します（最大根が重解の場合は多項式因子がつきます）。このソルバーは、ρ < 1 なら 0 に収束、ρ = 1 なら有界、ρ > 1 なら幾何学的増加として数列を分類します。

このツールでフィボナッチ数列を解くことはできますか？

はい。漸化式 a(n) = a(n−1) + a(n−2) と初期値 0, 1 を入力してください。電卓は特性方程式 r² − r − 1 = 0 から根 φ = (1 + √5)/2 と ψ = (1 − √5)/2 を導き出し、ビネの公式 a(n) = (φⁿ − ψⁿ) / √5 を返します。入力フォームの上のフィボナッチの例をクリックすると、完全な解答例が表示されます。

a(n) = a(n−1) + n のような非同次漸化式には対応していますか？

いいえ、このツールは同次漸化式（強制項がないもの）のみを解きます。非同次漸化式の場合、一般解を「同次部分の解（ここで計算可能）」と「強制項に合わせた特解」の和として分解して解く必要があります。一般的な特解の仮定（アンザッツ）には、多項式項には同次数の多項式、指数項には C·rⁿ、三角関数項には A·cos(nθ) + B·sin(nθ) などがあります。

参考文献

• 漸化式 — Wikipedia
• 線形漸化式 — Wikipedia
• 特性方程式 — Wikipedia
• フィボナッチ数 — Wikipedia
• デュラン＝ケルナー法 — Wikipedia

シーケンスツール:

• 算術シーケンス電卓 高精度
• キューブ番号リスト
• 最初のn個の素数
• 等比数列電卓
• フィボナッチ数のリスト
• 素数のリスト
• 平方数リスト おすすめ
• コラッツ予想電卓
• ハッピーナンバー電卓
• 魔方陣ジェネレーター おすすめ
• カタラン数ジェネレーター
• シグマ記法電卓 総和 新しい
• 積の記号電卓 (Π パイ記法) 新しい
• パスカルの三角形ジェネレーター 新しい
• 双子素数ファインダー 新しい
• 分割数電卓 新しい
• 漸化式ソルバー 新しい