clear
lamda=2;Tmax=50;
delta_t=0.01;%時間精度
i=1;a=random（‘exponential‘lamda）;
T（1）=round（a*10）/10;
w（1）=T（1）;%初始化
%%%%%%%%%%%%%%泊松過程模擬%%%%%%%%%%%%%%%
while（w（i）T（i）=random（‘exponential‘lamda）;%構造服從指數分布的時間間隔序列Tn
???T（i）=round（T（i）*10）/10;
???w（i+1）=w（i）+T（i）;%計算等待時間
???i=i+1;
end
w=w‘;
x=zeros（w（1）/delta_t1）;
for?k=1:size（w1）-1
length=w（k+1）/delta_t-w（k）/delta_t;
x=[x;ones（length1）*k];%得到泊松分布X（t）序列
shuai=mean（x）;
end
%%%%%%%%%%%泊松過程檢驗%%%%%%%%%%%%%%%%%
alpha=0.05;
la