MM1、MMk排队系统的仿真和比较实验matlab

下面是一个 MM1、MMk 排队系统的仿真和比较实验MATLAB代码示例，你可以参考一下：

% 设置仿真参数
lambda = 4; % 平均到达率
mu = 5; % 平均服务率
K = 3; % MMK系统的服务员数

% 初始化仿真状态（MM1系统）
queue_mm1 = 0; % 初始队列长度为0
arrival_time_mm1 = exprnd(1/lambda); % 第一个顾客到达时间
service_time_mm1 = exprnd(1/mu); % 第一个顾客服务时间
departure_time_mm1 = arrival_time_mm1 + service_time_mm1; % 第一个顾客离开时间

% 初始化仿真状态（MMK系统）
queue_mmk = 0; % 初始队列长度为0
arrival_time_mmk = exprnd(1/lambda); % 第一个顾客到达时间
service_times_mmk = zeros(1, K); % 所有服务员的服务时间
departure_times_mmk = zeros(1, K); % 所有服务员的离开时间
for i = 1:K % 初始化所有服务员的状态
    service_times_mmk(i) = exprnd(1/mu);
    departure_times_mmk(i) = arrival_time_mmk + service_times_mmk(i);
end

% 运行仿真（MM1系统）
num_customers = 10000; % 定义仿真的客户数量
for i = 2:num_customers % 第一个顾客已经被初始化，从第二个顾客开始仿真
    % 客户到达
    interarrival_time = exprnd(1/lambda);
    arrival_time_mm1(i) = arrival_time_mm1(i-1) + interarrival_time;
    % 客户进入队列
    queue_mm1(i) = queue_mm1(i-1) + 1;
    % 如果服务员空闲，则开始服务
    if queue_mm1(i) == 1
        service_time_mm1(i) = exprnd(1/mu);
        departure_time_mm1(i) = arrival_time_mm1(i) + service_time_mm1(i);
    else % 如果服务员忙碌，则继续等待
        service_time_mm1(i) = 0;
        departure_time_mm1(i) = departure_time_mm1(i-1) + exprnd(1/mu);
    end
    % 客户离开
    if queue_mm1(i) > 0
        queue_mm1(i) = queue_mm1(i) - 1;
    end
end

% 运行仿真（MMK系统）
for i = K+1:num_customers % 第一个顾客已经被初始化，从第 K+1 个顾客开始仿真
    % 客户到达
    interarrival_time = exprnd(1/lambda);
    arrival_time_mmk(i) = arrival_time_mmk(i-1) + interarrival_time;
    % 客户进入队列
    queue_mmk(i) = queue_mmk(i-1) + 1;
    % 找到空闲的服务员
    idle_servers = find(departure_times_mmk <= arrival_time_mmk(i));
    if isempty(idle_servers) % 如果没有空闲的服务员，则客户继续等待
        service_times_mmk(i) = 0;
        departure_times_mmk(i) = max(departure_times_mmk) + exprnd(1/mu);
    else % 如果有空闲的服务员，则客户开始服务
        server = idle_servers(1); % 选择第一个空闲的服务员
        service_times_mmk(i) = exprnd(1/mu);
        departure_times_mmk(i) = arrival_time_mmk(i) + service_times_mmk(i);
    end
    % 客户离开
    if queue_mmk(i) > 0
        queue_mmk(i) = queue_mmk(i) - 1;
    end
end

% 绘制仿真结果（MM1系统）
subplot(2,2,1);
stairs(arrival_time_mm1, queue_mm1, 'b');
xlabel('时间');
ylabel('队列长度');
title('MM1系统：队列长度随时间变化');

subplot(2,2,2);
stairs([0, departure_time_mm1], [0, queue_mm1], 'r');
xlabel('时间');
ylabel('队列长度');
title('MM1系统：队列长度随时间变化');

% 绘制仿真结果（MMK系统）
subplot(2,2,3);
stairs(arrival_time_mmk, queue_mmk, 'b');
xlabel('时间');
ylabel('队列长度');
title('MMK系统：队列长度随时间变化');

subplot(2,2,4);
stairs([0, max(departure_times_mmk)], [0, queue_mmk(end)], 'r');
xlabel('时间');
ylabel('队列长度');
title('MMK系统：队列长度随时间变化');

这段代码模拟了一个 MM1、MMk 排队系统，其中包括了客户到达、服务员忙碌、客户离开等步骤，并绘制了队列长度随时间变化的图表。你可以根据自己的需求修改代码中的参数和绘图方式。