Matlab演示神奇的斐波那契数列

士兵突击许三多

于 2025-06-12 11:57:34 发布

阅读量1k

点赞数 48

CC 4.0 BY-SA版权

分类专栏： matlab基础文章标签： matlab

本文链接：https://blog.csdn.net/weixin_43470904/article/details/148563431

matlab基础专栏收录该内容

60 篇文章

订阅专栏

Matlab演示神奇的斐波那契数列

斐波那契数列是数学中最迷人且应用广泛的数列之一，它在自然界、艺术和金融等领域都有惊人的出现。本文将介绍斐波那契数列的定义、性质，并展示如何使用MATLAB来实现和可视化这一神奇的数列。

斐波那契数列的定义

斐波那契数列是一个递归定义的整数序列，通常以 $F_n$ 表示第 $n$ 项，其定义如下：

$F_0 = 0, \quad F_1 = 1$

$F_n = F_{n-1} + F_{n-2} \quad \text{对于} \quad n \geq 2$

因此，数列的前几项为：0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144,…

斐波那契数列的性质

斐波那契数列有许多有趣的性质：

黄金分割比：相邻两项的比值 $Fn+1Fn\frac{F_{n+1}}{F_n}$ 趋近于黄金比例 $ϕ=1+52≈1.618033988749895\phi = \frac{1+\sqrt{5}}{2} \approx 1.618033988749895$
比奈公式（闭合形式表达式）：
$F_n = \frac{\phi^n - (-\phi)^{-n}}{\sqrt{5}}$
加法性质： $F_{n+k} = F_k F_{n+1} + F_{k-1} F_n$

MATLAB实现

下面我们展示几种在MATLAB中生成斐波那契数列的方法。

1. 递归方法（最直观但效率低）

function f = fib_recursive(n)
    if n == 0
        f = 0;
    elseif n == 1
        f = 1;
    else
        f = fib_recursive(n-1) + fib_recursive(n-2);
    end
end

2. 迭代方法（高效实现）

function fibs = fib_iterative(n)
    fibs = zeros(1, n+1);
    fibs(1) = 0; % F(0)
    fibs(2) = 1; % F(1)
    for k = 3:n+1
        fibs(k) = fibs(k-1) + fibs(k-2);
    end
end

3. 矩阵方法（基于线性代数）

function f = fib_matrix(n)
    A = [1 1; 1 0];
    f = A^n;
    f = f(2,1);
end

4. 比奈公式法（闭合形式）

function f = fib_binet(n)
    phi = (1 + sqrt(5))/2;
    f = round((phi^n - (-phi)^(-n))/sqrt(5));
end

斐波那契数列的可视化

让我们用MATLAB绘制斐波那契数列和黄金比例收敛的情况：

% 生成前30个斐波那契数
n = 30;
fib = fib_iterative(n);

% 绘制斐波那契数列
figure;
subplot(2,1,1);
stem(0:n, fib, 'filled');
title('斐波那契数列 F_n');
xlabel('n'); ylabel('F_n');
grid on;

% 计算相邻项比值并绘制黄金比例收敛情况
ratios = fib(2:end)./fib(1:end-1);
phi = (1+sqrt(5))/2;

subplot(2,1,2);
plot(1:n, ratios, 'b-', 1:n, phi*ones(1,n), 'r--');
title('相邻项比值 F_{n+1}/F_n 趋近于黄金比例');
xlabel('n'); ylabel('F_{n+1}/F_n');
legend('比值', '黄金比例');
grid on;