数学建模差分模型——会生娃的兔子，生生不息

Cherngul

于 2024-10-20 14:48:01 发布

阅读量1.4k

点赞数 41

分类专栏：数学建模文章标签：数学建模 matlab 算法

本文链接：https://blog.csdn.net/m0_74968460/article/details/143092545

版权

数学建模专栏收录该内容

5 篇文章

订阅专栏

“会生娃的兔子，生生不息：揭秘莱斯利种群模型的神奇运作！”

你是否曾好奇过，世界上所有的兔子是怎么保持生生不息的？别着急，今天就让我们一起来了解一个特别有趣的数学模型——莱斯利（Leslie）种群模型。这个模型不仅帮助我们理解兔子种群是如何以不同的年龄段进行增长的，还能带领我们走入更加广泛的生态学与生物学世界。坐稳了，接下来的一段旅程充满了数学、图表、编程以及很多幽默。

1. 莱斯利模型初探：给兔子加个数学帽子

莱斯利种群模型（Leslie population model）是由帕特里克·莱斯利（Patrick Leslie）在1945年提出的一个离散时间的种群模型，主要用于描述一个种群随时间变化的年龄结构。这个模型是基于矩阵的，并且使用线性代数来解决种群动态问题。这个模型非常适用于那些生命周期可以划分为离散阶段的生物，比如兔子、昆虫甚至植物。

2. 莱斯利模型的构成：三要素

莱斯利种群模型主要由三个要素组成：繁殖率、存活率和年龄结构。这三个要素是理解莱斯利模型运作的关键。

繁殖率（Fecundity, $f_i$ ）: 指每个年龄段的个体在下一个时间步生成的新个体数量。
存活率（Survival rate, $s_i$ ）: 指每个年龄段的个体存活到下一个时间步的概率。
年龄结构（Age distribution, $N_i$ ): 指种群在不同年龄段的个体数。

莱斯利矩阵模型通过一个矩阵来描述这些要素，并通过矩阵乘法来预测种群的未来状态。这个矩阵叫做莱斯利矩阵。

莱斯利矩阵看起来就像这样：

$\mathbf{L} = \begin{pmatrix} f_1 & f_2 & \cdots & f_n \\ s_1 & 0 & \cdots & 0 \\ 0 & s_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & s_{n-1} \end{pmatrix}$

其中：

$f_i$ 是第 $i$ 个年龄组的繁殖率，
$s_i$ 是第 $i$ 个年龄组的存活率。

每次的种群状态可以通过将上次的种群状态与莱斯利矩阵相乘来获得：

$\mathbf{L} \cdot N(t)$

其中， $N (t)$ 是时间步 $t$ 时的年龄结构。

3. 实例：兔子种群的增长

现在，假设我们有一个简单的兔子种群模型。我们的兔子分为三个年龄段：0-1 岁、1-2 岁和2-3 岁。以下是我们假设的参数：

繁殖率:
- $f_1 = 0$ （小兔子不生小兔子）
- $f_2 = 1.2$ （成年兔子的繁殖能力）
- $f_3 = 1.8$ （老年兔子的繁殖能力）
存活率:
- $s_1 = 0.5$ （只有一半的小兔子能长大）
- $s_2 = 0.8$ （成年兔子的生存率较高）

我们的初始种群分布如下：

$\begin{pmatrix} 100 \\ 50 \\ 30 \end{pmatrix}$

这意味着我们一开始有100只小兔子、50只成年兔子和30只老年兔子。

4. 用莱斯利矩阵计算下一代兔子

我们可以构建莱斯利矩阵来模拟这个种群的动态变化：

$\mathbf{L} = \begin{pmatrix} 0 & 1.2 & 1.8 \\ 0.5 & 0 & 0 \\ 0 & 0.8 & 0 \end{pmatrix}$

接着我们用这个莱斯利矩阵去计算下一代种群：

$\mathbf{L} \cdot N(0) = \begin{pmatrix} 0 & 1.2 & 1.8 \\ 0.5 & 0 & 0 \\ 0 & 0.8 & 0 \end{pmatrix} \cdot \begin{pmatrix} 100 \\ 50 \\ 30 \end{pmatrix}=\begin{pmatrix} 1.2 \times 50 + 1.8 \times 30 \\ 0.5 \times 100 \\ 0.8 \times 50 \end{pmatrix}=\begin{pmatrix} 126 \\ 50 \\ 40 \end{pmatrix}$

经过一个时间步后，我们的种群变为126只小兔子、50只成年兔子和40只老年兔子。哇哦！兔子种群的扩张可真快。

5. MATLAB实现莱斯利模型

接下来我们用MATLAB来实现这个莱斯利模型，并进行可视化分析。以下是代码：

% 定义莱斯利矩阵
L = [0 1.2 1.8;
     0.5 0 0;
     0 0.8 0];

% 初始种群分布
N = [100; 50; 30];

% 时间步长
T = 20;

% 存储每个时间步的种群数量
populations = zeros(3, T);
populations(:, 1) = N;

% 计算每个时间步的种群
for t = 2:T
    N = L * N;
    populations(:, t) = N;
end

% 绘制结果
figure;
hold on;
plot(1:T, populations(1, :), 'r', 'LineWidth', 2);
plot(1:T, populations(2, :), 'g', 'LineWidth', 2);
plot(1:T, populations(3, :), 'b', 'LineWidth', 2);
legend('0-1岁兔子', '1-2岁兔子', '2-3岁兔子');
xlabel('时间步长');
ylabel('种群数量');
title('莱斯利种群模型：兔子的种群动态');
hold off;