本发明属于模拟电路故障诊断技术领域，更为具体地讲，涉及一种基于切比雪夫的模拟电路元件故障参数识别方法。
背景技术：
在模拟电路工作过程中，元件退化会导致性能退化，及时对元件参数进行估计能够预防功能故障的发生。模拟电路发生故障时，除故障元件外，无故障元件参数是容差范围内的随机数，即所有元件参数都是变量。模拟集成电路测点数目有限，独立测试量的数目往往远远小于元件数目c，因此通过测试量和元件参数只能建立欠定方程组，无法精确计算出故障元件参数值。但是能够根据电路结构(传输函数)和容差范围，获得故障参数一个可能的故障范围。为电路性能退化预测提供支持。假定传输函数h(x)＝x1x2，x1、x2表示两个元件的参数值，两个元件标称值为10，则标准输出h为100。电路发生故障，测得输出为120，且已知故障源为x1,则容易得到x1＝12。考虑到无故障元件x2容差(容差范围[10(1-0.05),10(1 0.05)])的影响，当x2为容差下限9.5时，产生120的故障输出的x1应为12.6；当x2为容差上限10.5时，产生120的故障输出的x1应为11.4。即在±5％的容差影响下，x1在闭区间[11.4,12.6]任意取值都可能得到120的故障输出。当电路结构变得复杂，此闭区间的解析将很难精确计算，即难以确定故障参数范围。
技术实现要素：
本发明的目的在于克服现有技术的不足，提供一种基于切比雪夫的模拟电路元件故障参数识别方法，在遗传算法确定故障参数上限或下限过程中，采用切比雪夫函数值作为个体适应度值，实现对于故障元件参数范围的精确确定。
为实现上述发明目的，本发明基于切比雪夫的模拟电路元件故障参数识别方法包括以下步骤：
s1：获取模拟电路中元件数量c和各个元件的参数标称值xin，i＝1,2,…,c，确定模拟电路在预设测点的传输函数，获取当前故障状态下模拟电路在预设测点的故障电压相量
s2：基于切比雪夫确定故障元件参数下限，具体步骤包括：
s2.1：将元件参数向量x＝[x1,x2,…,xc]作为遗传算法种群的个体，生成n个个体构成初始种群p，具体方法为：根据需要设置故障元件c的参数值xc的故障取值范围
s2.2：预设范围[-λ2,-λ1]，λ1、λ2均为正常数，且λ2＞λ1；设置n个权重向量，记第n个权重向量wn＝(wn,1,wn,2)t，其中
s2.3：初始化迭代次数t＝1；
s2.4：对种群p中的个体进行交叉和变异，得到新种群q，在交叉和变异过程中需要保证故障元件c的参数值xc在故障取值范围内取值，非故障元件的参数值在容差范围内取值；
s2.5：将种群p和种群q进行合并，得到合并种群s，即s＝p∪q；
s2.6：分别计算2n个个体中每个个体xk到各权重向量的切比雪夫函数值，其中k＝1,2,…,2n；切比雪夫函数值的具体计算方法为：
根据传输函数计算个体所对应元件参数向量在预设测点处的输出电压相量
将每个个体xk中故障元件c的参数值xc,k进行归一化得到归一化后的参数值：
其中，xc,max表示2n个个体故障元件c的参数值中的最大值；
令f1(xk)＝e(xk)，
s2.7：对于每个权重向量wn，从2n个个体中选择到该权重向量的切比雪夫函数值最小的个体，作为该权重向量的最优个体，将n个权重向量对应的最优个体构成下一代种群p′；
s2.8：判断迭代次数t是否达到预设的最大迭代次数tmax，如果未达到，则进入步骤s2.9，否则进入步骤s2.10；
s2.9：令种群p＝p′，t＝t 1，返回步骤s2.4；
s2.10：将当前种群p′中具有最小误差的个体中故障元件c的参数值作为故障元件c的参数范围下限xcl；
s3：采用与步骤s2中的相同算法流程确定故障元件参数上限xcu，在算法执行过程中采用以下公式将每个个体xk中故障元件c的参数值xc,k进行归一化得到归一化后的参数值
算法执行完毕后，将最后一代种群p′中具有最小误差的个体中故障元件c的参数值作为故障元件c的参数范围上限xcu。
本发明基于切比雪夫的模拟电路故障参数范围确定方法，将元件参数向量作为遗传算法种群的个体，在生成初始种群的时候，故障元件的参数值在预设的故障取值范围中取值，其余元件在容差范围内取值，遗传算法迭代过程中，基于切比雪夫函数值进行个体优选，采用两次遗传算法分别在迭代完成后最后一代种群提取出故障元件参数的上限和上限。本发明结合遗传算法和切比雪夫函数，实现对于故障元件参数范围的精确确定。
附图说明
图1是本发明基于切比雪夫的模拟电路元件故障参数识别方法的具体实施方式流程图；
图2是本发明中基于切比雪夫确定故障元件参数下限的流程图；
图3是本实施例中二阶托马斯模拟滤波电路的电路图；
图4是本实施例中基于切比雪夫确定故障元件参数下限和上限的最后一代种群示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。
为了更好地说明本发明的技术方案，首先对本发明的技术思路进行简要说明。
假定模拟电路在预设测点上的传输函数为h(s,x)，其中s＝jω，j表示虚数单位，ω表示角频率，x表示元件参数向量，x＝[x1,x2,…,xc]，xi表示第i个元件的参数，i＝1,2,…,c，c表示模拟电路的元件数量。那么传输函数h(s,x)可以表示为：
其中，
如果频率不变，那么传输函数h(s,x)仅由参数向量x确定。如果将输入电压相量
如果模拟电路实测电压为
那么确定最小故障元件参数值xl可以表达成：
其中，xin表示第i个元件的参数标称值，α表示容差参数，其取值范围一般为α∈(0,0.05]。
那么确定最大故障元件参数值xu可以表达成：
将约束优化问题转化为双目标优化问题，式(4)和式(5)可以重写为：
其中，
切比雪夫(tchebycheff)方法定义如下：
其中，
基于以上思路，提出本发明基于切比雪夫的模拟电路元件故障参数识别方法。图1是本发明基于切比雪夫的模拟电路元件故障参数识别方法的具体实施方式流程图。如图1所示，本发明基于切比雪夫的模拟电路元件故障参数识别方法的具体步骤包括：
s101：获取模拟电路故障数据：
获取模拟电路中元件数量c和各个元件的参数标称值xin，i＝1,2,…,c，确定模拟电路在预设测点的传输函数，获取当前故障状态下模拟电路在预设测点的故障电压相量
s102：基于切比雪夫确定故障元件参数下限：
接下来基于切比雪夫确定故障元件参数下限。图2是本发明中基于切比雪夫确定故障元件参数下限的流程图。如图2所示，本发明中基于切比雪夫确定故障元件参数下限的具体步骤包括：
s201：初始化遗传算法种群：
将元件参数向量x＝[x1,x2,…,xc]作为遗传算法种群的个体，生成n个个体构成初始种群p，具体方法为：根据需要设置故障元件c的参数值xc的故障取值范围
s202：初始化权重向量：
由于本发明中种群个体数目为n，那么权重向量数目也为n。考虑到要让误差e(x)接近于0，因此一般需要设置权重向量中与误差e(x)对应的权重值小于10-2以满足精度要求。根据式(7)可知，权重值不能为0，但可以尽量小，如10-16。如果按照传统方法将[10-16,10-2]均匀分成n-1份，则会丢失接近0的细节部分，很难聚焦高精度区域，而本发明的目的是令e(x)接近于0。因此本发明对指数进行划分，例如将[-16,-2]均匀分成n-1份，然后以10为底取指数，这样在保证覆盖范围的同时又聚焦于高精度区域。初始化权重向量的具体方法为：
预设范围[-λ2,-λ1]，λ1、λ2均为正常数，且λ2＞λ1。设置n个权重向量，记第n个权重向量wn＝(wn,1,wn,2)t，其中
例如预设范围为[-16,-2]，n取15，则指数参数γn＝2,3,4,…,16，一共有15个权重向量，可以构成如下15×2的权重矩阵：
s203：初始化迭代次数t＝1。
s204：交叉变异：
对种群p中的个体进行交叉和变异，得到新种群q，在交叉和变异过程中需要保证故障元件c的参数值xc在故障取值范围内取值，非故障元件的参数值在容差范围内取值。
本实施例中个体交叉采用模拟二进制交叉，变异采用多项式变异，交叉率和变异率根据实际需要设置即可。
s205：合并种群：
将种群p和种群q进行合并，得到合并种群s，即s＝p∪q。显然合并种群s中包含2n个个体。
s206：计算个体的切比雪夫函数值：
分别计算2n个个体中每个个体xk到各权重向量的切比雪夫函数值，其中k＝1,2,…,2n。切比雪夫函数值的具体计算方法为：
根据传输函数计算个体所对应元件参数向量在预设测点处的输出电压相量
将每个个体xk中故障元件c的参数值xc,k进行归一化得到归一化后的参数值
其中，xc,max表示2n个个体故障元件c的参数值中的最大值。
令f1(xk)＝e(xk)，
其中，m＝1,2。
切比雪夫函数值在本发明中相当于是个体的适应度值，每个个体有n个切比雪夫函数值。
s207：个体优选：
对于每个权重向量wn，从2n个个体中选择到该权重向量的切比雪夫函数值最小的个体，作为该权重向量的最优个体，将n个权重向量对应的最优个体构成下一代种群p′。
s208：判断迭代次数t是否达到预设的最大迭代次数tmax，如果未达到，则进入步骤s209，否则进入步骤s210；
s209：令种群p＝p′，t＝t 1，返回步骤s204；
s210：确定故障参数下界：
将当前种群p′中具有最小误差的个体中故障元件c的参数值作为故障元件c的参数范围下限xcl。
s103：基于切比雪夫确定故障元件参数上限：
基于切比雪夫确定故障元件参数上限，其算法流程与步骤s102相同，但是由于此时需要确定故障元件参数上限，需要对故障件参数值的归一化计算公式进行调整，即采用以下公式将每个个体xk中故障元件c的参数值xc,k进行归一化得到归一化后的参数值
算法执行完毕后，将最后一代种群p′中具有最小误差的个体中故障元件c的参数值作为故障元件c的参数范围上限xcu。
实施例
为了更好地说明本发明的技术方案和技术效果，采用一个具体模拟电路对本发明进行实验验证。图3是本实施例中二阶托马斯模拟滤波电路的电路图。如图3所示，本实施例中二阶托马斯模拟滤波电路包括6个电阻元件，2个电容以及3个放大器，各元件参数的标称值如图3中标示。本实施例中以vout作为测点，其传输函数如下式所示：
本实施例中设置故障元件为电阻r2，令其参数值为19kω，其它元件在容差范围(本实施例设置容差参数α＝0.05，则容差范围为(xin×95％,xin×105％))内随机取值。输入信号为有效值为1v，频率为1khz的正弦信号。仿真得到故障电压相量
将元件参数向量x＝[x1,x2,…,xc]作为遗传算法种群的个体，本实施例中c＝8，前6位为电阻，最后2位为电容。设置故障元件2(即电阻r2)的参数值x2的故障取值范围[1ω,100mω]，初始种群p中每个个体中故障元件2的参数值在该故障取值范围内随机取值，其余元件i′的参数xi′在容差范围(xi′n×95％,xi′n×105％)内随机取值。设置种群数目n＝200，最大迭代次数tmax＝400。设置权重向量时所使用的范围[-λ2,-λ1]为[-16,-2]，即保证所有个体的误差均小于10-2。
图4是本实施例中基于切比雪夫确定故障元件参数下限和上限的最后一代种群示意图。如图4所示，可以确定本实施例中故障元件参数下限x2l＝17407ω，对应的误差e(x)＝3.4612e-13，接近零，精度完全满足要求。同理确定故障元件参数上限x2u＝20804ω，对应的误差e(x)＝2.0005e-11，接近零，精度完全满足要求。从而得到故障参数范围为[17407ω,20804ω]，即所有故障元件2的参数值在闭区间[17407ω,20804ω]内的故障都可以产生故障电压相量
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
1.一种基于切比雪夫的模拟电路元件故障参数识别方法，其特征在于，包括以下步骤：
s1：获取模拟电路中元件数量c和各个元件的参数标称值xin，i＝1,2,…,c，确定模拟电路在预设测点的传输函数，获取当前故障状态下模拟电路在预设测点的故障电压相量
s2：基于切比雪夫确定故障元件参数下限，具体步骤包括：
s2.1：将元件参数向量x＝[x1,x2,…,xc]作为遗传算法种群的个体，生成n个个体构成初始种群p，具体方法为：根据需要设置故障元件c的参数值xc的故障取值范围
s2.2：预设范围[-λ2,-λ1]，λ1、λ2均为正常数，且λ2＞λ1；设置n个权重向量，记第n个权重向量wn＝(wn,1,wn,2)t，其中
s2.3：初始化迭代次数t＝1；
s2.4：对种群p中的个体进行交叉和变异，得到新种群q，在交叉和变异过程中需要保证故障元件c的参数值xc在故障取值范围内取值，非故障元件的参数值在容差范围内取值；
s2.5：将种群p和种群q进行合并，得到合并种群s，即s＝p∪q；
s2.6：分别计算2n个个体中每个个体xk到各权重向量的切比雪夫函数值，其中k＝1,2,…,2n；切比雪夫函数值的具体计算方法为：
根据传输函数计算个体所对应元件参数向量在预设测点处的输出电压相量
将每个个体xk中故障元件c的参数值xc,k进行归一化得到归一化后的参数值：
其中，xc,max表示2n个个体故障元件c的参数值中的最大值；
令f1(xk)＝e(xk)，
s2.7：对于每个权重向量wk，从2n个个体中选择到该权重向量的切比雪夫函数值最小的个体，作为该权重向量的最优个体，将n个权重向量对应的最优个体构成下一代种群p′；
s2.8：断迭代次数t是否达到预设的最大迭代次数tmax，如果未达到，则进入步骤s2.9，否则进入步骤s2.10；
s2.9：令种群p＝p′，t＝t 1，返回步骤s2.4；
s2.10：将当前种群p′中具有最小误差的个体中故障元件c的参数值作为故障元件c的参数范围下限xcl；
s3：采用与步骤s2中的相同算法流程确定故障元件参数上限xcu，在算法执行过程中采用以下公式将每个个体xk中故障元件c的参数值xc,k进行归一化得到归一化后的参数值
算法执行完毕后，将最后一代种群p′中具有最小误差的个体中故障元件c的参数值作为故障元件c的参数范围上限xcu。
技术总结