MATLAB读写excel中指定sheet行列中的数据

%% 参考陆传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');