随机模拟泊松流.m

来自「many matlab script ,i colect many things」· M 代码 · 共 25 行

for D=5000:5800
x(1)=round(D*0.4);
y(1)=round(D*0.6);
xx=1:60;
p1=poissrnd(30,x(1),1);
p2=poissrnd(15,y(1),1);
px=hist(p1,xx)/x(1);
py=hist(p2,xx)/y(1);
sx=x(1);
sy=y(1);
for i=2:30
    s=0;
    for j=1:i-1
        s=s+round(x(j)*px(i-j))+round(y(j)*py(i-j));
    end
    x(i)=s*0.4;
    y(i)=s*0.6;
    sx=sx+x(i);
    sy=sy+y(i);
end
if (sx>=4000&sy>=6000)
    break;
end
end
D