之前归纳了一波数列,这里补充一个用于求形如 a n = p a n − 1 + q a n − 2 a_n=p\ a_{n-1}+q\ a_{n-2} an=p an−1+q an−2 的通项的方法。
(上面这样的递推式可被称作“二阶线性齐次常系数递推式”的特殊情况)
例题
有一数列 { a n } \{a_{n}\} {an}, a 1 = x a_1=x a1=x, a 2 = y a_2=y a2=y,且 ∀ n ≥ 3 \forall\ n \ge 3 ∀ n≥3,都有 a n = p a n − 1 + q a n − 2 a_n=pa_{n-1}+qa_{n-2} an=pan−1+qan−2。求该数列的通项公式。
解
不难想象这种线性递推式的答案都是指数形式的。
我们知道一阶线性递推式推出来一定会有一个指数项。那二阶递推式出来是不是会有两个指数项呢?
答案是肯定的,递推式一定能表示为
a
n
=
c
1
r
1
n
−
1
+
c
2
r
2
n
−
1
a_n=c_1r_1^{n-1}+c_2r_2^{n-1}
an=c1r1n−1+c2r2n−1。这里给结论给的比较突兀,严格证明限于编者能力暂时没法给出 ,但是读者可以想象:对递推式的两个项不断迭代,一个迭代到
a
1
a_1
a1,一个迭代到
a
2
a_2
a2,这样这两个项一定齐次,且次数与
n
n
n 同阶。
现在问题就是求描述中的常数。我们利用类似待定系数的思想先将该式代入递推式:
c 1 r 1 n + c 2 r 2 n = p c 1 r 1 n − 1 + p c 2 r 2 n − 1 + q c 1 r 1 n − 2 + q c 2 r 2 n − 2 c_1r_1^n+c_2r_2^n=pc_1r_1^{n-1}+pc_2r_2^{n-1}+qc_1r_1^{n-2}+qc_2r_2^{n-2} c1r1n+c2r2n=pc1r1n−1+pc2r2n−1+qc1r1n−2+qc2r2n−2
c 1 r 1 n − 2 ( r 1 2 − p r 1 − q ) + c 2 r 2 n − 2 ( r 2 2 − p r 2 − q ) = 0 c_1r_1^{n-2}(r_1^2-pr_1-q)+c_2r_2^{n-2}(r_2^2-pr_2-q)=0 c1r1n−2(r12−pr1−q)+c2r2n−2(r22−pr2−q)=0 、
要使原式成立,必然要求
r
1
2
−
p
r
1
−
q
=
0
r_1^2-pr_1-q=0
r12−pr1−q=0 且
r
2
2
−
p
r
2
−
q
=
0
r_2^2-pr_2-q=0
r22−pr2−q=0,即
r
1
、
r
2
r_1、r_2
r1、r2 为方程
x
2
−
p
x
−
q
=
0
x^2-px-q=0
x2−px−q=0
的两根。
我们把上述方程称作原递推式的 特征方程 。
然后我们又有 a 1 = c 1 + c 2 = x a_1=c_1+c_2=x a1=c1+c2=x, a 2 = c 1 r 1 + c 2 r 2 = y a_2=c_1r_1+c_2r_2=y a2=c1r1+c2r2=y,
此时 r 1 、 r 2 r_1、r_2 r1、r2 已知,所以 c 1 、 c 2 c_1、c_2 c1、c2 可求,这样我们就求完了。
补充说明
我们可以从矩阵的观点来看,由递推式子我们可以得到转移矩阵:
(图片来自 http://tieba.baidu.com/p/3820580861 @lgy_hello_life,侵权请联系笔者删除)
这个矩阵的特征多项式1即为
x
2
=
p
x
+
q
x^2=px+q
x2=px+q。