2009高教社杯全国大学生数学建模竞赛(6)

72 白内障(双眼) 9月20日，9月21日 73 白内障 9月15日，9月16日，9月21日 74 视网膜疾病 9月19日，9月24日 75 白内障 9月15日，9月16日，9月20日，9月21日，9月22日 76 青光眼 9月19日，9月24日 77 青光眼 9月19日，9月24日 78 视网膜疾病 9月19日，9月24日 79 白内障 9月15日，9月16日，9月20日，9月21日，9月22日 80 白内障 9月15日，9月16日，9月20日，9月21日，9月22日 81 视网膜疾病 9月19日，9月24日 82 白内障 9月16日，9月20日，9月21日，9月22日 83 白内障(双眼) 9月20日，9月21日，9月22日 84 白内障 9月16日，9月20日，9月21日，9月22日 85 白内障 9月16日，9月21日，9月22日，9月23日 86 白内障(双眼) 9月20日，9月21日，9月27日 87 白内障 9月16日，9月21日，9月22日，9月23日 88 青光眼 9月19日，9月24日 89 白内障(双眼) 9月20日，9月21日，9月27日 90 视网膜疾病 9月19日，9月24日，9月25日 91 视网膜疾病 9月19日，9月24日，9月25日 92 青光眼 9月19日，9月24日，9月25日 93 白内障(双眼) 9月20日，9月21日，9月27日 94 白内障(双眼) 9月20日，9月21日，9月27日 95 青光眼 9月19日，9月24日，9月25日 96 白内障(双眼) 9月20日，9月21日，9月27日 97 外伤 9月12日 98 白内障(双眼) 9月21日，9月27日 99 视网膜疾病 9月19日，9月24日，9月25日，9月26日 100 白内障 9月16日，9月21日，9月22日，9月23日，9月28日 101 视网膜疾病 9月19日，9月24日，9月25日，9月26日 102 视网膜疾病 9月24日，9月25日，9月26日 附录五：模拟所用的程序

（1）函数M 文件jsdata.m 用来计算连个日期之间的天数

function f=jsdata(data1,data2)

if str2num(data1(6))-str2num(data2(6))<0

f1=31*(str2num(data2(6))-str2num(data1(6))); else f1=0; end

n1=str2num(data1(8:end)); n2=str2num(data2(8:end)); f=n2-n1+f1;

（2）模拟的程序： clear;clc

26

load t1,load t2,load t3;

%t1 t2 t3分别是给出的三部分的数据 %时间统一化处理 for i=1:size(t1,1)

t1{i,2}=jsdata('2008-7-12',t1{i,2}); t1{i,3}=jsdata('2008-7-12',t1{i,3}); t1{i,4}=jsdata('2008-7-12',t1{i,4}); t1{i,6}=jsdata('2008-7-12',t1{i,6}); if ~strcmp(t1{i,5},'/')

t1{i,5}=jsdata('2008-7-12',t1{i,5}); end end

for i=1:size(t2,1)

t2{i,2}=jsdata('2008-7-12',t2{i,2}); t2{i,3}=jsdata('2008-7-12',t2{i,3}); t2{i,4}=jsdata('2008-7-12',t2{i,4}); if ~strcmp(t2{i,5},'/')

t2{i,5}=jsdata('2008-7-12',t2{i,5}); end end

for i=1:size(t3,1)

t3{i,2}=jsdata('2008-7-12',t3{i,2}); end

%计算机模拟产生对于t2中的79个人的出院情况

%已知到第61天都还没有出院，所以最早的出院时间是第62天 for i=1:size(t2,1)

if strcmp(t2{i,1},'白内障') time=unidrnd(3)+1; t2{i,6}=t2{i,4}+time;

elseif strcmp(t2{i,1},'白内障(双眼)') time=unidrnd(3)+3; t2{i,6}=t2{i,4}+time;

elseif strcmp(t2{i,1},'青光眼') time=unidrnd(9)+3; t2{i,6}=t2{i,4}+time;

elseif strcmp(t2{i,1},'视网膜疾病') time=unidrnd(11)+4; t2{i,6}=t2{i,4}+time;

elseif strcmp(t2{i,1},'外伤') time=unidrnd(8)+2; t2{i,6}=t2{i,4}+time; end

if t2{i,6}<62 t2{i,6}=62;

27

end end

for i=1:size(t1,1)

t1{i,7}=mod(t1{i,2}-1,7); t1{i,8}=mod(t1{i,3}-1,7); t1{i,9}=mod(t1{i,4}-1,7); t1{i,10}=mod(t1{i,6}-1,7); t1{i,11}=11; end

for i=1:size(t2,1)

t2{i,7}=mod(t2{i,2}-1,7); t2{i,8}=mod(t2{i,3}-1,7); t2{i,9}=mod(t2{i,4}-1,7); t2{i,10}=mod(t2{i,6}-1,7); t2{i,11}=11; end

for i=1:size(t3,1)

t3{i,7}=mod(t3{i,2}-1,7); t3{i,8}=-1; t3{i,9}=-1; t3{i,10}=-1; t3{i,11}=11; end

tt=t2;%初始的序列；

k=1;in={};inn=[];long=102;lon=[]; f=62;%9月12号的序号，以后以此递增； while k<=30

index=[];%存储需要记录的下标号。 n=0;%统计每天调度时的可用的病床数

for i=1:size(tt,1); if tt{i,6}==f n=n+1; end end

long=long-n;

for i=1:size(t3,1)%首先安排伤的病人

if strcmp(t3{i,1},'外伤')&t3{i,11}~=0&n>0 index=[index,i]; t3{i,11}=0; t3{i,3}=f; t3{i,4}=f+1;

t3{i,6}=f+1+unidrnd(8)+2; t3{i,8}=mod(t3{i,3}-1,7);

28

t3{i,9}=mod(t3{i,4}-1,7); t3{i,10}=mod(t3{i,6}-1,7); n=n-1; end end

if mod(f-1,7)==2%这天是周二时：

for i=1:size(t3,1)%此时确定是安排白内障还是（青光眼、视网膜） if strcmp(t3{i,1},'白内障')&t3{i,11}~=0&n>0 index=[index,i]; t3{i,11}=0; t3{i,3}=f; t3{i,4}=f+1;

t3{i,6}=f+1+unidrnd(3)+1; t3{i,8}=mod(t3{i,3}-1,7); t3{i,9}=mod(t3{i,4}-1,7); t3{i,10}=mod(t3{i,6}-1,7); n=n-1; end end

elseif mod(f-1,7)==1%这天是周一时：

for i=1:size(t3,1)%此时确定是安排白内障还是（青光眼、视网膜） if strcmp(t3{i,1},'白内障')&t3{i,11}~=0&n>0 index=[index,i]; t3{i,11}=0; t3{i,3}=f; t3{i,4}=f+2;

t3{i,6}=f+2+unidrnd(3)+1; t3{i,8}=mod(t3{i,3}-1,7); t3{i,9}=mod(t3{i,4}-1,7); t3{i,10}=mod(t3{i,6}-1,7); n=n-1; end end

elseif mod(f-1,7)==0%这天是周日时：

for i=1:size(t3,1)%此时确定是安排白内障还是（青光眼、视网膜）

if (strcmp(t3{i,1},'白内障')|strcmp(t3{i,1},'白内障(双眼)'))&t3{i,11}~=0&n>0

index=[index,i]; t3{i,11}=0; t3{i,3}=f; t3{i,4}=f+1;

if strcmp(t3{i,1},'白内障') t3{i,6}=f+1+unidrnd(3)+1; else

29

t3{i,6}=f+1+unidrnd(3)+3; end

t3{i,8}=mod(t3{i,3}-1,7); t3{i,9}=mod(t3{i,4}-1,7); t3{i,10}=mod(t3{i,6}-1,7); n=n-1; end end

elseif mod(f-1,7)==6%这天是周六时：

for i=1:size(t3,1)%此时确定是安排白内障还是（青光眼、视网膜）

if (strcmp(t3{i,1},'白内障')|strcmp(t3{i,1},'白内障(双眼)'))&t3{i,11}~=0&n>0

index=[index,i]; t3{i,11}=0; t3{i,3}=f; t3{i,4}=f+2;

if strcmp(t3{i,1},'白内障') t3{i,6}=f+2+unidrnd(3)+1; else

t3{i,6}=f+2+unidrnd(3)+3; end

t3{i,8}=mod(t3{i,3}-1,7); t3{i,9}=mod(t3{i,4}-1,7); t3{i,10}=mod(t3{i,6}-1,7); n=n-1; end end

elseif mod(f-1+2,7)~=1&mod(f-1+2,7)~=3 for i=1:size(t3,1)

if (strcmp(t3{i,1},'青光眼')|strcmp(t3{i,1},'视网膜疾病'))&t3{i,11}~=0&n>0

index=[index,i]; t3{i,11}=0; t3{i,3}=f; t3{i,4}=f+2;

if strcmp(t3{i,1},'青光眼') t3{i,6}=f+2+unidrnd(9)+3; else

t3{i,6}=f+2+unidrnd(11)+4; end

t3{i,8}=mod(t3{i,3}-1,7); t3{i,9}=mod(t3{i,4}-1,7); t3{i,10}=mod(t3{i,6}-1,7); n=n-1;

30

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库2009高教社杯全国大学生数学建模竞赛(6)在线全文阅读。