用matlab写出MM1和MMK的伪代码

MM1伪代码：

//初始化
ARRIVAL_TIME = 0;
SERVICE_TIME = 0;
DEPARTURE_TIME = 0;
QUEUE_SIZE = 0;
TOTAL_WAIT_TIME = 0;
TOTAL_CUSTOMER_SERVED = 0;
SIMULATION_TIME = 10000;

//循环模拟
while ARRIVAL_TIME < SIMULATION_TIME
    if QUEUE_SIZE == 0
        //队列为空，顾客直接服务
        ARRIVAL_TIME = ARRIVAL_TIME + exprand(lambda); //根据到达率生成顾客到达时间
        SERVICE_TIME = exprand(mu); //根据服务率生成服务时间
        DEPARTURE_TIME = ARRIVAL_TIME + SERVICE_TIME; //计算离开时间
        TOTAL_WAIT_TIME = TOTAL_WAIT_TIME + 0; //等待时间为0
        TOTAL_CUSTOMER_SERVED = TOTAL_CUSTOMER_SERVED + 1; //累计服务顾客数量
    else
        //队列不为空，顾客排队等待服务
        ARRIVAL_TIME = ARRIVAL_TIME + exprand(lambda); //根据到达率生成顾客到达时间
        QUEUE_SIZE = QUEUE_SIZE + 1; //队列长度+1
    end
    if ARRIVAL_TIME >= DEPARTURE_TIME
        //顾客已经离开，可以开始服务下一个顾客
        QUEUE_SIZE = QUEUE_SIZE - 1; //队列长度-1
        SERVICE_TIME = exprand(mu); //根据服务率生成服务时间
        TOTAL_WAIT_TIME = TOTAL_WAIT_TIME + (DEPARTURE_TIME - ARRIVAL_TIME); //计算等待时间
        DEPARTURE_TIME = DEPARTURE_TIME + SERVICE_TIME; //计算离开时间
        TOTAL_CUSTOMER_SERVED = TOTAL_CUSTOMER_SERVED + 1; //累计服务顾客数量
    end
end

//输出结果
AVG_WAIT_TIME = TOTAL_WAIT_TIME / TOTAL_CUSTOMER_SERVED;
AVG_SYSTEM_TIME = AVG_WAIT_TIME + 1 / mu;

MMK伪代码：

//初始化
ARRIVAL_TIME = 0;
SERVICE_TIME = zeros(1, K);
DEPARTURE_TIME = zeros(1, K);
QUEUE_SIZE = 0;
TOTAL_WAIT_TIME = 0;
TOTAL_CUSTOMER_SERVED = 0;
SIMULATION_TIME = 10000;

//循环模拟
while ARRIVAL_TIME < SIMULATION_TIME
    if QUEUE_SIZE < K
        //队列长度小于K，顾客直接服务
        ARRIVAL_TIME = ARRIVAL_TIME + exprand(lambda); //根据到达率生成顾客到达时间
        SERVICE_TIME(QUEUE_SIZE+1) = exprand(mu); //根据服务率生成服务时间
        DEPARTURE_TIME(QUEUE_SIZE+1) = ARRIVAL_TIME + SERVICE_TIME(QUEUE_SIZE+1); //计算离开时间
        TOTAL_WAIT_TIME = TOTAL_WAIT_TIME + 0; //等待时间为0
        TOTAL_CUSTOMER_SERVED = TOTAL_CUSTOMER_SERVED + 1; //累计服务顾客数量
        QUEUE_SIZE = QUEUE_SIZE + 1; //队列长度+1
    else
        //队列长度大于等于K，顾客排队等待服务
        ARRIVAL_TIME = ARRIVAL_TIME + exprand(lambda); //根据到达率生成顾客到达时间
        QUEUE_SIZE = QUEUE_SIZE + 1; //队列长度+1
    end
    for i = 1:K
        if ARRIVAL_TIME >= DEPARTURE_TIME(i)
            //第i个顾客已经离开，可以开始服务下一个顾客
            QUEUE_SIZE = QUEUE_SIZE - 1; //队列长度-1
            TOTAL_WAIT_TIME = TOTAL_WAIT_TIME + (DEPARTURE_TIME(i) - ARRIVAL_TIME); //计算等待时间
            SERVICE_TIME(i) = exprand(mu); //根据服务率生成服务时间
            DEPARTURE_TIME(i) = DEPARTURE_TIME(i) + SERVICE_TIME(i); //计算离开时间
            TOTAL_CUSTOMER_SERVED = TOTAL_CUSTOMER_SERVED + 1; //累计服务顾客数量
        end
    end
end

//输出结果
AVG_WAIT_TIME = TOTAL_WAIT_TIME / TOTAL_CUSTOMER_SERVED;
AVG_SYSTEM_TIME = AVG_WAIT_TIME + 1 / mu;