import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import norm

data = [
86.08695652,
86.08695652,
86.08695652,
73.33333333,
69.75,
68.27586207,
66,
66,
61.875,
58.23529412,
56.57142857,
56.57142857,
50.4,
48.29268293,
48.29268293,
41.25,
36.66666667,
36,
34.875,
34.53061224,
34.2,
34.13793103,
33.81818182,
33.6,
33.6,
33.36585366,
30,
27.2195122,
27.2195122,
27.12328767,
26.52631579,
26.06896552,
25.89473684,
22.90909091,
21.29577465,
19.01886792,
18.66666667,
18.51428571,
18.32727273,
18,
17.52212389,
17.51351351,
17.37931034,
17.37931034,
17.37931034,
17.35714286,
17.2173913,
17.05263158,
16.77966102,
16.75862069,
16.61538462,
15.67741935,
14.97029703,
14.7761194,
14.45255474,
13.62162162,
13.62162162,
13.56164384,
13.31506849,
13.2885906,
13.2885906,
13.2,
13.11258278,
13.11258278,
13.11258278,
13.0619469,
13.02631579,
13.02631579,
12.94117647,
12.94117647,
12.85714286,
12.85714286,
12.61146497,
12.53164557,
12.53164557,
12.29813665,
12.22222222,
12.07317073,
11.72093023,
11.71597633,
11.71597633,
11.45454545,
10.7027027,
10.53191489,
10.3803681,
10.28571429,
10.18181818,
10.10958904,
10.05076142,
8,
7.675126904,
6.962264151,
6.904109589,
6.675496689,
10.28571429,
10.28571,
6
]

# 计算均值和标准差
mean = np.mean(data)
std_dev = np.std(data)

# 绘制直方图
plt.hist(data, bins=20, density=True, alpha=0.6, color='b', label='Histogram')

# 绘制正态曲线
xmin, xmax = plt.xlim()
x = np.linspace(xmin, xmax, 100)
p = norm.pdf(x, mean, std_dev)
plt.plot(x, p, 'k', linewidth=2, label='Normal Distribution')

plt.xlabel('Speed')
plt.ylabel('Probability Density')
plt.title('Velocity distribution and normal curve')
plt.legend()
plt.show()

print("均值:", mean)
print("标准差:", std_dev)