生成斐波那契数列至100项，分析其向黄金比例的收敛，并使用比内公式直接计算特定的第n项。
| n | F(n) | F(n+1) | Fn+1/Fn | 与黄金比例之差 |
|---|
斐波那契数列（Fibonacci Sequence）是意大利数学家莱昂纳多·斐波那契于1202年在其著作《计算之书（Liber Abaci）》中介绍的数列，其中每一项都被定义为前两项之和。数列以F(0) = 0、F(1) = 1开始，之后遵循F(n) = F(n-1) + F(n-2)的递推关系式。前10项为0, 1, 1, 2, 3, 5, 8, 13, 21, 34。例如F(5) = F(4) + F(3) = 3 + 2 = 5，F(6) = F(5) + F(4) = 5 + 3 = 8。这个数列起源于兔子繁殖问题，但实际上出现在数学、自然、艺术、计算机科学等各个领域。斐波那契数呈指数级增长，对于较大的n近似为φⁿ / √5（φ为黄金比例1.618...）。斐波那契数在数论、组合数学、图论等纯数学研究中也扮演重要角色，在分析最大公约数算法（欧几里得辗转相除法）的最坏情况时也会用到。
斐波那契数列中连续两项的比值F(n+1)/F(n)会收敛到黄金比例φ = (1 + √5) / 2 ≈ 1.618033988749...。例如F(6)/F(5) = 8/5 = 1.6，F(10)/F(9) = 55/34 = 1.6176...，F(20)/F(19) = 6765/4181 = 1.6180339...，逐渐接近黄金比例。当n趋向无穷大时，该比值精确收敛于φ。黄金比例自古希腊时代起就被视为最美的比例，被广泛用于帕特农神庙、金字塔以及莱昂纳多·达·芬奇的《蒙娜丽莎》《最后的晚餐》等众多艺术作品中。黄金比例具有φ² = φ + 1和1/φ = φ - 1这样独特的性质，是满足这些关系的唯一正数。从黄金矩形（宽:高 = φ:1）中切下一个正方形后，剩下的仍是黄金矩形，不断重复这一过程便形成黄金螺旋（golden spiral）。自然界中的鹦鹉螺外壳、星系旋臂、飓风结构都遵循这一黄金螺旋。
比内公式（Binet's Formula）是由法国数学家雅克·菲利普·马里·比内发现的闭式公式，无需递归即可直接计算斐波那契数列的第n项：F(n) = (φⁿ - ψⁿ) / √5。其中φ = (1 + √5) / 2 ≈ 1.618033988749（黄金比例），ψ = (1 - √5) / 2 ≈ -0.618033988749。例如计算F(10)：φ¹⁰ ≈ 122.9918，ψ¹⁰ ≈ 0.0081，F(10) = (122.9918 - 0.0081) / 2.236 ≈ 55。递归方法计算F(50)需要数十亿次运算（O(2ⁿ)时间），而比内公式只需一次计算（O(1)时间）即可立即得到结果。不过由于浮点运算精度的限制，当n非常大时（n > 70）可能产生舍入误差，此时使用矩阵幂法（O(log n)）或采用记忆化的动态规划会更精确。比内公式在分析斐波那契数的渐近性质时也很有用。
斐波那契数列惊人地出现在自然界的各个角落。在植物的叶序（phyllotaxis）中最为显著，向日葵的种子以顺时针螺旋34条和逆时针螺旋55条（或55条和89条）排列，全部都是斐波那契数。这种模式能最高效地排布种子，最大化空间利用并均匀地接收阳光。松果的鳞片形成8条和13条螺旋，菠萝则形成8、13、21条螺旋。花瓣的数量也倾向于遵循斐波那契数：百合3瓣、毛茛5瓣、波斯菊8瓣、金盏花13瓣、雏菊21、34瓣。这是因为在进化过程中，以黄金角（golden angle，约137.5°）排布花瓣更有利于生存。树枝的分叉模式、罗马花椰菜的分形结构、DNA双螺旋的尺寸比例（宽21Å，一圈34Å）中也出现了斐波那契数。
斐波那契数列具有各种有趣的数学性质。卡西尼恒等式（Cassini's Identity）：F(n-1) × F(n+1) - F(n)² = (-1)ⁿ。例如当n=5时，3 × 8 - 5² = 24 - 25 = -1 = (-1)⁵。求和公式：F(0) + F(1) + ... + F(n) = F(n+2) - 1。前10项之和0+1+1+2+3+5+8+13+21+34 = 88 = F(12) - 1 = 144 - 1。最大公约数性质：gcd(F(m), F(n)) = F(gcd(m, n))。例如gcd(F(12), F(8)) = gcd(144, 21) = 3 = F(gcd(12, 8)) = F(4)。平方恒等式：F(2n) = F(n) × [2F(n+1) - F(n)]。相邻项的平方和：F(n)² + F(n+1)² = F(2n+1)。例如5² + 8² = 25 + 64 = 89 = F(11)。斐波那契数列的倒数之和收敛，1/F(1) + 1/F(2) + 1/F(3) + ... ≈ 3.359885...。此外，F(n)每隔4项能被3整除（n=0,4,8...），每隔5项能被5整除（n=0,5,10...）。