1、1结核病传播的数学模型代胜吉 焦胜博 朱益飞摘要本文建立结核病传染的改进 K-M 仓室模型，并对广州市结核病传播情况进行了分析和预测；讨论了不同预防和控制措施对结核病传染的影响；最后考虑了人口流动对模型的影响，将模型进一步完善。首先，将人群分为 7 类并确定了他们之间的关系，建立了改进型 K-M 仓室模型。列出了传染病动力学微分方程组和差分方程，计算并绘出了广州市地区结核病传染相关指标 （至每一季度季度累计染病人数） 、 （至每一季度累计无法治愈人数）21I 3I的理论曲线，将其与实际曲形进行了比较，得到了比较满意的结果。根据 2008 年广州市第一季度初始值预测了 2008 至 2012 年
2、广州市结核病累计染病人数（ ） ，得出21I结论：到 2012 年，广州累计结核染病人数将达 168940 例，累计不可治愈人数将达173660 例。接着，分析了政府采取的措施对结核病传染的影响。根据各措施对模型造成的不同影响将措施分成四类并分别分析其对结核病传播的影响，提出了改进意见：（1）着重建立和完善结核病监控体系，切实实行好各项措施；（2）扩大对结核病治疗的范围，增加用于治疗结核病的经费，提高对耐药性结核病的治疗力度；（3）艾滋病的检测和治疗对结核病流行趋势的抑制效果不能立刻体现，因此在采取各项预防和控制结核病措施的同时还应常备不懈地坚持对艾滋病的检测和治疗；（4）在控制艾滋病增长的同
3、时，扩大艾滋病免费治疗的范围，对艾滋病和结核病交叉感染的治疗实施政策优惠。最后，将流动人口这一因素考虑到模型之中，对问题一中建立的模型进行进一步改进。经计算得出结论：（1）流动人口对结核病传染的影响是非常大的， 人口的流动会极大促进结核病的传播。在考虑外来流动人口的情况下，结核病流行程度比无外来流动人口的情况严重，2005 年末前一种情况的总染病人数（ 145660 人）已是后一种情况（43362 人）的约 3.36 倍，且这种倍数关系随着时间推移有增大的趋势。 （2）人口流动还会影响正常的结核病预防控制措施的实施效果。2关键词 K-M 仓室改进模型 结核病 预防控制 流动人口一、问题的重述结
4、核病是一种古老的经呼吸道传播的慢性传染病，在历史上流行猖獗，被称为白色瘟疫。仅从 1882 年著名的德国科学家郭霍氏发现结核菌以来，至少有约 2 亿人被结核病夺去了生命。尽管上世纪 50 年代有效的抗结核药物问世以后，形势发生了根本的改变，但是进入八十年代以后，由于各国政府的忽视、移民、多耐药结核病和艾滋病等因素的影响，在全球，特别是在发达国家又重新出现了上升的趋势。我国是全球 22 个结核病高负担国家之一，结核病人数位居世界第二位，仅次于印度。我国卫生部于 1979 年至 2000 年在全国先后开展了四次结核病调查，结核病发病数居法定传染病之首。2010 年我国进行第五次结核病流行病调查。2
5、000 年结核病流行病调查结果分析：我国结核病疫情还相当严重，部分地区有蔓延趋势， 表现为结核菌感染人数多、现患肺结核病人多、结核病死亡人数多、农村结核病人多和传染性肺结核病疫情居高不下。结核病有着独特的传播规律。它的流行与其它的传染病非常不同。它的病情发展很慢、在疾病的潜伏期和活化期采取完善的措施完全可以根治，但如果治疗不当或治疗不彻底而发展为耐药性结核就会成为不治之症，而且结核病病愈后不会获得免疫力。目前，在全球范围内，由于耐药结核菌的产生与扩展，结核菌与人体免疫缺陷病毒(HIV)的双重感染以及许多国家结核病控制规则的不完善，使得全球结核病商情呈明显上升趋势。在全球所有传染性疾病中，结核病
6、已成为成年人的首要死因。结核病流行的另一在特点就是其潜伏期很长，而且因人而异，有的感染者也可能终生都不会发病。查询相关资料解决以下四个问题：1、建立你们认为能反应我国或某一地区(不限中国国内)的结核病传播规律的传染病动力学模型，并对其结核病流行趋势进行中长期预测；2、通过模型分析，对现有的预防和控制措施进行评价，并分析不同的预防和控制措施对结核病传播的作用效果，也可提出自己的最佳预防控制策略；3、尝试将人口流动特点(如：我国的农民工在城乡间的流动，或不同地区国家之间的人口迁移)等因素对结核病传播和发作的影响予以定量化描述，并引入到所建模型3之中，寻求对结核病传播和预防控制的影响。4、 给当地报
7、刊写一篇通俗短文，说明对结核病进行预防控制的重要性。二、模型假设1、非耐药性结核病治疗周期极短，且治疗后无传染性；2、认为未接种疫苗人群即为易感染人群；3、若艾滋病与结核病交叉感染则无法治愈；4、不考虑人口流动对艾滋病增长的影响；5、不考虑所研究地区的出生死亡率；6、每年所研究地区的净流入人数为定值；7、选定地区的人口年龄结构在所研究时间段里基本不变；8、只有儿童必须依政策进行结核病免疫接种。三、符号说明： 易感人群数（未接种疫苗）)(1tS： HIV/AIDS 人数2： 未及时治疗或治疗不当的发病人数)(1tI： 发病人群中治疗治愈人数2： 不可治愈人数（耐药性、HIV 交叉感染）)(3tI
8、： 所研究地区总人数N： 未及时治疗及治疗不当者所占总感染细菌者的比率： 每个季度每个染病者平均作 次能引起传染的接触CC： 人口增加影响系数： 易感人群受传染的传染率r： HIV 感染者中感染肺结核细菌的比率2： 每季度 HIV 的死亡人数1R： 每季度肺结核的死亡人数： 单位季度 HIV 增长率4： 每季度净流入广州的外来人数A： 每季度净流入广州的外来人数的患病比率p四、问题 1 初步模型建立和分析4.1 初步模型的建立在建立初步模型时，考虑到结核病具有不同于一般传染病的特点：（1）未及时治疗或治疗不当的结核病患者会发展为耐药性结核无法治愈；（2）HIV 患者易感染结核病且无法治愈；基于
9、以上特点，我们选择对传统的 K-M“仓室”( compartment) 模型进行分析并及加以改进。4.1.1 人群的分类:（1）不易感人群（接种疫苗不发病人群） 0S（2）易感人群 )(1tS（3）HIV/AIDS 人群 2（4）未治疗和治疗不当者 )(1tI（5）无法治愈者（耐药性、AIDS 交叉感染等） )(3tI（6）治愈者 )(2tI（7）死亡者 R以上几类人群的关系如改进后的仓室模型框图 1 所示：5HIV/AIDS 人群 2S未治疗和治疗不当 )(1tI无法治愈（耐药性、AIDS交叉感染等） )(3tI易感人群 )(1tS 死亡 R治愈 )(2tIrC2图 1 改进后仓室模型框图4
10、.1.2 初步模型的建立根据图 1 所示模型框图建立结核病动力学微分方程组如下： RStItICrStINtdtIrS213 21211)()()()(及艾滋病人数 随季度增长的差分方程：)(2t )1()()212tt4.2 模型的求解（1）确定研究对象考虑到作为改革开放后首批发展迅速的沿海城市，广州经济发展十分迅速，但在结核病传染防治等方面也表现出了其非常脆弱的缺点：人口复杂、成分多且结核病病毒携带者病例年年增加。为此，选取具有代表性的广州市为研究对象。（2）求解微分、差分方程求解所建立的微分方程组其解析解是十分困难，我们根据搜集的数据，运用MATLAB 软件给出了微分方程组的数值解，并推
11、导了差分方程的解析解。（1）（2）（4）（3）（5）64.2.1 模型初值的确定查 01 年第一季度广州未接种疫苗人群数据知 ；983)1(I查 01 年第一季度广州登记肺结核病治愈人数 ；204.2.2 推导和计算（1）求出差分方程（5）的解析解 并112 )()()(1)( ttRSt代入微分方程组中。 根据参考资料【】查得 ；%80（2）利用 MATLAB 工具箱求解微分方程组（）得到 、 、 的数值解；1I23I4.3 模型结果分析（2001-2006 年）4.3.1 对结核病流行趋势评断指标的确定不失一般性，设定以下两指标对广州市结核病流行趋势进行分析和预测：（1） ：至该季度总染病
12、人数；2I（2） ：至该季度无法治愈人数；34.3.2 对 2001-2006 年广州市结核病流行趋势分析的计算结果经计算，得到结核病流行趋势按季度变化结果如表 1 所示：表 1 模型数值计算解年 季度 （人）2I（人）3I1 9834 10024 2 11816 12036 3 13706 13960 20014 14716 15013 1 15023 15370 2 15815 16223 3 17446 17927 20024 19065 19634 1 20223 20896 2 21377 22175 20033 22965 23915 74 24869 26000 1 26738
13、28087 2 28514 30124 3 30483 32406 20044 32778 35077 1 35247 37998 2 37777 41069 3 40442 44383 20054 43362 48083 1 46559 52215 2 49960 56739 3 53556 61682 20064 57426 67168 为了便于观察，我们将模型计算结果和根据国家公共信息网公布数据得到的流行趋势绘制成图 2 进行比对如下：图 2 计算数据与实际染病数对比图对比模型计算结果和实际结果可见，在第 1 至第 20 季度计算结果和实际结果基本吻合，第 20 季度以后误差开始显著增大
14、，即该模型可以较好的预测五年之内的结核病流行趋势。4.4 未来 5 年结核病传染预测4.4.1 预测时间长度的确定从 1953 年第一个国家五年计划至今，我国已经完成了 11 个五年计划。国家对8一个计划的周期进行了科学的选择，计划包含了医药卫生在内的许许多多重要领域，对我国的快速发展起到了巨大的促进作用。我们根据国家五年计划这一规划，将广州市的结核病传染中长期预测定位为 5 年，五年的医疗卫生计划既符合国家的发展计划，针对结核病传染也可以科学地采取相应的措施，对于政府的决策也用很大的指导意义。4.4.2 未来五年结核病传染预测利用广州 2008 年第一季度数据作为初值，带入初步模型中，计算结
15、果列于表 2：表 2 模型预测值（2008-2012）年 季度 21I3I1 45318 455082 48131 483363 51336 5146420084 55087 552121 59143 595352 63419 640013 68004 6849820094 72938 733811 78229 788512 83899 847643 89982 9098720104 96506 976041 103500 1047902 111010 1126103 119060 12101020114 127690 1299901 136940 1396602 146870 1501403
16、 157520 16146020124 168940 173660流行趋势预测曲线如图 3 所示9图 3 2008-2012 总染病人数预测根据初步模型计算结果和预测曲线可见，如果不采取有效措施，随着时间推移，结核病总染病人数将会迅速增长。4 .4 模型的评价在考虑艾滋病的影响和结核病存在耐药性两因素的基础上，给出了广州市结核病发病及传染的模型，并做出了从 2008 年到 2012 年的流行趋势预测，给出了结核病总染病人数 和无法治愈人数 随时间推移的发展趋势和具体计算结果，对政府及医21I3I疗机构的科学决策提供了依据。该模型没有考虑人口出生死亡对人口总数的影响、医疗水平提高对治愈人数的影响
17、等因素，且认为艾滋病按常指数增长，这些都影响了使用该模型中长期预测的准确程度。五、问题 2 预防和控制措施对结核病传播的影响结核病的传播对人们的身体健康、社会安定、经济的发展等方面有着深远的影响，采取科学的预防和控制措施对结核病的传播具有巨大的遏制作用。采取不同的预防和控制措施对结核病的传播具有不同程度的遏制，高效的措施将发挥其不可替代的独特作用，对控制结核病传染起到“事半功倍”的效果。10我国目前采取 DOTS 策略来预防和控制结核病的传播，主要包括以下措施：（1）政府对结核病规划的承诺；（2）通过痰涂片镜检发现病人；（3）在正确的管理下，给予标准的的短程化疗；（4）建立正规的药物供应系统；
18、（5）建立对规划执行的监督，评价系统；问题 1 中建立了分析预测结核病传播规律的初步模型，下面以初步模型为依据，对现有的预防和控制措施进行评价，并给出改进意见：5.1 措施评价指标的确立为便于评价所采取的预防控制措施对结核病流行情况的影响，仍将季度总染病人数 和季度无法治愈人数 作为评价指标，并与未采取措施的情况进行对比。)(21tI3I5.2 现有预防和控制措施对结核病传播的影响预防和控制措施对于初步模型的影响，主要体现在对初步模型中的参数、未及时治疗及治疗不当者所占总感染细菌者的比率 、单位季度 HIV 增)0(1SrC 长率 、HIV 感染者中感染肺结核细菌的比率 值的变化。故分别这四种
19、情况加以分2析：5.2.1 改变 的措施及其效果)0(1SrC经过调查，现阶段改变 的措施主要有以下 4 项：（1）通过痰涂片镜检发现和治疗染病人群以切断传染源以降低有效接触次数 ；C（2）加强对易感人群（艾滋、弱体制、特定年龄段）的保护、定期体检以降低易感人群受传染的传染率 ；r（3）加强宣传力度，督促养成良好卫生习惯以降低 和 ；Cr（4）建立特定人群（学生、工人、军人等）定时检疫和接种疫苗的制度以降低易感人群数量初值 ；)0(1S为研究 减小对结核病流行情况的影响效果，我们将该参数值缩小rC后代入初步模型中计算 2001 至 2006 年总染病人数，并将其与参数值缩小前的总%1011染病
20、人数对应的绘制在图 4 中：图 4 采取减小 措施前后总染病人数对比图)0(1SrC从图 4 中可见，与未采取措施的情况相比较，采取减少 的策略，可以)0(1SrC将总染病人数控制在约 左右，明显的抑制总染病人数急剧增加的趋势。403.15.2.2 改变 的措施及其效果改变 的措施主要如下：（1）对结核病人采取免费治疗的政策以降低 （2）对医院进行正规化管理和监督，提高治疗效果以降低 在验证 降低对结核病流行趋势影响时，我们将该值缩小 后代入初步模型中计 %1算 2001 至 2006 年总染病人数，并将其与参数值缩小前的总染病人数对应的绘制在图5 中：12图 5 采取降低 措施后总染病人数对
21、比图从图 5 中可以发现，5 年（20 个季度）后总染病人数由 左右降至 左右，4105410可见 值即使有微小的降低，也会较明显的抑制结核病流行的趋势。5.2.3 改变 的措施及其效果改变 的措施主要有：（1）免费进行艾滋病检测（2）对艾滋病进行免费治疗（3）加强对艾滋病危害的宣传以上三项措施都会降低 值，故我们将 减小 15%后代入初步模型中计算 2001 至2006 年总染病人数，并将其与参数值缩小前的总染病人数对应的绘制在图 6 中：图 6 采取降低 措施后总染病人数对比图13研究图 6 知，减小 可以在一定程度上减缓结核病流行趋势，但效果需要约 15 个季度后才较为明显，故对艾滋病的
22、控制应常抓不懈。5.2.4 改变 的措施及其效果2可以改变 的措施主要如下：（1）加强对艾滋病人的教育和指导（2）加强对艾滋病人的健康保护以上两项措施都会降低 值，故我们将 减小 70%后代入初步模型中计算 200122至 2006 年总染病人数，并将其与参数值缩小前的总染病人数对应的绘制在图 7 中：图 7 采取降低 措施后总染病人数对比图2由图 7 可见， 值即使减小较大，对结核病流行的趋势影响仍然较小。由式（）2知，这主要是由于艾滋病人数增长较快、基数较大抵消了 值减小的效果造成的。25.4 科学预防和控制措施的几点建议根据以上分析可以发现， 、 、 值的减小对抑制结核病流行趋势有)0(
23、1SrC较好的效果，为此，我们在现有预防和控制结核病措施的基础上提出了可以改进方案：（1）考虑到减小 对结核病流行趋势的抑制效果最为明显，故应着重建)0(1Sr14立和完善结核病监控体系的建立，切实实行好 5.2.1 所列 4 条措施；（2） 的减小程度主要依赖于政府政策对于传染性结核病人的扶住程度。考虑到现阶段政府对于结核病人的免费治疗还仅限于首次得病治疗和非耐药性结核病治疗，故建议政府加强政策支持，扩大对结核病治疗的范围，增加用于治疗结核病的经费，提高对耐药性结核病的治疗力度。这将对降低 值起到较好的效果；（3） 值对结核病流行趋势的抑制效果不能立刻体现，因此在采取各项预防和控制结核病措施
24、的同时还应长期坚持对艾滋病的检测和治疗；（4）在控制艾滋病增长的同时，扩大艾滋病免费治疗的范围，对艾滋病和结核病交叉感染的治疗实施政策优惠。六、问题 3 模型的改进6.1 初步建立的模型的不足初步模型将研究对象作为一个封闭的系统进行考虑，没有考虑到外来流动人口涌入的情况，而实际上流动人口尤其是农民工的涌入会对模型产生诸多影响：（1）流动人口无法接受规范治疗，从而增大染病总人数 值)(21tI根据有关资料，农民工已成为新发肺结核病的主要人群，因为农民工群体缺乏结核病防治知识，自我防范意识较弱，他们患病后难以纳入规范化的治疗过程，治愈率较低，一旦脱离治疗就会形成新的传染源在农民工内部和广州本地市民
25、之间传播。这部分流动人口涌入城市，将会对染病人群 的值产生较大影响；)(21tI（2）流动人口流动性强，会增大有效接触次数 的值C农民工群体一方面由于不是本地常住居民，具有极大的流动性；另一方面存在工作环境的通风、卫生条件普遍不高的缺点，这些必然带来有效接触次数 的增加。C而在初步模型中， 中没有考虑外来人口，且认为 是一个常数，如果)(21tI考虑外来流动人口，该模型就存在不足了。6.2 模型改进的几点建议根据以上分析，我们对初步模型进行如下改进：（1）在初步模型中的 项中引入农民工中的染病人数 ，其中 为每)(21tIpA季度净流入广州的外来人数的患病比率， 为每季度净流入广州的外来人数；
26、A15（2）根据参考资料，引入 将初步模型中的常数 改为与总人口相关变量 ，CNC其中 为人口增加影响系数。6.3 考虑人口流动的结核病传染模型6.3.1 人口流动结核病传染模型建立改进后的初步模型如下图 8 所示：图 8 改进后的初步模型框图根据改进后的初步模型框图 8，列写改进后的初步模型结核病动力学微分方程组如下： RStIpAtICrSttdtIIr213 21211)()()()(以及艾滋病人数 随季度增长的差分方程)(2t )1()( 21tS（6）（7）（8）（9）（10）)1(Cr166.3.2 模型求解6.3.2.1 初值的确定（1）查 01 年相关数据，在第一问的基础上考虑
27、外来人口染病数量，确定 01 年第一季度 ；其中 假设为 0.5。107599384)(pAI（2）查 01 年第一季度登记肺结核病治愈人数 ；3705)1(230)(pAI6.3.2.2 推导和计算（1）依初步模型的方法求出差分方程（10）的解析解并带入微分方程组；112 )()()()( ttRSt（2）利用 matlab 工具箱求解微分方程组得到 、 、 的数值解；1I23I6.3.3 模型结果分析6.3.3.1 对结核病传播趋势评断指标的确定依然采用初步模型中设定的两个指标 、 来分析广州市人口流动特点21I3对结核病传播和发作的影响。6.3.3.2 对 2001-2006 年广州市结
28、核病流行趋势分析的对比计算结果我们将计算得到的考虑人口流动时结核病流行趋势按季度变化结果和经初步模型计算的不考虑人口流动时的计算结果同时列于表 3，并将对比计算结果在坐标图 9 中表示出来：表 3 改进后初步模型的计算结果及比较I1+I2 I3年 季度 考虑流动染病人数 不考虑流动染病人数考虑流动染病人数不考虑流动染病人数1 12584 9834 10949 10024 2 17671 11816 16937 12036 3 22463 13706 22242 13960 20014 27164 14716 27205 15013 1 31308 15023 31012 15370 2 357
29、55 15815 35136 16223 3 41062 17446 40523 17927 20024 46858 19065 46521 19634 2003 1 52681 20223 52320 20896 172 58609 21377 58035 22175 3 65053 22965 64362 23915 4 72153 24869 71522 26000 1 79649 26738 79054 28087 2 87420 28514 86725 30124 3 95642 30483 94818 32406 20044 104490 32778 103610 35077 1
30、113940 35247 113060 37998 2 123930 37777 123000 41069 3 134460 40442 133420 44383 20054 145660 43362 144520 48083 1 157610 46559 156410 52215 2 170310 49960 169040 56739 3 183750 53556 182380 61682 20064 198020 57426 196550 67168 图 9 模型改进前后总染病数对比图6.3.3.3 由计算结果分析人口流动特点对结核病传播的影响由表 3 和图 9 可见，在考虑外来流动人口的
31、情况下，结核病流行程度比无外来流动人口的情况严重，如在第 20 个季度（第 5 年）时，前一种情况的总染病人数（145660 人）已是后一种情况（43362 人）的约 3.36 倍，且这种倍数关系随着时间推移有增大的趋势。根据模型 2 和计算数据说明，人口的流动会促进结核病的传播。6.4 人口流动特点对结核病预防控制的影响18在问题二中，我们分析并证明了现有的结核病预防控制措施的效果主要通过减小、 、 来抑制疾病的流行，而根据式 1 和式 6 可知人口流动特点的引入，)0(1SrC会增加 和 的值，从而削弱正常的预防控制措施。 6.5 改进后模型的评价改进后的模型用有效接触率 和总染病人数 的
32、变化来表示人口流动特C)(21tI点对模型的影响，计算了人口流动对结核病传播的影响，并分析了其对现有预防控制措施的影响，较问题一中建立的初步模型有更好的适应性，为进一步制定预防和控制措施提供较好的参考。给报刊的短文科学预防控制结核病传染的重要性结核病传染模型启示 自人类出现在地球上，就一直受病毒的侵扰，传染病肆虐成狂。结核病就是一个非常典型的例子，它是一种古老的经呼吸道传播的慢性传染病，在历史上流行猖獗，被称为白色瘟疫。尽管上世纪 50 年代有效的抗结核药物问世以后，形势发生了根本的改变，但是进入八十年代以后，由于各国政府的忽视、移民、多耐药结核病和艾滋病等因素的影响，在全球，特别是在发达国家
33、又重新出现了上升的趋势，为此，2010 年我国再一次进行了结核病普查。针对我国 2000 年结核病流行病调查结果分析：我国结核病疫情还相当严重，部分地区有蔓延趋势， 表现为结核菌感染人数多、现患肺结核病人多、结核病死亡人数多、农村结核病人多和传染性肺结核病疫情居高不下等。它对我国的经济、社会安定等方面起到了很大的冲击作用，严重影响了我国的国民身体健康，结核病已成为了我国人口非自然死亡的首要因素！传染病对于人们来说一直是一个难题。对于一种传染病，医学工作者会不断深入的探索与研究，然后研制出相关的疫苗或是制定相关的措施，最大程度地减小疾病的危害。当然，在研究传染病本身的同时，我们还迫切需要了解传染
34、病的传染形式、与其相关的措施对病况的影响及未来的传染趋势进行预测等等，以防它对国家人民的健19康、社会经济的发展产生更加深远而不可逆的影响。而这些就要依靠数学工具进行处理解决：数学建模在基于现实情况的基础之上对问题就做出了详细的解答，指导人们实践，对传染病传染做出了科学的解释与预测，有效遏制了传染病的进一步发展。我国结核病发展现状令人担忧，防治任务艰巨，结核病的预防控制压力空前之大。从相关的数据上可以很明显的看出来近三十年来我国的结核病患者是逐年增高，由于部分地方政府承诺落实不足、结核病防治经费投入仍缺口较大；流动人口的不断增多、青壮年体质的下降、结核菌/艾滋病毒双重感染和耐药结核病的挑战等因
35、素的影响，这种趋势也是相当明显的，所以，科学的分析与预测与防控就显得更为重要！假若没有了科学的理解，对重要的措施如外来人口的医学检测、当地医疗水平没有进行加强等，没有对传染趋势做出科学的预测，那这些因素将直接影响我国的结核病传染形势，我们可以想象一下现在结核病会对我们的国家社会造成多么大的影响！结核病的预防控制迫在眉睫！对于结核病，有几个主要因素造成了传染趋势上升？又是什么因素造成的结果最严重？政府和医疗机构应各采取多大的力度来处理结核病？这些问题的结果对于结核病的防治是相当重要的也是非常必要的。如果能够科学的回答这些问题，那么将对我国今后的防治结核病工作产生巨大的推进作用！由此可见像我国这样
36、的结核病传染大国，科学的预防控制非常重要。只有我们在平时就做好防范，对于没有根除的疫情，当其重新抬头时，我们就可以从容面对，战胜疫情，并最终将现在还没有根除的传染病根除，为我们社会、国家、世界做出应有的贡献！结核病给我们的正常生产生活造成了很大的影响，2010 年全国结核病普查再一次给我们敲响了警钟。我们应不断加强结核病的预防与控制，并以此来指导生产实践，因为它关系着我们每个人的健康！七、参考文献1 王正林、刘明， 精通 MATLAB7 ，北京：电子工业出版社，2008。2 余洋、胡志兴， 一类具有常数移民且带隔离项传染病模型的分析 ，石家庄学院学报，第 8 卷第 6 期：43-48，2006
37、。3 姜启源、谢金星、叶俊， 数学模型（第三版） ，北京：高等教育出版社，2003。4 杨秀香， 接触率依赖于总人数的传染病动力学模型 ，渭南师范学院学报，第 22卷第 2 期：17-19，2007。5 公共卫生数中心，20012008 年广州市结核病监测数据，http:/ 年。6 新浪网，我国结核病有六多， ，2010。7 赵介民，肺结核的基本知识，北京：科学技术出版社，2009附录1.1 解初步模型中的微分方程function dx=Q1equation(t,x)dx=1;1;dx(1)=1.01*x(1)-x(1)+1.01*x(2);dx(2)=0.06*x(2)-x(2)+0.06*
38、x(1);21endt=1:0.1:24;n=length(t);x0=9834 2330;opt1=odeset(RelTol,1e100,AbsTol,1e100);%解微分方程组在此t,x=ode23(Q1equation,t,x0,opt1);plot(t,x(:,1),.,t,x(:,2),-.,t,x(:,1)+x(:,2),-g)%绘制总染病人数图hold ont,x1.2 由已知染病人数数据绘图deltaren=2330 1034 1394 1302 1336 1274 1526 3180 13122009 1554 2087 1875 2091 3092 1962 1638
39、1535 17471738 1797 3085 1741 1608ren=ren(1)=9834.6for k=2:24ren(k)=ren(k-1)+deltaren(k)plot(k,ren(k),o)hold onend221.3 预测 2008-2012 年染病总人数t=1:0.1:20;n=length(t);x0=45318 2023;%取初值opt1=odeset(RelTol,1e100,AbsTol,1e100);%解微分方程组在此t,x=ode23(Q1equation,t,x0,opt1);plot(t,x(:,1),.,t,x(:,2),-,t,x(:,1)+x(:,2
40、),-)hold ont,x1.4 解改进后初步模型微分方程function dx=improveofQ3equation(t,x)s=x(1)+x(2)+2750;dx=1;1;dx(1)=1.03*s-x(1);dx(2)=0.03*s-x(2);endt=1:0.1:20;n=length(t);23x0=10759 3705;opt1=odeset(RelTol,1e100,AbsTol,1e100);%解微分方程组在此t,x=ode23(improveofQ3equation,t,x0,opt1);plot(t,x(:,1),.,t,x(:,2),-.,t,x(:,1)+x(:,2)
41、+2750,-)hold ont,x1.5 参数值及其来源：未及时治疗及治疗不当者所占总感染细菌者的比率查耐药性肺结核占感染者平均比率 0.28（广州日报 07 年 3 月 22 日）：所研究地区总人数7125979（2001 年广州户籍人口）N：极易感人群数（主要是 HIV）查 2001 年 HIV 的人数 6902S：易感人群数（未接种疫苗）查广州已接种疫苗人群的覆盖率十岁以下儿1童人数或比率 0.213（深圳特区报 06.03.28）：HIV 感染者中感染肺结核细菌的比率查平均值 0.8(百科),0.128（南方日报207.03）：每季度肺结核的死亡人数查平均值 26 人（广州日报 07 年 3 月 22 日）R：单位季度 HIV 增长率查平均值 0.8（百度百科）：每季度净流入广州的外来人数查平均值 299.76（2001 中国城市发展报告）A：每季度净流入广州的外来人数的患病比率查广州外来人口肺结核人数均值p275（网易新闻 09.03.23）
