@@ -19,27 +19,24 @@ kernelspec:
1919
2020本讲将运用基本工具来近似计算由多个关键部件组成的系统的年度故障率的概率分布。
2121
22- 我们将使用对数正态分布来近似关键组件部件的概率分布 。
22+ 我们将使用对数正态分布来近似关键部件的概率分布 。
2323
24- 为了近似描述整个系统故障率的n个对数正态概率分布之 ** 和** 的概率分布,我们将计算这n个对数正态概率分布的卷积 。
24+ 为了近似描述整个系统故障率的 $n$ 个对数正态概率分布之 ** 和** 的概率分布,我们将计算这 $n$ 个对数正态概率分布的卷积 。
2525
2626我们将使用以下概念和工具:
2727
2828* 对数正态分布
2929* 描述独立随机变量之和的概率分布的卷积定理
30-
3130* 用于近似多组件系统故障率的故障树分析
3231* 用于描述不确定概率的层次概率模型
33- * 傅里叶变换和逆傅里叶变换作为计算序列卷积的高效方法
32+ * 傅里叶变换和傅里叶逆变换作为计算序列卷积的高效方法
3433
3534关于傅里叶变换的更多信息,请参见这个 quantecon 讲座 [ 循环矩阵] ( https://python.quantecon.org/eig_circulant.html )
3635以及这些讲座 [ 协方差平稳过程] ( https://python-advanced.quantecon.org/arma.html ) 和 [ 谱估计] ( https://python-advanced.quantecon.org/estspec.html ) 。
3736
3837El-Shanawany, Ardron 和 Walker {cite}` Ardron_2018 ` 以及 Greenfield 和 Sargent {cite}` Greenfield_Sargent_1993 ` 使用了这里描述的一些方法来近似核设施安全系统的故障概率。
3938
40- 这些方法响应了 Apostolakis {cite}` apostolakis1990 ` 提出的关于构建量化程序的一些建议。
41-
42- 对安全系统可靠性的不确定性。
39+ 这些方法响应了 Apostolakis {cite}` apostolakis1990 ` 提出的关于构建用于量化安全系统可靠性的不确定性的程序的一些建议。
4340
4441我们先引入一些Python工具。
4542
@@ -77,7 +74,7 @@ np.set_printoptions(precision=3, suppress=True)
7774
7875 * $\mu$ 和 $\sigma^2$ 是 $x = \exp (y)$ 的均值和方差
7976 * 它们** 不是** $y$ 的均值和方差
80- * 相反,$y$ 的均值是 $e ^{\mu + \frac{1}{2} \sigma^2}$,$y$ 的方差是 $(e^{\sigma^2} - 1) e^{2 \mu + \sigma^2} $
77+ * 相反,$y$ 的均值是 $e ^{\mu + \frac{1}{2} \sigma^2}$,方差是 $(e^{\sigma^2} - 1) e^{2 \mu + \sigma^2} $
8178
8279对数正态随机变量 $y$ 是非负的。
8380
@@ -92,7 +89,6 @@ $$ f(y) = \frac{1}{y \sigma \sqrt{2 \pi}} \exp \left( \frac{- (\log y - \mu)^2
9289$$
9390\begin{aligned}
9491 \textrm{均值:} & \quad e ^{\mu + \frac{1}{2} \sigma^2} \cr
95-
9692\textrm{方差:} & \quad (e^{\sigma^2} - 1) e^{2 \mu + \sigma^2} \cr
9793 \textrm{中位数:} & \quad e^\mu \cr
9894 \textrm{众数:} & \quad e^{\mu - \sigma^2} \cr
10298$$
10399
104100
105- 回顾两个独立正态分布随机变量的以下_稳定性_性质 :
101+ 回顾两个独立正态分布随机变量的以下 * 稳定性 * 性质 :
106102
107- 如果$x_1$是均值为$\mu_1$、方差为$\sigma_1^2$的正态分布,且$x_2$独立于$x_1$并且是均值为$\mu_2$、方差为$\sigma_2^2$的正态分布,那么$x_1 + x_2$是均值为$\mu_1 + \mu_2$、方差为$\sigma_1^2 + \sigma_2^2$的正态分布。
103+ 如果 $x_1$ 是均值为 $\mu_1$、方差为 $\sigma_1^2$ 的正态分布,且 $x_2$ 独立于 $x_1$ 并且是均值为$\mu_2$、方差为 $\sigma_2^2$ 的正态分布,那么 $x_1 + x_2$ 是均值为 $\mu_1 + \mu_2$、方差为 $\sigma_1^2 + \sigma_2^2$ 的正态分布。
108104
109105
110- 独立的对数正态分布具有不同的_稳定性_性质 。
106+ 独立的对数正态分布具有不同的 * 稳定性 * 性质 。
111107
112108独立对数正态随机变量的** 乘积** 也是对数正态分布。
113109
114110
115- 特别地,如果$y_1$是参数为$(\mu_1, \sigma_1^2)$的对数正态分布,且
116-
117- $y_2$ 是对数正态分布,参数为 $(\mu_2, \sigma_2^2)$,那么乘积 $y_1 y_2$ 也是对数正态分布,其参数为 $(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2)$。
111+ 特别地,如果 $y_1$ 是参数为 $(\mu_1, \sigma_1^2)$ 的对数正态分布,且
112+ $y_2$ 是参数为 $(\mu_2, \sigma_2^2)$ 的对数正态分布,那么乘积 $y_1 y_2$ 也是对数正态分布,其参数为 $(\mu_1 + \mu_2, \sigma_1^2 + \sigma_2^2)$。
118113
119114``` {note}
120115虽然两个对数正态分布的乘积是对数正态分布,但两个对数正态分布的**和**却**不是**对数正态分布。
@@ -144,7 +139,7 @@ $$ h(z) = (f * g)(z) \equiv \int_{0}^\infty f (z) g(z - \tau) d \tau $$
144139
145140下面,我们将使用上述公式的离散化版本。
146141
147- 具体来说,我们将把 $f$ 和 $g$ 都替换为离散化的对应形式,并归一化使其和为1 ,这样它们就是概率分布。
142+ 具体来说,我们将把 $f$ 和 $g$ 都替换为离散化的对应形式,并归一化使其和为 $1$ ,这样它们就是概率分布。
148143
149144 * ** 离散化** 指的是等间隔采样的版本
150145
@@ -253,15 +248,15 @@ def pdf_seq(μ,σ,I,m):
253248```
254249
255250<!-- #region -->
256- 现在我们要为我们的离散化设置一个网格长度$I$和网格增量大小$m =1$。
251+ 现在我们要为我们的离散化设置一个网格长度 $I$ 和网格增量大小 $m =1$。
257252
258253``` {note}
259- 我们将$I$设置为2的幂 ,因为我们希望能够自由使用快速傅里叶变换来计算两个序列(离散分布)的卷积。
254+ 我们将 $I$ 设置为 $2$ 的幂 ,因为我们希望能够自由使用快速傅里叶变换来计算两个序列(离散分布)的卷积。
260255```
261256
262- 我们建议尝试2的不同幂值 $p$。
257+ 我们建议尝试 $2$ 的不同幂值 $p$。
263258
264- 例如,将其设置为15而不是12 ,可以改善离散化概率质量函数对所研究的原始连续概率密度函数的近似程度。
259+ 例如,将其设置为 $15$ 而不是 $12$ ,可以改善离散化概率质量函数对所研究的原始连续概率密度函数的近似程度。
265260
266261<!-- #endregion -->
267262
@@ -299,29 +294,29 @@ mean, meantheory
299294
300295## 概率质量函数的卷积
301296
302- 现在让我们使用卷积定理来计算上面参数化的两个对数正态随机变量之和的概率分布 。
297+ 现在我们使用卷积定理来计算上面参数化的两个对数正态随机变量之和的概率分布 。
303298
304299我们还将计算上面构造的三个对数正态分布之和的概率。
305300
306301在进行这些计算之前,我们需要解释我们选择的用于计算两个序列卷积的Python算法。
307302
308- 由于要进行卷积的序列很长,我们使用` scipy.signal.fftconvolve ` 函数而不是numpy.convolve函数 。
303+ 由于要进行卷积的序列很长,我们使用` scipy.signal.fftconvolve ` 函数而不是 ` numpy.convolve ` 函数 。
309304
310305这两个函数给出的结果实际上是等价的,但对于长序列来说,` scipy.signal.fftconvolve ` 要快得多。
311306
312307程序` scipy.signal.fftconvolve ` 使用快速傅里叶变换及其逆变换来计算卷积。
313308
314- 让我们定义傅里叶变换和逆傅里叶变换 。
309+ 让我们定义傅里叶变换和傅里叶逆变换 。
315310
316- 序列$\{ x_t\} _ {t=0}^{T-1}$的** 傅里叶变换** 是一个复数序列$\{ x(\omega_j)\} _ {j=0}^{T-1}$,由下式给出:
311+ 序列 $\{ x_t\} _ {t=0}^{T-1}$ 的** 傅里叶变换** 是一个复数序列 $\{ x(\omega_j)\} _ {j=0}^{T-1}$,由下式给出:
317312
318313$$
319314 x(\omega_j) = \sum_{t=0}^{T-1} x_t \exp(- i \omega_j t)
320315$$ (eq:ft1)
321316
322317其中 $\omega_j = \frac{2 \pi j}{T}$,$j=0, 1, \ldots, T-1$。
323318
324- 序列 $\{x(\omega_j)\}_{j=0}^{T-1}$ 的**逆傅里叶变换 **为
319+ 序列 $\{x(\omega_j)\}_{j=0}^{T-1}$ 的**傅里叶逆变换 **为
325320
326321$$
327322 x_t = T^{-1} \sum_ {j=0}^{T-1} x(\omega_j) \exp (i \omega_j t)
@@ -335,9 +330,9 @@ $$ (eq:ift1)
335330
336331- 计算序列 $\{f_k\}$ 和 $\{g_k\}$ 的傅里叶变换 $F(\omega)$、$G(\omega)$
337332- 形成乘积 $H (\omega) = F(\omega) G (\omega)$
338- - 卷积 $f * g$ 是 $H(\omega)$ 的逆傅里叶变换
333+ - 卷积 $f * g$ 是 $H(\omega)$ 的傅里叶逆变换
339334
340- **快速傅里叶变换**和相关的**逆快速傅里叶变换 **能够非常快速地执行这些计算。
335+ **快速傅里叶变换**和相关的**快速傅里叶逆变换 **能够非常快速地执行这些计算。
341336
342337这就是 `scipy.signal.fftconvolve` 使用的算法。
343338
@@ -427,11 +422,11 @@ mean, 3*meantheory
427422
428423在应用卷积定理之前,我们首先描述将组成事件与我们要量化其故障率的**顶端**事件连接起来的模型。
429424
430- 该模型是El -Shanawany、Ardron和Walker {cite}`Ardron_2018`所描述的广泛使用的 **故障树分析**的一个例子。
425+ 该模型是广泛使用的、El -Shanawany、Ardron和Walker {cite}`Ardron_2018`所描述的 **故障树分析**的一个例子。
431426
432427为了构建统计模型,我们反复使用所谓的**稀有事件近似**。
433428
434- 我们想要计算事件 $A \cup B$的概率。
429+ 假设我们要计算事件 $A \cup B$的概率。
435430
436431* 并集$A \cup B$是事件$A$或$B$发生的情况
437432
@@ -494,22 +489,22 @@ $$ (eq:probtop)
494489
495490 * $P(F)$ 的隐含概率分布的离散程度表征了分析师对系统失效概率的不确定性
496491
497- 这导致了有时被称为 **层次化**模型,其中分析师对概率$P(A_i)$本身也有概率估计。
492+ 这就是所谓的 **层次化**模型,其中分析师对概率$P(A_i)$本身也有概率估计。
498493
499494分析师通过以下假设来形式化他的不确定性:
500495
501- * 失效概率$P(A_i)$本身是一个对数正态随机变量,其参数为$(\mu_i, \sigma_i)$。
502- * 对于所有$i \neq j$的配对,失效率$P(A_i)$和 $P(A_j)$在统计上是相互独立的。
496+ * 失效概率 $P(A_i)$ 本身是一个对数正态随机变量,其参数为 $(\mu_i, \sigma_i)$。
497+ * 对于所有 $i \neq j$ 的配对,失效率 $P(A_i)$ 和 $P(A_j)$ 在统计上是相互独立的。
503498
504- 分析师通过阅读工程论文中的可靠性研究来校准失效事件$i = 1, \ldots, n$的参数$(\mu_i, \sigma_i)$,这些研究考察了与系统中使用的组件尽可能相似的组件的历史失效率。
499+ 分析师通过阅读工程论文中的可靠性研究来校准失效事件 $i = 1, \ldots, n$的参数$(\mu_i, \sigma_i)$,这些研究考察了与系统中使用的组件尽可能相似的组件的历史失效率。
505500
506501分析师假设,这些关于年度失效率或失效时间的观测分散性的信息,可以帮助他预测零件在其系统中的性能表现。
507502
508503分析师假设随机变量 $P(A_i)$ 在统计上是相互独立的。
509504
510505分析师想要近似系统失效概率 $P(F)$ 的概率质量函数和累积分布函数。
511506
512- * 我们说概率质量函数是因为我们之前描述的对每个随机变量的离散化方式 。
507+ * 我们说概率质量函数是因为我们对每个随机变量进行了离散化,正如前文描述的那样 。
513508
514509分析师通过重复应用卷积定理来计算**顶事件** $F$(即**系统失效**)的概率质量函数,以计算独立对数正态随机变量之和的概率分布,如方程 {eq}`eq:probtop` 所述。
515510
@@ -523,11 +518,11 @@ $$ (eq:probtop)
523518
524519这个例子是{cite}`Greenfield_Sargent_1993`第27页表10中描述的设计方案B-2(案例I)。
525520
526- 该表描述了十四个对数正态随机变量的参数$\mu_i, \sigma_i$,这些随机变量由**七对**独立同分布的随机变量组成。
521+ 该表描述了十四个对数正态随机变量的参数 $\mu_i, \sigma_i$,这些随机变量由**七对**独立同分布的随机变量组成。
527522
528523* 在每一对内,参数$\mu_i, \sigma_i$是相同的
529524
530- * 如{cite}`Greenfield_Sargent_1993`第27页表10所述,七个唯一概率$P(A_i)$的对数正态分布参数已被校准为以下Python代码中的值:
525+ * 如{cite}`Greenfield_Sargent_1993`第27页表10所述,七个唯一概率 $P(A_i)$ 的对数正态分布参数已被校准为以下Python代码中的值:
531526
532527```{code-cell} ipython3
533528mu1, sigma1 = 4.28, 1.1947
0 commit comments