邏輯斯諦函數
邏輯斯諦函數(英語:logistic function)是一種常見的S型函數,其函數圖像稱為邏輯斯諦曲線(英語:logistic curve)。簡單的邏輯斯諦函數可用下式表示:
其中:
- x0為S形曲線中點的x值;
- L為曲線的最大值
- k為邏輯斯諦增長率或曲線的陡度。[1]
當x趨向於正無窮時,f(x)的值逼近L,而x趨向於負無窮時,f(x)的值逼近0。
邏輯斯諦函數應用領域廣泛,包括生物學(特別是生態學)、數理生物學、化學、人口學、經濟學、地球科學、數學心理學、概率、社會學、政治學、語言學、統計學和人工神經網絡等。例如,廣義邏輯斯諦曲線可以模仿一些情況人口增長(P)的S形曲線。起初階段大致是指數增長;然後隨着開始變得飽和,增長變慢;最後,達到成熟時增長停止。
歷史
[編輯]邏輯斯諦函數是皮埃爾·弗朗索瓦·韋呂勒於1838年至1847年間發表的三篇論文中提出的,他在阿道夫·凱特勒的指導下,通過調整指數增長模型,將其設計為人口增長模型。[2]韋呂勒在1830年代中期設計了該函數,並在1838年發表了一個簡短的說明,[1]然後在1844年進一步分析並命名了這個函數(發表於1845年)[3]第三篇論文調整了比利時人口增長模型中的修正項。[4]
增長的初始階段近似於指數增長(幾何級數);然後,隨着增長逐漸飽和,曲線放緩至接近線性,在成熟階段,增長停止。原本選用「邏輯斯諦」(法語:logistique,英語:logistic)一詞時,韋呂勒沒有解釋其原由,但這可能是為了區別於對數曲線,[5][a]並與算術和幾何進行對比。在提出該增長模型前,他討論了算術增長和幾何增長(他稱之為「對數曲線」,其現代通稱是指數曲線),因此「邏輯斯諦增長」可能是通過類比命名的,「邏輯斯諦」來自古希臘語λογῐστῐκός(logistikós),是指古希臘數學的一個分支。[b]「邏輯斯諦函數」中的「邏輯」與邏輯學(logic)和軍隊後勤/物流(logistics,自法語logis)均沒有關係。
數學特性
[編輯]標準邏輯斯諦函數的參數設定為, , ,即
實際上,由於指數函數的特性,函數的取值很快會逼近極限,x在很小的實數範圍內(例如[−6, +6])的取值就足以計算標準邏輯斯諦函數的極限。
標準邏輯斯諦函數具有如下對稱性:
因此,是奇函數。
標準邏輯斯諦函數可視為雙曲正切函數的偏移和縮放:
或
推導過程如下:
導數
[編輯]邏輯斯諦分布的均值為x0,方差為π2/3k2。
積分
[編輯]標準邏輯斯諦函數的不定積分可用換元積分法求得,令,,去掉積分常數,得到其不定積分:
在人工神經網絡中,它稱作線性整流函數,(縮放後)可視為平滑近似的斜坡函數,類似於邏輯斯諦函數(縮放後)是平滑近似的單位階躍函數。
邏輯斯諦微分方程
[編輯]標準邏輯斯諦函數是簡單的一階非線性常微分方程的解:
邊界條件為。該方程是邏輯斯諦映射的連續版本。注意倒數邏輯斯諦函數是簡單的一階線性常微分方程的解。[6]
邏輯斯諦差分方程
[編輯]是混沌理論的一個模型。[7][8]這個函數對初始值和參數的變化很敏感,往往微小的變化會引起混沌。如圖所示,當x1=0.3,參數k從0.1變到4時,系統變化很大。
- 當k由0.1變到1時,曲線很快趨向於0
- 當k繼續增加,曲線由0.3上升到 一個穩定值
- k繼續增加,曲線出現擺動,有2個穩定值。
- k繼續增加, 曲線相繼出現4個、8個、16個、32個....穩定值
- k增加到一個臨界值,系統進入混沌狀態。
- k再增加,系統突然垮塌。
變化
[編輯]應用
[編輯]生態學:種群增長模型
[編輯]邏輯斯諦方程的一個典型應用是種群(或人口)增長的通用模型(另見種群動態),最初由皮埃爾·弗朗索瓦·韋呂勒在1838年提出,其中繁殖率與現狀種群數量和可用資源量成正比,其他一切都條件均等。韋呂勒方程是他在閱讀馬爾薩斯的論文《An Essay on the Principle of Population 》後發表的,該論文描述了簡單(無約束條件)指數增長的馬爾薩斯模型。韋呂勒推導出他的邏輯斯諦方程來描述生物種群的自限性增長。該方程於1911年被A. G. McKendrick用於描述肉湯中細菌的生長,他使用非線性參數估計的方法進行了實驗測試。[9]在約翰斯·霍普金斯大學的Raymond Pearl(1879–1940)和Lowell Reed(1888–1966)於1920年使用該方程後,這一方程有時也稱為Verhulst-Pearl方程。[10]另一位科學家阿弗雷德·洛特卡在1925年再次推導出該方程,稱其為種群增長律(law of population growth)。
令P為種群(人口)規模(生態學經常用N代替),t代表時間,該模型用以下微分方程表示:
方程中,早期的幾乎無阻力的增長率來自+rP。增長率r代表種群(人口)數量P在一個單位時間內的增長比例。後來,隨着人口的增長,第二項-rP2/K變得幾乎和第一項一樣大,種群P內的個體之間開始爭奪某些關鍵資源(例如食物或生存空間)而相互干擾。這種對抗效應稱為「瓶頸」,由參數K代表。競爭會降低總合增長率,直到P停止增長(種群/人口成熟)。方程的解(P0為初始種群/人口數量)為
其中:
可以說,K是P的極限值,即經過無限長時間後(或在有限時間內近似),種群(人口)規模所能達到的最大值。須注意,只要初始值,無論取值多少,種群數量都會漸近環境承載力的值,包括的情況下。
生態學中有時稱一個物種是r策略或K策略的,這是指它們在自然選擇過程形成的生物生命週期策略。選取變量的量綱,使n代表以環境承載力單位計的種群數量,代表以的單位計量的時間,得出無量綱微分方程:
時變承載力
[編輯]由於環境條件會影響環境承載力,因此它可能是隨時間變化的,,得出以下數學模型:
其中一種特別重要的情況是承載力隨時期以T為周期變化的情況:
可見,只要初始值,無論具體取值為多少,會逼近一個周期為的周期解。
的典型取值為1年,在此情況下,可表示天氣條件的周期性變化。
另一個有趣的一般化情形是考慮承載能力K(t)作為關於較早時間的種群數量的函數,以表示種群改變其所處環境的延遲。這就構成了一個邏輯斯諦時滯方程,[11]它具有非常豐富的行為,在某些參數範圍內呈現雙穩定,以及單調衰減至零、平滑指數增長、間斷無限增長(即多個S形)、間斷增長或交替到平穩水平、振盪接近穩定水平、持續振盪、有限時間奇異點以及有限時間死亡。
統計學和機器學習
[編輯]邏輯斯諦函數在統計學中有多種應用。例如,它們是邏輯斯諦分布的累積分布函數,它們可用於模擬國際象棋棋手在埃洛等級分系統下擊敗對手的概率。以下是一些更具體的案例。
邏輯迴歸
[編輯]邏輯斯諦迴歸使用邏輯斯諦函數來模擬一個事件的概率p如何可能會受到一個或多個解釋變量的影響:一個案例模型如下
其中x為解釋變量,a和b為欲擬合的模型參數,f為標準邏輯斯諦函數。
邏輯斯諦迴歸和其他對數線性模型也常用於機器學習。將邏輯斯諦函數推廣至多元輸入情景即為Softmax激活函數,用於多元羅吉斯回歸。
神經網絡
[編輯]醫學:腫瘤生長模型
[編輯]在醫學上,邏輯斯諦微分方程可用於腫瘤生長的建模。這一用法可視為上述的生態學/人口學模型的延伸。以表示腫瘤在時間的大小,其變化動態遵循
屬於以下類型:
其中為腫瘤增殖率。
如果採用化療產生對數殺傷效果,則等式修改為
其中為治療引起的腫瘤死亡率。在理想化的極長的治療下,可模型化為周期為的周期函數或(在持續的輸液治療下)常數函數,有
即,如果平均治療引起的腫瘤死亡率大於基線增殖率,則疾病能被根除。當然,這是一個過於簡化的生長和治療模型(例如沒有考慮克隆抗性現象)。
醫學:傳染病模型
[編輯]在人群中未被免疫的新型傳染性病原體,通常會在早期呈指數級傳播,有大量易感個體尚未被感染。例如2020年初,導致2019冠狀病毒病的SARS-CoV-2病毒在多國的感染過程中呈現出指數級增長。[12]此後,易感宿主減少(持續感染直到超過群體免疫閾值)或通過社交距離措施減少潛在宿主的被傳染概率等因素,可能使呈指數增長的傳染曲線首先線性化,然後趨緩,達到最大值。[13]
邏輯斯諦函數或相關的函數(例如龔珀茲函數)通常以描述性或現象學方的式使用,因為它們非常符合早期的指數上升,也符合隨着人群形成群體免疫而最終趨於平穩的趨勢。它與流行病的實際模型不同,後者試圖根據大流行的動態(例如接觸率、潛伏期、社交距離等)來描述感染狀態。不過,一些簡單的模型有邏輯斯諦解。[14][15][16]
早期COVID-19病例數建模
[編輯]廣義邏輯斯諦函數(又稱Richards增長曲線)已應用於對COVID-19爆發的早期階段建模。[17]研究者將廣義邏輯斯諦函數擬合到累計感染病例數(稱為傳染軌跡)。文獻中對廣義邏輯斯諦函數有不同的參數化。一種常用的形式是:
其中取實數,為正實數。曲線的靈活性由賦予:(i)若,則曲線衰減為邏輯斯諦函數,(ii)若收斂至0,則曲線收斂至龔珀茲函數。在傳染病模型中,, 和分別代表傳染病最終的規模、感染率和滯後期。見右側的範例的傳染軌跡,其中設定為。
在流行病學建模中,使用類似廣義邏輯斯諦函數的增長函數的好處之一是它相對容易應用於多級模型框架,其中來自不同地理區域的信息可以匯總在一起。
化學:反應模型
[編輯]自催化反應中,反應物和產物的濃度遵循邏輯斯諦函數。例如燃料電池陰極中不含鉑族金屬(PGM-free)的氧還原反應催化劑的劣化遵循邏輯斯諦衰減函數,[18]表明這是一種自催化分解機制。
物理:費米-狄拉克分布
[編輯]費米子在熱平衡系統的能量狀態上的統計分布遵循邏輯斯諦函數。特別地,根據費米-狄拉克統計,它是每個可能的能級被一個費米子占據的概率分布。
語言學:語言變化
[編輯]語言學中,邏輯斯諦函數可用於對語言變化進行建模:[19]一種最初處於邊緣地位的新詞隨着時間的推移開始傳播得更快,然後傳播速度隨着其普及而減慢。
經濟學和社會學:創新擴散
[編輯]邏輯斯諦函數可用於描繪一項發明創新在其生命周期內擴散的過程.
參見
[編輯]注釋
[編輯]參考文獻
[編輯]- ^ 1.0 1.1 Verhulst, Pierre-François. Notice sur la loi que la population poursuit dans son accroissement (PDF). Correspondance Mathématique et Physique. 1838, 10: 113–121 [3 December 2014].
- ^ Cramer, J. S. The origins of logistic regression (PDF) (技術報告) 119. Tinbergen Institute: 167–178. 2002 [2022-05-05]. doi:10.2139/ssrn.360300. (原始內容存檔 (PDF)於2022-07-21).
|issue=
被忽略 (幫助) - ^ Verhulst, Pierre-François. Recherches mathématiques sur la loi d'accroissement de la population [Mathematical Researches into the Law of Population Growth Increase]. Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Bruxelles. 1845, 18: 8 [18 February 2013].
Nous donnerons le nom de logistique à la courbe [We will give the name logistic to the curve]
- ^ Verhulst, Pierre-François. Deuxième mémoire sur la loi d'accroissement de la population. Mémoires de l'Académie Royale des Sciences, des Lettres et des Beaux-Arts de Belgique. 1847, 20: 1–32 [18 February 2013]. (原始內容存檔於2016-05-18).
- ^ Shulman, Bonnie. Math-alive! using original sources to teach mathematics in social context. PRIMUS. 1998, 8 (March): 1–14. doi:10.1080/10511979808965879.
The diagram clinched it for me: there two curves labeled "Logistique" and "Logarithmique" are drawn on the same axes, and one can see that there is a region where they match almost exactly, and then diverge.
I concluded that Verhulst's intention in naming the curve was indeed to suggest this comparison, and that "logistic" was meant to convey the curve's "log-like" quality. - ^ Kocian, Alexander; Carmassi, Giulia; Cela, Fatjon; Incrocci, Luca; Milazzo, Paolo; Chessa, Stefano. Bayesian Sigmoid-Type Time Series Forecasting with Missing Data for Greenhouse Crops. Sensors. 7 June 2020, 20 (11): 3246. Bibcode:2020Senso..20.3246K. PMC 7309099 . PMID 32517314. doi:10.3390/s20113246 .
- ^ Garnett P. Williams Chaos Theory Tamed chapter 10
- ^ Edgar Peters Chaos and Order in the Capital Market p7
- ^ A. G. McKendricka; M. Kesava Paia1. XLV.—The Rate of Multiplication of Micro-organisms: A Mathematical Study. Proceedings of the Royal Society of Edinburgh. January 1912, 31: 649–653 [2022-05-05]. doi:10.1017/S0370164600025426. (原始內容存檔於2021-11-20).
- ^ Raymond Pearl & Lowell Reed. On the Rate of Growth of the Population of the United States (PDF) 6 (6). June 1920: 275 [2022-05-05]. (原始內容存檔 (PDF)於2021-07-10).
|journal=
被忽略 (幫助) - ^ Yukalov, V. I.; Yukalova, E. P.; Sornette, D. Punctuated evolution due to delayed carrying capacity. Physica D: Nonlinear Phenomena. 2009, 238 (17): 1752–1767. Bibcode:2009PhyD..238.1752Y. S2CID 14456352. arXiv:0901.4714 . doi:10.1016/j.physd.2009.05.011.
- ^ Worldometer: COVID-19 CORONAVIRUS PANDEMIC. [2022-05-05]. (原始內容存檔於2020-03-22).
- ^ Villalobos-Arias, Mario. Using generalized logistics regression to forecast population infected by Covid-19. 2020. arXiv:2004.02406 [q-bio.PE].
- ^ Postnikov, Eugene B. Estimation of COVID-19 dynamics "on a back-of-envelope": Does the simplest SIR model provide quantitative parameters and predictions?. Chaos, Solitons & Fractals. June 2020, 135: 109841. Bibcode:2020CSF...13509841P. PMC 7252058 . PMID 32501369. doi:10.1016/j.chaos.2020.109841.
- ^ Saito, Takesi. A Logistic Curve in the SIR Model and Its Application to Deaths by COVID-19 in Japan. medRxiv. June 2020 [July 20, 2020]. S2CID 220068969. doi:10.1101/2020.06.25.20139865. (原始內容存檔於2022-04-30).
- ^ Reiser, Paul A. Modified SIR Model Yielding a Logistic Solution. 2020. arXiv:2006.01550 [q-bio.PE].
- ^ Lee, Se Yoon; Lei, Bowen; Mallick, Bani. Estimation of COVID-19 spread curves integrating global data and borrowing information. PLOS ONE. 2020, 15 (7): e0236860. Bibcode:2020PLoSO..1536860L. PMC 7390340 . PMID 32726361. arXiv:2005.00662 . doi:10.1371/journal.pone.0236860 .
- ^ Yin, Xi; Zelenay, Piotr. Kinetic Models for the Degradation Mechanisms of PGM-Free ORR Catalysts. ECS Transactions. 13 July 2018, 85 (13): 1239–1250 [2022-05-05]. OSTI 1471365. S2CID 103125742. doi:10.1149/08513.1239ecst. (原始內容存檔於2022-05-05).
- ^ Bod, Hay, Jennedy (eds.) 2003, pp. 147–156
- Jannedy, Stefanie; Bod, Rens; Hay, Jennifer. Probabilistic Linguistics. Cambridge, Massachusetts: MIT Press. 2003. ISBN 0-262-52338-8.
- Gershenfeld, Neil A. The Nature of Mathematical Modeling. Cambridge, UK: Cambridge University Press. 1999. ISBN 978-0-521-57095-4.
- Kingsland, Sharon E. Modeling nature: episodes in the history of population ecology. Chicago: University of Chicago Press. 1995. ISBN 0-226-43728-0.
- 埃里克·韋斯坦因. Logistic Equation. MathWorld.
外部連結
[編輯]- L.J. Linacre, Why logistic ogive and not autocatalytic curve? (頁面存檔備份,存於網際網路檔案館), accessed 2009-09-12.
- https://web.archive.org/web/20060914155939/http://luna.cas.usf.edu/~mbrannic/files/regression/Logistic.html
- Modeling Market Adoption in Excel with a simplified s-curve (頁面存檔備份,存於網際網路檔案館)
- 埃里克·韋斯坦因. Sigmoid Function. MathWorld.
- Online experiments with JSXGraph (頁面存檔備份,存於網際網路檔案館)