【組合數學】遞推方程 ( 有重根遞推方程求解問題 | 問題提出 )

韓曙亮發表於2020-10-24





一、有重根遞推方程求解問題



有些 遞推方程特徵方程特徵根重根 的情況 , 特徵方程解出來的 特徵根有一部分是相等的 , 這樣就使得 通解中的常數無法獲取唯一的值 ;


參考 : 【組合數學】遞推方程 ( 通解定義 | 無重根下遞推方程通解結構定理 ) 二、無重根下遞推方程通解結構定理

在 “無重根下遞推方程通解結構定理” 章節中 , 通解要求 方程組中的 係數行列式不等於 0 0 0 , ∏ 1 ≤ i < j ≤ k ( q i − q k ) ≠ 0 \prod\limits_{1 \leq i < j \leq k} ( q_i - q_k ) \not= 0 1i<jk(qiqk)=0 , 如果有兩個特徵根 q i , q k q_i , q_k qi,qk 相等 , 則上面的 "係數行列式不等於 0 0 0" 便無法實現 ;


如果特徵方程有重根 , 就不能使用 “無重根下遞推方程公式求法” 進行遞推方程的求解 ;


針對有重根的遞推方程 , 需要將其 線性無關的元素 都找到 , 線性組合在一起 , 才能得到通解 ;


線性組合 : 將一個解乘以 c 1 c_1 c1 , 另一個解乘以 c 2 c_2 c2 , 相加之後的組合 ;





二、有重根遞推方程示例



遞推方程 : H ( n ) − 4 H ( n − 1 ) + 4 H ( n − 2 ) = 0 H(n) - 4H(n-1) + 4H(n-2) = 0 H(n)4H(n1)+4H(n2)=0

初值 : H ( 0 ) = 0 , H ( 1 ) = 1 H(0) = 0 , H(1) = 1 H(0)=0,H(1)=1



無重根下遞推方程求解完整過程 :

  • 1 . 寫出特徵方程 :
    • ( 1 ) 遞推方程標準形式 : 寫出遞推方程 標準形式 , 所有項都在等號左邊 , 右邊是 0 0 0 ;
    • ( 2 ) 特徵方程項數 : 確定 特徵方程項數 , 與 遞推方程項數相同 ;
    • ( 3 ) 特徵方程次冪數 : 最高次冪是 特徵方程項數 − 1 -1 1 , 最低次冪 0 0 0 ;
    • ( 4 ) 寫出 沒有係數 的特徵方程 ;
    • ( 5 ) 逐位將遞推方程的係數 抄寫 到特徵方程中 ;
  • 2 . 解特徵根 :特徵方程的特徵根解出來 , x = − b ± b 2 − 4 a c 2 a x = \cfrac{-b \pm \sqrt{b^2 - 4ac}}{2a} x=2ab±b24ac
  • 3 . 構造遞推方程的通解 : 構造 c 1 q 1 n + c 2 q 2 n + ⋯ + c k q k n c_1q_1^n + c_2q_2^n + \cdots + c_kq_k^n c1q1n+c2q2n++ckqkn 形式的線性組合 , 該線性組合就是遞推方程的解 ;
  • 4 . 求通解中的常數 : 將遞推方程初值代入通解 , 得到 k k k k k k 元方程組 , 通過解該方程組 , 得到通解中的常數 ;
    • ( 1 ) 常數代入通解 : 得到最終的遞推方程的解 ;


      遞推方程 -> 特徵方程 -> 特徵根 -> 通解 -> 代入初值求通解常數


根據上述求解過程進行求解 :


1 . 特徵方程 :

( 1 ) 遞推方程標準形式 : 遞推方程已經是標準形式 ;

( 2 ) 特徵方程項數 : 與 遞推方程項數 相同 , 3 3 3 項 ;

( 3 ) 特徵方程次冪數 : 最高次冪是 特徵方程項數減一 , 3 − 1 = 2 3-1=2 31=2 , 最低次冪 0 0 0 ;

( 4 ) 寫出 沒有係數 的特徵方程 : x 2 + x + 1 = 0 x^2 + x + 1 = 0 x2+x+1=0

( 5 ) 逐位將遞推方程的係數 抄寫 到特徵方程中 ;

1 x 2 + ( − 4 ) x + ( 4 ) 1 = 0 1x^2 + (-4)x + (4)1 = 0 1x2+(4)x+(4)1=0

x 2 − 4 x + 4 = 0 x^2 - 4x + 4 = 0 x24x+4=0


2 . 解特徵根 :特徵方程的特徵根解出來 , x = − b ± b 2 − 4 a c 2 a x = \cfrac{-b \pm \sqrt{b^2 - 4ac}}{2a} x=2ab±b24ac

x = 4 ± 16 − 16 2 = 2 x=\cfrac{4 \pm \sqrt{16 - 16}}{2} = 2 x=24±1616 =2

兩個特徵根都是 2 2 2 , q 1 = 2 , q 2 = 2 q_1=2, q_2 = 2 q1=2,q2=2 ;


3 . 構造遞推方程的通解 : 構造 c 1 q 1 n + c 2 q 2 n + ⋯ + c k q k n c_1q_1^n + c_2q_2^n + \cdots + c_kq_k^n c1q1n+c2q2n++ckqkn 形式的線性組合 , 該線性組合就是遞推方程的解 ;

通解是 : H ( n ) = c 1 2 n + c 2 2 n = c 2 n H(n) = c_12^n + c_22^n = c2^n H(n)=c12n+c22n=c2n


4 . 求通解中的常數 : 將遞推方程初值代入通解 , 得到 k k k k k k 元方程組 , 通過解該方程組 , 得到通解中的常數 ;

c 2 n c2^n c2n 代入到 x 2 − 4 x + 4 = 0 x^2 - 4x + 4 = 0 x24x+4=0 特徵方程中 , c c c 是無解的 ;



如果 兩個特徵根 都是 2 2 2 , 線性相關 , 此時就 無法確定通解中的 c 1 , c 2 c_1, c_2 c1,c2 待定常數 ;

觀察 n 2 n n2^n n2n 是解 , 該解與 2 n 2^n 2n 線性無關 , 將上述兩個解進行線性組合 ,

c 1 n 2 n + c 2 2 n c_1n2^n + c_22^n c1n2n+c22n 線性組合 , 是遞推方程的解 ,

將初值代入 , 可以解出 c 1 , c 2 c_1, c_2 c1,c2 常數的值 ;

相關文章