下载文档原格式
二阶常系数递推关系求解方法一、递推关系的定义与性质在数学中，递推关系是指通过递推公式来描述数列中各项之间的关系。
常系数递推关系是指递推关系中各项的系数都是常数。
设有一个序列 {an}，其中 n 表示序列中的项数。
如果序列满足递推关系 an = c1an-1+ c2an-2 + ... + ck an-k ，其中ci (1 ≤ i ≤ k) 为常数，那么我们称该序列满足一个 k 阶常系数递推关系。
常系数递推关系的性质：1. 齐次性：如果一个递推关系的非齐次项为0，即对于所有的 i，ci = 0，则该递推关系称为齐次线性递推关系。
2. 非齐次性：如果一个递推关系的非齐次项不为0，即存在一些 i，ci ≠ 0，则该递推关系称为非齐次线性递推关系。
3.初值条件：对于一个k阶线性递推关系，需要给出前k项的初值条件才能确定整个序列。
二、求解齐次线性递推关系的通解对于线性递推关系 an = c1an-1+ c2an-2 + ... + ck an-k ，其中ci (1 ≤ i ≤ k) 为常数，我们可以采用特征根法求解其通解。
1. 假设通解为an = λn ，将其代入递推关系，得到λ^n = c1λ^(n-1)+ c2λ^(n-2) + ... + ck λ^(n-k)2.将等式左边的λ^n移至等式右边，得到λ^n - c1λ^(n-1) - c2λ^(n-2) - ... - ck λ^(n-k) = 03.将该齐次方程转化为特征方程，即λ^k - c1λ^(k-1) - c2λ^(k-2) - ... - ck = 04.解特征方程，得到k个实数或复数根λ1,λ2,...,λk。
5.得到齐次线性递推关系的通解为an = A1λ1^n + A2λ2^n + ... + Akλk^n其中A1,A2,...,Ak为待定系数。
通过给定的初值条件，可以使用线性方程组求解方法来确定待定系数A1,A2,...,Ak。
三、求解非齐次线性递推关系的通解对于非齐次线性递推关系 an = c1an-1+ c2an-2 + ... + ck an-k + f(n)，其中 f(n) 为一个关于 n 的函数，我们可以采用常数变易法求解其通解。
零化多项式特征多项式最⼩多项式常系数线性齐次递推零化多项式/特征多项式/最⼩多项式/常系数线性齐次递推约定:I n是n阶单位矩阵，即主对⾓线是1的n阶矩阵⼀个矩阵A的|A|是A的⾏列式默认A是⼀个n×n的矩阵定义零化多项式:对于⼀个矩阵A，它的⼀个零化多项式f(λ)是满⾜f(A)=0的多项式，定义域包含矩阵最⼩多项式:次数最低的零化多项式特征多项式对于⼀个n阶的矩阵A，它的特征多项式p(λ)=|λI n−A|λ定义域不⽌是R，还可以是矩阵p(λ)是关于λ的⼀个不超过n+1次的多项式即p(λ)=∑n0a i x iCayley-Hamilton定理:矩阵的特征多项式也是它的零化多项式求解特征多项式带⼊n个数，求出得|xI n−A|,得到n个矩阵，通过⾼斯消元可以O(n3)地求出⾏列式然后可O(n2)拉格朗⽇插值求出原来的多项式，总复杂度受限于⾼斯消元，为O(n4)求解最⼩多项式构造矩阵序列a i=A i求出它的⼀个线性递推r i，即m∑j=0r j a i−j=m∑j=0r j A i−j=(m∑j=0r m−j A j)⋅A i−m=0∴m∑j=0r m−j A j=0所以可以由r i翻转得到f(λ)求解a i前n项的复杂度受限于矩阵乘法为O(n4)，求解递推式的复杂度为O(n3)考虑到实际求解递推式时，随机⽣成了两个向量u,v实际是计算标量序列{uA i v}的递推式，所以实际每次求出uA i复杂度应为O(n2)求这个递推式需要⽤到a i前2n项，求解复杂度为O(n3)因此总复杂度为O(n3)(但是如果只是求出来并没有什么⽤，因为求解⽅法是随机的，甚⾄连检查⼀次保证正确都需要O(n2(n+e))的时间(e为矩阵⾮0位置个数))求解稀疏⽅程组设⽅程系数⽤矩阵A表⽰，右侧每个⽅程的常数⽤向量b表⽰，答案⽤向量x表⽰，则满⾜关系式Ax=b，即x=A−1b求出{A i b}线性递推式，反推出A−1b即可反推⽅法:带⼊线性递推的m项，则∑m i=0A m−i b⋅r i=0A m−i br i=0两边同乘A−1，得到A−1b⋅r m+∑m−1i=0求解矩阵k次幂我们要求解A k，常规做法是直接⽤快速幂设矩阵A的⼀个零化多项式是f(λ)显然，A k可以⽤⼀个多项式表⽰A k=∑k0w i A i{w i}构成了⼀个k+1次多项式F k(x)存在⼀种合法的表⽰是F k(x)=x k∵f(A)=0∴∀i,f(A)A i=0也就是相当于我们要求出x k对于f(x)这个n+1多项式取模显然可以通过类似快速幂的⽅式倍增求解这个多项式，每次对f(x)取模复杂度是O(n log n)就能在O(n log m log n)时间得求出F(x)最后得到的F(x)是⼀个n次多项式那么带⼊就可以快速求出A k可以认为这个复杂度是受限于求解A0,A1,⋯,A n−1的O(n4)对于元矩阵A为稀疏矩阵的情况，设其包含e个⾮零位置那么求解B⋅A的过程是O(n⋅e)的，求解A0,A1,⋯,A n−1的过程，是O(n2e)的求解零化多项式的复杂度也是O(n2(n+e))的，因此总复杂度为O(n2(n+e))⽽⼀般的矩阵快速幂是O(n3log k)的，这种⽅法适⽤情况⾮常特殊另外，对于并不需要知道整个矩阵的答案，并且A0,A1,⋯,A n−1特殊的具体问题，这个⽅法也⼗分有效求解常系数线性齐次递推问题是要求数列f i=∑n j=1a j⋅f i−j给出f0,f1,⋯,f n−1，求第k项的值线性递推显然可以⽤初始向量列与转移矩阵的幂次的乘积表⽰，即f i=(S⋅A i)n，其中A为转移矩阵，S为初始向量列，我们求的是第n项对于n=4的情况，我们的转移矩阵A是12341a 421a 331a 241a 1鉴于它的特殊性，我们可以直接求出它的特征多项式表达式由λI n −A =12341λ−a 42−1λ−a 33−1λ−a 24−1λ−a 1带⼊⾏列式最暴⼒的求法枚举⼀个排列p i ，设排列p 的逆序对为f (p )，|A |=∑(−1)f (p )ΠA i ,pi 实际上合法的排列只有n 个，就是枚举p i =n那么p j =jj <i n j =i j −1j >i当i =n 时，(−1)f (p )ΠA i ,p i=λn −a 1λn −1当i >1时，f (p )=n −iΠA i ,p i=(−1)n −i +1λi ⋅a n −i +1(−1)f (p )ΠA i ,p i=−λi a n −i +1综上,转移矩阵A 的特征多项式有简单的表达p (λ)=|λI n −A |=λn −a 1λn −1−a 2λn −2−⋯−a n假设有f 0这⼀项(不需要知道是多少)，那么认为初始向量列为S =(f −(n −1),f −(n −2),⋯,f 0)这个问题，我们要求的是S ⋅A k 的第n 项，不需要知道整个矩阵类似求出A k 的过程，求出F_k(x)\mod p(\lambda)我们要求解(S\cdot A^k)_n=\sum_1^{n}[x^i]{F(x)}(S\cdot A^i)_n⽽(S\cdot A^i)_n=f_i 已知，求出F(x)后直接带⼊即可需要⽤到多项式取模，求解这个表达式是O(n\log n\log k)的，求完直接带⼊即可{Loading [MathJax]/extensions/TeX/mathchoice.js。
2阶常系数线性齐次递推关系Linear Homogeneous Relation of Degree 2•如果an的递推关系满足a n+C1a n-1+C2a n-2+…+C k a n-k=0，且初值为a0=d0, a1=d1, …a k-1=d k-1，则称这个等式为k阶常系数线性齐次递推关系（linear homogeneous relation of degree k）•多项式x k +C1x k-1+r2x k-2+…+r k=0称为它的特征多项式或特征方程（characteristic equation），其根称为特征根（characteristic root）。
是否常系数线性齐次递推关系？a n = a n-1⋅n Xa n = a n-1+2n-1 Xa n = 3a n-1+1 Xa n = 4a n-1- 4a n-2 √f n = f n-1 + f n-2 √T n = 2T n-1+1 Xx n+1 = x n(1-x n) X•例1•f= 1, f2 = 1, f n = f n-1+f n-21•特征方程是x2-x-1=0•例2•a= 1, a2 = 3, a n = 4a n-1-4a n-21•特征方程是x2-4+4=0•下面我们给出2阶常系数线性齐次递推关系的解法：•假设α, β是a n=c1a n-1+c2a n-2 的特征方程x2-c1x-c2=0 的两个根•a=(α+β)a n-1-(α⋅β)a n-2n•容易验证有a-α⋅a n-1=β⋅(a n-1-α⋅a n-2)n•递推可以得到：a n-αa n-1=β(a n-1-αa n-2)=β2(a n-2-αa n-3)=... =βn-1(a1-αa0)•由此倒推得到：a n-αn a0=(βn-1+αβn-2+α2βn-3+...+αn-1)(a1-αa0)•下面我们给出2阶常系数线性齐次递推关系的解法：•假设α, β是a n=c1a n-1+c2a n-2 的特征方程x2-c1x-c2=0 的两个根•若α≠β，则a n=uαn+vβn，其中u, v由初值决定•若α =β，则a n=a0⋅αn+(a1-αa0)⋅n⋅αn-1•例1•f 1 = 1, f 2 = 1, f n = f n -1+f n -2•特征方程是 x 2-x -1=0•两个特征根分别是 和 •于是 f n = us 1n + vs 2n•由 f 1 = 1和 f 2 = 1 得到1=us 1+vs 2 及1 = us 12 + vs 22 •解出1152+=s 2152-=s 15=u 15=-v•于是，斐波那契数列的第 n 项是1151152255n ⎛⎫⎛⎫+-=- ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭n n f•例2•a= 1, a2 = 3, a n = 4a n-1-4a n-21•特征方程是x2-4+4=0•特征根是α= 2，重数为二•可以求出a= 1/4•于是a n = a0⋅αn+(a1-αa0)⋅n⋅αn-1= (1/4)⋅2n+(1-2⋅1/4)⋅n⋅2n-1= (1+n)⋅2n-2E nd。
递推关系递归公式是用它自身来定义的一个公式，我们习惯称之为递推关系或递推式。
如正奇数序列可以用递推式描述为：f(n)=f(n-1)+2, n>1 且f(1)=1当n为很大的值时，直接用递推来计算f(n)会很麻烦，所以希望能够用一种封闭的式子来描述这个序列，从它入手可以直接计算f(n)。
如果找到这样一种封闭的式子，则称递推式已经解出。
下面的内容给出了求解基本的递推式的一些方法。
递推关系如果具有如下这种形式，则称为常系数线性齐次递推式：f(n)=a1f(n-1)+a2f(n-2)+…+a k f(n-k)这里f(n)称为k次的。
当一个附加项包括常数或者n的函数出现在递推中，那么它就称为非齐次的。
一、线性齐次递推式的求解令f(n)=a1f(n-1)+a2f(n-2)+…+a k f(n-k)的一般解含有f(n)=x n形式的特解的和。
用x n来代替上式中的f(n)，得到：x n =a1x n-1+a2 x n-2 +…+a k x n-k两边同时除以x n-k得到：x k =a1x k-1+a2 x k-2 +…+a k或者写成x k -a1x k-1-a2 x k-2 -…-a k =0以上两等式都称为原递推关系的特征方程。
下面我们只限于一阶和二阶的线性递推关系。
一阶齐次递推方程的解可以直接得到，令f(n)=af(n-1)，假定递推序列从f(0)开始，由于f(n)=af(n-1)=a2f(n-2)=…=a n f(0)所以f(n)=a n f(0)是递推的解。
如果递推的次数是2，那么特征方程变为x2-a1x-a2=0,令这个二次方程的根是r1和r2，递推的解是：f(n)=c1r1n+c2r2n(r1≠r2)f(n)=c1r n+c2nr n(r1=r2)代入序列初始的值f(n0)和f(n0+1)解方程得到c1和c2的值。
例1序列1,4,64,256,…可以用递推关系表示为f(n)=3f(n-1)+4f(n-2)，且f(0)=1,f(1)=4，求此递推式的解。
常系数齐次线性递推数列的周期性探究摘要: 本文借助初等数列通项推导法与高等代数的矩阵法分别对复数系二阶常系数递推数列周期的充要条件进行详尽、系统的探究.关键词: 递推数列 周期性 通项推导法 矩阵法二阶常系数齐次线性递推数列是高中数学中常见的数列,笔者发现有不少的题目考察了其数列的周期性,那么此类数列是周期数列的充分必要条件是什么呢？这引起了笔者的兴趣,故对其进行研究,并查阅相关文献,有了一些感悟,现与大家分享.在文[1]中,对于二阶常系数线性递推数列{}n a :()3,*21≥∈+=--n N n qa pa a n n n ,以及初值12,a a ()22120+a a≠其中,特征方程为2x px q =+,24p q ∆=+,按照0,0,0<=>ΔΔΔ共3种情况,分别进行}{n a 为周期数列的充分必要条件探讨,并给出了周期的判断方法.深究发现,文[1]存在两个漏洞,一是忽视初值对数列周期的影响,故0∆>所得结论欠缺严谨性.二是未限定基本条件12,,,p q a a R ∈,对0∆<的情况,却允许其特征根为复数.基于文[1]的探讨,本文将借助初等数列通项推导法与高等代数的矩阵法两种方式分别作详尽、系统的探讨:当12,,,p q a a C ∈时,二阶常系数递推数列周期的充要条件. 1 二阶常系数齐次线性递推数列的周期性1.1 二阶常系数齐次线性递推数列的概念满足递推关系21n n n a pa qa ++=+①,以及给定初值12,a a (其中12,,,p q a a C ∈是常数)的数列{}n a 称为二阶常系数齐次线性递推数列(下文简称“二阶递推数列”).把方程2x px q =+称为递推数列的特征方程,特征方程2x px q =+的两个复数根12,λλ称为递推数列的特征根.注:由于笔者的习惯,在下文中递推关系用21n n n a pa qa ++=+,而非文[1]的()3,*21≥∈+=--n N n qa pa a n n n .1.2二阶递推关系的周期分类由于二阶递推数列的周期性是受初值影响的,详例如下: (1)递推关系:2n n a a +=对任意的12,a a C ∈,数列{}n a 周期数列;(2)递推关系:213122n n n a a a ++=--当122,2a a ==-时,数列{}n a 的通项公式为2n a =,是周期数列;当121,3a a ==时,213+2n n a -⎛⎫=- ⎪⎝⎭并不是周期数列;(3)递推关系:21n n na a a ++=+对任意不全为零的初值,数列{}n a 均无界,故{}n a 不是周期数列（证明略）; 因此,二阶递推数列的周期定义—弱递推关系()21,n n n a pa qa p q C ++=+∈满足: (1)对任意的初值12,a a C ∈,数列{}n a 都是周期数列,则称递推关系①为周期递推关系; (2)存在的不全为零的初值12,a a C ∈,使得数列{}n a 是周期数列,也存在初''12,a a C ∈,使得数列{}n a 不是周期数列,则称递推关系①为弱周期递推数列;(3)对任意不全为零的初值12,a a C ∈,数列{}n a 都不是是周期数列,则称递推关①为非周期递推数列;2 复数系二阶常系数齐次线性递推数列周期的充要条件探究 2.1 初等数列通项推导法2.1.1理论基础引理1:复数列{}n a 是等比数列,其公比为q .则“数列{}n a 是以()*T T N ∈为周期数列”的充要条件是“1Tq =”.证明:因为{}n a 是等比数列,其公比为q ,且对任意的*,0n n N a ∈≠. 故1T Tn T n n n a a a q a q +=⇔⋅=⇔=.引理2:复数列{}n a 是以T 为周期数列,则对任意的数列C λ∈,{}1n n a a λ+-也是以T 为周期数列.证明:1111,n T n n T n n T n T n n a a a a a a a a λλ++++++++==⇒-=- ,故{}1n n a a λ+-也是以T 为周期数列.(逆命题不一定成立)2.1.2探究过程由韦达定理1212,p q λλλλ+==,于是将递推关系改写为:()()()211211212112221112n n n n n n n n n n n a a a a a a a a a a a λλλλλλλλλλ++++++++-=-⎧⎪∴=+-⇔⎨-=-⎪⎩②③然而并不能就此断定{}11n n a a λ+-与{}12n n a a λ+-都是等比数列,故分情况讨论: 情形(1) 211=0a a λ-,211=0a a λ-中至少一个为0时,不妨设211=0a a λ-,此时, 恒有011=-+n n a a λ,11n n a a λ+=.至此依然不能断定{}n a 为等比数列.所以继续分类讨论:(i)若10a =,此时0=n a ,即{}n a 为零数列;(ii) 若110,0a λ≠=,此时{}n a 除了首项之外,其余各项均为0,此数列不是周期数列;(iii)若 110,0a λ≠≠,此时{}n a 是等比数列,111n n a a λ-=,由引理1可知,数列{}n a是以T 为周期数列当且仅当1=1Tλ.情形(2) 211221,a a a a λλ--均不为0,12,λλ中至少有一为0时,不妨设10λ=. 此时,数列{}12n n a a λ+-除了首项之外,其余各项均为0.此时,{}12n n a a λ+-不是周期数列,从而,由引理2可知,数列{}n a 也不可能是周期数列.情形(3) 21122112,,,a a a a λλλλ--均不为0,此时,数列{}11n n a a λ+-以及{}12n n a a λ+-是等比数列.由②式得,()1112112n n n a a a a λλλ-+-=-⋅ ②由②式得,()1122211n n n a a a a λλλ-+-=-⋅ ⑤由引理2,若复数列{}n a 是以T 为周期数列,则{}11n n a a λ+-以及{}12n n a a λ+-也是以T为周期数列.由引理1,须满足121,1T Tλλ==,然而这只是情形(3)下数列为周期数列的必要条件,对于它的充分性,仍需继续依据12,λλ是否相等进行分类讨论:(i) 当12λλ≠时,则由-③④,得()()()112121122211n n n a a a a a λλλλλλ---=-⋅--⋅11211221212112n n n a a a a a λλλλλλλλ--⎛⎫⎛⎫--⇒=⋅+⋅ ⎪ ⎪--⎝⎭⎝⎭.此时,检验可知n T n a a +=成立,即数列{}n a 为周期数列;(ii)当12λλλ==时,则②和⑤式可以并为:()1121n n n a a a a λλλ-+-=-⋅12112n nn n a a a a λλλ+--⇒-=-⇒数列2n n aλ-⎧⎫⎨⎬⎩⎭是等差数列,公差为21a a λ-. 于是()()()21212112212n nn n a a a a n a a a n a a λλλλλλ--=+--⇒=-+-⎡⎤⎣⎦记2112,2k a a b a a λλ=-=-,则()2n n a kn b λ-=+.因为211,20T a a λλ=-≠,此时()()22n T n n T a k n T b k n T b λλ+--+=++=++⎡⎤⎡⎤⎣⎦⎣⎦220n n n a a k T λ-+-=⋅⋅≠,故此时{}n a 不是T 周期数列.2.1.3结论综述对于二阶递推复数列{}n a :()22,,n n n a pa qa p q C ++=+∈,其初值为12,a a C ∈,其特征根为12λλ，,有以下4个结论成立:结论1 {}n a 是()*T T N ∈周期数列当且仅当满足如下3个条件之一: (1) 初值120a a ==;(2) 存在某个特征根λ,使得=1T λ,且210a a λ-=(3) 有两个不相等的特征根12λλ，,且12==1T Tλλ.结论2 (i)“递推关系①为周期递推关系”当且仅当“其特征根12,λλ不相等,且存在*T N ∈使得12==1T Tλλ”; (ii)“递推关系①为弱周期递推关系”当且仅当“有且只有一个特征根(两个相等的特征根视为一个)λ,使得存在*T N ∈使得=1Tλ”;(iii)“递推关系①为非周期递推关系”当且仅当“任意的*T N ∈,121,1T Tλλ≠≠对于以上结论进行实例应用,如下:例1 递推关系:212n n n a a a ++=-,其特征根为12±,而1=≠,所以,对任意的*T N ∈,于是,121,1T T λλ≠≠.故212n n n a a a ++=-为非周期递推关系.例 2 ()211i n n n a a a ++=-+,121,i a a ==.其对应的特征根为: ()12i,1i λλ==-,满足:421i 0,i 1a a -⋅==可知数列以4为周期的周期数列.例 3 其对应的特征根为:()12i,1i λλ==-,21λ=≠,故对任意的*2,1T T N λ∈≠,而411λ=,但21i 0a a -≠,故此数列不是周期数列.例 4 21n n n a a a ++=-,其特征根为:12λλ==,6612=1λλ=.故以6T =为周期的递推关系.进一步,考虑1Tλ=的等价条件,可知λ在单位内,设其辐角为θ,则i e θλ=,21,=T k k N T πλθ=⇔∃∈使得.记2kt Q T=∈,则得到二阶递推数列{}n a 为周期数列的一个必要条件:结论 3 二阶递推数列{}n a 为周期数列,则必存在特征根λ,使得其辐角为,t t Q θπ=∈.由于实数域包含于复数域内,结合文[1]以及本文中的结论1,结论2以及结论3,可得到实数范围内相应的结论（4）,其中,情形(1)(2)(3)对应与0∆≥,情形(4)对应于0∆<:结论4 递推关系n n n qa pa a +=++12,以及给定初值21,a a (其中221212,,,,0p q a a R a a ∈+≠),数列{}n a 是周期数列当且仅当满足如下四个情形之一:(1) 1p q +=,21a a =.此时数列{}n a 是1为最小正周期的数列; (2) 1,q p -=12a a =-.此时,数列{}n a 是2为最小正周期的数列;(3) 0,1p q ==,对任意的12,a a ,数列{}n a 是1或2为最小正周期的数列;(4) 2cos ,1p t q π==-(其中*,,,st s r N s r s r r=∈<与互质,),对任意的12,a a ,数列{}n a 均为周期数列.若s 为奇数,则数列{}n a 是2r 为周期的数列;若r 为偶数,则数列{}n a 是r为周期的数列.2.2 高等代数的矩阵法由于初等的方法较复杂,受文[2]运用矩阵法求常系数递推数列的启发,所以尝试引入线性代数的相关知识,将数列的问题转化为矩阵的问题,实现高效推导.2.2.1 常系数递推数列与递推矩阵一般地,称满足下述递推公式的数列{}n a 为常系数齐次递推数列:递推关系1122n k n k n k k k a c a c a c a ++-+-=+++,以及初值:12,,,k a a a其中,1212,,,,,,,k k c c c a a a C ∈为常数.设1122111n k k n k n n k kk a c c c a X A a +-+-⨯⨯⎛⎫⎛⎫⎪ ⎪⎪ ⎪== ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭，,则常系数线性递推数列可以转化为 向量数列{}n X :满足递推关系1n n X AX +=,初始条件1X .称A 为递推矩阵,通过反复迭代,可得11n n X A X -=.2.2 .2 递推向量列的周期性研究若初始条件1X 为零向量,则{}n X 为恒为0向量,结论平凡,不予讨论.再者,由于{}n a 是周期数列等价于{}n X 是周期向量数列.故只需研究非零初值的递推向量列{}n X 的周期即可.(一) 递推数列为周期数列的充要条件设{}n X 是以T 为周期的向量数列.则1111TT X X A X X +=⇔=(注:这里表示A 的T 次幂,而不是矩阵的转置,下文同),即1是矩阵T A 的特征值,1X 是矩阵TA 的属于1的特征向量.另一方面,若11TA X X =成立,则()111111n T n T n n T n X A X A A X A X X +---+====,即有:{}n X 是以T 为周期的向量数列⇔1是矩阵TA 的特征值,1X 是矩阵TA 的属于1的特征向量.进一步,设矩阵A 的特征值为12,,,k λλλ,则矩阵T A 的全部特征值为12,,,T T T k λλλ,于是,得到:已知递推数列{}n a :1122n k n k n k k k a c a c a c a ++-+-=+++,在非零初值:12,,,k a a a .结论 5 {}n a 是以T 为周期的数列,当且仅当,存在存在特征值λ,使得1T λ=,且121n a a X a ⎛⎫ ⎪ ⎪= ⎪⎪⎝⎭是矩阵T A 对应于1的特征向量,其中A =1211k c c c ⎛⎫⎪⎪ ⎪ ⎪⎝⎭. （二）周期递推关系与弱周期递推关系 递推关系:1122n k n k n k k k a c a c a c a ++-+-=+++为周期递推关系⇔存在*T N ∈,使得任意的1k X C ⨯∈,都有T A X X ⋅=⇔()0T A E X -⋅=恒成立⇔TA E =.另一方面,对1211k c c c E A λλλλ---⎛⎫⎪-⎪-= ⎪ ⎪-⎝⎭实行初等矩阵变换, 得矩阵A 的Smith 标准型:1111k k k c c λλ-⎛⎫ ⎪⎪⎪ ⎪---⎝⎭由此可见,矩阵A 的最小多项式等于其特征多项式,即:()()11det k k k m E A c c λλλλ-=-=---于是T A E =()1T fλλ⇔=-是矩阵A 的化零多项式()()|m f λλ⇔()()111k k T k c c λλλ-⇔----⇔12,,,k λλλ是互不相等相等的T 次单位根. 结论6 递推关系:1122n k n k n k k k a c a c a c a ++-+-=+++为周期递推关系与下列条件等价:(1)存在正整数T ,使得()()111k k T k c c λλλ-----(2)存在正整数T ,使得12,,,k λλλ是互不相等相等的T 次单位根.(3)12,,,k λλλ均在单位圆上,且其辐角主值与π的比值均为有理数.同样,以上结论举例应用如下: 例5 3212n n n n a a a a +++=++,对应的递推矩阵112100010A ⎛⎫⎪= ⎪ ⎪⎝⎭,其特征值多项式为()()322212λλλλλλ---=++-,其特征值1232λλλ===,满足:33121λλ==.于是考虑3544232112A ⎛⎫⎪= ⎪ ⎪⎝⎭,3A 特征值为1,1,8.其中1对应的特征向量为12100,111X X ⎛⎫⎛⎫⎪ ⎪== ⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭.所以,当123,0,a a a a a ===-或1230,,a a a a a ===-时,数列{}n a 为周期数列;对于其他的初值,数列{}n a 均不为周期数列.所以,递推关系3212n n n n a a a a +++=++为弱周期递推关系.例6321n n n n a a a a +++=-+,对应的特征方程为:()()322111λλλλλ-+-=+-,其特征根为123,,1i i λλλ==-=,均是4次单位根.故递推关系321n n n n a a a a +++=-+为以4为周期递推关系.参考文献[1] 朱 勇 对一类递推数列周期性的再探究[J].中学数学教学,2014(4), 24-25[2] 李信明 常系数线性递推数列通项公式的矩阵求法[J] 高等数学研究,1999(4),23-25.。
