histogramdd 生成多维直方图数据并绘制3-D图

第一：函数说明：

histogramdd：可以生成一个多维度的直方图的数据。
Compute the multidimensional histogram of some data

函数说明

第二：参数解释：

2.1 sample： 用于生成统计直方图的输入数据，可以是多个维度的。

np.random.seed(0)                                 
sample1 = np.random.randn(10000) #1-D             
sample2 = np.random.randn(10000,2) # 2-D          
sample1.shape                                     
(10000,)
sample2.shape                                     
(10000, 2)

22 bins : 每个维度上的柱子个数和plt.hist参数一致。可以是个sequence，比如list，也可以是个int
如果是sequence：表示每个维度上每个bin的边界值。
如果是多个int，则代表每个维度上bin的个数。
如果是一个int，则代表所有维度上bin的个数，且相等。

2.3 range ：对每个维度上的数据进行显示的范围，默认为每个维度上的min到max和plt.hist参数一致
2.4 density：返回每个bin上的频率（density=True）或者频数（density= False， 默认） 和plt.hist参数一致，对于新版本，density取代了normed参数
2.5 normed：和density参数一样的作用，在新版本被density取代。
2.6 weights：每个数据的权重，和plt.hist参数一致

第三：返回值

返回值是一个2个元素组成的tuple，和plt.hist返回值类似，只是plt.hist返回值为3个元素组成的元组
第一个值：返回数据（可以是多维）在每个bin上的频率（density= True）或频数（density = False）。
第二个值：返回数据（可以是多维）上每个维度上bins的边界。

返回值

第四：其结果可以用于绘制图像

    N = 30
    density, edges = np.histogramdd(sample2, bins=[N, N])
    print('样本总数：', np.sum(density))
    density /= density.max()
    x = y = np.arange(N)
    print('x = ', x)
    print('y = ', y)
    t = np.meshgrid(x, y)
    print(t)
    fig = plt.figure(facecolor='w')
    # ax = fig.gca(projection='3d')
    # ax = Axes3D(fig)
    ax = fig.add_subplot(111, projection='3d')
    # s 参数是指点的大小, 值越大，点约大
    ax.scatter(t[0], t[1], density, c='r', s=50*density, marker='o', depthshade=True)
    # t[0], t[1], density分别代表X Y Z rstride（row）指定行的跨度  cstride(column)指定列的跨度 cmap = cm.Accent为颜色映射。
    #rstride = int1是指在行上的步长或跨度，即每次跨int1行，变现为Y轴上，因为Y轴上每个点，代表该行上Y的值。所以Y轴上就有bin_num/int1个值。
    #cstride = int2是指在列上的步长或跨度，即每次跨int2列，表现为X轴上，
因为X轴上每个店，代表带列上的X的值。所以X轴上就有bin_num/int2个值。
    ax.plot_surface(t[0], t[1], density, cmap=cm.Accent, rstride=1, cstride=1, alpha=0.9, lw=0.75, edgecolor='k')
    ax.set_xlabel('X')
    ax.set_ylabel('Y')
    ax.set_zlabel('Z')
    plt.title('二元高斯分布，样本个数：%d' % d.shape[0], fontsize=15)
    plt.tight_layout(0.1)
    plt.show()

ax.scatter s=100 *density

ax.scatter s=10 *density

s=50*density

s=1*density

s=100*density

rstride=10, cstride=10

rstride=3, cstride=5

rstride = int1是指在行上的步长或跨度，即每次跨int1行，变现为Y轴上，因为Y轴上每个点，代表该行上Y的值。所以Y轴上就有bin_num/int1个值。
cstride = int2是指在列上的步长或跨度，即每次跨int2列，表现为X轴上，
因为X轴上每个店，代表带列上的X的值。所以X轴上就有bin_num/int2个值。