贾占强,梁保卫,杨建永
(中国电子科技集团公司第五十四研究所 石家庄 050081)
文献[1-5]从参数是否具有相关性方面入手,对多性能参数退化型产品的可靠性评估方法进行了研究,有效促进了相关理论和评估方法的研究。但在实际工程中,研究的对象通常具有结构复杂、信号特征复杂等特点,上述方法存在过程复杂、计算量大、建模困难、针对性和适用性较差等问题。文献[6-8]提出了基于Copula 函数的多性能参数可靠性评估方法,但当参数较多时,则会遇到计算量大和模型建立难等问题。因此,对于具有多性能参数的复杂电子产品来讲,有必要研究一种相对简便易用的可靠性评估方法,以促进可靠性理论的进一步工程化。
基于小波分析技术建立的小波熵理论被广泛应用于机械产品、模拟电路等机电产品的故障诊断[9-11],该方法的核心思想是:产品从正常到故障状态,其输出信号必然会发生变化,同时其信息熵也会发生相应的变化,通过对其输出信号信息熵的监测便可判断产品是否故障。本文借鉴这一思想,将小波熵理论引入到复杂电子产品的可靠性评估中,设计出了基于小波熵理论的多性能参数产品可靠性评估方法及步骤,并以具体实例验证了方法的有效性。
在信息论中,熵表示每个符号所提供的平均信息量和信源的平均不确定性,它能提供关于信号潜在的动态过程的有用信息。复杂信号经过小波变换后,映射到时间-尺度平面上,在任意时间间隔里,可以在多个尺度(不同分辨率)下观察信号的变化。假设每一个尺度为一个信号源,那么,每个尺度上的小波系数相当于一个信号源发出的信息。这样,根据小波变换系数,可以计算信号的小波熵。
根据小波变换的标架理论[12],当小波基函数是一组正交基函数时,变换具有能量守恒的性质:
从而可以定义小波能量为小波系数的平方和
由此可知总能量的表达式为
若将信号总能量归一化,则各子带信号的相对能量为
依据Shannon 熵理论,小波熵[13](wavelet energy entropy,记为HW)即定义为
上述方法适用于随机信号的处理,对于一个单一频率的周期信号来说,除了包含典型信号频率的小波尺度,所有的其它小波系数几乎都是零,而此时信号的熵值将接近于0 或者是一个很小的值。在这种情况下可采用如下方法求取信号的小波熵。由于在实际应用中,信号波形都是由示波器采集得到的离散信号,这里设采集到的信号为x(k)(k=1,2,···,N),即采样点数为N,则采样点的对数能量熵为
该信号的总能量为
基于小波能量熵的多性能参数退化数据的可靠性评估算法,主要用于解决具有多输出参数的复杂电子产品可靠性评估,其基本步骤为:
Step 1.假设有n个试验样品,收集和采集试验样品在时间t1,t2,···,tm的性能监测信号,按照小波熵或小波对数能量熵(如式(7))的计算方法,求出n个样品在各个时刻的HW(ti)或E(ti),i=1,2,···,m;
Step 2.失效阈值的确定,主要通过两种方法实现:①在实际使用中监测到故障信号;② 采用仿真的方法,当组成该产品的关键单元或器件达到其失效阈值时,仿真产生故障信号;
Step 3.按Step 2的方法对故障信号进行特征提取,得到失效阈值HWDf或EDf;
Step 4.在得到受试样品各监测时刻的HW(ti)或E(ti)(i=1,2,···,m)及失效阈值HWDf或EDf后,记为(ti,HW(ti))或(ti,E(ti))。由于随机变量通常服从正态分布和韦布尔分布,且韦布尔分布可适应多种分布形式(当m=1 时为指数分布,当m>3 时可近似为正态分布),因此,在此基础上,通常选择正态分布或威布尔分布建立基于退化量分布的可靠性评估模型:
①正态分布
当y≤Df时产品失效:
当y>Df时产品失效:
② 威布尔分布
当y≤Df时产品失效:
当y>Df时产品失效:
式(8)至式(11)中,y为各监测时刻的HW(ti)或E(ti),Df为失效阈值HWDf或EDf;
Step 5.对各时刻退化量分布函数的参数进行估计,估计退化量样本均值与样本方差或形状参数与尺度参数
Step 6.画出分布函数参数随时间变化的轨迹,选择恰当的曲线模型,并对曲线方程系数进行估计;
Step 7.给定失效阈值Df,根据求得的分布参数随时间变化的函数,利用产品可靠性与性能退化量分布参数的关系(即式(8)~式(11))对产品进行可靠性评估。
线性调频信号是雷达装备中常用的信号,线性调频信号的定义为
其中,A为信号幅度,β为扫频带宽,τ为脉宽。某型雷达接收机的线性调频信号是在展宽脉冲形成单元产生的,它由窄脉冲电路产生的δ脉冲冲击声表面滤波器,产生线性调频脉冲,经放大与选通电路对线性调频信号进行脉冲选通,再送到限幅放大单元,经幅度整形和放大后,送给上变频器;在上变频器中将中频线性调频信号系统与本振信号系统变成二者之和的微波信号,经功率放大,送给发射机,再经放大后向空间辐射出去。本文以线性调频信号形成单元为研究对象,对其可靠性进行分析。由于涉及该型雷达的关键指标,这里假设主要技术参数的标准值为:A=2 V,β=3 MHz,τ=20 μs。各参数的技术要求分别为β=3±0.5 MHz,τ=20±5 μs。试验前初始状态的标准信号波形如图1 中实线所示。
Step 1.安排10 个样本进行性能仿真试验,每100 h 测试一次数据,到3 000 h 停止试验;
Step 2.按照本文方法求取各样本各监测时刻的小波能量熵,结果如图2 所示;
Step 3.试验过程中,由于β和τ呈递增趋势变化,因此将β和τ技术要求的上限作为失效阈值,即β=3.5 MHz,τ=25 μs,对应的信号波形如图1 中点划线所示,此时信号波形的小波熵为0.624 9、0.868 0 和0.617 7,即选择小波熵最小值作为失效阈值,即HWDf=0.617 7;
Step 4.采用本文给出的可靠性评估算法来分析线性调频信号形成单元的可靠性;
①这里假设各监测时刻的小波熵服从正态分布,而后求取各监测时刻小波熵的均值和标准差,散点图如图3 所示。从图3 中可以看出,均值和标准差都是时间的线性函数,因此可以求出它们随时间变化的曲线方程:
② 而后将式(13)、式(14)及HWDf=0.621 77代入式(9)得到产品在给定时间t的可靠度函数(如式(15)),其对应的可靠度曲线如图4 中实线所示。
Step 5.图4 中实线为基于小波熵方法计算得到的多变量退化产品的可靠性评估结果,虚线为采用文献[5]中基于多元正态分布求解联合概率密度函数得到的计算结果,两者的评估结果比较接近,间接验证了本文方法的有效性。
本文从信号处理角度出发,基于小波熵理论研究了具有复杂特性参数的多性能参数产品的可靠性评估方法,给出了详细步骤与具体算法,并通过具体实例对有效性进行了验证。该方法克服了现有方法计算量大、求解相关函数难、建模复杂等实际问题。复杂电子产品的可靠性评估方法研究是一个具有实际工程意义但难度极大的课题,具体实现过程中还有诸多技术难题,如复杂系统的竞争失效问题、小子样和极小子样条件下的系统可靠性评估方法等,需要在今后工作中针对实际问题做进一步研究和探索。