盛沛 许爱强 谷民亮

摘要：提出了一种基于局部特征尺度分解和基本尺度熵算法的模拟电路特征提取方法。首先采用LCD方法对被测信号进行分解，得到若干具有物理意义的分量;其次结合基本尺度熵算法得到特征向量。运用距离可分性指标的仿真结果表明，该方法提取的特征向量能够很好地描述电路的故障状态。

Abstract： A new method of feature vector extraction in analog circuit is proposed. Firstly， the measured signal is decomposed by using the LCD to obtain some signal components of physical significance. Secondly， combined with the base-scale entropy method， the feature vector are extracted. Simulation results show that the feature vectors extracted by this method can describe the fault state of the circuit well.

关键词：性能退化;模拟电路;局部特征尺度分解;基本尺度熵

Key words： performance degradation;analog circuit;LCD;base-scale entropy

中图分类号：N911;E917                                  文献标识码：A                                  文章编号：1006-4311（2019）25-0280-04

0  引言

传统的设备性能退化评估与预测流程主要由特征提取、退化评估、退化预测三大部分组成[1]。就目前电子装、设备发展现状[2]来看，往往是新型装设备在生产定型伊始即具备完善的测试手段，而老旧设备可测试性普遍不高，可达测试点不多、特征提取困难。因此，研究如何从有限的测试点中提取出丰富的特征数据并将其应用于诊断，就显得尤为重要。

局部特征尺度分解方法[3]-[5]（Local Characteristic-scale Decomposition，LCD）由湖南大学的程军圣团队提出。它可以把一个复杂的原始信号分解为n个ISC分量和一个剩余信号之和。该方法在一定程度上克服了传统时频分析方法的许多不足之处，并且在克服端点效应、计算时间等方面，较EMD[3]、LMD[4]等方法具有一定優势。

熵值能刻画时间序列的复杂程度，一般来说，时间序列复杂程度的增加会使熵值增加，当序列完全随机时熵值达到最大。基本尺度熵（Base-scale Entropy，BE）是文献[6]提出的一种复杂性度量方法，克服了近似熵、样本熵及排列熵的缺点，具有一致性好、适用于短时间序列、计算简便、抑制噪声能力强等优点。

基于此，本文根据局部特征尺度分解和基本尺度熵两个算法的特性，将二者相结合来提取模拟电路故障特征。更进一步的，引入距离可分性指标验证了该方法的有效性。

1  算法基本原理

1.1 LCD理论

LCD方法的原理是基于单分量信号“局部关于零均值对称”的结论，其基础在于构造均值曲线，并将均值曲线不断从原始信号中分离，直至信号为内禀尺度分量[5]。与传统的HHT方法不同的是，在LCD方法中，局部均值曲线定义为“基线信号”[8]。如图 1所示，基线信号的构造方法为：

2  基于LCD-基本尺度熵的故障特征提取方法

基于LCD-基本尺度熵的故障特征提取方法如下：

步骤1：首先，对原始信号进行LCD，得到多个ISC分量;

步骤2：选取包含主要故障特征的ISC分量作为有用分量;

步骤3：计算每个ISC分量的基本尺度熵，构成特征向量。

从文献[11]-[13]的有用分量选取过程中可以看出：有用分量往往是前几个ISC分量，这是由于故障信息集中在高频区，前几个分量无论是幅值、频率还是相关性指标均比较大。为了简化计算，本文借鉴文献[14]推荐的方法，选取了前4个ISC分量作为有用分量。

3  仿真实验

选用某型雷达发射机模拟电路中开关电源部分[15]作为验证对象。仿真平台ORCAD CIS和MATLAB16.2。电路图如图 2所示，节点OUT为输出测试点。状态集如表 1所示。

验证LCD-基本尺度熵的故障特征提取方法。选取正常状态20组数据（标称值5%容差范围内）、R1↑（超过标称值50%，下同。）故障状态40组数据、C1↓状态40组数据、L1↑状态40组数据，共140组数据。经LCD后的某组典型仿真数据分解结果如图 3所示。分别计算原始数据的基本尺度熵和前4个ISC的基本尺度熵，其均值见表2。

从表 2可以看出：

四种状态下各分量熵值发生一定的变化，总体趋势下降。这是因为随着时间尺度的变化，信号复杂程度变小。另外可以看出，在每一个时间尺度下四种状态的熵值大小关系与原始信号的熵值大小关系基本一致。因此，通过LCD将数据多尺度化后，将基本尺度熵扩维，拓展了熵值的作用区域，即可使故障特征更为突出。

为了进一步验证LCD-基本尺度熵的性能，参照文献[16]的近似熵算法，并引入距离可分性指标：

式中：Sw为样本的类内矩阵，表征了类内散度，其值越小，类别的可分性越好;Sb为样本类间散布矩阵，表征了类间散度，其值越大，类别的可分性越好。Ji为 的特征值，综合了特征向量的敏感度，Jj值越大，说明类别的可分性越好。

图 4给出了原始信号及ISC1～ISC5分量基本尺度熵和近似熵的距离可分性指标。从图中可以看出：

①基本尺度熵的各特征量的Ji值都要高于近似熵，因此，通过LCD对基本尺度熵升维的特征向量能够更好的描述电路的故障状态。

②随着各分量的变化，距离可分性指标在逐步变小，也就是说，虽然LCD可得到较多的分量，但是后面的分量对分类的贡献率不够明显，选取过多的分量计算熵值不合理。

4  结束语

本文给出了LCD-基本尺度熵的故障特征提取方法。在通过仿真实验验证了其分类性能优于LCD-近似熵的同时，说明了分量选取的方法，验证了所提出的特征向量提取方法适用于性能退化型模拟电路故障诊断。

参考文献：

[1]郭磊，陈进.设备性能退化评估与预测研究综述[J].振动与冲击，2008，27（S）：139-141.

[2]唐玉志，赵生伟.军工自动测试技术和设备现状及发展趋势[J].数据采集与处理，2009，24（S）：200-204.

[3]Huang N E， Shen Z， Long S R， et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of the Royal Society of London. Series A： Mathematical， Physical and Engineering Sciences， 1998， 454（1971）： 903-995.

[4]Smith J S. The local mean decomposition and its application to EEG perception data[J]. Journal of the Royal Society Interface， 2005， 2（5）： 443-454.

[5]程军圣，郑近德，杨宇.一种新的非平稳信号分析方法—局部特征尺度分解法[J].振动工程学报，2012，25（2）：215-220.

[6]李锦，宁新宝.短时心率变异性信号的基本尺度熵分析[J]. 科学通报，2005，50（14）：1438-1441.

[7]杨宇，曾鸣，程军圣.一种新的时频分析方法—局部特征尺度分解[J]. 湖南大学学报（自然科学版），2012，39（6）：35-39.

[8]杨宇，曾鸣，局部特征尺度分解方法及其分解能力研究[J]. 振動工程学报，2012，25（5）：602-609.

[9]程军圣，李海龙，杨宇.基于内禀尺度分量的自适应时频分析方法[J].中南大学学报，2013，44（4）：1425-1430.

[10]刘铁兵，姚文坡，宁新宝.功能磁共振成像的基本尺度嫡分析[J].物理学报，2013，62（21）：1-6.

[11]张俊红，李林洁，马文朋，等.EMD-ICA联合降噪在滚动轴承故障诊断中的应用[J].中国机械工程，2013，24（11）：1468-1472.

[12]赵阳，陈捷，洪荣晶，等.EMD和独立分量分析在转盘轴承故障诊断中的应用[J].轴承，2015，7：54-59.

[13]崔玲丽，吴春光，邬娜.基于EMD与ICA的滚动轴承复合故障诊断[J].北京工业大学学报，2014，40（10）：1459-1464.

[14]郑近德.局部特征尺度分解及其在机械诊断中的应用研究[D].长沙：湖南大学，2014.

[15]王刘军，刘昌锦.基于Hilbert-Huang变换的雷达发射机故障特征提取[J].兵工自动化，2013，32（8）：68-70.

[16]单剑锋，万国发.一种LMD和近似熵算法的模拟电路特征提取方法[J].机械科学与技术，2018，5（9）：78-83.