举个例子：
lambda = 2;
r = poissrnd(lambda,10000,1);
mean(r) % 均值
var(r) % 方差
y = poisspdf(r,lambda); % 概率密度
...功率谱应该可以用psd函数,自相关用xcorr,画图用plot...
依次往下算就是了.