捕食者猎物模型python

以下是一个简单的捕食者猎物模型的Python实现：

import matplotlib.pyplot as plt

# 设置模型参数
alpha = 0.01  # 猎物增长率
beta = 0.001  # 猎物被捕后的死亡率
gamma = 0.1  # 捕食者被养活的概率
delta = 0.002  # 捕食者自然死亡率

# 设置模型初值
prey = 1000
predator = 100

# 设置时间步长和模拟时间
dt = 0.01
t = 100

# 模拟模型
prey_list = [prey]
predator_list = [predator]
time_list = [0]

for i in range(int(t/dt)):
    prey_new = prey + alpha*prey*dt - beta*prey*predator*dt
    predator_new = predator + gamma*beta*prey*predator*dt - delta*predator*dt
    prey = prey_new if prey_new > 0 else 0
    predator = predator_new if predator_new > 0 else 0
    prey_list.append(prey)
    predator_list.append(predator)
    time_list.append((i+1)*dt)

# 绘制图像
plt.plot(time_list, prey_list, label='prey')
plt.plot(time_list, predator_list, label='predator')
plt.legend()
plt.xlabel('time')
plt.ylabel('population')
plt.show()

该代码使用Euler方法对模型进行了数值模拟，并绘制了猎物和捕食者的数量随时间变化的图像。