【数学建模集训系列】眼科病床安排问题——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