一种基于主元分析与多分类相关向量机的多电平逆变器故障诊断策略的制作方法
【专利摘要】本发明公开了一种基于主元分析与多分类相关向量机(PCA-mRVM)的多电平逆变器故障诊断策略。该诊断策略首先通过主元分析将原始样本进行降维，提取出多个带有故障特征的主成分，以便形成训练样本；其次通过多分类相关向量机对训练样本进行故障诊断，并输出各故障类别的概率，以概率最大的故障类别作为诊断结果。本发明的故障诊断策略用于具有较大样本空间以及较多类别的故障诊断中有明显的优势，具有模型稀疏度高，计算复杂度低等优点，更重要的是mRVM可以输出各类别成员概率，输出具有概率统计意义，便于分析不确定性问题。
【专利说明】—种基于主元分析与多分类相关向量机的多电平逆变器故障诊断策略
【技术领域】:
[0001]本发明涉及电力电子领域中的多电平逆变器的故障诊断，尤其涉及一种基于PCA-mRVM的故障诊断方法。
【背景技术】:
[0002]随着电力电子技术的发展和电力电子器件生产成本的降低，高压大功率变换器被广泛应用于各种电气设备中，例如大功率交流电机传动、有源电力滤波和新能源并网等。而为了满足日益提高的电力系统发展需求，多电平逆变器作为一种新型的变换器，因其主电路结构简单以及控制电路非常灵活等优点，作为电力电子领域中高压大功率变换方面最活跃的一个分支，在工业生产中得到了迅速的推广。其中级联H桥型逆变器因其不需要大量箝位二极管或电容，无需均衡电容电压，结构上易于模块化和扩展，工作效率高，因此得到了较为广泛的应用。虽然级联H桥型逆变器是多电平逆变器中输出同样数量电平而所需器件最少的一种，但它所需要的开关器件数量还是随着系统电平数的增加而成倍增加。大量使用功率半导体器件增加了故障发生的可能性(开关开路或短路)，而且随着电压的升高，故障出现的概率也变大。级联H桥型多电平逆变器的产生虽然为电力电子技术在高压、大功率场合的应用提供了很多便利，但一旦发生故障，轻则造成企业停产，重则会造成灾难性事故，给社会带来巨大的损失。研究表明，逆变器供电的变频调速系统中开关器件的故障占整个驱动系统故障的82.5%，是驱动系统中最容易发生故障的环节。
[0003]目前多电平逆变器开路故障诊断方法主要有=Rothenhagen提出了IGBT开路的故障诊断方法，该方法需要大量的检测点和传感器，这样会使电路变得更加复杂，成本更高，且容易受到外部干扰而导致误诊断；汤清泉等人提出了根据检测逆变器输出侧PWM电压和输出电流极性来诊断功率管开路故障的故障识别方案，该方法对逆变器的工作模式进行分区，不仅可用于单个功率管发生故障，还可诊断多个功率管发生故障的情况，但由于输出电流受负载影响，导致该方法泛化能力较差，且随着逆变器电平数的增多，对工作模式进行分区也会变得极为复杂；另外杨忠林等人提出了通过检测逆变器直流侧电流的频率成分就可以实现逆变器开关管的故障诊断，当各开关管工作正常时，直流侧电流中没有谐波成分，而当有开关管发生故障时，直流侧中则会出现谐波成分。虽然该方法可以对开关管故障进行检测，但却无法具体定位到故障开关管的位置，且当电路较复杂时，容易受到外部干扰而导致误诊断。由以上分析可知，目前针对多电平逆变器的故障诊断问题，主要存在以下困难:(I)故障特征样本较难提取；(2)样本数据维数较大，处理较难且计算量大；(3)目前的诊断方法仅能输出诊断精度，不能输出各个类别的具体结果，不便分析不确定性问题等。
【发明内容】
:
[0004]本发明所要解决的技术问题在于针对上述现有技术的技术问题而提出的一种适用于级联H桥多电平逆变器故障诊断策略，其目的就是为了克服已有技术中存在的泛化能力差、成本高、诊断速度慢等缺点，通过利用FFT将样本数据从时域转换到频域，从而便于故障特征样本的提取，其次通过PCA进行降维，减少计算量等，最后采用mRVM，可以输出各类别成员概率，输出具有概率统计意义，便于分析不确定性问题。
[0005]为了实现上述目的，本发明采用了如下技术方案:
[0006]1、构造一个多电平逆变器仿真模型
[0007]采用Matlab/Simulink搭建一个单相级联H桥五电平逆变电路系统模型，然后加入故障环节，依次进行不同位置的单个开关管开路故障测试，提取一个周期的故障特征信号，建立故障特征集；
[0008]2、数据预处理
[0009]利用FFT将信号转换到故障特征较为明显的频域并提取故障特征信号；
[0010]3、提取故障主元
[0011]由于故障特征样本维数较大，直接处理的计算复杂度较高，且处理速度较慢，因此此处采用PCA对样本进行降维，提取出对分类有影响的主成分，从而显著提高处理速度；
[0012]4、模型训练阶段
[0013]初始化mRVM参数，本发明采用RBF核函数作为mRVM的核函数，并使用交叉验证法来确定核参。然后将简化后的故障特征样本代入mRVM模型进行训练；训练结束后则用测试样本对该诊断策略进行测试，若训练得到的分类精度达到设定的精度要求，则训练完成；若训练结果未达到既定要求，则跳转至特征提取步骤并更改mRVM的核参。训练的精度要求可以根据工作条件自行设定。
[0014]5、实时诊断阶段
[0015]在该阶段，直接将待分类的样本数据导入训练好的模型中，即可得到各个样本的概率输出结果，并将其中概率最大的样本作为分类诊断结果。
[0016]本发明的技术效果在于:本发明所采用的PCA-mRVM策略，首先，在传统RVM分类诊断算法的基础上加入了降维环节，以满足具有较大样本的情况，且本发明所采用的PCA算法，所提取的特征样本更加准确且更具有代表性；其次，该诊断策略可以直接输出各类别成员概率，输出具有概率统计意义，便于分析不确定性问题。由于该诊断策略仅采用输出端电压信号作为其输入信号，因此避免了传感器的大量使用，并降低了系统的复杂度以及构建成本。
【专利附图】
【附图说明】
[0017]图1为本发明多电平逆变器实时故障诊断策略结构图，其中:
[0018]1:三相异步电机
[0019]2:级联H桥多电平逆变器
[0020]3:电压采样电路
[0021]图2为本发明故障诊断策略的模型训练过程流程图；
[0022]图3为本发明中逆变器单个开关管故障的输出电压波形图；
[0023]图4为对输出电压波形进行FFT之后的频谱图；
[0024]图5为分层贝叶斯模型结构图；
[0025]图6为Fast Type-1I ML的程序流程图。【具体实施方式】:
[0026]为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体图示，进一步阐释本发明。
[0027]—种多电平逆变器实时故障诊断策略结构图如图1所示。诊断策略主要由四部分组成:数据预处理以及特征提取、分类诊断模型、输出诊断结果、开关模式计算。本发明主要研究前三部分，开关模式计算这里暂不考虑。该策略工作原理如下:首先对逆变器的输出电压信号进行预处理并对其特征提取，从而降低样本的维数；其次，将处理后的数据带入已训练好的mRVM模型(该模型应用于陌生样本时只需在初始时训练一次即可)进行诊断分类；然后输出各类别成员概率，并将其中概率最大的样本作为分类诊断结果。
[0028]Surin Khomifoi等人分析了级联多电平其输出电压的波形是和开关的通断紧密相关的，以单相级联H桥为例，如果发生一个短路故障，某一个开关输出电压将一直为零。而如果发生一个开路故障，输出电压为(一个周期内均值)。当开关发生故障(开路或短路)时不能产生期望的输出电压波形。分析可知，不同H桥单元发生故障时，其故障特征不同，输出相电压信号与故障类型和位置是紧密相关的，且输出相电压不随负载变化。因此本发明采用输出电压作为故障诊断信号。
[0029]下面主要介绍本发明提出的PCA-mRVM多电平逆变器诊断策略，主要包括模型构建、数据预处理、提取故障主元、模型训练阶段以及测试阶段等。图2为所提出故障诊断策略的训练过程流程图，具体的步骤如下。
[0030]Stepl:模型构建
[0031]采用Matlab/Simulink搭建一个单相级联H桥五电平逆变电路系统模型，利用SHPWM (Subharmonic PWM)进行控制，调制波频率为50Hz，载波频率为4kHz，采样频率为40kHz。然后加入故障环节，依次进行不同位置的单个开关管开路故障测试，提取逆变器每种故障一个周期的输出电压信号，如图3所示，建立故障特征集；
[0032]St印2:数据预处理
[0033]采样得到的一个周期的离散输出电压信号
【权利要求】
1.一种基于主元分析与多分类相关向量机的多电平逆变器故障诊断策略，其特征在于，依次含有以下步骤: 步骤1:构造一个多电平逆变器仿真模型 采用Matlab/Simulink搭建一个单相级联H桥五电平逆变电路系统模型，利用SHPWM进行控制，调制波频率为50Hz，载波频率为4kHz，采样频率为40kHz，然后加入故障环节，依次进行不同位置的单个开关管开路故障测试，提取逆变器每种故障一个周期的输出电压信号，建立故障特征集； 步骤2:数据预处理 步骤I采样得到的一个周期的离散的所述输出电压信号，经过FFT变换后，得到数据集Xaxb，其中a为观测值个数，b为谐波个数(对应O到b-Ι次谐波)，FFT变换公式:Fk =Gk+ W1)Hk' H-WH，这里 Wb=e_j2〃b，k=0，1，...，b/2-l，将预处理后的故障样本数据Xaxb与其分类标签向量I 一一对应，作为输入样本数据[X，I]； 步骤3:提取故障主元 (1)由步骤2得预处理后的数据样本Xe Raxb,其中a为对应数据样本点数，b为对应样本维数； (2)计算协方差矩阵，求解协方差矩阵Rx的特征值λ与相对应的特征值向量P:
协方差矩阵=Rx=E {[X-E (X) ] [X-E (X) ]τ}； 通过求解方程I入1-Rx |=0和I λ J-RxIpi=O, i=l, 2，…，b求得特征值λ和特征值向量P，其中，XiS Rx的第i个特征值，满足A1S入^仏是相应于特征值Xi的特征向量，P=[pl,p2,…，pb]T ； (3)根据累计贡献率选择相对主元个数m:
【文档编号】G06F17/50GK103761372SQ201410005022
【公开日】2014年4月30日 申请日期:2014年1月6日 优先权日:2014年1月6日
【发明者】王天真, 徐浩, 吴昊, 张健, 钱坤 申请人:上海海事大学