浅谈特征方程及相关证明

前言

特征方程应该是大学里的内容,但最近做题的时候遇到了,就想把我的一点心得和大家分享一下。
但由于鄙人水平有限,故以下只讨论二阶常系数线性齐次递推式。

问题

已知f(n)=c1∗f(n−1)+c2∗f(n−2)f(n)=c1*f(n-1)+c2*f(n-2)f(n)=c1f(n1)+c2f(n2)c1,c2c1,c2c1,c2是常数),已知f(0)f(0)f(0)f(1)f(1)f(1),求f(n)f(n)f(n)的通项公式。

结论

先求出上面递推式的特征方程:x2−c1∗x−c2=0x^2-c1*x-c2=0x2c1xc2=0。设两根分别为x1,x2x1,x2x1,x2
x1≠x2x1\neq x2x1̸=x2,则f(n)=A∗x1n+B∗x2nf(n)=A*x1^n+B*x2^nf(n)=Ax1n+Bx2n
x1=x2x1=x2x1=x2,则f(n)=(A+B∗n)∗x1nf(n)=(A+B*n)*x1^nf(n)=(A+Bn)x1n。(A和B可通过f(0)f(0)f(0)f(1)f(1)f(1)求出)

例题

已知f(n)=4f(n−1)−3f(n)f(n)=4f(n-1)-3f(n)f(n)=4f(n1)3f(n)f(0)=3f(0)=3f(0)=3f(1)=5f(1)=5f(1)=5,求f(n)f(n)f(n)的通项公式。
解:
特征方程为:x2−4x+3=0x^2-4x+3=0x24x+3=0
x1=1,x2=3x1=1,x2=3x1=1,x2=3
∵x1≠x2\because x1\neq x2x1̸=x2
∴f(n)=A+B∗3n\therefore f(n)=A+B*3^nf(n)=A+B3n
当n=0时,3=A+B3=A+B3=A+B;当n=1时,5=A+3B5=A+3B5=A+3B
解得A=2,B=1A=2,B=1A=2,B=1
∴f(n)=3n+2\therefore f(n)=3^n+2f(n)=3n+2

证明

我们可以把递推式转化成一个类似等比数列的东西。
f(n)−r∗f(n−1)=s(f(n−1)−r∗f(n−2))f(n)-r*f(n-1)=s(f(n-1)-r*f(n-2))f(n)rf(n1)=s(f(n1)rf(n2)),则f(n)=(s+r)∗f(n−1)−r∗s∗f(n−2)f(n)=(s+r)*f(n-1)-r*s*f(n-2)f(n)=(s+r)f(n1)rsf(n2)
可得s+r=c1,s∗r=−c2s+r=c1,s*r=-c2s+r=c1,sr=c2
根据韦达定理,s和r是x2−c1∗x−c2=0x^2-c1*x-c2=0x2c1xc2=0的两根
x2−c1∗x−c2=0x^2-c1*x-c2=0x2c1xc2=0称为该递推式的特征方程,两根分别为x1,x2x1,x2x1,x2
不妨设x1=r,x2=sx1=r,x2=sx1=r,x2=s,则f(n)−x1∗f(n−1)f(n−1)−x1∗f(n−2)=x2\frac{f(n)-x1*f(n-1)}{f(n-1)-x1*f(n-2)}=x2f(n1)x1f(n2)f(n)x1f(n1)=x2
f(1)−x1∗f(0)=af(1)-x1*f(0)=af(1)x1f(0)=a
f(2)−x1∗f(1)=a∗x2f(2)-x1*f(1)=a*x2f(2)x1f(1)=ax2
……
f(n−2)−x1∗f(n−3)=a∗x2n−3f(n-2)-x1*f(n-3)=a*x2^{n-3}f(n2)x1f(n3)=ax2n3
f(n−1)−x1∗f(n−2)=a∗x2n−2f(n-1)-x1*f(n-2)=a*x2^{n-2}f(n1)x1f(n2)=ax2n2
f(n)−x1∗f(n−1)=a∗x2n−1f(n)-x1*f(n-1)=a*x2^{n-1}f(n)x1f(n1)=ax2n1
③+x1∗x1*x1②得:f(n)−x12∗f(n−2)=a∗x2n−1+a∗x1∗x2n−2f(n)-x1^2*f(n-2)=a*x2^{n-1}+a*x1*x2^{n-2}f(n)x12f(n2)=ax2n1+ax1x2n2
④+x12∗x1^2*x12①得:f(n)−x13∗f(n−3)=a∗x2n−1+a∗x1∗x2n−2+a∗x12∗x2n−3f(n)-x1^3*f(n-3)=a*x2^{n-1}+a*x1*x2^{n-2}+a*x1^2*x2^{n-3}f(n)x13f(n3)=ax2n1+ax1x2n2+ax12x2n3
发现规律了吗?

f(n)−x1n∗f(0)=a∗(x2n−1+x1∗x2n−2+x12∗x2n−3+……+x1n−2∗x2+x1n−1)f(n)-x1^n*f(0)=a*(x2^{n-1}+x1*x2^{n-2}+x1^2*x2^{n-3}+……+x1^{n-2}*x2+x1^{n-1})f(n)x1nf(0)=a(x2n1+x1x2n2+x12x2n3++x1n2x2+x1n1)
=a∗∑i=1n(x1n−1∗(x2x1)i−1)a*\sum_{i=1}^n({x1^{n-1}*(\frac{x2}{x1})^{i-1}})ai=1n(x1n1(x1x2)i1)
∑i=1n(x1n−1∗(x2x1)i−1)\sum_{i=1}^n({x1^{n-1}*(\frac{x2}{x1})^{i-1}})i=1n(x1n1(x1x2)i1)可以看成是以x1n−1x1^{n-1}x1n1为首项,x2x1\frac{x2}{x1}x1x2为公比的等比数列的前n项的和

在运用等比数列求和公式之前一定要讨论公比是否为1,接下来开始讨论:
1.当x1≠x2x1\neq x2x1̸=x2时:∑i=1n(x1n−1∗(x2x1)i−1)\sum_{i=1}^n({x1^{n-1}*(\frac{x2}{x1})^{i-1}})i=1n(x1n1(x1x2)i1)=x1n−1∗(x2x1)n−1x2x1−1=x2n−x1nx2−x1x1^{n-1}*\frac{(\frac{x2}{x1})^n-1}{\frac{x2}{x1}-1}=\frac{x2^n-x1^n}{x2-x1}x1n1x1x21(x1x2)n1=x2x1x2nx1n
f(n)=x1n∗f(0)+a∗∑i=1n(x1n−1∗(x2x1)i−1)=(f(0)−ax2−x1)∗x1n+ax2−x1∗x2nf(n)=x1^n*f(0)+a*\sum_{i=1}^n({x1^{n-1}*(\frac{x2}{x1})^{i-1}})=(f(0)-\frac{a}{x2-x1})*x1^n+\frac{a}{x2-x1}*x2^nf(n)=x1nf(0)+ai=1n(x1n1(x1x2)i1)=(f(0)x2x1a)x1n+x2x1ax2n
A=f(0)−ax2−x1,B=ax2−x1A=f(0)-\frac{a}{x2-x1},B=\frac{a}{x2-x1}A=f(0)x2x1a,B=x2x1a,则有f(n)=A∗x1n+B∗x2nf(n)=A*x1^n+B*x2^nf(n)=Ax1n+Bx2n
该情况证明完毕
2.当x1=x2x1=x2x1=x2时,∑i=1n(x1n−1∗(x2x1)i−1)=a∗n∗x1n−1\sum_{i=1}^n({x1^{n-1}*(\frac{x2}{x1})^{i-1}})=a*n*x1^{n-1}i=1n(x1n1(x1x2)i1)=anx1n1
f(n)=a∗n∗x1n−1+x1n∗f(0)=(a∗nx1+f(0))∗x1nf(n)=a*n*x1^{n-1}+x1^n*f(0)=(\frac{a*n}{x1}+f(0))*x1^nf(n)=anx1n1+x1nf(0)=(x1an+f(0))x1n
A=f(0),B=ax1A=f(0),B=\frac{a}{x1}A=f(0),B=x1a,则有f(n)=(A+B∗n)∗x1nf(n)=(A+B*n)*x1^nf(n)=(A+Bn)x1n
至此,命题证明完毕

应用

求斐波那契数列通项公式
f(n)=f(n−1)+f(n−2),f(0)=0,f(1)=1f(n)=f(n-1)+f(n-2),f(0)=0,f(1)=1f(n)=f(n1)+f(n2),f(0)=0,f(1)=1
根据特征方程:x2−x−1=0x^2-x-1=0x2x1=0
x1=1+52,x2=1−52x1=\frac{1+\sqrt{5}}{2},x2=\frac{1-\sqrt{5}}{2}x1=21+5,x2=215
那么,f(n)=A∗x1n+B∗x2nf(n)=A*x1^n+B*x2^nf(n)=Ax1n+Bx2n
代入f(0)=0,f(1)=1f(0)=0,f(1)=1f(0)=0,f(1)=1,得:A=55,B=−55A=\frac{\sqrt{5}}{5},B=-\frac{\sqrt{5}}{5}A=55,B=55
整理得:f(n)=(1+52)n−(1−52)n5f(n)=\frac{(\frac{1+\sqrt{5}}{2})^n-(1-\frac{\sqrt{5}}{2})^n}{\sqrt{5}}f(n)=5(21+5)n(125)n

总结

对于这个问题来说,其实结论比证明更加重要。但我想要在网上找证明,却很难找到,无奈,只好自己证明。不知是这个问题太过简单,各位高手不愿证明,还是别的证明比较复杂,难度较高。总之,这些证明只是我的一点小想法,希望能抛砖引玉。

评论 13
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值