Leslie矩阵人口模型的matlab代码

以下是一个简单的 Leslie 矩阵人口模型的 MATLAB 代码示例：

% 参数设置
n = 5; % 人口年龄组数
b = [0 0.75 0.8 0.6 0.3]; % 出生率
s = [0.1 0.4 0.6 0.6 0.4]; % 存活率
L = diag(s(1:n-1), -1); % 存活率矩阵
L(n, :) = [b(1:n-1) 0]; % 出生率行

% 初始人口分布
x = [1000;500;300;200;100]; % 每个年龄组的初始人数

% 模拟时间段
t = 1:20;

% 仿真人口数量变化
y = zeros(n, length(t));
y(:, 1) = x;
for i = 2:length(t)
    y(:, i) = L * y(:, i-1);
end

% 可视化结果
plot(t, y(1,:), 'r', t, y(2,:), 'g', t, y(3,:), 'b', t, y(4,:), 'c', t, y(5,:), 'm');
title('Leslie Matrix population model');
xlabel('Time (years)');
ylabel('Population');
legend('0-4 years', '5-9 years', '10-14 years', '15-19 years', '20-24 years');

这段代码实现了一个假设有5个年龄段的人口模型，其中每个年龄段的出生率和存活率已经预先设置好，使用 Leslie 矩阵来模拟人口数量的变化。运行代码将会生成一个简单的人口数量变化图表。