【数学建模集训系列】眼科病床安排问题——FCFS

%function sec_problem
%*********************************************************
% ※第二问代码：FCFS
% ※时间：2011/8/31
%*********************************************************
clear,clc
%**************根据术后恢复满足的正态分布随机产生术后恢复时间*****************
miu=[2.90 4.96 10.16 8.07 6.03];%各个正态分布均值
sigm=[0.69 0.59 1.56 2.36 1.82];%方差
num_type=[7 22 9 33 8];%各类病人数
for i=1:1:5
    shuhou{i,1}=(fix(normrnd(miu(i),sigm(i),num_type(i),1)))';%产生随机数,5*1单元阵列
    shuhou{i,1};
end
%**************计算可以入院病人的优先级，安排入院****************************
cha=693960;%matlab和excel日期之差
out1_time_num = xlsread('筛选_数据', 'FCFS', 'M2:N28');%前一组出院时间
in_order_type_time=xlsread('筛选_数据', 'FCFS', 'A2:C103');%门诊时间
%------------------计算优先级矩阵------------------------------
% 说明：值越小表示优先级越高,第二项为手术
%       白内障（单眼） 白内障（双眼） 青光眼  视网膜疾病  外伤
% 周日       3,1         2,1         4,2       4,2      1,1
% 周一       2,2         4,2         3,3       3,3      1,1
% 周二       2,1         4,1         3,2       3,2      1,1
% 周三       4,5         3,5         2,2       2,2      1,1
% 周四       4,4         3,4         2,2       2,2      1,1
% 周五       4,3         3,3         2,2       2,2      1,1
% 周六       3,2         2,2         4,3       4,3      1,1
%-------------------------------------------------------------
prio_mat=[3 2 4 4 1;2 4 3 3 1;2 4 3 3 1;4 3 2 2 1;4 3 2 2 1;4 3 2 2 1;3 2 4 4 1];
ssjg=    [1 1 2 2 1;2 2 3 3 1;1 1 2 2 1;5 5 2 2 1;4 4 2 2 1;3 3 2 2 1;2 2 3 3 1];
flag_visited=zeros(102,1);
in_order_type_time_flag=[in_order_type_time flag_visited];
ruyuan_mat=nan(102,1);
shoushu1_mat=nan(102,1);
shoushu2_mat=nan(102,1);
chuyuan_mat=nan(102,1);
tic
  for i=1:size(out1_time_num,1)
% for i=1:3
    to_day=out1_time_num(i,1);%当天excel时间
    m=find(in_order_type_time_flag(:,3)<to_day&in_order_type_time_flag(:,4)==0);
    if isempty(m)
       continue;
    else
        temp_in=in_order_type_time(m,:);%当天有机会入院的病人
        %判断当天星期几y=weekday(x)，返回值y：周日-周六分别为1-7
        to_week=weekday(to_day+cha);
        %附加优先级:序号+类型+门诊时间+优先级
        for j=1:size(temp_in,1)
            if temp_in(j,2)==5
                prio(j,1)=0;
            else
                prio(j,1)= temp_in(j,1);
            end
            prio(j+1:end)=[];
        end
        temp_in=[temp_in prio];
        %按优先级选取与病床相同数量的病人:按优先级排序，优先级相同的按门诊顺序时间排
        temp_in_sort=sortrows(temp_in,4);
        if size(out1_time_num,1)<=size(temp_in_sort,1)
            num_2=out1_time_num(i,2);        
        else
            num_2=size(temp_in_sort,1);
        end
        temp_write=temp_in_sort(1:num_2,:);
        %计算excel表格写入范围,逐个入院
        for j=1:num_2%out1_time_num(i,2)        
            ruyuan_mat(temp_write(j,1))=to_day;
            %write_str=['D' num2str(temp_write(j,1)+1)];
            %xlswrite('筛选_数据',to_day,'FCFS',write_str);    %入院
            %------------------------------------------------
            %安排手术时间
            shoushu1_mat(temp_write(j,1))=to_day+ssjg(to_week,temp_write(j,2));
            if temp_write(j,2)==2
                shoushu2_mat(temp_write(j,1))=to_day+ssjg(to_week,temp_write(j,2))+2;
            end
            %随机产生出院时间
            chuyuan_mat(temp_write(j,1))=to_day+ssjg(to_week,temp_write(j,2))+...
                                    fix(normrnd(miu(temp_write(j,2)),sigm(temp_write(j,2)),1));
            %更新病床数
            index=find(out1_time_num(:,1)==chuyuan_mat(temp_write(j,1)));
            out1_time_num(index,2)=out1_time_num(index,2)+1;
            %------------------------------------------------
        end
        %标记已经入院的数据
        in_order_type_time_flag(temp_write(:,1),:)=1;    
    end    
 end
%write_str='D2:D103';
%xlswrite('筛选_数据',ruyuan_mat,'FCFS',write_str);    %入院
disp('安排完成')
toc
write_mat=[ruyuan_mat shoushu1_mat shoushu2_mat chuyuan_mat]
write_str='D2:G103';
xlswrite('筛选_数据',write_mat,'FCFS',write_str);
out1_time_num