python拟合威布尔分布

Python中可以使用SciPy库来进行威布尔分布的拟合威布尔分布是一种连续概率分布，常用于描述可靠性分析和寿命数据分析。下面是使用Python进行威布尔分布拟合的步骤：

1. 导入所需的库：

import numpy as np
from scipy.stats import weibull_min
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt

2. 准备数据：
首先，你需要准备一组威布尔分布的样本数据。假设你的数据存储在一个NumPy数组中，可以使用以下代码创建一个示例数据集：

data = np.array([10, 12, 15, 18, 20, 22, 25, 28, 30])

3. 定义威布尔分布函数：
在进行拟合之前，需要定义一个威布尔分布的概率密度函数。可以使用SciPy库中的`weibull_min`函数来定义威布尔分布函数：

def weibull(x, shape, loc, scale):
    return weibull_min.pdf(x, shape, loc, scale)

其中，`shape`、`loc`和`scale`是威布尔分布的参数。

4. 进行拟合：
使用`curve_fit`函数进行拟合，该函数会返回拟合后的参数值和协方差矩阵。以下是拟合的代码：

params, cov = curve_fit(weibull, data, bins, p0=[1, 0, 1])

其中，`data`是样本数据，`bins`是数据的分组数，`p0`是初始参数的猜测值。

5. 绘制拟合曲线：
可以使用以下代码绘制拟合后的威布尔分布曲线：

x = np.linspace(min(data), max(data), 100)
y = weibull(x, *params)

plt.plot(x, y, label='Fitted Weibull Distribution')
plt.hist(data, bins=bins, density=True, alpha=0.5, label='Histogram')
plt.legend()
plt.show()

以上就是使用Python进行威布尔分布拟合的基本步骤。