数字商品
教育 / 知识
2024年美国大学生数学建模竞赛（A题） 建模解析，小鹿学长带队指引全代码文章与思路
小鹿学长与Tech · 2 · 0
商品介绍
我是小鹿学长，就读于上海交通大学，截至目前已经帮200+人完成了建模与思路的构建的处理了～
这回带大家体验一下2024年美国大学生数学建模竞赛呀！
我是小鹿学长，就读于上海交通大学，截至目前已经帮200+人完成了建模与思路的构建的处理了～
## 问题重述
2024 MCM Problem A: 资源可用性与性别比例
1. **性别比例的变异和生态系统影响：**
- 当七鳃鳗的种群具有改变性别比例的能力时，这对更大的生态系统会产生何种影响？
- 探讨性别比例的变化如何在生态系统中引起涟漪效应，是否对其他物种和生态平衡造成影响。
2. **七鳃鳗种群的优势和劣势：**
- 就七鳃鳗种群能够调整性别比例的能力而言，这一特性带来了哪些优势和劣势？
- 分析种群在不同环境条件下性别比例变化对其繁衍、生存和生态交互的影响。
3. **性别比例变化对生态系统稳定性的影响：**
- 随着七鳃鳗性别比例的变化，这对生态系统的稳定性产生何种影响？
- 考察性别比例变异是否可能引起生态系统中的动荡，以及这种动荡可能如何影响其他生物的生存条件。
4. **可变性别比例对其他生物的影响：**
- 具有可变性别比例的七鳃鳗种群是否能够为生态系统中的其他生物提供优势，尤其是对寄生虫这类生物？
- 分析性别比例变异是否导致了一种互惠关系，其中七鳃鳗的性别比例变化对其他生物的繁衍和生存产生积极影响。
## 问题一
--综合生态系统动力学模型和Agent-Based Modeling (ABM)
**生态系统动力学模型介绍：**
生态系统动力学模型是一种基于差分方程或微分方程的数学模型，用于描述生态系统中物种之间的相互作用和种群动态的变化。这种模型通常涉及到种群的增长、死亡、竞争、捕食等基本生态过程。生态系统动力学模型可以帮助我们理解生物群体在特定环境条件下是如何演变的，以及它们之间的相互关系。
**Agent-Based Modeling (ABM)介绍：**
代理模型是一种计算模型，它通过对系统中个体的行为进行建模，来模拟整个系统的行为。在代理模型中，每个个体被视为一个代理，拥有自己的属性、状态和行为规则。这些个体代理根据环境中的变化做出决策，与其他代理进行交互，并在模拟中演化。ABM能够更细致地考虑个体层面的异质性和相互作用。
**为什么要结合使用这两种模型：**
1. **捕捉个体层面的行为：**
- ABM能够更好地捕捉个体层面的异质性和复杂行为，而生态系统动力学模型则更适合描述整体种群动态。结合这两者可以在模型中兼顾个体和整体的行为。
2. **考虑环境的影响：**
- ABM能够很好地模拟个体对环境的感知和响应，而生态系统动力学模型可以更全面地考虑环境因素对整体种群的影响。综合考虑这两者可以更全面地理解环境与个体、种群之间的复杂关系。
3. **探究动态的传播效应：**
- 生态系统动力学模型能够很好地捕捉群体之间的动态传播效应，而ABM则能够更详细地揭示个体行为对这种传播效应的贡献。综合使用有助于深入研究这些效应如何在个体行为层面传播。
4. **提高模型的可解释性：**
- ABM通常更容易解释，因为它涉及到具体的个体行为。将ABM与生态系统动力学模型结合使用，可以提高整个模型的可解释性，使研究者更容易理解模型的输出。
5. **更全面的理解生态系统的复杂性：**
- 生态系统通常涉及到多个层面的复杂相互作用，包括个体行为、物种之间的相互作用以及环境的影响。结合使用ABM和生态系统动力学模型，可以更全面地理解这些层面之间的交互作用，提供更深刻的见解。
**综合生态系统动力学模型和Agent-Based Modeling (ABM)的解决思路：**
**1. 生态系统动力学模型的基础：**
- **差分方程模型：**
$$\frac{dM}{dt} = r_M \cdot M \cdot (1 - \frac{M}{K_M}) - \alpha_{MF} \cdot M \cdot F$$
$$\frac{dF}{dt} = r_F \cdot F \cdot (1 - \frac{F}{K_F}) - \alpha_{FM} \cdot F \cdot M$$
- 这里，$M$和$F$分别是雄性和雌性的种群密度，$r_M$和$r_F$是其增长率，$K_M$和$K_F$是种群的饱和容量，$\alpha_{MF}$和$\alpha_{FM}$是雄性和雌性之间的相互作用系数。
**2. 引入Agent-Based Modeling (ABM)来模拟个体层面的交互：**
- **代理模型中的个体行为：**
- 个体代理具有属性，如性别、年龄、繁殖状态等。
- 个体代理能够感知环境中的食物可用性、温度等因素。
- **个体代理的决策行为：**
- 个体代理基于感知到的环境因素作出决策，包括繁殖、寻找伴侣等行为。
- 繁殖决策可能受到性别比例、食物可用性等因素的影响。
- **与生态系统动力学模型的耦合：**
- 个体代理的繁殖行为和生态系统动力学模型中的种群动态相互影响。
- **性别比例的变异嵌入到个体行为中：**
- 个体代理的繁殖决策可能受到性别比例的变异的影响。
- 可以使用以下公式：
$$\text{Sex Ratio} = \frac{M}{M + F}$$
**3. 个体与种群的关系：**
- **生态系统动力学模型中的其他生物：**
$$\frac{dP}{dt} = r_P \cdot P \cdot (1 - \frac{P}{K_P}) - \beta_{PM} \cdot P \cdot M$$
- 这里，$P$是其他生物的种群密度，$r_P$和$K_P$是其增长率和饱和容量，$\beta_{PM}$是七鳃鳗对其他生物的影响系数。
- **ABM中的个体与其他生物的交互：**
- 个体代理的行为可能影响其他生物的种群动态，如寄生虫在七鳃鳗身上的寄生等。
**4. 环境因素的耦合：**
- **ABM中的环境因素：**
- 模拟食物可用性、温度等环境因素，这些因素会影响个体代理的行为，如寻找食物、选择繁殖地点等。
- **环境因素与生态系统动力学模型的耦合：**
- 环境因素在生态系统动力学模型中影响七鳃鳗的种群动态。
```python
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp
from mesa import Agent, Model
from mesa.time import RandomActivation
from mesa.space import MultiGrid
from mesa.datacollection import DataCollector
# 定义生态系统动力学模型
def ecosystem_dynamics(t, y, r_M, r_F, alpha_MF, alpha_FM):
M, F = y
dMdt = r_M * M * (1 - M / (M + F)) - alpha_MF * M * F
dFdt = r_F * F * (1 - F / (M + F)) - alpha_FM * F * M
return [dMdt, dFdt]
# 定义代理模型中的个体代理
class LampreyAgent(Agent):
def __init__(self, unique_id, model):
super().__init__(unique_id, model)
self.sex = np.random.choice(['Male', 'Female'])
def step(self):
# 个体行为规则，例如繁殖、移动等
pass
# 定义代理模型
#见完整版
# 运行生态系统动力学模型
r_M = 0.1
r_F = 0.1
alpha_MF = 0.02
alpha_FM = 0.02
y0 = [100, 100] # 初始种群密度
sol = solve_ivp(
ecosystem_dynamics,
[0, 100],
y0,
args=(r_M, r_F, alpha_MF, alpha_FM),
dense_output=True
)
# 绘制种群动态图
time_points = np.linspace(0, 100, 1000)
population_evolution = sol.sol(time_points)
plt.figure(figsize=(10, 6))
plt.plot(time_points, population_evolution[0], label='Male Population')
plt.plot(time_points, population_evolution[1], label='Female Population')
plt.xlabel('Time')
plt.ylabel('Population')
plt.legend()
plt.title('Population Dynamics of Sea Lampreys')
plt.show()
# 运行代理模型
lamprey_model = LampreyModel(N=100, width=10, height=10)
for i in range(100):
lamprey_model.step()
# 数据收集
agent_sex_data = lamprey_model.datacollector.get_agent_vars_dataframe()
print(agent_sex_data.head())
```
在上述示例中，使用了 `matplotlib` 库进行数据可视化。以下是如何使用 `matplotlib` 可视化代理模型中个体的性别分布的代码。
```python
# 在 LampreyModel 类中添加一个可视化方法
class LampreyModel(Model):
# ... (之前的代码)
def visualize_agents(self):
agent_counts = self.datacollector.get_agent_counts()
male_count = agent_counts['Male']
female_count = agent_counts['Female']
labels = ['Male', 'Female']
counts = [male_count, female_count]
plt.figure(figsize=(6, 4))
plt.bar(labels, counts, color=['blue', 'pink'])
plt.xlabel('Sex')
plt.ylabel('Count')
plt.title('Distribution of Lamprey Agents')
plt.show()
# 运行代理模型和可视化
lamprey_model = LampreyModel(N=100, width=10, height=10)
for i in range(100):
lamprey_model.step()
# 可视化个体性别分布
lamprey_model.visualize_agents()
```
上述代码添加了 `visualize_agents` 方法，该方法使用 `matplotlib` 绘制了一个柱状图，显示了代理模型中雄性和雌性个体的数量分布。
## 问题二
--结合蒙特卡洛模拟和Agent-Based Modeling (ABM)
问题. **七鳃鳗种群的优势和劣势：**
- 就七鳃鳗种群能够调整性别比例的能力而言，这一特性带来了哪些优势和劣势？
- 分析种群在不同环境条件下性别比例变化对其繁衍、生存和生态交互的影响。
方法介绍：
**Agent-Based Modeling (ABM)：**
Agent-Based Modeling (ABM) 是一种建模和仿真方法，其中个体被建模为独立的代理，这些代理具有自己的属性、状态和行为规则。这些代理通过模拟它们之间的相互作用，以及它们对环境的响应，从而形成整体的系统行为。ABM 是一种强大的方法，特别适用于模拟个体行为对整体系统行为的影响。在问题二中，使用ABM的原因包括：
- **个体行为模拟：** ABM 允许模拟每个七鳃鳗个体的行为，包括繁殖、移动、食物获取等。这对于研究性别比例变化对个体层面的生态交互影响至关重要。
- **动态适应性：** 七鳃鳗种群的性别比例调整可能是一种动态的适应策略。ABM 可以帮助模拟种群在不同环境条件下如何动态调整性别比例以适应外部变化。
- **模拟不确定性：** ABM 允许引入随机性，模拟不同个体之间的差异和环境的变化。这有助于考虑到不确定性，使模型更贴近实际情况。
**Monte Carlo Simulation：**
蒙特卡洛模拟是一种通过随机抽样来近似数值结果的方法。在蒙特卡洛模拟中，通过多次模拟系统的随机状态变化，以获得对系统性质的概率性估计。在问题二中，使用蒙特卡洛模拟的原因包括：
- **不确定性的建模：** 蒙特卡洛模拟允许引入不确定性，通过多次运行模拟来模拟参数的分布和可能的情景，从而更全面地考虑不确定性。
- **参数范围的探索：** 蒙特卡洛模拟能够在给定的参数范围内对模型进行大量的随机抽样。这对于研究不同环境条件下性别比例变化对生态系统的影响至关重要。
- **统计分析：** 蒙特卡洛模拟提供了对模拟结果的统计分析能力，如均值、方差等。这有助于确定结果的可靠性和置信区间。
综合使用ABM和蒙特卡洛模拟，可以更全面地考虑个体层面的动态行为和系统整体性质之间的相互作用，同时能够处理模型中的不确定性，提高对性别比例变化对生态系统影响的理解。
**思路展开：**
1. **建立代理模型（ABM）：**
- 定义代理模型中个体代理的行为规则，包括繁殖、移动、食物获取等。
- 引入性别比例的变化作为模型的可调参数。
- 使用概率分布来模拟不同环境条件下个体的性别比例。
2. **确定评估指标：**
- 定义用于评估生态系统影响的指标，例如生态系统的稳定性、其他生物种群的变化、种群繁衍率等。
3. **设定实验参数范围：**
- 确定需要研究的参数范围，包括性别比例的初始值、变化速率等。
- 对每个参数组合运行一定数量的模拟实验，引入蒙特卡洛模拟的随机性。
4. **蒙特卡洛模拟（Monte Carlo Simulation）：**
- 针对每个参数组合，多次运行代理模型，每次运行引入一定程度的随机性。
- 收集每次模拟的结果，包括生态系统指标和性别比例的变化。
5. **分析结果：**
- 对模拟结果进行统计分析，计算生态系统指标的平均值、方差等。
- 观察不同性别比例变化对生态系统的影响，特别是在不同的参数组合下。
6. **性别比例变化的影响分析：**
- 比较不同性别比例下种群的繁衍率和生存率。
- 观察性别比例变化对生态系统中其他生物种群的影响，例如掠食者和被掠食者之间的相互作用。
7. **制作可视化图表：**
- 使用可视化工具绘制图表，展示性别比例变化和生态系统指标之间的关系。
- 绘制散点图、折线图等，以直观地展示模拟结果。
8. **敏感性分析：**
- 进行敏感性分析，评估不同参数对模拟结果的影响。
- 确定哪些参数对于性别比例变化和生态系统影响最为敏感。
9. **提出结论：**
- 根据模拟结果和分析，提出关于性别比例变化对生态系统的优劣势的结论。
- 考虑不确定性和模型的局限性，明确研究的结论范围。
**回答问题二的子问题：**
**优势：**
- *繁衍适应性：* 通过调整性别比例，种群能够更好地适应环境变化，提高繁衍的适应性。
- *生存优势：* 在不同环境条件下，调整性别比例可能提供生存优势，使得种群更容易应对资源的变化。
**劣势：**
- *生态平衡破坏：* 过度的性别比例调整可能破坏生态平衡，影响生态系统的稳定性。
- *与其他物种互动变化：* 种群性别比例的变化可能对与其他物种的生态交互产生不可预测的变化，可能对掠食者-被掠食者关系等产生负面影响。
```python
import numpy as np
from mesa import Agent, Model
from mesa.time import RandomActivation
from mesa.space import MultiGrid
from mesa.datacollection import DataCollector
import matplotlib.pyplot as plt
class LampreyAgent(Agent):
def __init__(self, unique_id, model, sex):
super().__init__(unique_id, model)
self.sex = sex
def step(self):
# 个体行为规则，例如繁殖、移动等
pass
class LampreyModel(Model):
def __init__(self, N, width, height, initial_male_ratio):
self.num_agents = N
self.grid = MultiGrid(width, height, True)
self.schedule = RandomActivation(self)
# 创建代理，初始时按照给定的性别比例创建
num_initial_males = int(N * initial_male_ratio)
for i in range(N):
sex = 'Male' if i < num_initial_males else 'Female'
a = LampreyAgent(i, self, sex)
self.schedule.add(a)
x = self.random.randrange(self.grid.width)
y = self.random.randrange(self.grid.height)
self.grid.place_agent(a, (x, y))
#见完整版
# 数据收集
agent_sex_data = lamprey_model.datacollector.get_agent_vars_dataframe()
return agent_sex_data
# 蒙特卡洛模拟，运行多次模型并统计结果
num_simulations = 10
results = []
for _ in range(num_simulations):
initial_male_ratio = np.random.uniform(0.4, 0.6) # 随机初始化不同的性别比例
simulation_result = run_model(initial_male_ratio)
results.append(simulation_result)
# 绘制结果的柱状图
plt.figure(figsize=(12, 6))
for idx, result in enumerate(results):
plt.bar(idx, result['Sex'].value_counts()['Male'], color='blue', alpha=0.5)
plt.bar(idx, result['Sex'].value_counts()['Female'], bottom=result['Sex'].value_counts()['Male'], color='pink', alpha=0.5)
plt.xlabel('Simulation')
plt.ylabel('Population Count')
plt.title('Population Distribution in Different Simulations')
plt.show()
```
当涉及到生态系统的模拟和性别比例变化时，可以考虑绘制多种图表来更全面地呈现模拟结果。
1. **时间序列图：** 绘制种群性别比例随时间的变化，以观察其动态演化。
2. **散点图：** 在二维平面上绘制个体的分布，可以用不同颜色或形状表示不同性别。
3. **箱线图：** 展示不同模拟中性别比例的分布情况，包括中位数、四分位数等统计信息。
4. **生态网络图：** 如果模型考虑到了七鳃鳗与其他生物的相互作用，可以绘制生态网络图展示物种之间的关系。
5. **密度图：** 在空间上绘制个体的密度分布，观察不同区域内性别比例的差异。
6. **参数敏感性分析图：** 绘制不同参数对模型输出的敏感性，以观察模型对参数变化的响应。
7. **生存曲线：** 如果模型中考虑到了个体的寿命和生存率，可以绘制生存曲线来展示个体的寿命分布。
## 问题三
--ABM结合LampreyAgent模型
**建模思路：**
1. **建立 Agent-Based Model (ABM)：**
- 定义 LampreyAgent 代理，模拟七鳃鳗个体的行为，包括繁殖、移动等。
- 引入其他生物代理，模拟它们的行为，如被掠食者等。
2. **构建生态网络模型：**
- 将七鳃鳗代理和其他生物代理表示为网络节点。
- 引入边表示生态交互关系，包括食物链、掠食关系等。
3. **考虑性别比例的动态变化：**
- 在 LampreyAgent 中添加性别比例的动态变化规则，考虑食物、环境等因素的影响。
- 在网络模型中，将七鳃鳗性别比例的变化作为节点属性。
4. **生态网络动态更新：**
- 使用 Lotka-Volterra 模型描述捕食者（七鳃鳗）和被捕食者（其他生物）的数量变化。
- 引入性别比例变化的方程，考虑七鳃鳗个体繁殖和性别比例的动态变化。
5. **性别比例变化公式：**
- 使用基于环境因素和个体行为的公式，计算七鳃鳗性别比例的变化率。例如：
$$\frac{dM}{dt} = \gamma \cdot M \cdot \left(1 - \frac{M}{K_M}\right) + \delta \cdot H \cdot P$$
- 这个函数 $f$ 应考虑到环境的影响，例如食物充足时性别比例趋于平衡，食物稀缺时可能导致雄性或雌性占优势。
6. **生态系统稳定性评估：**
- 定义生态系统稳定性的评估指标，考虑到七鳃鳗群体和其他生物群体的相互作用。
- 使用节点度的方差、网络连通性等指标来评估生态系统的稳定性。
7. **模拟实验和统计分析：**
- 运行多次模拟实验，观察七鳃鳗性别比例变化对生态系统稳定性的影响。
- 使用统计方法对模拟结果进行分析，包括平均值、方差等。
8. **敏感性分析：**
- 进行敏感性分析，评估不同参数对模型输出和生态系统稳定性的影响。
- 调整模型参数，观察模型响应的变化。
以上建模思路是一个综合考虑了 ABM 和 Lotka-Volterra 模型的框架，可以进一步细化和扩展，根据具体需求调整参数和模型结构。
**使用 Lotka-Volterra 模型的思路：**
Lotka-Volterra 模型是一种常见的捕食者-被捕食者模型，适用于描述生态系统中两个种群之间的相互作用。在这里，我们将其扩展以考虑七鳃鳗种群（捕食者）和其他生物种群（被捕食者），并加入性别比例的动态变化。
**基本方程：**
1. **捕食者种群（七鳃鳗）：**
$$\frac{dP}{dt} = r_P \cdot P \cdot \left(1 - \frac{P}{K_P}\right) - \alpha \cdot P \cdot H$$
2. **被捕食者种群（其他生物）：**
$$\frac{dH}{dt} = r_H \cdot H \cdot \left(1 - \frac{H}{K_H}\right) - \beta \cdot H \cdot P$$
**扩展考虑性别比例的动态变化：**
在上述方程的基础上，引入性别比例的动态变化：
3. **性别比例变化：**
$$\frac{dM}{dt} = \gamma \cdot M \cdot \left(1 - \frac{M}{K_M}\right) + \delta \cdot H \cdot P$$
其中：
- $P$ 是七鳃鳗种群数量，$H$ 是其他生物种群数量，$M$ 是雄性七鳃鳗数量。
- $r_P$, $r_H$ 是种群的固有增长率。
- $K_P$, $K_H$ 是种群的环境容量。
- $\alpha$, $\beta$ 是相互作用的强度。
- $\gamma$, $\delta$ 是性别比例变化的影响参数。
**模型的伪代码：**
```python
class LampreyAgent(Agent):
def __init__(self, unique_id, model, sex):
super().__init__(unique_id, model)
self.sex = sex
def step(self):
# 个体行为规则，包括繁殖、移动等
# 计算性别比例的变化率并更新
class EcologyModel(Model):
def __init__(self, N, K_P, K_H, alpha, beta, gamma, delta):
# 初始化种群数量和参数
self.num_lampreys = N
self.K_P = K_P
self.K_H = K_H
self.alpha = alpha
self.beta = beta
self.gamma = gamma
self.delta = delta
# 初始化 LampreyAgent 代理，并将其添加到模型中
for i in range(N):
sex = 'Male' if np.random.rand() < 0.5 else 'Female'
a = LampreyAgent(i, self, sex)
self.schedule.add(a)
def step(self):
# 计算捕食者和被捕食者的数量变化率
dP = self.r_P * self.num_lampreys * (1 - self.num_lampreys / self.K_P) - self.alpha * self.num_lampreys * self.num_preys
dH = self.r_H * self.num_preys * (1 - self.num_preys / self.K_H) - self.beta * self.num_preys * self.num_lampreys
# 计算性别比例的变化率
dM = self.gamma * self.num_males * (1 - self.num_males / self.K_M) + self.delta * self.num_preys * self.num_lampreys
# 更新数量
self.num_lampreys += dP
self.num_preys += dH
self.num_males += dM
# 更新每个代理的行为
self.schedule.step()
# 运行模型
ecology_model = EcologyModel(N=100, K_P=your_defined_value, K_H=your_defined_value,
alpha=your_defined_value, beta=your_defined_value,
gamma=your_defined_value, delta=your_defined_value)
for i in range(100):
ecology_model.step()
# 收集数据并分析
# ...见完整版
```
可以考虑以下一些图表来更全面地了解模型的行为和结果：
1. **网络结构图：**
- 使用 NetworkX 库将生态系统的网络结构可视化，显示七鳃鳗和其他生物之间的关系。
- 这可以帮助你理解生态系统中不同种群之间的交互。
```python
import networkx as nx
import matplotlib.pyplot as plt
# 获取模型中的网络结构
G = ecology_model.grid.G
# 绘制网络结构图
pos = nx.spring_layout(G)
nx.draw(G, pos, with_labels=True)
plt.show()
```
2. **种群分布图：**
- 绘制七鳃鳗和其他生物在生态系统中的空间分布图。
- 这有助于观察个体在环境中的分布情况。
```python
import matplotlib.pyplot as plt
# 获取模型中的代理
agents = ecology_model.schedule.agents
# 绘制空间分布图
plt.figure(figsize=(8, 8))
for agent in agents:
x, y = agent.pos
color = 'blue' if agent.sex == 'Male' else 'red'
plt.scatter(x, y, c=color)
plt.xlabel('X Coordinate')
plt.ylabel('Y Coordinate')
plt.title('Population Distribution')
plt.show()
```
## 问题四
-ABM模型
解决问题四需要考虑七鳃鳗种群的可变性别比例对生态系统中其他生物的影响。一种适合解决这类问题的方法是 Agent-Based Modeling (ABM)，结合寄生虫的动态模型。以下是一些实现思路：
**第一步：定义代理类**
在 `LampreyAgent` 类中，我们已经考虑了性别比例的动态变化。现在，我们需要扩展这个类，以包含对寄生虫的交互规则。
```python
class LampreyAgent(Agent):
def __init__(self, unique_id, model, sex):
super().__init__(unique_id, model)
self.sex = sex
def step(self):
# 个体行为规则，包括繁殖、移动等
self.reproduce()
self.move()
# 计算性别比例的变化率并更新
self.calculate_sex_ratio()
# 寄生虫交互规则
parasites = self.model.grid.get_neighbors(self.pos, moore=True, include_center=False, radius=1)
for parasite in parasites:
self.interact_with_parasite(parasite)
def reproduce(self):
# 繁殖规则
# ...
def move(self):
# 移动规则
# ...
def calculate_sex_ratio(self):
# 计算性别比例的变化率并更新
# ...
def interact_with_parasite(self, parasite):
# 寄生虫相互作用规则
# ...
```
**第二步：创建寄生虫代理类**
在上述代码中，我们已经引入了 `ParasiteAgent` 类，现在需要在这个类中定义寄生虫的行为规则。
```python
class ParasiteAgent(Agent):
def __init__(self, unique_id, model):
super().__init__(unique_id, model)
def step(self):
# 寄生虫的行为规则，包括繁殖、寄生等
self.reproduce()
self.parasitize()
def reproduce(self):
# 寄生虫繁殖规则
# ...
def parasitize(self):
# 寄生规则
# ...
```
**第三步：扩展生态系统模型**
在 `EcologyModel` 类中，我们需要更新模型的步进规则，考虑寄生虫数量和性别比例的动态变化。我们将采用 Lotka-Volterra 模型的扩展来描述七鳃鳗和寄生虫之间的相互作用。
```python
class EcologyModel(Model):
def __init__(self, N, N_parasites, K_P, K_H, alpha, beta, gamma, delta, epsilon, zeta):
# ...
# 添加寄生虫代理
for i in range(N_parasites):
parasite = ParasiteAgent(i, self)
self.schedule.add(parasite)
x, y = np.random.random(), np.random.random()
self.grid.place_agent(parasite, (x, y))
def step(self):
# 计算捕食者和被捕食者的数量变化率
dP = self.alpha * self.num_lampreys * (1 - self.num_lampreys / self.K_P) - self.beta * self.num_lampreys * self.num_preys
# 计算性别比例的变化率
dM = self.gamma * self.num_males * (1 - self.num_males / self.K_M) + self.delta * self.num_preys * self.num_lampreys
# 计算寄生虫数量的变化率
dParasites = self.epsilon * self.num_parasites * (1 - self.num_parasites / self.K_P) - self.zeta * self.num_parasites * self.num_lampreys
# 更新数量
self.num_lampreys += dP
self.num_males += dM
self.num_parasites += dParasites
# 更新每个代理的行为
self.schedule.step()
# 数据收集
self.datacollector.collect(self)
```
**第四步：定义寄生虫与七鳃鳗的相互作用规则**
在 `LampreyAgent` 类中，我们添加了 `interact_with_parasite` 方法，该方法描述了七鳃鳗与寄生虫之间的相互作用规则。
```python
class LampreyAgent(Agent):
# ...
def interact_with_parasite(self, parasite):
# 模拟七鳃鳗与寄生虫的相互作用
# 可以根据具体情况定义交互规则，如寄生虫对性别比例的影响
# 以下是一个简单的示例
if self.sex == 'Male':
self.model.num_preys += 1 # 七鳃鳗捕食寄生虫，增加被捕食者数量
self.model.num_parasites -= 1 # 寄生虫数量减少
```
对于问题四，即“可变性别比例对其他生物的影响”，可以通过模拟实验和数据分析绘制多种图表来展示结果。以下是一些可能的图表：
1. **性别比例变化趋势图：**
- 绘制七鳃鳗种群的总体性别比例随时间的变化趋势，以展示模型中性别比例的动态变化。
2. **寄生虫数量变化趋势图：**
- 绘制寄生虫数量随时间的变化趋势，以展示七鳃鳗种群可变性别比例对寄生虫数量的影响。
3. **其他生物数量变化趋势图：**
- 如果模型中包含其他生物，例如食物链中的其他物种，可以绘制它们的数量随时间的变化趋势图，以展示可变性别比例对整个生态系统的影响。
4. **七鳃鳗与寄生虫相互作用图：**
- 绘制七鳃鳗与寄生虫之间相互作用的图表，以展示不同性别比例下它们之间的动态关系。
5. **生态系统稳定性分析图：**
- 分析模型中生态系统的稳定性，可以绘制相对稳定性指标随时间的变化图表，以了解可变性别比例对生态系统整体稳定性的影响。
6. **参数敏感性分析图：**
- 进行参数敏感性分析，绘制不同模型参数下性别比例、寄生虫数量等指标的变化图表，以确定模型中哪些参数对结果影响较大。
7. **可变性别比例优势和劣势分析图：**
- 绘制不同环境条件下，可变性别比例对七鳃鳗种群的优势和劣势的分析图表，以了解在不同情境下可变性别比例的影响。
8. **性别比例与其他生物种群相互作用图：**
- 如果模型中包含其他生物，可以绘制它们与七鳃鳗性别比例之间的相互作用图表，以揭示可变性别比例对生态系统中其他生物的影响。
这回带大家体验一下2024年美国大学生数学建模竞赛呀！
我是小鹿学长，就读于上海交通大学，截至目前已经帮200+人完成了建模与思路的构建的处理了～
## 问题重述
2024 MCM Problem A: 资源可用性与性别比例
1. **性别比例的变异和生态系统影响：**
- 当七鳃鳗的种群具有改变性别比例的能力时，这对更大的生态系统会产生何种影响？
- 探讨性别比例的变化如何在生态系统中引起涟漪效应，是否对其他物种和生态平衡造成影响。
2. **七鳃鳗种群的优势和劣势：**
- 就七鳃鳗种群能够调整性别比例的能力而言，这一特性带来了哪些优势和劣势？
- 分析种群在不同环境条件下性别比例变化对其繁衍、生存和生态交互的影响。
3. **性别比例变化对生态系统稳定性的影响：**
- 随着七鳃鳗性别比例的变化，这对生态系统的稳定性产生何种影响？
- 考察性别比例变异是否可能引起生态系统中的动荡，以及这种动荡可能如何影响其他生物的生存条件。
4. **可变性别比例对其他生物的影响：**
- 具有可变性别比例的七鳃鳗种群是否能够为生态系统中的其他生物提供优势，尤其是对寄生虫这类生物？
- 分析性别比例变异是否导致了一种互惠关系，其中七鳃鳗的性别比例变化对其他生物的繁衍和生存产生积极影响。
## 问题一
--综合生态系统动力学模型和Agent-Based Modeling (ABM)
**生态系统动力学模型介绍：**
生态系统动力学模型是一种基于差分方程或微分方程的数学模型，用于描述生态系统中物种之间的相互作用和种群动态的变化。这种模型通常涉及到种群的增长、死亡、竞争、捕食等基本生态过程。生态系统动力学模型可以帮助我们理解生物群体在特定环境条件下是如何演变的，以及它们之间的相互关系。
**Agent-Based Modeling (ABM)介绍：**
代理模型是一种计算模型，它通过对系统中个体的行为进行建模，来模拟整个系统的行为。在代理模型中，每个个体被视为一个代理，拥有自己的属性、状态和行为规则。这些个体代理根据环境中的变化做出决策，与其他代理进行交互，并在模拟中演化。ABM能够更细致地考虑个体层面的异质性和相互作用。
**为什么要结合使用这两种模型：**
1. **捕捉个体层面的行为：**
- ABM能够更好地捕捉个体层面的异质性和复杂行为，而生态系统动力学模型则更适合描述整体种群动态。结合这两者可以在模型中兼顾个体和整体的行为。
2. **考虑环境的影响：**
- ABM能够很好地模拟个体对环境的感知和响应，而生态系统动力学模型可以更全面地考虑环境因素对整体种群的影响。综合考虑这两者可以更全面地理解环境与个体、种群之间的复杂关系。
3. **探究动态的传播效应：**
- 生态系统动力学模型能够很好地捕捉群体之间的动态传播效应，而ABM则能够更详细地揭示个体行为对这种传播效应的贡献。综合使用有助于深入研究这些效应如何在个体行为层面传播。
4. **提高模型的可解释性：**
- ABM通常更容易解释，因为它涉及到具体的个体行为。将ABM与生态系统动力学模型结合使用，可以提高整个模型的可解释性，使研究者更容易理解模型的输出。
5. **更全面的理解生态系统的复杂性：**
- 生态系统通常涉及到多个层面的复杂相互作用，包括个体行为、物种之间的相互作用以及环境的影响。结合使用ABM和生态系统动力学模型，可以更全面地理解这些层面之间的交互作用，提供更深刻的见解。
**综合生态系统动力学模型和Agent-Based Modeling (ABM)的解决思路：**
**1. 生态系统动力学模型的基础：**
- **差分方程模型：**
$$\frac{dM}{dt} = r_M \cdot M \cdot (1 - \frac{M}{K_M}) - \alpha_{MF} \cdot M \cdot F$$
$$\frac{dF}{dt} = r_F \cdot F \cdot (1 - \frac{F}{K_F}) - \alpha_{FM} \cdot F \cdot M$$
- 这里，$M$和$F$分别是雄性和雌性的种群密度，$r_M$和$r_F$是其增长率，$K_M$和$K_F$是种群的饱和容量，$\alpha_{MF}$和$\alpha_{FM}$是雄性和雌性之间的相互作用系数。
**2. 引入Agent-Based Modeling (ABM)来模拟个体层面的交互：**
- **代理模型中的个体行为：**
- 个体代理具有属性，如性别、年龄、繁殖状态等。
- 个体代理能够感知环境中的食物可用性、温度等因素。
- **个体代理的决策行为：**
- 个体代理基于感知到的环境因素作出决策，包括繁殖、寻找伴侣等行为。
- 繁殖决策可能受到性别比例、食物可用性等因素的影响。
- **与生态系统动力学模型的耦合：**
- 个体代理的繁殖行为和生态系统动力学模型中的种群动态相互影响。
- **性别比例的变异嵌入到个体行为中：**
- 个体代理的繁殖决策可能受到性别比例的变异的影响。
- 可以使用以下公式：
$$\text{Sex Ratio} = \frac{M}{M + F}$$
**3. 个体与种群的关系：**
- **生态系统动力学模型中的其他生物：**
$$\frac{dP}{dt} = r_P \cdot P \cdot (1 - \frac{P}{K_P}) - \beta_{PM} \cdot P \cdot M$$
- 这里，$P$是其他生物的种群密度，$r_P$和$K_P$是其增长率和饱和容量，$\beta_{PM}$是七鳃鳗对其他生物的影响系数。
- **ABM中的个体与其他生物的交互：**
- 个体代理的行为可能影响其他生物的种群动态，如寄生虫在七鳃鳗身上的寄生等。
**4. 环境因素的耦合：**
- **ABM中的环境因素：**
- 模拟食物可用性、温度等环境因素，这些因素会影响个体代理的行为，如寻找食物、选择繁殖地点等。
- **环境因素与生态系统动力学模型的耦合：**
- 环境因素在生态系统动力学模型中影响七鳃鳗的种群动态。
```python
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp
from mesa import Agent, Model
from mesa.time import RandomActivation
from mesa.space import MultiGrid
from mesa.datacollection import DataCollector
# 定义生态系统动力学模型
def ecosystem_dynamics(t, y, r_M, r_F, alpha_MF, alpha_FM):
M, F = y
dMdt = r_M * M * (1 - M / (M + F)) - alpha_MF * M * F
dFdt = r_F * F * (1 - F / (M + F)) - alpha_FM * F * M
return [dMdt, dFdt]
# 定义代理模型中的个体代理
class LampreyAgent(Agent):
def __init__(self, unique_id, model):
super().__init__(unique_id, model)
self.sex = np.random.choice(['Male', 'Female'])
def step(self):
# 个体行为规则，例如繁殖、移动等
pass
# 定义代理模型
#见完整版
# 运行生态系统动力学模型
r_M = 0.1
r_F = 0.1
alpha_MF = 0.02
alpha_FM = 0.02
y0 = [100, 100] # 初始种群密度
sol = solve_ivp(
ecosystem_dynamics,
[0, 100],
y0,
args=(r_M, r_F, alpha_MF, alpha_FM),
dense_output=True
)
# 绘制种群动态图
time_points = np.linspace(0, 100, 1000)
population_evolution = sol.sol(time_points)
plt.figure(figsize=(10, 6))
plt.plot(time_points, population_evolution[0], label='Male Population')
plt.plot(time_points, population_evolution[1], label='Female Population')
plt.xlabel('Time')
plt.ylabel('Population')
plt.legend()
plt.title('Population Dynamics of Sea Lampreys')
plt.show()
# 运行代理模型
lamprey_model = LampreyModel(N=100, width=10, height=10)
for i in range(100):
lamprey_model.step()
# 数据收集
agent_sex_data = lamprey_model.datacollector.get_agent_vars_dataframe()
print(agent_sex_data.head())
```
在上述示例中，使用了 `matplotlib` 库进行数据可视化。以下是如何使用 `matplotlib` 可视化代理模型中个体的性别分布的代码。
```python
# 在 LampreyModel 类中添加一个可视化方法
class LampreyModel(Model):
# ... (之前的代码)
def visualize_agents(self):
agent_counts = self.datacollector.get_agent_counts()
male_count = agent_counts['Male']
female_count = agent_counts['Female']
labels = ['Male', 'Female']
counts = [male_count, female_count]
plt.figure(figsize=(6, 4))
plt.bar(labels, counts, color=['blue', 'pink'])
plt.xlabel('Sex')
plt.ylabel('Count')
plt.title('Distribution of Lamprey Agents')
plt.show()
# 运行代理模型和可视化
lamprey_model = LampreyModel(N=100, width=10, height=10)
for i in range(100):
lamprey_model.step()
# 可视化个体性别分布
lamprey_model.visualize_agents()
```
上述代码添加了 `visualize_agents` 方法，该方法使用 `matplotlib` 绘制了一个柱状图，显示了代理模型中雄性和雌性个体的数量分布。
## 问题二
--结合蒙特卡洛模拟和Agent-Based Modeling (ABM)
问题. **七鳃鳗种群的优势和劣势：**
- 就七鳃鳗种群能够调整性别比例的能力而言，这一特性带来了哪些优势和劣势？
- 分析种群在不同环境条件下性别比例变化对其繁衍、生存和生态交互的影响。
方法介绍：
**Agent-Based Modeling (ABM)：**
Agent-Based Modeling (ABM) 是一种建模和仿真方法，其中个体被建模为独立的代理，这些代理具有自己的属性、状态和行为规则。这些代理通过模拟它们之间的相互作用，以及它们对环境的响应，从而形成整体的系统行为。ABM 是一种强大的方法，特别适用于模拟个体行为对整体系统行为的影响。在问题二中，使用ABM的原因包括：
- **个体行为模拟：** ABM 允许模拟每个七鳃鳗个体的行为，包括繁殖、移动、食物获取等。这对于研究性别比例变化对个体层面的生态交互影响至关重要。
- **动态适应性：** 七鳃鳗种群的性别比例调整可能是一种动态的适应策略。ABM 可以帮助模拟种群在不同环境条件下如何动态调整性别比例以适应外部变化。
- **模拟不确定性：** ABM 允许引入随机性，模拟不同个体之间的差异和环境的变化。这有助于考虑到不确定性，使模型更贴近实际情况。
**Monte Carlo Simulation：**
蒙特卡洛模拟是一种通过随机抽样来近似数值结果的方法。在蒙特卡洛模拟中，通过多次模拟系统的随机状态变化，以获得对系统性质的概率性估计。在问题二中，使用蒙特卡洛模拟的原因包括：
- **不确定性的建模：** 蒙特卡洛模拟允许引入不确定性，通过多次运行模拟来模拟参数的分布和可能的情景，从而更全面地考虑不确定性。
- **参数范围的探索：** 蒙特卡洛模拟能够在给定的参数范围内对模型进行大量的随机抽样。这对于研究不同环境条件下性别比例变化对生态系统的影响至关重要。
- **统计分析：** 蒙特卡洛模拟提供了对模拟结果的统计分析能力，如均值、方差等。这有助于确定结果的可靠性和置信区间。
综合使用ABM和蒙特卡洛模拟，可以更全面地考虑个体层面的动态行为和系统整体性质之间的相互作用，同时能够处理模型中的不确定性，提高对性别比例变化对生态系统影响的理解。
**思路展开：**
1. **建立代理模型（ABM）：**
- 定义代理模型中个体代理的行为规则，包括繁殖、移动、食物获取等。
- 引入性别比例的变化作为模型的可调参数。
- 使用概率分布来模拟不同环境条件下个体的性别比例。
2. **确定评估指标：**
- 定义用于评估生态系统影响的指标，例如生态系统的稳定性、其他生物种群的变化、种群繁衍率等。
3. **设定实验参数范围：**
- 确定需要研究的参数范围，包括性别比例的初始值、变化速率等。
- 对每个参数组合运行一定数量的模拟实验，引入蒙特卡洛模拟的随机性。
4. **蒙特卡洛模拟（Monte Carlo Simulation）：**
- 针对每个参数组合，多次运行代理模型，每次运行引入一定程度的随机性。
- 收集每次模拟的结果，包括生态系统指标和性别比例的变化。
5. **分析结果：**
- 对模拟结果进行统计分析，计算生态系统指标的平均值、方差等。
- 观察不同性别比例变化对生态系统的影响，特别是在不同的参数组合下。
6. **性别比例变化的影响分析：**
- 比较不同性别比例下种群的繁衍率和生存率。
- 观察性别比例变化对生态系统中其他生物种群的影响，例如掠食者和被掠食者之间的相互作用。
7. **制作可视化图表：**
- 使用可视化工具绘制图表，展示性别比例变化和生态系统指标之间的关系。
- 绘制散点图、折线图等，以直观地展示模拟结果。
8. **敏感性分析：**
- 进行敏感性分析，评估不同参数对模拟结果的影响。
- 确定哪些参数对于性别比例变化和生态系统影响最为敏感。
9. **提出结论：**
- 根据模拟结果和分析，提出关于性别比例变化对生态系统的优劣势的结论。
- 考虑不确定性和模型的局限性，明确研究的结论范围。
**回答问题二的子问题：**
**优势：**
- *繁衍适应性：* 通过调整性别比例，种群能够更好地适应环境变化，提高繁衍的适应性。
- *生存优势：* 在不同环境条件下，调整性别比例可能提供生存优势，使得种群更容易应对资源的变化。
**劣势：**
- *生态平衡破坏：* 过度的性别比例调整可能破坏生态平衡，影响生态系统的稳定性。
- *与其他物种互动变化：* 种群性别比例的变化可能对与其他物种的生态交互产生不可预测的变化，可能对掠食者-被掠食者关系等产生负面影响。
```python
import numpy as np
from mesa import Agent, Model
from mesa.time import RandomActivation
from mesa.space import MultiGrid
from mesa.datacollection import DataCollector
import matplotlib.pyplot as plt
class LampreyAgent(Agent):
def __init__(self, unique_id, model, sex):
super().__init__(unique_id, model)
self.sex = sex
def step(self):
# 个体行为规则，例如繁殖、移动等
pass
class LampreyModel(Model):
def __init__(self, N, width, height, initial_male_ratio):
self.num_agents = N
self.grid = MultiGrid(width, height, True)
self.schedule = RandomActivation(self)
# 创建代理，初始时按照给定的性别比例创建
num_initial_males = int(N * initial_male_ratio)
for i in range(N):
sex = 'Male' if i < num_initial_males else 'Female'
a = LampreyAgent(i, self, sex)
self.schedule.add(a)
x = self.random.randrange(self.grid.width)
y = self.random.randrange(self.grid.height)
self.grid.place_agent(a, (x, y))
#见完整版
# 数据收集
agent_sex_data = lamprey_model.datacollector.get_agent_vars_dataframe()
return agent_sex_data
# 蒙特卡洛模拟，运行多次模型并统计结果
num_simulations = 10
results = []
for _ in range(num_simulations):
initial_male_ratio = np.random.uniform(0.4, 0.6) # 随机初始化不同的性别比例
simulation_result = run_model(initial_male_ratio)
results.append(simulation_result)
# 绘制结果的柱状图
plt.figure(figsize=(12, 6))
for idx, result in enumerate(results):
plt.bar(idx, result['Sex'].value_counts()['Male'], color='blue', alpha=0.5)
plt.bar(idx, result['Sex'].value_counts()['Female'], bottom=result['Sex'].value_counts()['Male'], color='pink', alpha=0.5)
plt.xlabel('Simulation')
plt.ylabel('Population Count')
plt.title('Population Distribution in Different Simulations')
plt.show()
```
当涉及到生态系统的模拟和性别比例变化时，可以考虑绘制多种图表来更全面地呈现模拟结果。
1. **时间序列图：** 绘制种群性别比例随时间的变化，以观察其动态演化。
2. **散点图：** 在二维平面上绘制个体的分布，可以用不同颜色或形状表示不同性别。
3. **箱线图：** 展示不同模拟中性别比例的分布情况，包括中位数、四分位数等统计信息。
4. **生态网络图：** 如果模型考虑到了七鳃鳗与其他生物的相互作用，可以绘制生态网络图展示物种之间的关系。
5. **密度图：** 在空间上绘制个体的密度分布，观察不同区域内性别比例的差异。
6. **参数敏感性分析图：** 绘制不同参数对模型输出的敏感性，以观察模型对参数变化的响应。
7. **生存曲线：** 如果模型中考虑到了个体的寿命和生存率，可以绘制生存曲线来展示个体的寿命分布。
## 问题三
--ABM结合LampreyAgent模型
**建模思路：**
1. **建立 Agent-Based Model (ABM)：**
- 定义 LampreyAgent 代理，模拟七鳃鳗个体的行为，包括繁殖、移动等。
- 引入其他生物代理，模拟它们的行为，如被掠食者等。
2. **构建生态网络模型：**
- 将七鳃鳗代理和其他生物代理表示为网络节点。
- 引入边表示生态交互关系，包括食物链、掠食关系等。
3. **考虑性别比例的动态变化：**
- 在 LampreyAgent 中添加性别比例的动态变化规则，考虑食物、环境等因素的影响。
- 在网络模型中，将七鳃鳗性别比例的变化作为节点属性。
4. **生态网络动态更新：**
- 使用 Lotka-Volterra 模型描述捕食者（七鳃鳗）和被捕食者（其他生物）的数量变化。
- 引入性别比例变化的方程，考虑七鳃鳗个体繁殖和性别比例的动态变化。
5. **性别比例变化公式：**
- 使用基于环境因素和个体行为的公式，计算七鳃鳗性别比例的变化率。例如：
$$\frac{dM}{dt} = \gamma \cdot M \cdot \left(1 - \frac{M}{K_M}\right) + \delta \cdot H \cdot P$$
- 这个函数 $f$ 应考虑到环境的影响，例如食物充足时性别比例趋于平衡，食物稀缺时可能导致雄性或雌性占优势。
6. **生态系统稳定性评估：**
- 定义生态系统稳定性的评估指标，考虑到七鳃鳗群体和其他生物群体的相互作用。
- 使用节点度的方差、网络连通性等指标来评估生态系统的稳定性。
7. **模拟实验和统计分析：**
- 运行多次模拟实验，观察七鳃鳗性别比例变化对生态系统稳定性的影响。
- 使用统计方法对模拟结果进行分析，包括平均值、方差等。
8. **敏感性分析：**
- 进行敏感性分析，评估不同参数对模型输出和生态系统稳定性的影响。
- 调整模型参数，观察模型响应的变化。
以上建模思路是一个综合考虑了 ABM 和 Lotka-Volterra 模型的框架，可以进一步细化和扩展，根据具体需求调整参数和模型结构。
**使用 Lotka-Volterra 模型的思路：**
Lotka-Volterra 模型是一种常见的捕食者-被捕食者模型，适用于描述生态系统中两个种群之间的相互作用。在这里，我们将其扩展以考虑七鳃鳗种群（捕食者）和其他生物种群（被捕食者），并加入性别比例的动态变化。
**基本方程：**
1. **捕食者种群（七鳃鳗）：**
$$\frac{dP}{dt} = r_P \cdot P \cdot \left(1 - \frac{P}{K_P}\right) - \alpha \cdot P \cdot H$$
2. **被捕食者种群（其他生物）：**
$$\frac{dH}{dt} = r_H \cdot H \cdot \left(1 - \frac{H}{K_H}\right) - \beta \cdot H \cdot P$$
**扩展考虑性别比例的动态变化：**
在上述方程的基础上，引入性别比例的动态变化：
3. **性别比例变化：**
$$\frac{dM}{dt} = \gamma \cdot M \cdot \left(1 - \frac{M}{K_M}\right) + \delta \cdot H \cdot P$$
其中：
- $P$ 是七鳃鳗种群数量，$H$ 是其他生物种群数量，$M$ 是雄性七鳃鳗数量。
- $r_P$, $r_H$ 是种群的固有增长率。
- $K_P$, $K_H$ 是种群的环境容量。
- $\alpha$, $\beta$ 是相互作用的强度。
- $\gamma$, $\delta$ 是性别比例变化的影响参数。
**模型的伪代码：**
```python
class LampreyAgent(Agent):
def __init__(self, unique_id, model, sex):
super().__init__(unique_id, model)
self.sex = sex
def step(self):
# 个体行为规则，包括繁殖、移动等
# 计算性别比例的变化率并更新
class EcologyModel(Model):
def __init__(self, N, K_P, K_H, alpha, beta, gamma, delta):
# 初始化种群数量和参数
self.num_lampreys = N
self.K_P = K_P
self.K_H = K_H
self.alpha = alpha
self.beta = beta
self.gamma = gamma
self.delta = delta
# 初始化 LampreyAgent 代理，并将其添加到模型中
for i in range(N):
sex = 'Male' if np.random.rand() < 0.5 else 'Female'
a = LampreyAgent(i, self, sex)
self.schedule.add(a)
def step(self):
# 计算捕食者和被捕食者的数量变化率
dP = self.r_P * self.num_lampreys * (1 - self.num_lampreys / self.K_P) - self.alpha * self.num_lampreys * self.num_preys
dH = self.r_H * self.num_preys * (1 - self.num_preys / self.K_H) - self.beta * self.num_preys * self.num_lampreys
# 计算性别比例的变化率
dM = self.gamma * self.num_males * (1 - self.num_males / self.K_M) + self.delta * self.num_preys * self.num_lampreys
# 更新数量
self.num_lampreys += dP
self.num_preys += dH
self.num_males += dM
# 更新每个代理的行为
self.schedule.step()
# 运行模型
ecology_model = EcologyModel(N=100, K_P=your_defined_value, K_H=your_defined_value,
alpha=your_defined_value, beta=your_defined_value,
gamma=your_defined_value, delta=your_defined_value)
for i in range(100):
ecology_model.step()
# 收集数据并分析
# ...见完整版
```
可以考虑以下一些图表来更全面地了解模型的行为和结果：
1. **网络结构图：**
- 使用 NetworkX 库将生态系统的网络结构可视化，显示七鳃鳗和其他生物之间的关系。
- 这可以帮助你理解生态系统中不同种群之间的交互。
```python
import networkx as nx
import matplotlib.pyplot as plt
# 获取模型中的网络结构
G = ecology_model.grid.G
# 绘制网络结构图
pos = nx.spring_layout(G)
nx.draw(G, pos, with_labels=True)
plt.show()
```
2. **种群分布图：**
- 绘制七鳃鳗和其他生物在生态系统中的空间分布图。
- 这有助于观察个体在环境中的分布情况。
```python
import matplotlib.pyplot as plt
# 获取模型中的代理
agents = ecology_model.schedule.agents
# 绘制空间分布图
plt.figure(figsize=(8, 8))
for agent in agents:
x, y = agent.pos
color = 'blue' if agent.sex == 'Male' else 'red'
plt.scatter(x, y, c=color)
plt.xlabel('X Coordinate')
plt.ylabel('Y Coordinate')
plt.title('Population Distribution')
plt.show()
```
## 问题四
-ABM模型
解决问题四需要考虑七鳃鳗种群的可变性别比例对生态系统中其他生物的影响。一种适合解决这类问题的方法是 Agent-Based Modeling (ABM)，结合寄生虫的动态模型。以下是一些实现思路：
**第一步：定义代理类**
在 `LampreyAgent` 类中，我们已经考虑了性别比例的动态变化。现在，我们需要扩展这个类，以包含对寄生虫的交互规则。
```python
class LampreyAgent(Agent):
def __init__(self, unique_id, model, sex):
super().__init__(unique_id, model)
self.sex = sex
def step(self):
# 个体行为规则，包括繁殖、移动等
self.reproduce()
self.move()
# 计算性别比例的变化率并更新
self.calculate_sex_ratio()
# 寄生虫交互规则
parasites = self.model.grid.get_neighbors(self.pos, moore=True, include_center=False, radius=1)
for parasite in parasites:
self.interact_with_parasite(parasite)
def reproduce(self):
# 繁殖规则
# ...
def move(self):
# 移动规则
# ...
def calculate_sex_ratio(self):
# 计算性别比例的变化率并更新
# ...
def interact_with_parasite(self, parasite):
# 寄生虫相互作用规则
# ...
```
**第二步：创建寄生虫代理类**
在上述代码中，我们已经引入了 `ParasiteAgent` 类，现在需要在这个类中定义寄生虫的行为规则。
```python
class ParasiteAgent(Agent):
def __init__(self, unique_id, model):
super().__init__(unique_id, model)
def step(self):
# 寄生虫的行为规则，包括繁殖、寄生等
self.reproduce()
self.parasitize()
def reproduce(self):
# 寄生虫繁殖规则
# ...
def parasitize(self):
# 寄生规则
# ...
```
**第三步：扩展生态系统模型**
在 `EcologyModel` 类中，我们需要更新模型的步进规则，考虑寄生虫数量和性别比例的动态变化。我们将采用 Lotka-Volterra 模型的扩展来描述七鳃鳗和寄生虫之间的相互作用。
```python
class EcologyModel(Model):
def __init__(self, N, N_parasites, K_P, K_H, alpha, beta, gamma, delta, epsilon, zeta):
# ...
# 添加寄生虫代理
for i in range(N_parasites):
parasite = ParasiteAgent(i, self)
self.schedule.add(parasite)
x, y = np.random.random(), np.random.random()
self.grid.place_agent(parasite, (x, y))
def step(self):
# 计算捕食者和被捕食者的数量变化率
dP = self.alpha * self.num_lampreys * (1 - self.num_lampreys / self.K_P) - self.beta * self.num_lampreys * self.num_preys
# 计算性别比例的变化率
dM = self.gamma * self.num_males * (1 - self.num_males / self.K_M) + self.delta * self.num_preys * self.num_lampreys
# 计算寄生虫数量的变化率
dParasites = self.epsilon * self.num_parasites * (1 - self.num_parasites / self.K_P) - self.zeta * self.num_parasites * self.num_lampreys
# 更新数量
self.num_lampreys += dP
self.num_males += dM
self.num_parasites += dParasites
# 更新每个代理的行为
self.schedule.step()
# 数据收集
self.datacollector.collect(self)
```
**第四步：定义寄生虫与七鳃鳗的相互作用规则**
在 `LampreyAgent` 类中，我们添加了 `interact_with_parasite` 方法，该方法描述了七鳃鳗与寄生虫之间的相互作用规则。
```python
class LampreyAgent(Agent):
# ...
def interact_with_parasite(self, parasite):
# 模拟七鳃鳗与寄生虫的相互作用
# 可以根据具体情况定义交互规则，如寄生虫对性别比例的影响
# 以下是一个简单的示例
if self.sex == 'Male':
self.model.num_preys += 1 # 七鳃鳗捕食寄生虫，增加被捕食者数量
self.model.num_parasites -= 1 # 寄生虫数量减少
```
对于问题四，即“可变性别比例对其他生物的影响”，可以通过模拟实验和数据分析绘制多种图表来展示结果。以下是一些可能的图表：
1. **性别比例变化趋势图：**
- 绘制七鳃鳗种群的总体性别比例随时间的变化趋势，以展示模型中性别比例的动态变化。
2. **寄生虫数量变化趋势图：**
- 绘制寄生虫数量随时间的变化趋势，以展示七鳃鳗种群可变性别比例对寄生虫数量的影响。
3. **其他生物数量变化趋势图：**
- 如果模型中包含其他生物，例如食物链中的其他物种，可以绘制它们的数量随时间的变化趋势图，以展示可变性别比例对整个生态系统的影响。
4. **七鳃鳗与寄生虫相互作用图：**
- 绘制七鳃鳗与寄生虫之间相互作用的图表，以展示不同性别比例下它们之间的动态关系。
5. **生态系统稳定性分析图：**
- 分析模型中生态系统的稳定性，可以绘制相对稳定性指标随时间的变化图表，以了解可变性别比例对生态系统整体稳定性的影响。
6. **参数敏感性分析图：**
- 进行参数敏感性分析，绘制不同模型参数下性别比例、寄生虫数量等指标的变化图表，以确定模型中哪些参数对结果影响较大。
7. **可变性别比例优势和劣势分析图：**
- 绘制不同环境条件下，可变性别比例对七鳃鳗种群的优势和劣势的分析图表，以了解在不同情境下可变性别比例的影响。
8. **性别比例与其他生物种群相互作用图：**
- 如果模型中包含其他生物，可以绘制它们与七鳃鳗性别比例之间的相互作用图表，以揭示可变性别比例对生态系统中其他生物的影响。
同店推荐
包含文件 · 6 · 341.6 KB
- 01_美赛A题鹿鹿学长思路+代码.docx 30.3 KB
- 02_问题二.py 2.4 KB
- 03_问题三.py 2.6 KB
- 04_问题四.py 3.8 KB
- 05_问题一.py 2.3 KB
- 06_美赛A题鹿鹿学长思路+代码.pdf 300.2 KB
支付金额
¥17.20