一种石英挠性加速度计参数稳定性建模方法技术

一种石英挠性加速度计参数稳定性建模方法，即一种基于随机点退化模型的石英挠性加速度计的参数稳定性建模方法，其步骤如下：一：建立含随机变化点的退化模型；二：观测数据的似然函数估计；三：从观测数据中分离缺失数；四：EM算法进行参数估计；通过以上步骤，本发明专利技术用产品主要性能参数来建立退化模型，对于高造价、高可靠性、长寿命产品来讲，节约试验成本和试验时间，有现实意义；本发明专利技术主要针对石英挠性加速度计的参数重复性差的问题，建立加速度计参数稳定性模型和随机变化点的估计，以确定使得加速度计参数进入稳定的剖面循环数，对于实际工程应用有着很大的现实意义。

【技术实现步骤摘要】
一种石英挠性加速度计参数稳定性建模方法
本专利技术提供一种石英挠性加速度计参数稳定性建模方法，它涉及一种基于随机点退化模型的石英挠性加速度计参数稳定性建模方法，它是基于含随机变化点的退化模型，针对石英挠性加速度计的参数重复性较差的问题，设计加速稳定性试验，对加速度计参数稳定性试验数据进行建模。这属于可靠性
中的加速退化试验数据处理。
技术介绍
石英挠性加速度计是一种机电产品，其主要功能是检测载体的直线运动信息。它广泛应用于导弹、卫星、火箭等飞行器上。在产品的整体制造和装配后，由于温度应力和残余应力的影响，加速度计的参数在时间和环境的影响下仍处于漂移状态。加速稳定性试验，是通过设计试验剖面，使得加速度计释放内部应力，以确保加速度计快速进入稳定状态。根据已设计的试验剖面，每一次循环试验之后离线测试一次，得到试验数据，当达到一定的试验循环数时，试验停止。根据收集到的试验退化数据建立退化模型，最终找到能使得加速度计进入稳定状态时的试验剖面循环数。在已有的退化模型研究工作中，出现了两大类方向：1.大多研究的退化模型是单一的退化轨迹。而实际工程中发现，产品的退化轨迹分为多阶段退化而并非单一。2.在研究中，产品的退化分为多阶段，但是将拐点视为固定不变的。而由于样本间的个体差异，使得变化点并非固定的，而是呈现出了随机性。因此，将退化轨迹多阶段和变化点随机性结合起来，利用加速度计稳定性试验的数据，对加速度计进行稳定性模型的参数估计，确定加速稳定试验剖面，是亟待解决的一个问题。
技术实现思路
本专利技术的目的在于提供一种石英挠性加速度计参数稳定性建模方法，它解决石英挠性加速度计的参数重复性较差(加速度计稳定性差)的问题，设计了加速稳定性试验，基于含随机变化点的多阶段退化模型，建立加速度计参数稳定性模型和随机变化点的估计，以确定使得加速度计参数进入稳定的剖面循环数。本专利技术是采用以下技术方案实现的，本专利技术是一种石英挠性加速度计参数稳定性建模方法，即一种基于随机点退化模型的石英挠性加速度计的参数稳定性建模方法，其步骤如下：步骤一：建立含随机变化点的退化模型试验通常采用了多个样本，因此退化轨迹由于样本间个体差异性，从一个阶段到另一个阶段的变化点，并非固定而考虑随机性；1.D(t)为t时刻的退化量，产品性能的退化速率在τ时刻发生非平滑的变化，τ为变化点，且假设产品的τ值服从某种随机分布，记τ的概率密度函数为其中θτ记为变化点τ待估参数的集合；认为D(t)为一个独立增量的随机过程，分阶段退化，其内的所有待估参数的集合记为θ；2.假设寿命试验中测试了I个样本，对于第i个样本，第j个观测时刻，分别在ti,j记录退化量Di,j，共获得ni个观测值；令τ1,τ2,...,τI分别为I个样本的变化点，记退化量为ΔDi,j＝Di,j+1-Di,j，时间增量为Δti,j＝ti,j+1-ti,j，从而对样本可获得一组退化增量且ΔDi,j的概率密度函数记为步骤二：观测数据的似然函数估计根据步骤一的退化模型构建观测数据的似然函数为：式中：为观测数据ΔDi,j的概率密度函数，为缺失数据变化点τi的概率密度函数；步骤三：从观测数据中分离缺失数据由步骤二中,极大似然估计的式子可以明显地看出，观测数据的似然函数中不仅包含观测数据Di,j，还包含缺失数据τi，因此必须将两者分离开来；1.对Lobs取对数可得完全数据对数似然函数为：其中，式中：表示变化点τi的对数似然函数，表示观测数据Di,j的对数似然函数，和表示概率密度函数取对数；2.上步中，由式子可以看出，对极大似然估计取对数之后，仍然没有完全分离观测数据Di,j和缺失数据τi；因此，需要进一步分离式子3.可将和写成如下形式：和其中例如：那么可得即可解决分离观测数据Di,j和缺失数据τi的问题；和分别表示对数似然函数Lobs的第一部分和第二部分；步骤四：EM算法进行参数估计EM算法包含E步(求期望)和M步(对期望极大化)；1.E步(求期望)：基于观测数据计算完全数据对数似然函数的期望；其中，E(v|ΔD)＝[f1,f2,f3]；上步中可知假定变化点τi的取值范围为：ti,j＜τi＜ti,j+1；因此，表达式可写为：式中：Q为完全数据对数似然函数的期望，Q1为缺失数据τi的对数似然函数期望，Q2为观测数据ΔDi,j的对数似然函数期望，和分别为式和中的系数矩阵；2.M步(对期望极大化)：1)选取一组参数估计初值θ0；θ0表示所有待估参数的初值集合，初始迭代次数q＝0；2)从E步期望公式获得Q(θ)；3)更新迭代步数q＝q+1，根据Q(θ)计算极大似然估计θq；4)如果||Q(θq+1|θq，D)-Q(θq|θq,D)||小于某一给定误差时，该程序终止；否则转到步骤2)；通过以上步骤，首先，可以看出该退化模型考虑到了样本个体间的差异性，使得退化轨迹中的变化点并非固定，而是随机分布的；其次，用产品主要性能参数来建立退化模型，对于高造价、高可靠性、长寿命产品来讲，节约试验成本和试验时间，有现实意义；最后，本专利技术主要针对石英挠性加速度计的参数重复性差的问题，建立加速度计参数稳定性模型和随机变化点的估计，以确定使得加速度计参数进入稳定的剖面循环数；基于以上三点明显的优势，本专利技术对于实际工程应用有着很大的现实意义。本专利技术具有以下优点：1.本专利技术利用加速度计主要性能参数，建立退化模型。在试验时间和试验成本的制约下，非常适用；2.本专利技术利用加速稳定试验方法，根据设计的试验剖面，最终确定加速度计达到稳定状态时的试验剖面循环数。解决了加速度计参数重复性较差的问题。对实际工程应用有着实际意义。附图说明图1本专利技术所述方法流程图。图2功率谱密度示意图。图3加速稳定性试验的剖面图；图4一号、三号、五号加速度计零偏压K0退化量的变化趋势图。图5一号、三号、五号加速度计K1退化量的变化趋势图。图6零偏压K0的变化点的置信区间图。图7K1的变化点的置信区间图。图中符号代号说明如下：图1中，K0表示石英挠性加速度计零偏压退化值,K1表示石英挠性加速度计退化值，每个虚线框表示一个步骤；图2中，该图为石英挠性加速度计稳定性试验一个循环周期下振动的功率谱密度图；图3中，gRMS表示振动幅值的单位；具体实施方式下面将结合附图和实施例对本专利技术作进一步的详细说明。以下实施例是以石英挠性加速度计为研究对象按照如图1所示的流程进行实施的，加速稳定性试验施加两种类型的应力：快速温度变化循环和随机振动。综合起来的应力细分为以下七个试验因素：高温、低温、温度加热/冷却变化率、高温保持时间、低温保持时间、振动幅度和振动时间。对于温度和振动幅度具体值的确定见表1。具体的试验剖面图见图3，图2是振动幅度的功率谱密度。本试验采用了8个加速度计样本进行试验。由图4,图5试验数据分析可知，该退化呈现两阶段，因此建立两阶段退化物理模型。最终确定石英挠性加速度计参数估计结果，得到加速度计关键性能参数稳定性模型，以确定加速度计达到稳定状态时的试验循环数。本专利技术一种石英挠性加速度计参数稳定性建模方法，即一种基于随机点退化模型的石英挠性加速度计的参数稳定性建模方法，见图1所示，其具体实施步骤如下：步骤一：建立含随机变化点的退化模型1.考虑到退化过程为一个随机过程，因此本专利技术选择Wie本文档来自技高网...

【技术保护点】
一种石英挠性加速度计参数稳定性建模方法，即一种基于随机点退化模型的石英挠性加速度计的参数稳定性建模方法，其特征在于：其步骤如下：步骤一：建立含随机变化点的退化模型试验通常采用了复数个样本，因此退化轨迹由于样本间个体差异性，从一个阶段到另一个阶段的变化点，并非固定而考虑随机性；1).D(t)为t时刻的退化量，产品性能的退化速率在τ时刻发生非平滑的变化，τ为变化点，且假设产品的τ值服从某种随机分布，记τ的概率密度函数为f

【技术特征摘要】
2016.10.19 CN 20161091066171.一种石英挠性加速度计参数稳定性建模方法，即一种基于随机点退化模型的石英挠性加速度计的参数稳定性建模方法，其特征在于：其步骤如下：步骤一：建立含随机变化点的退化模型试验通常采用了复数个样本，因此退化轨迹由于样本间个体差异性，从一个阶段到另一个阶段的变化点，并非固定而考虑随机性；1).D(t)为t时刻的退化量，产品性能的退化速率在τ时刻发生非平滑的变化，τ为变化点，且假设产品的τ值服从某种随机分布，记τ的概率密度函数为fτ(·；θτ)，其中θτ记为变化点τ待估参数的集合；认为D(t)为一个独立增量的随机过程，分阶段退化，其内的所有待估参数的集合记为θ；2).假设寿命试验中测试了I个样本，对于第i个样本，第j个观测时刻，分别在ti,j记录退化量Di,j，共获得ni个观测值；令τ1,τ2,...,τI分别为I个样本的变化点，记退化量为ΔDi,j＝Di,j+1-Di,j，时间增量为Δti,j＝ti,j+1-ti,j，从而对样本能获得一组退化增量且ΔDi,j的概率密度函数记为步骤二：观测数据的似然函数估计根据步骤一的退化模型构建观测数据的似然函数为：式中：为观测数据ΔDi,j的概率密度函数，为缺失数据变化点τi的概率密度函数；步骤三：从观测数据中分离缺失数据由步骤二中,极大似然估计的式子明显地看出，观测数据的似然函数中不仅包含观测数据Di,j，还包含缺失数据τi，因此必须将两者分离开来；1).对Lobs取对数得完全数据对数似然函数为：其中，式中：表示变化点τi的对数似然函数，表示观测数据Di,j的对数似然函数，和表示概率密度函数取...

【专利技术属性】
技术研发人员：许丹，何娇兰，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11