zoukankan      html  css  js  c++  java
  • 矩阵的特征值 => 矩阵的幂 => 广义斐波拉契数列的通项公式

           众所周知,计算机里动辄涉及到斐波拉契数列,本文主要是运用线性代数的方法求出广义斐波拉契数列的通项公式。

    广义斐波拉契数列的定义(自定义,可能还不够严谨,欢迎指教)如下:

    且a,b满足a^2 + 4b > 0(原因后面有,即保证分母Δ>0)

     

    观察递推公式可知:

        (1)

     

     

     

     

     

     

     

     

     

     

    OK!写这篇博客有两点启示:

    1.亦可用差分方程的方法求通项公式,不过我还没学,我发现把递推关系用矩阵的形式写成再求通项也是可以的,不过要构造出这个关系比较困难,比如(1)式。

    2.在一般的算法中求矩阵的幂几乎都是用二分法求,殊不知若λi下有ki个线性无关的特征向量(ki是λi的重根数),即可用相似对角形的方法求,这在高等数学里可是小儿科。虽然前面的步骤有点繁,但作用是“磨刀不误砍柴工”,直接手算矩阵的幂。

     

     2012/4/7更新:

    今天发现递推公式可简化为:

    这样算起来更方便!

     

    2012/4/15更新:

    想要更刺激的吗?好吧,直接来个n阶的!

    From: Fibonacci数计算中的两个思维盲点及其扩展数列的通用高效解法

     

     

     持续更新中 ing

     

    附:上一张关于菲波拉契数列的图:Fib与指数爆炸

     

     

     

     

     

     

     

  • 相关阅读:
    化了妆的祝福 4
    桥牌感悟 2
    关于送礼
    东京disney sea流水账 1
    东京disney sea流水账 3
    晕倒的候任日本驻华大使
    一饮一啄
    就算神游 之二:行路 4
    关于送礼续
    桥牌感悟 3
  • 原文地址:https://www.cnblogs.com/jjtx/p/2533504.html
Copyright © 2011-2022 走看看