宋丹丹，王晶囡，楊光旭，李佳思，楊嘉欣

（哈爾濱理工大學(xué) 應(yīng)用科學(xué)學(xué)院，黑龍江 哈爾濱 150080）

?

基于數(shù)值積分的傳染病SIRS模型參數(shù)估計

宋丹丹，王晶囡，楊光旭，李佳思，楊嘉欣

（哈爾濱理工大學(xué) 應(yīng)用科學(xué)學(xué)院，黑龍江 哈爾濱 150080）

摘要：針對具有免疫的傳染病SIRS模型，利用三次Hermite插值函數(shù)及數(shù)值積分公式，基于患病的各個種群人數(shù)估計值的誤差最小原則，將參數(shù)估計問題轉(zhuǎn)化為非約束優(yōu)化問題.將數(shù)據(jù)帶入后可得關(guān)于模型參數(shù)的多項式，為求得該式最小值，將其分別對各個參數(shù)進行微分，得到關(guān)于模型參數(shù)的非線性方程組.使用最速下降法獲得較為合理與精確的初值，在該初值的基礎(chǔ)上利用牛頓法對非線性方程組進行求解，得到了該模型的高精度參數(shù)估計值.并對計算結(jié)果進行數(shù)值仿真，數(shù)值仿真實驗表明，所給出的參數(shù)估計方法能夠較為精確地估計出相應(yīng)參數(shù)值.

關(guān)鍵詞：SIRS模型；數(shù)值積分；牛頓法；數(shù)值仿真；參數(shù)估計

對于各類新型突發(fā)傳染病，人類對它的防治還處于初級階段，如何有效地從宏觀上了解和掌握這些流行病的傳播規(guī)律，控制傳染病的蔓延就顯得尤為重要［1-8］.本文研究現(xiàn)實復(fù)雜情形下（包含非線性傳染率，有隔離措施，群外個體遷入，生育與死亡以及疾病可以水平或者垂直傳播等）SIRS傳染病模型的參數(shù)估計問題.采用三次Hermite插值函數(shù)，得到具有五階精度的數(shù)值積分公式，利用該公式對SIRS傳染病模型進行離散化，將參數(shù)估計問題轉(zhuǎn)化為無約束的優(yōu)化問題，使用Matlab軟件進行操作，最終得到估計誤差最小的參數(shù)估計值.

1　SIRS傳染病模型

在文獻［1］的基礎(chǔ)上，本文考慮了感染者具有免疫作用并對該文獻中傳染病模型進行適當(dāng)刪減，建立傳染病模型其中：S= S（t），I= I（t），R= R（t）分別表示易感個體、染病個體（又分為未被隔離染病個體I1和隔離染病個體I2，即I= I1+ I2）和康復(fù)個體的數(shù)量；ε表示隔離率；β（N）表示傳染病通過有效接觸進行傳播的傳播系數(shù)；a表示接種比例；As表示易感個體的輸入率；b表示自然生育率；d表示自然死亡率；d表示染病個體的死亡率；γ表示治愈率.初始條件為S（ t0）= S0，I1（ t0）= I10，I2（ t0）= I20，R（t0）= R0.As，a，b，c，d，β（N），ε，γ為待求參數(shù).

2　三次Hermite插值多項式與數(shù)值積分

一般地，求三次Hermite插值函數(shù)要求給定x0，x1和相應(yīng)的函數(shù)值F（x0），F(xiàn)（x1）以及微商值F'（x0），F(xiàn)'（x1） .

其中：ξ∈［x0，x1］.

3　SIRS傳染病模型的參數(shù)估計

方程（1）為非線性微分方程組，很難求出其解析解，本文采用數(shù)值方法進行參數(shù)估算.對于S（t），I1（t），I2（t），R（t）的一組測量數(shù)據(jù)S=｛S（ti）｝，I1=｛ I1（ti）｝，I2=｛ I2（ti）｝，R=｛ R（ ti）｝（i =1，2，…，n），令則參數(shù)估計問題轉(zhuǎn)化為無約束非線性優(yōu)化問題，求得f的最小值minf即為目標(biāo)函數(shù).其中：S，I1，I2， R為觀測值；為相應(yīng)的估計值；表示泛函中的2-范數(shù)的平方.為得到函數(shù)（3）的具體表達式，利用數(shù)值積分公式可得

其中：A，B，C，D分別為式（4）～（7）所表示的估計值.

將數(shù)據(jù)代入式（8）中即可得到關(guān)于8個參數(shù)的多項式，記為f（a，b，d，d，As，ε，β（N），γ），為求得f（ a，b，d，d， As，ε，β（N），γ）的最小值，將它分別對8個參數(shù)進行求導(dǎo)，可得8個非線性方程所組成的非線性方程組，記為

采用牛頓法求解非線性方程組（9）.牛頓法求解非線性方程組就是將非線性問題線性化，為減少計算次數(shù)和得到較為精確的結(jié)果，選取合理初值就顯得尤為重要，這里采用最速下降法選取合理初值，記為

4　數(shù)值實驗與結(jié)果分析

以2004—2013年江蘇省相關(guān)麻風(fēng)病數(shù)據(jù)（數(shù)據(jù)來源于公共衛(wèi)生科學(xué)數(shù)據(jù)中心，相關(guān)網(wǎng)站為http：//www.phsciencedata.cn/）為例，對本文所建模型的未知參數(shù)進行估計，時間區(qū)間為［0，9］，采樣步長為Δt = 1 .利用軟件得到的參數(shù)計算結(jié)果為a=0.31，b =0.87，d =0.08，ε=0.23，γ=0.93，β（N）=0.65，As=0.26，d =0.05.為驗證本文方法的有效性與精確性，將參數(shù)估計值代入模型（1）中，采用Runge-Kutta法求得其數(shù)值并進行數(shù)據(jù)仿真，得到I1（t），I2（t），R（t）的擬合曲線（見圖1）.由圖1可以看出，該方法所得計算結(jié)果與實際觀測數(shù)據(jù)吻合程度相對較高.

圖1 I1（t），I2（t），R（t）仿真結(jié)果

由于查找的數(shù)據(jù)十分有限且時間跨度較大，部分缺失值的填充參入過多主觀因素，使得數(shù)據(jù)擬合效果與預(yù)想有一定的偏差，對本文計算方法的精確性的驗證產(chǎn)生了影響，但并不能否認該方法的正確性與實用性.若有相對較為精確的數(shù)據(jù)，本文計算方法將會得到較好的驗證.

5　結(jié)論

對于SIRS傳染病模型的參數(shù)估計問題，本文運用了基于數(shù)值算法的參數(shù)估計，首先結(jié)合三次Hermite插值多項式得到一個高精度的數(shù)值積分公式，利用該公式對傳染病模型進行離散化，進而將參數(shù)估計問題轉(zhuǎn)化為無約束優(yōu)化問題，將數(shù)據(jù)代入后得到關(guān)于各個參數(shù)的多項式.為求得該式的最小值，將其分別對各個參數(shù)進行求導(dǎo)，進而得到關(guān)于各個參數(shù)的非線性方程組，利用最速下降法取得參數(shù)的合理初始值，最后利用牛頓法求得各個參數(shù)值.并針對江蘇省患麻風(fēng)病情況進行了仿真對比，仿真結(jié)果表明，該算法所得結(jié)果與實際數(shù)據(jù)統(tǒng)計結(jié)果具有較好的擬合度.因此，本文的參數(shù)估計算法在估計類似傳染病模型參數(shù)上是可行和可信的.

參考文獻：

［1］ 劉利斌，劉翠萍，張永，等.基于粒子群算法的捕食者-食餌模型的參數(shù)估計［J］.純粹數(shù)學(xué)與應(yīng)用數(shù)學(xué)，2016，32（1）：19-25

［2］ 孫有發(fā)，郭旭沖.現(xiàn)實復(fù)雜情形下的SIRS型傳染病模型及其控制策略［J］.系統(tǒng)仿真學(xué)報，2010，22（1）：195-200

［3］ 張海峰，傅新楚.含有免疫作用的SIR傳染病模型在復(fù)雜網(wǎng)絡(luò)上的動力學(xué)行為［J］.上海大學(xué)學(xué)報：自然科學(xué)版，2007，13（2）：189-192

［4］ 燕振剛，胡賀年，李廣.基于粒子群算法的Richards模型參數(shù)估計和算法有效性［J］.計算機應(yīng)用，2014，34（10）：2827-2830

［5］ 曹懷火，李海燕，伏升茂.基于混合加速粒子群算法的捕食者-食餌模型參數(shù)估計［J］.生物數(shù)學(xué)學(xué)報，2013，28（3）：553-557

［6］ 馬新生，王來群，胡文玉.基于潛周期模型的兩種群食餌-捕食者模型的參數(shù)估計［J］.南昌大學(xué)學(xué)報，2008，30（2）：134-137

［7］ 柳寅，馬良，黃鈺.基于模糊粒子群算法的非線性函數(shù)優(yōu)化［J］.上海理工大學(xué)學(xué)報，2012，34（4）：314-317

［8］ 李斌，李文峰.基于仿真的優(yōu)化的粒子群算法參數(shù)選取研究［J］.計算機工程與應(yīng)用，2011，47（33）：23-24

Parameter estimation of the SIRS infectious disease model based on numerical integration

SONG Dan-dan，WANG Jing-nan，YANG Guang-xu，LI Jia-si，YANG Jia-xin
（School of Applied Science，Harbin University of Science and Technology，Harbin 150080，China）

Abstract：The third-order Hermite interpolation function and numerical integral formula are applied to an SIRS infectious disease model.Based on the principle of minimum error， the parameter estimation problem in the model is transformed to the unconstrained optimization problem and a polynomial can be obtained after the data is substituted into the formula.In order to get the minimum value of the polynomial， the partial derivatives of the polynomial with respect to each parameter are obtained after which a system of nonlinear equations related to the model is obtained. The steepest descent method is used to get a set of accurate and reasonable initial values，based on which Newton method is used to solve the nonlinear equations and then the specific parameter values are obtained.Finally，the equilibriums of SIRS model and the effectiveness of general parameter estimation are validated by numerical simulations.

Key words：SIRS model；numerical integration；Newton method；numerical simulations；parameter estimation

中圖分類號：O29

文獻標(biāo)識碼：A

doi：10.3969/j.issn.1007-9831.2016.05.004

文章編號：1007-9831（2016）05-0010-05

收稿日期：2016-04-15

基金項目：哈爾濱理工大學(xué)大學(xué)生創(chuàng)新項目（201512）；哈爾濱理工大學(xué)教學(xué)改革項目（320150023）

作者簡介：宋丹丹（1991-），女，黑龍江齊齊哈爾人，在讀本科生.

通信作者；王晶囡（1978-），女，山東平度人，副教授，博士，從事微分方程研究.E-mail：wangjingnan@hrbust.edu.cn