X=-0:0.01:30;
Y1=poisspdf(X,20);
subplot(2,2,1);
plot(X,Y1);
title('概率密度');
Y2=poisscdf(X,20);
subplot(2,2,2);
plot(X,Y2);
title('概率分布');
R=xcorr(Y1,'coeff');
subplot(2,2,3);
plot(R);
title('相關函數');
F1=fft(R);         %由相關函數通過快速傅里葉變換得到功率譜密度
Fx=(0:length(F1)-1)*1000/length(F1);
subplot(2,2,4)
plot(Fx,10*log10(abs(F1)));
title('功率譜密度');