基于系统搜索的容差模拟电路多参数故障诊断

董海迪 刘刚 何兵 李永强 王世涛

摘要：针对容差条件下非线性模拟电路的多参数故障诊断问题，提出一种系统搜索算法求解故障参数，实现故障元件的定位和参数辨识。首先将非线性测试方程的多解求解转化为差分方程的初值问题;然后采用牛顿法迭代生成简单曲线簇，沿曲面方向进行测试方程的解搜索，并通过四阶龙格库塔法计算初值问题的数值解，显著提高方程组解的求解效率，降低漏解的风险;最后将故障集验证转化为容差约束下的线性规划问题求解，实现故障元件定位，并通过实际电路实验验证该算法的有效性。该算法适用于含BJT和CMOS等非线性电路的故障诊断，具有较高的诊断准确度和参数辨识准确度。

关键词：模拟电路;参数故障;系统搜索;非线性

中图分类号：TN06 文献标志码：A 文章编号：1674-5124（2019）07-0134-06

收稿日期：2018-01-18;收到修改稿日期：2018-05-20

基金项目：国家自然科学基金资助项目（61403399）

作者简介：董海迪（1988-），男，湖北武汉市人，博士研究生，研究方向为模拟电路测试与故障诊断。

0 引言

模拟电路故障诊断是电子设备有效性验证的重要问题[1-2]，模拟电路故障诊断方法依据仿真进行的时间划分为测前仿真方法和测后仿真方法[3-4]。参数故障诊断通常采用测后仿真方法，通过求解测试方程得到待测元件的参数值，与标称值比对进行故障隔离定位。参数辨识法分为伴随网络法和多频诊断法[5-6]。伴随网络法是根据Tellegen定理计算得到测点处各元件的灵敏度值，采用牛顿迭代方法和区间法等求解被测电路构成的多元非线性方程组，该方法计算量大。多频诊断法是指通过增加测试频率，得到更多的故障方程数，据此来确定元件值，该方法收敛效果差。近几年来参数辨识法取得了新的进展，He等[7]将L1最优范数引入到非线性电路的参数估计，采用神经网络方法求解;Guo和Savir[8]采用自回归模型和系统辨识工具进行模拟电路故障参数辨识;LONG等[9]采用改进的灵敏度计算方法进行模拟电路的故障参数辨识。

对于复杂线性电路或非线性电路，电路代数方程很难求取或不存在、求解效率低等问题仍然存在，困扰着参数辨识方法的进一步发展。此外，由于元件容差影响，存在故障掩盖现象。同时，考虑方程的非线性，测试方程可能存在多解情况。综上分析，有必要寻求新技术途径来实现对模拟电路的故障参数辨识，本文针对容差条件下模拟电路的多软故障诊断问题开展研究，提出一种基于系统搜索的故障诊断和参数辨识方法。

1 算法原理

假设某被测电路中可能存在n个参数故障，给该被测电路施加激励信号，同时在r个测点处测量输出电压。对于m个激励电压信号集合，测量得到mr个输出电压，其中mr≥n。从中选取n个测量电压，作为输出电压向量u=[U₁，…，u_n]^T的元素，每個测量电压u_i是电路可能故障参数x₁，…，x_n的函数，表示为u_i=f_i（x），i=1，…，n，其中x=[x₁，…，x_n]^T。得到的测试方程为

f（x）=0（1）式中：f（x）=f（x）-u;

f（x）=[f₁（x），…，f_n（x）]^T。

实际电路中，往往无法得到非线性函数f（x）的准确数学表达式，同时由于可用测点和激励点有限，测试方程中包含的信息往往不足以找到实际的故障参数集合，在数学上意味着式（1）包含多个解，因此需要找到一种能得到测试方程多解的方法。参照文献[10]，将代数问题转化为差分方程求解是一种有效寻找多解的方法。Chao等[11]将该思想拓展为寻找非线性电阻电路多直流工作点的系统搜索方法。本文进一步将系统搜索的概念拓展应用于模拟电路故障诊断过程中测试方程的多解求解。

为了找到式（1）的多个解，将该问题转化为初值问题。设计差分方程为

由式（2）得到：

（f_i（x（t））=f_i（x（0））e^-t=0.i=1，…，n-1

f_n（x（t））=f_n0e^±t（3）

初值条件x（0）的选择满足如下条件：f_i（x（0））=0，i=1，…，n-1，即x（0）位于空间曲线L上，该空间曲线为曲面f_i（x）=0，i=1，…，n-1的交叉部分。解式（3）中前n-1个方程得到：对于任意x（0）∈l，跟踪点x（t）保留在曲线l上。初值问题式（3）的解x（t）能够通过数值积分方法得到。任何位于曲线l上的点x（t_k），满足f_n（x（t_k））=0都是式（1）的解。

由于df/dt=（df/dx）（dx/dt），式（2）在x空间具有如下形式：

dx/dt=J^-1（x）w（x）（4）式中：J（x）=df/dx是f（x）的雅克比矩阵;

w（x）=[-f₁（x），…，-f_n-1（x），±f_n（x）]^T。

在文献[11]中，已经证明了f_n（x）在沿曲线l方向的导数满足：式中：△_ni（x）（i=1，…，n）由J（x）中最后一行元素组成。

由式（3）可知，式（2）中f_n（x）取负号时，f_n（x）沿曲线l吸引到0;式（2）中f_n（x）取正号时，f_n（x）沿曲线l远离0。文献[11]已证明在搜索过程中，式（2）中f_n（x）的符号改变发生在雅克比矩阵J（x）的行列式改变符号时或找到了方程的解时。

2 算法求解

兼顾截断精度和计算量，采用四阶显式龙格库塔法[12]求解初值问题式（4）的数值解。因为隐式积分方法要求对电路进行二阶灵敏度分析，其无法对式（4）进行求解。考虑到测试式（1）的代数方程无法准确得到，测点电压f_i（x）（i=1，…，n）以及f_i（x）对参数x_j（j=1，…，n）的偏导数需要通过对给定参数X₁，…，x_n的电路进行仿真，通过直流和灵敏度分析得到，即得到了w（x）和J（x）。

为了找到测试方程的所有可行解，本文提出采用系统搜索方法，关键是找到简单的空间曲线l。否则，如果曲线l是多分枝的，可能会失去部分解。因此，本文提出了一种产生多个简单曲线的方法，生成的每条曲线都可能找到多个解。考虑含n-1个方程的系统：

f₁（x₁，x₂，…，x_n-1，x_n）=0

f₂（x₁，x₂，…，x_n-1，x_n）=0

…

f_n-1（x₁，x₂，…，x_n-1，x_n）=0

选择某参数取标称值，其余参数为变量，对该方程组进行n次求解。例如，首先设置，针对参数x₁，x₂，…x_n-1采用牛顿法解式（6），得到解x₁^*，x₂^*，…，x_n-1^*与x_n^nom组成初值向量x（0），该向量在空间曲线l（标记为l^（1），见式（2））上定义了起始点，使用第1节描述的搜索方法，产生空间曲线l^（1），得到式（1）的一个或多个解;为找到其他可能的解，设置x_n-1=，针对参数x₁，x₂，…，x_n-2，x_n，解式（6），得到新的向量，采用上述過程产生空间曲线l（2）。最后是设置x₁=x₁^nom，考虑参数x₂，…，x_n。上述n次求解方程组的过程中，如果遇到牛顿法发散，则舍弃该结果。

为简化说明，式（6）的系统用以下形式描述：p（x）= 0，其中p（x）=[f₁（x），…，f_n-1（x）^T。此外，用x_-k表示向量[x₁，…，x_k^nom，x_k+1，…，x_n]^T，由解x_-k^*得到的空间曲线表示为l^（n-k+1）。沿着曲线l（n-k+1），k∈{n，n-1，…，1}，采用系统搜索方法能够找到方程f（x）=0的解。该方法从空间曲线的初始点开始，分别从两个方向进行跟踪。首先，将负号赋予f_n（x），求解式（4），跟踪曲线，该曲线标记为N部分;之后，将正号赋予f_n（x），选择同样的起始点，求解式（4），跟踪曲线，该曲线标记为P部分。当获得一个解时，搜索过程暂时锁定，检验是否满足物理限制条件。如果解不被接受，过程解锁并继续。如果计算步长m达到最大值M，该过程从起始点开始搜索曲线的P部分。当总步长m=M=2M时，方法锁定。因此，模块的输出量是搜索次数m和方程解（m

实际应用中测量精度有限，且考虑元件存在容差，因此对上述方法进行修正。对于曲线上的任何点x，检测是否满足f_n（x）∈[u_n△u_n，u_n+△u_n]，其中△u_n依赖于测量精度，而不是检测是否满足f_n（x）=f_n（x）-u_n=0。基于大量数值仿真实验，△u_n设置为5ε，其中ε为测量精度。通常，曲线上的若干点满足这一条件，而且所有点都必须视为测试方程的合适解。因此，不同于理想情况，得到参数范围集合而不是一系列具体值。计算中值元=（l/2）（万+对），i=1'...，。，考虑x=巨1，"..xn1T为方程的合适解。

3 故障辨识

诊断的目的是找到电路中满足测试方程的可能发生故障的参数x₁，…x_n的实际值，它们是非线性代数式（1）的解。由于方程可能存在多解，因此需考虑如何从中找到真实解。采用上述广义系统搜索方法寻找到方程解后，利用文献[13]中方法进行故障验证。

考虑电路中所有元件参数在容差范围内扰动，由文献[13]得到：式中：y_i——第i个激励信号下测量值与标称参数下测量值的偏差;

x_j——第j个元件的参数偏差;

a_ij——第i个激励信号下电压e_i对第j个元件参数的电压灵敏度;

m——激励信号数量。

容差模拟电路正常状态下的元件参数偏差满足-ε_j^-≤x_j≤ε_j⁺（j=1，…，n），其中ε_j⁺和-ε_j^-分别表示参数容差的上下限值。判断式（7）解中所有参数偏差是否在容差范围内即可定位到故障元件。考虑将电路故障状态验证转化为容差约束下的线性规划问题求解。

引入变量x_j=x_j+ε_j^-，带入式（7），得到：

为判断式（8）解的存在性，采用单纯形法的第一阶段方法进行求解。引入人工变量z_i，对疑似故障元件“和‘的参数偏差限值（-ε_k^-，ε_k⁺）和（-ε_l^-，ε_l⁺）进行修正，得到修正后的偏差限值为（-γ_k^-，γ_k⁺）和（-γ_l^-，γ_l⁺），得到辅助线性规划方程：

若修正后的式（9）存在最小目标值ZO，则表示元件k和l确实发生故障。电路中存在更多故障的诊断过程与此类似。

4 实例验证

本文提出的算法使用大量BJT和CMOS电路进行测试。为了说明算法的有效性，考虑两个数值例子。在所有例子中，BJT电路的步长为0.01，CMOS的步长为0.05，M=20000。

4.1 案例一

在实验室条件下搭建的BJT电路[14]如图1所示。该电路含5个BC108B晶体管，测试电压精度为1mV。考虑到电阻Rz.R5和R6的参数变化对输出电压的影响很小，故障无法检测，因此仅考虑电阻R₁、R₃和R₄可能发生故障情况。实际测量电阻值为R₂=5.024kΩ，R₅=10.12kΩ，R₆=99.8kΩ。采用算法仿真计算时，采用它们的标称值。晶体管采用G-P模型，参数设置如下：IS=7.59×10^-15A，BF=480，NF=1，VAF=19.3V，IKF=0.071A，ISE=7.477×10^-15A，NE=1.2，BR=5，NR=1，VAR=12.5V，ISC=2.00×10^-13A，NC=1.3808，RB=0Ω，RE=0Ω，RC=0Ω。选择△u_i=5mV。测试激励信号如表1所示，考虑两种故障情况，实验测量电压见表1。采用本文提出算法，计算得到结果如表2所示。可以看出，计算电阻范围基本涵盖了故障电阻实值，计算得最大误差为4.35%。

如果提高电压测量精度到0.01mV，更多数量的参数故障能够被诊断出来。考虑R₁、R₂、R₃、R₄和R₅参数，激励信号、测量电压及计算结果如表3所示。可以看出，計算电阻范围基本涵盖了故障电阻实值，计算得最大误差为1.38%。在上述讨论的例子中，计算得到的结果都在曲线l^（1）上。

仿真实验包括30个不同的故障模式，实验结果为：计算得到19个故障模式的参数值位于曲线l^（1）上，11个故障模式的参数值位于曲线l^（2）上。

4.2 案例二

考虑如图2所示的含CMOS电压跟随器电路[15]。MOS晶体管仿真采用PSPICE中的LEVEL 1模型，标称参数设置：pMOS-GAMMA=0.53 V^0.5，IS=10_-16A，KP=28.3 μA/V²，PHI=0.58V，RD=RS=94Ω，VTO=-0.84 V：nMOS-GAMMA=0.38V^0.5，IS=10^-16A，KP=79.7μA/V²，PHI=0.53V，RD=RS=63Ω，VTO=0.79V。

晶体管标称每米的通道宽度和长度如图2所示。考虑沟道的宽度与长度的比值（W/L）偏差。晶体管的标称比值（W/L）为：3.750，8.750，50.000，15.000，15.000，6.875，6.875，75.000。

假设测点为1和2，测量精度取0.1mV，测试安排如表4所示。假定晶体管T1、T5、T6、T7、T8可能发生故障，晶体管T2、T3、T4的宽度与长度的比值（W/L）取标称值，故障参数值位于曲线l^（2）上，结果如表5所示。可以看出，W/L，计算范围基本涵盖了实值，计算得最大误差为1.08%。

将测量精度提高到0.001mV，测试安排如表6所示。假设所有的晶体管可能发生故障，该方法得到参数值范围（位于曲线l^（2）上）如表7所示。可以看出，W/L计算范围基本涵盖了实值，计算得最大误差为1.16%。

仿真实验包括30个不同故障模式，实验结果为：计算得到24个故障模式的参数值位于曲线l^（1）上，6个故障模式的参数值位于曲线l（2）上。

5 结束语

模拟电路参数故障诊断的关键是求解可能存在多解的非线性测试方程组，并从所有解中辨识出真解。本文研究将非线性测试方程的多解求解转化为差分方程的初值问题，提出一种系统搜索算法求解故障元件参数，将故障验证转化为线性规划问题最优解的存在性判断，BJT和CMOS电路实验验证了算法的有效性。本方法能够有效处理非线性模拟电路的软故障诊断问题，后续将进一步完善工程化的试验研究和对比分析。

参考文献

[1]TADEUSIEWICZ M，HALGAS S，KUCZYNSKI A.Newaspects of fault diagnosis of nonlinear analog circuits[J].International Journal of Electronics and Telecommunications，2015，61（1）：83-93.

[2]唐小峰，许爱强，李文海.多层次模拟电路随机变异生成方法研究[J].中国测试，2016，42（1）：12-15，64.

[3]周启忠，谢永乐，毕东杰，等一种模拟电路参数型故障诊断新方法[J].电子科技大学学报，2017（2）：386-391.

[4]XIE Y，LI X，XIE S，et al.Soft fault diagnosis of analogcircuits via frequency response function measurements[J].Journal of Electronic Testing，2014，30（2）：243-249.

[5]RAPISARDA L A，CARLO R D.Analog multi -frequencyfault diagnosis[J].IEEE Transactions on Circuits andSystems，1983，30（4）：223-234.

[6]BANDLER J W，BIERNACKI R M，SALAMA A E，et al.Fault isolation in linear analog circuits using the Llnorm[C]//Proc.IEEE Int.Symp.Circuits and Systems，1982：1140-1143.

[7]HE Y，SUN Y C.Neural network-based l1-norm optimisationapproach for fault diagnosis of nonlinear circuits withtolerance[J].IEE Proceedings-Circuits Devices and Systems，2001，148（4）：223-228.

[8]GUO Z，SAVIR J.Coefficient-based test of parametric faultsin analog circuits[J].IEEE Transactions on Instrumentationand Measurement，2006，55（1）：150-157.

[9]LONG T，WANG H，LONG B.A classical parameteridentification method and a modern test generationalgorithm[J].Circuits Systems and Signal Processing，2011，30（2）：391-412.

[10]BRANIN F H.Widely convergent method for finding multiplesolutions of simultaneous nonlinear equations[J].IBMJournal of Research and Development，1972，16（5）：504-522.

[11]CHAO K，LIU D，PAN C.A systematic search method forobtaining multiple solutions of simultaneous nonlinearequations[J].IEEE Transactions on Circuits and Systems，1975，22（9）：748-753.

[12]CHUA L O，LIN P M.Computer-aided analysis of electroniccircuits：algorithms and computational techniques[M].UpperSaddle River：Prentice Hall，1975.

[13]董海迪，劉刚，何兵，等.容差模拟电路参数故障诊断[J].中国测试，2017，43（9）：128-133.

[14]TADEUSIEWICZ M，HALGAS S.A new approach tomultiple soft fault diagnosis of analog BJT and CMOScircuits[J].IEEE Transactions on Instrumentation andMeasurement，2015，64（10）：2688-2695.

[15]TADEUSIEWICZ M，HALGAS S.Multiple soft faultdiagnosis of nonlinear circuits using the continuationmethod[J].Journal of Electronic Testing，2012，28（4）：487-493.

（编辑：李刚）