邏輯斯諦函數

數學函數

邏輯斯諦函數(英語:logistic function)是一種常見的S型函數,其函數圖像稱為邏輯斯諦曲線(英語:logistic curve)。簡單的邏輯斯諦函數可用下式表示:

標準邏輯斯諦函數,其中

其中:

x0為S形曲線中點的x值;
L為曲線的最大值
k為邏輯斯諦增長率或曲線的陡度。[1]

x趨向於正無窮時,f(x)的值逼近L,而x趨向於負無窮時,f(x)的值逼近0。

邏輯斯諦函數應用領域廣泛,包括生物學(特別是生態學)、數理生物學化學人口學經濟學地球科學數學心理學機率社會學政治學語言學統計學人工神經網絡等。例如,廣義邏輯斯諦曲線英語generalized logistic curve可以模仿一些情況人口增長(P)的S形曲線。起初階段大致是指數增長;然後隨著開始變得飽和,增長變慢;最後,達到成熟時增長停止。

歷史

編輯
 
邏輯斯諦曲線的原始圖像,與指數曲線對比

邏輯斯諦函數是皮埃爾·弗朗索瓦·韋呂勒英語Pierre François Verhulst於1838年至1847年間發表的三篇論文中提出的,他在阿道夫·凱特勒的指導下,通過調整指數增長模型,將其設計為人口增長模型。[2]韋呂勒在1830年代中期設計了該函數,並在1838年發表了一個簡短的說明,[1]然後在1844年進一步分析並命名了這個函數(發表於1845年)[3]第三篇論文調整了比利時人口增長模型中的修正項。[4]

增長的初始階段近似於指數增長幾何級數);然後,隨著增長逐漸飽和,曲線放緩至接近線性,在成熟階段,增長停止。原本選用「邏輯斯諦」(法語:logistique,英語:logistic)一詞時,韋呂勒沒有解釋其原由,但這可能是為了區別於對數曲線英語logarithmic curve[5][a]並與算術和幾何進行對比。在提出該增長模型前,他討論了算術增長和幾何增長(他稱之為「對數曲線」,其現代通稱是指數曲線),因此「邏輯斯諦增長」可能是通過類比命名的,「邏輯斯諦」來自古希臘語λογῐστῐκόςlogistikós),是指古希臘數學的一個分支。[b]「邏輯斯諦函數」中的「邏輯」與邏輯學(logic)和軍隊後勤/物流(logistics,自法語logis)均沒有關係。

數學特性

編輯

標準邏輯斯諦函數的參數設定為 ,  ,  ,即

 

實際上,由於指數函數 的特性,函數的取值很快會逼近極限x在很小的實數範圍內(例如[−6, +6])的取值就足以計算標準邏輯斯諦函數的極限。

標準邏輯斯諦函數具有如下對稱性:

 

因此, 奇函數

標準邏輯斯諦函數可視為雙曲正切函數的偏移和縮放:

 

 

推導過程如下:

 

導數

編輯

標準邏輯斯諦函數的導數稱為邏輯斯諦分布英語logistic distribution密度,公式如下:

 
 

邏輯斯諦分布的均值為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再增加,系統突然垮塌。

變化

編輯
 
 

應用

編輯

生態學:族群增長模型

編輯
 
皮埃爾·弗朗索瓦·韋呂勒(1804–1849)

邏輯斯諦方程式的一個典型應用是族群(或人口)增長的通用模型(另見族群動態英語population dynamics),最初由皮埃爾·弗朗索瓦·韋呂勒英語Pierre François Verhulst在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代表時間,該模型用以下微分方程式表示:

 

其中常數r族群(人口)增長率K環境承載力

方程式中,早期的幾乎無阻力的增長率來自+rP。增長率r代表族群(人口)數量P在一個單位時間內的增長比例。後來,隨著人口的增長,第二項-rP2/K變得幾乎和第一項一樣大,族群P內的個體之間開始爭奪某些關鍵資源(例如食物或生存空間)而相互干擾。這種對抗效應稱為「瓶頸」,由參數K代表。競爭會降低總合增長率,直到P停止增長(族群/人口成熟)。方程式的解(P0為初始族群/人口數量)為

 

其中:

 

可以說,KP的極限值,即經過無限長時間後(或在有限時間內近似),族群(人口)規模所能達到的最大值。須注意,只要初始值 ,無論取值多少,族群數量都會漸近環境承載力的值,包括 的情況下。

生態學中有時稱一個物種r策略或K策略的,這是指它們在自然選擇過程形成的生物生命週期策略。選取變量的量綱,使n代表以環境承載力單位計的族群數量, 代表以 的單位計量的時間,得出無量綱微分方程式:

 

時變承載力

編輯

由於環境條件會影響環境承載力,因此它可能是隨時間變化的, ,得出以下數學模型:

 

其中一種特別重要的情況是承載力隨時期以T為週期變化的情況:

 

可見,只要初始值 ,無論具體取值為多少, 會逼近一個週期為 的週期解 

 的典型取值為1年,在此情況下, 可表示天氣條件的週期性變化。

另一個有趣的一般化情形是考慮承載能力K(t)作為關於較早時間的族群數量的函數,以表示族群改變其所處環境的延遲。這就構成了一個邏輯斯諦時滯方程式,[11]它具有非常豐富的行為,在某些參數範圍內呈現雙穩定,以及單調衰減至零、平滑指數增長、間斷無限增長(即多個S形)、間斷增長或交替到平穩水平、振盪接近穩定水平、持續振盪、有限時間奇異點以及有限時間死亡。

統計學和機器學習

編輯

邏輯斯諦函數在統計學中有多種應用。例如,它們是邏輯斯諦分布英語Logistic distribution累積分布函數,它們可用於模擬西洋棋棋手在埃洛等級分系統下擊敗對手的機率。以下是一些更具體的案例。

邏輯迴歸

編輯

邏輯斯諦迴歸使用邏輯斯諦函數來模擬一個事件的機率p如何可能會受到一個或多個解釋變量的影響:一個案例模型如下

 

其中x為解釋變量,ab為欲擬合的模型參數,f為標準邏輯斯諦函數。

邏輯斯諦迴歸和其他對數線性模型英語log-linear model也常用於機器學習。將邏輯斯諦函數推廣至多元輸入情景即為Softmax激活函數,用於多元邏輯迴歸英語multinomial logistic regression

神經網絡

編輯

醫學:腫瘤生長模型

編輯

在醫學上,邏輯斯諦微分方程式可用於腫瘤生長的建模。這一用法可視為上述的生態學/人口學模型的延伸。以 表示腫瘤在時間 的大小,其變化動態遵循

 

屬於以下類型:

 

其中 為腫瘤增殖率。

如果採用化療產生對數殺傷效果,則等式修改為

 

其中 為治療引起的腫瘤死亡率。在理想化的極長的治療下, 可模型化為週期為 的週期函數或(在持續的輸液治療下)常數函數,有

 

即,如果平均治療引起的腫瘤死亡率大於基線增殖率,則疾病能被根除。當然,這是一個過於簡化的生長和治療模型(例如沒有考慮克隆抗性現象)。

醫學:傳染病模型

編輯

在人群中未被免疫的新型傳染性病原體,通常會在早期呈指數級傳播,有大量易感個體尚未被感染。例如2020年初,導致2019冠狀病毒病SARS-CoV-2病毒在多國的感染過程中呈現出指數級增長。[12]此後,易感宿主減少(持續感染直到超過群體免疫閾值)或通過社交距離措施減少潛在宿主的被傳染機率等因素,可能使呈指數增長的傳染曲線首先線性化,然後趨緩,達到最大值。[13]

邏輯斯諦函數或相關的函數(例如龔珀茲函數英語Gompertz function)通常以描述性或現象學方的式使用,因為它們非常符合早期的指數上升,也符合隨著人群形成群體免疫而最終趨於平穩的趨勢。它與流行病的實際模型不同,後者試圖根據大流行的動態(例如接觸率、潛伏期、社交距離等)來描述感染狀態。不過,一些簡單的模型有邏輯斯諦解。[14][15][16]

早期COVID-19病例數建模

編輯
 
流行病模型中的廣義邏輯斯諦曲線英語Generalised logistic function(Richards增長曲線)

廣義邏輯斯諦函數英語Generalised logistic function(又稱Richards增長曲線)已應用於對COVID-19爆發的早期階段建模。[17]研究者將廣義邏輯斯諦函數擬合到累計感染病例數(稱為傳染軌跡)。文獻中對廣義邏輯斯諦函數有不同的參數化。一種常用的形式是:

 

其中 取實數, 為正實數。曲線 的靈活性由 賦予:(i)若 ,則曲線衰減為邏輯斯諦函數,(ii)若 收斂至0,則曲線收斂至龔珀茲函數。在傳染病模型中, ,   分別代表傳染病最終的規模、感染率和滯後期。見右側的範例的傳染軌跡,其中 設定為 

 
受COVID-19嚴重影響的40個國家的傳染曲線外推和截至5月14日的母體(人口)平均值

在流行病學建模中,使用類似廣義邏輯斯諦函數的增長函數的好處之一是它相對容易應用於多級模型框架,其中來自不同地理區域的資訊可以匯總在一起​​。

化學:反應模型

編輯

自催化反應中,反應物和產物的濃度遵循邏輯斯諦函數。例如燃料電池陰極中不含鉑族金屬(PGM-free)的氧還原反應催化劑的劣化遵循邏輯斯諦衰減函數,[18]表明這是一種自催化分解機制。

物理:費米-狄拉克分布

編輯

費米子在熱平衡系統的能量狀態上的統計分布遵循邏輯斯諦函數。特別地,根據費米-狄拉克統計,它是每個可能的能級被一個費米子占據的機率分布。

語言學:語言變化

編輯

語言學中,邏輯斯諦函數可用於對語言變化進行建模:[19]一種最初處於邊際地位的新詞隨著時間的推移開始傳播得更快,然後傳播速度隨著其普及而減慢。

經濟學和社會學:創新擴散

編輯

邏輯斯諦函數可用於描繪一項發明創新在其生命週期內擴散的過程.

參見

編輯

注釋

編輯
  1. ^ 韋呂勒首先提到算術級數和幾何級數,並將幾何增長曲線稱為「對數曲線 / logarithmic curve」(但須注意,現代術語稱之為「指數曲線 / exponential curve」)。然後他把他提出的曲線叫做邏輯斯諦(logistic),以區別於對數(logarithmic),並在他論文的圖中比較了「對數」曲線和邏輯斯諦曲線。
  2. ^ 在古希臘,λογῐστῐκός指實務的計算與會計,而ἀριθμητικήarithmētikḗ)指對於數的理論或哲學研究。然而在英語中,表示算術的詞反而是源自ἀριθμητική的arithmetic。

參考文獻

編輯
  1. ^ 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]. 
  2. ^ 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=被忽略 (幫助)
  3. ^ 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] 
  4. ^ 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). 
  5. ^ 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.
     
  6. ^ 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 . 
  7. ^ Garnett P. Williams Chaos Theory Tamed chapter 10
  8. ^ Edgar Peters Chaos and Order in the Capital Market p7
  9. ^ 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). 
  10. ^ 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=被忽略 (幫助)
  11. ^ 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. 
  12. ^ Worldometer: COVID-19 CORONAVIRUS PANDEMIC. [2022-05-05]. (原始內容存檔於2020-03-22). 
  13. ^ Villalobos-Arias, Mario. Using generalized logistics regression to forecast population infected by Covid-19. 2020. arXiv:2004.02406  [q-bio.PE]. 
  14. ^ 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. 
  15. ^ 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). 
  16. ^ Reiser, Paul A. Modified SIR Model Yielding a Logistic Solution. 2020. arXiv:2006.01550  [q-bio.PE]. 
  17. ^ 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 . 
  18. ^ 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). 
  19. ^ Bod, Hay, Jennedy (eds.) 2003, pp. 147–156

外部連結

編輯