2014年第三届数学建模国际赛小美赛C题埃博拉病毒:为什么不能遏制传播解题全过程文档及程序

2014年第三届数学建模国际赛小美赛

C题 埃博拉病毒:为什么不能遏制传播？

原题再现：

  致命的埃博拉出血热最早发现于1976年，自从理查德·普雷斯顿的《热区》出版以来，它就一直困扰着公众20年。然而，在这段时间里，没有一种药物被批准用于治疗这种疾病。据世界卫生组织称，自2013年12月以来，西非持续爆发的埃博拉病毒已在几内亚、塞拉利昂和利比里亚感染至少567人，其中350人死亡。这次爆发似乎是历史上规模最大的一次，超过了2000年乌干达埃博拉爆发时发生的425例病例。埃博拉病毒感染者只接受旨在支持患者的一般治疗。根据疾病控制和预防中心的说法，他们可能需要补充液体（埃博拉患者经常脱水），或者进行旨在维持血压和氧气水平的治疗，如果出现感染，还可以进行治疗。（向两名美国患者提供的实验药物已经用完）埃博拉病毒治疗的缺乏令人不安。但是，考虑到药物开发的资金来源，这也是可以预测的。
  当制药公司决定将研发资金投向何处时，他们自然会评估候选药物的潜在市场。这意味着他们有动机针对那些影响更富有的人（首先是发达国家的人）的疾病，这些人有能力支付很多费用。他们有动机去制造许多人会服用的药物。他们有一个动机，使药物，人们将定期服用很长一段时间的药物，如他汀类药物。
  这个系统在让西方人得到他们想要的药物（尽管通常价格很高）方面做得很合理。但这也导致了对某些疾病和某些药物种类的巨大投资不足。主要影响穷国穷人的疾病并不是研究的重点，因为这些市场不太可能提供可观的回报。因此，像疟疾和肺结核这样每年夺去200万人生命的疾病，在制药公司受到的关注比高胆固醇要少。世界卫生组织称之为“被忽视的热带疾病”，如查加斯病和登革热；它们影响到10亿多人，每年造成50万人死亡。一项研究发现，在1975年至2004年间上市的1500多种药物中，只有10种是针对这些疾病的。当一种疾病的受害者既穷又不多的时候，这是双重打击。从这两个方面来看，一种治疗埃博拉病毒的药物看起来都是一项糟糕的投资：到目前为止，这种疾病只出现在贫穷国家，影响到的人数相对较少。
  然而，这一体系损害的不仅仅是发展中国家。近年来，耐药微生物的增多使得我们使用的抗生素的效果降低，增加了传染病失控的风险。卫生官员同意，西方人需要的是我们可以保留的新药，以防常规抗生素无法控制的疾病爆发。然而，在过去的三十年里，新抗生素的供应已经放缓到涓涓细流。波士顿大学卫生法项目的联合主任凯文·奥特森（Kevin Outterson）是美国疾病预防控制中心（C.D.C.）抗生素耐药性工作组的创始成员，他告诉我：“抗生素耐药性确实有可能使我们生活方式的方方面面发生变化。因此，我们需要推动这条管道。”
  问题还是在于商业模式。如果一家制药公司真的发明了一种功能强大的新抗生素，我们就不会希望它被广泛使用，因为它的目标是延缓耐药性。奥特森说：“公共卫生官员会尽可能适当地限制这种药物的销售。”；投资前景不佳。
  那么，我们如何才能获得我们所需要的药物，而不神奇地改变开发药物的行业呢？关键是奖励那些创造了大量公共卫生福利的公司。最简单的方法就是为新药提供奖励。Outterson描述了一个场景：“政府将向公司支付一笔或一系列款项，作为交换，公司将放弃销售该产品的权利。”制药公司将获得报酬，并将避免尝试推出新产品的所有费用（无论如何，你不希望使用最后手段的抗生素）。社会将得到一种新药，公共卫生官员将能够控制它的推广和使用方式。
  任务#1：选择您认为在非洲埃博拉疫情中最严重的国家。建立一个模型，在没有任何额外药物的情况下，估算2006年至2025年该国埃博拉感染人数的预期变化率。充分解释您的模型和模型的基础假设。此外，估计我们应对危机的成本。
  任务#2：建立一个数学模型，确定在此期间该地区是否有任何新疾病需要开发药物。如果是，估计药物开发过程的截止日期和产量。
  任务#3：制定制药企业“奖励”政策，支持药品研发工作。您制定的方法应利用公共卫生福利和美国预算。
  任务#4：多公司联合开展药物研究，可有效降低成本和风险。如何选择一些企业建立合作关系是一个棘手的问题。您的团队可以帮助政府做出选择吗？

整体求解过程概述(摘要)

  针对这一问题，我们从各个方面对埃博拉病毒进行了研究。我们查阅了互联网上有关埃博拉病毒在非洲爆发的数据。比较不同国家的病例、死亡率等影响最严重的因素。找到问题的内部信息和特征对我们来说不是很明显，也就是说，有一部分是未知的，然后找到不同的数学模型来求解，利用收集到的数据绘制出速率曲线，并预测其未来的变化。
  围绕任务1，我们从不同方面找到了美国成本的财务状况。预测埃博拉病毒治疗药物的经济发展。并预测美国将花费多少来应对这场危机。收集有关非洲病毒的信息。联系经济、病毒等诸多因素，寻找合适的国家。利用灰色预测模型和GM（1,1）模型在MATLAB中进行编程，绘制各因素的曲线。
  任务#2是多项式回归曲线的拟合，我们用Logistic回归模型进行求解，将不同的数据进行拟合。通过这两个模型，我们预测了药物开发的可能性。最后对埃博拉病毒检测结果进行了分析。研究药物研究的影响因素。发现主要影响来源是病毒的发育速度、死亡率、变化的脱色率。将三个因素拟合到变化曲线中。
  针对任务#3，建立了“奖励”策略模型，描述了“奖励”策略各方面的客观性。在此基础上，提出了创新的“奖励”政策。激励公司创造更大的社会效益。完善“奖励”政策。
  在任务4中，按照各公司的方案处理药品开发风险和不同成本。通过任务分配模型，降低药物开发的成本和风险。优化作业，提高工作效率。

模型假设：

  (1） 在没有其他新的病毒感染的情况下，我们只讨论埃博拉病毒的发展。
  （2） 摆脱环境和温度的巨大变化。
  （3） 确保病毒不会有很多变态。然后，它保持其主要的增长曲线。
  （4） 该新药尚未开发，因此不会影响埃博拉病毒的传播速度。
  （5） 不同的病毒有不同的发展速度。在我们的模型中，我们关注的是会被感染的病毒。
  （6） 我们的模型不考虑非感染病毒。当我们确认研究的合理性时，我们可以在合理的时间内完成。
  （7） 我们研究的人口问题是根据人口的平均密度。

问题分析：

  问题1分析
  我们需要选择一个非洲国家作为埃博拉病毒感染的主要研究目标。根据世界卫生组织的数据，我们关注乌干达。根据2000年、2007年、2011年、2012年和2013年埃博拉病毒感染的现有数据，通过线性插值得到2000年至2014年的相关数据。然后利用灰色预测方法对2014-2025年的感染风险人数进行估计，最后计算出感染人数的变化率。对于美国对这种病毒的财政援助，我们的目标是医药投资的成本。根据相关文献提供的药品不同阶段的机会成本，我们以现金和资本两种形式评估支出。结合美国的贴现率和成本消耗的增长，我们可以计算出2010年每种药物的平均成本。

  问题2分析
  这个问题要求我们确定2006年至2025年乌干达是否需要新药。我们选择该区域更重要的疾病类型：结核病、艾滋病和疟疾。如果方差检验有显著性差异，则新药的需求是明显的。在确定新药需求后，我们提取每种疾病每年的异常值感染数。药品的总产量是根据一定时期内感染的确切人数来估计的。

  问题3分析
  这个问题要求我们设计一个方案，在药物开发过程中为制药公司提供资金支持。该方案旨在降低药品价格，压缩药品销售，使非洲较贫穷的国家能够负担费用，使药品不会因广泛使用而降低药效。我们将对制药企业的补贴分为两个阶段：药品研发阶段和销售阶段。政府补贴来自政府自己的预算和公共卫生福利。为了便于计算，我们将公共卫生福利资本化，从而获得最佳的援助方案。

  问题4的分析
  第四个问题要求政府设计一个选择公司合作药品的方案，以降低成本和风险。为了有效发挥不同企业的优势，发挥其在医药行业的优势，根据公司的相关数据，我们定义了四个指标来评价其药物研发能力，根据各公司的各种参数值，选择各公司的最大优势进行混搭。最后，对于可持续性的考虑，我们考虑了公司数量增加所带来的影响。

模型的建立与求解整体论文缩略图

全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可

部分程序代码：(代码和文档not free)

clc,clear 
Number=input('please enter the codes of problem(1 min risks 2lowest cost）Number=：');
if Number==1
 F=[5 4 5 1 2 3 2 4 5 5 1 5 1 3 1 2 3 2 4 2 1 3 4 2 3];
elseif Number==2
 F=[20 19 20 28 17 18 24 27 20 20 26 16 15 18 15 17 20 24 19 16 15 18 21 17 21];
end
A=[];B=[];
Aeq=[ones(1,5) zeros(1,20);zeros(1,5) ones(1,5) zeros(1,15);zeros(1,10) ones(1,5) 
zeros(1,10);...
 zeros(1,15) ones(1,5) zeros(1,5);zeros(1,20) ones(1,5);...
 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4);...
 zeros(1,1) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,3);...
 zeros(1,2) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,2);...
 zeros(1,3) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,1);...
 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,4) 1 zeros(1,0)];
Beq=ones(1,10);
if Number==1
 [x,min_fval]=bintprog(F,A,B,Aeq,Beq);
 min_fval2=min_fval/5;
 fprintf('(1)the match function of company：\n')

fprintf(' (A1) the match function of company A（1 replace match 0 replace not match）：
\n')
 disp(' A-1 A-2 A-3 A-4 A-5'),disp(x(1:5)')
 fprintf(' (B1) the match function of company B（1 replace match 0 replace not match）：
\n')
 disp(' B-1 B-2 B-3 B-4 B-5'),disp(x(6:10)')
 fprintf(' (C1) the match function of company C（1 replace match 0 replace not match）：
\n')
 disp(' C-1 C-2 C-3 C-4 C-5'),disp(x(11:15)')
 fprintf(' (D1) the match function of company D（1 replace match 0 replace not match）：
\n')
 disp(' D-1 D-2 D-3 D-4 D-5'),disp(x(16:20)')
 fprintf(' (E1) the match function of company E（1 replace match 0 replace not match）：
\n')
 disp(' E-1 E-2 E-3 E-4 E-5'),disp(x(21:25)')
 fprintf('min risks：'),disp(min_fval2)
elseif Number==2
 [x,min_fval]=bintprog(F,A,B,Aeq,Beq);
 min_fval2=min_fval/5;
 fprintf('(2)the match function of company：\n')
 fprintf(' (A2) the match function of company A（1 replace match 0 replace not match）：
\n')
 disp(' A-1 A-2 A-3 A-4 A-5'),disp(x(1:5)')
 fprintf(' (B2) the match function of company B（1 replace match 0 replace not match）：
\n')
 disp(' B-1 B-2 B-3 B-4 B-5'),disp(x(6:10)')
 fprintf(' (C2) the match function of company C（1 replace match 0 replace not match）：
\n')
 disp(' C-1 C-2 C-3 C-4 C-5'),disp(x(11:15)')
 fprintf(' (D2) the match function of company D（1 replace match 0 replace not match）：
\n')
 disp(' D-1 D-2 D-3 D-4 D-5'),disp(x(16:20)')
 fprintf(' (E2) the match function of company E（1 replace match 0 replace not match）：
\n')
 disp(' E-1 E-2 E-3 E-4 E-5'),disp(x(21:25)')
 fprintf('lowest costs：'),disp(min_fval2)
end
%The result of these as follows：
Please enter the codes of problem（1min risks 2 lowest cost）Number=：1
(1）the match function of company：
(A1) the match function of company A（1 replace match 0 replace not match）：
 A-1 A-2 A-3 A-4 A-5
0 0 0 1 0
(B1) the match function of company B（1 replace match 0 replace not match）：
 B-1 B-2 B-3 B-4 B-5
 0 1 0 0 0
(C1) the match function of company C（1 replace match 0 replace not match）：
 C-1 C-2 C-3 C-4 C-5
 0 0 0 0 1
(D1) the match function of company D（1 replace match 0 replace not match）：
 D-1 D-2 D-3 D-4 D-5
 0 0 1 0 0
(E1) the match function of company E（1 replace match 0 replace not match）：
 E-1 E-2 E-3 E-4 E-5
 1 0 0 0 0
min risks： 1.4000
Please enter the codes of problem（1min risks 2 lowest cost）Number=：2
(1）the match function of company：
(A2) the match function of company A（1 replace match 0 replace not match）：
 A-1 A-2 A-3 A-4 A-5
 0 1 0 0 0
(B2) the match function of company B（1 replace match 0 replace not match）：
 B-1 B-2 B-3 B-4 B-5
 0 0 0 1 0
(C2) the match function of company C（1 replace match 0 replace not match）：
 C-1 C-2 C-3 C-4 C-5
 0 0 1 0 0
(D2) the match function of company D（1 replace match 0 replace not match）：
 D-1 D-2 D-3 D-4 D-5
 0 0 0 0 1
(E2) the match function of company E（1 replace match 0 replace not match）：
 E-1 E-2 E-3 E-4 E-5
 1 0 0 0 0
lowest costs： 17

全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可