Possion分布的车辆到达随机数

交叉口八个车道，每个进口两条车道的 车辆到达随机数，服从Possion分布

function a=GENTIME(N,rate)
y1=rand(1,N);
y2=rand(1,N);
y3=rand(1,N);
y4=rand(1,N);
y5=rand(1,N);
y6=rand(1,N);
y7=rand(1,N);
y8=rand(1,N);%生成5个范围为（0,1）的随机数，、上同，八条车道，因此有8行.
lamda=1/((3600/rate)-1)%到达率为6s(3.5s+1s)/辆，即4.5辆/s.
originy1=-log(1-y1)/lamda+1;
originy2=-log(1-y2)/lamda+1;
originy3=-log(1-y3)/lamda+1;
originy4=-log(1-y4)/lamda+1;
originy5=-log(1-y5)/lamda+1;
originy6=-log(1-y6)/lamda+1;
originy7=-log(1-y7)/lamda+1;
originy8=-log(1-y8)/lamda+1;%到达率为4.5辆/s的泊松分布
a=zeros(8,3);
sumtime1=0;
sumtime2=0;
sumtime3=0;
sumtime4=0;
sumtime5=0;
sumtime6=0;
sumtime7=0;
sumtime8=0;
for j=1:1:N
sumtime1=sumtime1+originy1(1,j)
sumtime2=sumtime2+originy2(1,j)
sumtime3=sumtime3+originy3(1,j)
sumtime4=sumtime4+originy4(1,j)
sumtime5=sumtime5+originy5(1,j)
sumtime6=sumtime6+originy6(1,j)
sumtime7=sumtime7+originy7(1,j)
sumtime8=sumtime8+originy8(1,j)
a(1,j)=sumtime1
a(2,j)=sumtime2
a(3,j)=sumtime3
a(4,j)=sumtime4
a(5,j)=sumtime5
a(6,j)=sumtime6
a(7,j)=sumtime7
a(8,j)=sumtime8
end
%上面为生成服从到达率为800辆/h的tij-,即车辆生成的时间，八条车道，独立随机生成