[发明专利]一种恒温器寿命DOE重要因子的识别方法

权利要求书

1.一种恒温器寿命DOE重要因子的识别方法，其特征在于，所述方法包括以下步骤：

S1、根据服从威布尔分布的恒温器寿命数据建立线性回归模型，其中线性回归模型包含尺度参数λ和形状参数ν、待估参数β和待估参数θ；

S2、利用扩散先验分布确定线性回归模型中的待估参数β和待估参数θ的先验分布：

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

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

其中，a₁，b₁，c₁，d₁，a₂，b₂，c₂，d₂为超参数；

S3、根据贝叶斯定理得到线性回归模型中的待估参数β和待估参数θ的后验分布，并从待估参数β和待估参数θ的后验分布中识别出重要因子；

所述步骤S1中利用恒温器寿命数据建立线性回归模型的具体步骤为：

X1、恒温器寿命数据服从威布尔分布，即y～dweibull(λ,v)

威布尔分布概率密度函数：f(y|λ,v)＝λv(y)^v-1exp[-λ(y)^v](1)威布尔分布累积分布函数：F(y|λ,v)＝1-exp[-λ(y)^v](2)

非截尾数据对似然函数的贡献为：

右截尾数据对似然函数的贡献为：

其中，m为试验因子水平组合的数量，共有n个样品，则y_ij表示为第i组第j个样品的寿命数据，且i＝1,2,...,m；j＝1,2,...,n；影响恒温器可靠性主要是由于腐蚀引起的隔膜穿孔，所述试验因子包括：隔膜镀层漂洗、电流密度、硫酸清洗、隔膜电解法净化时、铍铜合金晶粒尺寸、应力方向、湿度、热处理、焊接机、动力原件电解法净化程度和动力原件镀层清洗；

X2、对第i组试验因子x水平组合进行分析，威布尔分布下寿命数据的似然函数为：

X3、建立威布尔分布下有关尺度参数λ和形状参数ν与试验因子x之间的线性回归模型，即：

其中，k表示DOE中有k个试验因子x，x_ik为协变量。

2.如权利要求1所述的一种恒温器寿命DOE重要因子的识别方法，其特征在于，所述步骤S2中的扩散先验分布包括均匀分布和正态分布。

3.如权利要求1所述的一种恒温器寿命DOE重要因子的识别方法，其特征在于，所述步骤S3具体包括如下步骤：

Y1、利用贝叶斯定理，得到线性回归模型参数：

则边缘概率密度函数m(y)＝∫L(data|θ)f(θ)dθ为常数；

Y2、对待估参数β和待估参数θ进行收敛性判断；

Y3、从待估参数β和待估参数θ的后验分布中识别出重要因子，得到包含尺度参数λ和形状参数ν的线性回归具体模型。