%lamda為λ。Tmax為時間長度，hao是時間間隔，Num是產生的樣本函數的個數；
lamda=2;Tmax=300;hao=0.1;?Num=100;
sum=zeros（1Tmax）;
N=zeros（NumTmax）;
for?j=1:Num
????i=1;
????N（j1）=?0;
????while（i????????N（ji+1）=N（ji）+poissrnd（lamda*hao）;
????????i=i+1;
????end
????%subplot（22j）;
????%stairs（N（j:））;
????sum=?sum+N（j:）;
end
%title（‘泊松隨機過程的樣本函數‘）;

sum=sum/Num;%產生均值函數；
figure（）;
plot（sum）;%顯示均值函數。
title（‘均值函數‘）;

var=zeros（1Tmax）;
for?j=1:Tmax
????for?i=1:Num
????????var（j）=（（N（ij）-sum（j））^2）/（Num-1）;
????end
end
figure（）;
plot（var）;%顯示方差函數
title（‘方差函數‘）;
n=zeros（NumTmax）;
lam=0;
for?i=1:Tmax
????if?i==1
????????n（:i）=N（:i）;
????else
????????n（:i）=N（:i）-N（:i-1）;
????end
????
????Lambdahat=poissfit（n（:i））;
????lam=lam+Lambdahat;
end
lam=lam/Tmax;
lam=lam*hao;%得到對速率的估計

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-12-22?11:43??possion\
?????文件?????????883??2012-12-22?11:11??possion\possion.m
?????文件?????????239??2012-12-21?21:05??possion\possion2.m