lotka-volterra种间竞争模型代码

下面是一个简单的Lotka-Volterra种间竞争模型的Python代码：

import numpy as np
import matplotlib.pyplot as plt

# 定义Lotka-Volterra模型：
def lotka_volterra(x, t, a, b, c, d):
    dx = x[0] * (a - b * x[1])
    dy = -x[1] * (c - d * x[0])
    return np.array([dx, dy])

# 定义初始条件：
x0 = np.array([10, 5])

# 定义模型参数：
a = 1.5
b = 1
c = 3
d = 1

# 定义时间步长和时间范围：
t = np.linspace(0, 20, 1000)

# 使用odeint求解微分方程：
from scipy.integrate import odeint
x = odeint(lotka_volterra, x0, t, args=(a, b, c, d))

# 绘制结果：
plt.plot(t, x[:,0], label='Prey')
plt.plot(t, x[:,1], label='Predator')
plt.xlabel('Time')
plt.ylabel('Population')
plt.title('Lotka-Volterra Model')
plt.legend()
plt.show()

在这个代码中，我们使用SciPy库中的odeint函数来求解微分方程，绘制了模型中捕食者和猎物的数量随时间变化的图像。