在Python中拟合Maxwell-Boltzman分布

Question

是否有可能在matplotlib或python中的类似模块中拟合Maxwell-Boltzmann数据？

Answer 1

scipy.stats支持maxwell distribution。

---

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

maxwell = stats.maxwell 
data = maxwell.rvs(loc=0, scale=5, size=10000) 

params = maxwell.fit(data, floc=0) 
print(params) 
# (0, 4.9808603062591041) 

plt.hist(data, bins=20, normed=True) 
x = np.linspace(0, 25, 100) 
plt.plot(x, maxwell.pdf(x, *params), lw=3) 
plt.show() 

的第一个参数是位置或从零移位远离。 第二个参数是缩放参数，在the wikipedia page上由a表示。

为了生成随机变元（随机数据）与该分布，使用其rvs方法：

newdata = maxwell.rvs(*params, size=100)