小波滤波器组在模拟电路故障特征提取中的应用.pdf

第23卷第3期2023年6月Vol.23 No.3Jun.2023湖 南 工 业 职 业 技 术 学 院 学 报JOURNAL OF HUNAN INDUSTRY POLYTECHNIC小波滤波器组在模拟电路故障特征提取中的应用张倩（三门峡职业技术学院智能制造学院，河南 三门峡，472000）摘要 针对模拟电路伪随机信号测试输出响应的特点，提出基于小波滤波器组技术的模拟电路故障特征提取，利用小波滤波器组将测试响应信号分解成多个子频带后分析其时域和功率谱特征，比较正常电路响应各子频带的时域、功率谱特征量，从中选择与无故障特征差异最明显的特征量作为测试与诊断依据。实验结果表明，该故障特征提取方法有效减少故障混叠，提高故障测试分辨率和测试精度。关键词 故障诊断；伪随机激励信号；小波滤波器组；仿真中图分类号TN713 文献标识码A文章编号1671-5004（2023）03-0026-08DOI：10.13787/ki.43-1374/z.2023.03.006Application of Wavelet Filter Banks in Fault Feature Extractionof Analog CircuitsZHANG Qian(Sanmenxia Polytechnic,Sanmenxia 472000,Henan)AbstractAnalog circuit excitation response technique is an important method for analog circuit fault diagnosis.According to the characteristics of pseudo-random signals in analog circuits,the fault feature extraction of analog circuitsbased on wavelet filter banks is proposed,and the test response signals are decomposed into several sub-bands by usingwavelet filter banks,the time domain and power spectrum of the signal in each sub-band are analyzed,and then comparedwith the time domain and power spectrum characteristic of the normal circuit response to each sub-band signal,the mostobvious difference between fault feature and non-fault feature is selected as the basis of test and diagnosis.Experimentalresults show that the fault feature extraction method can effectively reduce fault aliasing,improve fault resolution and testprecision.Key wordsfault diagnosis;pseudo-random excitation;wavelet filter banks;simulation引言目前数字电路故障诊断技术已经达到了较高水平，但由于模拟电路系统自身复杂性，如电参数和元件参数的连续性、元件参数的容差性、电路状态的非线性等，模拟电路故障诊断技术的研究进展一直比较缓慢。1-3模拟电路激励响应技术是模拟电路故障诊断的关键技术环节，其在时域内对模拟电路动态特性进行描述，能够完整表征电路本身的固有特性，可有效区分电路的不同故障状态。4获取电路脉冲响应的方法主要有两种，即直接将单位脉冲信号作为激收稿日期作者简介2022-10-05张倩（1984），女，三门峡职业技术学院智能制造学院副教授，研究方向：电力系统。26第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报励信号，或采用基于伪随机信号的相关辨识方法。5-7基于伪随机信号的相关辨识原理，可获得电路的脉冲响应，进而有效完成模拟电路的故障定位。以伪随机激励信号作为模拟电路测试信号时，响应信号是含有多个频率成分的复杂周期信号，而故障的主要特征可能只集中于某个频带范围内，需要对测试响应信号进行频带分解，实现模拟电路伪随机测试的子带特征提取。子带滤波器组技术将输入的全带信号通过子带分解、子带处理、子带重构等操作，提取模拟电路故障特征。8,9本文针对伪随机信号下模拟电路输出响应的特点，采用小波滤波器组技术在频域内对测试响应信号做频带分解，将不同故障响应之间差异最大的子带信号作为故障特征提取的依据，以有效减少故障混叠，提高故障测试分辨率和测试准确度。1 小波分析的滤波器组实现小波分析作为常用的时频分析方法，可自动调节时频窗，灵活性强，在信号处理中应用广泛。10,111.1 小波分析基本原理信号x(t)的连续小波变换定义为：Wx(a,b)-+x(t)*a,b(t)dt（1）式中，a,b(t)=()1a()t-ba，称为小波基函数，a，b分别为尺度因子、平移因子。对于离散小波变换，可将a按幂级数离散，b在尺度内均匀离散，即a=aj0，b=nb0aj0，其中a0 1，b0 0，j，n为整数。取a0=2，并将 t轴用b0归一化，有：j,n(t)=2-j/2()2-jt-n（2）信号x(t)的离散小波变换为：Wx(j,n)-+x(t)*j,n(t)dt（3）小 波 变 换 窗 口 形 状 为 矩 形b-a,b+a、()0-a,()0+a，窗口中心为()b,0a，时窗和频窗分别为a和a，其中是母小波函数(t)的半时宽，是(t)的半频宽。由于b仅影响窗口在相平面时间轴的位置，而a影响窗口在频率轴上的位置及窗口形状，因此小波变换对不同频率在时域上的取样步长具有调节性，这种特性符合低频信号变化缓慢而高频信号变化迅速的特点。在常规滤波器方面，利用小波的多尺度分解特性能够方便、快捷实现低通滤波、高通滤波、带通滤波等常见滤波。1.2 多分辨分析的滤波器组实现多分辨分析可将信号分解成不同函数空间。设子空间集Vm:m Z为尺度函数(t)的多分辨分析，该子空间集满足时移不变性、单调性、唯一性、伸缩性、Riesz基存在性等性质。基于Octave分解结构的Haar小波滤波器组以小波变换中的多分辨分析为基础，由图1给出的三级Octave分解结构可知，信号总频带多分辨分析仅对低频信号部分进行进一步分解，因此，当需要获取低频段无故障信号与故障信号特征时，可以采用基于Octave分解结构的滤波器组。0()Hz()x n1()H z2200()Hz01()Hz22000()Hz001()Hz220()x n1()x n2()x n3()x n图 1三级 Octave 分解结构基于Octave分解结构的Haar小波滤波器组对信号分解的过程为：首先，寻找一个合适的多分辨分析Vm:m Z，根据尺度函数与小波函数的尺度平移所获得的函数族，对信号进行分解，而 后 由 分 解 系 数 得 到 低 通 滤 波 器 响 应27第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报h(n):n Z和高通滤波器响应g(n):n Z，最后对信号做双通道滤波，将原信号分解为低频和高频部分。设与初始输入信号 x（t）的离散信号序列cm,n:m,n Z相 匹 配 的 尺 度 空 间 为 Vm，把cm,n:m,n Z作为信号在尺度空间Vm中的尺度系数的近似。设h0(n)为低通滤波器，h1(n)为高通滤波器，其输出分别对应原信号的低频和高频部分，称为双通道滤波器组，其输出序列长度为输入序列和滤波器长度之和。对输出序列进行二抽取，对应的cm+1,n和dm+1,n数据长度缩减一半，使得各输出序列长度与输入序列长度保持一致。每次分解，都可将该次输入信号分解成低频和高频部分，由于输出采样率再次缩减一半，总输出序列长度不变，由此产生低频序列和高频序列实现多分辨分解。将上述分解输出的实际信号与空间理论概念联系起来，过渡到空间理论，分解过程示意图如图 2 所示。设 x(t)在总频带(0，）上占据空间Vm，一级分解后，Vm分为两个子空间，频带为0，2的低频带 Vm+1和频带为()2,的高频带Wm+1；二 级 分 解 后，Vm+1又 进 一 步 分 解 为()0,20/4频带的低频Vm+2和()4,2频带的高频Wm+2，如此类推。图 2频率空间（0，）的逐级分解1.3 小波包分析的滤波器组实现小波包分析对尺度空间域小波空间的频带进行同时划分，使得多分辨分析无法细分的高频部分进一步分解。12,13该方法能够根据被分析信号特征，自适应选择对应的频带，与信号频谱匹配，提高时频分辨率。14给出双尺度方程：12()2mt=m Zh()n()2m+1t-n12()2mt=m Zg()n()2m+1t-n（4）式中，h(n):n Z与g(n):n Z分别为多分辨分析中定义的低通和高通共轭滤波器，即g(n)=(-1)nh(1-n)。双尺度方程的递推关系可做进一步推广，记：u0()t=(t)u1()t=(t)（5）递推关系为：u2n()t=2k Zh()k un()2t-ku2n+1()t=2k Zg()k un()2t-k确定的函数集un(t):n N称为由正交尺度(t)生成的正交小波包，N为非负整数集。通常采用金字塔结构小波滤波器组实现小波分析，如图3所示，该结构将频带进行多层次划分，在进一步细分高频部分后自适应选择频带并使之与信号频谱匹配。15图中，第0层为信号占据带宽，二分后得到第1层，依此类推，各层滤波器子带均需二分并传递至下一层。这样每层的滤波器子带都能覆盖信号所占据的带宽。0()Hz()x n1()Hz2200()Hz01()Hz22000()Hz001()Hz220()x n1()x n2()x n3()x n010()Hz011()Hz2210()Hz11()Hz22100()Hz101()Hz224()x n5()x n6()x n7()x n110()Hz111()Hz22图 3三级金字塔分解结构小波包频带多级分解如图4所示。图中，信号x(t)占据频带（0，）的空间为Vm，第一级分解后划分为低频子空间Vm+1和高频子空间Wm+1；第28第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报二级分解后进一步划分为 V2m+1、W2m+1、V2m+2、W2m+2。通过上述方案，对频率（0，）区间进行多次分解，得到相互不包含的多个频域区间。这样原始信号被分解为若干大大小小的包，根据需要分析的信号在各频段的特性，可以适当选取不同大小的包来组装原始信号，使其对信号有一定的适应性。图 4频率空间（0，）的小波包分解2 提取方法及实例分析2.1 基于小波滤波器组的模拟电路故障特征提取方法基于子带滤波的模拟电路伪随机测试故障仿真过程如图5所示。模拟电路在伪随机激励信号下的响应经采样量化为时间序列xn，响应序列经子带滤波器组 H0(z),H1(z),Hn(z)后的输出序列记为 v0n,v1n,vnn，对子带输出序列进行时域与频域的计算，完成故障特征提取。图 5基于小波滤波器组的模拟电路故障仿真过程模拟电路测试与故障诊断时，利用小波滤波器组将测试信号分解成多个子带，针对各个子带进行时域和功率谱分析，而后与无故障对应输出进行对比，从中选择与无故障特征差异最为明显的特征量作为测试与诊断依据，可大大减少故障混叠现象，提升测试分辨率和测试精度。2.2 实例分析仿真电路如图 6所示，选用的是 ITC 97国际标准电路中的状态变量滤波器电路16，各元件参数分别为R1=R2=R3=R4=R5=10 k，R6=3 k，R7=7 k，C1=C2=20 nF。图 6状态变量滤波器电路电 路 仿 真 软 件 选 用 Cadence PSD 15.0，Pspice仿真类型为时域瞬态仿真，温度为27，其他条件均采用默认值。实验中设置PRN信号带宽为fp=5.5 kHz，PRN生成过程中间参数P=3，m 序列的信号长度 N=31，经计算可得 PRN信号的周期为1.67410-3s，生成的PRN信号经过截止频率为 5.5 kHz的低通滤波器后得到的信号如图7（a）所示，经过滤波后的PRN信号的功率谱密度如图7（b）所示。（a）（b）图 7伪随机噪声信号及功率谱密度根据图6所示模拟电路可生成PRGN信号，带宽为fpg=6.4 kHz，信号生成过程中移位叠加次数S=5，生成信号所需的信号源时钟周期为1.110-5s，m序列信号长度取N=127，则生成的29第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报伪高斯噪声信号周期为1.39710-3s。图8（a）为伪高斯噪声信号经截止频率为6.5KHz的低通滤波器后得到的波形，对应的功率谱密度为图 8（b）所示。（a）（b）图 8伪高斯噪声信号及功率谱密度首先选择电路中两个参数型故障，仿真在PRN信号激励下小波-子带滤波方法。故障1：电容C1参数正漂移10%；故障2：电容C1参数正漂移20%。利用电路仿真分析的方法得到电路在正常状态与故障1、故障2状态下的测试响应信号，分别如图9（a）、（b）及（c）所示。（a）（b）（c）图 9PRN 激励情况下被测电路不同状态的输出响应由图9得出，电路在正常状态与故障1状态下的输出响应比较接近，很难看出其中差异，计算未采用滤波器组时几种电路状态测试响应之间的关联度，故障1与正常状态下的输出响应的互相关系数为0.9962，故障2与正常状态响应的互相关系数为0.9922，故障1和故障2之间的互相关系数为 0.9948，在未采用滤波器组进行滤波时，3种电路状态的时域输出响应极为近似，难以区分。利用金字塔结构的Haar小波滤波器组进行故障特征提取，可得到输出响应信号8个不同子带。这里给出了正常状态、故障1状态下从第1-8个子带的输出响应信号，分别如图10（a）-（h）所示，其中左图为正常状态，右图为故障1状态。（a）（b）（c）30第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报（d）（e）（f）（g）（h）图 10不同测试响应各个子带的信号实验结果表明，经小波滤波器组滤波后，故障1和故障2在第1子带时域波形难以区分；而在第4子带，故障1和正常状态输出响应之间的差异较大。为了从数学上说明各个响应信号在不同子带的差异程度，分别计算电路在正常状态、故障1状态和故障2状态下输出响应信号在各个子带的相关系数，考察各个子带之间的相关性，不同子带的相关系数如表1所示。表 1故障 1、故障 2 与正常状态输出信号各子带的相关性（PRN）子带第1子带第2子带第3子带第4子带第5子带第6子带第7子带第8子带故障1&正常0.99670.40580.54370.10160.58130.16910.11890.2946故障2&正常0.99260.36370.44780.01780.47900.06300.28380.0867故障2&故障10.99520.31020.39460.07700.43270.03780.11270.0083分析表1，电路三种不同状态输出的响应信号与各个不同子带的相关系数是不相同的，其中第1子带相关系数较大，说明三种状态的输出响应信号在低频段比较近似，即三种状态响应的轮廓大致相同。第4子带两种故障状态的输出响应信号与正常响应之间的相关系数较小，说明在该子带故障1、故障2两种状态与正常状态之间的差异较大。而故障1和故障2之间差异最大的子带为第8子带，因此对故障1和故障2的诊断宜在第8子带进行。由此可以看出，直接利用测试响应信号区分故障比较困难，对所有测试响应进行子带滤波，即可得到不同电路状态下差异明显的故障特征，说明采用基于子带滤波的方法可以得到差异明显的故障特征。以PRGN信号为测试激励信号，验证基于小波滤波器组技术的故障特征提取方法。故障状态选择故障1和故障2，图11给出了电路在故障1、故障2、正常状态下的测试响应。图中，实线表示故障状态测试响应，虚线表示正常状态测试响应。未采用滤波器组滤波时，首先计算三种电路状态测试响应的相关性，故障1与正常状态输出响应之间的相关系数为0.9978，故障2与正常31第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报状态输出响应之间的相关系数为0.9914，故障1与故障2输出响应之间的相关系数为0.9956。将该参数对比于PRN信号诊断的结果，相关系数相似，说明采用 PRN诊断与采用 PRGN 信号诊断分辨率近似。从图11可以看出，三种电路状态的时域输出响应难以区分。（a）故障 1 与正常状态输出响应的比较（b）故障 2 与正常状态输出响应的比较图 11故障 1、故障 2 与正常输出响应的比较对三种状态下的响应进行子带滤波，分别计算电路在正常状态、故障1状态和故障2状态下输出响应信号各个子带间的相关系数，考察各个子带之间的相关性，不同子带的相关系数如表2所示。表 2故障 1、故障 2 与正常状态输出信号各子带的相关性（PRGN）子带第1子带第2子带第3子带第4子带第5子带第6子带第7子带第8子带故障1&正常0.99780.47760.53050.04860.56200.06010.06270.2077故障2&正常0.99140.32640.39360.02650.40270.15860.04230.1413故障2&故障10.99780.66510.70560.19690.73310.20300.11400.0924分析表2，电路三种状态响应信号在第1子带相关系数较大，说明三种状态输出的响应信号在低频段比较近似。第4子带两种故障状态输出的响应信号与正常响应之间的相关系数较小，说明在该子带故障1、故障2两种状态与正常状态之间的差异较大。而故障1和故障2之间差异最大的子带为第8子带，因此对故障1和故障2的诊断宜在第8子带进行。分析上述实验，测试响应经过滤波器后，故障1、故障2输出响应信号特征向量相对差异变大，说明采用小波滤波器组技术对测试响应进行滤波能够提取差异明显的特征向量。在实际进行诊断时，测前对电路响应的各个子带进行仿真，得到故障分辨率最高的子带（即相关性最小的子带），将不同故障响应之间差异较大的子带作为故障特征提取的依据，可有效提高电路诊断的分辨率，使诊断结果具有较高的可靠性。3 结论针对模拟电路伪随机测试中输出响应的特点，提出了一种基于信号子带分解技术的模拟电路故障特征提取方法，采用小波滤波器组技术将信号的频谱均匀或非均匀地分解成若干部分，对不同电路状态响应的各个子带进行比较，从中找出不同电路状态响应之间差异较大的子带，并将其作为故障定位的依据。将基于小波滤波器组的信号子带分解技术应用于模拟电路故障诊断，可有效减少故障混叠，提高故障测试分辨率和测试准确度。参考文献1李光升,欧博,石海滨.现代模拟电路故障诊断技术发展综述J.计算机测量与控制,2019,27(5):5-8.2郭珂,伞冶,朱奕.现代模拟电路智能故障诊断方法研究与发展J.电子设计工程,2012,20(2):177-180.3阳景,潘强,潘红兵.基于时间卷积网络的模拟电路故障诊断方法J.电子测量技术,2019,42(5):128-132.4张朝龙,何怡刚,杜博伦,等.一种基于多核相关向量机的模拟电路故障预测方法J.电子测量与仪器学报,2019,33(11):96-101.5薛国强,底青云,王若,等.多通道瞬变电磁法资料处理方法技术综述J.地球物理学进展,2020,35(1):211-215.6刘恒,孙晋,邰凡彬,等.基于相关分析法的脉冲响应辨识实验设计J.实验技术与管理,2018,35(5):53-59.32第3期2023年湖 南 工 业 职 业 技 术 学 院 学 报7王学伟,王艳君.m序列调制的正弦离散伪随机动态测试信号的完备性分析J.中国电机工程学报,2018,38(12):3529-3537,12.8陈勇,王玉文,洪文军,等.面向5G三大应用场景的F-OFDM系统PAPR抑制算法研究J.通信技术,2020,53(4):816-823.9赵艳,万坚.一种单载波宽带信号非线性均衡技术J.计算机工程与应用,2012,48(5):135-137,142.10童耀南,曹鹂晨,赵舜楠.天牛须搜索算法实现小波滤波器直接频域逼近J.电子测量与仪器学报,2021,35(3):56-63.11童耀南,何怡刚,李宏民,等.频域逼近优化的开关电流电路小波变换方法J.电子测量与仪器学报,2014,28(3):262-271.12李默.压电耦合智能板结构中微损伤识别方法研究J.电子科技,2017,30(4):5-10.13宋固全,涂飞,张纯.基于小波包降噪与有限元模型修正的结构损伤识别方法J.力学季刊,2010,31(4):528-533.14颉潭成,马君达,袁子皓,等.基于多元传感器的数控刀具磨损状态的实验研究J.组合机床与自动化加工技术,2016(9):28-31.15杨琳,王从庆,缪鹏,等.基于Huffman最优二叉树支持向量机的舱音记录器背景信号识别J.宇航学报,2011,32(6):1428-1434.16杨晓朋,陈伟,王鹏展,等.基于DCCA-IWO-MKSVM的模拟电路故障诊断方法J.计算机应用与软件,2020,37(1):271-276.责任编辑/校对：周哲民压力曲线自动跟踪精准斜率控制与恒压控制图 6实验曲线5 结论本文根据压力变送器的发展现状以及其检测需求，设计了基于NI cRIO-9073 PAC（作为下位机）与LabVIEW的高精度动态测试系统，针对系统的硬件设计和软件设计进行了详述。该动态测试系统实现了压力变送器特性参数高精度全自动测试，以及数据采集和对比分析，自动生成测试报告，有效提高了工作效率。参考文献1罗扬锋.压力变送器检定、使用中的问题及其相应解决J.化工管理,2014(17):187-192.2武文凯,李明辉,巩强令.基于S71200PLC与LabVIEW的泵轮检测设备控制系统J.仪表技术与传感器,2020(12):82-85.3Matias J,Garay J,Toledo N,et al.Toward an SDN-EnabledNFVArchitectureJ.IEEECommunicationsMagazine:2015,4(4):187-193.4王华忠,刘涛,姚俊,等.基于云平台的远程控制系统虚拟实验J.建模与仿真,2019,8(3):95-101.责任编辑/校对：周哲民（上接第25页）33