泊松过程模拟

两个程序：1.采用增量迭加法产生泊松过程。根据样本函数产生均值函数和方差函数，并对速率进行估计。2.采用点间间距迭加法产生泊松过程。

%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