如何为双变量数据生成3D联合概率分布(例如，绘图)？

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# 准备双变量数据
x = np.random.normal(0, 1, 1000)
y = np.random.normal(0, 1, 1000)

# 统计分析，计算联合概率分布
hist, xedges, yedges = np.histogram2d(x, y, bins=10, range=[[-3, 3], [-3, 3]])
xpos, ypos = np.meshgrid(xedges[:-1] + 0.25, yedges[:-1] + 0.25)
xpos = xpos.flatten()
ypos = ypos.flatten()
zpos = np.zeros_like(xpos)
dx = dy = 0.5
dz = hist.flatten()

# 绘制3D图形
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.bar3d(xpos, ypos, zpos, dx, dy, dz, color='b', alpha=0.8)

# 设置图形参数
ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Probability')

# 显示图形
plt.show()