論壇 › 可靠度工程與管理 › 壽命或失效時間分布的表示法

hlperng 發表於 2013-3-7 10:38:22

壽命或失效時間分布的表示法

本帖最後由 hlperng 於 2023-4-7 09:43 編輯 <br /><br />在可靠度領域裡，有很多機率分布可以作為描述失效發生時間 (time to failure) 或壽命 (life) 的隨機特性，常見的有指數 (exponential) 分布、韋伯 (Weibull) 分布、對數常態分布 (log-normal)、伽瑪 (Gamma) 分布等四種，這些機率分布可以應用的領域為大於等於零的隨機變數。&nbsp;<div><br></div><div>指數分布與伽瑪分布屬於計時或打點型 (pointing) 過程隨機變數的波桑 (Poisson) 家族機率分布，均值波桑隨機過程 (Homogeneous Poisson Process, HPP)。波桑家族機率分布，包括指數分布、歐郎分布 (Erlang)、伽瑪分布、波桑分布等，此一家族討論在選定時間內有興趣的對象事件（例如失效）會發生幾次，當事件發生強度會隨時間變化時，則延伸到韋伯分布或非均質波桑過程 (Non-homogeneous Poisson Process, NHPP）。&nbsp;
</div><div><br></div><div>每一種機率分布各有其用途與目的，早期美國國防部的電子設備可靠度諮詢小組 (Advisory Group on Reliability of Electronic Equipment, AGREE) 在處理電子產品可靠度問題時，以被動零件為主，假設壽命或失效時間為常態分布，後來轉向討論系統設備產品，多以指數分布為主說明失效時間，壽命或可靠度試驗數據分析時，涉及抽樣統計分析，使用伽瑪分布(或是實際計算時常用的卡方分布)，半導體產業初期構造較簡單假設積體電路的壽命為對數常態分布，發展至更複雜的微電路與積體電路時，漸漸使用韋伯分布，而機械業者一直都認為壽命或失效時間可以用韋伯分布表示。
&nbsp;</div><div><br></div><div>從工程觀點，壽命或失效發生時間都是大於等於零的變數，常態分布的領域由負無窮大到正無窮大，在應用上難免有罣礙，因此對數常態、韋伯及指數三種變數領域為大於等於零的機率分布最符合工程原則。標準常態分布與伽瑪分布，則是當推論機率分布參數（例如：失效率 λ）時所需使用到的抽樣分布。&nbsp;
</div><div><br></div><div>四個常用的壽命或失效時間分布中，從數學之間的關係，加上由物理科學（零件）、工程技術（組裝）、工程管理（系統設備產品）、驗證確認（管理運維）的系統工程及生命週期觀點，依次為對數常態分布、韋伯分布、指數分布、伽瑪分布。&nbsp;<span style="line-height: 1.5;">這四個機率分布之間關係的知識地圖說明如下表。&nbsp;</span></div><div><br></div><div><div><div><table cellspacing="0" cellpadding="0" width="40%" class="t_table" style="border: 1px solid rgb(227, 237, 245); table-layout: auto; width: 755.556px;"><tbody><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;"></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;">對數常態分布<br>lognormal distribution</td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;">韋伯分布<br>Weibull distribution</td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;">指數分布<br>exponential distribution</td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;">伽瑪分布<br>gamma distribution</td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">機率密度函數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="text-align: -webkit-center;">見公式 (1)</span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;">見公式 (2)</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="text-align: -webkit-center;"><img src="http://latex.codecogs.com/gif.latex?f(t)={\lambda}exp(-\lambda%20t)" title="f(t)={\lambda}exp(-\lambda t)"></span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;">見公式 (3)</td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">參數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="line-height: 1.5;">μ：尺度參數</span><br>σ：形狀參數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">η：尺度參數<br>β：形狀參數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">λ：尺度參數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">λ：尺度參數<br>κ：形狀參數</td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">失效率函數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="text-align: center;">(複雜公式)</span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="text-align: -webkit-center;"><img src="http://latex.codecogs.com/gif.latex?h(t)%20=%20\frac%20{\beta}%20{\eta}%20(\frac%20{t}%20{\eta})^{\beta%20-%201}" title="h(t) = \frac {\beta} {\eta} (\frac {t} {\eta})^{\beta - 1}"></span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?h(t)%20=%20\lambda" title="h(t) = \lambda  "></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?h(t)%20=%20\frac%20{\lambda%20(\lambda%20t)^{\kappa%20-%201}%20exp(-%20\lambda%20t)}%20{\Gamma(\lambda%20t,%20\kappa)}" title="h(t) = \frac {\lambda (\lambda t)^{\kappa - 1} exp(- \lambda t)} {\Gamma(\lambda t, \kappa)}   "></td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">平均數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="text-align: center;"><img src="http://latex.codecogs.com/gif.latex?\theta%20=%20exp(\mu%20+%20{\frac%20{\sigma^2}%20{2}})" title="\theta = exp(\mu + {\frac {\sigma^2} {2}})"></span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><span style="text-align: -webkit-center;"><img src="http://latex.codecogs.com/gif.latex?\theta%20=%20\eta%20\Gamma%20(\frac%20{1}%20{\beta}%20+%201%20)" title="\theta = \eta \Gamma (\frac {1} {\beta} + 1 )"></span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?\theta%20=%20{{1}%20\over%20{\lambda}}" title="\theta = {{1} \over {\lambda}}"></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?\theta%20=%20\frac%20{\kappa}%20{\lambda}" title="\theta = \frac {\kappa} {\lambda}"></td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">特徵數</td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?t_{0.63}" title="t_{0.63}"></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?t_{0.63}" title="t_{0.63} "></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?t_{0.63}" title="t_{0.63}"></td><td style="border-color: rgb(227, 237, 245); overflow: hidden; text-align: center;"><img src="http://latex.codecogs.com/gif.latex?t_{0.50}" title="t_{0.50}"></td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">機率密度函數圖</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><ignore_js_op><img id="aimg_386" src="http://tw-redi.com/data/attachment/forum/201412/25/160259tvkrk452v23tc6h4.png" zoomfile="data/attachment/forum/201412/25/160259tvkrk452v23tc6h4.png" file="data/attachment/forum/201412/25/160259tvkrk452v23tc6h4.png" class="zoom" width="120" inpost="1" alt="pns_lognor-pdf-02" title="pns_lognor-pdf-02" style="cursor: pointer;"></ignore_js_op></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><ignore_js_op><img id="aimg_387" src="http://tw-redi.com/data/attachment/forum/201412/25/160327o5n4vetioqz62v7q.png" zoomfile="data/attachment/forum/201412/25/160327o5n4vetioqz62v7q.png" file="data/attachment/forum/201412/25/160327o5n4vetioqz62v7q.png" class="zoom" width="120" inpost="1" alt="pns_wei-pdf-02" title="pns_wei-pdf-02" initialized="true" style="cursor: pointer;"></ignore_js_op></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><div style="text-align: center;"></div><br></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;"><br></td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">特性<br></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">形狀參數 σ 愈大、分布愈寬。<br><img src="http://latex.codecogs.com/gif.latex?\rho%20=%20\frac%20{\sigma_T}%20{\mu_T}" title="\rho = \frac {\sigma_T} {\mu_T}"><br>ρ = 變異係數<br>σ ≈ 0.1，對數常態近似常態。<br>隨機變數有乘法性（對數常態分布隨機變數的乘或除結果所得到的隨機變數也是對數常態分布）。</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">形狀參數 β 愈大、分佈愈窄。<br>令 t'=t^{\beta}，則 t' 為指數分布。<br>β = 1 時為指數分布。<br>β ≈ 4.0，韋伯近似常態。<br>一般建議，當形狀參數遠大於 4.0 時 (β &gt;&gt; 4.0)，要特別注意資料的物理合理性。</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">尺度參數（表示<span style="line-height: 1.5;">失效率時）有加法性。</span></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">形狀參數 κ 愈大、分布愈寬。<br>κ ≈ 30，伽瑪近似常態。<br>隨機變數有加法性（伽瑪分布隨機變數的相加結果所得到的隨機變數也是伽瑪分布）。</td></tr><tr><td style="border-color: rgb(227, 237, 245); overflow: hidden;">應用<br></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">物理科學<br>零件<br>單一失效模式<br>半導體產品</td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">工程技術<br>組裝件<br>實體失效模式（多重，最弱鏈環）&nbsp;<br>機械產品<br>電子元件<br><br></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">營運管理<br>系統<br>功能失效模式<br>電子產品(設備)<br><br></td><td style="border-color: rgb(227, 237, 245); overflow: hidden;">驗證保證<br>系統設備<br>壽命試驗<br>現場操作<br>抽樣推論（推定與檢定）</td></tr></tbody></table></div><div><br></div><div><div><span style="line-height: 1.5;">公式 1：</span>對數常態分布機率密度函數有下列兩種寫法</div><div><span style="line-height: 1.5;"> f(t) = \frac {1}{t \sigma \sqrt {2 \pi}} exp \left( -\frac{1}{2} \left(\frac {t - \mu}{\sigma} \right)^2 \right) </span></div><div><span style="line-height: 1.5;">或</span></div><div><span style="line-height: 1.5;">f(t) = \frac {1}{t \sigma \sqrt{2\pi}} exp\left(-{{1} \over {2}} \left(\frac {ln(\frac {t} {t_{0.50}}) } {\sigma}\right)^2\right) &nbsp;</span></div><div><br></div><div><span style="line-height: 1.5;">公式 2：</span>韋伯分布機率密度函數</div><div><span style="line-height: 1.5;">f(t) = \frac {\beta} {\eta} \left(\frac {t} {\eta} \right)^{\beta - 1} exp \left( - \left(\frac {t} {\eta} \right)^ {\beta} \right) &nbsp;</span></div><div><span style="line-height: 1.5;">或</span></div><div><span style="line-height: 1.5;">f(t) =\frac {\beta} {t_{0.63}} \left(\frac {t} {t_{0.63}} \right)^{\beta -1} exp \left(- \left(\frac {t} {t_{0.63}} \right)^{\beta} \right)&nbsp;</span></div><div><br></div><div><span style="line-height: 1.5;">公式 3：</span>伽瑪分布機率密度函數</div><div><span style="line-height: 1.5;">f(t) = \frac {\lambda} {\Gamma (\kappa)} (\lambda t)^{\kappa - 1} exp( - \lambda t) &nbsp;</span></div></div><div><br></div><div><span style="line-height: 1.5;">韋伯分布係因為研究機械產品失效發生時間或壽命而發展得到，從機率數學觀點，本身又是有界 (bounded) 極小值分布（剛柏 III 型極小值，Gumbel Type III Minimum）的特性，可用於解釋系統中最小鏈環、最弱環節 (weakest linkage) 的問題，加上數學式子簡潔，機率分布的尺度參數 η 代表特徵壽命、形狀參數 β 與失效機制的工程關聯意義合理，加上形狀參數 β &lt; 1、β = 1、β &gt; 1 可以說明物品生命週期早夭、有用、及磨耗三個階段的現象與特質，其應用空間極為寬廣，是描述可靠度進階的機率分布。&nbsp;</span></div></div></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">指數分布是韋伯分布的特例，當韋伯分佈的形狀參數 β = 1 時，韋伯分布變成指數分布。美軍可靠度試驗標準 MIL-STD-781，就是以指數分布作為軍用系統設備的壽命或失效時間分佈，主要是因為指數分佈簡單、易懂，且頗符合電子產品的失效時間或壽命的隨機發生特性（波桑過程）。
&nbsp;</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">伽瑪分布是 κ 個指數分布隨機變數加總之後的新隨機變數的機率分布，伽瑪分布的隨機變數有加法性，當 κ=1，伽瑪分布變成指數分佈，指數分布為伽瑪分布的特例。伽瑪分布適合於描述壽命試驗累積失效時間的機率分布。</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">卡方分布為</span><span style="line-height: 1.5;">伽瑪分布經過適當的數學轉換處理後的機率分布，當伽瑪分布的形狀參數 κ 為整數時，可以利用以兩倍形狀參數 κ 為形狀參數（自由度）(ν = 2 κ) 的卡方分布。</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">一般在處理壽命或失效時間的統計分析問題，例如：判斷試驗數據是否有太大或太小的離群數、計算尺度參數失效率 λ 的區間推定值、或檢定樣本失效率與假設失效率真值是否相同，大多應用卡方分布處理，就是因為指數分布、伽瑪分布與卡方分布之間有密切關係的緣故。&nbsp;</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">伽瑪分布可以用下述兩種數學表達方式：</span></div><div><div>f(t) = \frac {\lambda} {\Gamma (\kappa)} (\lambda t)^{\kappa - 1} exp( - \lambda t) </div><div>或者</div><div>f(t) = \frac {1} {\eta {\Gamma (\kappa)}} \left(\frac {t} {\eta} \right)^{\kappa-1} exp \left(- \frac {t} {\eta} \right)&nbsp;</div><div>其中：&nbsp;</div><div>\eta = \frac {1} {\lambda} 。&nbsp;</div><div><br></div><div>重點是如何解釋與應用這兩種數學式子及所要描述或解釋的問題。當 t 為失效時間或壽命時，單位為小時 (hr)，第一種寫法的尺度參數 λ 為失效強度 (failure intensity)，單位為小時-1 (fr/hr)；第二種寫法的尺度參數 η 為特徵壽命，單位為小時 (hr)。失效強度表示造成物品失效的量化根源，特徵壽命則為表示 63.2 % 物品失效的時間點。</div><div><br></div><div>失效率是從底層、微觀、根源、工程技術的角度，描述、觀察、與處理物品的可靠度，屬於物品的內在或本質 (intrinsic) 特性；特徵壽命是從表面、宏觀、結果、驗證保證的角度，描述、觀察、及處理物品的可靠度。

失效率與特徵壽命兩者之間有互為因果的關係，但是在應用上有其個別的目的與特色。</div><div><br></div><div>失效率可以根據樣本的壽命試驗結果，配合假設的機率理論推論出來的，因此又稱之為表觀失效率或視在失效率 (apparent failure rate)。視在失效率也是隨機變數，因此 2006 年 Telcordia 發佈 SR-332 《電子設備可靠度預估程序》時，就假設失效率為伽瑪分佈，每一種電子零件的失效率都提供平均值與標準差兩組數據。表面觀察資料雖然接近現實，但又受限於觀察樣本，追根究底就是想要瞭解事物的根源，發現原貌。只有瞭解或掌握根源，才能解決問題。&nbsp;
</div><div><br></div><div>指數分布既是伽瑪分布的特例（伽瑪分布為 κ 個指數分布之和的分布，κ 為伽瑪分布的形狀參數），也是韋伯分布的特例（當韋伯分布的形狀參數 β =1 時），那麼韋伯分布與伽瑪分布有沒有關係，若有的話又是什麼關係？關鍵在於韋伯分布與伽瑪分布應用兩個數學的特性，也就是說 κ 與 β 這兩個形狀參數的物理意義，它們在可靠度領域時對可靠度的解釋有何幫助。&nbsp;</div><div><br></div><div>美國主導半導體裝置研發與製造的電子產業協會 (EIA)（後來獨立成為JEDEC），早期認定半導體的壽命分布為對數常態分布，最近則有建議採韋伯分布的趨勢（參閱 Section 5.19 Reliability Distributions,<span style="line-height: 1.5;">&nbsp;JEDEC JEP-122G, 2011）。電容器的壽命也假設符合韋伯分布，可靠度等級為 B 級、C 級、D 級等，與傳統假設為指數分佈的 L 級、M 級、S 級有所區別。
&nbsp;</span></div><div><span style="line-height: 1.5;"><br></span></div><div><span style="line-height: 1.5;">數學是一種工具，可靠度工程就是工程概念在探討物品可靠度的應用，只要符合定義與邏輯、簡單易懂，最能協助描述事物，就是最佳解答。&nbsp;</span></div><div><br></div><div>&nbsp;(備註：本文因2013年3月27日系統故障及操作不當產生亂碼，2013年6月22日已重新編輯。)</div></div><div><br></div><div><br></div><div><br></div>

a69071183 發表於 2014-9-18 13:08:53

上面所述敘的公式，可以藉由那些軟體取得

hlperng 發表於 2014-9-18 23:06:22

a69071183 發表於 2014-9-18 13:08 static/image/common/back.gif
上面所述敘的公式，可以藉由那些軟體取得

Microsoft Excel !

liaojenyi 發表於 2014-9-19 19:43:54

http://office.microsoft.com/en-us/excel-help/gammadist-HP005209101.aspx

查看完整版本: 壽命或失效時間分布的表示法