阻塞率matlab仿真

%本程序时间单位是小时 clear; k=0; %设定初始阻塞次数为0 s=3; %信道数为3 a=1/100; %服务时间 num=1000; %呼叫次数 c=zeros(1,s); %信道初始化，全部为空闲 y=zeros(1,s); %y为某一信道的剩余时间 g=rand(1,num); pr=zeros(600,1); for q=1:num %随机按指数分布产生服务时间 ser(q)=-a*log(1-g(q)); end for r=1:600 %r为信号到达率 k=0; x=rand(1,num); for n=1:num %产生信号时间间隔 arrive(n)=-(1/r)*log(1-x(n)); end for i=1:num for j=1:s %依次分析信道占用情况 if c(j)==0 c(j)=1; y(j)=ser(i); break end end if c(3)==1 k=k+1; %如果信道全部占用则阻塞次数加1 end for t=1:s %分析信号到达间隔期间的信道情况 y(t)=y(t)-arrive(i); if y(t)>0 c(t)=1; else c(t)=0; end end end pr(r)=k/num; %输出数据 pr(r) end save simulate.txt pr -ASCII