Lotka-Volterra捕食者-猎物模型R代码片段逻辑解读咨询

1. 两层for循环逻辑解释 #

你提问里贴的循环代码存在一个笔误，n2的seq参数多了一个多余右括号，实际生效的正确逻辑如下：

# 外层循环：遍历物种1的初始种群数量
for (n1 in seq((0.1*k[1]),(1.2*k[1]),100)){
  # 内层循环：遍历物种2的初始种群数量
  for (n2 in seq(.1*k[2],1.2*k[2],.1*k[2])){

• k[1]、k[2]是预先定义的两个物种的环境容纳量，代码里初始赋值为1000
• 物种1的初始数量取值范围是容纳量的10%到120%，步长固定为100，原参数下取值为100、200...1200
• 物种2的初始数量取值范围同样是容纳量的10%到120%，步长为容纳量的10%，原参数下取值为100、200...1200
• 两层循环的作用是遍历不同初始种群组合，逐一运行LV模型模拟，最终生成相图的轨线簇。

2. if判断语句逻辑解释 #

if(n1==100&n2==100){plot(y[,2],y[,3],type="l",xlim=c(0,1000),...)}

• 100是原参数下两个物种初始数量为容纳量10%的取值，这个判断的作用是第一次模拟时初始化画布，设置坐标轴范围、标签、基础绘图样式
• 后续其他初始值的模拟结果不需要重新创建画布，直接用后面的lines()函数把轨线叠加到同一张图上，避免重复绘图
• 后续的4个if判断是为了只绘制边界初始值的轨线（两个物种初始值取上下限的情况），避免图中轨线过多杂乱，保持可视化效果清晰。

3. 适配其他尺度数据的修改方案 #

原代码大量用了基于k=1000的硬编码数值，修改时按以下步骤调整即可：

• 第一步修改核心参数：直接修改开头的r（种群增长率）、k（环境容纳量）、alpha（种间竞争系数）三个矩阵即可，不需要动后续逻辑
• 第二步修正循环步长：把外层n1的固定步长100改成和n2一致的相对步长，改为seq(0.1*k[1], 1.2*k[1], 0.1*k[1])，避免k修改后步长不合理
• 第三步修正画布坐标范围：把plot里写死的xlim=c(0,1000), ylim=c(0,1000)改成xlim=c(0, 1.2*k[1]), ylim=c(0, 1.2*k[2])，自动适配新的容纳量尺度
• 第四步修正判断条件的硬编码：把所有if判断里的固定数值替换为相对值：
  • n1==100 改为 n1 == 0.1*k[1]
  • n2==100 改为 n2 == 0.1*k[2]
  • n1==1200 改为 n1 == 1.2*k[1]
  • n2==1200 改为 n2 == 1.2*k[2]
• 按需调整t_stop（模拟时长）即可，种群收敛慢的话可以调大该数值。

内容的提问来源于stack exchange，提问作者springle