%% 参考陆传X《排队论》第二版[M].北京邮电大学出版社,2009,p72-76
%%=================初始化并赋初值================
clc;clear all
numda=xlsread('data1.xlsx','sheet1','b2:b12');
mu=xlsread('data1.xlsx','sheet1','c2:c12');
n=xlsread('data1.xlsx','sheet1','d2:d12');
m=xlsread('data1.xlsx','sheet1','e2:e12');
%Lqq=AA('sheet1','f2:f12')';
%Wqq=AA('sheet1','g2:g12')';
%numda=2.5392; %numda表示排队模型M/M/n/m中,产品到达速率
%mu=0.94; %mu表示工站中某一台设备的加工速率
%n=3; %n表示某个工站中并联加工设备的数量
%m=15; %m表示缓冲区大小(和ARENA中“NQ(xxx.Queue)<N”的关系为N=m-n,或者:m=N+n)
ro=numda./(n.*mu); %ro表示工站负荷水平或强度(ro=ro1/n)。
%============计算ro===========
if ro~=1
sum = 0;
for k = 0:n-1 % 循环变量根据需要改吧
sum = sum+((n.*ro).^k)/(factorial(k));
end
P0=(sum+((n.*ro).^n).*(1-(ro.^(m-n+1)))./(factorial(n).*(1-ro))).^(-1); %计算p0表示饥饿概率
Pm=((n.^n).*(ro.^m).*P0)./factorial(n); %pm表示缓冲区为m时的阻塞概率
wq=((((n.^n).*(ro.^(n+1)).*P0)./( factorial(n).*((1-ro).^2).*numda.*(1-Pm))).*(1-(m-n+1).*(ro.^(m-n))+(m-n).*(ro.^(m-n+1)))); %平均排队等待时间
p_loss=Pm; %(1)系统损失率,与上面的pm相同。
Q= 1-Pm; %(2)系统相对通过能力
numda1=numda.*Pm; %(3)单位时间平均损失的产品数量
numdae=numda.*(1-Pm); % (3)平均进入系统的产品数量
n_busy=numdae./mu; %(4)平均忙着的设备数量
Lq=((n.^n).*(ro.^(n+1)).*P0.*(1-(m-n+1).*(ro.^(m-n))+(m-n).*(ro.^(m-n+1))))./(factorial(n).*((1-ro).^2)); %(5)平均排队等待数量(排队等待队长)
Ls=Lq+(numdae./mu); %(6)平均队长(即平均产品数量)
Wq=Lq./numdae; %(7)平均排队等待时间
Ws=Wq+(1./mu); %(7)产品在系统内平均逗留时间
%===========================显示结果================================
%disp(['P0(饥饿概率/系统空闲率)=' num2str(P0)]);
% disp(['Q(系统相对通过能力)=' num2str(Q)]);
% disp(['numda1(单位时间平均损失的产品数量)=' num2str(numda1)]);
% disp(['numdae(平均进入系统的产品数量)=' num2str(numdae)]);
% disp(['Pm(缓冲区为m时的阻塞概率)=' num2str(Pm)]);
% disp(['n_busy(平均忙着的设备数量)=' num2str(n_busy')]);
% disp(['Ls(平均队长)=' num2str(Ls')]);
% disp(['Lq(平均等待队长)=' num2str(Lq')]);
% disp(['Ws(平均逗留时间)=' num2str(Ws')]);
else
sum = 0;
for k = 0:n-1 % 循环变量根据需要改吧
sum = sum + (n^k)/(factorial(k));
end
P0=(sum+(n^n)*(m-n+1))/(factorial(n))^(-1) %p0表示排队模型M/M/n/m中,饥饿率
end
%Lqerror=abs((Lqq-Lq))./Lqq;
%Wqerror=abs((Wqq-Wq))./Wqq;
%disp(['Lq(平均排队等待数量)=' num2str(Lq)])
% disp(['Lqerror(平均排队等待数量误差)=' num2str(Lqerror)])
% disp(['Wq(平均排队等待时间)=' num2str(Wq)])
%disp(['Wqerror(平均排队等待时间误差)=' num2str(Wqerror)])
xlswrite('data1.xlsx',Wq,'Sheet1','f2:f12');
xlswrite('data1.xlsx',Lq,'Sheet1','g2:g12');
xlswrite('data1.xlsx',ro,'Sheet1','h2:h12');
xlswrite('data1.xlsx',n_busy,'Sheet1','i2:i12');
WWq=xlsread('data1.xlsx','sheet1','j2:j12'); %读取excel表中arena仿真得到的数据,区别与matlab计算获得的Wq.
LLq=xlsread('data1.xlsx','sheet1','k2:k12');
rro=xlsread('data1.xlsx','sheet1','l2:l12');
nn_busy=xlsread('data1.xlsx','sheet1','m2:m12');
error_Wq=abs(100*(WWq-Wq))./WWq;
error_Lq=abs(100*(LLq-Lq))./LLq;
error_ro=abs(100*(rro-ro))./rro;
error_n_busy=abs(100*(nn_busy-n_busy))./nn_busy;
%将误差数据写入excel表中,以方便word直接复制。
xlswrite('data1.xlsx',error_Wq,'Sheet1','n2:n12');
xlswrite('data1.xlsx',error_Lq,'Sheet1','o2:o12');
xlswrite('data1.xlsx',error_ro,'Sheet1','p2:p12');
xlswrite('data1.xlsx',error_n_busy,'Sheet1','q2:q12');
MATLAB读写excel中指定sheet行列中的数据
最新推荐文章于 2024-07-11 18:11:12 发布