模型参数
解析解: $N(t) = \dfrac{K}{1 + \left(\dfrac{K-N_0}{N_0}\right)e^{-rt}}$
当N≪K时 (1-N/K)≈1 呈指数增长,N→K时dN/dt→0(增长停止)。增长速度dN/dt在 N=K/2 时达最大值(变曲点)。
Lotka-Volterra(捕食者-被食者)方程
捕食者: $\dfrac{dP}{dt} = \delta NP - \gamma P$
解为周期性振荡(位相错开的循环)
通过滑块操控增长率r、环保容纳量K、初始个体数N₀,实时计算个体数的S型增长曲线。对比指数增长或体验捕食者-被食者系(Lotka-Volterra)。
当N≪K时 (1-N/K)≈1 呈指数增长,N→K时dN/dt→0(增长停止)。增长速度dN/dt在 N=K/2 时达最大值(变曲点)。
逻辑斯谛增长模型描述个体群增加受环保容纳量 \( K \) 限制的现象。基本方程是微分方程 \( \frac{dN}{dt} = rN \left(1 - \frac{N}{K}\right) \),其中 \( r \) 是内禀自然增长率,\( N \) 是当前个体数。初始个体数 \( N_0 \) 越小,初期增长越接近指数函数,但随着个体数接近 \( K \),增速减缓,最终在 \( N = K \) 处达到平衡。这种S型(Sigmoid)曲线反映了资源有限性。相比之下,指数增长模型 \( \frac{dN}{dt} = rN \) 无限制,显示不可持续的爆炸性增长。此外,捕食者-被食者系统(Lotka-Volterra模型)用 \( \frac{dx}{dt} = \alpha x - \beta xy \)、\( \frac{dy}{dt} = \delta xy - \gamma y \) 表达被食者 \( x \) 和捕食者 \( y \) 的相互作用,可观察到两者个体数的周期振荡。本模拟器基于这些方程,通过滑块调整参数并实时更新图表,让用户直观理解生态系统动态。
$\frac{dN}{dt} = rN\left(1 - \frac{N}{K}\right)$产业实际应用案例
在水产养殖业中,本模拟器应用于三文鱼、黄鰤等养殖规划。例如,日本丸红生鲜公司的养殖事业中,将环保容纳量K设定为养殖池容量或氧供应量,通过调整初始个体数N₀和增长率r,计算出避免过度拥挤导致增长停滞的最优放养密度。这使得从鱼苗到出货的增长曲线得以预测,提高了饲料投喂量和收获时期的计划效率。
研究与教育应用
在大学生态学和资源管理学讲座中,这是理解种群动态的必要工具。在东京大学的演习中,学生通过滑块改变r和K,视觉对比指数增长和逻辑斯谛增长的差异。在Lotka-Volterra模式下重现捕食者(狼)和被食者(鹿)的振荡,讨论生态系统的稳定性和灭绝风险,作为教学素材被广泛使用。
与CAE分析的联动及实务定位
本模拟器用于CAE前期处理阶段的生物学参数敏感性分析。例如,与水槽流体分析(CFD)联动,从氧浓度分布推估局部环保容纳量,作为K值输入,计算增长曲线以在设计阶段优化养殖效率。在实务中,通过实验数据拟合进行参数辨识,或作为设备投资决策支持工具。
容易误认为"个体数超过K会自动减少",但实际上逻辑斯谛增长模型中,个体数超过K时增长率为负,理论上开始减少。然而现实生态系统中,K超额导致的资源枯竭或环境恶化会引发复杂变化,简单S曲线无法完全再现。
容易误认为"增长率r越高,到达K越快",但高r会导致个体数大幅超过K后急剧下降,随后产生振荡。特别是捕食者-被食者系统模拟中,r的微小变化会剧烈改变周期和稳定性,参数设置需谨慎。
容易误认为"初始个体数N₀为零时永远为零",但逻辑斯谛模型中N₀=0是平衡点,在无外部移入或随机事件情况下不会增加。在进行人口预测时,初始条件的微小差异可能对长期结果产生巨大影响,需充分理解这一点。
日本梅花鹿5年预测:初始个体数N0=80只,内禀自然增长率r=0.28/年,环保容纳量K=1200只。设定这些参数后,第1年约103只、第3年约286只、第5年约621只,呈阶段性增长并描绘S型曲线。若用同一参数应用Lotka-Volterra方程(捕食-被食系)模拟天敌狼导入,则观察约2年周期的个体数振荡,最终在平衡点附近稳定。