【学习笔记】数列特征方程与特征根

觉得有意思就稍微微写写, 内容大多摘自某本书.


1.不动点求数列通项

对于函数 f(x) , 若存在实数 x0 使得 f(x0)=x0 , 则称 x0 是函数 f(x) 的一阶不动点, 类似的, 若 f(f(x0))=x0 则称 x0 是函数 f(x) 的二阶不动点, 容易发现 f(f(x0))=f(x0)=x0 , 所以一阶不动点一定是二阶不动点.

几何意义上是与直线 y=x 相交的点.

一般的, 数列 {an} 可以由公式 an+1=f(an) 给出, 若满足此条件且存在实数 x0 使得 f(x0)=x0 , 则称 x0 是数列 {an} 的不动点, 不动点可能不存在, 也就是为复数.

可以发现, 数列的不动点性质即为从某一项开始 ak=x0 且再也不变.

有时候, 数列 {an} 中的值可能无法取到 x0 , 但会收敛 x0 , 即 limnan=x0 , 带入递推式可得到 f(x0)=x0 .

2.一阶线性递推数列

一阶线性递推数列通常指以下这种数列, 数列 {an} 满足 an+1=pan+q , 其中 p,q 是给定实数.

这里我们只考虑 p1 的情况, p=1 就是等差数列了.

p1 时我们可以通过待定系数法构造一个公比为 p 的数列, 假设存在实数 x0 , 使得 an+1x0=p(anx0) , 得 x0=q1p .

因此数列 {anq1p} 为等比数列, 通项即为 an=q1p+pn1(a1q1p) .

事实上我们解出的 x0 非常特殊, 是数列 {an} 的不动点.

所以也可以直接在等式两边同时减去不动点解通项.

一阶线性递推通项是 "常数+等比" 的形式.

3.分式递推数列

不动点在分式递推数列的应用才更为优秀, 分式递推数列指一下数列, 数列 {xn} 满足 xn+1=axn+bcxn+d , 其中 a,b,c,d 为给定实数.

这时候求数列的不动点, cx2+(da)xb=0 , 得到一个一元二次, 得到两个根 x1,x2 , 当 x1=x2 时还是减去不动点, 在构造等差数列, 构造的时候要聪明点.

有两个不同的根时, 分别相减再将两式相除, 总得来说, 若只有一个不动点 x0 , 则 {1xx0} 是等差数列, 若有两个不动点, 则 {xnαxnβ} 是等比数列.

4.没有不动点

知乎上看到的, 抄一下.

看一个例子, 设数列 {an} 满足 a1=2 , an+1=12(an1an) , 求通项.

这时候考虑方程 x2+1=0 , 不存在实数不动点, 将不动点扩展到复数域上, 可以得到 ii 是两个不动点, 所以有

an+1i=an212ani=(ani)22anan+1+i=an212an+i=(an+i)22an

两式相除得 an+1+ian+1i=(an+iani)2 , 又 a1+ia1i=2+i2i=35+45i , 迭代得 an+1+ian+1i=(35+45i)2n , 解得 an=(3+4i)2n1+g2n1(3+4i)2n152n1i , 并且根据递推式知道, 虽然通项公式由复数给出, 但每一项都是实数.

仅供此题也有一个漂亮的做法, 注意到 tan2θ=2tanθ1tan2θ , 根据 cotθ=1tanθ , 变形即可得到 cot2θ=12(cotθ1cotθ) .

因此作换元 an=cotθn , 整理得到 θn+1=2θn , 因此 {θn} 是等比数列, 又根据 a1=2θ1=arctan12 , 故 θn=2n1arctan12 , 因此 an=cot(2n1arctan12) .

不动点是复数的时候, 数列可以是周期的, 例如上面这个例子, 以及另一个例子.

设数列 {an} 满足 a1=12 , an+1=1+an1an , 求通项.

其中

an+1=1+an1an=1+1+an11an111+an11an1=1an1=an3

因此数列周期为 4 , 此时数列不动点为复数.

但值得注意的是就算有复数解数列也不一定是周期的.

5.二阶线性递推数列

二阶线性递推通常指以下数列, 数列 {an} 满足 an+1=pan+qan1 , 其中 p,q 是给定常数.

类比一阶线性递推数列, 尝试通过待定系数法构造等比数列, 假设存在实数 a,b , 使得 xn+1axn=b(xnaxn1) .

整理得到 xn+1=(a+b)xnabxn1 , 令 {p=a+bq=ab , 由韦达定理可知实数 a,b 满足方程 x2pxq=0 .

因此 {xn+1axn} 是公比为 b 的等比数列, 累成得 xn+1axn=(x2ax1)bn1 .

同理 xn+1bxn=a(xnbxn1) , xn+1bxn=(x2bxa)an1 , 两式作差解得 xn=x2bx1aban1+x2ax1babn1 .

二阶线性递推通项是 "等比+等比" 的形式.

为了更方便的得到该形式, 我们可以再次用待定系数法, 设 xn=αan1+βbn1 , 列出如下方程:

{x1=α+βx2=αa+βb,{α=x2bx1abβ=x2ax1ba

6.特征根法

终于回归主题了.

设数列 {xn} 的前两项 x1,x2 已知, 且 xn+1=pxn+qxn1 , 则称方程 x2pxq=0 为该数列的特征方程, 该方程有两个根 a,b 则称它们为该数列的特征根.

因此设数列 xn=αan1+βbn1 , 由 {x1=α+βx2=αa+βb , 解得 {α=x2bx1abβ=x2ax1ba .

因此通项为 xn=x2bx1aban1+x2ax1babn1 .

所以特征根其实挺简单且好用的.

当特征根只有一个的时候也可以不用生成函数, 设重根为 r , 只得到一个等比数列 xn+1rxn=(x2rx1)rn1 , 两边同时除以 rn+1 , 变成 xn+1rn+1xnrn=x2rx1r2 , 变成了一个 {xnrn} 为一个公差为 xnr2 的等差数列.

提一嘴, 实际上 xn=Apn+Bqn , xn=αan1+βbn1 , xn=T(Apn+Bqn) 都是一样的.

posted @   Lonely923  阅读(2408)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 使用C#创建一个MCP客户端
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
点击右上角即可分享
微信分享提示