1、2009 高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D 中选择一项填写）： B 我们的参赛报名号为（如果
2、赛区设置报名号的话）： 所属学校（请填写完整的全名）： 三峡大学 参赛队员 (打印并签名) ：1. 王 恒 2. 湛敏峰 3. 童 圣 指导教师或指导教师组负责人 (打印并签名)： 指导教师组 日期：2009 年 9 月 14 日赛区评阅编号（由赛区组委会评阅前进行编号）：2009 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：1眼科病床合理安排摘要本文讨论了病床的合理安排问题，属于优化问题中的排队问题。我们根据
3、始数据利用 EXCEL 软件进行了统计分析，得出各类眼科病人的平均等待时间等相关数据信息。对于问题一，我们综合考虑医院与病人的利益，提出了平均病床周转次数、病人住院平均等待时间 、等待住院病人队列长度 、等待住院病人队列ABC变化趋势这四项评价指标，用以对病床安排模型的优劣进行评价。并利用该评价指标体系对医院当前的病床安排模型进行了评价。对于问题二，我们基于医院的当前情况，以平均病床周转次数 为优化目A标，以改进后的优先非抢占排队思想为依据，采用优先级随时间变化的规则来进行病床安排，并根据五类眼科病人的平均住院时间设置了初始优先度值，建立起单目标优化模型一。我们利用模型一对前来门诊的病人重新进
4、行病床安排，得出了相关结果。由结果我们可以看出，模型一可以较好的解决医院的等待住院病人队列越来越长的问题。我们利用问题一里确定的评价指标体系对模型一进行了评价，并将其与医院当前采用的模型进行了对比分析，突显出模型一的优势。对于问题三，我们根据问题二里得出的病人信息，统计出了各类病人的平均等待时间和等待队列长度，发现在模型一的病床分配方案下，每天门诊总病人数与出院总人数大致平衡。于是，我们可以根据各类病人的等待时间分布来给出门诊病人的入院时间区间：外伤：1 天；视网膜疾病：（10，15）天；青光眼：（7，12）天；白内障单眼：（4，8）天。白内障双眼病人需视门诊时间而定。对于问题四，在周六、周日
5、不安排手术的情况下，利用模型一重新对病人进行入院安排，并用评价指标体系对结果进行了评价，发现分配结果并不理想，等待队列长度很长，且等待入院的病人队列会越来越长。因此，我们认为医院手术时间应该调整，我们建议将白内障双眼病人的手术时间由原来的每周一、周三调整到每周三、周五。对于问题五，我们利用多服务台排队系统 来进行求解。我们以平cM/均逗留时间最短为目标函数，建立起病床优化分配模型，并通过 MATLAB7.0 软件进行编程求解，得出病人在系统内的平均逗留时间最短为 13.92 天，五类疾病病床的分配比例如下表所示：疾病类型 外伤 视网膜疾病 青光眼 白内障单眼 白内障双眼分配病床数 8 37 1
6、1 8 15我们在模型改进里对模型一中可能出现的多种可行解提出了一种评价方法，通过定义医生工作强度值 来对不同可行解进行评判。Q2关键词：病床安排 排队论 改进的优先非抢占排队思想 多服务台排队系统cM/1.问题重述1.1 问题背景到医院就医时，一般要经历以下过程：到门诊就诊、到收费处划价、到药房取药、到注射室打针、等待住院等。而在这些过程里，往往需要排队等待以接受服务。现有某医院眼科门诊住院部的大致情况如下：该医院眼科门诊每天开放，住院部共有 79 张病床，当前对全体非急症病人是按照 FCFS（First come, First serve）的规则来安排住院的，但是等待住院的病人队列却越来越
7、长。因此，解决该住院部的病床合理安排问题，提高对医院资源的有效利用，以减少等待住院的病人队列就显得很有必要了。1.2 该医院眼科门诊相关信息该医院眼科手术主要分为四大类：白内障、视网膜疾病、青光眼和外伤。白内障手术比较简单，没有急症。目前该院是每周一、三做白内障手术，此类病人的术前准备时间为 1-2 天。其中，做两只眼的病人比做一只眼的要多，大约占到 60%。如果要做双眼，便安排周一先做一只，周三再做另一只。外伤疾病通常属于急症，有空的病床时就立即安排住院，住院后第二天便会安排手术。其他眼科疾病比较复杂，有各种不同情况，但大致住院以后 2-3 天内就可以接受手术，主要是手术后的观察时间较长。这
8、类疾病手术时间可根据需要安排，一般不安排在周一、周三。此类疾病急症数量较少，可不予考虑。该医院眼科手术条件比较充分，但考虑到手术医生的安排问题，通常情况下白内障手术与其他眼科手术（急症除外）不安排在同一天做。1.3 需要解决的问题题目附录中给出了 2008 年 7 月 13 日至 2008 年 9 月 11 日这段时间里各类病人的情况。根据这些病人情况和该医院眼科门诊的相关信息，我们需要通过采用数学建模的方法来帮助解决以下问题：问题一：分析并确定出合理的评价指标体系，用以评价该问题的病床安排模型的优劣。问题二：试就该住院部当前的情况，建立合理的病床安排模型，以根据已知的第二天拟出院病人数来确定
9、第二天应该安排哪些病人住院。并对我们的模型利用问题一中的指标体系作出评价。问题三：作为病人，自然希望尽早知道自己大约何时能住院。能否根据当时住院病人及等待住院病人的统计情况，在病人门诊时即告知其大致入住时间区间。 问题四：在该住院部周六、周日不安排手术的情况下，对问题二重新回答，并讨论医院的手术时间安排是否应作出相应调整?问题五：有人从便于管理的角度提出建议，在一般情形下，医院病床安排可采取使各类病人占用病床的比例大致固定的方案，试就此方案，建立使得所3有病人在系统内的平均逗留时间（含等待入院及住院时间）最短的病床比例分配模型。2.问题分析医院就医排队是大家都非常熟悉的现象，题目中该医院眼科门
10、诊对全体非急症病人的住院安排所采取的先到先服务规则，是排队规则中应用最为普遍的一种情形，这种规则能够很好的体现出公平性。然而，从实际情况看来，在这种规则下，等待住院的病人队列却越来越长。因此，改进病床安排规则，更加合理的利用床位资源，对于解决实际的问题是非常重要的。2.1 问题一的分析问题一要求我们给出合理的评价指标体系，用以对病床安排模型的优劣进行评价。我们首先找出了与病床工作相关的指标，从中筛选出一些指标来作为评价指标体系，用以对病床工作效率进行评价。另外，我们还可以考虑从病人的角度来确定评价指标，使病床分配模型被更多的病人所接受。最终，我们确定了四个评价指标：平均病床周转次数、病人住院平
11、均等待时间、等待住院病人队列长度、等待住院病人队列变化趋势。2.2 问题二的分析问题二要求我们给出一种合理的病床安排模型，根据已知的第二天拟出院人数来安排病人住院。我们已经知道，医院目前所采用的先到先服务规则无法解决等待人数过多的实际情况。因此，我们重新设定住院病人的排队规则，根据各类病人的平均住院时间对各类病人设定初始优先度值，优先度值越低，其优先级越高，优先级高的病人先住院。同时，为了不使优先级低的病人长期排队而无法入院，我们规定病人每等待一天，其优先度值便减 1。根据这种规则，我们以平均周转次数为目标函数建立优化模型，然后用问题一所确定的评价指标体系对模型进行评价，突显模型优势。2.3
12、问题三的分析问题三希望可以让病人在门诊时便得知大致住院时间区间，这就需要对当时已住院人数和等待住院人数进行统计，并对已住院病人的出院日期进行预测，然后考虑在具体的病床分配模型下，根据每天门诊人数和病人具体情况和出院人数来对未来的病床使用进行大致安排。2.4 问题四的分析问题四假定该医院住院部在周六、周日不安排手术，则原先可以安排在周六、周日的手术必须安排到周一至周五去，这就使得等待住院病人队列更长，也必将导致病床的合理安排更为复杂。在这种情况下，重新对问题二进行回答，然后根据评价结果来确定医院的手术安排是否应该作出相应调整。2.5 问题五的分析如果按照建议将病床按一定比例分配给各类病人，则对于
13、不同的病人，其所对应能够入住的病床数便会有限制，问题可以等效为多服务台排队问题，我4们可以利用排队论中的多服务台排队系统 模型来解决这一问题。我们/Mc以系统内所有病人的平均逗留时间为优化目标建立优化模型，利用 MATLAB7.0软件进行编程求解，得出病床分配比例与平均逗留时间。3.模型假设1. 假定该医院眼科手术条件充分，不必考虑手术条件对病床安排的限制；2. 假定不需要考虑白内障、青光眼、视网膜疾病的急诊情况；3. 假定对于外伤病人，只要一有空病床，便安排其住院，并在住院第二天后安排手术；4. 白内障双眼手术时间必须间隔一天；5. 假定题目附录中给定的数据真实可靠，具有较好的代表性；6.
14、我们用序号 1-61 来对应表示 2008 年 7 月 13 日至 2008 年 9 月 11 日。4.符号说明平均病床周转次数；A:住院平均等待时间；B等待住院病人队列长度；C：统计时间内的出院人数；cR：统计时间；T：平均开放床位数；w：眼科疾病类型， 时，分别对应外伤、视网膜疾病、青光眼、白i 1,2345i内障单眼、白内障双眼；：病人的门诊时间；Tm：病人的入院时间；r：病人的手术时间；s：病人手术后的留院观察时间g：求余运算式；Re(,)n：各类病人每天前来就诊的人数；i：各类病人的平均住院时间的倒数；i：各类病人在这段时间里前来就诊的总人数；ia：种类病人对应分配床位数；ic：队列
15、长度期望值；sL：期望等待时间；qW5：期望逗留时间。sW5.模型的建立与求解5.1 数据处理题目附录中给出了 2008 年 7 月 13 日至 2008 年 9 月 11 日这段时间里各类病人的情况，我们利用 EXCEL 软件对这些数据进行了相关统计分析，得到了相关信息。5.1.1 各类病人就医平均时间根据统计分析可得，五类病人在就医过程中所经历的平均等待时间、平均准备时间、平均恢复时间和平均住院时间如下表 1 所示。表 1：各类病人就医平均时间疾病类型 平均等待时间 平均准备时间 平均恢复时间 平均住院时间外伤 1.00 1.00 6.04 7.05视网膜疾病 12.54 2.38 10.
16、17 12.54青光眼 12.26 2.41 8.08 10.49白内障单眼 12.67 2.33 2.90 5.28白内障双眼 12.51 3.60 2.96 8.56其中，等待时间为病人入院时间与门诊时间的差值；准备时间为病人做第一次手术时间与入院时间的差值；恢复时间为出院时间与出院前最后一次手术时间的差值；住院时间为出院时间与入院时间的差值。从上表 1 中，我们可以看出，除了外伤病人的平均等待时间为 1 天外，其他眼科病人从接受门诊到住进医院的平均时间都至少需要 12 天。这五种病人的手术前准备时间与题目中提到的白内障 1-2 天、外伤 1 天、其他疾病 2-3 天大致吻合，只有做白内障
17、双眼手术的病人平均术前准备时间略有出入，至少需要3 天。各类病人的平均恢复时间，白内障最短为 3 天，视网膜疾病最长，至少需要 10 天。五类病人的平均住院时间长短依次为白内障单眼、外伤、白内障双眼、青光眼、视网膜疾病，其中，视网膜疾病病人平均住院时间最长，至少需要 12 天。5.1.2 各类病人每天门诊人数根据统计分析，我们得出五类病人在 2008 年 7 月 13 日至 2008 年 9 月 11日这段时间里每门诊的人数，见附录表 2 所示。根据表 2 中的数据，我们利用 EXCEL 软件进行统计分析。我们采用 6SQ 统计插件中的假设检验下的卡方拟和优度来对其进行泊松分布的检验，在显著性
18、水平 时，发现五类病人每天门诊的人数和总门诊人数都服从泊松分布。0.1表 3：泊松分布检验白内障 白内障（双眼） 视网膜统计量 统计量 统计量数据个数 61 数据个数 61 数据个数 61总和 100 总和 133 总和 1706最大值 5 最大值 7 最大值 7平均值 1.6 平均值 2.2 平均值 2.8假设检验 假设检验 假设检验零假设服从泊松分布 零假设服从泊松分布 零假设服从泊松分布自由度 3 自由度 4 自由度 4卡方统计量 1.8 卡方统计量 1.9 卡方统计量 3.7p 值 0.6 p 值 0.8 p 值 0.5显著性水平 0.1 显著性水平 0.1 显著性水平 0.1结果 接
19、受零假设 结果 接受零假设 结果 接受零假设外伤： 青光眼： 所有病人：统计量 统计量 统计量数据个数 61 数据个数 61 数据个数 61总和 64 总和 63 总和 530最大值 3 最大值 4 最大值 16平均值 1 平均值 1 平均值 8.7假设检验 假设检验 假设检验零假设服从泊松分布 零假设服从泊松分布 零假设服从泊松分布自由度 2 自由度 3 自由度 7卡方统计量 1 卡方统计量 4 卡方统计量 5.1p 值 1 p 值 0 p 值 0.6显著性水平 0.1 显著性水平 0.1 显著性水平 0.1结果 接受零假设 结果 接受零假设 结果 接受零假设5.1.3 其他统计数据我们通过
20、对题目附录中的病人数据信息进行分析，发现从 2008 年 7 月 13日到 2008 年 8 月 7 日之间的 26 天里，入院人数总和（126 人）与出院人数总和（47 人）的差值为 79，刚好等于该住院部的病床总数。由此，我们可以知道，在 2008 年 8 月 7 日的时候，该住院部的病床已经全部被占用。具体每天入院出院人数见附录表 4。从表 4 中，我们可以看到，从 2008 年 8 月 8 日开始（在表3 中序号为 27） ，每天入院和住院的人数都相等。另外，我们可以看到，2008 年 7 月 13 日前来门诊（外伤除外）的病人（外伤病人除外） ，直到 2008 年 7 月 25 日才
21、被安排入院，由此我们也可以知道，在这段时间内，医院的病床基本上也是全部被占用了。5.2 评价指标体系（针对问题一）在本题里，选定评价指标时，我们既需要考虑这些指标对于医院利益的影响，也应该从病人的角度出发，提高病人对于医院服务的满意度。从医院角度7考虑，应该使每张病床的工作效率提高，而与工作效率相关的指标 1有平均病床工作日和平均病床周转次数。从病人角度来说，自然是希望等待住院的时间越短越好。于是，我们可以尝试从上述两个方面来设定评价指标，用以反映病床安排模型的优劣。通常情况下，医院为了自身的利益，会希望有更多的人前来做手术，而不是有病人长期占用病床，因为相对住院费用来说，手术费用更能增加医院
22、的利益。因此，从医院角度来说，单位时间里越多的病人住院，会使其利益越高。这可以反映到平均病床周转次数上来，平均病床周转次数越高，说明单位时间内住院做手术的病人越多。对于病人来说，把平均等待时间作为评价指标便可以很好的反映出其利益，病人等待住院的时间越短，他们对医院的满意度就会超高。另外，根据题目里该医院眼科门诊的具体情况，我们主要应该解决的是等待队列越来越长的问题，因此，我们设想将等待住院病人的队列长度也作为评价指标，这样就可以结合本题目的具体情况更加直接的反映出病床安排的合理性。考虑到等待住院病人的队列长度是一个不断变化的值，我们需要另外再设定一个指标来反映其长期的发展状况，我们取定用病人等
23、待住院队列的变化趋势来作为指标反映这一情形。我们一共得到四个评价指标：1.平均病床周转次数 ：指“出院人数”与“平均开放床位数”之比。A平均病床周转次数=出院人数/平均开放床位数，在本题中，我们用单位时间内每张病床平均每天接待的病人数量来表示平均病床周转次数。则平均病床周转次数=出院人数/（统计时间*平均开放床位数） 。2. 病人住院平均等待时间 ：指所有病人从门诊到入院的等待时间的平均B值。病人住院平均等待时间=（入院时间总和门诊时间总和）/总人数。3.等待住院病人队列长度 ：每天等待住院的病人人数。C4. 等待住院病人队列变化趋势：等待住院病人人数的变化趋势，我们用图表来表示这一变化趋势。
24、综上所述，我们可以将上述四个评价指标用如下表 5 来表示。表 5：评价指标体系评价指标平均病床周转次数 A病人住院平均 等待时间 B等待住院病人 队列长度 C等待住院病人队列变化趋势含义每张病床平均每天接待的病人数量病人从门诊到入院的平均等待时间每天等待住院的病人人数等待住院病人人数的变化趋势（图表表示）其中，因为题目给定的数据是从 2008 年 7 月 13 日到 2008 年 8 月 7 日这一段时间里的病人信息，因此，我们将平均病床周转次数平均到每一天上。因为图表可以更加明确的反映出变化趋势，所以我们要求用图表形式来表示等待住院病人的队列变化趋势。5.2.1 对医院当前病床安排模型的评价
25、现在，我们来利用确定出的评价指标体系对医院目前的病床安排模型进行评价。根据统计数据，我们可以得出，在 2008 年 7 月 13 日到 2008 年 8 月 7 日这861 天时间里，一共有 349 个病人出院，由此我们可以算出平均病床周转次数。在前面的数据处理中，我们已经得出了各类病人的平均等待时间。对于等待住院的病人队列，我们只取 2008 年 9 月 11 日时的人数来表示这个队列长度，即为 102 人。我们将每天的等待住院病人队列长度值在 MATLAB7.0 软件上进行三次曲线拟和，得到其变化趋势图形。各指标值如下：平均病床周转次数： 0724.9*613A病人平均等待时间，见下表 6
26、：表 6：各类病人平均等待时间疾病类型 外伤 视网膜疾病 青光眼 白内障单眼 白内障双眼平均等待时间 B1.00 12.54 12.26 12.67 12.51等待住院病人队列长度： 102C等待住院病人队列变化趋势，见下图 1：图 1：当前模型下等待住院病人队列变化趋势由评价结果，我们可以看出，在医院当前所采用的病床安排模型下，等待住院的病人队列一直呈上升趋势，住院病人等待队列会越来越长。如果按照这个趋势发展下去，将不利于医院自身的发展，因此，必须寻找更为合理的病床安排模型来解决这一问题。5.3 病床安排模型的建立与求解（针对问题二）当平均病床周转次数增大，也就是说出院病人的数目在增加，医院
27、就可以进行更多的手术，医院的利益将会增大。因此，我们可以将平均病床周转次数9作为优化目标来建立模型。当出院人数增多，就可以有更多的病人住院，这样就可以减少等待住院的病人数量，从而使住院病人的等待队列逐渐减小。因此，在安排病人住院的时候，我们可以尽量先安排住院时间短的病人住院，使单位时间内出院人数增加，以使等待入院的病人数量减少，从而减少等待住院的病人队列长度。因此，我们可以按照各类病人的平均住院时间来对各类病人设定优先权，让平均住院时间短的病人先住院，这样就可以增大单位时间内平均病床周转次数。然而，如果单纯的根据各类病人的平均住院时间长短来安排病人入院时间，就可能导致平均住院时间长的病人始终停
28、留在等待入院的队列中，这肯定是不合理的。从病人的角度来看，随着病人在等待入院的队列中等待时间的不断增加，病人对于医院的不满意度也会随之增长，这就可能导致病人转投其他医院。对于医院来说，当然是不希望看到病人流失的。为了均衡上述矛盾，在提高平均病床周转次数的同时，又保证病人不流失，我们提出一种基于优先非抢占排队思想的改进排队模型。5.3.1 基于优先非抢占排队思想的改进排队模型优先非抢占排队系统 2是指各类顾客有着不同的优先服务权，同类顾客仍然按照先到先服务的规则进行服务，当优先级高的顾客到达系统时，不中断正在服务的优先级低的顾客的服务，而是等到该优先级低的顾客服务完离开后才使用服务台。在本题目的
29、具体情况下，我们不能直接采用优先非抢占排队思想来建立模型，因为这就相当于只考虑了平均住院时间短的病人优先住院，虽然可以使平均病床周转次数增大，却可能导致优先级低的病人长期无法住院。因此，我们考虑对这种非抢占排队思想进行改进，以解决这一问题。改进后的排队思想：我们令顾客的优先级不断变化，即优先级较低的顾客在排队等待的过程中，随着其排队时间的不断增加，我们使其优先级也跟着上升，这样就可以保证优先级较低的顾客在经过一定时间的排队等待后优先级上升到较大值，从而可以接受到服务，而不是一直排队等待下去。5.3.2 模型一建立根据上述改进后的排队思想模型，我们可以采用这种思想建立病床分配模型。我们根据统计出
30、来的五类病人的平均住院时间，运用“去尾法”处理之后，将得到的整数值作为各类病人在门诊后等待住院时所对应的初始优先度。表 7：各类病人初始优先度值疾病类型 外伤 视网膜疾病 青光眼 白内障单眼 白内障双眼初始优先度值 12 10 5 5其中，我们将外伤一律看作是急诊，赋其优先度值为 。对于白内障病人，为方便计算，我们将做单眼手术和双眼手术的病人同等看待，取做单眼手术的平均住院时间，经“去尾法”处理后作为它们的初始优先度值。排队规则：我们规定，所有病人在门诊之后，等待的时间每增加一天，其优先度降低101（优先度可以是负数） 。医院安排病人入院的时候，在满足术前准备时间（白内障 1-2 天，其他病科
31、（不含急诊）2-3 天）的情况下，按照病人的优先度值从小到大的顺序安排病人入院。同时，考虑到急诊病人需要当病床有空时立即安排住院，我们规定急诊病人的优先度为 ，即保证无论急诊病人何时门诊，都能在病床有空时立即安排住院。则我们可以根据上述排队规则建立如下病床分配模型：目标函数： ()*cwRMinATC约束条件： 1, 124,5,3230,6Re(,7)2,4513,2,3.1634,5183102iii iTsriiimsistTrigi其中， 表示平均病床周转次数。我们用序号 1-61 来对应表示 2008 年 7A月 13 日至 2008 年 9 月 11 日这段时间里的每一天。 表示统
32、计时间内的出院cR人数， ，表示统计时间， ，表示平均开放床位数， 表示病人61T79wCTm的门诊时间， 表示病人的入院时间， 表示病人的手术时间， 表示病人rTsg的手术后的留院观察时间。 表示求余运算，因为 2008 年 7 月 13 日Re(,)mn是星期日，对应序号为 1，我们便可以用 的值来表示各天时间所对(,7)is应的星期值，其余数可为 0、1、2、3、4、5、6，分别对应星期六、星期日、星期一、星期二、星期三、星期四、星期五。5.3.3 模型一求解我们根据题目中给定的 2008 年 7 月 13 日至 2008 年 9 月 11 前来门诊的各类病人信息，利用我们建立的模型一，
33、重新安排病人住院，得到了一份类似题11目附录中给定的病人信息，详细信息见附录表 8。根据表 8 中的相关信息，我们同样对其进行统计分析，得到如下数据：表 9：各类病人就医平均时间疾病类型 平均等待时间 平均准备时间 平均恢复时间 平均住院时间外伤 1.00 1.00 6.00 7.00视网膜疾病 13.94 2.09 9.94 12.03青光眼 11.16 2.16 7.86 10.02白内障单眼 8.07 1.70 3.00 4.70白内障双眼 8.78 1.83 3.00 6.83从表 8 中，我们可以看到，利用模型一来进行病床分配后，截止 2008 年 9月 11 日，这 61 天里前来
34、门诊的病人已经有 382 人出院，等待住院的病人人数为 69 人。5.3.4 对模型一的评价现在，我们来利用评价指标体系对模型一的结果进行评价。根据统计数据，我们可以得出，在 2008 年 7 月 13 日到 2008 年 8 月 7 日这61 天时间里，一共有 382 个病人出院，由此我们可以算出平均病床周转次数。对于等待住院的病人队列，我们同样只取 2008 年 9 月 11 日时的人数来表示这个队列长度，即为 69 人。我们将每天的等待住院病人队列长度值在 MATLAB7.0软件上进行三次曲线拟和，得到其变化趋势图形。各指标值如下：平均病床周转次数： 0793.*61382A病人平均等待
35、时间，见下表 6：表 10：各类病人平均等待时间疾病类型 外伤 视网膜疾病 青光眼 白内障单眼 白内障双眼平均等待时间 B1.00 13.94 11.16 8.07 8.78等待住院病人队列长度： 9C等待住院病人队列变化趋势，见下图 2：图 2：模型二下等待住院病人队列变化趋势12由评价结果，我们可以看出，如果医院采用我们的的病床安排模型一对病床进行分配，等待住院的病人队列变化会呈现下降趋势，住院病人的等待队列长度会相应减小。按照这个趋势发展下去，医院的病人等待队列越来越长的问题将得到较好缓解乃至解决，2008 年 9 月 11 日这一天的病人等待队列长度为69 人。因此，我们认为模型二的病
36、床分配模型是比较合理的。5.4 两种模型的评价结果比较为了更好的说明我们所建立的病床分配模型的优越性，我们来对上述两种模型的评价结果进行比较分析，见下表 11。表 11：医院当前采用模型与模型一的评价结果比较病人住院平均等待时间 B评价指标平均病床周转次数 A外伤 视网膜疾病青光眼白内障单眼白内障双眼等待住院病人队列长度 C等待住院病人队列变化趋势当前模型0.0724 1.00 12.54 12.26 12.67 12.51 102 呈上升 趋势模型一0.0793 1.00 13.94 11.16 8.07 8.78 69 呈下降 趋势由上表 11 可以看出，相对于医院目前所采用的病床安排模型
37、而言，我们所建立的病床安排模型一，除了视网膜病人的平均等待时间略有上升外（由平1312.54 天变为 13.94 天） ，而其相差也只有一天左右，我们认为这是可以接受的。至于其他各项指标，模型二的评价结果均远远优于医院当前所采用的模型。平均病床周转次数由 0.0724 上升到 0.0793，青光眼病人的平均等待时间减少了一天左右，白内障病人（包括做单眼和做双眼手术的病人）的平均等待时间减少了 4 天左右，在 2008 年 9 月 11 日这一天，等待住院的病人队列长度由原先的 102 人减少到 69 人。另外，从等待住院的病人队列变化趋势来看，医院当前所采用的模型会导致病人等待住院的队列越来越
38、长，而由我们的模型一对病床进行安排后，病人等待住院的队列将逐渐减小。从上述两种模型的对比分析中，我们可以很清楚的看到模型一的优越性，同时，从模型一的评价结果看来，病人等待住院的队列在逐渐减小，从具体等待队列的长度来看，在 2008 年 9 月 11 日，这个队列长度为 69 人，与医院当前具体情况的 102 人相比，减少了 33 人，减少比率为 32.35%。通过以上具体的数据分析，我们的模型一可以较好的解决医院目前的问题，即减少病人住院等待队列的长度。5.5 根据门诊时间确定住院时间区间（针对问题三）我们利用表 8 中的病人信息，取后一个月的医院门诊人数和出院人数作统计。我们发现，在我们的模
39、型运行一段时间后，系统基本上趋于稳定。总体的门诊人数和出院人数大体上趋于平衡,如图 3 所示。所以，采用模型一来进行病床分配，一段时间后等待住院的病人队列将基本上保持在一个相对稳定的状态。图 3：门诊人数与出院人数关系我们再取表 8 中 8 月 11 号（序号为 30）到 8 月 31 号（序号为 50）的数据来进行统计分析，得出各类疾病每天对应的病人等待入院队列长度和该类病人14的平均等待时间，如表 12 所示。表 12：各类病人等待时间与列长外伤 视网膜 青光眼 单眼 双眼日期队列长度等待时间队列长度等待时间队列长度等待时间队列长度等待时间队列长度等待时间总队列长度30 0 37 13 8
40、 11 13 14 5 7231 0 40 14 6 10 5 17 4 7332 0 44 15 7 9 11 6 22 10 8433 1 1 41 15 6 12 11 25 9 8434 0 43 14 8 11 12 8 27 8 9035 1 1 45 13 7 6 7 14 7 7336 3 1 31 13 7 5 6 15 6 6137 2 1 17 12 8 11 9 7 19 6 5538 0 18 13 6 11 9 6 23 11 5639 2 1 42 13 6 11 5 24 10 8540 1 1 33 13 6 12 9 24 7641 1 1 29 12 4
41、13 8 26 8 7342 1 1 26 12 6 10 11 7 11 7 5543 0 29 11 4 10 9 11 6 5344 1 1 28 10 5 9 9 6 13 5 5645 2 1 27 9 3 11 5 14 4 5746 3 1 30 13 3 10 4 15 4 6147 0 31 13 4 7 11 5 18 9 6448 2 1 24 13 5 6 13 8 21 9 6549 0 26 12 5 5 16 7 8 8 5550 0 27 11 7 11 16 6 9 7 59其中，空格表示该天没有该类病人前来门诊。通过分析，我们发现外伤的等待时间都是 1 天，
42、而随着各类病人等待住院队列长度的不同和总队列的队长不同，视网膜疾病，青光眼，白内障单眼的等待时间基本上都趋于稳定，大致分布在区间（10，15） ， （7，12） ， （4，8）内。我们发现白内障双眼的等待时间波动性很大，因为该类病人的等待时间还与星期有关，主要原因是白内障双眼的手术时间只能是星期一，星期三。为了减少队列长度，依据我们的模型，在安排白内障双眼病人时优先在星期六进行安排，所以白内障双眼病人的住院时间基本上都安排在星期六。所以，在预知白内障双眼病人的大致入住时间区让时，不但要依据等待住院病人的统计情况，还需要根据当时的星期值来大致判断其入住时间，我们用公式表示如下：1566,2013
43、,27abc其中， 表示门诊当天星期值， 表示白内障双眼的等待人数， 表示等待abc时间。则我们可以根据当时住院病人及等待住院病人的统计情况，告知其大致的入住时间区间：1. 外伤：1 天；2. 视网膜疾病：（10，15）天；3. 青光眼：（7，12）天；4. 白内障单眼：（4，8）天；5. 白内障双眼：66,2013,27abc其中， 表示门诊当天星期值（例 时表示星期一） ， 表示白内障双眼a1ab的等待人数， 表示等待时间。c5.6 手术时间安排的调整（针对问题四）在该住院部周六、周日不安排手术的情况下，我们根据模型一来重新对病床进行分配，完成了对病人入院时间的重新安排。我们将每天的等待住
44、院病人队列长度值在 MATLAB7.0 软件上进行三次曲线拟和，得到其变化趋势图形。图 4：问题四中等待入院病人队列变化趋势16与问题二中所得到的结果相比较，在该住院部周六、周日不安排手术的情况下，病人等待队列长度变大，且有逐渐上升的趋势。我们用评价指标体系来对这种情况下的病床分配模型进行评价，各指标值如下：平均病床周转次数： 3450.7261*9A病人平均等待时间，见下表 13：表 13：各类病人平均等待时间疾病类型 外伤 视网膜疾病 青光眼 白内障单眼 白内障双眼平均等待时间 B1.00 16.07 13.82 8.99 9.75等待住院病人队列长度： 。127C我们发现，当住院部周六、
45、周日不安排手术时，模型一的分配结果并不理想，在 2008 年 9 月 11 日，等待队列为 127 人，等待入院的病人队列同样样会越来越长。分析出现该种情况的原因，当住院部周六、周日不安排手术时，由于术前准备时间的约束，导致每周四、周五的时候，即使住院部有空病床，除了外伤病人外无法安排其他病人入院，造成医院病床资源的浪费，这是应该避免的情况。而到了每周六、周日的时候，又会出现入院“拥挤”的状况，使得部分病人的等待时间相对于问题二中要长。由此可以判定，在我们建立的病床安排模型一下，如果该住院部周六、周日不安排手术，需要对医院的手术时间安排作出相应调整。为避免病床资源浪费以及入院“拥挤”的现象，我
46、们认为最简单方便的手术时间调整方案，是将白内障病人的手术时间由原来的每周一、周三调整到每周三、周五，即如果要做白内障双眼手术，则周三先做一只，周五再做另一只。17这样便能避免在某一天出现大量空病床，却无法安排外伤病人之外的病人入院。调整方案：将白内障病人的手术时间由原来的每周一、周三调整到每周三、周五。5.7 平均逗留时间最短的病床比例分配模型的建立与求解（针对问题五）依据建议，医院可以将总数为 79 的病床，按照各类病人占用病床的比例来大致进行固定分配。在这种方案下，我们可以把病床安排问题等效看作是多服务台排队问题。其中病床可以表示为服务台 ，每天前来就诊的病人数可以表c示为顾客到达率 ，各
47、类病人平均住院时间的倒数可以表示为每个服务台的平均服务率 。这样，我们就可以利用多服务台排队系统 来对问题五进 cM/行求解。5.7.1 多服务台排队系统 相关知识 2cM/在多服务台排队系统中，服务台用 表示，系统内顾客数为 ，顾客的平均n到达率为常数 ，每个服务台的平均服务率 是相同的，且规定各服务台的工作是相互独立的。对于整个服务机构来说，其平均服务率与系统状态有关，可表示如下： cncn, 1,2为了使系统不会排成无限队列，需要使服务强度 。在分析这个排队系统时，仍然从状态间的转移关系开始，其状态平衡方程为 cnPcPcnnn,)(1)1(0由以上差分方程可求得， 1100()()!(
48、)*,1!1,ccnnnncPnc系统的队列长度期望值可以求得，如下所示： 020()()*!1cqnnLPP由李特尔公式 3即可求出队长期望值 ，期望等待时间 和期望逗留时间sLqW。sW18综合上述结果，我们可以把 系统的主要绩效测度指标归纳如下：cM/101002 ()*,1()!(),! ,()*,!1ncnnn ncc qsq sqPncPcLLPLW 5.7.2 模型二建立我们可以采用上述多服务台排队系统 模型来解决问题五。cM/我们假定五类病人的服务台数分别为 ，分别对应其分配到的54321,病床数；我们用顾客到达率 来表示五类病人每天前来就诊的人数；54321,用五类病人的平均住院时间的倒数来表示每个服务台的平均服务率。54321,则期望逗留时间可以表示如下： 102 120()*()()1!1*()!cc cq ns nPLW 根据题目中给定的 2008 年 7 月 13 日至 2008 年 9 月 11 日这段时间里各类病人的情况，我们可以得出下列数据，见表 12。表 14：各类病人服务系统值外伤 视网膜疾病 青光眼 白内障单眼 白内障双眼到达率 i1.0492 2.7869 1.0328 1.6393 2.1803服务率 0.1420 0.0822 0.0953 0.2500 0.1667总人数 ia55 101 39 72 82其中， 表示
