[發明專利]一種利用特征數的威布爾型備件需求量的近似計算方法

有效

• 申請號: 201710343600.1
• 申請日: 2017-05-16
• 公開（公告）號: CN107220216B
• 公開（公告）日: 2020-07-28
• 發明（設計）人: 邵松世;翟亞利;李華;王俊龍
• 申請（專利權）人: 中國人民解放軍海軍工程大學
• 主分類號: G06F17/18
• 分類號: G06F17/18
• 代理公司: 武漢東喻專利代理事務所(普通合伙) 42224
• 代理人: 張英

1.一種利用特征數的威布爾型備件需求量的近似計算方法，所述威布爾型備件的壽命服從威布爾分布W(α,b)，α、b為威布爾分布參數，α為尺度參數，b為形狀參數；所述特征數為均值、方差、偏度和峰度，其特征在于，該計算方法包括如下步驟：

步驟一：利用所述威布爾分布參數α、b計算伽瑪分布的參數α_g、λ，

由威布爾分布的參數α、b，可得其均值為方差為其中Γ為伽瑪函數；當伽瑪分布的參數為α_g、λ時，其均值為方差為按照所求伽瑪分布的均值和方差，與威布爾分布的均值和方差相等的原則，計算出α_g、λ：

計算正態分布的參數μ、σ，

當正態分布的參數為μ、σ時，其均值為μ，方差為σ²；按照所求正態分布的均值和方差，與威布爾分布的均值和方差相等的原則，計算出μ、σ，

步驟二：計算偏度和峰度，

依據所述伽瑪分布的參數α_g、λ及所述正態分布的參數μ、σ，按下式計算三種分布的所述特征數中的偏度和峰度：

伽瑪分布，偏度為峰度為

正態分布，偏度為0，峰度為0；

威布爾分布，偏度為

峰度為

步驟三：按如下規則比較所述特征數中的偏度和峰度，

判斷所述步驟二中的伽瑪分布與正態分布中偏度和峰度與所述威布爾分布的偏度和峰度的絕對差值情況；

步驟三：計算備件需求量，

以所述偏度特征數的絕對差值為首要比較條件，相比正態分布，若所述步驟三所得伽瑪分布與所述威布爾分布的偏度絕對差值較小，則按照式(1)計算備件保障概率；否則，按照式(2)計算備件保障概率P_s：

式(1)、(2)中，T_w為保障任務時間，所述保障任務時間為所述備件完成任務的預期累積工作時間；

設置所述備件保障概率閾值，令j從0開始逐一遞增，使得所述保障概率P_s大于或等于所述概率閾值的j值即為計算出的備件需求量。