泊松过程的模拟及检验 matlab程序

%泊松过程的模拟、检验及参数估计 
syms Un T W; 
n=10;%生成n*n个随机数 
r=1;%参数 
temp=0; 
tem=0; 
Un=rand（n1）;%共产生n*1个随机数 
for i=1:1:n 
    T（i）=-log（Un（i））/r; 
end 
T=subs（T）; 
for i=1:1:n 
    for j=1:1:i 
        temp=temp+T（j）;     
    end 
    W（i）=temp; 
    temp=0; 
end 
W=subs（W）; 
%检验泊松过程使用第四条 
for i=1:1:n 
    tem=tem+W（i）; 
end 
sigmaN=tem; 
Time=W（n）; 
alpha=0.05;%置信水平 
p=sigmaN/Time; 
p1=（1/2）*（n-1.96*（n/3）^（1/2））; 
p2=（1/2）*（n+1.96*（n/3）^（1/2））; 
c1=subs（p-p1） 
c2=su