基于红外热像采集的BPSO-WD-ESN电路板故障预测

贾彦荣,崔展博,吴倩倩,侯昌宇,李兴丽,高云红

(1.石家庄海山实业发展总公司， 石家庄 050200；2.石家庄海山航空电子科技有限公司 电子室，石家庄 050200；3.沈阳航空航天大学 自动化学院，沈阳 110136)



基于红外热像采集的BPSO-WD-ESN电路板故障预测

贾彦荣1,2,崔展博1,2,吴倩倩1,2,侯昌宇1,李兴丽1,2,高云红3

(1.石家庄海山实业发展总公司， 石家庄 050200；2.石家庄海山航空电子科技有限公司 电子室，石家庄 050200；3.沈阳航空航天大学 自动化学院，沈阳 110136)

针对电路板集成度高、故障规律波动无序、伪周期性等难以准确预测的问题，借助红外热成像仪应用于电路板故障诊断中，提出了一种基于Bootstrap采样技术的粒子群优化(BPSO)小波回声状态网络(WD-ESN)模型，运用红外热像故障分析模式，对功能电路板进行动态分析，获取相关模块的热成像变化数据。利用WD解析出表征不同元器件变化区间，根据不同热量变化规律，建立相应的储备池网络预测模型，采用Bootstrap粒子群算法对WD-ESN网络的组合权值矩阵进行优化。分别采用三种不同算法和经过BPSO优化后的WD-ESN模型对某型航电产品功能电路温度图谱进行分析，仿真结果表明，经过PSO优化后的WD-ESN模型能够更加准确预测温度变化趋势(误差低于15%)，满足了高速、实时、准确的红外温度预测的要求。

红外热像技术;小波分解解析；WD-ESN网络；SCADA系统；整合重构

随着对功能电路板集成性和低能耗的要求越来越高，快速准确地对电路板的各种故障进行分析和定位[1]，已被作为衡量航电产品可靠性技术的重要特性。传统的电路板检测技术是一种接触式检测方法[2]，首先需要测量电路中功能模块电流、电压等元器件必要的参数，然后由检测人员结合电路原理和故障现象进行分析[3]，才能诊断出电路板上的故障部位或元器件[4]，不仅整个检测过程耗费时间较长，而且要求检测人员具有较强的专业技能[5]。 航空电子产品在通电情况下，各类电子器件(电阻、电容、光耦、集成电路芯片等)均会发热[6]，其工作过程会伴随功耗变化，进而影响电路板的性能变化。有效分析功能电路板的温度变化规律，对故障模块预测显得尤为重要。

GontarZ于2001年提出了改进型径向基神经网络(RBF)预测模型[7]，可以有效反映输入输出间的非线性关系，在学习过程中动态确定隐层神经元个数，但很难保证所建立网络精确预测。SVM采用最小化系统泛化误差可有效获得全局最优解[8]，但是求解过程多是随机确定模型参数，将问题转化为求解高维空间下的规划问题且不具备动态记忆能力，计算量较大。针对这一问题，李明等学者利用动态储备池理论建立了减少预测训练权值计算量的回声状态网络(ESN)[9]，通过对红外热像仪采集的功能电路温度变化的仿真分析表明，利用常规的ESN算法对功能电路的温度变化进行预测时，训练过程和学习速率方面得到一定改善，能够减少误差，但因储备池结构固定，仅能解决特定温度变化范围的预测，对波动较大的复杂非线性时间序列处理过程较差。

针对以上问题，本文将红外热像仪技术结合小波分解应用到时序波动，无序的电路板故障预测中，提出了基于Bootstrap采样技术的粒子群优化小波回声状态网络(WD-ESN)模型，对网络权值进行优化，改善了预测效果，减少了权值训练的计算过程。红外热像仪即是基于这一热量变化过程来记录电路板各观测点的温度变化，并将其记录数据与数据库中存储的正常电路状态进行对比分析，预测出热量变化“失常单元”，以此为电子产品设计人员及电路维修人员提供排故依据。

1　红外热像系统构建

1.1红外热像系统功能

本应用系统借助红外热成像仪采集某航电产品正常测试数据和故障测试数据等多方面信息，采用数据处理、知识挖掘、专家经验等多方面手段，通过融合数据挖掘技术、优化模型等多种信息，给出电路故障分析方案。系统功能模块组成如图1所示。

图1　系统功能模块组成

系统遵照模块化方式进行设计，核心功能分为原始电路红外采集数据处理与建库、模型优化及故障定位三个部分。具体流程如下：鉴于现场实时数据和历史数据形成的EMS原始数据噪声高、动态特性强等[10]问题，首先将原始数据通过除噪、滤波、缺失数据填补等处理环节，形成模型使用的应用数据库。另一方面，实时数据、历史数据和计划数据将一并送入实时预测模块。在此过程中，系统将依据SCADA系统实时监控数据来完成对电路板故障分析环节的在线校正，将误差减到最小，完成在线实时预测。红外热像技术应用于电路板故障预测的功能流程图如图2所示。

系统输入来源为EMS/MES系统数据库的计划数据、实时数据和历史数据，还包括现场定修等专家预案和人工干预等。此类信息将为系统四大功能模块提供服务。

工程实践中，技术人员有时更加关注诸多不确定性因素影响下的故障变化趋势，并能有效衡量各类不确定性因素对预测结果的影响。可通过实施区间定位给用户展示出相应阶段电路温度曲线的最高/最低波动区间，一般设定其区间置信率为80%。

图2　系统功能流程图

1.2数据库基本设计

本红外热像系统关系数据库管理系统采用DB2，表的设计大致分为四部分，即主要处理计划数据的计划类表，主要处理历史数据的数据表，主要处理电耗类数据的数据表以及负责帐号管理等辅助功能的数据表。基本表设计包括一些主外键连接以方便联合查询，但也未设置过多此类连接，以免造成数据库调试修改困难。

数据库与系统的每个模块都具备输入输出接口[11]，方便各模块依功能需要从数据库取得相应数据，或存储一些关键参数、操作信息等。

2　小波分解改进ESN网络模型(WD-ESN)

2.1ESN模型算法

回声状态网络预测(echostatenetwork,ESN)算法是一种动态递归网络算法[12]，其基本思路是通过内部的动态储备池DR(DynamicalReservoir，DR)产生一个随输入信号u(t)变化生成的动态非线性系统。该系统时序预测方面表现出较好的性能。

ESN预测模型建模过程如下：

Step1 设定input(输入单元)、output(输出单元)，可得预测系统如公式(1)、(2)所示：

UcDi(n)=[cDi(n),cDi(n-V1),…,cDi(n-VM-1)]

(1)

YcDi(n)=[cDi(n+VM)]

(2)

式中，V1,V2,…VM-1,VM分别表示温度变化延迟时间，M表征输入维数。

Step2 储备池更新 确定储备池(DR)大小NcDi，缩放矩阵使谱半径ρcDi<1，生成连接矩阵，依据生成的输入连接和输出反馈权值，结合ESN建模理论，得到式(3)的DR状态更新推导公式：

(3)

Step3 输出权值计算

令xcDi(0)=0，当实时功能电路温度变化数据动态导入ESN预测模型，DR(储备池)随(4)式不断更新，输出权值矩阵如式(4)所示：

(4)

(4)式中 DcDi=[ucDi(n)xcDi(n)]，TcDi=[ycDi(n)]

Step4 温度变化预测

温度变化预测方程如式(5)所示：

(5)

2.2小波分解(wavelet decomposition,WD)基本算法

小波分解是建立在基函数基础上的时间～尺度(时间～频率)分析方法，适合将复杂非线性信号分解成若干不同频带的细节和基本信号序列[13]，针对不同的尺度信号分别建模，然后再还原信号得到更为精确的预测信号，分解过程如式(6)所示：

(6)

(7)

分解步骤如下：

步骤1：设置小波基函数，将其与待分解温度变化信号初始时刻对齐；

步骤2：选择初始小波变换尺度a0，计算初始时刻待分解的温度变化信号与小波基函数的相似程度C。C值的大小反应出温度变化信号与所选小波函数逼近效果；

步骤3：以温度变化信号为参照，令小波函数沿着时间轴平行移动b0，利用步骤1～步骤2，计算平移后的小波变换系数C1，重复这一过程直到完整覆盖温度变化信号；

步骤4：重新设置小波变换尺度a1，然后重复步骤1～步骤3；

步骤5：针对待分解信号重复步骤1～步骤4(尺度伸缩)，最后，将得到各电路板功能模块温度变化在这些小波函数的投影效果。

2.3WD-ESN预测模型

由于电路板各功能模块温度变化波动无序、幅度大、易受外界干扰，这就增加了温度变化预测的难度，一旦发生波动，或产品性能下降，都将严重影响其算法模型的应用推广。

基于以上分析，为了解决电路板时序逻辑芯片温度变化信号非平稳和多尺度变化的问题，本文利用小波分解算法改进ESN网络泛化效果。通过WD将温度变化信号进行分频解析，经平移及伸缩后分解为1个低频基本分量和多个高频细节分量，希望利用分频预测来提高ESN网络的预测效果，最终将分解预测效果重构，得到电路板功能模块温度变化整体预测效果。

WD-ESN预测建模步骤如下：

(1)确定时序信号的分解层次m，初始化伸缩变量a0，通过式(7)，依次平移b0个单位。可将温度变化初始信号分解为基本分量cA1，高频分量cD1；

(2)改变伸缩变量，对细节分量cA1进行小波变换，可将其分解为基本分量cA2和细节变量cD2；

(3)根据分解层次m，可依次将原始温度变化序列分解为cAm，…，cDm，cD2,cD1；

通过小波重构整合理论，可整合出电路板温度变化终值，如式(8)所示：

(8)

为便于准确预测温度变化趋势，以真实值y实(n)和预测值y(n)的偏差最小值为目标，借助线性回归来确定各分支系数，如式(9)所示：

O=(JTJ)-1JTT

(9)

其中，O=[w1w2w3…wm+1] ，

T=[y(1)y(2)…y(λ)],

3　基于BPSO的WD-ESN预测模型优化

本文温度变化预测用户选择某航空产品的5个主要功能模块，各预测环节的运行状态不同，其用电特征和温度变化波动情况各异，故应针对各模块分别进行预测，模型参数的选取就显得尤为重要。对多参数、过拟合、非线性、多尺度、多约束等耦合系统寻优问题，BPSO优化效果较为理想。

文化是民族的血脉，是人民的精神家园，人类的发展史亦可称为一部文化的继承与发展史。而高校是人才培养的摇篮，是滋生文化的沃土，是传播文化的重要领地，理应是文化建设中的重大抓手。在日新月异、瞬息万变的时代里，高校文化建设更是应该把握时代发展脉络，不断调整发展方向、丰富发展内涵，积极推进文化大发展与大繁荣。

3.1基于BPSO的WD-ESN参数优化

PSO算法用于WD-ESN网络组合权值优化。一方面采用数值迭代的求解方式，不需要进行矩阵求逆操作，可以避免因病态矩阵而导致权值偏移的现象；另一方面编码方式采取实数编码，即用粒子编码来表征待求的网络权值，将网络权值寻优过程转化为粒子寻优问题。为提高模型实用性，本文在建立WD-ESN算法模型基础上，提出基于Bootstrap采样技术的粒子群优WD-ESN模型。

在利用BPSO算法进行组合权值优化时，将粒子进行归一化操作，并在“评价标准(适应度函数)”中引入罚函数项，以弱化噪声和共线性的干扰。因此，对于给定好的m个回声网络，cD1,cD2,cD3,…,cDm,其优化组合权值w(w1,w2,…,wm)步骤如下：

Step2：将初始化粒子群按照式(10)进行处理，种群粒子代入问题“评价标准”之中。初始20组粒子，通过MATLAB2012中的HeavySine函数进行小波分解层数的自适应分析，对于本文预测对象，适应度函数选取如下：

φ(x)=1/L

(10)

给定输入参数作为寻优变量，WD-ESN的动态优化目标函数选取如式(11)所示：

(11)

(12)

FP表征罚函数项，用以约束权值和为1，如式(13)所示：

(13)

Step3：当t

(14)

(15)

为限制PSO粒子的盲目搜素，将位置和速度分别限制在[-Xmax,+Xmax]和[-Vmax,+Vmax]区间，k是当前迭代次数，Vid粒子的速度是非负加速度学习因子,σ是惯性权重。为了改善算法的收敛性能通常取值如式(16)所示：

(16)

其中σmax、σmin分别为最大、最小加权系数，t为当前迭代次数，tmax为算法迭代总次数。

每执行一次，需对粒子约束条件进行验证，如超过速度或位置约束，则排除该粒子，并对其他粒子进行更新操作，直到达到最大迭代次数为止；

(17)

(18)

Step6：误差对比分析，将得到最优位置或者达到迭代进化次数上限(本文中取80) 的寻优粒子代入WD-ESN网络，如果达到预期温度变化预测精度需求，则寻优停止，否则转至Step2。

3.2模型预测流程

Step1 针对输入的电路板温度变化数据，借助WD分解得到cAm，…,cDm,cD2,cD1；

Step2 结合红外成像SCADA系统，对WD分解后不同尺度的电路板温度变化分量进行单支重构，得到对应尺度的电路板温度变化序列，分析其对应尺度变化情况，寻找其规律，为后续ESN建模提供指导；

Step3 依据不同尺度温度变化序列的特征，借助BPSO优化选择合适的ESN网络权值参数，分别预测对应频率温度变化序列ycDi(n)。

其BPSO算法优化权重过程如下：

(1)输入待处理的组合权值矩阵数据改变目标优化函数L；

(2)约定最大寻优次数tmax=50，初始值t=1；

(5)t=t+1 ；

(7)对20个过程(6)所示最优组合求均值如式(19)所示：

(19)

(8) 相应的ESN网络分支输出如式(20)所示：

(20)

Step4 通过线性回归模型整合各分支预测结果(温度变化预测流程如图3所示)。

图3　温度变化预测流程

3.3预测仿真分析

为验证本文所建立红外热像技术应用于电路板温度变化预测模型的实用性，分别采用经典SVM，径向基预测模型及ESN网络进行对比预测分析，预测对象为CPU模块、电源模块、驱动电路模块、通信模块，选择2015年5月份某航电产品通电测试段内的1 920个温度变化数据(每隔15秒的电路板温度变化值为一个采样时刻数据)作为训练样本，预测电路板温度变化趋势，仿真效果拟选用三种评价准则进行综合分析。

(1)平均绝对百分比误差如式(21)所示：

(21)

(2)归一化均方根误差如式(22)所示：

nrmse=

(22)

(3)均方误差如式(23)所示：

(23)

式中，yreal(i)表征温度变化验证真实数据，ypre(i)为预测点，q表征温度变化预测点数，不同模型预测效果如表1所示。

由预测对比结果可知，本文提出的基于红外热像技术组合模型对某航电产品电路板系统主要模块进行温度变化预测时，不同的评价准则的预测效果均优于其他几种模型。ESN虽然输出波动不大，但预测精度不高,原因可能是由于其网络结构固定，性能易受储备池神经元不同选择的影响。径向基网络因为其学习过程易形成局部极值点，因此导致模型预测效果不够理想。SVM模型输出波动最大，性能易受核函数不同选择的影响,预测精度最低，但运算速度上具有明显优势(0.309 11秒)，且适合于小样本处理(数据点数小于4 000)。而本文提出的BPSO优化WD-ESN预测模型，无论输出精度(误差均低于10%)，还是预测效果的稳定性都远优于其他几种算法。

表1　各主要电路板用户温度变化预测误差率对比

4　结束语

(1)电路板标准信息库(SCADA系统)的建立是电路板检测的基础工作。要完善故障的相关信息，如温度特征信息、温度特征曲线等。可根据相关的温度特征信息将故障分级，以方便对电路板的维修。本文以红外热像技术应用于电路板温度变化预测为背景，提出基于BPSO和WD-ESN组合模型，将电路板功能模块波动无序和非线性的温度变化预测问题转化为不同尺度的温度变化分量，针对不同的尺度信号分别建立ESN网络，更有针对性；

(2)基于Bootstrap采样技术的粒子群优化算法有效解决了多参数、过拟合、非线性、多尺度、多约束等耦合系统寻优问题；

(3)针对不同温度变化预测对象的仿真实例表明，经BPSO优化后的WD-ESN模型的预测精度更高，适应性更好，可以有效处理电路板温度变化难以准确预测的问题，为解决时序电路故障定位提供了一个良好的思路。

[1]张聪.PCB外观检测仪中基于FPGA的CMOS图像采集[D].成都：电子科技大学,2013.

[2]贾延玺.加载电路板检测系统研制[D].哈尔滨：哈尔滨工业大学,2009.

[3]崔昊杨,许永鹏,杨俊杰,等.基于MIV和BRBP神经网络的电路板红外诊断方法[J].激光与红外,2014,44(4):382-386.

[4]FARROKHRAHIMI,ALIIPAKCHI.Overviewofdemandresponseunderthesmartgridandmarketparadigms[C].OpenAccessTechnol.Int.Inc.(OATI),Minneapolls,MN,USAinGaithersburg,MD,2010:1-7.

[5]龚镇，李运祯.基于红外温度变化规律诊断电路故障[J].电子技术，2009，8(25)：37-39.

[6]周孟特.基于ARM11的便携式红外电路故障检测硬件系统[D].成都：电子科技大学,2013.

[7]田裕鹏.红外检测与诊断技术[M].北京：化学工业出版社，2006.

[8]崔伟.电路板故障红外热像检测关键技术研究[D].南京：南京航空航天大学,2011.

[9]MENG LI,WEI DONG.,Model of multi-sensor data fusion and trajectory prediction based on echo state network[C].IEEE Industrial Electronics Society Beijing(Shenzhen) Chapter、Changchun University of Technology,China、Intelligent Information Technology Application Research Association(IITA Association),Hong Kong,2010:4.

[10]周孟特,叶玉堂,王智芳,等.便携式红外电路故障检测系统[J].光电工程,2011,38(8):101-106.

[11]曾佩佩.基于主动轮廓模型的PCB红外图像分割[D].天津：中国民航大学,2015.

[12]景华.基于ARM11的嵌入式碳氢分析仪的研制[D].南京：南京信息工程大学,2015.

[13]李加超.基于NiosⅡ的多功能图像信号源的设计与实现[D].太原：中北大学,2015.

[14]刘宁.便携式红外电路故障检测仪底层软件设计[D].成都：电子科技大学,2013.

[15]雷冬梅.非侵入式电路故障诊断及家用电器用电状况在线监测算法研究[D].重庆：重庆大学,2013.

[16]刘宁,周孟特,王智芳,等.便携式红外电路故障检测仪视频采集驱动设计[J].计算机工程与科学,2011,33(12):189-193.

(责任编辑：刘划英文审校：赵欢)

FaultpredictionofBPSO-WD-ESNcircuitboardbasedonthermalinfraredimaging

JIAYan-rong1,2,CUIZhan-bo1,2,WUQian-qian1,2,HOUChang-yu1,LIXing-li1,2,GAOYun-hong3

(1.HaishanCorporationofIndustrialDevelopmentinShijiazhuang,Shijiazhuang050200,China;2.TheElectronicChamber，ShijiazhuangHaishanAvionicsCo.,Ltd.,Shijiazhuang050200,China;3.CollegeofAutomation，ShenyangAerospaceUniversity,Shenyang110136,China)

Somefaultsofcircuitboardaredifficulttopredictaccuratelyduetohighintegrationoftheboard,disorderlyfluctuationsoftimingmodulefailure,andpseudo-periodicity.Awaveletechostatenetwork(WD-ESN)modebasedonbootstrapsamplingtechniqueandparticleswarmoptimization(BPSO)wasproposedusingthermalinfraredimagingtodynamicallyanalyzethefunctionoftheboardandobtainthechangesinthermalimageofrelativemodule.ThevariationrangeofdifferentcomponentscharacterizedbyWDanalysis,correspondingreservepoolnetworkpredictionmodelwasestablishedbasedontheheatvariation,andcombinationweightmatrixofWD-ESNnetworkwasoptimizedusingBootstrapPSOalgorithms.ThreedifferentalgorithmsandWD-ESNmodeloptimizedbyBPSOwereusedrespectivelytoanalyzethefeaturescircuitofacertainavionicsproduct.SimulativeresultsshowthatWD-ESNmodeloptimizedbyPSOcanmoreaccuratelypredictthechangeoftemperature(errorLowto15%)andmeetthepredictionrequirementssuchashigh-speed,real-time,andaccurate.

hermalinfraredimagingtechnology;waveletanalysis;WD-ESNnetwork;SCADAsystem;integrationandreconstruction

2015-11-08

国家自然科学基金(项目编号：51277126);航空科学基金(项目编号：2014ZB54010)

贾彦荣(1982-)，男，工程师，主要研究方向：航空电子产品可靠性维修，E-mail:2363022568@qq.com。

2095-1248(2016)04-0066-07

N945.12

A

10.3969/j.issn.2095-1248.2016.04.012