CN112485650B - 基于pbi的模拟电路故障参数范围识别方法 - Google Patents

Abstract

本发明公开了一种基于PBI的模拟电路故障参数范围确定方法，将元件参数向量作为遗传算法种群的个体，在生成初始种群的时候，故障元件的参数值在预设的故障取值范围中取值，其余元件在容差范围内取值，遗传算法迭代过程中，采用PBI计算个体适应度值，采用两次遗传算法分别在迭代完成后最后一代种群提取出故障元件参数的上限和上限。本发明结合遗传算法和PBI方法，实现对于故障元件参数范围的精确确定。

Description

基于PBI的模拟电路故障参数范围识别方法

技术领域

本发明属于模拟电路故障诊断技术领域，更为具体地讲，涉及一种基于PBI的模拟电路故障参数范围识别方法。

背景技术

在模拟电路工作过程中，元件退化会导致性能退化，及时对元件参数进行估计能够预防功能故障的发生。模拟电路发生故障时，除故障元件外，无故障元件参数是容差范围内的随机数，即所有元件参数都是变量。模拟集成电路测点数目有限，独立测试量的数目往往远远小于元件数目C，因此通过测试量和元件参数只能建立欠定方程组，无法精确计算出故障元件参数值。但是能够根据电路结构(传输函数)和容差范围，获得故障参数一个可能的故障范围。为电路性能退化预测提供支持。假定传输函数H(X)＝x₁x₂，x₁、x₂表示两个元件的参数值，两个元件标称值为10,则标准输出H为100。电路发生故障，测得输出为120，且已知故障源为x₁,则容易得到x₁＝12。考虑到无故障元件x₂容差(容差范围[10(1-0.05),10(1+0.05)])的影响，当x₂为容差下限9.5时，产生120的故障输出的x₁应为12.6；当x₂为容差上限10.5时，产生120的故障输出的x₁应为11.4。即在±5％的容差影响下,x₁在闭区间[11.4,12.6]任意取值都可能得到120的故障输出。当电路结构变得复杂，此闭区间的解析将很难精确计算，即难以确定故障参数范围。

发明内容

本发明的目的在于克服现有技术的不足，提供一种基于PBI的模拟电路故障参数范围确定方法，在遗传算法确定故障参数上限或下限过程中，采用PBI计算个体适应度，实现对于故障元件参数范围的精确确定。

为实现上述发明目的，本发明基于PBI的模拟电路故障参数范围确定方法包括以下步骤：

S1：获取模拟电路中元件数量C和各个元件的参数标称值x_iN，i＝1,2,…,C，确定模拟电路在预设测点的传输函数，获取当前故障状态下模拟电路在预设测

点的故障电压相量

以及所检测出的故障元件序号c；

S2：基于PBI确定故障元件参数下限，具体步骤包括：

S2.1：将元件参数向量X＝[x₁,x₂,…,x_C]作为遗传算法种群的个体，生成N个个体构成初始种群P，具体方法为：根据需要设置故障元件c的参数值x_c的故障取值范围

初始种群P中每个个体中故障元件c的参数值x_c在故障取值范围内取值，其余元件i′的参数x_i′在容差范围[x_i′N×(1-α),x_i′N×(1+α)]内取值，其中x_i′N表示元件i′的参数标称值，i′＝1,2,…,C&i′≠c；

S2.2：初始化迭代次数t＝1；

S2.3：对种群P中的个体进行交叉和变异，得到新种群Q，在交叉和变异过程中需要保证故障元件c的参数值x_c在故障取值范围内取值，非故障元件的参数值在容差范围内取值；

S2.4：将种群P和种群Q进行合并，得到合并种群S，即S＝P∪Q；

S2.5：分别计算2N个个体中每个个体X_k所对应的适应度值F_k，其中k＝1,2,…,2N，适应度值越小个体越优；适应度值F_k的具体计算方法为：

根据传输函数计算个体所对应元件参数向量在预设测点处的输出电压相量

然后计算该输出电压相量

与故障电压相量

之间的误差

将误差进行归一化得到归一化后的误差

其中，E_max、E_min分别表示2N个个体所对应误差中的最大值和最小值。

将每个个体X_k中故障元件c的参数值x_c,k进行归一化得到归一化后的参数值

其中，x_c,max、x_c,min分别表示2N个个体故障元件c的参数值中的最大值和最小值。

采用以下公式计算得到每个个体X_k的适应度值F_k：

S2.6：根据适应度值从2N个个体中优选出N个个体作为下一代种群P′；

S2.7：判断迭代次数t是否达到预设的最大迭代次数t_max，如果未达到，则进入步骤S2.8，否则进入步骤S2.9；

S2.8：令种群P＝P′，t＝t+1，返回步骤S2.3；

S2.9：将当前种群P′中误差最小个体中故障元件c的参数值作为故障元件c的参数范围下限x_cL；

S3：采用与步骤S2中的相同算法流程确定故障元件参数上限x_cU，在算法执行过程中采用以下公式计算个体X_k的适应度值F_k：

适应度值越小个体越优；

算法执行完毕后，将最后一代种群P′中误差最小个体中故障元件c的参数值作为故障元件c的参数范围下限x_cU。

本发明基于PBI的模拟电路故障参数范围确定方法，将元件参数向量作为遗传算法种群的个体，在生成初始种群的时候，故障元件的参数值在预设的故障取值范围中取值，其余元件在容差范围内取值，遗传算法迭代过程中，采用PBI计算个体适应度值，采用两次遗传算法分别在迭代完成后最后一代种群提取出故障元件参数的上限和上限。本发明结合遗传算法和PBI方法，实现对于故障元件参数范围的精确确定。

附图说明

图1是PBI优化示意图；

图2是本发明基于PBI的模拟电路故障参数范围确定方法的具体实施方式流程图；

图3是本发明中基于PBI确定故障元件参数下限的流程图；

图4是本实施例中二阶托马斯模拟滤波电路的电路图；

图5是本实施例中基于PBI确定故障元件参数下限的最后一代种群示意图。

具体实施方式

下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。

为了更好地说明本发明的技术方案，首先对本发明的技术思路进行简要说明。

假定模拟电路在预设测点上的传输函数为h(s,X)，其中s＝jω，j表示虚数单位，ω表示角频率，X表示元件参数向量，X＝[x₁,x₂,…,x_C]，x_i表示第i个元件的参数，i＝1,2,…,C，C表示模拟电路的元件数量。那么传输函数h(s,X)可以表示为：

其中，

表示输入电压相量，

表示测点的输出电压相量，a_n,a_n-1,…,a₀、b_m,b_m-1,…,b₀表示以元件参数值为变量的函数。

如果频率不变，那么传输函数h(s,X)仅由参数向量X确定。如果将输入电压相量

作为参考向量，即

那么：

如果模拟电路实测电压为

故障编号为f，其故障元件参数用x_f表示，则故障参数范围识别就是找到满足以下约束的最大故障元件参数值x_U和最小故障元件参数值x_L：

那么确定最小故障元件参数值x_L可以表达成：

其中，x_iN表示第i个元件的参数标称值，α表示容差参数，其取值范围一般为α∈(0,0.05]。

那么确定最大故障元件参数值x_U可以表达成：

将约束优化问题转化为双目标优化问题，式(4)和式(5)可以重写为：

其中，

E(X)表示输出电压相量和故障电压相量的误差，当x_f取相反数-x_f时，最小化-x_f就能实现求取最大故障元件参数值x_U。

基于惩罚的边界交叉法(penalty-based boundary intersection，PBI)定义如下：

其中，

为理想点，

表示第j个目标处的极小值。X为自变量，W为参考向量，θ为惩罚因子。d₁是作为X对PF面的收敛性指标，d₂作为度量种群多样性的指标。g(X|W,Z^*)＝d₁+θd₂作为收敛性和多样性的X的复合度量。PBI方法的目标是尽可能地降低目标函数，使其尽可能的达到PF面。

图1是PBI优化示意图。如图1所示，就本发明所要解决的问题而言，式(6)中需要E(X)接近0，而x_f尽量小。在E(X)-x_f平面上可以表达为沿着权W＝[0,1]最小化。此时，式(7)中d₁就是归一化之后的故障元件参数x_f，d₂就是归一化之后的误差E(X)。

因此，最小参数识别问题就可以表达为：

其中，

分别表示故障元件参数x_f、误差E(X)的归一化值。

基于以上思路，提出本发明基于PBI的模拟电路故障参数范围确定方法。图2是本发明基于PBI的模拟电路故障参数范围确定方法的具体实施方式流程图。如图2所示，本发明基于PBI的模拟电路故障参数范围确定方法的具体步骤包括：

S201：获取模拟电路故障数据：

获取模拟电路中元件数量C和各个元件的参数标称值x_iN，i＝1,2,…,C，确定模拟电路在预设测点的传输函数，获取当前故障状态下模拟电路在预设测点的故障电压相量

以及所检测出的故障元件序号c。

S202：基于PBI确定故障元件参数下限：

接下来基于PBI确定故障元件参数下限。图3是本发明中基于PBI确定故障元件参数下限的流程图。如图3所示，本发明中基于PBI确定故障元件参数下限的具体步骤包括：

S301：初始化遗传算法种群：

将元件参数向量X＝[x₁,x₂,…,x_C]作为遗传算法种群的个体，生成N个个体构成初始种群P，具体方法为：根据需要设置故障元件c的参数值x_c的故障取值范围

初始种群P中每个个体中故障元件c的参数值x_c在故障取值范围内取值，其余元件i′的参数x_i′在容差范围[x_i′N×(1-α),x_i′N×(1+α)]内取值，其中x_i′N表示元件i′的参数标称值，i′＝1,2,…,C&i′≠c。本实施例中故障元件c的参数值x_c的故障取值范围为[x_cN×10^-3,x_cN×10³]。

S302：初始化迭代次数t＝1。

S303：交叉变异：

对种群P中的个体进行交叉和变异，得到新种群Q，在交叉和变异过程中需要保证故障元件c的参数值x_c在故障取值范围内取值，非故障元件的参数值在容差范围内取值。

本实施例中个体交叉采用模拟二进制交叉，变异采用多项式变异，交叉率和变异率根据实际需要设置即可。

S304：合并种群：

将种群P和种群Q进行合并，得到合并种群S，即S＝P∪Q。显然合并种群S中包含2N个个体。

S305：基于PBI计算个体适应度值：

分别计算2N个个体中每个个体X_k所对应的适应度值F_k，其中k＝1,2,…,2N，适应度值F_k的具体计算方法为：

根据传输函数计算个体所对应元件参数向量在预设测点处的输出电压相量

然后计算该输出电压相量

与故障电压相量

之间的误差

将误差进行归一化得到归一化后的误差

其中，E_max、E_min分别表示2N个个体所对应误差中的最大值和最小值。

将每个个体X_k中故障元件c的参数值x_c,k进行归一化得到归一化后的参数值：

其中，x_c,max、x_c,min分别表示2N个个体故障元件c的参数值中的最大值和最小值。

采用以下公式计算得到每个个体X_k的适应度值F_k：

其中，θ为预设的惩罚因子。

显然适应度值越小个体越优。

S306：个体优选：

根据适应度值从2N个个体中优选出N个个体作为下一代种群P′。本实施例中采用二选一锦标赛优选算法进行个体优选。

S307：判断迭代次数t是否达到预设的最大迭代次数t_max，如果未达到，则进入步骤S308，否则进入步骤S309；

S308：令种群P＝P′，t＝t+1，返回步骤S303；

S309：确定故障参数下界：

将当前种群P′中误差最小个体中故障元件c的参数值作为故障元件c的参数范围下限x_cL。

S203：基于PBI确定故障元件参数上限：

基于PBI确定故障元件参数上限，其算法流程与步骤S202相同，但是由于此时需要确定故障元件参数上限，需要对适应度计算公式进行调整，即在算法执行过程中采用以下公式计算得到每个个体X_k的适应度值F_k：

同样地，适应度值越小个体越优。

算法执行完毕后，将最后一代种群P′中误差最小个体中故障元件c的参数值作为故障元件c的参数范围上限x_cU。

实施例

为了更好地说明本发明的技术方案和技术效果，采用一个具体模拟电路对本发明进行实验验证。图4是本实施例中二阶托马斯模拟滤波电路的电路图。如图4所示，本实施例中二阶托马斯模拟滤波电路包括6个电阻元件，2个电容以及3个放大器，各元件参数的标称值如图4中标示。本实施例中以V_out作为测点，其传输函数如下式所示：

本实施例中设置故障元件为电阻R₂，令其参数值为19kΩ，其它元件在容差范围(本实施例设置容差参数α＝0.05，则容差范围为(x_iN×95％,x_iN×105％))内随机取值。输入信号为有效值为1V，频率为1KHz的正弦信号。仿真得到故障电压相量

将元件参数向量X＝[x₁,x₂,…,x_C]作为遗传算法种群的个体，本实施例中C＝8，前6位为电阻，最后2位为电容。设置故障元件2(即电阻R₂)的参数值x₂的故障取值范围[1Ω,100MΩ]，初始种群P中每个个体中故障元件2的参数值在该故障取值范围内随机取值，其余元件i′的参数x_i′在容差范围(x_i′N×95％,x_i′N×105％)内随机取值。设置种群数目N＝200，最大迭代次数t_max＝400。

图5是本实施例中基于PBI确定故障元件参数下限的最后一代种群示意图。如图5所示，可以确定本实施例中故障元件参数下限x_2L＝17410Ω，此时其他元件的参数值分别为：R₁＝10087Ω，R₃＝9900Ω，R₄＝9900Ω，R₅＝9900Ω，R₆＝11000Ω，C₁＝9.9nF，C₂＝9.9nF。对应的误差E(X)＝1.3515e-06，接近零，精度完全满足要求。同理基于PBI确定故障元件参数上限x_2U＝20808Ω，此时其他元件的参数值分别为：R₁＝9073Ω，R₃＝11000Ω，R₄＝11000Ω，R₅＝11000Ω，R₆＝9900Ω，C₁＝11nF，C₂＝11nF。对应的误差E(X)＝1.1471e-07，接近零，精度完全满足要求。从而得到故障参数范围为[17410Ω,20808Ω]，即所有故障元件2的参数值在闭区间[17410Ω,20808Ω]内的故障都可以产生故障电压相量

显然，设定的故障R₂＝19kΩ也在此范围内。

尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

Claims (1)

1.一种基于PBI的模拟电路故障参数范围确定方法，其特征在于，包括以下步骤：

S1：获取模拟电路中元件数量C和各个元件的参数标称值x_iN，i＝1,2,…,C，确定模拟电路在预设测点的传输函数，获取当前故障状态下模拟电路在预设测点的故障电压相量

以及所检测出的故障元件序号c；

S2：基于PBI确定故障元件参数下限，具体步骤包括：

S2.1：将元件参数向量X＝[x₁,x₂,…,x_C]作为遗传算法种群的个体，生成N个个体构成初始种群P，具体方法为：根据需要设置故障元件c的参数值x_c的故障取值范围

初始种群P中每个个体中故障元件c的参数值x_c在故障取值范围内取值，其余元件i′的参数x_i′在容差范围[x_i′N×(1-α),x_i′N×(1+α)]内取值，其中x_i′N表示元件i′的参数标称值，i′＝1,2,…,C&i′≠c，α表示容差参数，取值范围为α∈(0,0.05]；

S2.2：初始化迭代次数t＝1；

S2.3：对种群P中的个体进行交叉和变异，得到新种群Q，在交叉和变异过程中需要保证故障元件c的参数值x_c在故障取值范围内取值，非故障元件的参数值在容差范围内取值；

S2.4：将种群P和种群Q进行合并，得到合并种群S，即S＝P∪Q；

S2.5：分别计算2N个个体中每个个体X_k所对应的适应度值F_k，其中k＝1,2,…,2N，适应度值越小个体越优；适应度值F_k的具体计算方法为：

根据传输函数计算个体所对应元件参数向量在预设测点处的输出电压相量

然后计算该输出电压相量

与故障电压相量

之间的误差

将误差进行归一化得到归一化后的误差

其中，E_max、E_min分别表示2N个个体所对应误差中的最大值和最小值；

将每个个体X_k中故障元件c的参数值x_c,k进行归一化得到归一化后的参数值

其中，x_c,max、x_c,min分别表示2N个个体故障元件c的参数值中的最大值和最小值；

采用以下公式计算得到每个个体X_k的适应度值F_k：

其中，θ为预设的惩罚因子；

S2.6：根据适应度值从2N个个体中优选出N个个体作为下一代种群P′；

S2.7：判断迭代次数t是否达到预设的最大迭代次数t_max，如果未达到，则进入步骤S2.8，否则进入步骤S2.9；

S2.8：令种群P＝P′，t＝t+1，返回步骤S2.3；

S2.9：将当前种群P′中误差最小个体中故障元件c的参数值作为故障元件c的参数范围下限x_cL；

S3：采用与步骤S2中的相同算法流程确定故障元件参数上限x_cU，在算法执行过程中采用以下公式计算个体X_k的适应度值F_k：

适应度值越小个体越优；

算法执行完毕后，将最后一代种群P′中误差最小个体中故障元件c的参数值作为故障元件c的参数范围上限x_cU。

Images