几类特殊递推数列的矩阵算法(前篇)
前言
搁置了许久,终于回来更新以兑现前周立下的flag了。
在学习数列时,我们会遇到这种递推式,其中一种算法是用特征根,而特征根之前做过一期专栏写了个人的奇妙理解方法,感兴趣的可以参考:
而这期专栏,我们可以换种视角,从线性代数的角度来解决。
另外,在做一些数列的题中难免会遇到像,
这种恶心的数列。答案是通过"暴力"反复迭代得出的周期,不禁让人摸不着头脑!
"这谁能想到啊?"[恼]
那么命题人的心机何在呢?今儿就来一探究竟~
另外,此篇专栏会涉及到不少矩阵运算的知识,萌新可以先参考3b1b的视频讲解:
此篇文章不针对应试,而是给数学爱好者们拓展下思路

整式的线性递推
先来讲讲其概念,所谓线性递推,也即(~表系数)形式,这里的
次数均为1
比如为线性递推,而像
这种有次数不为1的就不是线性递推了
再来引入矩阵表示式以便良好地衔接
对于线性方程组,如:
我们想更直观地表示与
的映射关系,可将其化为矩阵形式:
几何意义上,即向量由向量
作线性变换
得到,其中矩阵
就类似于函数(function),只不过函数
的功能是将一个x映射为一个y,而矩阵的功能则是将一个向量映射为另一个向量。由于向量起点均默认于原点,故也可视为将一个坐标(向量终点)映射为另一个坐标

先来看一道例题:已知,求通项公式?
由于笔者能力暂有限,想不到更好地衔接铺垫了,所以提前剧透(
有了上述的铺垫,如果我们能将递推式化为该多好!
为什么好呢?因为这样就有:,
我们便惊奇地发现,这运算不就和等比数列很相近么?
于是得到
因此,要求,只需对初值构成的向量
作n次矩阵A的变换即可
而我们要求,那么就需作n-1次矩阵A的变换

回到题目,先选取已知递推式,再选取另一条式子
写成矩阵形式,即:
于是可形成递推:
则
下面计算
对角化后,得:
代入化简得:
即得:

如果是3阶递推又怎样呢?
比如①
方法是一样的,则考虑构建3阶矩阵:
则还需选取:
②
③
由①②③,写成矩阵形式,即:
递推即得:
后面就也是对角化处理了
不过求特征值时,特征多项式的次数跟矩阵的阶数是对应的,因此上述矩阵特征多项式就是3次方程。而上面的系数是随便选取,因此解析解理论可求,只是可能不平凡。这也是出题时最多只出到二阶线性递推的主要原因。

推广到任意阶线性递推
其中表系数
考虑构建n阶矩阵:
则还需选取:
写成矩阵形式,即:
观察系数矩阵可知,其有如下特点:
第一行的数字即递推式中对应的由到
的系数
从第二行开始,从左往右依次写1,其余写0。也即从第1列第2行起,沿"左上--右下"走向处数字为1,其余为0
举个例子当练习
对应矩阵形式递推,即:

有了以上背景,我们便可以从更高观点来证明这个奇葩式子了
对应矩阵形式递推,即:
递推得:
下面求解
求得矩阵特征根为一对共轭复数:
ps:特征值为复数则对应伸缩+旋转(模长对应伸缩,辐角对应旋转)矩阵运算在复数域成立已有严格证明,这里就先不作拓展了
那么对角化后,则有:
其中
于是有
说明对向量作6次J的变换后回到原位,因此周期为6
到此悟性好的读者脑海中已经掀起了滔天巨浪!
这正是特殊的辐角具有旋转周期性造成的!!!

原来这令人头大的结论背后有如此绝妙的背景!
这便是数学!纯真美妙的数学!不被应试名利铜臭玷污的数学!

总结
这期专栏主要讲解了用矩阵方法求解整式线性递推数列,顺带发掘了一些奇葩结论背后美妙的命题背景。
其中矩阵次方的处理采用了矩阵对角化的方法。而这种方法对于二阶矩阵而言需要满足特征根为2不等实根或一对共轭复根时使用。换而言之,还存在特征根为重根的情况,这种情况矩阵不可对角化。由于个人尚未了解充分,因此先把尚未解决的这种特殊情况先遗留与此。后续掌握后再补充讨论。
整式线性递推写完也花了不少篇幅了,因此一阶分式线性递推就放后续再讲了~
另外,感觉缺少例题的讲解,有空翻翻好像还没丢的一轮书找找例题附上[滑稽]