注:在产生随机数过程中注意相关包的导入。
step1:创建3*4的矩阵
step2:基于包numpy下的泊松分布方法random.poisson(lam, size)生成所需的随机数。
x1=np.random.poisson(lam=5,size=4)
x2=np.random.poisson(lam=40,size=4)
x3=np.random.poisson(lam=10,size=4)
step3:对其进行可视化展示得到相应的效果图
a=plt.hist(x1,bins=pillar,normed=True,range=[0,pillar],color='g',alpha=0.5)
plt.plot(a[1][0:pillar],a[0],'r')
plt.grid()
plt.show()
step4:生成的随机数形成矩阵
C_zd=np.mat((x1,x2,x3))
得到的结果如下所示
C_zd= [[ 2 6 8 3]
[36 46 37 47]
[10 10 10 16]]
step5:对矩阵进行相关的计算
y=C_zd.sum(axis=1)#计算矩阵C_zd的边缘累积和
T=C_zd.sum(axis=0)
c=y.sum(axis=0)
M=y*T/c
得到的结果为
y= [[ 19]
[166]
[ 46]]
T= [[48 62 55 66]]
M= [[ 3.94805195 5.0995671 4.52380952 5.42857143]
[34.49350649 44.55411255 39.52380952 47.42857143]
[ 9.55844156 12.34632035 10.95238095 13.14285714]]