邏輯斯諦函數 (英語:logistic function )是一種常見的S型函數 ,其函數圖像稱為邏輯斯諦曲線 (英語:logistic curve )。簡單的邏輯斯諦函數可用下式表示:
f
(
x
)
=
L
1
+
e
−
k
(
x
−
x
0
)
{\displaystyle f(x)={\frac {L}{1+e^{-k(x-x_{0})}}}}
標準邏輯斯諦函數,其中
L
=
1
,
k
=
1
,
x
0
=
0
{\displaystyle L=1,k=1,x_{0}=0}
其中:
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 )均沒有關係。
標準邏輯斯諦函數的參數設定為
k
=
1
{\displaystyle k=1}
,
x
0
=
0
{\displaystyle x_{0}=0}
,
L
=
1
{\displaystyle L=1}
,即
f
(
x
)
=
1
1
+
e
−
x
=
e
x
e
x
+
1
=
1
2
+
1
2
tanh
(
x
2
)
{\displaystyle f(x)={\frac {1}{1+e^{-x}}}={\frac {e^{x}}{e^{x}+1}}={\frac {1}{2}}+{\frac {1}{2}}\tanh \left({\frac {x}{2}}\right)}
實際上,由於指數函數
e
−
x
{\displaystyle e^{-x}}
的特性,函數的取值很快會逼近極限 ,x 在很小的實數範圍內(例如[−6, +6])的取值就足以計算標準邏輯斯諦函數的極限。
標準邏輯斯諦函數具有如下對稱性:
1
−
f
(
x
)
=
f
(
−
x
)
{\displaystyle 1-f(x)=f(-x)}
因此,
x
↦
f
(
x
)
−
1
/
2
{\displaystyle x\mapsto f(x)-1/2}
是奇函數 。
標準邏輯斯諦函數可視為雙曲正切函數的偏移和縮放:
f
(
x
)
=
1
2
+
1
2
tanh
(
x
2
)
{\displaystyle f(x)={\frac {1}{2}}+{\frac {1}{2}}\tanh \left({\frac {x}{2}}\right)}
或
tanh
(
x
)
=
2
f
(
2
x
)
−
1.
{\displaystyle \tanh(x)=2f(2x)-1.}
推導過程如下:
tanh
(
x
)
=
e
x
−
e
−
x
e
x
+
e
−
x
=
e
x
⋅
(
1
−
e
−
2
x
)
e
x
⋅
(
1
+
e
−
2
x
)
=
f
(
2
x
)
−
e
−
2
x
1
+
e
−
2
x
=
f
(
2
x
)
−
e
−
2
x
+
1
−
1
1
+
e
−
2
x
=
2
f
(
2
x
)
−
1.
{\displaystyle {\begin{aligned}\tanh(x)&={\frac {e^{x}-e^{-x}}{e^{x}+e^{-x}}}={\frac {e^{x}\cdot \left(1-e^{-2x}\right)}{e^{x}\cdot \left(1+e^{-2x}\right)}}\\&=f(2x)-{\frac {e^{-2x}}{1+e^{-2x}}}=f(2x)-{\frac {e^{-2x}+1-1}{1+e^{-2x}}}=2f(2x)-1.\end{aligned}}}
標準邏輯斯諦函數的導數 稱為邏輯斯諦分布 密度,公式如下:
f
(
x
)
=
1
1
+
e
−
x
=
e
x
1
+
e
x
,
{\displaystyle f(x)={\frac {1}{1+e^{-x}}}={\frac {e^{x}}{1+e^{x}}},}
d
d
x
f
(
x
)
=
e
x
⋅
(
1
+
e
x
)
−
e
x
⋅
e
x
(
1
+
e
x
)
2
=
e
x
(
1
+
e
x
)
2
=
f
(
x
)
(
1
−
f
(
x
)
)
{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} x}}f(x)={\frac {e^{x}\cdot (1+e^{x})-e^{x}\cdot e^{x}}{(1+e^{x})^{2}}}={\frac {e^{x}}{(1+e^{x})^{2}}}=f(x){\big (}1-f(x){\big )}}
邏輯斯諦分布的均值為x0 ,變異數為π2 /3k 2 。
標準邏輯斯諦函數的不定積分 可用換元積分法 求得,令
u
=
1
+
e
x
{\displaystyle u=1+e^{x}}
,
f
(
x
)
=
e
x
1
+
e
x
=
u
′
u
{\displaystyle f(x)={\frac {e^{x}}{1+e^{x}}}={\frac {u'}{u}}}
,去掉積分常數 ,得到其不定積分:
∫
e
x
1
+
e
x
d
x
=
∫
1
u
d
u
=
ln
u
=
ln
(
1
+
e
x
)
.
{\displaystyle \int {\frac {e^{x}}{1+e^{x}}}\,dx=\int {\frac {1}{u}}\,du=\ln u=\ln(1+e^{x}).}
在人工神經網絡 中,它稱作線性整流函數 ,(縮放後)可視為平滑近似的斜坡函數 ,類似於邏輯斯諦函數(縮放後)是平滑近似的單位階躍函數 。
標準邏輯斯諦函數是簡單的一階非線性常微分方程式 的解:
d
d
x
f
(
x
)
=
f
(
x
)
(
1
−
f
(
x
)
)
{\displaystyle {\frac {d}{dx}}f(x)=f(x){\big (}1-f(x){\big )}}
邊界條件 為
f
(
0
)
=
1
/
2
{\displaystyle f(0)=1/2}
。該方程式是邏輯斯諦映射 的連續版本。注意倒數邏輯斯諦函數是簡單的一階線性常微分方程式的解。[ 6]
x
n
+
1
=
k
x
n
(
1
−
x
n
)
{\displaystyle x_{n+1}=kx_{n}(1-x_{n})}
是混沌理論 的一個模型。[ 7] [ 8] 這個函數對初始值和參數的變化很敏感,往往微小的變化會引起混沌。如圖所示,當x 1 =0.3 ,參數k 從0.1變到4時,系統變化很大。
當k 由0.1變到1時,曲線很快趨向於0
當k 繼續增加,曲線由0.3上升到 一個穩定值
k 繼續增加,曲線出現擺動,有2個穩定值。
k 繼續增加, 曲線相繼出現4個、8個、16個、32個....穩定值
k 增加到一個臨界值,系統進入混沌狀態。
k 再增加,系統突然垮塌。
x
n
+
1
=
k
x
n
(
1
−
x
n
2
)
{\displaystyle x_{n+1}=kx_{n}(1-x_{n}^{2})}
皮埃爾·弗朗索瓦·韋呂勒(1804–1849)
邏輯斯諦方程式的一個典型應用是族群(或人口)增長的通用模型(另見族群動態 ),最初由皮埃爾·弗朗索瓦·韋呂勒 在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 代表時間,該模型用以下微分方程式 表示:
d
P
d
t
=
r
P
(
1
−
P
K
)
,
{\displaystyle {\frac {dP}{dt}}=rP\left(1-{\frac {P}{K}}\right),}
其中常數r 為族群(人口)增長率 ,K 為環境承載力 。
方程式中,早期的幾乎無阻力的增長率來自+rP 。增長率r 代表族群(人口)數量P 在一個單位時間內的增長比例。後來,隨著人口的增長,第二項-rP2 /K 變得幾乎和第一項一樣大,族群P 內的個體之間開始爭奪某些關鍵資源(例如食物或生存空間)而相互干擾。這種對抗效應稱為「瓶頸」,由參數K 代表。競爭會降低總合增長率,直到P 停止增長(族群/人口成熟)。方程式的解(P0 為初始族群/人口數量)為
P
(
t
)
=
K
P
0
e
r
t
K
+
P
0
(
e
r
t
−
1
)
=
K
1
+
(
K
−
P
0
P
0
)
e
−
r
t
,
{\displaystyle P(t)={\frac {KP_{0}e^{rt}}{K+P_{0}\left(e^{rt}-1\right)}}={\frac {K}{1+\left({\frac {K-P_{0}}{P_{0}}}\right)e^{-rt}}},}
其中:
lim
t
→
∞
P
(
t
)
=
K
.
{\displaystyle \lim _{t\to \infty }P(t)=K.}
可以說,K 是P 的極限值,即經過無限長時間後(或在有限時間內近似),族群(人口)規模所能達到的最大值。須注意,只要初始值
P
(
0
)
>
0
{\displaystyle P(0)>0}
,無論取值多少,族群數量都會漸近環境承載力的值,包括
P
(
0
)
>
K
{\displaystyle P(0)>K}
的情況下。
生態學中有時稱一個物種 是r 策略或K 策略 的,這是指它們在自然選擇 過程形成的生物生命週期 策略。選取變量的量綱 ,使n 代表以環境承載力單位計的族群數量,
τ
{\displaystyle \tau }
代表以
1
/
r
{\displaystyle 1/r}
的單位計量的時間,得出無量綱微分方程式:
d
n
d
τ
=
n
(
1
−
n
)
.
{\displaystyle {\frac {dn}{d\tau }}=n(1-n).}
邏輯斯諦函數在統計學中有多種應用。例如,它們是邏輯斯諦分布 的累積分布函數 ,它們可用於模擬西洋棋 棋手在埃洛等級分系統 下擊敗對手的機率。以下是一些更具體的案例。
邏輯斯諦迴歸 使用邏輯斯諦函數來模擬一個事件的機率p 如何可能會受到一個或多個解釋變量 的影響:一個案例模型如下
p
=
f
(
a
+
b
x
)
,
{\displaystyle p=f(a+bx),}
其中x 為解釋變量,a 和b 為欲擬合的模型參數,f 為標準邏輯斯諦函數。
邏輯斯諦迴歸和其他對數線性模型 也常用於機器學習 。將邏輯斯諦函數推廣至多元輸入情景即為Softmax激活函數 ,用於多元邏輯迴歸 。
在醫學上,邏輯斯諦微分方程式可用於腫瘤生長的建模。這一用法可視為上述的生態學/人口學模型的延伸。以
X
(
t
)
{\displaystyle X(t)}
表示腫瘤在時間
t
{\displaystyle t}
的大小,其變化動態遵循
X
′
=
r
(
1
−
X
K
)
X
,
{\displaystyle X'=r\left(1-{\frac {X}{K}}\right)X,}
屬於以下類型:
X
′
=
F
(
X
)
X
,
F
′
(
X
)
≤
0
,
{\displaystyle X'=F(X)X,\quad F'(X)\leq 0,}
其中
F
(
X
)
{\displaystyle F(X)}
為腫瘤增殖率。
如果採用化療產生對數殺傷效果,則等式修改為
X
′
=
r
(
1
−
X
K
)
X
−
c
(
t
)
X
,
{\displaystyle X'=r\left(1-{\frac {X}{K}}\right)X-c(t)X,}
其中
c
(
t
)
{\displaystyle c(t)}
為治療引起的腫瘤死亡率。在理想化的極長的治療下,
c
(
t
)
{\displaystyle c(t)}
可模型化為週期為
T
{\displaystyle T}
的週期函數或(在持續的輸液治療下)常數函數,有
1
T
∫
0
T
c
(
t
)
d
t
>
r
→
lim
t
→
+
∞
x
(
t
)
=
0
,
{\displaystyle {\frac {1}{T}}\int _{0}^{T}c(t)\,dt>r\to \lim _{t\to +\infty }x(t)=0,}
即,如果平均治療引起的腫瘤死亡率大於基線增殖率,則疾病能被根除。當然,這是一個過於簡化的生長和治療模型(例如沒有考慮克隆抗性現象)。
在人群中未被免疫的新型傳染性病原體,通常會在早期呈指數級傳播,有大量易感個體尚未被感染。例如2020年初,導致2019冠狀病毒病 的SARS-CoV-2 病毒在多國的感染過程中呈現出指數級增長。[ 12] 此後,易感宿主減少(持續感染直到超過群體免疫閾值)或通過社交距離措施減少潛在宿主的被傳染機率等因素,可能使呈指數增長的傳染曲線首先線性化,然後趨緩,達到最大值。[ 13]
邏輯斯諦函數或相關的函數(例如龔珀茲函數 )通常以描述性或現象學方的式使用,因為它們非常符合早期的指數上升,也符合隨著人群形成群體免疫而最終趨於平穩的趨勢。它與流行病的實際模型不同,後者試圖根據大流行的動態(例如接觸率、潛伏期、社交距離等)來描述感染狀態。不過,一些簡單的模型有邏輯斯諦解。[ 14] [ 15] [ 16]
流行病模型中的廣義邏輯斯諦曲線 (Richards增長曲線)
廣義邏輯斯諦函數 (又稱Richards增長曲線)已應用於對COVID-19爆發的早期階段建模。[ 17] 研究者將廣義邏輯斯諦函數擬合到累計感染病例數(稱為傳染軌跡)。文獻中對廣義邏輯斯諦函數有不同的參數化。一種常用的形式是:
f
(
t
;
θ
1
,
θ
2
,
θ
3
,
ξ
)
=
θ
1
[
1
+
ξ
exp
(
−
θ
2
⋅
(
t
−
θ
3
)
)
]
1
/
ξ
{\displaystyle f(t;\theta _{1},\theta _{2},\theta _{3},\xi )={\frac {\theta _{1}}{[1+\xi \exp(-\theta _{2}\cdot (t-\theta _{3}))]^{1/\xi }}}}
其中
θ
1
,
θ
2
,
θ
3
{\displaystyle \theta _{1},\theta _{2},\theta _{3}}
取實數,
ξ
{\displaystyle \xi }
為正實數。曲線
f
{\displaystyle f}
的靈活性由
ξ
{\displaystyle \xi }
賦予:(i)若
ξ
=
1
{\displaystyle \xi =1}
,則曲線衰減為邏輯斯諦函數,(ii)若
ξ
{\displaystyle \xi }
收斂至0,則曲線收斂至龔珀茲函數。在傳染病模型中,
θ
1
{\displaystyle \theta _{1}}
,
θ
2
{\displaystyle \theta _{2}}
和
θ
3
{\displaystyle \theta _{3}}
分別代表傳染病最終的規模、感染率和滯後期。見右側的範例的傳染軌跡,其中
(
θ
1
,
θ
2
,
θ
3
)
{\displaystyle (\theta _{1},\theta _{2},\theta _{3})}
設定為
(
10000
,
0.2
,
40
)
{\displaystyle (10000,0.2,40)}
。
受COVID-19嚴重影響的40個國家的傳染曲線外推和截至5月14日的母體(人口)平均值
在流行病學建模中,使用類似廣義邏輯斯諦函數的增長函數的好處之一是它相對容易應用於多級模型 框架,其中來自不同地理區域的資訊可以匯總在一起。
自催化反應 中,反應物和產物的濃度遵循邏輯斯諦函數。例如燃料電池 陰極中不含鉑族金屬 (PGM-free)的氧還原反應催化劑的劣化遵循邏輯斯諦衰減函數,[ 18] 表明這是一種自催化分解機制。
費米子 在熱平衡系統的能量狀態上的統計分布遵循邏輯斯諦函數。特別地,根據費米-狄拉克統計 ,它是每個可能的能級被一個費米子占據的機率分布。
語言學中,邏輯斯諦函數可用於對語言變化 進行建模:[ 19] 一種最初處於邊際地位的新詞隨著時間的推移開始傳播得更快,然後傳播速度隨著其普及而減慢。
邏輯斯諦函數可用於描繪一項發明創新在其生命週期內擴散的過程 .
韋呂勒首先提到算術級數和幾何級數,並將幾何增長曲線稱為「對數曲線 / logarithmic curve」(但須注意,現代術語稱之為「指數 曲線 / exponential curve」)。然後他把他提出的曲線叫做邏輯斯諦(logistic),以區別於對數(logarithmic),並在他論文的圖中比較了「對數」曲線和邏輯斯諦曲線。
在古希臘,λογῐστῐκός 指實務的計算與會計,而ἀριθμητική (arithmētikḗ )指對於數的理論或哲學研究。然而在英語中,表示算術的詞反而是源自ἀριθμητική 的arithmetic。
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.
Garnett P. Williams Chaos Theory Tamed chapter 10
Edgar Peters Chaos and Order in the Capital Market p7
Villalobos-Arias, Mario. Using generalized logistics regression to forecast population infected by Covid-19. 2020. arXiv:2004.02406 [q-bio.PE ].
Bod, Hay, Jennedy (eds.) 2003, pp. 147–156