一种IGBT性能退化特征参数的提取方法与流程

本发明属于电力电子器件与装置可靠性
技术领域：
，涉及一种igbt性能退化特征参数的提取方法，具体涉及一种基于核主成分分析和加权马氏距离的igbt性能退化特征参数的提取方法，可应用于对igbt性能退化进行在线监测。
背景技术：
：绝缘栅双极型晶体管igbt(insulatedgatebipolartransistor)是由双极型三极管和绝缘栅型场效应管组成的复合全控型电压驱动式功率半导体器件。igbt具有输入阻抗高、快速开关、电流密度大、导通压降低等优点，被广泛应用于逆变器、电机传动控制、汽车、轨道交通等民用领域，而且在航空、航天等军事领域也拥有很好的发展前景。但是igbt在面对过高电压和高温等恶劣环境时容易发生失效，如果igbt突然发生失效，轻则造成电子设备出现故障停机，重则造成严重事故。因此，需要提取表征igbt性能退化的特征参数并根据性能退化特征参数对igbt性能退化进行在线监测，实现提前更换维护、降低维护成本和提高系统稳定性。性能退化特征参数是一种表征系统退化状态的参数，能够衡量系统的健康状态，也称为特征参数或退化因子或性能参数或退化指标或退化参数。目前国内外表征igbt的性能退化状态多采用单一退化特征建立性能退化特征参数，如饱和开态电压、集电极电流、栅极电压、栅极阈值电压、开通/关断时间、结温、导通电阻、拖尾电流拟合系数。然而由于不同性能退化特征参数包含的igbt性能退化信息不同，另外性能退化特征参数对igbt性能退化机制的敏感度不同，因此单一退化特征不能全面的表征igbt的性能退化过程。有效地建立融合多退化特征的性能退化特征参数能够提高性能退化特征参数对igbt性能退化过程表征的准确性，进而提高了igbt性能退化在线监测的精度，具有重要意义。例如，申请公布号为cn106124957a，名称为“一种绝缘栅双极型晶体管退化的在线监测方法”的专利申请，公开了一种基于结温优化的igbt退化参数提取方法。该方法首先检测igbt在正常情况下表征结温水平的不变压降和表征性能退化程度的导通电阻，并计算相同结温下导通电阻变化量与结温变化量的比例系数；其次在线检测igbt在性能退化情况下的不变压降和导通电阻，结合正常情况下得到的比例系数对导通电阻进行去结温影响的优化。该方法考虑了结温对导通电阻的影响，得到了更能代表igbt性能退化状态的性能退化特征参数，实现了对igbt性能退化的在线监测。其存在的不足之处是，首先该方法采用的igbt性能退化特征参数只有一种，不能全面地包含igbt性能退化信息，不利于实现高精度的igbt性能退化的在线监测。技术实现要素：本发明的目的在于克服上述现有技术存在的缺陷，提供了一种igbt性能退化特征参数的提取方法，旨在提高对igbt性能退化在线监测的精度。为了实现上述目的，本发明采用的技术方案包括如下步骤：(1)采集igbt性能退化的检测数据：采集igbt在每一个开通/关断周期的开通状态集电极-发射极的饱和开态电压vce-oni、关断状态t时间段内集电极的拖尾电流和栅极的漏电流并将m个开通/关断周期内所采集的检测数据组成检测数据矩阵其中i＝1,2…m，m≥2，ti表示第i个开通/关断周期采集数据的时间点；(2)计算的拟合系数和的拟合系数：(2a)构建指数模型f1i(t)＝exp(p1(i)t3+p2(i)t2+p3(i)t+p4(i))和指数模型其中p1(i)、p2(i)、p3(i)和p4(i)分别为f1i(t)的系数，g1(i)、g2(i)、g3(i)和g4(i)为的系数，exp(·)为以自然常数e为底的指数；(2b)通过f1i(t)对进行拟合，得到p1(i)、p2(i)、p3(i)和p4(i)的值，并将p1(i)作为拖尾电流的拟合系数，同时通过对进行拟合，得到g1(i)、g2(i)、g3(i)和g4(i)的值，并将g4(i)作为漏电流的拟合系数，然后对p1(i)和g4(i)进行备份；(3)构建igbt的特征矩阵xm×3：对检测数据矩阵sm×3内m个开通/关断周期的饱和开态电压，以及拖尾电流的拟合系数的备份数据和栅极的漏电流的拟合系数的备份数据进行合并，得到igbt的特征矩阵xm×3：(4)对特征矩阵xm×3进行核主成分分析：采用核主成分分析算法对特征矩阵xm×3进行核主成分分析，得到核主成分矩阵zm×p及贡献率矩阵l1×p为：zm×p＝[z1,…,zj,…,zp]l1×p＝[l1,…,lj,…,lp]其中，为第j个核主成分，lj为第j个核主成分的贡献率，j＝1,2…p，p为核主成分个数；(5)构建样本集和健康样本集；构建样本集同时构建健康样本集其中为核主成分矩阵zm×p的第i行元素构成的样本，为第j个核主成分的前n个元素构成的健康样本，n≥1；(6)获取igbt的性能退化特征参数hi：计算每个核主成分zj的马氏距离加权系数得到加权矩阵w＝diag(w1,…,wj,…,wp)，并采用加权马氏距离公式，通过w计算与z*的偏离度di(i)，并将m个开通/关断周期的偏离度hi＝[di(1)，…,di(i),…,di(m)]作为igbt的性能退化特征参数，其中diag(w1,…,wj,…,wp)为以w1,…,wj,…,wp为对角元素的对角阵。本发明与现有技术相比，具有如下优点：1.本发明采用核主成分分析对三种退化特征进行核主成分分析，剔除了三种退化特征间的冗余信息，再采用核主成分分析得到的贡献率作为加权马氏距离的权重参考，通过加权马氏距离得到了igbt性能退化特征参数，合理地区分了核主成分之间的重要性区别，提高了igbt性能退化特征参数准确性，有效提高了igbt性能退化在线监测的精度。2.本发明使用指数模型从栅极漏电流和拖尾电流中提取出包含igbt性能退化信息的漏电流拟合系数和拖尾电流拟合系数，并在漏电流拟合系数、拖尾电流拟合系数和饱和开态电压三种退化特征分别代表igbt栅极、集电极和集电极-发射极输出的退化信息，相比于单一退化特征，从三种退化特征中提取的性能退化特征参数更全面地包含igbt性能退化信息，进一步提高了igbt性能退化在线监测的精度。附图说明图1为本发明的实现流程图；图2为本发明集电极拖尾电流拟合系数p1提取结果图；图3为本发明栅极漏电流拟合系数g4提取结果图；图4为本发明开态饱和电压vce-on提取结果图；图5为本发明性能退化特征参数hi提取结果图。具体实施方式下面结合附图和具体实施例，对本发明作进一步详细描述：参照图1，本发明包括如下步骤：步骤1)采集igbt器件性能退化的检测数据：保持壳温在(268℃,270℃)的区间，在型号为irg4bc30k的igbt的栅极施加方波信号，实现对igbt寿命的加速，直到igbt出现闩锁效应时截止。在对igbt进行寿命加速期间，采集igbt在每一个开通/关断周期的开通状态集电极-发射极的饱和开态电压vce-oni、关断状态t时间段内集电极的拖尾电流和栅极的漏电流并将m个开通/关断周期内所采集的检测数据组成检测数据矩阵其中i＝1,2…m，m≥2，ti表示第i个开通/关断周期采集数据的时间，本实施例m＝710；步骤2)计算的拟合系数和的拟合系数：(2a)构建指数模型f1i(t)＝exp(p1(i)t3+p2(i)t2+p3(i)t+p4(i))和指数模型其中p1(i)、p2(i)、p3(i)和p4(i)分别为f1i(t)的系数，g1(i)、g2(i)、g3(i)和g4(i)为的系数，exp(·)为以自然常数e为底的指数；(2b)通过f1i(t)对进行拟合，得到p1(i)、p2(i)、p3(i)和p4(i)的值，并将p1(i)作为拖尾电流的拟合系数，同时通过对进行拟合，得到g1(i)、g2(i)、g3(i)和g4(i)的值，并将g4(i)作为漏电流的拟合系数，然后对p1(i)和g4(i)进行备份，得到如图2所示的拖尾电流的拟合系数集合从图2中可以看出p1在igbt性能退化初期平缓退化，在igbt性能退化后期快速退化，符合igbt性能退化特点，如图3所示的漏电流的拟合系数集合从图3中可以看出g4在igbt性能退化整个阶段都为快速退化。步骤3)构建igbt的特征矩阵x710×3：对检测数据矩阵s710×3内710个开通/关断周期的饱和开态电压，以及拖尾电流的拟合系数的备份数据和栅极的漏电流的拟合系数的备份数据进行合并，得到igbt器件的特征矩阵x710×3：其中饱和开态电压集合如图4所示，从图4中可以看出g4在igbt性能退化初期没有明显退化趋势，在igbt性能退化后期快速退化；步骤4)对特征矩阵x710×3进行核主成分分析：(4a)为了消除量纲的影响，对特征矩阵x710×3进行标准化，得到标准特征矩阵p1、g4和vce-on之间存在很多冗余信息，如果直接将性能退化特征参数进行融合会降低性能退化特征参数的准确性，而核主成分分析能够很好的剔除冗余信息并且能够适用于非线性数据，因此采用核主成分分析对特征矩阵x710×3提取核主成分。采用高斯径向基核函数将x*710×3映射到710维特征空间，得到大小为710×710的核矩阵k，其中为特征矩阵x710×3的第j列元素的平均值；为第j列元素的标准差，i＝1,2,…710，j＝1,2,3；(4b)计算核矩阵k的特征值λ1,…,λd,…,λ710及其对应的特征向量v1,…,vd,…,v710，并按照降序的方式分别对710个特征值和710个特征向量进行排序，得到特征值λ′1,…,λ′d,…,λ′710和特征向量v′1,…,v′d,…,v′710；(4c)设提取效率为并计算特征向量v′1,…,v′d,…,v′710的贡献率l1,…,ld,…,l710及累积贡献率b1,…,bd,…,b710，如表1所示；表1贡献率及累计贡献率结果特征向量v′d贡献率ld/％累计贡献率bd/％172.61172.611213.61886.229310.73896.96741.76098.72750.67499.40160.23099.63170.16899.799(4d)按照从前到后的顺序逐个对b1,…,bd,…,b710与进行比较，得到当d＝2时bd大于或等于时，令p＝2，得到特征向量矩阵v＝[v′1,v′2]及贡献率矩阵l1×p＝[l1,l2]；(4e)计算核矩阵k在特征向量矩阵v上的投影，得到核主成分矩阵z710×2＝k·v。步骤5)构建样本集和健康样本集；构建样本集同时构建健康样本集其中为核主成分矩阵zm×p的第i行元素构成的样本，为第j个核主成分的前100个元素构成的健康样本；(6)获取igbt的性能退化特征参数hi：加权马氏距离相对于马氏距离考虑到了变量之间不同的重要性，以重要性为加权系数的参考，是一种有效的距离融合方法。使用贡献率来评价核主成分的重要性，计算每个核主成分zj的马氏距离加权系数得到加权矩阵并采用加权马氏距离算法，通过w计算与z*的偏离度并将710个开通/关断周期的偏离度hi＝[di(1)，…,di(i),…,di(710)]作为igbt的性能退化特征参数，性能退化特征参数的提取结果如图5所示，从图中可以看出hi在igbt性能退化初期，为缓慢退化状态，在igbt性能退化后期，为快速退化状态。计算hi、vce-on、g4、p1及开通/关断周期t的相关系数矩阵：从相关矩阵中可以看出，hi、vce-on、g4、p1都与t具有显著的相关性，证明hi、vce-on、g4、p1都包含了igbt的健康状态信息，另外hi与vce-on、g4、p1都具有显著相关性，证明了hi有效保留了vce-on、g4和p1的退化状态信息，能够有效地代表igbt的退化状态。另外，在igbt器件性能退化初期时，g4处于快速退化状态，p1和vce-on处于缓慢退化状态，此时hi退化趋势相对平缓，而在igbt性能退化中后期，vce-on、g4和p1都处于快速退化状态，hi也处于快速退化状态。由此也可以看出经过一系列的数据处理，本发明提取的性能退化特征参数相对于单一性能退化特征参数，更全面地保留了igbt的性能退化信息，能够有效地表征igbt的性能退化过程，提高了igbt性能退化在线监测的精确度。综上所述，本发明能有效地提取igbt器件的性能退化特征参数，提高了对igbt性能退化在线监测的精确度。当前第1页12