常系数线性递推方程解的一般公式.docVIP

2001年赣南师范学院学报No.6　　　　　　　　　　　　　　　　第6期JournalofGannanTeachersCollegeDec.2001 常系数线性递推方程解的一般公式 吴顺唐 (镇江师范高等专科学校数学系,江苏镇江　212003)Ξ 摘　要:讨论了常系数线性非齐次递推方程初值问题的解,并利用Z变换方法,得到了解的一个新的求解公式。 关键词:递推方程;Z变换 中图分类号:O174.41　　　文献标识码:A-(: y(n1k-)(nk-2)+…+aky(n)+f(n) (0)0(1)b1…,y(k-1)=bk-1,(1) 其中a12,,ak;0,b1,…,bk-1为待定常数,f:N→R为给定函数.在许多实际问题中,特别是在对离散系数作定性分析时,求出方程(1)解的准确表达式,是十分必要的.关于方程(1)的求解,一般都采用特征值法,将(1)化归为一高次代数方程的求根问题.而有些作者则讨论了解的一般公式.[1]和[2],给出了齐次方程初值问题解的一般公式,在[2]中证明了下面的定理: 定理A　当f(n)=0时,常系数线性齐次递推方程初值问题(1)的解为 k-1k-1 yA(n)= j=0a∑(∑i=jk-i+jbi)D(n-k-j)(2) 其中D(r)=　　0,　　　　　　　　　r0;　　i,　　　　　　　　　r=0; kk r+2r+…+kr=ri=112k∑(∑r)ii=1i　r0;ri!r(3) 且当n-k-j0时,D(n-k-j)=0. 公式(2)比较简洁,但[2]是用组合分析的方法,即通过考察赋权图的路径的权和而得到的,而不是用直接的分析方法得到的.另外,在求和公式中还包含了当n-k-j0时的那些多余的项.对于f(n)≠0时的非齐次递推方程,[3]给出了一个一般公式,但其中包含了方程的特征根;[4]和[5]则分别利用不同的方法,得到了一个相同的一般公式,即他们证明了下面的定理: 定理B　常系数线非齐次递推方程初值问题(1)的解为 k-1k-1n-k yB(n)= j=0a∑(∑i=jk-i+jbi)D(n-k-j)+i=0∑D(i)f(n-k-i)　(n≥k)(4) Ξ收稿日期:2001-05-14 　作者简介:吴顺唐,男,江苏常熟人,教授,长期从事逼近论方面的研究工作。 第6期　　　　　　　　吴顺唐　常系数线性递推方程解的一般公式　　　　　　　　　　11其中D(i)如(3)式所示. 我们看到(4)式中也包含了和(2)式相同的多余项.本文则利用与[1],[2]不同的方法,即利用Z变换方法,很方便地给出了初值问题(1)的用给定常数ai,bi与f(n)直接表示且形式比较简单的一个新的一般公式,其结构比[4]和[5]中的结果更为紧凑. 先介绍一些预备知识.设y=y(n)为N到R的映射.如果存在常数M与λ,使得Πn∈ ∞ λ,则称y(n)有限增长的.当y(n)为有限增长时,如果令F(Z)=N,|y(n)|≤Mn n=0n,Z则F(Z)为∞点某邻域内的解析函数,并称F(z)为函数(数列)f(n)的Z变换式,记为 -1Z{f(n)}=F(z),而F(z)到f(n)的相应变换便称为Z逆变换,记为Z.显然,Z变换是线 性变换,且以下结论成立: 引理1　设Z{f(n)}=F(z),则Πk∈N, Z{f(n+k)}=zkF(z)-f(0)kf)-…-()　　引理2　设Z{f1(n)}=F(z,Z{)((z,f(n)为f1(n),f2(n)的 n 卷积:f(n)(f2=0()f2(n-i)则Z{f(n)}=F1(z)F2(z)关于Z变 ,[63　设h(x)=1-a1x-a2x2-…-akxk,则在x=0的邻域内有展开式h(x)∞= n-0∑D(n)xn,其中 kk 　　　　　D(r)=r+2r+…+kr=ri=112k∑(∑r)!ii=1i,　r0;ri! (5)r　　　i,　　　　　　　　　r=0; 　　　0,　　　　　　　　　r 证　当x充分小,使|a1x+a2x2+…+akxk|1时, h(x) ∞ m=0∞=n=0∑(a1x2n+a2x+…+akx)n=∑(s+s+…+s=m12k∑ ∞()ssss+2s+…Ksk)a11a22…akkX12s1!s2!…sk!()ssa11a22…askk)=s1!s2!…sk!将级数中满足条件s1+2s2+…Ksk=n的项集中在一起,当n1时,它们的系数总和为s+2s+…ks=n12k∑(m=0∑　 ∑s+s+…+s=m12k∑ [ s+2s+…ks=n12k(s+s+s)!sssa11a22…akk)]s1!s2!…sk! 而x0的系数显然为1,故(5)式成立.现在,我们考虑非齐次线性递推方程零初值问题 y(n+k)=a1y(n+k-1)+a2y(n+k-2)+