python的leslie人口模型代码

下面是一个简单的 Python 实现 Leslie 人口模型的代码：

import numpy as np

# 初始化模型参数
b = 0.05   # 出生率
l = np.array([0.2, 0.4, 0.3, 0.1])   # 年龄组存活率
m = np.array([0, 0.01, 0.05, 0.1])   # 年龄组生育率
n0 = np.array([1000, 800, 600, 400])   # 初始人口数

# 定义模型函数
def leslie(b, l, m, n):
    n_next = np.zeros_like(n)
    n_next[0] = b * n[3]   # 第一年新生儿数量
    for i in range(1, len(n)):
        n_next[i] = l[i-1] * n[i-1] + m[i-1] * n[i-1]   # 计算每个年龄组的人口数
    return n_next

# 运行模型
n = n0
num_years = 10
populations = [n]
for i in range(num_years):
    n = leslie(b, l, m, n)
    populations.append(n)

# 输出结果
for i, pop in enumerate(populations):
    print("Year {}: {}".format(i, pop))

这个代码中，我们定义了一个 `leslie()` 函数来计算每一年的人口数量。在主函数中，我们使用这个函数来模拟人口数量在 10 年内的变化，并将结果打印出来。