数学建模基于matlab单列多服务台排队系统仿真含Matlab源码 1698期

Posted 2022-02-17 紫极神光

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【数学建模】基于matlab单列多服务台排队系统仿真【含Matlab源码 1698期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

function varargout = MMc(varargin)
% MMC M-file for MMc.fig
%      MMC, by itself, creates a new MMC or raises the existing
%      singleton*.
%	
%      copy of xuwentao
%
%      H = MMC returns the handle to a new MMC or the handle to
%      the existing singleton*.
%
%      MMC('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in MMC.M with the given input arguments.
%
%      MMC('Property','Value',...) creates a new MMC or raises the
%      existing singleton*.  Starting from the left, property value pairs are
%      applied to the GUI before MMc_OpeningFcn gets called.  An
%      unrecognized property name or invalid value makes property application
%      stop.  All inputs are passed to MMc_OpeningFcn via varargin.
%
%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one
%      instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help MMc

% Last Modified by GUIDE v2.5 03-Nov-2021 23:20:51

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
    'gui_Singleton',  gui_Singleton, ...
    'gui_OpeningFcn', @MMc_OpeningFcn, ...
    'gui_OutputFcn',  @MMc_OutputFcn, ...
    'gui_LayoutFcn',  [] , ...
    'gui_Callback',   []);
if nargin && ischar(varargin1)
    gui_State.gui_Callback = str2func(varargin1);
end

if nargout
    [varargout1:nargout] = gui_mainfcn(gui_State, varargin:);
else
    gui_mainfcn(gui_State, varargin:);
end
% End initialization code - DO NOT EDIT


% --- Executes just before MMc is made visible.
function MMc_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
% varargin   command line arguments to MMc (see VARARGIN)

% Choose default command line output for MMc
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% UIWAIT makes MMc wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = MMc_OutputFcn(hObject, eventdata, handles)
% varargout  cell array for returning output args (see VARARGOUT);
% hObject    handle to figure
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout1 = handles.output;



function num_server_EDIT_Callback(hObject, eventdata, handles)
% hObject    handle to num_server_EDIT (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of num_server_EDIT as text
%        str2double(get(hObject,'String')) returns contents of num_server_EDIT as a double


% --- Executes during object creation, after setting all properties.
function num_server_EDIT_CreateFcn(hObject, eventdata, handles)
% hObject    handle to num_server_EDIT (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
%       See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
    set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in OK.
function OK_Callback(hObject, eventdata, handles)
% hObject    handle to OK (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
num_serverTemp= str2double(get(handles.num_server_EDIT,'String'))  ;

ST_Idle=0;%%%服务台为空闲
ST_Busy=1;%%%服务台为繁忙
ST_All_Busy=-1;%%%所有服务台为繁忙
EV_NULL=0;
EV_Arrive=1;
EV_Depart=2;
EV_LEN=3;
Q_LIMIT=str2double(get(handles.Q_LIMIT_EDIT,'String'))  ;                %%%%%%%%%%排队系统容量
Q_LIMIT=1e10;

time_arrival=[];                 %到达时刻
time_next_event=zeros(1,EV_LEN);
%仿真参数
num_events=EV_LEN-1;
num_serveCounter=num_serverTemp;                  %%%%M/M/m  m为服务台数量
mean_interarrival=str2double(get(handles.mean_interarrival_EDIT,'String'))  ;      %%%%%顾客平均到达间隔时间
mean_service=str2double(get(handles.mean_service_EDIT,'String'))  ;             %%%%%为顾客服务平均服务时间
num_peo_required=str2double(get(handles.num_delays_required_EDIT,'String'))  ;           %%%%%%%顾客来源总量

serv_time=zeros(1,num_peo_required);
outfile=fopen('MMc.txt','w');
fprintf(outfile, 'MMc多服务台排队仿真系统：\\r\\n');
fprintf(outfile, '平均到达间隔时间为%11.3f minutes\\r\\n',mean_interarrival);
fprintf(outfile, '平均服务时间为%16.3f minutes\\r\\n', mean_service);
fprintf(outfile, '服务台数为%20d\\r\\n', num_serveCounter);
fprintf(outfile, '顾客总数为%14d\\r\\n', num_peo_required);

%part1 初始化
sim_time=0.0;
server_status   =zeros(1,num_serveCounter);   %状态空闲
num_in_q        = 0;%队列人数
time_last_event = 0.0;%事件持续时间
num_custs_simulated  = 0;%%%已经完成仿真的顾客数量
total_of_delays    = 0.0;%%%%%%队列中顾客等待总的时间
total_of_time    = 0.0;%%%%%%系统中顾客逗留总时间
area_num_in_q      = 0.0;
area_server_status = 0.0;

%/* Initialize event list. 初始事件表，因为没有顾客到达接受服务，离开事件时刻初始为无穷大*/
time_next_event(EV_Arrive) = sim_time + exprnd(mean_interarrival);%%记录时间事件表下一到达事件时间
time_next_event(EV_Depart) = 1.0e+230;%%记录时间事件表下一离开事件时间

time_serveDepart=zeros(1,num_serveCounter);%%%记录每个服务台实体中的顾客，服务完毕后离开此服务台的时刻
ordinal_serveCounter_depart=0;
    
%part2  
while (num_custs_simulated < num_peo_required)
%/*决定下次发生什么事件 */
min_time_next_event = 1.0e+290;
%  找出time_serveDepart(1,:)最小值所对应的下标i赋给ordinal_serveCounter_depart，确定当前时刻第几个服务台中的顾客最先离开
min_time_serveDepart=1e290;
ordinal_serveCounter_depart=0;
for i=1:num_serveCounter%%%%%%% 也可用[val,pos]=min(time_serveDepart)代替，得出val和pos
    if(server_status(i)==1 && time_serveDepart(i)<min_time_serveDepart)
        min_time_serveDepart=time_serveDepart(i);
        ordinal_serveCounter_depart=i;
    end
end
time_next_event(2)=min_time_serveDepart;

%/* 确定下一最早发生事件的类型： 1代表Arrive，2代表Departure */ 
next_event_type = 0;
for i = 1: num_events
    if (time_next_event(i) < min_time_next_event)
        min_time_next_event = time_next_event(i);
        next_event_type     = i;
    end
end

%/* 事件列表是否为空. */
if (next_event_type == 0)
    % 空事件表，终止仿真
    fprintf(outfile, '\\r\\nEvent list empty at time %f', sim_time);
    exit(1);
end
%/* 事件表不空，推进仿真时钟到下一最早时刻*/
sim_time = min_time_next_event;
double time_since_last_event; %%%刚发生事件持续时间
time_since_last_event = sim_time - time_last_event;
time_last_event       = sim_time;
%/* 更新队列中的（人数*排队时间），为了计算平均队列长即队列中的平均顾客数/
area_num_in_q=area_num_in_q +  num_in_q * time_since_last_event;
%/* 为了计算平均队长即系统中的平均顾客数 */
for i=1:num_serveCounter
    area_server_status =area_server_status + server_status(i) * time_since_last_event;
end

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]李昕.MATLAB数学建模[M].清华大学出版社.2017
[2]王健,赵国生.MATLAB数学建模与仿真[M].清华大学出版社.2016
[3]余胜威.MATLAB数学建模经典案例实战[M].清华大学出版社.2015