• 二阶常系数线性齐次递推式的特征方程


    参照liuzibujian的博客。

    问题

    已知(f(n)=c_1∗f(n−1)+c_2∗f(n−2))(c_1,c_2) 是常数),已知(f(0))(f(1)),求(f(n))的通项公式。

    结论

    先求出上面递推式的特征方程:(x^2-c_1x-c_2=0)(式子有点像解(n)次方程)。设两根分别为(x_1,x_2)
    (x_1≠x_2),则(f(n)=A*x_1^n+B*x_2^n)
    (x1=x2),则(f(n)=(A+B∗n)∗x_1^n) 。((A)(B)可通过(f(0))(f(1))求出)

    例题

    已知(f(n)=4f(n-1)-3f(n),f(0)=3,f(1)=5),求(f(n))的通项公式。
    解:
    特征方程为:(x^2-4x+3=0)
    (x_1=1,x_2=3)
    (ecause x_1 e x_2)
    ( herefore f(n)=A+B*3^n)
    (n=0)时,(3=A+B);当(n=1)时,(5=A+3B)
    解得(A=2,B=1)
    ( herefore f(n)=3^n+2)

    证明

    我们可以把递推式转化成一个类似等比数列的东西。

    [f(n)-r*f(n-1)=s(f(n-1)-r*f(n-2)) ]

    [f(n)=(s+r)*f(n-1)-r*s*f(n-2) ]

    可得(s+r=c_1,s*r=-c_2)
    根据韦达定理,(s)(r)(x^2-c_1x-c_2=0)的两根
    (x^2-c_1*x-c_2=0)称为该递推式的特征方程,两根分别为(x_1,x_2)
    不妨设(x_1=r,x_2=s),则(frac{f(n)-x_1*f(n-1)}{f(n-1)-x_1*f(n-2)}=x_2)

    [f(1)-x_1*f(0)=a ]

    [f(2)-x_1*f(1)=a*x_2 \ dots \ f(n-2)-x_1*f(n-3)=a*x_2^{n-3}① \ f(n-1)-x_1*f(n-2)=a*x_2^{n-2}②\ f(n)-x_1*f(n-1)=a*x_2^{n-1}③\ ]

    (③+x_1*②)得:

    [f(n)-x_1^2*f(n-2)=a*x_2^{n-1}+a*x_1*x_2^{n-2}④ ]

    (④+x_1^2*①)得:

    [f(n)-x_1^3*f(n-3)=a*x_2^{n-1}+a*x_1*x_2^{n-2}+a*x_1^2*x_2^{n-3} ]

    发现规律了吗?

    [f(n)-x_1^n*f(0)=a*(x_2^{n-1}+x_1*x_2^{n-2}+x_1^2*x_2^{n*3}+dots+x_1^{n-2}*x_2+x_1^{n-1}) \ =a*sum_{i=1}^n(x_1^{n-1}*(frac{x_2}{x_1})^{i-1}) ]

    (sum_{i=1}^n(x_1^{n-1}*(frac{x_2}{x_1})^{i-1}))可以看成是以(x_1^{n-1})为首项,(frac{x_2}{x_1})为公比的等比数列的前(n)项的和。

    在运用等比数列求和公式之前一定要讨论公比是否为(1),接下来开始讨论:

    1. (x_1 e x_2)时:

    [sum_{i=1}^n(x_1^{n-1}*(frac{x_2}{x_1})^{i-1}) \ = x_1^{n-1}*frac{(frac{x_2}{x_1})^n-1}{frac{x_2}{x_1}-1} = frac{x_2^n-x_1^n}{x_2-x_1} ]

    所以有

    [f(n)=x_1^n*f(0)+a*sum_{i=1}^n(x_1^{n-1}*(frac{x_2}{x_1})^{i-1}) \ =(f(0)-frac{a}{x_2-x_1})*x_1^n+frac{a}{x_2-x_1}*x_2^n ]

    (A=f(0)-frac{a}{x_2-x_1},B=frac{a}{x_2-x_1}),则有

    [f(n)=A*x_1^n+B*x_2^n ]

    该情况证明完毕。
    2. 当(x_1=x_2)

    [sum_{i=1}^n(x_1^{n-1}*(frac{x_2}{x_1})^{i-1}) \ = a*n*x_1^{n-1} ]

    所以有

    [f(n)=a*n*x_1^{n-1}+x_1^n*f(0)=(frac{a*n}{x_1}+f(0))*x_1^n ]

    (A=f(0),B=frac{a}{x_1}),则有

    [f(n)=(A+B*n)*x_1^n ]

    至此,命题证明完毕。

    应用

    求斐波那契数列的通项公式
    (f(n)=f(n-1)+f(n-2),f(0)=1,f(1)=1)
    根据特征方程:(x^2-x-1=0)

    [x_1=frac{sqrt{5}+1}{2},x_2=frac{sqrt{5}-1}{2} ]

    那么,(f(n)=A*x_1^n+B*x_2^n)
    带入(f(0)=1,f(1)=1),得:(A=frac{sqrt{5}}{5},B=-frac{sqrt{5}}{5})
    整理得:(f(n)=frac{1}{sqrt{5}}left((frac{sqrt{5}+1}{2})^n-(frac{sqrt{5}-1}{2})^n ight))

  • 相关阅读:
    使用sshfs挂载远程服务器目录
    yum install --downloadonly 下载依赖包到本地 但不安装
    CentOS 7安装SSHFS 实现远程主机目录 挂载为本地目录
    CentOS7下FTP的安装与配置
    标签传播算法(Label Propagation Algorithm, LPA)初探
    线性模型(linear model)基本定义及参数求解数学本质、损失函数的选择与评估数学原理、及其基于线性模型衍生的其他机器学习模型相关原理讨论
    贝叶斯A/B测试
    谨慎选择我们的先验
    关于年会抢红包游戏的一个思考
    WMI攻击技术研究
  • 原文地址:https://www.cnblogs.com/autoint/p/10422212.html
Copyright © 2020-2023  润新知