clc
clear
load P1.txt;
load P2.txt;
load P3.txt;
load P4.txt;
a=P1;
b=P2;
c=P3;
d=P4;
n=-5:1:30;
%绘制Poisson(0.2)的分布统计直方图和理论曲线
figure(1);
subplot(2,2,1);
hist(a,n);%绘制分布直方图
hold on;
u=10000*poisspdf(n,0.2);
plot(n,u,'g');%绘制理论曲线
title('Poisson(0.2)');
xlabel('x');
ylabel('PDF');
legend('simulation result','theory curve');
hold off
%绘制Poisson(1)的分布的统计直方图和理论曲线
subplot(2,2,2);
hist(b,n);%绘制分布直方图
hold on;
v=10000*poisspdf(n,1);
plot(n,v,'g');%绘制理论曲线
title('Poisson(1)');
xlabel('x');
ylabel('PDF');
hold off
%绘制Poisson(6)的分布的统计直方图和理论曲线
subplot(2,2,3);
hist(c,n);%绘制分布直方图
hold on;
w=10000*poisspdf(n,6);
plot(n,w,'g');%绘制理论曲线
title('Poisson(6)');
xlabel('x');
ylabel('PDF');
hold off
%绘制Poisson(10)的分布的统计直方图和理论曲线
subplot(2,2,4);
hist(d,n);%绘制分布直方图
hold on;
s=10000*poisspdf(n,10);
plot(n,s,'g');%绘制理论曲线
title('Poisson(10)');
xlabel('x');
ylabel('PDF');
hold off
评论1