CN110308386B - 基于多项式拟合与状态监测的模拟电路故障诊断方法 - Google Patents
基于多项式拟合与状态监测的模拟电路故障诊断方法 Download PDFInfo
- Publication number
- CN110308386B CN110308386B CN201910233050.7A CN201910233050A CN110308386B CN 110308386 B CN110308386 B CN 110308386B CN 201910233050 A CN201910233050 A CN 201910233050A CN 110308386 B CN110308386 B CN 110308386B
- Authority
- CN
- China
- Prior art keywords
- fault
- matrix
- representing
- analog circuit
- representative
- 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.)
- Expired - Fee Related
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000012544 monitoring process Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 43
- 239000013598 vector Substances 0.000 claims abstract description 37
- 230000005540 biological transmission Effects 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000011156 evaluation Methods 0.000 claims description 13
- 230000005284 excitation Effects 0.000 claims description 11
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/28—Testing of electronic circuits, e.g. by signal tracer
- G01R31/316—Testing of analog circuits
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tests Of Electronic Circuits (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
本发明公开了一种基于多项式拟合与状态监测的模拟电路故障诊断方法,获取测点的传输函数,分析模拟电路的模糊组信息,确定每个模糊组的代表故障元件,基于传输函数得到各个代表故障元件的特征矩阵,构建多项式拟合的超定方程组,计算得到各个代表故障元件对应的系数向量,当模拟电路发生故障时,对测点进行状态监测,获得测试矩阵和常数项矩阵,求使得测试矩阵与各个代表故障元件特征向量确定的方程具有最小二乘距离的特征向量,对应的代表故障元件即为故障诊断结果。采用本发明可以有效实现对模拟电路的故障诊断。
Description
技术领域
本发明属于模拟电路故障诊断技术领域,更为具体地讲,涉及一种基于多项式拟合与状态监测的模拟电路故障诊断方法。
背景技术
随着集成电路的快速发展,为了提高产品性能、降低芯片面积和费用,需将数字和模拟元件集成在同一块芯片上。据资料报道,虽然模拟部分仅占芯片面积的5%,但其故障诊断成本却占总诊断成本的95%,模拟电路故障诊断一直是集成电路工业中的一个“瓶颈”问题。现阶段已经有一些发展的比较完善的模拟电路故障诊断理论应用到实际中了,例如:测前模拟诊断法中的故障字典法测后模拟诊断法中的元件参数辨识法和故障验证法。但这些方法仅能处理离散的参数故障和硬故障,不能完整诊断模拟元件的连续参数故障。复数域圆模型能够完整建模模拟元件的所有参数漂移故障,是一种软硬统一的故障诊断模型。
图1是二阶托马斯模拟滤波电路的拓扑图。如图1所示,Vout为测点,该电路模糊组情况为:{R1}、{R2}、{R3,C1}、{R4,R5,R6,C2}。模糊组内部元件的故障不可区分。模糊组之间的故障理论上都能被区分。每个模糊组选出一个元件作为代表性故障元件。图2是图1所示电路中代表性故障元件的特征曲线。他们的公共交点为无故障点。除去此点,各曲线无交点,即所有故障的所有参数故障产生的输出电压都不同,是可以被区分的。但是如果考虑无故障元件存在±5%的容差,则每个代表性故障的特征曲线应该在一定范围内变化。图3是图2中代表性故障元件含容差的特征区域。如图3所示,此时特征曲线变成了特征区域,各特征区域存在重叠区域。图4是图3中特征区域的重叠区域标识图。图4中黑色部分表示特征区域的重叠区域。重叠区域说明原本能够区分的不同故障源的故障,由于容差影响,导致这些故障具有相同的特征,从而不可分。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于多项式拟合与状态监测的模拟电路故障诊断方法,对各个模糊组中代表故障元件的数据进行多项式拟合得到系数向量作为特征向量,在模拟电路故障时通过状态监测得到测试矩阵,通过与各个系数向量的最小二乘距离得到故障诊断结果。
为实现上述发明目的,本发明基于多项式拟合与状态监测的模拟电路故障诊断方法包括以下步骤:
S1:获取模拟电路在测点t处的传输函数;
S2:分析得到模拟电路经测点t输出电压进行故障诊断的模糊组信息,将得到的模糊组数量记为N,每个模糊组选取一个元件作为代表故障元件;
S3:对于每个代表故障元件,分别采用以下方法获取代表故障元件的特征矩阵:
令第i个代表故障元件的参数值pi在可能取值范围[pi min,pi max]中随机取M个值,pi min、pi max分别表示pi的最小可能取值和最大可能取值,i=1,2,…,N,其余故障元件参数值在容差范围内取值,根据传输函数得到预设激励信号下第i个代表故障元件第m次取值对应的传输函数值根据每个代表故障元件的M个传输函数值构建大小为M×2的代表故障元件特征矩阵Hi:
S4:在每个代表故障元件特征矩阵Hi的最后一列***一列M维单位列向量,得到大小为M×3的扩展特征矩阵Ai:
如果Ai的秩r(Ai)=3,则令第i个代表故障元件的系数向量Ki:
Ki=[ki,1,ki,2,ki,3]T
令第i个代表故障元件的常数项矩阵bi:
求解超定方程组AiKi=bi的系数向量Ki的最小二乘解:
如果矩阵Ai的秩r(Ai)=2,则构建如下矩阵A′i:
令第i个代表故障元件的系数向量Ki:
K′i=[ki,2,ki,3]T
令第i个代表故障元件的常数项矩阵bi:
求解超定方程组A′iK′i=bi的系数向量K′i的最小二乘解:
得到系数向量Ki=[0,ki,2,ki,3]T;
S5:当模拟电路发生故障时,在与步骤S3相同的激励信号下对测点t处的输出电压进行状态监测,测量得到D个输出电压其中d=1,2,…,D,D>3;计算得到 表示激励信号电压,构成大小为D×3的测试矩阵以及大小为D×2的常数项矩阵
本发明基于多项式拟合与状态监测的模拟电路故障诊断方法,获取测点的传输函数,分析模拟电路的模糊组信息,确定每个模糊组的代表故障元件,基于传输函数得到各个代表故障元件的特征矩阵,构建多项式拟合的超定方程组,计算得到各个代表故障元件对应的系数向量,当模拟电路发生故障时,对测点进行状态监测,获得测试矩阵和常数项矩阵,求使得测试矩阵与各个代表故障元件特征向量确定的方程具有最小二乘距离的特征向量,对应的代表故障元件即为故障诊断结果。采用本发明可以有效实现对模拟电路的故障诊断。
附图说明
图1是二阶托马斯模拟滤波电路的拓扑图;
图2是图1所示电路中代表性故障元件的特征曲线;
图3是图2中代表性故障元件含容差的特征区域;
图4是图3中特征区域的重叠区域标识图;
图5是本发明基于多项式拟合与状态监测的模拟电路故障诊断方法的具体实施方式流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
为了更好地说明本发明的技术方案,首先对本发明的技术原理进行说明。
众所周知,模拟电路任意测点相对输入点的传输函数的表达式为:
其中,s=jω,j表示虚数单位,ω表示角频率,an、bm都是以元件参数值为变量的函数。
假如频率固定,单故障假设下,只有故障源元件的参数值x未知,其它元件参数值都是标称值,此时,传输函数变成x的函数,可以用双线性方程表达为:
其中,A、B、C、D都是常复数,显然H(x)也为复数,其实部hr、虚部hj均为以故障源元件参数值x为变量的函数:
hr=f1(x) (3)
hj=f2(x) (4)
将H(x)作为故障特征,在复平面上画出特征曲线,就是要表达出实部hr和虚部hj的关系函数,即联立(3)和(4),消去x,获得hr和hj的关系方程:
F(hr,hj)=0 (5)
众所周知公式(2)所示的双线性变换将直线或者圆方程变换为圆方程(如公式(6)所示)和直线方程(如公式(7)所示):
因此,只要确定了故障特征参数K=[k1,k2,k3],则故障源特征就唯一确定。
图5是本发明基于多项式拟合与状态监测的模拟电路故障诊断方法的具体实施方式流程图。如图5所示,本发明基于多项式拟合与状态监测的模拟电路故障诊断方法的具体步骤包括:
S501:获取传输函数:
获取模拟电路在测点t处的传输函数。
S502:模糊组分析:
分析得到模拟电路经测点t输出电压进行故障诊断的模糊组信息,将得到的模糊组数量记为N,每个模糊组选取一个元件作为代表故障元件。
S503:获取代表故障元件特征矩阵:
对于每个代表故障元件,分别采用以下方法获取代表故障元件的特征矩阵:
令第i个代表故障元件的参数值pi在可能取值范围[pi min,pi max]中随机取M个值,pi min、pi max分别表示pi的最小可能取值和最大可能取值,i=1,2,…,N,可能取值范围中包含正常取值和故障取值,一般为电路实际运行过程时元件参数值的取值范围,以电阻为例,可能会发生短路或断路,那么其可能取值范围即为[0,+∞)。令其余故障元件参数值在容差范围内取值,根据传输函数得到预设激励信号下第i个代表故障元件第m次取值对应的传输函数值m=1,2,…,M。根据每个代表故障元件的M个传输函数值构建大小为M×2的代表故障元件特征矩阵Hi:
S504:确定代表故障元件特征向量:
在每个代表故障元件特征矩阵Hi的最后一列***一列M维单位列向量,得到大小为M×3的扩展特征矩阵Ai:
那么公式(6)和(7)可以统一表达为多项式拟合的超定方程组:
AiKi=bi (10)
其中,Ki表示第i个代表故障元件的系数向量,bi表示第i个代表故障元件的常数项矩阵。
如果Ai的秩r(Ai)=3(对应圆方程),则令第i个代表故障元件的系数向量Ki:
Ki=[ki,1,ki,2,ki,3]T (11)
令第i个代表故障元件的常数项矩阵bi:
求解超定方程组AiKi=bi的系数向量Ki的最小二乘解:
Ki=(Ai TAi)(-1)Ai Tbi (13)
如果矩阵Ai的秩r(Ai)=2(对应直线方程),则构建如下矩阵A′i:
同时令第i个代表故障元件的系数向量K′i:
K′i=[ki,2,ki,3]T (15)
令第i个代表故障元件的常数项矩阵bi:
求解超定方程组A′iK′i=bi的系数向量K′i的最小二乘解:
得到系数向量Ki=[0,ki,2,ki,3]T。
每个系数向量Ki即为对应代表故障元件的特征向量。
S505:基于状态监测进行故障诊断:
当模拟电路发生故障时,在与步骤S503相同的激励信号下对测点t处的输出电压进行状态监测,测量得到D个输出电压其中d=1,2,…,D,D需要大于系数向量的维数,即D>3,其值越大越好,可以根据实际需要设置。计算得到 表示激励信号电压,构成大小为D×3的测试矩阵以及大小为D×2的常数项矩阵
从所有N个代表故障元件的特征向量Ki中,找出与测试矩阵最为相符的一个,即求使得测试矩阵与特征向量确定的方程具有最小二乘距离的特征向量,具体方法为:计算评价参数选取N个评价参数Wi中最小值所对应的系数向量Ki,其所对应的代表故障元件即为故障诊断结果,也就是说,诊断得到的代表故障元件序号可以用如下公式表示:
实施例
为了更好地说明本发明的技术效果,采用如图1所示的二阶托马斯模拟滤波电路为例对本发明进行说明。本实施例中二阶托马斯模拟滤波电路以Vout作为测点,该测点下的模糊组情况为:{R1}、{R2}、{R3,C1}、{R4,R5,R6,C2}。模糊组是由电路结构确定的,与激励信号无关,只与测点选取相关。标称情况下,直流电源为5V,在1V,1kHZ正弦信号激励下,采用本发明方法获得代表故障元件的特征向量。表1是本实施例中代表故障元件的特征向量。
模糊组 | k<sub>1</sub> | k<sub>2</sub> | k<sub>3</sub> | |
{R<sub>1</sub>} | K<sub>1</sub> | 0 | 0.9632 | 0 |
{R<sub>2</sub>} | K<sub>2</sub> | 1.6523 | 0 | 0 |
{R<sub>3</sub>,C<sub>1</sub>} | K<sub>3</sub> | 0 | -1.5915 | 0 |
{R<sub>4</sub>,R<sub>5</sub>,R<sub>6</sub>,C<sub>2</sub>} | K<sub>4</sub> | 1 | -0.6283 | 0 |
表1
从表1中可知,模糊组{R1}、{R3,C1}所对应的模型为直线,模糊组{R2}、{R4,R5,R6,C2}
为对本发明进行验证,本实施例中模拟两个故障:元件R1故障和元件R5故障。
·元件R1故障
设置容差范围为±5%,元件R1的故障取值范围为[10-3×10kΩ,0.95×10kΩ)∪(1.05×10kΩ,103×10kΩ],在故障取值范围内随机取值50次,构建大小为50×3的测试矩阵和大小为50×2的矩阵然后计算得到4个评价参数Wi。表2是本实施例中元件R1故障的评价参数列表。
K<sub>1</sub> | K<sub>2</sub> | K<sub>3</sub> | K<sub>4</sub> |
162.4689 | 4.2379e+03 | 3.8724e+03 | 4.0916e+03 |
表2
如表2所示,根据评价参数可以得到特征向量K1所对应的代表故障元件为诊断结果,即故障元件为模糊组{R1}中的元件,与实际相符。
·元件R5故障
同样地,设置容差范围为±5%,元件R5的故障取值范围为[10-3×10kΩ,0.95×10kΩ)∪(1.05×10kΩ,103×10kΩ],在故障取值范围内随机取值50次,构建大小为50×3的测试矩阵和大小为50×2的矩阵然后计算得到4个评价参数Wi。表3是本实施例中元件R5故障的评价参数列表。
K<sub>1</sub> | K<sub>2</sub> | K<sub>3</sub> | K<sub>4</sub> |
4.4358 | 10.0151 | 26.2568 | 0.0148 |
表3
如表3所示,根据评价参数可以得到特征向量K4所对应的代表故障元件为诊断结果,即故障元件为模糊组{R4,R5,R6,C2}中的元件,与实际相符。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种基于多项式拟合与状态监测的模拟电路故障诊断方法,其特征在于,包括以下步骤:
S1:获取模拟电路在测点t处的传输函数;
S2:分析得到模拟电路经测点t输出电压进行故障诊断的模糊组信息,将得到的模糊组数量记为N,每个模糊组选取一个元件作为代表故障元件;
S3:对于每个代表故障元件,分别采用以下方法获取代表故障元件的特征矩阵:
令第i个代表故障元件的参数值pi在可能取值范围[pimin,pimax]中随机取M个值,pimin、pimax分别表示pi的最小可能取值和最大可能取值,i=1,2,…,N,其余故障元件参数值在容差范围内取值,根据传输函数得到预设激励信号下第i个代表故障元件第m次取值对应的传输函数值j表示虚数单位;根据每个代表故障元件的M个传输函数值构建大小为M×2的代表故障元件特征矩阵Hi:
S4:在每个代表故障元件特征矩阵Hi的最后一列***一列M维单位列向量,得到大小为M×3的扩展特征矩阵Ai:
如果Ai的秩r(Ai)=3,则令第i个代表故障元件的系数向量Ki:
Ki=[ki,1,ki,2,ki,3]T
令第i个代表故障元件的常数项矩阵bi:
求解超定方程组AiKi=bi的系数向量Ki的最小二乘解:
Ki=(Ai TAi)(-1)Ai Tbi
如果矩阵Ai的秩r(Ai)=2,则构建如下矩阵A′i:
令第i个代表故障元件的系数向量K′i:
K′i=[ki,2,ki,3]T
令第i个代表故障元件的常数项矩阵bi:
求解超定方程组A′iK′i=bi的系数向量K′i的最小二乘解:
得到系数向量Ki=[0,ki,2,ki,3]T;
S5:当模拟电路发生故障时,在与步骤S3相同的激励信号下对测点t处的输出电压进行状态监测,测量得到D个输出电压其中d=1,2,…,D,D>3;计算得到 表示激励信号电压,构成大小为D×3的测试矩阵以及大小为D×2的常数项矩阵
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910233050.7A CN110308386B (zh) | 2019-03-26 | 2019-03-26 | 基于多项式拟合与状态监测的模拟电路故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910233050.7A CN110308386B (zh) | 2019-03-26 | 2019-03-26 | 基于多项式拟合与状态监测的模拟电路故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110308386A CN110308386A (zh) | 2019-10-08 |
CN110308386B true CN110308386B (zh) | 2020-09-18 |
Family
ID=68074776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910233050.7A Expired - Fee Related CN110308386B (zh) | 2019-03-26 | 2019-03-26 | 基于多项式拟合与状态监测的模拟电路故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110308386B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113156303A (zh) * | 2021-05-16 | 2021-07-23 | 桂林电子科技大学 | 基于输出响应矩阵特性分析的模拟电路故障诊断和定位法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3670553B2 (ja) * | 2000-03-27 | 2005-07-13 | 株式会社東芝 | 半導体集積回路解析装置、半導体集積回路解析方法及び半導体集積回路解析方法を実行するためのプログラムを記録した記録媒体 |
JP2004061415A (ja) * | 2002-07-31 | 2004-02-26 | Agilent Technologies Japan Ltd | デバイスの特性試験方法 |
CN101900789B (zh) * | 2010-07-07 | 2012-05-02 | 湖南大学 | 基于小波变换与分形维数的容差模拟电路故障诊断方法 |
CN105548862B (zh) * | 2016-01-25 | 2019-02-05 | 合肥工业大学 | 一种基于广义多核支持向量机的模拟电路故障诊断方法 |
-
2019
- 2019-03-26 CN CN201910233050.7A patent/CN110308386B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110308386A (zh) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110308384B (zh) | 基于圆模型和神经网络的模拟电路故障诊断方法 | |
CN109581203B (zh) | 基于遗传算法的测后仿真故障诊断方法 | |
CN110907810B (zh) | 基于粒子群算法的模拟电路单故障诊断方法 | |
JP4783394B2 (ja) | センサの故障適応方法 | |
US9322871B2 (en) | Current measurement circuit and method of diagnosing faults in same | |
CN110470979B (zh) | 基于故障特征区域的模拟电路故障诊断方法 | |
CN109839583B (zh) | 基于改进遗传算法的模拟电路多故障诊断方法 | |
US20160061349A1 (en) | Device and Method for Valve Signature Testing | |
CN112485652B (zh) | 基于改进正余弦算法的模拟电路单故障诊断方法 | |
CN111260063A (zh) | 基于遗传算法的模拟电路故障定位与参数辨识方法 | |
EP3690658B1 (en) | Method for detecting repair-necessary motherboards and device using the method | |
CN110673017B (zh) | 基于遗传算法的模拟电路故障元件参数辨识方法 | |
CN104237770B (zh) | 一种模拟电路故障诊断方法 | |
CN110308386B (zh) | 基于多项式拟合与状态监测的模拟电路故障诊断方法 | |
CN111308327B (zh) | 模拟电路故障定位与故障元件参数辨识方法 | |
EP0992807B1 (en) | Method and apparatus for board model correction | |
CN112084910A (zh) | 故障诊断方法及*** | |
CN104090228B (zh) | 一种模拟电路模糊组识别方法 | |
CN112505532A (zh) | 基于改进粒子群算法的模拟电路单故障诊断方法 | |
CN107576904B (zh) | 基于锯齿波的模拟电路故障诊断方法 | |
CN109948267B (zh) | 基于圆模型参数的线性模拟电路故障诊断方法 | |
US20060155734A1 (en) | Apparatus and methods for evaluating a dynamic system | |
CN113189513B (zh) | 一种基于纹波的冗余电源均流状态识别方法 | |
CN113156303A (zh) | 基于输出响应矩阵特性分析的模拟电路故障诊断和定位法 | |
CN108664936B (zh) | 一种基于机器故障的诊断方法及*** |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200918 |