[發(fā)明專利]一種恒溫器壽命DOE重要因子的識別方法

有效

• 申請?zhí)枺?/td>: 202111111679.8
• 申請日: 2021-09-23
• 公開（公告）號: CN113705112B
• 公開（公告）日: 2023-06-27
• 發(fā)明（設(shè)計）人: 王國東;邵夢穎;孔美迪;薛麗;陳洪根
• 申請（專利權(quán)）人: 鄭州航空工業(yè)管理學(xué)院
• 主分類號: G06F30/27
• 分類號: G06F30/27;G06F17/18;G06N7/01;G06F111/08;G06F119/02
• 代理公司: 鄭州博派知識產(chǎn)權(quán)代理事務(wù)所(特殊普通合伙) 41137
• 代理人: 榮永輝

1.一種恒溫器壽命DOE重要因子的識別方法，其特征在于，所述方法包括以下步驟：

S1、根據(jù)服從威布爾分布的恒溫器壽命數(shù)據(jù)建立線性回歸模型，其中線性回歸模型包含尺度參數(shù)λ和形狀參數(shù)ν、待估參數(shù)β和待估參數(shù)θ；

S2、利用擴(kuò)散先驗分布確定線性回歸模型中的待估參數(shù)β和待估參數(shù)θ的先驗分布：

β_i～dunif(a₁,b₁)；θ_i～dunif(c₁,d₁)

β_i～dnorm(a₂,b₂)；θ_i～dnorm(c₂,d₂)

其中，a₁，b₁，c₁，d₁，a₂，b₂，c₂，d₂為超參數(shù)；

S3、根據(jù)貝葉斯定理得到線性回歸模型中的待估參數(shù)β和待估參數(shù)θ的后驗分布，并從待估參數(shù)β和待估參數(shù)θ的后驗分布中識別出重要因子；

所述步驟S1中利用恒溫器壽命數(shù)據(jù)建立線性回歸模型的具體步驟為：

X1、恒溫器壽命數(shù)據(jù)服從威布爾分布，即y～dweibull(λ,v)

威布爾分布概率密度函數(shù)：f(y|λ,v)＝λv(y)^v-1exp[-λ(y)^v](1)威布爾分布累積分布函數(shù)：F(y|λ,v)＝1-exp[-λ(y)^v](2)

非截尾數(shù)據(jù)對似然函數(shù)的貢獻(xiàn)為：

右截尾數(shù)據(jù)對似然函數(shù)的貢獻(xiàn)為：

其中，m為試驗因子水平組合的數(shù)量，共有n個樣品，則y_ij表示為第i組第j個樣品的壽命數(shù)據(jù)，且i＝1,2,...,m；j＝1,2,...,n；影響恒溫器可靠性主要是由于腐蝕引起的隔膜穿孔，所述試驗因子包括：隔膜鍍層漂洗、電流密度、硫酸清洗、隔膜電解法凈化時、鈹銅合金晶粒尺寸、應(yīng)力方向、濕度、熱處理、焊接機(jī)、動力原件電解法凈化程度和動力原件鍍層清洗；

X2、對第i組試驗因子x水平組合進(jìn)行分析，威布爾分布下壽命數(shù)據(jù)的似然函數(shù)為：

X3、建立威布爾分布下有關(guān)尺度參數(shù)λ和形狀參數(shù)ν與試驗因子x之間的線性回歸模型，即：

其中，k表示DOE中有k個試驗因子x，x_ik為協(xié)變量。

2.如權(quán)利要求1所述的一種恒溫器壽命DOE重要因子的識別方法，其特征在于，所述步驟S2中的擴(kuò)散先驗分布包括均勻分布和正態(tài)分布。

3.如權(quán)利要求1所述的一種恒溫器壽命DOE重要因子的識別方法，其特征在于，所述步驟S3具體包括如下步驟：

Y1、利用貝葉斯定理，得到線性回歸模型參數(shù)：

則邊緣概率密度函數(shù)m(y)＝∫L(data|θ)f(θ)dθ為常數(shù)；

Y2、對待估參數(shù)β和待估參數(shù)θ進(jìn)行收斂性判斷；

Y3、從待估參數(shù)β和待估參數(shù)θ的后驗分布中識別出重要因子，得到包含尺度參數(shù)λ和形狀參數(shù)ν的線性回歸具體模型。