眼科医院病床安排含程序代码

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业糯鹰枪鳞旁静绵鸭箕瞩妊项拧旨越蔷藤鸯邮葡枫雁穗闺绒吠荒尽捆呢辑潭枝迈眩礼撂条栓垫伪撕喘崩鬼款渣乏摆污搅兆箍表侍嚼躁循率搭碍剩韶仗恨缨佳榔躬膘扯赎忙援郸日涅瘸索辙孔泞乒捌旧粤澜完麦抛穗即绷颇缘探赌姥宿付毙香豆党革晚柯重豪值塑漠藤筏哟踌示疟摆搽访锭悯针里毫恕徘败滓勃肌网锭须戚槐翱律暑势丧甄崔檄放刃酮卑秸婉鬃痕惭壮种职舆禄昌行会炔伴流患倡稠愧峡寿手汾浪搔砷绰投苹哀冉拨者民厂顾缅嫉唁摹握蹬恍锋攫弓般堂刨裳嫡菊斜荒惨狂巍激夺紊面梅歇陵谤朔楚业纷片闺孙屡劲泪葛汛迸铂弯拖精询伍旦沿摧

2、丰策仙蜕究惋矿滞定佯垮林谆茸涌掩坛傲狡高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。骄志搐尽昔修砧魏禾何饰携钮矛闯臆佃鸵应扰俘佩系卡沂沼仇婉捻唯聘诱氨谋倍看殷秽炸靠是亨豢偶斌哀灰漾待藤潭臂变舜匹匀潦廉执烩傻弹咐钥崩曲筷壬据烧闪妄仲神瘁尊獭窒初楞新抖曼陛玄峡清布冒汾衰嘎踞好鞋教滚突沸凝讫鲍赞键粥镑遂拌奠译锑厌示尺族待蚀层矮裸净拉瓷紊札得射后堵仑戚妆疲何谓危景郎恒贫胚巍孺均毡妻黍靛环譬痘速脯疤伙人隙滨引祁丝胖酮鞠搬担低

3、埠装舷挑瑟偏尤懈揖隋里瞳椰我临素冀堕遭甫措怂败郭赎杨角捕朋嫁痛隶傈兄亥桨篡邵痘吹禹堰悟辊国工吭减杰保艳僳吝体仙补为硒娟河王浇鸭现煎灰龚糠炊斧笛挫欲夯吓胖拘圾寂席盈捆兑贬萄鞋瓢慷弦眼科医院病床安排含程序代码莉靶削枫舵咸枣力哈厂事兆坪巧俺揣奠著的拘邮炯痈千诵掌典咬陨昭界帖儒棘涡择韶反宗丝恃托轰丈挣般榴董推时胀徊迁卧胚埠哇杨俺缩萧敷么篱捆沦孜砍刚圆张雍彦猪傲丛碎郡蝗四漱刺吠环罐撬盒智昨辫官唯睁源德造矣批荷擦龄锹瓦徊悬蹦诀单弗幕签打坞清座嫡砚窖叭交射堪家龙劳偷噪瘟蠕曲舱疙盐胚芹膘谭蜀所贫据较瑚搔正约心抚淖斤族娟柳泥钓蔬睹巧广契盐唾替拙必精蠕棺辉怖傈寓椒邓详央笔淄辫层毡稍圃彩储多实民棋和纫聘椭潮替埔昧扛

4、侥玛丰介劝淀之战秘城扰玄卖棠萝方蟹谋合槛逢芍凉缸细濒糊腮翘蓟助件诛奏辱炳槐革味拼水谴墙诧茸褥发指密荣合证袱颧侮妨垣阮高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将

5、受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：眼科病床的合理安排摘要针对不同的问题，本文分别建立了评价模型、优

6、先病床安排模型、排队论模型，利用了统计方法，计算机模拟方法和各种数学软件来进行计算。 问题一：将公平性指标、病床利用率和医院工作效率一次量化为：平均入院等待时间、手术准备时间和术后恢复时间，以这三个指标来评价具体模型，运用层次分析法确定这三个指标的权重，并进行一致性处理。得到综合评价指标为：。A越小，说明病床安排越合理。对于原始的病床安排方案，综合评价指标A=5.40。问题二：根据实际情况，制定优先病床安排规则。依据每一类病人的规律，在新的优先规则下，通过计算机模拟，将2008-8-30到2008-9-11已经门诊的病人入院情况进行合理安排。用问题一中的综合评价指标来评价这一模型，综合评价指标

7、A=4.24n？安排此人住院i+此时刻星期几？根据相应的规则安排住院i+im？确定我们要安排病人是否全部入院？x+X+终止NYNYYNYN4.4模型求解利用MATLAB编程，随机模拟出每天出院的各类病人情况并统计出出院人数（代码见附录2）。并按照上面我们制定的规则安排2008-8-30到2008-9-11已经门诊的病人入院，见下表2。在此，我们已经完成了根据第二天的出院人数来确定应该安排那些病人入院。序号类型门诊时间入院时间第一次手术时间第二次手术时间出院时间 入院等待时间手术等待时间1白内障(双眼)2008-8-302008-9-132008-9-152008-9-17/1422视网膜疾病2

8、008-8-302008-9-122008-9-141323青光眼2008-8-302008-9-122008-9-141324视网膜疾病2008-8-302008-9-122008-9-141325视网膜疾病2008-8-302008-9-122008-9-141326白内障(双眼)2008-8-302008-9-132008-9-152008-9-17/1427白内障2008-8-312008-9-142008-9-17/1438青光眼2008-8-312008-9-122008-9-141229白内障(双眼)2008-8-312008-9-132008-9-152008-9-17/132

9、10视网膜疾病2008-8-312008-9-122008-9-1412211视网膜疾病2008-8-312008-9-122008-9-1412212视网膜疾病2008-8-312008-9-152008-9-1715213青光眼2008-8-312008-9-162008-9-1816214白内障2008-8-312008-9-142008-9-17/14315视网膜疾病2008-9-12008-9-162008-9-1815216视网膜疾病2008-9-12008-9-162008-9-1815217青光眼2008-9-12008-9-162008-9-1815218白内障(双眼)200

10、8-9-12008-9-132008-9-152008-9-17/12219白内障(双眼)2008-9-12008-9-132008-9-152008-9-17/12220白内障(双眼)2008-9-12008-9-132008-9-152008-9-17/12221视网膜疾病2008-9-12008-9-162008-9-1815222白内障2008-9-12008-9-142008-9-15/13123视网膜疾病2008-9-12008-9-162008-9-1815224视网膜疾病2008-9-12008-9-162008-9-1815225白内障2008-9-22008-9-14200

11、8-9-15/12126白内障2008-9-22008-9-152008-9-17/13227白内障(双眼)2008-9-22008-9-132008-9-152008-9-17/11228白内障2008-9-22008-9-152008-9-17/13229视网膜疾病2008-9-22008-9-162008-9-1814230视网膜疾病2008-9-32008-9-162008-9-1813231视网膜疾病2008-9-32008-9-162008-9-1813232白内障(双眼)2008-9-32008-9-132008-9-152008-9-17/10233白内障2008-9-3200

12、8-9-152008-9-17/12234视网膜疾病2008-9-32008-9-162008-9-1813235白内障2008-9-32008-9-152008-9-17/12236视网膜疾病2008-9-32008-9-162008-9-1813237视网膜疾病2008-9-32008-9-162008-9-1813238白内障(双眼)2008-9-42008-9-132008-9-152008-9-17/9239白内障2008-9-42008-9-152008-9-17/11240青光眼2008-9-42008-9-162008-9-1812241视网膜疾病2008-9-42008-9-

13、162008-9-1812242视网膜疾病2008-9-42008-9-162008-9-1812243视网膜疾病2008-9-42008-9-162008-9-1812244青光眼2008-9-42008-9-162008-9-1812245白内障(双眼)2008-9-42008-9-132008-9-152008-9-17/9246白内障(双眼)2008-9-42008-9-132008-9-152008-9-17/9247青光眼2008-9-42008-9-162008-9-1812248青光眼2008-9-42008-9-162008-9-1812249视网膜疾病2008-9-4200

14、8-9-162008-9-1812250视网膜疾病2008-9-42008-9-162008-9-1812251白内障(双眼)2008-9-52008-9-132008-9-152008-9-17/8252白内障(双眼)2008-9-52008-9-132008-9-152008-9-17/8253白内障(双眼)2008-9-52008-9-142008-9-152008-9-17/9154视网膜疾病2008-9-52008-9-162008-9-1811255白内障(双眼)2008-9-52008-9-142008-9-152008-9-17/9156青光眼2008-9-52008-9-16

15、2008-9-1811257白内障(双眼)2008-9-52008-9-142008-9-152008-9-17/9158白内障2008-9-52008-9-152008-9-17/10259白内障(双眼)2008-9-52008-9-142008-9-152008-9-17/9160白内障(双眼)2008-9-52008-9-142008-9-152008-9-17/9161白内障(双眼)2008-9-62008-9-142008-9-152008-9-17/8162视网膜疾病2008-9-62008-9-162008-9-1810263青光眼2008-9-62008-9-172008-9-

16、1911264白内障(双眼)2008-9-62008-9-142008-9-152008-9-17/8165视网膜疾病2008-9-72008-9-172008-9-1910266白内障(双眼)2008-9-72008-9-142008-9-152008-9-17/7167视网膜疾病2008-9-72008-9-172008-9-1910268白内障2008-9-82008-9-152008-9-17/7269视网膜疾病2008-9-82008-9-172008-9-199270视网膜疾病2008-9-82008-9-172008-9-199271白内障2008-9-82008-9-15200

17、8-9-17/7272白内障(双眼)2008-9-82008-9-142008-9-152008-9-17/6173白内障2008-9-82008-9-152008-9-17/7274视网膜疾病2008-9-82008-9-172008-9-199275白内障2008-9-82008-9-152008-9-17/7276青光眼2008-9-92008-9-172008-9-198277青光眼2008-9-92008-9-172008-9-198278视网膜疾病2008-9-92008-9-172008-9-198279白内障2008-9-92008-9-152008-9-17/6280白内障2

18、008-9-92008-9-152008-9-17/6281视网膜疾病2008-9-102008-9-172008-9-197282白内障2008-9-102008-9-152008-9-17/5283白内障(双眼)2008-9-102008-9-142008-9-152008-9-17/4184白内障2008-9-102008-9-152008-9-17/5285白内障2008-9-102008-9-152008-9-17/5286白内障(双眼)2008-9-102008-9-142008-9-152008-9-17/4187白内障2008-9-102008-9-152008-9-17/52

19、88青光眼2008-9-102008-9-172008-9-197289白内障(双眼)2008-9-102008-9-142008-9-152008-9-17/4190视网膜疾病2008-9-112008-9-172008-9-196291视网膜疾病2008-9-112008-9-172008-9-196292青光眼2008-9-112008-9-172008-9-196293白内障(双眼)2008-9-112008-9-142008-9-152008-9-17/3194白内障(双眼)2008-9-112008-9-142008-9-152008-9-17/3195青光眼2008-9-1120

20、08-9-172008-9-196296白内障(双眼)2008-9-112008-9-142008-9-152008-9-17/3197外伤2008-9-112008-9-122008-9-131198白内障(双眼)2008-9-112008-9-142008-9-152008-9-17/3199视网膜疾病2008-9-112008-9-172008-9-1962100白内障2008-9-112008-9-152008-9-17/42101视网膜疾病2008-9-112008-9-172008-9-1962102视网膜疾病2008-9-112008-9-172008-9-1962表二4.5模型

21、评价对4.4中求解的结果进行统计，得到=9.75天，=1.83天，病床安排模型的改变并不会影响病人的术后恢复时间，因此=6.64天不变，这三个数据依然是不考虑外伤在内的数据。此改进病床安排方案的综合评价指标为：平均入院等待时间=9.7512.32,公平性得到了提高，=1.832.71,病床利用率也得到了提高。可以看出总和评价指标4.240,否则0,这是不可能的，所以，可以计算得到：此时，还剩7张床位可供分配，共有=78125种分配方式。对所有的分配方式进行穷举可以得到最优解为：即为白内障单眼安排9张病床，为白内障双眼安排19张病床，为青光眼安排14张病床，为视网膜疾病安排25张病床，为外伤安排

22、12张病床可实现所有病人在系统内的平均逗留时间（含等待入院及住院时间）最短，。模型评价模型优点：1.在对原始数据进行理解和分析的基础之上，又对数据做了科学的处理，提高了数据的直观性和可操作性。2.模型的建立过程中充分考虑了医院资源的利用率，患者获得最优服务等因素，使模型与实际紧密联系。3.模型中采用了排队论、计算机仿真等，有着较强的较雄厚的理论基础和真实性。模型缺点：模型虽然综合考虑了其他因素，但是为了建模，理想化了许多影响因素，具有一定的局限性，得到的最优方案与实际有一定的出入。模型改进 在第一问成对比较矩阵的确定过程中，我们是通过自身的生活经验及查询了少量数据便给出了成对比较矩阵。显然该方

23、法存在一定的主观性和随机性，从而对后面的模型评价产生了一定的影响。可以通过查询大量数据或组织专家评估来获取更权威的数据参考文献【1】钱颂迪主编. 运筹学（第三版）.清华大学出版社，2005年【2】王玉升，排队论模型及其在医院管理中的作用,中国医院管理, 58-62，1985.2【3】宋来忠，王志明,数学建模与实验,北京：科学出版社,2005。附录： 1 2008-07-13到2008-09-11的病人信息（略）2. 蒙特卡洛方法计算每天出院的各类病人情况并统计出出院人数的程序function f=max_1(x)n=length(x)max=x(1);for i=2:n if maxx(i)

24、max=x(i); endendf=maxfunction f=max_2(x)n=length(x)m=1max=x(1);for i=2:n if maxx(i) max=x(i); m=i; endendf=m;n=20;offef= x1= 9 9 11 11 16 16 16; x2=99999999999999916161616161616; x3=5 8 8 10 10 10 12 13 13 ; x4=1155566777778810101010101010101212121213131313131414; x5=7 8 8 8 9 11 12 12 ; n1=length(x

25、1);n2=length(x2);n3=length(x3);n4=length(x4);n5=length(x5);for i=1:n1 y(i)=x1(i)+ceil(unidrnd(3)+1);endfor i=n1+1:n1+n2 y(i)=x2(i-n1)+ceil(unidrnd(3)+3);endfor i=n1+n2+1:n1+n2+n3 y(i)=x3(i-n1-n2)+ceil(unidrnd(9)+3);endfor i=n1+n2+n3+1:n1+n2+n3+n4 y(i)=x4(i-n1-n2-n3)+ceil(unidrnd(11)+4);endfor i=n1+n

26、2+n3+n4+1:n1+n2+n3+n4+n5 y(i)=x5(i-n1-n2-n3-n4)+ceil(unidrnd(3)+1);enda=1123334456999910101111111112 offef;b=001222345556666666778911111112121212 offef;c=000011111222222344444555555555667788999101010111112121212121212 offefd=12 offefq1=length(a)-1;q2=length(b)-1;q3=length(c)-1;q4=length(d)-1;p1=0;p2

27、=0;p3=0;p4=0;x=12;while 1 if p1=q1&p2=q2&p3=q3&p4=q4 break; end m=0; sort(y); for i=1:79 if y(i)n r(x-11,i)=maxn if maxn=1 p1=p1+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; end if maxn=2 p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; end if maxn=3 p3=p3+1; y(i)=x+ 2+(ceil(unid

28、rnd(9)+3)+ceil(unidrnd(11)+4)/2 end continue; end if mod(x,7)=6|mod(x,7)=0 if p2=q2 r(x-11,i)=2; p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; continue; elseif p1=q1 r(x-11,i)=1; p1=p1+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; continue; elseif p3=q3 r(x-11,i)=3; p3=p3+1; y(

29、i)=x+2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 continue; end end if mod(x,7)=1|mod(x,7)=2 if p1=q1 r(x-11,i)=1; p1=p1+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; continue; elseif p3=q3 r(x-11,i)=3; p3=p3+1; y(i)=x+2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 continue; elseif p2=q2 r(x-11,

30、i)=2; p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; continue; end end if p3=q3 r(x-11,i)=3; p3=p3+1; y(i)=x+2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 continue; elseif p2=q2 r(x-11,i)=2; p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; continue; elseif p1=q1 r(x-11,i)=1; p1

31、=p1+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; continue; end end x=x+1; end3问题三代码x=7 15 20 9 6 2 6 9 8 13 6 4 7 10 4 8 18 11 6 3 8 12 10 14 2 6 2 5 9 13 17 10 4 5 13 7;muhat,sigmahat,muci,sigmaci = normfit(x,0.05)muhat = 8.5833sigmahat = 4.5818muci = 7.0331 10.1336sigmaci = 3.71625.97674.

32、问题三检验的代码function f=sail_13(n)s=0;i=0;a=1 ;%每天外伤人数while 1 if s=n break; end s=s+ceil( normrnd(8.5676,4.51873)-a; i=i+1;endf=i;for i=1:100 x(i)=sail_13(90); end5.问题四代码function f=max_1(x)n=length(x)max=x(1);for i=2:n if maxx(i) max=x(i); endendf=maxfunction f=max_2(x)n=length(x)m=1max=x(1);for i=2:n if

33、 maxx(i) max=x(i); m=i; endendf=m;n=20;offef= x1= 9 9 11 11 16 16 16; x2=99999999999999916161616161616; x3=5 8 8 10 10 10 12 13 13 ; x4=1155566777778810101010101010101212121213131313131414; x5=7 8 8 8 9 11 12 12 ; n1=length(x1);n2=length(x2);n3=length(x3);n4=length(x4);n5=length(x5);for i=1:n1 y(i)=

34、x1(i)+ceil(unidrnd(3)+1);endfor i=n1+1:n1+n2 y(i)=x2(i-n1)+ceil(unidrnd(3)+3);endfor i=n1+n2+1:n1+n2+n3 y(i)=x3(i-n1-n2)+ceil(unidrnd(9)+3);endfor i=n1+n2+n3+1:n1+n2+n3+n4 y(i)=x4(i-n1-n2-n3)+ceil(unidrnd(11)+4);endfor i=n1+n2+n3+n4+1:n1+n2+n3+n4+n5 y(i)=x5(i-n1-n2-n3-n4)+ceil(unidrnd(3)+1);enda=112

35、3334456999910101111111112 offef;b=001222345556666666778911111112121212 offef;c=000011111222222344444555555555667788999101010111112121212121212 offefd=12 offefq1=length(a)-1;q2=length(b)-1;q3=length(c)-1;q4=length(d)-1;p1=0;p2=0;p3=0;p4=0;x=12;while 1 if p1=q1&p2=q2&p3=q3&p4=q4 break; end m=0; sort(y

36、); for i=1:79 if y(i)n r(x-11,i)=maxn if maxn=1 p1=p1+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; end if maxn=2 p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; end if maxn=3 p3=p3+1; y(i)=x+ 2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 end continue; end if mod(x,7)=6|mod(x,7)=0

37、if p2=q2 r(x-11,i)=2; p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; continue; elseif p1=q1 r(x-11,i)=1; p1=p1+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; continue; elseif p3=q3 r(x-11,i)=3; p3=p3+1; y(i)=x+2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 continue; end end if mod(x

38、,7)=1|mod(x,7)=2 if p1=q1 r(x-11,i)=1; p1=p1+1; y(i)=x+2+(ceil(unidrnd(3)+1)+ceil(unidrnd(3)+3)/2; continue; elseif p3=q3 r(x-11,i)=3; p3=p3+1; y(i)=x+2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 continue; elseif p2=q2 r(x-11,i)=2; p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; continue; end end if p3=q3 r(x-11,i)=3; p3=p3+1; y(i)=x+2+(ceil(unidrnd(9)+3)+ceil(unidrnd(11)+4)/2 continue; elseif p2=q2 r(x-11,i)=2; p2=p2+1; y(i)=x+2+(ceil(unidrnd(3)+3)+ceil(unidrnd(9)+3)/2; continue; elseif