实验一 EXCEL规划求解

实验一EXCEL规划求解实验目的和任务训练建模能力.应用EXCEL建模及求解的方法应用；1-3通过实验进一步掌握运筹学有关方法原理、求解过程，提高学生分析问题和解决问题能力。实验仪器、设备及材料计算机、Excel实验内容JtM*Jl一个人口不到15万人的小**根章妁级英叶植情K..务N住花.*(A)、(B)W(C)过三全我的清工5上011鼐将JL*达，与收的再由a柱内上分送到金■的8)基**. 珞情丸、冬雄以乐舄(单位：100m)及各收的晨、赛事小①.….®的具体位JLXS3.2.按翕年情义.A、B、C三个4f.♦天收的量分刷为200、170160(单位：100k«>9各款审场的■•天才求量及女生柒&缰候时♦皋的擅大(尢/100kg)JL4.3-48.设从收成点重分装・塔或某斗1速食用力1>t/(100kg*100m).・3-4S»««每天■学(100k«)*9♦央(%/1009①7510⑦608(3)903④7010⑤10010SS8⑦903808(a)为《率设计一个从各收0点£方最事*的;t*轶墟才童.使用于建东调述及衣期的姓金君失为最小,(b)冬蝮之各基卡场缰段量一律不曜过禽生量的20%.★新谩计文晟供膜方金：(c)为遇及**后艮的■.菜供&・比明♦的e导版则瑞加・.冢“状面机，二问Jt户衲天原分则向A、B、C三个采购上备供及多少最缓济会履建模解：该问题的目标是求运费最小，则设如表：收购①②③④⑤⑥⑦⑧量EAA1A2A3A4A5A6A7A8200BB1B2B3B4B5B6B7B8170CC1C2C3C4C5C6C7C8160每天需求75608070100559080610短缺1512161420111816122必需量6048645680447264488£530设增产的蔬菜每天应分别向A、B、C三个采购点各供应M、N、。最经济合理。短缺量为xl,x2,x3,x4,x5,x6,x7,x8omi为每天的需求量，ni为每天的必需量。目标函数为：MaxZ=4A1+8A2+8A3+19A4+11A5+6A6+22A7+20A8+14B1+7B2+7B3+16B4+12B5+16B6+23B7+17B8+20C1+19C2+11C3+14C4+6C5+15C64-5C74-10C8+10X1+8X24-5X3+10X4+10X54-8X6+5X74-8X8三个小题目标函数一样，但约束条件分别为：问题(一):⑴约束条件：Ai+Bi+Ci<mi(i=l,2,…，8;mi=75,60,80,70,l00,55,90,80)Ai+Bi+Ci+xi=mi(i=l,2,…，8^ni=75,60,80,70,100,55,90,80)Al+A2+A3+A4+A5+A6+A7+A8=200Bl+B2+B3+B4+B5+B6+B7+B8=170Cl+C2+C3+C4+C5+C6+C7+C8=160问题(-)：(2)约束条件：Ai+Bi+Bi<mi(i=l,2,…，8;mi=75,60,80,70,l00,55,90,80)Ai+Bi+Ci+xi=mi(i=12,…,8;mi=75,60,80,70,l00,55,90,80)Al+A2+A3+A4+A5+A6+A7+A8=200Bl+B2+B3+B4+B5+B6+B7+B8=170Cl+C2+C3+C4+C5+C6+C7+C8=160Ai+Bi+Bi>ni(i=1,2,...,8;ni=60,48,64,56,80,44,72,64)问题(三)：(3)约束条件：Ai+Bi+Bi<mi(i=l,2,…,8;mi=75,60,80,70,100,55,90,80)Al+A2+A3+A4+A5+A6+A7+A8・A=200Bl+B2+B3+B4+B5+B6+B7+B8-B=170C1+C2+C3+C4+C5+C6+C7+C8・C=160(所有变量均大于零)求解过程

求解的方法步骤和案例3.1的方法步骤一样，简要抓图如下问题(―)：licrosoftExcel11.0敏感性报告工作表[复件菜篮子l.*ls]Sheetl报告的建立：2007-12-0415:59:34可变单元格 终通社目标武允许的允许的TOC\o"1-5"\h\z单元格名字值成本系数增量减量$A$3 M 75 0 4 n 1的30$B$3 A2 0 0 8 1E+30 0$C$3 A3 0 0 8 1E+30 0$D$3 A4 0 2 19 1E+30 2$E$3 A5 70 0 Ho 9炉$ A； 55 0 6 9 1E+30%$3 A7 0 12 22 1E+30 12Hr A3 0 5 IR” F|I$3Bl I'11 14!=+:!? 11目标式单元格名字值$AG$1 4610目标式

结果75 4610<J 461。0 4610。 461。70 的历55 46100 4610目标式单元格名字值$AG$1 4610目标式

结果75 4610<J 461。0 4610。 461。70 的历55 46100 46100 46160 461060 4610±M 目标式极网 结果TOC\o"1-5"\h\z75 46100 _46100 4610―0—她g70 4；1「55 4610r>4”。。的00的060 4610问题(-):

工作表［复件菜篮子l.xls］（2）报告的建立：2007-12-0416:20:12可变单元格终递减目标式允许的允许的单元格名字值成本系数增量减量$A$3Al750471E+30相$3A2100802$C$3A30081E+300$D$3A402191E+302$E$3A56001121$F$3A6550631E+30和$3A7012221E+3012$H$3A805201E+305$I$3Bl011141E+3011licrosoftExcel11.0极限值ft（告工作我［艮件1报告的建立：2007-12F416:20:14单元格目标式名字佰$AG$l4806单元馅费量各牛侑下眼做吊目标式结果目标式结果加375754806754806加3A210104806104806|C|3A30T.68434E-144806-5.68434E-144806|D13M0-5.68434E-144806-5.68434B-1448061E13A560604妫604806林$3A555554806554806$G$3A70-5.68434E-144806-5.68434E-144806呻A8Q-5.68434E-144806-5.68434E-M4806IH3BI004806-Q4806问题（三）:licrosoftExcel11.0敏感性报告工作表【复件艮件菜筌子（2）收告的建立：2007-12-1317:45:39可变单元格终且出L各生值一SM3 A1 75A2 可变单元格终且出L各生值一SM3 A1 75A2 0$03 微 加口标式

累数,8允许的

增量允许的11 1E+301E<j 001.999999997$g3 A4 0 2.000000004 1； 1E+30 2.OOOOQOOO43BI3 A5 30 0 13 L999999WT 21XX>000003“13M56Q$11 1龄30103 W 0 12 22 1E+30 12$H13 AS 0 4.999999%519.999^9999 1E+30 4.91SI$3Bl0n 1E+30n匚$3 E二 60 0 7 0 1F-CrlicrosottExcel11.UIftMUIAtWi［作表［fcftftft菜IS子3」】・】极用侦报告I程分的■立：200T12-1317:45:40尊元餐川$1名字・4770里冗余名名—(A卜限ffFF日磨式上用砧用口标式结果75754770754770皿3A202.36966E-1147702.36966E-114770$53A340404770404770晔A400477004770IBI3A53030477。30477。"13然55554770554770IGI3A7__CL .0477004770$咫3A800477004770$【13B1一&JL4770JL4770IJS3B260的4TT0的4770结果分析问题（一）：调运及预期短缺损失最小为4610元。在这情况下，由A点供应3、4、7、8市场的供应两分别为8、56、72、64,由B点供应1、5、6三市场的供应量分别为60、80、30,由C点供应2、3、6市场的供应量分别为90、56、14,A、B、C点的蔬菜皆运出去，1市场短缺15,3市场短缺16,4市场短缺14,5市场短缺20,6市场短缺11,7市场短缺18,8市场短缺16,所以为使调运及预期短缺损失最小2市场多出30.问题（二）：各市场短缺量一律不小于需求量20%时，调运及预期短缺损失最小为4806元。在这样的情况下，由A点供应1、2、5、6市场的量分别为75、10、60、55,由B点供应2、3、4市场的量分别为50、64、56,由C点供应5、7、8市场的量分别为24、72、64,3市场缺量为16,4市场短缺14,5市场缺量为16,7市场短缺18,8市场短缺16.问题（三）：保证蔬菜供应情况下增产蔬菜向C点供应80即可满足要求，此时最小费用为4770元。由A点向1、3、5、6市场的供应量分别为75、40、30、55,由B点向2、3、4市场的供应量分别为60、40、70,由C点向5、7、8市场的供应量分别为70、90、80.实验二：最小生成树与最短路一、实验目的：通过题目编写最小生成树与最短路的Matlab程序，掌握最小生成树的K算法和P算法，以及最短路的D算法和F算法。二、实验内容：(1)最小生成树问题：用Matlab编程分别求解该题的K算法和P算法程序。⑵ 最短路问题：分别用D算法和F算法求解下题的最短路。(1)最小生成树K算法：Kruskal算法适用于一般非负权的无向网络求最小树。编程思想：输入加权连通图的边权矩阵g(i,力]““3，定点数n。1)整理边权矩阵将仍(1，J)]mx3按第三行由小到大的顺序重新排列，得到新的边权矩阵j)]mx3;2)初始化j- — —0;对所有3)更新T,c,t(i),若1(5(1,/))=7(8(2")),则转4);否则，若则对所有i,若=max{r(B(l,j)),/(B(2,；))},则r(i)<-min{r(B(l,J)),r(B(2,j))};4)若卜=①1或)=11,则终止，输出T,c；否则返回3)。程序代码：%k算法b=[l1122334;23435455;17364852];%边权矩阵[B,i]=sortrows(b\3);B=B'; %按边权由小到大重新排列矩阵b的列m=size(b,2);n=5;t=1:n;k=O;T=[];c=0;fori=l:mift(B(13))-=t(B(2,i))%判断第i条边是否与树中的边形成圈k=k+l;T(k,l:2)=B(l:23);c=c+B(3,i)；tmin=min(t(B(l3)),t(B(2,i)));tmax=max(t(B(lj)),t(B(2j)));forj=l:nift(j)==tmaxt(j)=tmin;endendendifk==n-lbreak;endendT,c运行结果：T=TOC\o"1-5"\h\z1 24 545c=11因此，图一的最小生成树的边集合为{(1,2),(4,5),(1,4),(3,5)},总费用为11。得到的最小生成树如下图三：① 1 3 ③④ 2 @— (图三)P算法：编程思想：输入加权图的带权邻接矩阵j)]nxno1）建立初始候选边表，T<-（/>■2）从候选边中选出最短边（u,v）,T<—TU（m,v）;3)调整候选边集；4)重复2),3)直到T含有n-1条边。程序代码：a=[0173inf;l06inf4;76085;3inf802ynf4520];T=[];c=0;v=1;n=5;sb=2:n;forj=2:nb(l,j-l)=l;b(2,)-l)=j;b(3,j-l)=a(l,j)；endwhilesizc(T,2)<n-l[m,i]=min(b(3,:));T(:,sizeCT^)+l)=b(:,i)；c=c4-b(3,i)；v=b(2J);temp=find(sb==b(2,i))；sb(tcmp)=[];b(:j)=0；forj=l:length(sb)d=a(v,b(2,j));ifd<b(3,j)b(l,j)=v;b(3,j)=d;endendendT,c运行结果：TOC\o"1-5"\h\zT=114 52 4 5 313 2 5c=11因此,图一的最小生成树的边集合为{(1,2),(1,4),(4,5),(5,3)},总费用为11。得到的最小生成树同图三。(2)最短路D算法：Dijkstra算法适用于非负权有向网络求最短有向路，当与•任A时，W--=ooo程序代码：%SP_d.mn=input(*Entertheverticesnumberofthegraph:n=*)W=inputfEntertheweightadjacentmatrixofthegraph:[W(l,l),..,W(l,n);..;VC'(n,l),..,\K'(n,n)]=l)%stcpO:preparationT=l:n;T(l)=0;P=zcros(l,n);i=l；U=oncs(l,n)*inf;U(l)=0;%mainstepsforr=l:n%thestop-conditionofthealgorithmisS=0.%stcp3forj=l:nifU(j)>W(i,j)+U(i)U(j)=W(i,j)+U(i);P(j)=i;endend%U(j)=min{Ug),W(k,j)}%step1min二inf;forj=l:nif(T(j)>0)&(U(j)<=min)min=U(j);i=j;endend%U(i)=min{U(j)}%stcp2T(i)=0;%choosev(i)endinputfthepre-venixlabellingoftheshortestpath(from1totheother)is:)P运行结果：Entertheverticesnumberofthegraph:n=6n=6Entertheweightadjacentmatrixofthegraph:p^'(l,l),..,W(l,n);..;W(n,l),.-,W(n5n)]=[inf,l,3,inf,inf,inf;inf,inf,inf,4,inf,inf;inf,2,inf,inf,5,inf;inf,inf,2,inf,inf,3;inf,inf,inf,5,inf,2;inf,inf,inf,inf,inf,inf)W=Inf13InfInfInfInfInfInf4InfInfInf2InfInf5InfInfInf2InfInf3InfInfInf5Inf2InfInfInfInfInfInfthepre-vertixlabellingoftheshortestpath(from1totheother)is:P=TOC\o"1-5"\h\z0 112 3 4因此，由MATLAB程序的运行结果可知网络G中自点1到其他各点的最短路为：②一4 ⑷>二 J①/⑥、-5 L (图四)F算法：编程思想：d(i,J):喈；path(ij):对应于畤)的路径上i的后继结点，最终的取值为i到j的最短路径上i的后继点。输入带权邻接矩阵A=[a(i,，)]“*"o1)赋初值对所有i,j,d(/,j)=a(i,j);当a(i,j)=8时，path(i,j)=0,否则path(i,j)=j;k=l;2)更新d(i,j),path。,j)对所有ij,若d(i,k)+d(k,j)Nd(i,j),则转3);否则，d(i,j)=d(i,k)+d(k,j),path®j)=path(i,k),k=k+l,继续执行3);3)重复2)直到k=n+1。程序代码：%Floydfunction[D,path]=flodl(a)n=size(a,l);D=a;path=zeros(n,n);fori=l:nforj=l:nifD(i,j)〜=infpath(i,j)=j;endendendfork=l:nfori=l:nforj=l:nifD(i,k)+D(k,j)<D(i,j)D(i,j)=D(i,k)+D(k,j);path(i,j)=path(i,k);endend

endend运行结果：a=[0,l,3,inf,inf,inf;inf,0,inf,4,inf,inf;inf,2,0,inf,5,inf;inf,inf,2,0,inf,3；inf,inf,inf,5,0,2;inf,inf,inf,inf,inf,0];[D,path]=flodl(a)D=013588Inf064117Inf20657Inf42073Inf97502InfInfInfInfInf0path=123232024444023255033436044456000006因此，由MATLAB程序的运行结果可知网络G中任意两点的最短路径及其长度。其最短路见图四。四、实验心懵：通过这次实验，复习并巩固了最小生成树和最短路的算法，且熟悉了MATLAB的编程思想，为用MATLAB软件解决实际问题打好基础。实验三数学建模实验目的学习运用运筹学的方法解决实际问题。实验要求.任选一个问题，建立数学模型。.用LIND。软件求解模型。.对论文格式的统一要求如下：.摘要、2.问题的提出3.基本假设4.建立数学模型、5.求解模型6.模型评价7.问题分析实验类型综合性实验需用仪器设备：pc486微机、mindows环境、LIND。软件、mathcad软件。实验步骤.建立工业区发展的数学模型.用LINDQ软件求解模型.写出工业区发展规划报告实验内容问题一：工业区发展规划某小城市的政府当局拟将一块300英亩的土地发展成工业区以增加税收，有5家工厂与政府当局签定了初步协议，市府允诺对这些工厂所需的资金，施工力量及长期出租有关厂房及设备方面加以支持。市府现有一切公用设施，且最近完成了一项为这些工厂所需的发展，以便需求能有节制地增长，可用资金最多六百万元。对该地劳动力情况的研究表明，可得到的本地临时工与外地工不超过表1所列数字：劳动力表1-1种类普通工半技工技工职员技术员管理员最大数40018012015012033平均年收入500070001200050001400030000研究表明，在缴纳州与国家的税收后，净收在本地消费的情况（计算是以每种收水平的家庭的平均大小为基础进行的）列于表1-2中。劳动力消费形式表1—2种类普通工半技工技工职员技术员管理员州，国税占总收入的百分比10%12%16%10%18%24%年净收入450061601010045001150022800流入市区的现金4050524070703830740012520还估算了城市从财产税与销售收入及公用设施、燃料、卫生设施及其它方面所得的税收，平均为每人在市内消费额的2.5%,市府从各种工人所得税收汇总于表1-3中：城市得自每个劳动力的税收（表1-3） 单位：元工人种类普通工半技工技工职员技术员管理员

平均税收10113117796187313城市拥有的公共设施可承担下列负荷：.电力；14.5百万度/年.水；75百万加仑/年..煤气：40百万立方英尺/年4.下水：40百万加仑/年品处理；2万吨/年城市的公共服务所得的单位利润；.电力：0.0007元/度 2.水：7元/百万加仑 3.下水：1元/百万加仑4.煤气：0.004元/立方英尺 5.废品；0.15元/吨市政府与之签订初步协定的五家拟建工厂，将主要使用本地原材料，加上工厂引起的其它利润虽将提高社会的收入水平，但因不能确定的加以定量，故在最优化对策中未于考虑。然而,在决定实施此计划时，将这些看作为有利的因素。表1-4给出了各工厂所需的最大与最小厂房面积。工 厂12345主要产品电子器件压力罐铸件塑料薄膜再生轮胎最大面积（平方英尺）120,00080,00070,000100,000150,000最小面积（平方英尺）40,00020,00026,00028,00035,000表1-5给出的是各工厂所需人力资源的数每千平方英尺厂房所需劳动力的工日数（每班8小时）工厂主要产品普通工半技工技工职员技术人员管理人员1电子器件0.750.230.380.300.150.052压力罐0.901.350.300.600.300.153铸件1.300.460.460.580.230.154塑料薄膜0.430.640.320.750.320.075再生轮胎1.370.390.11表1-6给出每千平方英尺厂房所需的场地及服务设施厂名主要产品场地（英亩）电力（度/年）水（百万加仑/年）煤气（立方英尺/年）下水（百万加仑/年）废品处理（吨/年）1电子器件0.5040,0000.05060,0000.040202压力罐1.0032,0000.10080,0000.058303铸件1.1562,0000.300120,0000.200504塑料薄膜1.4348,0000.20075,0000.150455再生轮胎1.4358,0000.100150,0000.06025表1-7给出每千平方英尺厂房的年产值厂名12345产品电子器件压力罐铸件塑料薄膜再生轮胎年产值（元/千平方英尺）50007000120001500010000市府将对各厂征收占其年总产值1.5%的税。该厂区的地价为：2000元/每英亩，各厂厂房的平均造价为:1厂（电子器件）一一15元/平方英尺 2厂（压力罐）——17元/平方英尺3厂（铸件）——19元/平方英尺 4厂■料薄膜）——22元/平方英尺5厂（再生轮胎）——16元/平方英尺市府已同意将厂房（包括场地）按每年L75元/平方英尺建筑面积租给工厂，维修费估计为0.05元/平方英尺建筑面积。市府怎样确定对各厂的空间分配，以使政府的税收为最多？问题二SARS的传播SARS（SevereAcuteRespiratorSyndrome,严重急性呼吸道综合症，俗称：非典型肺炎）是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到了许多重要的经验和教训，认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。请你们对SARS的传播建立数学模型，具体要求如下：（1）对附件1所提供的一个早期的模型，评价其合理性和实用性。（2）建立你们自己的模型，说明为什么优于附件1中的模型；特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型，这样做的困难在哪里？对于卫生部门所采取的措施做出评论，如：提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计。附件2提供的数据供参考。（3）收集SARS对经济某个方面影响的数据，建立相应的数学模型并进行预测。附件3提供的数据供参考。（4）给当地报刊写一篇通俗短文，说明建立传染病数学模型的重要性。附件1：SARS疫情分析及对北京疫情走势的预测2003年5月8日在病例数比较多的地区，用数理模型作分析有一定意义。前几天，XXX老师用解析公式分析了北京SARS疫情前期的走势。在此基础上，我们加入了每个病人可以传染他人的期限（由于被严格隔离、治愈、死亡等），并考虑在不同阶段社会条件下传染概率的变化，然后先分析香港和广东的情况以获得比较合理的参数，最后初步预测北京的疫情走势。希望这种分析能对认识疫情，安排后续的工作生活有帮助。1模型与参数假定初始时刻的病例数为平均每病人每天可传染K个人（人一般为小数），平均每个病人可以直接感染他人的时间为心天。则在原之内，病例数目的增长随时间《单位天）的关系是：N（t）=NOQ+W如果不考虑对传染期的限制，则病例数将按照指数规律增长。考虑传染期限4勺作用后，变化将显著偏离指数律，增长速度会放慢。我们采用半模拟循环计算的办法，把到达乙天的病例从可以引发直接传染的基数中去掉。参数人和原有比较明显的实际意义。同理解为平均每个病人在被发现前后可以造成直接传染的期限，在此期限后他失去传染作用，可能的原因是被严格隔离、病愈不再传染或死去等等。从原理上讲，这个参数主要与医疗机构隔离病人的时机和隔离的严格程度有关，只有医疗机构能有效缩短这个参数。但我们分析广东、香港、北京现有的数据后发现，不论对于疫情的爆发阶段，还是疫情的控制阶段，这个参数都不能用得太小，否则无法描写好各阶段的数据。该参数放在15-25之间比较好，为了简单我们把它固定在20（天）上这个值有一定统计上的意义，至于有没有医学上的解释，需要其他专家分析。参数K显然代表某种社会环境下一个病人传染他人的平均概率，与全社会的警觉程度、政府和公众采取的各种措施有关。在疾病初发期，社会来不及防备，此时K值比较大。为了简单起见，我们从开始至到高峰期间均采用同样的充值（从拟合这一阶段的数据定出），即假定这阶段社会的防范程度都比较低，感染率比较高。到达高峰期后，我们在10天的范围内逐步调整人值到比较小，然后保持不变，拟合其后在控制阶段的全部数据，即认为社会在经过短期的剧烈调整之后，进入一个对疫情控制较好的常态。显然，如果疫情出现失控或反复的状态，则人值需要做更多的调整。2计算结果对香港疫情的计算和分析。香港的数据相对比较完整准确。但在初期，由于诊断标准等不确切，在3月17日之前，没有找到严格公布的数据。我们以报道的2月15日作为发现第一例病人的起点，2月27日从报道推断为7例。3月17日后则都是正式公布的数据。累积病例数在图1中用三角形表示。我们然后用上述方法计算。4月1日前后（从起点起45天左右）是疫情高峰时期，在此之前我们取K=0.16204。此后的10天，根据数据的变化将K逐步调到0.0273,然后保持0.0273算出后面控制期的结果。短期内“调整的幅度很大，反映社会的变化比较大。图中实心方黑点是计算的累积病例数。从计算累积病例数，很容易算出每天新增病例数（当然只反映走向，实际状况有很大涨落）。可以看出，香港疫情从起始到高峰大约45天，从高峰回落到1/10以下（每天几个病例）大约40天（5月上中旬），到基本没有病例还要再经过近一个月（到6月上中旬）。对广东疫情的计算和分析。广东的起点是02年11月16日，到今年2月下旬达到高峰，经过了约100天。在今年2月10日以前的数据查不到，分析比较困难。总体上看，广东持续的时间比香港长得多，但累积的总病例数却少一些，这反映出广东的爆发和高峰都不强烈。但广东的回落也比较慢。从2月下旬高峰期到现在经过了约70天，还维持着每天10来个新增病例，而同样过程香港只用了约40天。这种缓慢上升和下降的过程也反映到尺值上。比较好的拟合结果是，在高峰期之前（t<101天），K=0.0892;在随后的10天逐步调整到0.031。用这组参数算出的后期日增病例数比实际公布的偏小，说明实际上降低得更慢。这种情况与疫情的社会控制状况有没有什么关系，需要更仔细的分析。对北京疫情的分析与预测。北京的病例起点定在3月1日，经过大约59天在4月29日左右达到高峰。我们通过拟合起点和4月20日以后的数据定出高峰期以前的弥=0.13913。这个值比香港的0.16204来得低，说明北京初期的爆发程度不如香港，但遗憾的是上升时间持续了近60天，而香港是45天，这就造成了累积病例数大大超过香港。从图2中还看出4月20日以前公布的数据大大低于计算值。而我们从对香港、广东情况的计算中，知道疫情前期我们的计算还是比较可行的。从而可以大致判断出北京前期实际的病例数。图中的公布数据截止到5月7日（从起点起67天），其后的计算采用的是香港情况下获得的参数。按这种估算，北京最终累积病例数将达到3100多。

图2对北京疫情的分析图3是计算的日增病例数。后期下降得较快的实心方黑点是采用香港参数获得的。这就是说，如果北京的疫情控制与香港相当或更好的话，就可以在高峰期后的40天（从起点起100天）左右，即6月上中旬下降到日增几例。然后再经过约一个月，即7月上中旬达到日增0病例。但如果北京的新病例下降速度与广东类似的话，则要再多花至少一个月，才能达到上述的效果，且累积总病例数会到3800左右。至于什么原因造成香港下降速度快而广东下降速度慢，需要有关方面作具体分析。

■采用齐港缢O采用广东缴图3北京日增病例走势分析3结论每个病人可以造成直接感染他人的期限平均在20天左右，这个值在不同地区和不同疫情阶段似乎变化不大。病人的平均每天感染率与社会状况有关，在疫情爆发期较大，在疫情控制期要小很多。香港的初期爆发情况比广东和北京都剧烈，但控制效果明显比较好。北京后期如果控制在香港后期的感染率水平上，则有望在6月上中旬下降到日增几例。然后再经过约一个月，即7月上中旬达到日增0病例。而累积总病例数将达到3100多。但如果北京的新病例下降速度与广东类似的话，则要再多花至少一个月，才能达到上述的效果，且累积总病例数会到3800左右。附件2:北京市疫情的数据(据：/Resource/Detail.asp?ResourcerD=66070)日期已确诊病例累计现有疑似病例死亡累计治愈出院累计4月20日33940218334月21日48261025434月22日58866628464月23日69378235554月24日7748633964

4月25日87795442734月26日988109348764月27日1114125556784月28日1199127559784月29日1347135866834月30日1440140875905月01日15531415821005月02日16361468911095月03日17411493961155月04日180315371001185月05日189715101031215月06日196015231071345月07日204915141101415月08日213614861121525月09日217714251141685月10日222713971161755月11日226514111201865月12日230413781292085月13日234713381342445月14日237013081392525月15日238813171402575月16日240512651412735月17日242012501453075月18日243412501473325月19日243712491503495月20日244412251543955月21日244412211564475月22日245612051585285月23日246511791605825月24日249011341636675月25日249911051677045月26日250410691687475月27日251210051728285月28日25149411758665月29日25178031769285月30日252076017710065月31日252174718110876月01日252273918111246月02日252273418111576月03日252272418111896月04日252271818112636月05日252271618113216月06日252271318314036月07日252366818314466月08日252255018415436月09日252245118416536月10日252235118617476月11日252325718618216月12日252315518718766月13日25227118719446月14日2522418919946月15日2522318920156月16日2521319020536月17日2521519021206月18日2521419121546月19日2521319121716月20日2521319121896月21日2521219122316月22日2521219122576月23日252121912277附件3:北京市接待海外旅游人数（单位：万人）年1月2月3月4月5月6月7月8月9月10月11月12月19979.411.316.819.820.318.820.924.924.724.319.4199818.619999.611.715.819.919.517.817.823.321.424.520.1200015.9200110.112.917.721.021.020.421.925.829.329.823.6200216.5200311.426.019.625.927.624.323.027.827.328.532.818.511.526.420.426.128.928.025.230.828.713.729.723.128.929.027.426.032.231.432.62