[论文解读] Exact Failure Frequency Calculations for Extended Systems

Annie Druault-Vicard, Christian Tanguy|Dec 28, 2006

Software Reliability and Analysis Research参考文献 27被引用 5

一句话总结

本文提出一种统一的单次遍历方法，通过矩阵乘积表达式和线性微分算子，计算大型可修复系统在稳态下的可用性与故障频率。关键贡献在于证明：对于 k-out-of-n:G、线性连续 k-out-of-n:F 及梯形网络等系统，故障率在渐近极限下与系统规模呈线性关系，从而可通过传递矩阵的特征值导数实现高效且精确的计算。

ABSTRACT

This paper shows how the steady-state availability and failure frequency can be calculated in a single pass for very large systems, when the availability is expressed as a product of matrices. We apply the general procedure to $k$-out-of-$n$:G and linear consecutive $k$-out-of-$n$:F systems, and to a simple ladder network in which each edge and node may fail. We also give the associated generating functions when the components have identical availabilities and failure rates. For large systems, the failure rate of the whole system is asymptotically proportional to its size. This paves the way to ready-to-use formulae for various architectures, as well as proof that the differential operator approach to failure frequency calculations is very useful and straightforward.

研究动机与目标

• 开发一种统一且高效的方法，用于精确计算大型可修复系统中的故障频率与稳态可用性。
• 展示微分算子方法在 k-out-of-n:G 和线性连续 k-out-of-n:F 等复杂系统架构中的有效性。
• 建立在系统组件完全相同时，故障率在渐近区域与系统规模呈线性关系的理论基础。
• 通过推导生成函数与基于特征值的近似方法，为常见网络拓扑提供即用型公式。

提出的方法

• 将系统可用性表示为传递矩阵的乘积，其中每个矩阵元素均为组件可用性的多重线性多项式。
• 对矩阵乘积应用线性微分算子 ∑λᵢpᵢ∂/∂pᵢ，以单次遍历方式计算故障频率。
• 故障频率通过将算子作用于乘积链中每个矩阵后所得矩阵的迹之和来计算。
• 对于相同组件，推导生成函数以分析渐近行为，故障率由传递矩阵主特征值的导数决定。
• 在 k-out-of-n:G、线性连续 k-out-of-n:F 及两终端梯形网络（含边与节点故障）上验证该方法。
• 渐近分析表明，故障率 ∝ n（当 n 较大时），比例常数由 ∂lnζ₊/∂lnp 和 ∂lnα₊/∂lnp 确定。

实验结果

研究问题

• RQ1 当可用性以矩阵乘积形式表达时，是否能以单次遍历方式精确且高效地计算大型系统的故障频率？
• RQ2 在大型 k-out-of-n:G 与线性连续 k-out-of-n:F 系统中，故障率的渐近行为如何？
• RQ3 对于具有相同组件的梯形网络，故障率如何随系统规模变化？
• RQ4 微分算子方法能否系统性地应用于同时包含边与节点故障的递归网络拓扑？
• RQ5 在组件可用性相同时，可用性与故障频率的精确生成函数是什么？

主要发现

• 通过在矩阵乘积上应用微分算子，可单次遍历精确计算故障频率，避免了昂贵的包含-排除展开。
• 对于大型系统，故障率在渐近下与系统规模 n 呈正比，比例常数由传递矩阵主特征值的导数决定。
• 在节点完全可靠的梯形网络中，故障率满足 ∂lnζ₊/∂lnp ≈ 1.138（在最优链路可靠性 p ≈ 0.2516 时），表明对 n 具有强烈的线性依赖性。
• 对于非常可靠的组件，故障率在首阶近似下为 (2n + 4)λq，与最小割集枚举结果一致。
• k-out-of-n:G 系统的生成函数为 Gₖ(z) = [pz / (1 - (1-p)z)]ᵏ / (1 - z)，支持闭式分析。
• 该方法可推广至基于不可用性的表达式，并在组件可靠性替代可用性时适用于时变系统。