逻辑斯谛方程,又称逻辑斯蒂增长模型,是一种描述在有限资源环境下种群规模增长的S形曲线数学模型。该模型由比利时数学家皮埃尔·弗朗索瓦·韦尔霍斯(Pierre François Verhulst)特在19世纪提出,如今已成为生态学、社会学、医学和机器学习等多个领域中的核心模型之一。
历史
逻辑斯谛方程的历史发展跨越了两个世纪,融合了人口学、生物学和数学的交叉研究:
- 理论奠基:逻辑斯蒂函数由比利时数学家皮埃尔·弗朗索瓦·韦尔霍斯特在1838年至1847年间发表的三篇论文中系统性地提出并发展。他在阿道夫·凯特勒的指导下,将其作为调整指数增长模型后的人口增长模型。韦尔霍斯特在1838年的论文中首次提出了逻辑斯蒂方程的微分形式,其中包含种群平衡点对应环境承载能力[1]。1845年,他正式将方程的解曲线命名为“逻辑斯蒂曲线”[2]。1847年,他的第三篇论文调整了对比利时人口增长模型中的修正项[3]。
- 重新发现与推广:韦尔霍斯特的工作在当时并未引起广泛关注。直到1920年,美国生物学家雷蒙德·珀尔和洛厄尔·里德在研究果蝇种群增长和人口问题时独立重新发现了这一方程,并提出了现在广为流传的形式。他们的工作使得该方程在生物学界广为人知,因此该方程有时也被称为韦尔霍斯特-珀尔阻碍方程。
- 混沌理论发展:随着20世纪系统科学和混沌理论的发展,该方程的离散形式——逻辑斯谛映射被深入研究。罗伯特·梅在1976年发表于《自然》杂志的论文《表现非常复杂的动力学的简单数学模型》中,系统性地阐述了这一简单非线性方程如何产生复杂的混沌行为,从而使该模型成为混沌研究的经典范例。
简介
问题背景
在自然界中,种群的增长并非总是无限制的。马尔萨斯模型描述了在资源无限情况下的指数增长,但这在现实中几乎不存在。大多数增长过程会受到环境承载力(如食物、空间、天敌等)的限制。韦尔霍斯特意识到人口增长不但和现有人口相关,还和可用资源有关,即存在一个人口的承载量。
核心思想
逻辑斯谛方程的核心思想是:种群的增长率不仅与当前的种群规模成正比,还与剩余的可利用资源空间成正比。随着种群规模接近环境承载力,其增长率将逐渐下降至零,从而实现种群的自我调节。这种增长模式表现为典型的S形曲线,包含开始期、加速期、转折期、减速期和饱和期五个阶段。
主要内容
微分形式与推导
逻辑斯谛方程的标准微分形式为:
-
[math]\displaystyle{ \begin{aligned} \frac{dN}{dt} &= rN(1-\frac{N}{K}) \end{aligned} }[/math]
(1)
其中:
- N(t): 在时间t时的种群规模;
- r: 内禀增长率,表示在理想条件下种群的最大增长潜力;
- K: 环境承载力,即环境所能支持的最大种群规模;
- [math]\displaystyle{ \frac{dN}{dt} }[/math]: 种群在时刻t的瞬时增长率。
该方程也可以从原始形式推导而来。韦尔霍斯最初提出的形式是[math]\displaystyle{ \frac{dN}{dt} = rN - \alpha N^2 }[/math] ,其中[math]\displaystyle{ \alpha }[/math]表示密度依赖的拥挤效应(种内竞争),通过变量替换[math]\displaystyle{ K = \frac{r}{\alpha} }[/math]可得标准形式。
[[文件:Schematic-diagram-of-a-simple-logistic-S-curve-defined-by-three-parameters-1 W640.jpg|缩略图|图1. 简单逻辑S曲线的示意图。增长率函数在同一图上以钟形曲线表示。引自[4]
解析解与S形曲线
逻辑斯谛方程是一个可分离变量的微分方程,其解析解为:
-
[math]\displaystyle{ \begin{aligned} N(t) = \frac{K}{1+(\frac{K-N_0}{N_0})e^{-rt}} \end{aligned} }[/math]
(2)
或者等价形式:[math]\displaystyle{ N(t) = \frac{K}{1 + C K e^{-rt}} }[/math]
其中[math]\displaystyle{ C = \frac{1}{N_0} - \frac{1}{K} }[/math]是由初始条件决定的常数。
该解描述了一条S形曲线(又称sigmoid曲线),其增长过程可分为三个典型阶段:
- 初始阶段(N ≪ K):增长近似指数增长,曲线呈凹形上升;
- 过渡阶段(N ≈ K/2):增长速度达到最大值,曲线出现拐点;
- 饱和阶段(N → K):增长逐渐减缓并趋近于零,种群规模稳定在承载能力K附近。
方程的性质
不动点与稳定性
逻辑斯谛方程有两个不动点(平衡点):
- [math]\displaystyle{ N = 0 }[/math]:不稳定平衡点,表示种群灭绝;
- [math]\displaystyle{ N = K }[/math]:稳定平衡点,表示种群达到环境承载能。
当 [math]\displaystyle{ N \lt K }[/math]时,逻辑斯蒂系数为正值,种群数量上升;当[math]\displaystyle{ N \gt K }[/math] 时,逻辑斯蒂系数为负值,种群数量下降;当[math]\displaystyle{ N = K }[/math] 时,逻辑斯蒂系数等于零,种群数量保持不变。
拐点
逻辑斯谛曲线的拐点出现在 [math]\displaystyle{ N = \frac{K}{2} }[/math]处,此时种群增长率达到最大值[math]\displaystyle{ \frac{rK}{4} }[/math] 。这一性质在实际应用中具有重要意义,例如在资源管理和流行病控制中,拐点可以帮助确定干预措施的最佳时机。
离散形式的逻辑斯谛方程
逻辑斯谛方程的离散形式称为逻辑斯谛映射。
推广
逻辑斯谛方程有多种推广形式,包括:
- 广义逻辑斯蒂曲线:引入更多参数以增强拟合灵活性。
- 冈珀茨函数:常用于描述生长速度不对称的生长过程。
- 超博拉斯特函数:用于更复杂的生长模式建模。
在统计学中,逻辑斯蒂函数的多元推广是softmax函数,用于处理多类别分类问题。
应用
生态学与人口学
逻辑斯谛方程最初且最直接的应用领域是生态学与人口学:
- 种群动态:描述有限资源下生物种群的增长规律,如果蝇、细菌等种群的生长。
- 人口预测:用于模拟和预测人类人口增长,如珀尔和里德对美国人口的研究。
- 资源管理:指导渔业、林业等可再生资源的可持续利用,避免过度开发。
医学:流行病传播模型
在医学领域,逻辑斯谛方程被广泛应用于流行病学建模:
- 疾病传播:在经典的SIR模型中,逻辑斯谛方程可以近似描述在疾病爆发初期,易感人群数量庞大时,感染者人数的增长模式。
- 疫情预测:在COVID-19疫情期间,广义逻辑斯蒂函数(理查兹增长曲线)被用于拟合感染轨迹,预测疫情发展趋势。
- 肿瘤生长:用于描述肿瘤在营养限制条件下的生长动力学。
化学与物理学
逻辑斯谛方程在化学和物理学中也有重要应用:
- 自催化反应:描述自催化化学反应中反应物和产物浓度的变化规律。
- 费米-狄拉克分布:在统计物理中,逻辑斯蒂函数决定了费米子在热平衡系统中各能态的占据概率分布。
- 相变模型:用于描述材料科学中的某些相变过程。
经济学与技术创新扩散
在经济学和社会科学中,逻辑斯谛方程用于描述各种扩散过程:
- 创新扩散**:巴斯扩散模型的核心就是逻辑斯谛方程,用于描述新产品或新技术在人群中采纳的过程:从早期采用者到晚期大众,最后市场饱和。
- 经济发展**:用于分析经济增长的长期周期和基础设施建设的扩散过程。
- 市场预测**:预测新产品市场占有率的增长和饱和趋势。
机器学习:逻辑斯谛回归
在机器学习领域,逻辑斯谛函数是逻辑斯谛回归的基础:
- 二分类问题:通过逻辑斯蒂函数将线性组合的输出映射到(0,1)区间,解释为概率。
- 多分类扩展:通过softmax函数推广到多类别分类问题。
- 深度学习:作为神经网络的激活函数,引入非线性特性。
逻辑斯谛回归源于1944年约瑟夫·伯克森基于逻辑斯蒂函数提出的logit函数,后来由大卫·考克斯在1958年正式命名为逻辑斯蒂回归。如今,它已成为分类问题中最广泛使用的算法之一。
参考
- ↑ Verhulst, Pierre-Francois (1838). Notice sur la loi que la population suit dans son accroissement. Correspondance mathématique et Physigue, publiée par A. Quetelet, 10, 113-120.
- ↑ Verhulst, Pierre-Francois (1845). Recherches mathématiques sur la loi d'accroissement dela population. Nouveaur Mémoires de l'Académie Royale des Sciences, des Lettres et desBeaux-Arts de Belgique, 18,1-38.
- ↑ Verhulst, Pierre-Francois (1847). Deuxième Mémoire sur la loi d'accroissement de la popula.tion. Nouveaur Mémoires de l'Académie Rouale des Sciences, des Lettres et des Beaur-Artsde Belgique, 20,1-32.
- ↑ Kucharavy, D., & De Guio, R. (2015). Application of logistic growth curve. Procedia engineering, 131, 280-290. https://doi.org/10.1016/j.proeng.2015.12.390