隨機(jī)數(shù)組統(tǒng)計(jì)偏差對Weibull分布參數(shù)估計(jì)的影響

王桂金

(原鋼鐵研究總院，北京 100081)

1 均布于(0，1)隨機(jī)數(shù)的統(tǒng)計(jì)量

理想均勻分布于(0，1)的隨機(jī)數(shù)組具有如下統(tǒng)計(jì)量[1]：

均值(MEAN)μ0=0.5；

標(biāo)準(zhǔn)差(STDEV)σ0=0.288 68;

斜度 (SKEW)γ10=0.0;

過剩峭度(KURT)γ20=-1.2。

然而，在(0，1)準(zhǔn)均勻分布的1組隨機(jī)數(shù)組的統(tǒng)計(jì)量很難嚴(yán)格符合上述要求，存在如下偏差：

均值偏差δ(μ)=μ1-μ0;

標(biāo)準(zhǔn)差偏差δ(σ)=σ1-σ0;

斜度偏差δ(γ1)=γ11-γ10；

過剩峭度偏差δ(γ2)=γ21-γ20,

其中，下角標(biāo)1為準(zhǔn)隨機(jī)數(shù)組的統(tǒng)計(jì)量。根據(jù)文獻(xiàn)[2],隨機(jī)數(shù)α對應(yīng)的Weibull隨機(jī)壽命L為

L=λ0{ln(α)}1/κ0，

(1)

下文中模擬形狀參數(shù)κ0和尺寸參數(shù)λ0均設(shè)為1。

2 參數(shù)κ和λ的計(jì)算

按(1)式產(chǎn)生的隨機(jī)壽命組L，可用極大似然法(MLE)通過下式迭代計(jì)算獲得該數(shù)組的實(shí)際Weibull分布形狀參數(shù)κ和尺寸參數(shù)λ。

(2)

(3)

式中：N為樣本數(shù)。為保證結(jié)果的可靠性，計(jì)算進(jìn)行到(2)式左側(cè)絕對值小于10-5后停止。從文獻(xiàn)[3]中取得28組(0，1)準(zhǔn)均勻隨機(jī)數(shù)組，每組100個數(shù)據(jù)，各隨機(jī)數(shù)組的統(tǒng)計(jì)量及其隨機(jī)壽命的實(shí)際極大似然法Weibull參數(shù)κ和λ見表1。

表1 28組(0,1) 準(zhǔn)均勻隨機(jī)數(shù)組的統(tǒng)計(jì)量及其參數(shù)κ和λ

3 κ和λ的最小二乘法擬合

令隨機(jī)壽命數(shù)組的實(shí)際形狀參數(shù)κ和尺寸參數(shù)λ偏離期望值為

δ(κ)=κ-κ0;

δ(λ)=λ-λ0。

(4)

假定擬合的δ(κ1)和δ(λ1)分別為準(zhǔn)隨機(jī)數(shù)組4個統(tǒng)計(jì)量偏差的線性函數(shù),則

δ(m1)=a0+a1δ(μ)+a2δ(σ)+a3δ(γ1)+a4δ(γ2);m1=κ1,λ1，

(5)

于是應(yīng)用最小二乘法得

min{∑i[δ(m1)i-δ(m)i]2};i=1～28,m=κ,λ,

(6)

式中：δ(m)為隨機(jī)壽命實(shí)際Weibull參數(shù)偏離值，可分別得到δ(m1)的最佳系數(shù)a0,a1,a2,a3和a4。實(shí)際上，δ(μ),δ(σ),δ(γ1)和δ(γ2)分別與隨機(jī)參數(shù)分布的1～4次矩有關(guān)。

3.1 形狀參數(shù)κ的最小二乘法(LSD)擬合

經(jīng)過運(yùn)算，這28組隨機(jī)壽命擬合δ(κ1)的結(jié)果為

δ(κ1)=0.012 8-0.345δ(μ)-6.95δ(σ)+0.276δ(γ1)-0.287δ(γ2),

(7)

min{∑i[δ(κ1)i-δ(κ)i]2}=0.017 329,i=1～28。

最大和最小的絕對偏差為

max|δ(κ1)i-δ(κ)|=0.077 964，

min|δ(κ1)i-δ(κ)|=0.001 320。

最小二乘法擬合的形狀參數(shù)κ1和隨機(jī)壽命數(shù)組實(shí)際參數(shù)κ相差的絕對值平均只有0.004 7，而單個數(shù)據(jù)相差的絕對值不超過0.08。κ相對期望值1的最大偏離是+0.250 39和-0.125 95(表1)。因此用準(zhǔn)隨機(jī)數(shù)組的4個統(tǒng)計(jì)量偏差可以恰當(dāng)描述κ圍繞期望值的分布規(guī)律。把κ的偏差δ(κ)由小到大排列，再把擬合曲線的結(jié)果，即最小二乘法擬合形狀參數(shù)偏離值δ(κ1)畫在圖1中，可以看出擬合結(jié)果令人滿意。

圖1 28組δ(κ)與δ(κ1)值的比較

3.2 尺寸參數(shù)λ的最小二乘法擬合

經(jīng)過運(yùn)算，這28組隨機(jī)壽命擬合δ(λ1)的結(jié)果為

δ(λ1)=-0.002 4-3.485δ(μ)+

2.23δ(σ)-0.033δ(γ1)+0.165 5δ(γ2),

(8)

min{∑i[δ(λ1)i-δ(λ)i]2}=0.005 108,i=1～28。

最大和最小的絕對偏差為

max|δ(λ1)i-δ(λ)|=0.035 44,

min|δ(λ1)i-δ(λ)|=2.78×10-5。

最小二乘法得到的尺寸參數(shù)λ1和隨機(jī)數(shù)組實(shí)際參數(shù)λ平均僅相差0.002 6，而單個數(shù)據(jù)和隨機(jī)壽命的實(shí)際值相差不超過0.036。λ偏離期望值1的范圍是0.276 640～-0.273 886，因此，用準(zhǔn)隨機(jī)數(shù)組的4個統(tǒng)計(jì)量偏差同樣可以恰當(dāng)描述λ圍繞期望值的分布規(guī)律。把λ的偏差δ(λ)由小到大排列，再把擬合曲線的結(jié)果δ(λ1)畫在圖2中，可以看出擬合結(jié)果更令人滿意。

圖2 28組δ(λ)與δ(λ1)的比較

4 由試驗(yàn)數(shù)據(jù)產(chǎn)生的(0，1)隨機(jī)數(shù)組

(1)式的另一個重要應(yīng)用是可以從疲勞壽命試驗(yàn)數(shù)據(jù)產(chǎn)生與之相關(guān)的(0，1)準(zhǔn)隨機(jī)數(shù)組。為此對文獻(xiàn) [4] 中提到的 1208K+H208，6104, 7208和6307軸承[5-7]進(jìn)行了計(jì)算。作為例子，圖3給出H208軸承37個疲勞數(shù)據(jù)對應(yīng)的準(zhǔn)隨機(jī)數(shù)分布，其曲線接近于1條直線, 即接近于均勻分布(0，1)。為便于比較，表3同時列出這4種軸承隨機(jī)數(shù)分布統(tǒng)計(jì)量偏差，并代入(6)和(7) 式估計(jì)Weibull參數(shù)κ和λ的修正量。結(jié)果表明，修正后的額定壽命L10均降低，可以接受。

表3 4種軸承額定疲勞壽命L10的修正

圖3 1208K+H208軸承疲勞壽命對應(yīng)的(0， 1)準(zhǔn)隨機(jī)數(shù)序列

5 結(jié)論

(1)當(dāng)Weibull分布形狀參數(shù)κ和尺寸參數(shù)λ已知，用最大似然法重新計(jì)算28組隨機(jī)數(shù)組(每組100個數(shù)據(jù))的κ和λ，所得平均值很接近設(shè)定值，同時與文獻(xiàn)[8-9]的結(jié)果相吻合。

(2)各隨機(jī)壽命數(shù)組的Weibull參數(shù)κ和λ實(shí)際值相對設(shè)定值的偏差取決于隨機(jī)數(shù)組統(tǒng)計(jì)量(均值、標(biāo)準(zhǔn)差、斜度和過剩峭度)偏離期望值的程度，并且可以合理表達(dá)為這4個統(tǒng)計(jì)量偏差δ(μ),δ(σ),δ(γ1)和δ(γ2)的四元一次線性方程。這4個系數(shù)可由最小二乘法迭代計(jì)算求得。

(3) 從疲勞試驗(yàn)壽命數(shù)據(jù)組可以計(jì)算出其相應(yīng)(0，1)準(zhǔn)均勻分布的隨機(jī)數(shù)組，該數(shù)組的統(tǒng)計(jì)量可以反映出壽命數(shù)據(jù)的隨機(jī)特性，并可以對最大似然法計(jì)算出的Weibull分布參數(shù)κ和λ進(jìn)行修正。