有没有人有使用scipy.stats.distributions的示例代码？

import numpy as np
from scipy.stats import norm

# 生成正态分布的随机数
mean = 0
std_dev = 1
num_samples = 1000
random_samples = np.random.normal(mean, std_dev, num_samples)

# 计算正态分布的概率密度函数
x = np.linspace(norm.ppf(0.01, mean, std_dev), norm.ppf(0.99, mean, std_dev), 100)
y = norm.pdf(x, mean, std_dev)

# 绘制正态分布的概率密度函数图形
import matplotlib.pyplot as plt
plt.plot(x, y)
plt.hist(random_samples, bins=50, density=True, alpha=0.6, color='g')
plt.show()