本发明涉及一种模拟电路健康预测方法。
背景技术：
随着科学技术的快速发展，电子设备广泛应用于生活中的各行各业，模拟电路作为电子设备中的重要组成部分，人们对其可靠性的要求越来越高，尤其是军工、航天等领域，即使是微小的误差都有可能造成重大的事故，因此，提高模拟电路的可靠性具有十分重要的意义。近几年来故障预测一直是国内外研究的热点问题，例如著名的美国航空航天局(NASA)，作为美国联邦政府的重要政府机构负责太空计划，一直致力于研究飞行器的故障预测技术。同时，中国也将故障预测技术作为研究重点领域之一。虽然在电子设备中模拟电路的应用没有数字电路广泛，但是有统计资料显示电子设备总故障中的80％来自于模拟电路，可见对模拟电路进行健康预测具有十分重要的意义和必要性。
目前,国内外对模拟电路的故障与预测的研究还不成熟，处于基础理论研究阶段。模拟电路故障健康预测主要存在的问题有以下四个方面。
(1)电路故障的标准不同。电路发生故障的标准是进行故障预测的基础，而国内外发表的很多学术论文中，对故障的定义则不同。
(2)电路各元件容限不同。一般来说，模拟电路中包含大量的电子元件，越复杂的电路，元件越多，也越容易发生故障。对模拟电路进行故障预测，分析各元件失效机理时，必须考虑各元件的容限，而各元件的容限不同则造成了预测电路故障的困难。
(3)复杂的失效机理。电路中通常包含电阻、电容、电感、运放等元件，而各元件的退化、失效机理不同，造成了故障预测的又一困难。
(4)人为操作和环境因素的影响。人为操作对电路的故障预测存在不确定因素，不恰当的操作可能会加快电路发生故障，因此在分析上存在难点。同时，包括温度、湿度、压强等在内的环境因素也是影响电路故障预测的重要因素。
技术实现要素：
本发明要解决的技术问题是，克服现有技术存在的上述缺陷，提供一种模拟电路故障定期检修到智能预测检修过渡的模拟电路健康预测方法。
本发明之模拟电路健康预测方法的实质在于利用经验小波技术把电路样本容差信号和待测时刻的偏差信号作经验小波分解，分别得到其基频信号，然后根据统计学知识求取各自基频信号的能量熵值，以电路无故障模式下的容差电压信号的基频能量熵与电路待测时刻电压增量信号的基频能量熵的比值作为电路健康指数(HI),并以此来判断电路的健康状况。
本发明解决其技术问题所采用的技术方案是：
一种模拟电路健康预测方法，包括如下步骤：
(1)采集电路容差、无故障模式下的测试样本集：利用数据采集器获取待测模拟电路容差、无故障模式F0下的M组电压样本向量其中为采集到的无故障模式模式F0的第i组样本向量；
(2)计算电路容差、无故障模式F0下的M组电压样本向量的统计平均值
(3)计算电路容差、无故障模式F0下模拟电路测试样本容差电压增量信号集△V0：i＝1,2,…M，令则i＝1,2,…M，其中为第i个样本的容差电压增量信号；
(4)计算电路容差、无故障模式F0下的容差电压增量信号i＝1,2,…M的基频能量熵
(5)计算无故障容差电压增量信号的最大能量熵：其中max(·)代表取最大值；
(6)利用数据采集器采集电路待测时刻电压信号v；
(7)计算电路待测时刻电压增量信号△v：其中与步骤(2)中相同，即为无故障模式下的样本向量统计平均值；
(8)计算电路待测时刻电压增量信号△v的基频能量熵HT；
(9)计算待测电路健康指数HI(Health Index)：
(10)按如下规则判断电路健康情况：
若HI≥1,则判为健康；
若HI＜1,则判为电路出现故障，需要维修。
进一步，所述步骤(4)和(8)中，基频能量熵计算方法如下：
(a)将待处理信号赋给变量f，形成信号f：步骤(4)中，待处理信号为无故障模式F0下的容差电压增量信号i＝1,2,…M；步骤(8)中，待处理信号为电路待测时刻电压增量信号△v；
(b)确定信号f的基频信号的个数N：
(b-1)计算信号f的傅立叶变换f(ω)，找出f(ω)的J个幅度局部最大值及其所对应的频率点
(b-2)按照幅度递减的规律重排及其所对应的频率点即使L1≥L2≥…≥LJ；
(b-3)对数据集进行标准化处理：其中，
(b-4)设阈值若数据集中大于阈值Th的个数为Mth,则N＝Mth+1；
(c)采用经验小波变换技术求取信号f的基频信号fn,n＝0,1,2，…N-1,N为信号f的基频信号的个数；
(d)计算各基频信号fn的能量En：En＝<fn,fn>,n＝0,1,2,…N-1，其中<·>代表求内积运算；
(e)计算基频信号fn的能量概率密度pn：
(f)计算信号f的能量熵H：
进一步，所述步骤(c)中，采用经验小波变换技术求取信号f的基频信号fn,n＝0,1,2，…N-1的具体步骤为：
(c-1)对信号f的傅立叶变换f(ω)在频带[0,π]内的前N-1个局部幅度最大值所对应的频率ωn，n＝1,2,…,N-1按从小到大的顺序排列，令ω0＝0,ωN＝π，则以ωn，n＝0，1,2,…,N为边界,将频带[0,π]分割成N个连续的频谱段，第n个频谱段表示为Λn＝[ωn-1,ωn],n＝1,2,…N，且有其中∪表示求“并”运算；
(c-2)对每个频谱段边界，以ωn,n＝1,2,…,N-1为中心，划分过渡带区间(ωn-τn,ωn+τn),n＝1,2,…,N-1，则过渡带的宽带Tn＝2τn,n＝1,2,…,N-1，其中τn＝γωn,n＝1,2,…,N-1,γ取实数区间n＝1,2,…N-1上的任一常数,其中min(·)表示取最小值；
(c-3)确定经验小波变换的经验尺度函数和经验小波函数n＝1,2,…N-1：
式中：β(x)＝x 4(35-84x+70x2-20x3)；和为经验尺度函数频谱形式和经验小波函数频谱形式，即时域经验尺度函数的傅里叶变换后的形式和经验小波函数的傅里叶变换后的形式，此时ω为其频率，|·|为绝对值；
(c-4)计算经验小波变换系数：n＝1,2,…,N-1；
计算n＝0时的经验小波变换的近似系数
式中：f(ω)表示信号f的傅立叶变换，表示求横线下面变量的复共轭，F-1[·]表示求括号内信号的傅立叶逆变换；
(c-5)计算信号基波fn(t)，n＝0，1,2,…,N-1：
式中符号“*”表示卷积，和ψn(t)分别为经验尺度函数的时域形式和经验小波函数的时域形式。
本发明能在不丢失原始测量信息的前提下，利用EWT分解出测试信号基频信号，通过分别计算基频信号的能量熵与电路健康参数，能够及时对电路健康状况进行判断，进而对故障电路进行及时维修，在线实时性好，此外本发明通过计算基频信号的能量熵与电路健康参数来判断电路健康状况，计算方法简单，易于实现。
附图说明
图1为模拟电路健康预测方法系统框图；
图2为基频信号能量熵计算框图；
图3为基频信号个数确定流程。
具体实施方式
以下结合附图对本发明进行详细的说明。
参照图1，模拟电路健康预测方法，包括如下步骤：
(1)采集电路容差、无故障模式下的测试样本集：利用数据采集器获取待测模拟电路容差、无故障模式F0下的M组电压样本向量i＝1,2,…M,其中为采集到的无故障模式模式F0的第i组样本向量；
(2)计算电路容差、无故障模式F0下的M组电压样本向量的统计平均值i＝1,2,…M；
(3)计算电路容差、无故障模式F0下模拟电路测试样本容差电压增量信号集△V0：i＝1,2,…M，令则i＝1,2,…M，其中为第i个样本的容差电压增量信号；
(4)计算电路容差、无故障模式F0下的容差电压增量信号i＝1,2,…M的基频能量熵
(5)计算无故障容差电压增量信号的最大能量熵：其中max(·)代表取最大值；
(6)利用数据采集器采集电路待测时刻电压信号v；
(7)计算电路待测时刻电压增量信号△v：其中与步骤(2)中相同，即为无故障模式下的样本向量统计平均值；
(8)计算电路待测时刻电压增量信号△v的基频能量熵HT；
(9)计算待测电路健康指数HI：
(10)按如下规则判断电路健康情况：
若HI≥1,则判为健康；
若HI＜1,则判为电路出现故障，需要维修。
参照图2、图3，所述步骤(4)和(8)中，基频能量熵计算方法如下：
(a)将待处理信号赋给变量f，形成信号f：步骤(4)中，待处理信号为无故障模式F0下的容差电压增量信号i＝1,2,…M；步骤(8)中，待处理信号为电路待测时刻电压增量信号△v；
(b)确定信号f的基频信号的个数N：
(b-1)计算信号f的傅立叶变换f(ω)，找出f(ω)的J个幅度局部最大值及其所对应的频率点
(b-2)按照幅度递减的规律重排及其所对应的频率点即使L1≥L2≥…≥LJ；
(b-3)对数据集进行标准化处理：其中，
(b-4)设阈值若数据集中大于阈值Th的个数为Mth,则N＝Mth+1；
(c)采用经验小波变换技术求取信号f的基频信号fn,n＝0,1,2，…N-1,N为信号f的基频信号的个数；
(d)计算各基频信号fn的能量En：En＝<fn,fn>,n＝0,1,2,…N-1，其中<·>代表求内积运算；
(e)计算基频信号fn的能量概率密度pn：
(f)计算信号f的能量熵H：
进一步，所述步骤(c)中，采用经验小波变换技术求取信号f的基频信号fn,n＝0,1,2，…N-1的具体步骤为：
(c-1)对信号f的傅立叶变换f(ω)在频带[0,π]内的前N-1个局部幅度最大值所对应的频率ωn，n＝1,2,…,N-1按从小到大的顺序排列，令ω0＝0,ωN＝π，则以ωn，n＝0，1,2,…,N为边界,将频带[0,π]分割成N个连续的频谱段，第n个频谱段表示为Λn＝[ωn-1,ωn],n＝1,2,…N，且有其中∪表示求“并”运算；
(c-2)对每个频谱段边界，以ωn,n＝1,2,…,N-1为中心，划分过渡带区间(ωn-τn,ωn+τn),n＝1,2,…,N-1，则过渡带的宽带Tn＝2τn,n＝1,2,…,N-1，其中τn＝γωn,n＝1,2,…,N-1,γ取实数区间n＝1,2,…N-1上的任一常数,其中min(·)表示取最小值；
(c-3)确定经验小波变换的经验尺度函数和经验小波函数n＝1,2,…N-1：
式中：β(x)＝x4(35-84x+70x2-20x3)；和为经验尺度函数频谱形式和经验小波函数频谱形式，即时域经验尺度函数的傅里叶变换后的形式和经验小波函数的傅里叶变换后的形式，此时ω为其频率，|·|为绝对值；
(c-4)计算经验小波变换系数：n＝1,2,…,N-1；
计算n＝0时的经验小波变换的近似系数
式中：f(ω)表示信号f的傅立叶变换，表示求横线下面变量的复共轭，F-1[·]表示求括号内信号的傅立叶逆变换；
(c-5)计算信号基波fn(t)，n＝0，1,2,…,N-1：
式中符号“*”表示卷积，和ψn(t)分别为经验尺度函数的时域形式和经验小波函数的时域形式。