基于遗传算法的模拟电路故障定位与参数辨识方法与流程

本发明属于模拟电路故障诊断
技术领域：
，更为具体地讲，涉及一种基于遗传算法的模拟电路故障定位与参数辨识方法。
背景技术：
：随着电子技术的飞速发展，各种电子设备、芯片等的体积快速减小，而故障诊断的技术尤其是电子设备和芯片中的模拟部分的故障诊断技术未能跟上发展的速度。传统的模拟电路故障诊断可以分为测前诊断(simulation-before-test，sbt)和测后诊断(simulationafterthetest，sat)。sbt即在测试前对模拟电路进行大量仿真提取模拟电路的故障特征建立故障字典，或者通过神经网络、支持向量机等分类算法对故障特征进行分类，以分类的结果定位故障元件。而sat是对测试得到的数据进行分析得到故障诊断结果，本发明的故障定位是基于遗传算法的测后诊断。现阶段的故障诊断，无论是sat还是sbt，主要是针对故障进行定位。少部分也对参数识别进行了研究，而目前的参数识别也只是给出一个具体的诊断结果，该结果在一定范围内接近元件真实参数。而模拟电路在容差影响下，不同的元件参数组合能够得到相同的输出，识别得到的该具体诊断结果仅能体现一种可能性，而不能体现出故障元件参数的可能范围。技术实现要素：本发明的目的在于克服现有技术的不足，提供一种基于遗传算法的模拟电路故障定位与参数辨识方法，通过遗传算法进行故障定位并确定故障定位后的故障元件的参数值范围。为实现上述发明目的，本发明基于遗传算法的模拟电路故障定位与参数辨识方法包括以下方法：s1：获取模拟电路在测点t处的传输函数；s2：分析得到模拟电路经测点t输出电压进行故障诊断的模糊组信息，每个模糊组选取一个元件作为代表故障元件，记代表性故障元件的数量为n，记其他非代表性故障元件的数量为m；s3：当模拟电路出现故障时，在预设的激励信号下测量得到测点t处的输出电压分别表示输出电压的实部和虚部，j为虚数单位；s4：以x＝{x1,…,xn,x′1，…,x′m}作为遗传算法中的个体，其中xn表示第n个代表性故障元件的参数值，n＝1,2,…,n，x′m表示第m个非代表性故障元件的参数值，m＝1,2,…,m，每个个体中有一个代表性故障元件的参数值随机取值，其余故障元件的参数值在其容差范围内随机取值，生成g个个体构成种群p；s5：判断是否达到故障定位遗传算法的迭代结束条件，如果是，进入步骤s10，否则进入步骤s6；s6：对当前种群进行交叉和变异操作，得到子种群q，在进行交叉和变异操作时，需要保证每个个体中参数值位于故障范围内的代表性故障元件的数量小于等于1，非代表性故障元件的参数值在容差范围内取值；s7：将种群p和种群q进行合并，构成种群s；s8：将种群s中的每个个体分别代入传输函数，得到预设的激励信号下在测点t处的输出电压ug＝αg+jβg，αg、βg分别表示输出电压的实部和虚部，g＝1,2,…,2g，然后采用以下公式计算第g个个体输出电压与当前模拟电路的输出电压之间的欧式距离dg，计算公式如下：s9：根据欧式距离dg从种群s中优选欧式距离较小的g个个体作为下一代种群p，返回步骤s5；s10：从当前种群中选择欧式距离最小的个体，该个体中参数值位于故障范围内的代表性故障元件即为故障诊断结果；s11：采用以下方法进行故障元件参数辨识：s10.1：记故障定位所得到的代表性故障元件对应的序号为在所定位个体中第个代表性故障元件的参数值为以作为遗传算法中的个体，表示第个代表性故障元件的参数值，种群中所有个体的参数值范围内均匀取值，λ是预设的[0,1]范围内的常数，其余故障元件的参数值在其容差范围内随机取值，生成r个个体构成种群w；s10.2：判断是否达到预设的最大迭代次数，如果是，进入步骤s10.7，否则进入步骤s10.3；s10.3：对种群w进行交叉和变异操作，在进行交叉和变异操作过程中第个代表性故障元件的参数值不变，其余故障元件的参数值在容差范围内取值，得到子种群v；s10.4：将种群w和种群v进行合并，构成种群h；s10.5：将种群h中的每个个体分别代入传输函数，得到预设的激励信号下在测点t处的输出电压u′r＝α′r+jβ′r，α′r、β′r分别表示输出电压的实部和虚部，r＝1,2,…,2r，然后采用以下公式计算第r个个体输出电压与当前模拟电路的输出电压之间的欧式距离d′r，计算公式如下：s10.6：将种群h中2r个个体中第个代表性故障元件参数值相同的两个个体分为一对，将每对个体中欧式距离较小的个体作为下一代种群w中的个体，返回步骤s10.2；s10.7：将种群w中的所有个体以第个代表性故障元件参数值作为横坐标，对应欧式距离作为纵坐标作图，在得到的曲线中搜索出欧式距离在[0,a]范围内的连续曲线，a为预设的正参数，将连续曲线对应的第个代表性故障元件参数值范围作为模拟电路的参数识别结果。本发明基于遗传算法的模拟电路故障定位与参数辨识方法，获取测点的传输函数，分析模拟电路的模糊组信息，确定每个模糊组的代表故障元件，当模拟电路出现故障时首先获取当前输出，再利用遗传算法确定当前输出对应的故障元件参数值向量，从而定位出发生故障的代表性故障元件，然后再基于遗传算法确定该代表性故障元件的参数值范围。采用本发明可以有效实现对模拟电路的故障诊断和故障元件参数辨识。附图说明图1是本发明基于遗传算法的模拟电路故障定位与参数辨识方法的具体实施方式流程图；图2是本发明中基于遗传算法进行故障元件参数辨识的流程图；图3是本实施例中二阶托马斯模拟滤波电路的结构图；图4是本实施例中参数辨识的曲线图；图5是本实施例中4种故障的参数辨识曲线图。具体实施方式下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。图1是本发明基于遗传算法的模拟电路故障定位与参数辨识方法的具体实施方式流程图。如图1所示，本发明基于遗传算法的模拟电路故障定位与参数辨识方法的具体步骤包括：s101：获取传输函数：获取模拟电路在测点t处的传输函数。s102：模糊组分析：分析得到模拟电路经测点t输出电压进行故障诊断的模糊组信息，每个模糊组选取一个元件作为代表故障元件，记代表性故障元件的数量为n，显然n也表示模糊组数量，记其他非代表性故障元件的数量为m。s103：确定模拟电路当前输出：当模拟电路出现故障时，在预设的激励信号下测量得到测点t处的输出电压分别表示输出电压的实部和虚部，j为虚数单位。为了使故障状态下的输出电压更加准确，可以在多次测量输出电压后进行平均，从而得到输出电压s104：初始化故障定位遗传算法种群：以x＝{x1,…,xn,x′1，…,x′m}作为遗传算法中的个体，其中xn表示第n个代表性故障元件的参数值，n＝1,2,…,n，x′m表示第m个非代表性故障元件的参数值，m＝1,2,…,m，每个个体中有一个代表性故障元件的参数值随机取值，其余故障元件的参数值在其容差范围内随机取值，生成g个个体构成种群p。s105：判断是否达到故障定位遗传算法的迭代结束条件，如果是，进入步骤s110，否则进入步骤s106。遗传算法的迭代结束条件一般有两种，一是达到最大迭代次数，一是目标函数值达到预设阈值，可以根据实际需要进行设置。s106：生成子种群：对当前种群进行交叉和变异操作，得到子种群q。在进行交叉和变异操作时，需要保证每个个体中参数值位于故障范围内的代表性故障元件的数量小于等于1，非代表性故障元件的参数值在容差范围内取值。本实施例中个体交叉采用模拟二进制交叉(sbx)和多项式变异(pol)。s107：合并种群：将种群p和种群q进行合并，构成种群s，即s＝p∪q，显然合并种群中个体数量为2g。s108：计算个体适应度值：接下来需要对种群s中的每个个体分别计算个体适应度值，对于本发明而言，是采用每个个体在预设的激励信号下得到的输出电压与当前模拟电路的输出电压之间的欧式距离来作为适应度值的，因此具体计算方法如下：将种群s中的每个个体分别代入传输函数，得到预设的激励信号下在测点t处的输出电压ug＝αg+jβg，αg、βg分别表示输出电压的实部和虚部，g＝1,2,…,2g，然后采用以下公式计算第g个个体输出电压与当前模拟电路的输出电压之间的欧式距离dg，计算公式如下：显然就故障诊断而言，应当是距离越小，表示输出电压与当前模拟电路和输出电压越接近，个体越优。s109：生成下一代种群：根据欧式距离dg从种群s中优选欧式距离较小的g个个体作为下一代种群p，返回步骤s105。本实施例中采用锦标赛选择，即打乱种群s中个体的顺序，两两比较适应度值，保留欧式距离较小的个体，淘汰欧式距离较大的个体，从而得到下一代种群。s110：得到故障定位结果：从当前种群中选择欧式距离最小的个体，该个体中参数值位于故障范围内的代表性故障元件即为故障诊断结果。s111：基于遗传算法进行故障元件参数辨识：接下来需要对代表性故障元件的参数进行辨识，得到当前的故障参数值。图2是本发明中基于遗传算法进行故障元件参数辨识的流程图。如图2所示，本发明中基于遗传算法进行故障元件参数辨识的具体步骤包括：s201：初始化参数辨识遗传算法种群：记故障定位所得到的代表性故障元件对应的序号为在所定位个体中第个代表性故障元件的参数值为以作为遗传算法中的个体，表示第个代表性故障元件的参数值，种群中所有个体的参数值在范围内均匀取值，λ是预设的[0,1]范围内的常数，其余故障元件的参数值在其容差范围内随机取值，生成r个个体构成种群w。s202：判断是否达到预设的最大迭代次数，如果是，进入步骤s207，否则进入步骤s203。s203：生成子种群：对种群w进行交叉和变异操作，在进行交叉和变异操作过程中第个代表性故障元件的参数值不变，其余故障元件的参数值在容差范围内取值，得到子种群v。本实施例中个体交叉采用仿二进制交叉(sbx)，变异操作为多项式变异(pol)。s204：合并种群：将种群w和种群v进行合并，构成种群h，即h＝w∪v，显然合并种群中个体数量为2r，种群中有r对个体中第个代表性故障元件的参数值相同。s205：计算个体适应度值：接下来需要对种群h中的每个个体分别计算适应度值，其计算方法与故障定位时相同，具体方法为：将种群h中的每个个体分别代入传输函数，得到预设的激励信号下在测点t处的输出电压u′r＝α′r+jβ′r，α′r、β′r分别表示输出电压的实部和虚部，r＝1,2,…,2r，然后采用以下公式计算第r个个体输出电压与当前模拟电路的输出电压之间的欧式距离d′r，计算公式如下：s206：生成下一代种群：将种群h中2r个个体中第个代表性故障元件参数值相同的两个个体分为一对，将每对个体中欧式距离较小的个体作为下一代种群w中的个体，返回步骤s202。s207：确定模拟电路参数识别结果：将种群w中的所有个体以第个代表性故障元件参数值作为横坐标，对应欧式距离作为纵坐标作图，在得到的曲线中搜索出欧式距离在[0,a]范围内的连续曲线，a为预设的正数，将连续曲线对应的第个代表性故障元件参数值范围作为模拟电路的参数识别结果，该参数值范围内的第个代表性故障元件参数值均可以达到当前模拟电路的实际电压输出。实施例为了更好地说明本发明的技术效果，采用二阶托马斯模拟滤波电路为例对本发明进行说明。图3是本实施例中二阶托马斯模拟滤波电路的结构图。如图3所示，本实施例中二阶托马斯模拟滤波电路以vout作为测点，其传输函数为：本实施例采用
专利名称：为“一种模拟电路模糊组识别方法”、专利号为“20141033627.7”的专利中的方法，基于圆模型进行模糊组分析，得到该测点下的模糊组情况为：{r1}、{r2}、{r3,c1}、{r4,r5,r6,c2}，模糊组内部元件的故障不可区分，模糊组之间的故障理论上都能被区分。因此取r1、r2、r3、r4故障作为故障类型，其余元件为非故障类型元件。首先故障定位遗传算法种群初始化数目为50，其中每种故障类型初始化10个个体，最大迭代次数为100，交叉概率为1，变异概率为0.1，限定r5、r6、c1、c2非故障类型元件始终在容差范围内，各元件容差为5％，且r1、r2、r3、r4始终只允许随机一个元件参数超出容差范围。模拟电路输入激励频率为1khz，模拟电路各元件参数为x＝[r1,r2,r3,r4,r5,r6,c1,c2]。表1是本实施例中托马斯电路元件标称值和故障值。表1以元件参数x＝[13000ω,9800ω,9600ω,10140ω,9850ω,10400ω,9.7960nf,9.9650nf]为例，得到测点vout处输出电压为[-0.7361,0.2339j]。表2是本实施例的故障定位结果表。故障类型诊断结果欧式距离正确率r1r18.3210e-06100％表2得到的故障定位最优个体为：x＝[13354ω,10181ω,9598ω,10500ω,9542ω,10021ω,10.158nf,9.6804nf]接下来进行参数辨识，初始化种群数1600，设置参数λ＝50％，种群个体中r1的参数值在13354±(50％)内均匀取1600个数，其余故障元件参数均始终限定在容差范围±5％内。交叉概率为1，变异概率为0.1，最大迭代次数为800。图4是本实施例中参数辨识的曲线图。表3是本实施例中参数辨识结果。故障类型及故障值参数识别结果(区间)r1(13kω)[12.37kω,13.93kω]表3为了进一步测试本发明的方法，对二阶托马斯电路的五种故障类型均做了测试，并对结果进行统计。表4是本实施例中4种故障分别进行50次故障定位结果统计表。故障分组r1r2r3r4r150000r204900r300500r401050诊断正确率100％98％100％100％表4图5是本实施例中4种故障的参数辨识曲线图。表5是本实施中故障辨识结果表。表5根据表4、表5和图5所示，可以看出采用本发明可以准确定位代表性故障元件，并且能够准确识别代表性故障元件参数得到其参数值范围。尽管上面对本发明说明性的具体实施方式进行了描述，以便于本
技术领域：
的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本
技术领域：
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。当前第1页12