Leslie矩阵人口模型matlab 代码

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

% 定义 Leslie 矩阵
L = [0.2, 2.0, 0.0, 0.0, 0.0;
     0.6, 0.0, 1.0, 0.0, 0.0;
     0.0, 0.5, 0.0, 1.0, 0.0;
     0.0, 0.0, 0.3, 0.0, 1.0;
     0.0, 0.0, 0.0, 0.2, 0.0];

% 初始人口分布向量
P0 = [100; 50; 30; 20; 10];

% 模拟时间步长
dt = 1;

% 模拟时间范围
T = 20;

% 初始化人口分布矩阵
P = zeros(5, T);

% 设置初始人口分布
P(:, 1) = P0;

% 模拟人口分布随时间的变化
for i = 2:T
    P(:, i) = L * P(:, i-1);
end

% 绘制人口分布随时间的变化图像
plot(1:T, P(1,:), 'r', 1:T, P(2,:), 'g', 1:T, P(3,:), 'b', 1:T, P(4,:), 'm', 1:T, P(5,:), 'c');
legend('Age 0-4', 'Age 5-9', 'Age 10-14', 'Age 15-19', 'Age 20-24');
xlabel('Time (years)');
ylabel('Population size');
title('Leslie Matrix Population Model');

该代码定义了一个 Leslie 矩阵，设置了初始人口分布向量，并使用循环模拟了人口分布随时间的变化。最后，使用 MATLAB 的 plot 函数绘制了人口分布随时间的变化图像。