数列的递推公式求通项(特征方程)

本文详细介绍了线性常系数递推关系的解法,包括齐次与非齐次递推关系,并通过实例演示了如何利用特征方程求解递推序列。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

如: a n + c 1 a n − 1 + c 2 a n − 2 + ⋯ c k a n − k = b n a_{n}+c_{1}a_{n-1}+c_{2}a_{n-2}+\cdots c_{k}a_{n-k}=b_{n} an+c1an1+c2an2+ckank=bn(n>k-1)
( c 1 , c 2 ⋯ c k c_{1},c_{2}\cdots c_{k} c1,c2ck为常数,k为正整数)的递推关系为k阶线性常系数递推关系
a 0 = d 0 , a 1 = d 1 , ⋯   , a k − 1 = d k − 1 a_{0}=d_{0},a_{1}=d_{1},\cdots ,a_{k-1}=d_{k-1} a0=d0,a1=d1,,ak1=dk1称为初值条件

1. 当 b n b_{n} bn=0时##

b n b_{n} bn=0,时 a n + c 1 a n − 1 + c 2 a n − 2 + ⋯ c k a n − k = 0 a_{n}+c_{1}a_{n-1}+c_{2}a_{n-2}+\cdots c_{k}a_{n-k}=0 an+c1an1+c2an2+ckank=0称为k阶线性常系数齐次递推关系
如斐波那契数列 F n − F n − 1 − F n − 2 = 0 F_{n}-F_{n-1}-F_{n-2}=0 FnFn1Fn2=0(n>2)称为二阶线性常系数齐次递推关系

引入特征方程的定义
给定 a n + c 1 a n − 1 + c 2 a n − 2 + ⋯ c k a n − k = 0 a_{n}+c_{1}a_{n-1}+c_{2}a_{n-2}+\cdots c_{k}a_{n-k}=0 an+c1an1+c2an2+ckank=0
C ( x ) = x n + c 1 x n − 1 + c 2 x n − 2 + ⋯ + c k C(x)=x^{n}+c_{1}x^{n-1}+c_{2}x^{n-2}+\cdots +c_{k} C(x)=xn+c1xn1+c2xn2++ck称为性常系数齐次递推关系的特征多项式,而称C(x)=0为特征方程,方程的根称为特征根
如斐波那契数列 F n − F n − 1 − F n − 2 = 0 F_{n}-F_{n-1}-F_{n-2}=0 FnFn1Fn2=0,(n>2)的特征方程为 x 2 − x − 1 = 0 x^{2}-x-1=0 x2x1=0

1. 为实数根时
a n − a n − 1 − 12 a n − 2 = 0 , a 0 = 3 , a 1 = 26 a_{n}-a_{n-1}-12a_{n-2}=0,a_{0}=3,a_{1}=26 anan112an2=0,a0=3,a1=26的解
: 特征方程 x^{2}-x-1=0
根为: x=4 x=-3
故通解为 a n = A 1 4 n + A 2 ( − 3 ) n , n ≥ 0 a_{n}=A_{1}4^{n}+A_{2}(-3)^{n} ,n\geq 0 an=A14n+A2(3)n,n0
由初值条件 3 = A 1 + A 2 3=A_{1}+A_{2} 3=A1+A2 26 = A 1 4 + A 2 ( − 3 ) 26=A_{1}4+A_{2}(-3) 26=A14+A2(3)
解得 A 1 = 5 A_{1}=5 A1=5 A 2 = − 2 A_{2}=-2 A2=2
所求得解 a n = 5 ⋅ 4 n − 2 ( − 3 ) n a_{n}=5 \cdot 4^{n}-2(-3)^{n} an=54n2(3)n

2. 有共轭复根时
设 $ \alpha _{1}$和 $ \alpha _{2}$是一对共轭复根则有
$ \alpha _{1}=\rho (cos\theta +isin\theta)$ α 2 = α 1 ˉ = ρ ( c o s θ − i s i n θ ) \alpha _{2}=\bar{\alpha _{1}}=\rho (cos\theta -isin\theta) α2=α1ˉ=ρ(cosθisinθ)
A 1 α 1 n + A 2 α 2 n A_{1}\alpha _{1}^{n}+A_{2}\alpha _{2}^{n} A1α1n+A2α2n
= A 1 ρ n ( c o s θ + i s i n θ ) n + A 2 ρ n ( c o s θ − i s i n θ ) n A_{1}\rho^{n}(cos\theta +isin\theta)^{n}+A_{2}\rho^{n} (cos\theta -isin\theta)^{n} A1ρn(cosθ+isinθ)n+A2ρn(cosθisinθ)n
= A 1 ρ n ( c o s n θ + i s i n n θ ) + A 2 ρ n ( c o s n θ − i s i n n θ ) A_{1}\rho^{n}(cosn\theta +isinn\theta)+A_{2}\rho^{n} (cosn\theta -isinn\theta) A1ρn(cosnθ+isinnθ)+A2ρn(cosnθisinnθ)
= ( A 1 + A 2 ) ρ n c o s n θ + i ( A 1 − A 2 ) ρ n s i n n θ (A_{1}+A_{2})\rho^{n}cosn\theta +i(A_{1}-A_{2})\rho^{n}sinn\theta (A1+A2)ρncosnθ+i(A1A2)ρnsinnθ
= A ρ n c o s n θ + B ρ n s i n n θ A\rho^{n}cosn\theta +B\rho^{n}sinn\theta Aρncosnθ+Bρnsinnθ
其中 A = A 1 + A 2 A=A_{1}+A_{2} A=A1+A2 B = i ( A 1 − A 2 ) B=i(A_{1}-A_{2}) B=i(A1A2)
计算时,可先求出各对共轭复根,再求A,B

例如 求 a n − a n − 1 + a n − 2 = 0 , a 0 = 1 , a 1 = 1 a_{n}-a_{n-1}+a_{n-2}=0,a_{0}=1,a_{1}=1 anan1+an2=0,a0=1,a1=1的解
: 特征方程 α 2 − α + 1 = 0 \alpha^{2}-\alpha+1=0 α2α+1=0

根为: α = 1 2 ± 3 2 i = c o s π 3 + i s i n π 3 \alpha =\frac{1}{2}\pm \frac{\sqrt{3}}{2}i=cos \frac{\pi}{3}+isin\frac{\pi }{3} α=21±23 i=cos3π+isin3π

故通解为 a n = A 1 α 1 n + A 2 α 2 n = A c o s n π 3 + B s i n n π 3 a_{n}=A_{1}\alpha_{1}^{n}+A_{2}\alpha_{2}^{n}=Acos \frac{n\pi}{3}+Bsin\frac{n\pi }{3} an=A1α1n+A2α2n=Acos3nπ+Bsin3nπ

由初值条件 1=A 1 = A 1 2 + B 3 2 1=A\frac{1}{2}+B\frac{\sqrt{3}}{2} 1=A21+B23

解得 A=1 B = 3 3 B=\frac{\sqrt{3}}{3} B=33

所求解为 a n = c o s n π 3 + 3 3 s i n n π 3 a_{n}=cos \frac{n\pi}{3}+\frac{\sqrt{3}}{3}sin\frac{n\pi }{3} an=cos3nπ+33 sin3nπ

3. 有重根时

a n = ( B 0 + B 1 n + ⋯ + B k − 1 n k − 1 ) α n a_{n}=(B_{0}+B_{1}n+\cdots +B_{k-1}n^{k-1})\alpha ^{n} an=(B0+B1n++Bk1nk1)αn

例如求 a n − 2 a n − 1 + a n − 2 = 0 , a 1 = 2 , a 2 = 3 a_{n}-2a_{n-1}+a_{n-2}=0,a_{1}=2,a_{2}=3 an2an1+an2=0,a1=2,a2=3的解
:特征方程为 x^{2}-2x+1=0 解得x=1为二重根
故通解为 a n = ( A 1 + A 2 n ) ⋅ 1 n = A 1 + A 2 n a_{n}=(A_{1}+A_{2}n)\cdot 1^{n}=A_{1}+A_{2}n an=(A1+A2n)1n=A1+A2n
由初值条件 2 = A 1 + A 2 2=A_{1}+A_{2} 2=A1+A2 3 = A 1 + A 2 2 3=A_{1}+A_{2}2 3=A1+A22
解得 A 1 = 1 A_{1}=1 A1=1 A 2 = 1 A_{2}=1 A2=1
所求解为 a n = 1 + n a_{n}=1+n an=1+n, n ≥ 1 n\geq 1 n1

总之
若$\alpha $是特征方程的单根,则递推关系的解中含有项
a n = A a n a_{n}=Aa^{n} an=Aan
若$\alpha $是特征方程的k重根,则递推关系的解中含有项
a n = ( B 0 + B 1 n + ⋯ + B k − 1 n k − 1 ) α n a_{n}=(B_{0}+B_{1}n+\cdots +B_{k-1}n^{k-1})\alpha ^{n} an=(B0+B1n++Bk1nk1)αn
若$\alpha_{1},\alpha _{2} 是 一 对 k 重 共 轭 复 根 , 且 是一对k重共轭复根,且 k\alpha _{2}=\bar{\alpha _{1}}=\rho (cos\theta -isin\theta)$则递推关系的解中含有项
a n = ( B 0 + B 1 n + ⋯ + B k − 1 n k − 1 ) ρ n c o s n θ + i ( C 0 + C 1 n + ⋯ + C k − 1 n k − 1 ) ρ n s i n n θ a_{n}=(B_{0}+B_{1}n+\cdots +B_{k-1}n^{k-1})\rho^{n}cosn\theta+i(C_{0}+C_{1}n+\cdots +C_{k-1}n^{k-1})\rho^{n}sinn\theta an=(B0+B1n++Bk1nk1)ρncosnθ+i(C0+C1n++Ck1nk1)ρnsinnθ

2. 当 b n ≠ 0 b_{n}\neq 0 bn=0

b n ≠ 0 b_{n}\neq 0 bn=0,时 a n + c 1 a n − 1 + c 2 a n − 2 + ⋯ c k a n − k = b n a_{n}+c_{1}a_{n-1}+c_{2}a_{n-2}+\cdots c_{k}a_{n-k}=b_{n} an+c1an1+c2an2+ckank=bn称为k阶线性常系数非齐次递推关系
a n a_{n} an的解可表示为齐次关系的解郁非齐次关系的解之和
而齐次递推关系的解渴由特征方程特征根的方法求得,所以我们只要求解非齐次递推关系的一个特解即可
非齐次的解法和非齐次微分方程的通解解法较为类似可以结合二者进行理解
1. 待定系数法1
b n = r n P m ( n ) b_{n}=r^{n}P_{m}(n) bn=rnPm(n),其中 P m ( n ) P_{m}(n) Pm(n)为n的m次多项式(如a+bn为一次多项式),则可设非齐次递推关系的解为
β n ∗ = r n n k q m ( n ) \beta _{n}^{*}=r^{n}n^{k}q_{m}(n) βn=rnnkqm(n)
其中 q m ( n ) q_{m}(n) qm(n)为系数待定的 q m ( n ) = B m n m + B m − 1 n m − 1 + ⋯ + B 1 n 1 + B 0 q_{m}(n)=B_{m}n^{m}+B_{m-1}n^{m-1}+\cdots +B_{1}n^{1}+B_{0} qm(n)=Bmnm+Bm1nm1++B1n1+B0
其中k为r作为特征根的个数,即
\begin{matrix}
k=0 & r不是特征根\
k=k& r是k重特征根
\end{matrix}
如求解 a n − a n − 1 − 12 a n − 2 = 3 n , a 0 = 3 , a 1 = 26 a_{n}-a_{n-1}-12a_{n-2}=3^{n},a_{0}=3,a_{1}=26 anan112an2=3n,a0=3,a1=26的解
: 特征方程 x 2 − x − 12 = 0 x^{2}-x-12=0 x2x12=0
根为: x=4 x=-3
故齐次关系通解为 a n = A 1 4 n + A 2 ( − 3 ) n , n ≥ 0 a_{n}=A_{1}4^{n}+A_{2}(-3)^{n},n \geq 0 an=A14n+A2(3)n,n0
自由项为 b n = 3 n b_{n}=3^{n} bn=3n,所以可设特解为 a n ∗ = B 3 n a_{n}^{*}=B3^{n} an=B3n
那么代入递推关系式有 B 3 n − B 3 n − 1 − 12 B 3 n − 2 = 3 n B3^{n}-B3^{n-1}-12B3^{n-2}=3^{n} B3nB3n112B3n2=3n
解得 B = − 3 2 B=-\frac{3}{2} B=23
故非齐次关系通解为 a n = A 1 4 n + A 2 ( − 3 ) n − 3 2 3 n , n ≥ 0 a_{n}=A_{1}4^{n}+A_{2}(-3)^{n}-\frac{3}{2}3^{n},n \geq 0 an=A14n+A2(3)n233n,n0
由初值条件 3 = A 1 + A 2 − 3 2 3=A_{1}+A_{2}-\frac{3}{2} 3=A1+A223 26 = A 1 4 + A 2 ( − 3 ) − 3 2 ⋅ 3 26=A_{1}4+A_{2}(-3)-\frac{3}{2}\cdot3 26=A14+A2(3)233
解得 A 1 = 44 7 A_{1}=\frac{44}{7} A1=744 A 2 = − 25 14 A_{2}=-\frac{25}{14} A2=1425
所求解为
a n = 44 7 ⋅ 4 n − 25 14 ( − 3 ) n − 3 2 ⋅ 3 n , n ≥ 0 a_{n}=\frac{44}{7}\cdot 4^{n}-\frac{25}{14}(-3)^{n}-\frac{3}{2}\cdot 3^{n},n\geq 0 an=7444n1425(3)n233n,n0

再如求 a n + 3 a n − 1 − 10 a n − 2 = 2 n ( n + 5 ) 的 通 解 a_{n}+3a_{n-1}-10a_{n-2}=2^{n}(n+5)的通解 an+3an110an2=2n(n+5)
:特征方程 $ x^{2}+3x-10=0$
根为: x=-5 x=2
根据自由项:$ b_{n}=2^{2}(n+5)$,所以特解设为
a n ∗ = n 1 ( B 0 + B 1 n ) 2 n a_{n}^{*}=n^{1}(B_{0}+B_{1}n)2^{n} an=n1(B0+B1n)2n
代入递推关系式
解得
$ a_{n}=A_{1}2{n}+A_{2}(-5){n}-(\frac{87}{49}n+\frac{1}{7}n{2})2{n}$

2. 待定系数法2
b n = r n P m ( n ) s i n θ b_{n}=r^{n}P_{m}(n)sin\theta bn=rnPm(n)sinθ b n = r n P m ( n ) c o s θ b_{n}=r^{n}P_{m}(n)cos\theta bn=rnPm(n)cosθ,其中 P m ( n ) P_{m}(n) Pm(n)为n的m次多项式(如a+bn为一次多项式),则可设非齐次递推关系的解为
β n ∗ = r n n k [ q m ( n ) s i n n θ + h m ( n ) s i n n θ ] \beta _{n}^{*}=r^{n}n^{k}[q_{m}(n)sinn\theta+h_{m}(n)sinn\theta] βn=rnnk[qm(n)sinnθ+hm(n)sinnθ]
其中 q m ( n ) q_{m}(n) qm(n) h m ( n ) h_{m}(n) hm(n)为系数待定的
q m ( n ) = B m n m + B m − 1 n m − 1 + ⋯ + B 1 n 1 + B 0 q_{m}(n)=B_{m}n^{m}+B_{m-1}n^{m-1}+\cdots +B_{1}n^{1}+B_{0} qm(n)=Bmnm+Bm1nm1++B1n1+B0
h m ( n ) = C m n m + C m − 1 n m − 1 + ⋯ + C 1 n 1 + C 0 h_{m}(n)=C_{m}n^{m}+C_{m-1}n^{m-1}+\cdots +C_{1}n^{1}+C_{0} hm(n)=Cmnm+Cm1nm1++C1n1+C0
其中k为r作为特征根的个数,即
\begin{matrix}
k=0 & r(cos\theta+isin\theta)不是特征根\
k=k& r(cos\theta+isin\theta)是k重特征根
\end{matrix}

如求 a n − 3 a n − 1 + 2 a n − 2 = 3 s i n ( n π 2 ) , a 0 = 0 , a 1 = 1 a_{n}-3a_{n-1}+2a_{n-2}=3sin(n\frac{\pi}{2}),a_{0}=0,a_{1}=1 an3an1+2an2=3sin(n2π),a0=0,a1=1的解
: 特征方程 x 2 − 3 x + 2 = 0 x^{2}-3x+2=0 x23x+2=0
根为:x=1 x=2
故齐次关系的通解为 a n = A 1 1 n + A 2 2 2 , n ≥ 0 a_{n}=A_{1}1^n+A_{2}2^{2},n\geq0 an=A11n+A222,n0
由于自由项: $ b_{n}=3sin(n\frac{\pi}{2})$,所以特解为
a n ∗ = B c o s ( n π 2 ) + C s i n ( n π 2 ) a_{n}^{*}=Bcos(n\frac{\pi}{2})+Csin(n\frac{\pi}{2}) an=Bcos(n2π)+Csin(n2π)
代入递推关系式解得 B = − 9 10 B=-\frac{9}{10} B=109 C = − 3 10 C=-\frac{3}{10} C=103
故通解为
a n = A 1 + A 2 2 2 − 9 10 c o s ( n π 2 ) − 3 10 s i n ( n π 2 ) , n ≥ 0 a_{n}=A_{1}+A_{2}2^{2}-\frac{9}{10}cos(n\frac{\pi}{2})-\frac{3}{10}sin(n\frac{\pi}{2}),n\geq0 an=A1+A222109cos(n2π)103sin(n2π),n0
再由初值条件可得
a n = 15 10 − 6 10 2 2 − 9 10 c o s ( n π 2 ) − 3 10 s i n ( n π 2 ) , n ≥ 0 a_{n}=\frac{15}{10}-\frac{6}{10}2^{2}-\frac{9}{10}cos(n\frac{\pi}{2})-\frac{3}{10}sin(n\frac{\pi}{2}),n\geq0 an=101510622109cos(n2π)103sin(n2π),n0

3. 非齐次递推关系齐次化法
b n = r b_{n}=r bn=r,r为常数
a n + c 1 a n − 1 + c 2 a n − 2 + ⋯ c k a n − k = r a_{n}+c_{1}a_{n-1}+c_{2}a_{n-2}+\cdots c_{k}a_{n-k}=r an+c1an1+c2an2+ckank=r
a n − 1 + c 1 a n − 2 + c 2 a n − 3 + ⋯ c k a n − k − 1 = r a_{n-1}+c_{1}a_{n-2}+c_{2}a_{n-3}+\cdots c_{k}a_{n-k-1}=r an1+c1an2+c2an3+ckank1=r
两式相减则变为齐次递推关系用特征方程求解
b n = r n b_{n}=r^{n} bn=rn
a n + c 1 a n − 1 + c 2 a n − 2 + ⋯ c k a n − k = r n a_{n}+c_{1}a_{n-1}+c_{2}a_{n-2}+\cdots c_{k}a_{n-k}=r^{n} an+c1an1+c2an2+ckank=rn
a n − 1 + c 1 a n − 2 + c 2 a n − 3 + ⋯ c k a n − k − 1 = r n − 1 a_{n-1}+c_{1}a_{n-2}+c_{2}a_{n-3}+\cdots c_{k}a_{n-k-1}=r^{n-1} an1+c1an2+c2an3+ckank1=rn1
式子二乘以r,两式相减,变为齐次关系式,用特征方程求解
如求 a n − a n − 1 − 6 a n − 2 = 3 n , a 0 = 5 , a 1 = 2 a_{n}-a_{n-1}-6a_{n-2}=3^{n},a_{0}=5,a_{1}=2 anan16an2=3n,a0=5,a1=2的解
: 有 a n − a n − 1 − 6 a n − 2 = 3 n a_{n}-a_{n-1}-6a_{n-2}=3^{n} anan16an2=3n
a n − 1 − a n − 2 − 6 a n − 3 = 3 n − 1 a_{n-1}-a_{n-2}-6a_{n-3}=3^{n-1} an1an26an3=3n1
两式相减得
a n − 4 a n − 1 − 3 a n − 2 + 18 a n − 3 = 0 , a 2 = 41 a_{n}-4a_{n-1}-3a_{n-2}+18a_{n-3}=0,a_{2}=41 an4an13an2+18an3=0,a2=41
特征方程为 x{3}-x{2}-3x+18=0
解得 x=-2,x=3(为2重根)
a n = A ( − 2 ) n + ( B + C n ) 3 n a_{n}=A(-2)^{n}+(B+Cn)3^{n} an=A(2)n+(B+Cn)3n
由初值条件得
a n = 74 25 ( − 2 ) n + ( 51 25 + 3 5 n ) 3 n a_{n}=\frac{74}{25}(-2)^{n}+(\frac{51}{25}+\frac{3}{5}n)3^{n} an=2574(2)n+(2551+53n)3n

本文参考太原理工大学魏毅强教授组合数学课件

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值