譚景寶 李明
(1. 合肥幼兒師范高等??茖W(xué)?;A(chǔ)部, 合肥 230013; 2. 安徽文達信息工程學(xué)院, 合肥 231201)
近年來,我國政府對計劃生育政策逐步作了重大調(diào)整。自2015年開始,已全面實施二孩政策。人口問題是一個復(fù)雜的系統(tǒng)工程,生育政策的變化,必然會使人口數(shù)量、人口質(zhì)量和人口結(jié)構(gòu)發(fā)生變化。在這種新形勢下,原有的人口系統(tǒng)模型研究也需隨之調(diào)整,以適應(yīng)新的人口發(fā)展趨勢。本次研究將在現(xiàn)有人口系統(tǒng)模型的基礎(chǔ)上[1-4],探討二孩政策下的人口系統(tǒng)模型,建立和補充新的人口系統(tǒng)模型。
簡單描述生物種群模型:
(1)
在方程組(1)中,通過微分公式,對種群群落數(shù)量N(t)變化的速率進行分析,即分析離散狀態(tài)下的時間與種群數(shù)量關(guān)系。其中,b為種群群落的出生系數(shù);d為種群群落的死亡系數(shù);γ為種群群落數(shù)量N(t)的增長系數(shù)。
Malthus人口模型:
(2)
方程組(2)與方程組(1)在函數(shù)表達形式上有著一定的相似性。p(t)代表的是不同時間的人口增長系數(shù),即相關(guān)研究內(nèi)容方面的種群樣本總數(shù)。然而,方程組(1)中種群群落的增長系數(shù)在種群總量增加的情況下卻不斷降低,這不利于種群群落人口數(shù)量的研究。針對此問題,可應(yīng)用Logistic分析法構(gòu)建相關(guān)模型,使其與生態(tài)學(xué)中的人口模型研究相適應(yīng)。具體模型如方程組(3)所示:
(3)
在方程組(3)中,某一時刻t的種群群落總數(shù)為N(t),種群環(huán)境生存常數(shù)為k,種群生活環(huán)境的實際容量(也稱種群增長系數(shù)) 為γ。相對于方程組(1)(2)來說,方程組(3)的適用性較強。在此基礎(chǔ)上,推導(dǎo)出與我國人口發(fā)展相適應(yīng)的模型,如方程組(4)所示:
(4)
在方程組(4)中,p(t)表示在t時刻某一穩(wěn)定環(huán)境中的人口總數(shù)情況。其中,δ0為環(huán)境總量,ω0則為其人口增長系數(shù)。利用假設(shè)法,假設(shè)穩(wěn)定生態(tài)環(huán)境中人口出生、死亡與年齡均無相關(guān)性。方程組(4)能夠在一定程度上應(yīng)對傳統(tǒng)線性人口動力學(xué)模型對人口數(shù)量增長與自然資源的適應(yīng)性問題。由此,得到對應(yīng)的適應(yīng)模型:
假設(shè)增量h>0,用方程表示為:
(5)
dp(a,t)代表的是t p(a,0)=p0(a) (6) 這里忽略了遷移的影響,仿照方程組(4),假設(shè): d(a,t)=[μ(a,t)-kN(t)]p(a,t) (7) 在方程(7)中,μ(a,t)為死亡率相對函數(shù),k是環(huán)境容量,N(t)為在t時刻下社會人口總量。設(shè)p(a,t)可微,則有: (8) 結(jié)合式(6)、(7)、(8),得到式(9): =-[μ(a,t)-kN(t)]×p(a,t) (9) 考慮遷移因素,則有式(10): =-[μ(a,t)-kN(t)]×p(a,t)+f(a,t) (10) 于是,得到非線性人口動力學(xué)方程: (11) 其中:p0(a)、β(t)、h(a,t)、k(a,t)分別代表最大育齡總數(shù)、規(guī)格化條件、生育模式、總和生育率。方程組(11)有別于線性動力學(xué)方程,式中社會死亡率和總?cè)藬?shù)相關(guān),即人口的生存依賴于人口總?cè)萘俊?/p> 我們提出的非線性人口動力學(xué)方程,可以更準(zhǔn)確地反映社會人口的出生與死亡與其總?cè)萘康年P(guān)系。在方程(6)中,令 d(a,t)=μ(a,t,N(t))p(a,t) (12) 為了區(qū)分,稱μ(′, ″ ,′)為死亡率,設(shè)β(a,t,N(t))為t時刻年齡為a的個體在平均單位時間內(nèi)的平均生育量,即生育率。β依賴N(t),則邊界條件為: (13) 設(shè)人口初始密度函數(shù)為 p(a,0)=p0(a) (14) 結(jié)合式(6)(12)(13)(14),得到方程組(15): (15)