CN107894564A - 一种基于交叉小波特征的模拟电路故障诊断方法 - Google Patents

一种基于交叉小波特征的模拟电路故障诊断方法 Download PDF

Info

Publication number
CN107894564A
CN107894564A CN201711097153.2A CN201711097153A CN107894564A CN 107894564 A CN107894564 A CN 107894564A CN 201711097153 A CN201711097153 A CN 201711097153A CN 107894564 A CN107894564 A CN 107894564A
Authority
CN
China
Prior art keywords
mrow
msub
training sample
mover
sample set
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201711097153.2A
Other languages
English (en)
Other versions
CN107894564B (zh
Inventor
何怡刚
何威
李志刚
佐磊
李兵
何鎏璐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201711097153.2A priority Critical patent/CN107894564B/zh
Priority to PCT/CN2017/115011 priority patent/WO2019090879A1/zh
Priority to US16/762,926 priority patent/US11409933B2/en
Publication of CN107894564A publication Critical patent/CN107894564A/zh
Application granted granted Critical
Publication of CN107894564B publication Critical patent/CN107894564B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/28Testing of electronic circuits, e.g. by signal tracer
    • G01R31/316Testing of analog circuits

Abstract

一种基于交叉小波特征的模拟电路故障诊断方法,包括以下步骤:向被测模拟电路输入激励信号,采集时域响应输出信号,构成原始数据样本集;将原始数据样本集分成训练样本集和测试样本集;将训练样本集和测试样本集,均进行交叉小波分解,分别获取训练样本集和测试样本集的小波交叉谱;应用双向二维线性判别分析对训练样本集和测试样本集的小波交叉谱进行处理,提取训练样本集和测试样本集的故障特征向量;将训练样本集的故障特征向量提交给支持向量机用于训练SVM分类器,构建支持向量机故障诊断模型;向该模型输入测试样本集的故障特征向量,进行故障分类。使用本发明,能高效地识别模拟电路故障,并显著提高模拟电路故障诊断的精度。

Description

一种基于交叉小波特征的模拟电路故障诊断方法
技术领域
本发明属于电路测试技术领域,具体涉及一种基于交叉小波特征的模拟电路故障诊断方法。
背景技术
随着微电子技术的发展,大规模数模混合集成电路在电子产品、工业控制、通讯设备等领域中发挥着广泛的应用。模拟电路作为集成电路的重要组成部分,其对信号所进行的传输、滤波、放大、转换是系统正常运行的必要功能。因此,极有必要对模拟电路开展故障诊断技术的研究,以保障系统的可靠性。
然而,由于模拟电路故障模拟难以构建、电路元器件的非线性和容差的影响,模拟电路故障诊断仍然面临着诸多问题。例如,外界噪声的干扰、故障发生未知的不可预测性及参数偏离正常范围程度的不确定性,都将极大的影响故障诊断方法的准确性及效率。
针对上述难点,国内外诸多学者均对其展开了大量的研究工作。这些研究主要涉及故障特征提取方法。例如,Spina R,Upadhyaya S.Linear circuit fault diagnosisusing neuromorphic analyzers[J].Circuits&Systems II Analog & Digital SignalProcessing IEEE Transactions on,1997,44(3):188-196.Negnevitsky M,PavlovskyV.Neural networks approach to online identification of multiple failures ofprotection systems[J].IEEE Transactions on Power Delivery,2005,20(2):588-594.均直接将被测电路的故障响应信号不加任何预处理作为分类器的输入,导致分类器训练时间过长,结构规模过大。又如,何星,王宏力,陆敬辉,等.基于优选小波包和ELM的模拟电路故障诊断[J].仪器仪表学报,2013,34(11):2614-2619.肖玉飞,刘祖润,李目.基于小波包能量熵与SVM的模拟电路故障诊断[J].电子测量技术,2011,34(6):110-113.是通过计算小波系数能量值或能量熵,并将其作为故障特征输入到分类器以实现故障识别,但各节点特征值较小且特征区分不明显。又如,Yuan L,He Y,Huang J,et al.A New Neural-Network-Based Fault Diagnosis Approach for Analog Circuits by Using Kurtosis andEntropy as a Preprocessor[J].IEEE Transactions on Instrumentation&Measurement,2010,59(3):586-595.Long Y,Xiong Y,He Y,et al.A new switchedcurrent circuit fault diagnosis approach based on pseudorandom test andpreprocess by using entropy and Haar wavelet transform[J].Analog IntegratedCircuits&Signal Processing,2017(3):1-17.作者采用故障响应信号的高阶统计量作为故障特征,输入到分类器进行识别,并获得了很好的诊断效果。但是没有考虑多故障情况下即故障类别重叠性较大时,故障诊断方法的性能。
综上,有必要研究一种可以有效提高电路故障诊断效率,并且能够诊断多种故障类型的模拟电路故障诊断方法。
发明内容
本发明要解决的技术问题是,克服现有技术存在的上述缺陷,提供一种可以有效提高电路故障诊断效率,并且能够诊断多种故障类型的基于交叉小波特征的模拟电路故障诊断方法。
本发明解决其技术问题所采用的技术方案是:
一种基于交叉小波特征的模拟电路故障诊断方法,包括以下步骤:
步骤一:向被测模拟电路输入激励信号,采集正常状态下的时域响应输出信号和故障状态下的时域响应输出信号共同构成一个原始数据样本集。此处,原始数据样本集中的正常时域响应输出信号和故障时域响应输出信号,按照采集顺序依次排列在原始数据样本集中;
然后,从原始数据样本集内的正常时域响应输出信号和故障时域响应输出信号中,各选择50%构成训练样本集,将剩余的50%构成测试样本集。
该步骤可通过Multisim软件对被测模拟电路进行仿真来实现。
步骤二:将训练样本集和测试样本集,均进行交叉小波(Cross waveletTransform,XWT)分解,分别获取训练样本集的小波交叉谱(Wavelet Cross Spectrum,WCS)和测试样本集的小波交叉谱。其中:小波交叉谱以二维矩阵形式储存于仿真软件matlab内存中。
步骤三:应用双向二维线性判别(Bilateral Two-Dimensional LinearDiscriminant Analysis,B2DLDA)分析对训练样本集的小波交叉谱进行处理,提取训练样本集的故障特征向量;应用双向二维线性判别分析对测试样本集的小波交叉谱进行处理,提取测试样本集的故障特征向量。
步骤四:将训练样本集的故障特征向量提交给支持向量机(Support VectorMachine,SVM),用于训练SVM分类器,构建基于训练样本特征向量的支持向量机(supportvector machine,SVM)故障诊断模型。待模型训练完成后,向该模型输入测试样本集的故障特征向量,进行故障分类。
进一步,所述步骤二中,训练样本集的小波交叉谱的获取方法为:
首先,将训练样本集内的电路故障状态下的时域响应输出信号及电路正常状态下的时域响应输出信号均进行小波分解以获取各自的信号小波系数,可表示为:
式中:t为时间点;x(t)为训练样本集内的电路故障状态下的时域响应输出信号,y(t)为训练样本集内的电路正常状态下的时域响应输出信号;a为尺度参数;τ为时移参数;*表示复数共轭。ψ是小波基函数。本发明中选择Morlet函数为基函数,其表达式为
其次,对训练样本集内的电路正常状态下的时域响应输出信号及电路故障状态下的时域响应输出信号进行交叉处理,以获得训练样本集的小波交叉谱Sx,y(a,τ)。该过程可表示为:
同理,可获得测试样本集的小波交叉谱。
训练样本集的小波交叉谱和测试样本集的小波交叉谱均为矩阵形式。
应用双向二维线性判别分析对训练样本集的小波交叉谱和测试样本集的小波交叉谱进行降维处理,去除原始矩阵(即训练样本集的小波交叉谱和测试样本集的小波交叉谱)行列上的信息冗余,减少特征维数,提取最有利于分类的故障特征数据。
进一步,所述步骤三中,应用双向二维线性判别分析对训练样本集的小波交叉谱进行处理,提取训练样本集的故障特征向量,具体方法为:
首先对矩阵(即训练样本集的小波交叉谱)行方向上进行降维。设在训练样本集中有C类模式,M为训练样本集内所有训练样本的总个数(即训练样本集内所有的时域响应输出信号的总个数);Mi是训练样本集中第i类训练样本的个数(即训练样本集中第i故障类中时域响应输出信号的个数)设m×n矩阵Aij表示第i类训练样本中的第j个元素。则定义类间离散度矩阵Sb和类内离散度矩阵Sw为:
式中,是训练样本集中第i类训练样本的平均值;是训练样本集中所有训练样本的总平均值;
然后求解矩阵的最大的d个特征值和对应的d个正交特征向量:x1,x2,…,xd。正交特征向量x1,x2,…,xd为训练样本集的小波交叉谱在行方向上的最佳投影。此处,d≤(C-1)。d的具体值,由使用者在进行实例验证时,自行确定。
由正交特征向量x1,x2,…,xd构成的向量集即为一个n×d的矩阵X=[x1,x2,…,xd]。则给定一个m×n的样本矩阵A,有
B=AX;
其中,B为矩阵A的特征矩阵,是一个m×d的矩阵。
其次,对矩阵列方向上应用二维线性判别(2DLDA)方法进行降维,同理可得如下公式:
然后求解矩阵的最大的f个特征值和对应的f个正交特征向量:将正交特征向量作为训练样本集的小波交叉谱在列方向上的最佳投影。此处,f≤(C-1)。由正交特征向量构成的向量集即为一个m×f的矩阵f的具体值,由使用者在进行实例验证时,自行确定。
综上,则样本矩阵A同时在行方向上和列方向上投影:
其中,Z即是样本矩阵A利用双向2DLDA方法提取出的特征矩阵,维数为f×d。
将上式提取到的特征矩阵Z,转换为一维向量,从而构成训练样本集的故障特征向量。
同理,应用双向二维线性判别分析对测试样本集的小波交叉谱进行处理,可提取测试样本集的故障特征向量。
步骤四中,如文献“孙永奎,陈光,李辉.基于自适应小波分解和SVM的模拟电路故障诊断[J].仪器仪表学报,2008,29(10):2105-2109”中所述,首先构造基于训练样本特征向量的支持向量机(support vector machine,SVM)故障诊断模型。SVM模型中惩罚参数及核参数通过网格搜索法寻优得到。待模型构造完毕,输入测试样本集的故障特征向量,从而实现对故障类别进行识别。
本发明应用交叉小波处理被测电路输出响应信号,输出信号经交叉小波变换后,生成为小波交叉谱,以二维矩阵方式储存,小波交叉谱难以直接应用于训练及测试故障诊断模型。本发明提出应用二维线性判别分析方法对小波交叉谱进行特征向量提取;以便于后续训练及测试故障诊断模型。
采用本发明所获得的有益效果有:
本发明通过交叉小波变换较为全面地获取被测电路的故障信号与正常信号中的有用信息;相比其他信号处理技术,通过应用交叉小波变换所获得的小波交叉谱以图像的方式直观全面地表征信号的时间-频率-能量信息。B2DLDA技术可以直接对图像矩阵进行特征提取,实现各故障特征之间明显区分。通过SVM自动实现快速准确地故障分类。
对电路的软故障及多故障的仿真实验表明:本发明方法模拟电路的故障响应与正常响应特征之间的区别,以及各故障响应特征之间的区别较为明显,特征向量充分涵盖了各类故障本质信息;通过SVM所进行的故障诊断具备准确、迅速的效果。
附图说明
图1是本发明故障诊断方法的体系结构示意图;
图2是四运放低通滤波器电路图;
图3是四运放低通滤波器电路的故障特征图;
图4是四运放低通滤波器的故障分类图。
具体实施方式
下面结合附图和具体实例对本发明作进一步详细说明。
1.故障诊断方法
如图1所示,一种基于交叉小波特征的模拟电路故障诊断方法,具体步骤如下:
步骤一:向被测模拟电路输入激励信号,采集正常状态下的时域响应输出信号和故障状态下的时域响应输出信号共同构成一个原始数据样本集。此处,原始数据样本集中的正常时域响应输出信号和故障时域响应输出信号,按照采集顺序依次排列在原始数据样本集中;然后,从原始数据样本集内的正常时域响应输出信号和故障时域响应输出信号中,各选择50%构成训练样本集,将剩余的50%构成测试样本集。
也就是说,通过Multisim仿真,采集被测模拟电路输出响应信号。
步骤二:将训练样本集和测试样本集,均进行交叉小波(Cross waveletTransform,XWT)分解,分别获取训练样本集的小波交叉谱(Wavelet Cross Spectrum,WCS)和测试样本集的小波交叉谱。其中:小波交叉谱以二维矩阵形式储存于仿真软件matlab内存中。
步骤三:应用双向二维线性判别(Bilateral Two-Dimensional LinearDiscriminant Analysis,B2DLDA)分析对训练样本集的小波交叉谱进行处理,提取训练样本集的故障特征向量;应用双向二维线性判别分析对测试样本集的小波交叉谱进行处理,提取测试样本集的故障特征向量。
本实施例中,设置B2DLDA在列方向上保留维数为3,在行方向上保留维数为2;因此,对每个故障信号经过特征提取后,每个样本相应的对应一个6维特征向量。
步骤四:将训练样本集的故障特征向量提交给支持向量机(Support VectorMachine,SVM),用于训练SVM分类器,构建基于训练样本特征向量的支持向量机(supportvector machine,SVM)故障诊断模型。待模型训练完成后,向该模型输入测试样本集的故障特征向量,进行故障分类。
下面将本发明故障诊断方法中的核心技术:交叉小波变换、双向二维线性判别分析作进一步详细说明。
1.1交叉小波分析
对于能量有限信号x(t),其小波变换为
式中:a为尺度参数;τ为时移参数;*表示复数共轭。ψ是小波基函数。本发明中选择Morlet函数为基函数,其表达式为
对于2个能量有限信号x(t)和y(t),其交叉小波变换(cross-wavelet transform,XWT)为
在实际应用中,设置信号x(t)为电路故障信号即电路故障状态下的时域响应输出信号,y(t)为电路正常状态下响应信号即电路正常状态下的时域响应输出信号。
本发明中采用交叉小波谱进行对故障信号进行分析,其定义为:
小波变换将一维信号在时域和频域中展开,可以分析信号变换的任何细节。在小波变换基础之上建立的交叉小波变换得到2个信号相关性在时频域中分布状况,其变换系数表示这2个信号在时频域中存在相关性的大小,其值越大,说明相关性越密切。此处,本发明选择Morlet函数作为小波基函数。
1.2双向二维线性判别分析(B2DLDA)
由于经交叉小波变换得到的小波交叉谱为包含复值的矩阵,不能直接提交给SVM分类器,以实现故障诊断。因此还需对小波交叉谱进一步处理,即采用B2DLDA提取WCS数据集故障特征。
B2DLDA通过对矩阵进行2次降维处理,去除原始矩阵行列上的信息冗余,减少特征维数,提取最有利于分类的特征数据。
假设在训练样本集中有C类模式,M为训练样本集内采集的时域响应输出信号的总个数即训练样本集内所有训练样本的总个数;Mi是训练样本集中第i故障类中响应输出信号个数即训练样本集中第i类训练样本的个数设m×n矩阵Aij表示第i类训练样本中的第j个元素。则定义类间离散度矩阵Sb和类内离散度矩阵Sw为:
式中,是训练样本集中第i类训练样本的平均值;是训练样本集中所有训练样本的总平均值;
然后求解矩阵的最大的d个特征值和对应的d个特征向量:x1,x2,…,xd。正交特征向量x1,x2,…,xd为训练样本集的小波交叉谱在行方向上的最佳投影。此处,d≤(C-1)。d的具体值,由使用者在进行实例验证时,自行确定。
其中,由正交特征向量x1,x2,…,xd构成的向量集即为一个n×d的矩阵X=[x1,x2,…,xd]。则给定一个m×n的样本矩阵A,有
B=AX;
其中,B为m×d的矩阵,表示矩阵A的特征矩阵。
其次,利用2DLDA方法将矩阵在列方向进行投影降维,同理可得如下公式:
然后求解矩阵的最大的f个特征值和对应的f个特征向量:将正交特征向量作为训练样本集的小波交叉谱在列方向上的最佳投影。此处,f≤(C-1)。由正交特征向量构成的向量集即为一个m×f的矩阵f的具体值,由使用者在进行实例验证时,自行确定。
综上,则样本矩阵A在行方向和列方向上的最佳投影:
其中,Z即是矩阵A利用双向2DLDA方法提取出特征矩阵,维数仅为f×d。
通过上式提取到的特征矩阵,转换为一维向量,从而构成了训练样本集的电路故障特征向量。
同理,通过上述方法可提取测试样本集的故障特征向量。
2.案例电路与方法应用:
下面以两个典型故障电路为例,分别对具有硬故障的视频放大电路和具有多故障的四运放低通滤波器电路应用本发明方法进行故障诊断作说明。
2.1故障电路:
案例:本发明所采用的故障电路为如图2所示的四运放高通滤波器电路,各元件的标称值均于图上标出。在此,以此电路为例展示本发明提出的故障诊断方法的整个流程,激励源采用持续时间为10us,幅值为10V的脉冲波,故障时域响应信号在电路输出端采样获得。设定该电路中电阻与电容分别为5%。
选择R1,R2,R3,R4,C1和C2作为测试对象,表1给出了电路中各测试元件的故障码,故障类别,标称值和故障值,其中↑和↓分别表示高出和低于标称值;此处,正常状态也视为故障状态,正常状态的故障码为F0。为每种故障类别(正常状态及故障状态)分别采样60个数据,分成2部分,前30个用于建立基故障诊断模型,后30个数据用于测试该诊断模型的性能。
表1
故障码 故障类别 故障值 故障码 故障类别 故障值
F1 R1↓ 3000Ω F2 R1↑ 15000Ω
F3 R2↓ 2000Ω F4 R2↑ 18000Ω
F5 R3↓ 2700Ω F6 R3↑ 12000Ω
F7 R4↓ 500Ω F8 R4↑ 2500Ω
F9 C1↓ 2.5nF F10 C1↑ 10nF
F11 C2↓ 1.5nF F12 C2↑ 15nF
F13 C1↑C2↓ - F14 R1↑C2↓ -
F15 R1↑R2↓ - F16 R2↑R4↓ -
步骤一:
在四运放双二阶高通滤波器中,施加的激励响应为:幅值为10V,持续时间为10us的脉冲序列。电路在不同故障模式下的输出响应通过Multisim仿真。采集时域响应输出信号,以构成原始数据样本集,将原始数据样本集平均分成两部分,训练样本集及测试样本集。
步骤二:
下面采用交叉小波变换对电路输出响应进行分析以获得训练样本集的交叉小波谱矩阵和测试样本集的交叉小波谱矩阵。其中,交叉小波的小波基选用复morlet小波。鉴于交叉小波谱矩阵为复负矩阵,不便于进一步处理。本发明中,对矩阵中复数元素求幅值后,组成新的交叉谱矩阵。
步骤三:
对矩阵进行双向二维线性判别分析,提取训练样本集的故障特征向量和测试样本集的故障特征向量。分别对矩阵在行方向和列方向上进行降维。此处,设置B2DLDA在列方向上保留维数为f=3,在行方向上保留维数为d=2;故,每个样本对应的特征向量为6维。
不同故障之间或故障与正常状态之间的特征数值差别越大,不同故障之间或故障与正常状态之间的信号响应差别就越明显,这样的故障特征就越有利于故障诊断。从图3可知,按本发明方法所得到的电路故障响应特征与电路正常状态响应特征之间的数值差异,以及电路在不同故障模式下的特征数值差异明显,这充分说明本发明故障特征提取的有效性。
步骤四:
将训练样本集的的故障特征向量输入到SVM中以训练分类模型,待训练完成后,将测试样本集的故障特征向量输入到SVM分类器模型中,故障诊断的结果如图所示4。SVM模型成功地识别了绝大部分的故障,故障整体准确率达到99.22%。

Claims (3)

1.一种基于交叉小波特征的模拟电路故障诊断方法,其特征在于,包括以下步骤:
步骤一:向被测模拟电路输入激励信号,采集正常状态下的时域响应输出信号和故障状态下的时域响应输出信号共同构成一个原始数据样本集;此处,原始数据样本集中的正常时域响应输出信号和故障时域响应输出信号,按照采集顺序依次排列在原始数据样本集中;
然后,从原始数据样本集内的正常时域响应输出信号和故障时域响应输出信号中,各选择50%构成训练样本集,将剩余的50%构成测试样本集;
步骤二:将训练样本集和测试样本集,均进行交叉小波分解,分别获取训练样本集的小波交叉谱和测试样本集的小波交叉谱;其中:小波交叉谱以二维矩阵形式存在;
步骤三:应用双向二维线性判别分析对训练样本集的小波交叉谱进行处理,提取训练样本集的故障特征向量;应用双向二维线性判别分析对测试样本集的小波交叉谱进行处理,提取测试样本集的故障特征向量;
步骤四:将训练样本集的故障特征向量提交给支持向量机用于训练SVM分类器,构建基于训练样本特征向量的支持向量机故障诊断模型;待模型训练完成后,向该模型输入测试样本集的故障特征向量,进行故障分类。
2.根据权利要求1所述的基于交叉小波特征的模拟电路故障诊断方法,其特征在于,所述步骤二中,训练样本集的小波交叉谱的获取方法为:
首先,将训练样本集内的电路故障状态下的时域响应输出信号及电路正常状态下的时域响应输出信号均进行小波分解以获取各自的信号小波系数,表示为:
<mrow> <msub> <mi>W</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mi>a</mi> </msqrt> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </msubsup> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mi>&amp;psi;</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>&amp;tau;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>,</mo> <mi>a</mi> <mo>&gt;</mo> <mn>0</mn> <mo>;</mo> </mrow>
<mrow> <msub> <mi>W</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mi>a</mi> </msqrt> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </msubsup> <mi>y</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mi>&amp;psi;</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>&amp;tau;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>,</mo> <mi>a</mi> <mo>&gt;</mo> <mn>0</mn> <mo>;</mo> </mrow>
式中:t为时间点;x(t)为训练样本集内的电路故障状态下的时域响应输出信号,y(t)为训练样本集内的电路正常状态下的时域响应输出信号;a为尺度参数;τ为时移参数;*表示复数共轭;ψ是小波基函数;选择Morlet函数为基函数,其表达式为
其次,对训练样本集内的电路正常状态下的时域响应输出信号及电路故障状态下的时域响应输出信号进行交叉处理,以获得训练样本集的小波交叉谱Sx,y(a,τ);该过程表示为:
<mrow> <msub> <mi>C</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>W</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <msubsup> <mi>W</mi> <mi>y</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>S</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>C</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <msubsup> <mi>C</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>C</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>;</mo> </mrow>
同理,可获得测试样本集的小波交叉谱;
训练样本集的小波交叉谱和测试样本集的小波交叉谱均为矩阵形式。
3.根据权利要求1或2所述的基于交叉小波特征的模拟电路故障诊断方法,其特征在于,所述步骤三中,应用双向二维线性判别分析对训练样本集的小波交叉谱进行处理,提取训练样本集的故障特征向量,具体方法为:
首先对矩阵即训练样本集的小波交叉谱行方向上进行降维;设在训练样本集中有C类模式;M为训练样本集内所有训练样本的总个数即训练样本集内所有时域响应输出信号的总个数;Mi是训练样本集中第i类训练样本的个数即训练样本集中第i故障类中时域响应输出信号的个数设m×n矩阵Aij表示第i类训练样本中的第j个元素;则定义类间离散度矩阵Sb和类内离散度矩阵Sw为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mi>b</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <msub> <mi>M</mi> <mi>i</mi> </msub> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mi>w</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
式中,是训练样本集中第i类训练样本的平均值;是训练样本集中所有训练样本的总平均值;
然后求解矩阵的最大的d个特征值和对应的d个正交特征向量:x1,x2,…,xd;正交特征向量x1,x2,…,xd为训练样本集的小波交叉谱在行方向上的最佳投影;此处,d≤(C-1);
由正交特征向量x1,x2,…,xd构成的向量集即为一个n×d的矩阵X=[x1,x2,…,xd];则给定一个m×n的样本矩阵A,有
B=AX;
其中,B为矩阵A的特征矩阵,是一个m×d的矩阵;
其次,对矩阵列方向上应用二维线性判别(2DLDA)方法进行降维,同理得如下公式:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>S</mi> <mo>~</mo> </mover> <mi>b</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <msub> <mi>M</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mi>T</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>S</mi> <mo>~</mo> </mover> <mi>w</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>c</mi> </munderover> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
然后求解矩阵的最大的f个特征值和对应的f个正交特征向量:将正交特征向量作为训练样本集的小波交叉谱在列方向上的最佳投影;此处,f≤(C-1);由正交特征向量构成的向量集即为一个m×f的矩阵
综上,则样本矩阵A同时在行方向上和列方向上投影:
<mrow> <mi>Z</mi> <mo>=</mo> <msup> <mover> <mi>X</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mi>A</mi> <mi>X</mi> <mo>;</mo> </mrow>
其中,Z即是样本矩阵A利用双向2DLDA方法提取出的特征矩阵,维数为f×d;
将上式提取到的特征矩阵Z,转换为一维向量,从而构成训练样本集的故障特征向量;
同理,应用双向二维线性判别分析对测试样本集的小波交叉谱进行处理,可提取测试样本集的故障特征向量。
CN201711097153.2A 2017-11-09 2017-11-09 一种基于交叉小波特征的模拟电路故障诊断方法 Active CN107894564B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201711097153.2A CN107894564B (zh) 2017-11-09 2017-11-09 一种基于交叉小波特征的模拟电路故障诊断方法
PCT/CN2017/115011 WO2019090879A1 (zh) 2017-11-09 2017-12-07 一种基于交叉小波特征的模拟电路故障诊断方法
US16/762,926 US11409933B2 (en) 2017-11-09 2017-12-07 Method for diagnosing analog circuit fault based on cross wavelet features

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711097153.2A CN107894564B (zh) 2017-11-09 2017-11-09 一种基于交叉小波特征的模拟电路故障诊断方法

Publications (2)

Publication Number Publication Date
CN107894564A true CN107894564A (zh) 2018-04-10
CN107894564B CN107894564B (zh) 2020-02-18

Family

ID=61804719

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711097153.2A Active CN107894564B (zh) 2017-11-09 2017-11-09 一种基于交叉小波特征的模拟电路故障诊断方法

Country Status (3)

Country Link
US (1) US11409933B2 (zh)
CN (1) CN107894564B (zh)
WO (1) WO2019090879A1 (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108629348A (zh) * 2018-05-09 2018-10-09 江南大学 一种基于显著性检验的机械设备故障诊断方法
CN109239585A (zh) * 2018-09-06 2019-01-18 南京理工大学 一种基于改进优选小波包的故障诊断方法
CN109613610A (zh) * 2019-01-17 2019-04-12 中南大学 微震信号到时差的自动拾取方法
CN109782158A (zh) * 2019-02-19 2019-05-21 桂林电子科技大学 一种基于多级分类的模拟电路诊断方法
CN110991481A (zh) * 2019-10-31 2020-04-10 河海大学 一种基于交叉小波变换的高压并联电抗器内部松动故障诊断方法
CN111337819A (zh) * 2020-03-02 2020-06-26 武汉大学 基于优选小波基函数的模拟电路故障特征提取方法及系统
CN111626329A (zh) * 2020-04-20 2020-09-04 国网上海市电力公司 一种基于lda寻优多尺度纹理特征的绝缘管母故障诊断方法
CN112101116A (zh) * 2020-08-17 2020-12-18 北京无线电计量测试研究所 一种基于深度学习的模拟电路故障诊断方法
CN112819037A (zh) * 2021-01-12 2021-05-18 广东石油化工学院 基于交叉注意力和自注意力的分类参数分布的故障诊断方法
CN113985261A (zh) * 2021-07-07 2022-01-28 四川大学锦城学院 一种基于神经网络的模拟电路故障诊断方法

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110991240B (zh) * 2019-10-31 2023-06-09 湖南世优电力科技股份有限公司 一种低压中性点不接地系统的单相故障识别及分类方法
CN111339872A (zh) * 2020-02-18 2020-06-26 国网信通亿力科技有限责任公司 一种基于分类模型的电网故障分类方法
CN111898447B (zh) * 2020-06-30 2024-02-02 西安理工大学 基于辛几何模态分解的风电机组故障特征提取方法
US11182525B1 (en) * 2020-07-07 2021-11-23 Infineon Technologies Ag Fault aware analog model (FAAM)
CN112712046B (zh) * 2021-01-06 2023-06-16 浙江大学 一种基于设备硬件指纹的无线充电设备认证方法
CN113011084B (zh) * 2021-02-26 2023-08-11 华北电力大学 基于相关向量机的输电线路故障原因智能辨识方法
CN113255777B (zh) * 2021-05-28 2023-01-31 郑州轻工业大学 基于多模态敏感特征选取融合的设备故障预警方法及系统
CN113311407B (zh) * 2021-05-28 2024-01-19 西安电子科技大学 基于支持向量机的空中飞机目标分类在线学习方法
CN113434973B (zh) * 2021-06-21 2023-04-07 西安科技大学 基于改进pso-som-bpnn的机械故障诊断方法
CN113657556B (zh) * 2021-09-23 2023-12-26 华北电力大学 基于多元统计分析的燃气轮机进口导叶系统故障诊断方法
CN114062845A (zh) * 2021-11-11 2022-02-18 国网北京市电力公司 线路故障的检测方法、装置以及电子设备
CN114118292B (zh) * 2021-12-06 2024-04-16 东北大学 一种基于线性判别邻域保持嵌入的故障分类方法
CN114266401B (zh) * 2021-12-23 2022-08-05 北京石油化工学院 一种底盘发动机故障预测方法及系统
CN115688018B (zh) * 2023-01-04 2023-09-08 湖南大学 一种多工况下轴承的状态监测与故障诊断方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101221213A (zh) * 2008-01-25 2008-07-16 湖南大学 基于粒子群算法的模拟电路故障诊断神经网络方法
CN101329697A (zh) * 2008-06-11 2008-12-24 电子科技大学 一种基于混杂算法的模拟电路状态预测方法
CN102636742A (zh) * 2012-05-15 2012-08-15 长沙河野电气科技有限公司 基于小波神经网络的大规模模拟电路故障诊断方法
CN103824135A (zh) * 2014-03-11 2014-05-28 合肥工业大学 一种模拟电路故障预测方法
CN104198924A (zh) * 2014-09-11 2014-12-10 合肥工业大学 一种新颖的模拟电路早期故障诊断方法
CN105701470A (zh) * 2016-01-13 2016-06-22 合肥工业大学 一种基于最优小波包分解的模拟电路故障特征提取方法
CN107101715A (zh) * 2017-05-25 2017-08-29 贵州电网有限责任公司 一种基于交叉小波的变压器振动信号幅频特征量提取方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100589118C (zh) * 2007-11-29 2010-02-10 北京中星微电子有限公司 特征提取方法、装置与模式识别方法、装置
CN102982349B (zh) * 2012-11-09 2016-12-07 深圳市捷顺科技实业股份有限公司 一种图像识别方法及装置
CN105044590B (zh) * 2015-08-07 2017-08-25 合肥工业大学 一种单测量节点模拟电路故障诊断方法
CN106443297B (zh) * 2016-09-07 2018-05-22 江南大学 光伏二极管箝位型三电平逆变器的决策树svm故障诊断方法
CN106597260B (zh) * 2016-12-29 2020-04-03 合肥工业大学 基于连续小波分析和elm网络的模拟电路故障诊断方法
CN107884706B (zh) * 2017-11-09 2020-04-07 合肥工业大学 基于向量值正则核函数逼近的模拟电路故障诊断方法
CN108535636A (zh) * 2018-05-16 2018-09-14 武汉大学 一种模拟电路基于参数随机分布邻近嵌入胜者为王的故障特征提取方法
US11486925B2 (en) * 2020-05-09 2022-11-01 Hefei University Of Technology Method for diagnosing analog circuit fault based on vector-valued regularized kernel function approximation

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101221213A (zh) * 2008-01-25 2008-07-16 湖南大学 基于粒子群算法的模拟电路故障诊断神经网络方法
CN101329697A (zh) * 2008-06-11 2008-12-24 电子科技大学 一种基于混杂算法的模拟电路状态预测方法
CN102636742A (zh) * 2012-05-15 2012-08-15 长沙河野电气科技有限公司 基于小波神经网络的大规模模拟电路故障诊断方法
CN103824135A (zh) * 2014-03-11 2014-05-28 合肥工业大学 一种模拟电路故障预测方法
CN104198924A (zh) * 2014-09-11 2014-12-10 合肥工业大学 一种新颖的模拟电路早期故障诊断方法
CN105701470A (zh) * 2016-01-13 2016-06-22 合肥工业大学 一种基于最优小波包分解的模拟电路故障特征提取方法
CN107101715A (zh) * 2017-05-25 2017-08-29 贵州电网有限责任公司 一种基于交叉小波的变压器振动信号幅频特征量提取方法

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108629348A (zh) * 2018-05-09 2018-10-09 江南大学 一种基于显著性检验的机械设备故障诊断方法
CN109239585A (zh) * 2018-09-06 2019-01-18 南京理工大学 一种基于改进优选小波包的故障诊断方法
CN109613610A (zh) * 2019-01-17 2019-04-12 中南大学 微震信号到时差的自动拾取方法
CN109782158B (zh) * 2019-02-19 2020-11-06 桂林电子科技大学 一种基于多级分类的模拟电路诊断方法
CN109782158A (zh) * 2019-02-19 2019-05-21 桂林电子科技大学 一种基于多级分类的模拟电路诊断方法
CN110991481A (zh) * 2019-10-31 2020-04-10 河海大学 一种基于交叉小波变换的高压并联电抗器内部松动故障诊断方法
CN110991481B (zh) * 2019-10-31 2022-08-19 河海大学 一种基于交叉小波变换的高压并联电抗器内部松动故障诊断方法
CN111337819B (zh) * 2020-03-02 2021-02-19 武汉大学 基于优选小波基函数的模拟电路故障特征提取方法及系统
CN111337819A (zh) * 2020-03-02 2020-06-26 武汉大学 基于优选小波基函数的模拟电路故障特征提取方法及系统
CN111626329A (zh) * 2020-04-20 2020-09-04 国网上海市电力公司 一种基于lda寻优多尺度纹理特征的绝缘管母故障诊断方法
CN111626329B (zh) * 2020-04-20 2023-07-28 国网上海市电力公司 一种基于lda寻优多尺度纹理特征的绝缘管母故障诊断方法
CN112101116A (zh) * 2020-08-17 2020-12-18 北京无线电计量测试研究所 一种基于深度学习的模拟电路故障诊断方法
CN112819037A (zh) * 2021-01-12 2021-05-18 广东石油化工学院 基于交叉注意力和自注意力的分类参数分布的故障诊断方法
CN112819037B (zh) * 2021-01-12 2024-01-30 广东石油化工学院 基于交叉注意力和自注意力的分类参数分布的故障诊断方法
CN113985261A (zh) * 2021-07-07 2022-01-28 四川大学锦城学院 一种基于神经网络的模拟电路故障诊断方法

Also Published As

Publication number Publication date
WO2019090879A1 (zh) 2019-05-16
CN107894564B (zh) 2020-02-18
US11409933B2 (en) 2022-08-09
US20200394354A1 (en) 2020-12-17

Similar Documents

Publication Publication Date Title
CN107894564A (zh) 一种基于交叉小波特征的模拟电路故障诊断方法
CN105046279B (zh) 一种模拟电路故障模式分类方法
CN107884706B (zh) 基于向量值正则核函数逼近的模拟电路故障诊断方法
CN104793124B (zh) 基于小波变换和ica特征提取的开关电路故障诊断方法
WO2017128455A1 (zh) 一种基于广义多核支持向量机的模拟电路故障诊断方法
CN106483449B (zh) 基于深度学习与复数特征的模拟电路故障诊断方法
CN101819253A (zh) 一种基于概率神经网络的容差电路故障诊断方法
CN102279358B (zh) 一种基于mcskpca的神经网络模拟电路故障诊断方法
CN103995237A (zh) 一种卫星电源系统在线故障诊断方法
CN104198924A (zh) 一种新颖的模拟电路早期故障诊断方法
CN104714171A (zh) 基于小波变换和ica特征提取的开关电路故障分类方法
CN113297922B (zh) 一种高压开关柜故障诊断方法、装置及存储介质
CN104678288A (zh) 基于信息熵和小波变换的开关电流电路故障字典获取方法
CN103267941A (zh) 一种集成开关电流电路故障模式测试方法
CN103064009B (zh) 基于小波分析和有限高斯混合模型em方法的模拟电路故障诊断方法
CN108828437B (zh) 基于云关联系数矩阵的模拟电路故障特征提取方法
CN109782158A (zh) 一种基于多级分类的模拟电路诊断方法
CN105823634A (zh) 基于时频关联向量卷积玻尔兹曼机的轴承损伤识别方法
CN109948267B (zh) 基于圆模型参数的线性模拟电路故障诊断方法
CN101614799A (zh) 应用于电流传感器角差在线监测系统的信号分离择优方法
CN115356599B (zh) 一种多模态城市电网故障诊断方法及系统
CN110160778A (zh) 基于序贯假设检验的齿轮箱故障状态识别方法
CN114355173A (zh) 一种基于多输入残差网络的模拟滤波电路故障诊断方法
Wang et al. A novel approach of feature extraction for analog circuit fault diagnosis based on WPD-LLE-CSA
CN114297963A (zh) 基于K-Means小波神经网络的软电路故障诊断方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant