一類受媒體影響的傳染病模型的研究

劉玉英，肖燕妮

摘要：目的：在突發(fā)性傳染病爆發(fā)時，大眾媒體報道對疫情的有效報道極大的影響著群眾的心理和行為，這種行為改變（如采取防護措施）對傳染病的感染率有著重要的影響。在傳染病爆發(fā)初期，媒體報道對傳染病感染率的削弱作用尤為明顯。隨著疫情的擴散，這種在感染率上的削減作用會越來越弱。隨著疫情繼續(xù)蔓延，易感者人群已不能形成新的習慣來避免感染。這意味著感染者數目有一個臨界值，當達到這個臨界值時，易感者人群的行為變化達到最大，也即媒體報道的削減作用達到最大；隨著染病者人數的繼續(xù)增加，易感者人群仍然需要參加必要的社會活動與別人進行接觸，媒體影響因子函數達到飽和。本文構造了一個能刻畫媒體影響飽和性的函數，探究大眾媒體報道對傳染病傳播與控制的作用，特別是研究飽和型的媒體反映函數是如何影響傳染病動力學行為。方法：我們建立了感染率為非光滑函數的傳染病模型，該非光滑函數（即分段函數）能準確的刻畫媒體報道對感染率的飽和作用影響。引入了Lambert W 函數來明確計算出平衡態(tài)表達式，得到了模型的所有平衡態(tài)，并利用了Lambert W 函數在區(qū)間上的單調性來分析了平衡點局部漸進穩(wěn)定性的過程，根據基本再生數的大小分析了各平衡態(tài)的局部穩(wěn)定性。通過排除極限換的存在性證明了平衡態(tài)的全局漸近穩(wěn)定性，討論了媒體影響傳染病模型的動力學性態(tài)。結果：計算了系統(tǒng)的基本再生數，并得到當R0＜1 時，疾病消除平衡態(tài)是全局漸近穩(wěn)定的；當1＜R0且1＜Rc＜emIc+βIc/d 時，子系統(tǒng)S1的地方病平衡態(tài)E1是全局漸近穩(wěn)定的，當emIc+βIc/d＜Rc時，子系統(tǒng)S2的地方病平衡態(tài)E2是全局漸近穩(wěn)定的，并用圖形分別顯示了不同參數下地方病平衡態(tài)E1和E2的全局穩(wěn)定性。結論：較具有飽和性的媒體影響函數與不具飽和性的媒體影響函數對傳染病動力學性質的影響，具有飽和傳染率的系統(tǒng)的地方病平衡態(tài)中感染者數目要比采用光滑感染因子的模型的感染者數目多。這意味著如果忽視媒體影響的飽和效應，媒體對傳染病感染率的影響將會被高估。

來源出版物：應用數學和力學, 2013, 34(4): 399-407

入選年份：2017

基于非等時距加權灰色模型與神經網絡的組合預測算法

韓晉，楊岳，陳峰，等

摘要：目的：在實際生活和工程應用中，有很多數據是采用非等時間間距采樣，對于這一類數據的預測，和傳統(tǒng)的等間距數據預測方法不同，需要考慮到時間序列對數據發(fā)展的影響。針對這一問題，傳統(tǒng)的灰色系統(tǒng)理論將一切隨機變量都視為在一定范圍內變化的灰色量，用數據處理的方法找數據間的規(guī)律，將原始數據看成灰色量，對其進行一次累加生成處理后，所得到的數據序列將呈現指數規(guī)律變化，從而實現預測。但是灰色模型對波動性比較大的數據序列，難以把握其波動性，造成較大的預測誤差，而神經網絡模型在處理隨機性、非線性數據時有很大的優(yōu)越性，對數據的波動性預測結果較好?；谝陨戏治?，提出的基于非等時距加權灰色模型與神經網絡的組合預測模型可以避免單一模型丟失信息的缺憾，減小隨機性，提高非等時間間距采樣數據的預測精度。方法：針對非等時距的數據，設計的組合預測算法基本思想如下：首先對原始數據進行預檢驗，通過級比計算對原始數據進行可靠性判斷，決定是否進行平移處理，然后利用非等時距灰色預測模型進行初步預測，并且對該模型中累加序列進行初值和背景值優(yōu)化，利用最小二乘法估計灰色模型的常參數，進行加權處理，求出最終的時間響應函數，得到初步預測值。將初步得到的預測值與實測值比較，求出殘差，利用BP 神經網絡對殘差序列進行預測修正，并還原到原始序列，獲取最終的組合預測值。結果：本文針對不等時間間隔數據預測的特點，將灰色預測理論與神經網絡有機結合，建立了一種非等時距加權灰色模型與神經網絡的組合預測算法，對模型精度及后驗差比值進行計算，通過實例驗證本文算法的有效性。以某防洪大堤沉降監(jiān)測數據及不同溫度條件下的鈦合金疲勞強度預測為例，本文算法的模型精度分別為98.279%和99.530%，后驗差比值分別為0.0259 和0.0056，模型精度等級均為1 級，高于傳統(tǒng)預測算法。結論：本文提出的組合預測算法，采用原始數據序列的平均值作為累加序列初值從而對初值進行優(yōu)化，綜合利用了所有歷史數據的特點，使結果更加準確；在計算累積序列微分方程背景值時，將累積曲線看成指數函數的形式，利用積分面積代替背景值，而不是簡單地使用梯形公式，提高了預測精度；在傳統(tǒng)灰色預測模型的基礎上，引入了加權的概念，通過定義權矩陣來表征歷史數據的可靠性，對離預測時間越近的數據值賦予更大的權值，從而使預測結果更加可靠；利用神經網絡對灰色預測結果進行殘差修正，可以較好地解決數據發(fā)展趨勢中復雜的不確定性問題，克服了單一灰色模型預測的不足；既可應用于非等時距采樣數據的預測，也可應用于等時距采樣數據的預測；本文算法與灰色預測模型相比，計算精度高，且誤差可控；與神經網絡預測模型相比，計算量小，在小樣本情況下也可達到較高的預測精度。

來源出版物：應用數學和力學, 2013, 34(4): 408-419

入選年份：2017

基于PFC3D 的露天礦邊坡爆破過程模擬 及穩(wěn)定性研究

崔鐵軍，馬云東，王來貴

摘要：目的：為研究爆破在露天礦邊坡內發(fā)展過程，基于能量守恒理論，假設爆炸時刻產生的能量全部由爆點周邊一定范圍內的巖體承受，并部分轉化為動能；能量在碎裂的巖塊中傳遞、吸收，最終達到平衡，爆炸過程結束。為實現上述過程，使用顆粒流理論的PFC3D 作為模擬平臺，模擬在露天礦邊坡內，不同高度、埋深和裝藥量的單孔爆破過程，并對爆破后邊坡穩(wěn)定性進行論述。為采礦工程提供安全提供指導。方法：根據PFC3D爆破計算步驟圖實施爆破。從兩方面對PFC3D 構造的邊坡模型中顆粒進行爆破初始狀態(tài)的設置。一是設炸藥的能量轉化為顆粒動能。設R＝0.07 m，考慮爆炸區(qū)域劃分，壓縮區(qū)范圍較小，為[0.21 m, 0.59 m]；振動區(qū)主要吸收殘余能量起阻尼作用，不發(fā)生斷裂；破裂區(qū)主要承受爆炸能量，碎裂并飛濺，炸藥化學能轉化為動能。故設壓縮區(qū)、破裂區(qū)、振動區(qū)能量分配為10%、80%、10%。壓縮區(qū)集中在炸點的顆粒內不參與爆破。振動區(qū)不發(fā)生破碎所以所吸收的能量不轉化為動能。動能主要集中在破裂區(qū)。通過炸藥化學能轉化為動能，根據動能公式，將每個顆粒得到的動能轉化為顆粒的爆炸初始瞬間速度，實現爆炸模擬。二是考慮爆炸瞬間釋放的氣體等沖擊波對巖體造成的碎裂和劣化作用，對破裂區(qū)范圍內的顆粒之間鏈接力設為0，并調整顆粒的接觸連接（contact-bond）和平行連接（parallel-bond）使其只承受壓力，不承受拉力和剪力。對振動區(qū)只產生劣化作用，分配10%能量完成巖石結構破壞。故根據顆粒位置到破裂區(qū)邊緣距離線性改變顆粒的鏈接屬性，靠近破裂區(qū)的顆粒鏈接屬性減小且接近破裂區(qū)內顆粒屬性，遠離破裂區(qū)的顆粒屬性接近正常。規(guī)定超過300 R 的顆粒不受爆破影響。結果：以中國，遼寧，阜新海州礦邊坡為例進行爆破模擬研究。A3 點100 kg TNT 為例說明爆炸過程。動力計算中1 步約等于1 ms。0 步時是對邊坡內顆粒作初始能量分配。從爆炸點A3 向外輻射速度矢量，接近爆破點的顆粒速度矢量大，外圍顆粒的速度小。200 步時，靠近邊坡自由面的顆粒向外飛出，速度矢量只受重力影響減小不明顯；其余部分方向的顆粒速度矢量減小明顯。500 步時邊坡方向顆粒繼續(xù)向外移動，其余方向顆粒將動能轉化為彈性勢能。此過程出現唯一一次速度矢量的震蕩。從720 步到1340 步保持了顆粒的返回震蕩特征，顆粒向爆炸中心移動。前1340 步爆炸是主導能量，重力幾乎不起作用。從1340 開始，爆炸能量逐漸消散，爆點上方沒有飛散出去的顆粒開始由于重力作用塌方。到14090 步大范圍的顆粒塌落已經結束，顆粒開始進行局部調整。26590 步顯示了最后的滾石顆粒的模擬，砂巖坡腳部分已經坍塌。計算達到平衡，邊坡自由面重新變得平緩。下層砂巖缺失一部分顆粒，這些顆粒就是被炸除的。結論：使用PFC3D 結合爆炸區(qū)劃理論的爆炸模擬是合理的，應考慮兩個方面。爆炸過程分3 個階段，第1 階段主要是爆炸沖擊起主導作用，有速度矢量的回蕩現象；第2 階段是重力占優(yōu)勢的上覆巖層塌落過程；第3 階段顆粒下滑局部調整，最后平衡。在經歷時間上，后一個階段比前一個階段大一個數量級。上述各種爆破后上層砂巖是穩(wěn)定的，即坡頂是穩(wěn)定的。下層砂巖和砂質泥巖會受到一定程度的破壞，但是沒有發(fā)生大面積的滑坡，是可控范圍內的。

來源出版物：應用數學和力學, 2014, 35(7): 759-767

入選年份：2017