matlab关于商场人数的泊松随机过程仿真作业

X = random（‘exp‘1/2510001）; %构造服从指数分布的Xn序列，个数为1000个
S = [ ];
S（1） = 0; 
for i = 1:length（X）  %length（ ）测出向量长度  
    S（i+1） = S（i）+X（i）;
end %这个函数的作用是生成顾客进入商场的时间序列
N = [ ];M = [];W = [];
N（1）=0; M（1）=0;W（1）=0; %规定好初始值为0
b = rand（1length（S）-1）;
for i = 2:length（S）
    N（i）=i-1;  
    if b（i-1） <= 0.4
        M（i） = M（i-1）+1;
        W（i） = W（i-1）;
    else 
        M（i）=M（i-1）;
        W（i）=W（i-1）+1;