CN105548792B - 基于预测控制的矩阵变换器开关开路故障诊断方法 - Google Patents

基于预测控制的矩阵变换器开关开路故障诊断方法 Download PDF

Info

Publication number
CN105548792B
CN105548792B CN201510995957.9A CN201510995957A CN105548792B CN 105548792 B CN105548792 B CN 105548792B CN 201510995957 A CN201510995957 A CN 201510995957A CN 105548792 B CN105548792 B CN 105548792B
Authority
CN
China
Prior art keywords
mrow
msub
mtd
mtr
msubsup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510995957.9A
Other languages
English (en)
Other versions
CN105548792A (zh
Inventor
彭涛
但汉兵
邓慧
孙尧
粟梅
杨建�
韩华
朱奇
王春生
桂卫华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN201510995957.9A priority Critical patent/CN105548792B/zh
Publication of CN105548792A publication Critical patent/CN105548792A/zh
Priority to PCT/CN2016/105673 priority patent/WO2017114001A1/zh
Application granted granted Critical
Publication of CN105548792B publication Critical patent/CN105548792B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • G01R31/54—Testing for continuity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/327—Testing of circuit interrupters, switches or circuit-breakers
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2117/00—Details relating to the type or aim of the circuit design
    • G06F2117/02—Fault tolerance, e.g. for transient fault suppression

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Ac-Ac Conversion (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种基于预测控制的矩阵变换器开关开路故障诊断方法，包括步骤：在矩阵变换器正常工作的情况下，建立输出电压与输入电压、输入电流与负载电流之间的关系模型；确定矩阵变换器的所有开关组合状态的数目；建立状态空间模型，获得在下一采样周期负载电流、输入电流与输入电压的预测值；确定评价函数；采用有限集模型预测控制策略，在每一采样周期，在所有开关组合状态中选择使评价函数值最小的开关状态，作为下一采样周期的开关状态；建立状态空间模型，得到故障相负载电流在开关开路故障后的变化规律；对负载电流进行在线监测，根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置。本发明能及时诊断开路故障。

Description

基于预测控制的矩阵变换器开关开路故障诊断方法
技术领域
本发明涉及电路故障诊断领域，尤其涉及基于预测控制的矩阵变换器开关开路故障诊断方法。
背景技术
近年来，随着电力电子技术的迅猛发展，出现了许多新型的电源转换器，在诸多交-交变换器中，矩阵变换器可以实现所有交流参数(相数、相位、振幅、频率)的变换，并以其特有的优势吸引了大批学者的关注。与传统的变换器相比，矩阵变换器具有以下的优点：不需要中间直流储能环节、能够四象限运行、电源设计紧凑、具有优良的输入电流波形和输出电压波形、能量可双向流动、体积小、重量轻以及可自由控制输入功率因数等。另外，由于矩阵变换器没有电解电容器，所以寿命较长，可以广泛地应用于航空航天、军事和医疗等对体积、重量有较高要求的领域。
然而，矩阵变换器***的开关器件长期处于高频调制下，发热较为严重，开关损耗较大，也极易发生故障。当某个开关器件出现开路故障时，如不能及时发现并进行处理，很容易损坏整个矩阵变换器***。特别是长时间运行后，会严重影响***的运行性能和设备使用寿命。因此，对矩阵变换器***进行故障诊断与容错至关重要。
预测控制近年来被应用于电力电子变换器，因为该方法具有诸多优点，比如动态响应快速、思路清晰和约束简单等，可灵活地应用于诸多***。然而目前尚未见到用预测控制的思想进行矩阵变换器故障诊断与容错的对策。因此，需设计一种基于预测控制的矩阵变换器开关开路故障诊断方法，及时诊断***故障的发生并采取相应的容错措施，以提高矩阵变换器***的可靠性及使用寿命。
发明内容
本发明目的在于提供一种基于预测控制的矩阵变换器开关开路故障诊断方法，以解决矩阵变换器的开关器件出现开路故障时无法及时发现的技术问题。
为实现上述目的，本发明提供了一种基于预测控制的矩阵变换器开关开路故障诊断方法，包括以下步骤：
S1：在矩阵变换器正常工作的情况下，依据采样周期采集矩阵变换器的输出电压、输入电压、输入电流与负载电流，建立输出电压与输入电压、输入电流与负载电流之间的关系模型；
S2：在满足约束条件的情况下，确定矩阵变换器的所有开关组合状态的数目；
S3：建立矩阵变换器负载端状态空间模型，获得负载电流在下一采样周期的预测值；
S4：建立矩阵变换器输入端滤波器状态空间模型，离散化后，获得下一采样周期输入电流与输入电压的预测值；
S5：确定评价函数；
S6：采用有限集模型预测控制策略，在每一采样周期，在所有开关组合状态中选择使评价函数值最小的开关状态，作为下一采样周期的开关状态；
S7：建立开关开路故障条件下的等效电路的状态空间模型，得到故障相负载电流在开关开路故障后的变化规律；
S8：对负载电流进行在线监测，当负载电流持续为零时，根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置。
作为本发明的进一步改进：
优选地，步骤S8完成后，方法还包括：
S9：在识别故障开关的位置后，将所有开关组合状态中与故障开关的动作有关的开关状态去除，确定余下所有开关组合状态的数目；
S10：在余下的所有开关组合状态中，重新确定具体故障下的评价函数；
S11：采用有限集模型预测控制策略，在每一采样周期，在余下的所有开关组合状态中选择使评价函数值最小的开关组合状态，作为下一采样周期的开关状态。
优选地，输出电压与输入电压、输入电流与负载电流之间的关系模型如下：
其中，uoA、uoB和uoC分别为A、B和C三相的输出电压；uea,、ueb和uec分别为A、B和C三相的输入电压；iea、ieb和iec分别为A、B和C三相的输入电流；ioA、ioB和ioC分别为A、B和C三相的负载电流；S为开关状态矩阵；S中的元素SXY(X∈{A,B,C},Y∈{a,b,c})表示该开关的开关状态，SXY＝1表示开关SXY为闭合状态，SXY＝0表示开关SXY为断开状态。
优选地，约束条件包括：
S201：输入端不短路，输出端不断路；
S202：在任意时刻，连接到同一相输出的三个双向开关中，有且仅有一个开关导通，另外两个开关关断；
S203：满足各个开关占空比的非负性；
在上述约束条件下，开关状态矩阵S中的每个元素满足：
可得矩阵变换器正常工作时的所有开关组合状态有3*3*3＝27种。
优选地，步骤S3，包括以下步骤：
S301：建立矩阵变换器负载端状态空间模型如下：
其中，io和uo分别为负载电流与输出电压；L和R分别表示负载的电感和电阻；
S302：定义一个开关周期Ts，根据向前欧拉公式得到负载电流在下一采样周期的预测值io k+1为：
其中，io k和uo k分别为当前采样时刻的负载电流值和输出电压值；系数F1和F2分别为：
优选地，步骤S4，包括以下步骤：
S401：建立矩阵变换器输入端滤波器状态空间模型如下：
其中us、ue、is和ie分别表示电源电压、滤波电容电压、电源电流和输入电流；Ri、Li和Ci分别表示滤波电阻、电源和滤波电感以及滤波电容；
S402：结合(9)式和(10)式得到：
其中：
S403：将(11)式离散化，得到：
式中和分别表示在下一采样周期的滤波电容电压和电源电流值；和分别表示在当前采样时刻的电源电压、滤波电容电压、电源电流和输入电流值；且有：
S404：将(15)式代入(14)式，得到下一采样周期的输入电流以及滤波电容电压：
其中，E2、E3、D2和D3均为矩阵G的元素，E1、E4、D1和D4均为矩阵H的元素，由Ri、Li、Ci和采样周期Ts决定。
优选地，步骤S5，包括以下步骤：
S501：确定与矩阵变换器的输入电流相关的评价函数；
矩阵变换器输入电流的参考值为：
其中φ是输入电压的相角，为参考输入电流的幅值，并且有：
其中Usm为电源电压us的幅值；η为矩阵变换器的效率，由如下公式获得：
其中的Po和Pin分别表示输出功率和输入功率，分别由如下公式得到：
其中I* om表示参考负载电流i* o的幅值；期望输入电流is功率因数为1，所以式中θ为0；
S502：确定预测控制中关于输入电流is期望达到单位输入功率因数性能的评价函数g1为：
式中，△isα和△isβ分别表示两相静止坐标系下的输入电流在下一采样周期预测值is p与参考值is *偏差的α和β分量，由如下公式得到：
式中，Tabctoαβ是三相静止坐标系转换为两相旋转坐标系的变换矩阵，由下式得到：
S503：确定预测控制中与矩阵变换器负载电流相关的评价函数：
定义关于负载电流预测值期望达到跟踪参考负载电流性能的评价函数g2为：
式中，△ioα和△ioβ分别表示两相静止坐标系下的负载电流在下一采样周期预测值io p与参考值io *偏差的α和β分量，由如下公式得到：
S504：定义总评价函数为：
g＝λ*g1+g2 (28)
其中，λ是权重因子，表示输入电流的性能指标在整个矩阵变换器***的控制中的重要程度；
则，步骤S6，具体为：在矩阵变换器运行过程中，采用有限集模型预测控制策略，即在每个采样周期计算所有开关组合状态的总评价函数的g值，选择其中能使g值最小的一种开关状态，作为下一采样周期的开关状态，实现达到***期望的动态性能和稳态性能指标的运行。
优选地，步骤S7，包括以下步骤：
S701：假设如下条件S7011至S7013成立：
S7011：采样周期为Ts；
S7012：将输入电压周期分为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ和Ⅵ顺序排列的六个扇区，输入电压相角在扇区Ⅴ并且即将进入扇区Ⅵ；A输出相的参考电流值i* oA和实际值ioA都是正值；
S7013：SAa发生开路故障，在kth个采样周期Ts k的时间内，预测控制器为了缩小A输出相电流的实际值与参考值的偏差，始终给SAa导通信号；
则有，A输出相开路，A相负载电流ioA会流经箝位电路，其电路状态空间模型描述如下：
式中，uN是负载中性点与电源接地点之间的电压；
输出相的电压表示如下：
其中，Ucp是箝位电路的电容电压，umax为输入三相电压的最大值；
由于A相负载电流流过箝位电路，所以Ucp远大于umax，表示为：
其中，umin表示输入三相电压的最小值；
根据(29)式和(30)式，得到A相负载电流ioA的变化率为：
由(31)可得：
将(33)代入(32)，得到：
当ioA>0时 (34)
S703：构造能量函数v：
其变化率表示为：
于是，得到：
当ioA<0时 (37)
由(34)式和(37)式可知：当故障开关满足SAa＝1时，A相负载电流的幅值会越来越小，从而使其实际值与参考值之间的偏差越来越大；因为SAa在扇区VI和I连通A输出相对应的最大输入电压，预测控制器始终给开关SAa导通信号，然而，SAa因为故障无法导通，结果就是A相负载电流幅值减小并且连续几个采样周期内SAa＝1；在一定时间内，A相负载电流的幅值逐渐减小到零；
S704：当输入电压相角落在扇区III和IV时，A输出相的参考电流值i* oA和实际值ioA都是负值，此时的能量函数变化率表示为：
在一定的时间内，A相负载电流的幅值也会减小到零。
优选地，步骤S8，包括以下步骤：
S801：当负载电流保持为零的持续时间Tzero达到故障诊断所需时间时，再根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置，故障情况与输入电压周期和负载电流基波周期To有关，将故障持续时间表示为：
其中，Tin和To分别表示输入电压基波周期和负载电流基波周期；
定义故障相的负载电流通过箝位电路变到零所需的时间Tdischarge，故障相电流保持为零的时间Tzero表示为：
Tzero＝Tfault-Tdischarge (39)
故障诊断时间需要满足：
kfTs>Tε (40)
其中，kf是故障诊断的敏感系数；Ts是正常运行情况下输出相的参考电流的过零(输出相的参考电流小于阈值ε)时间；阈值ε是为了减少故障诊断的误报率而设置的一个小的正值；
综合(38)～(40)式，故障诊断所需时间应该满足：
S802：根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关位置，假设SAa发生开路故障，故障开关定位的逻辑可以表示为：
如果几个连续的采样周期内，满足|ioA|≤ε且SAa＝1，那么SAa发生开路故障，即：
优选地，步骤S9中当SAa开路时，矩阵变换器的运行的约束条件为：
从而余下的开关组合状态共有2*3*3＝18种；其余开关开路时，依此类推；
步骤S10中，在余下的所有开关组合状态中，不考虑与矩阵变换器输入电流相关的评价函数，即(28)式总评价函数中的λ取零值，重新确定具体故障下的评价函数为：
gF＝g2 (44)
步骤S11中，在矩阵变换器带故障运行过程中，采用有限集模型预测控制方法，即在每个采样周期计算余下的所有开关组合状态下的gF值，选择其中使gF最小值的开关组合状态作为下一采样周期的开关状态；实现故障后的满足***动、稳态性能指标的容错运行。
本发明具有以下有益效果：
1、本发明的基于预测控制的矩阵变换器开关开路故障诊断方法，实时监测三相负载电流，通过判断输出相电流保持为零，检测矩阵变换器***是否有故障发生，并且根据当时预测控制器的开关状态定位故障开关。本发明能及时诊断开路故障，以便控制器及时进行处理，从而延长使用寿命。
2、本发明的基于预测控制的矩阵变换器开关开路故障诊断方法，进一步还能对故障进行容错处理，当矩阵变换器的一个开关发生开路故障时，利用剩余的正常开关组成新的开关状态组合，然后继续使用预测控制思想选择最合适的开关状态，使矩阵变换器***在单开关开路故障的情况下继续运行。本发明可以避免因开路故障造成的器件损坏，进一步延长了器件的使用寿命。
除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照附图，对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是本发明优选实施例的基于预测控制的矩阵变换器开关开路故障诊断方法的流程示意图；
图2是本发明优选实施例的矩阵变换器的拓扑结构图；
图3是本发明优选实施例的矩阵变换器的模型预测控制策略的流程示意图；
图4是本发明优选实施例的矩阵变换器在正常工作条件下的输入/输出电压以及输入/负载电流实验波形图；
图5是本发明优选实施例的矩阵变换器的输入电压基波周期分扇区示意图；其中，Voltage为电压；Sector为扇区；
图6是本发明优选实施例的矩阵变换器开关SAa故障条件下的等效电路示意图；
图7是本发明优选实施例的矩阵变换器开关开路故障诊断的流程图；
图8是本发明优选实施例的矩阵变换器开关开路故障容错策略的流程图；
图9是发明优选实施例的矩阵变换器发生开路故障的条件下加入诊断与容错措施前后的负载电流波形对比示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
参见图1，本发明实施例的基于预测控制的矩阵变换器开关开路故障诊断方法，包括以下步骤：
S1：在矩阵变换器正常工作的情况下，依据采样周期采集矩阵变换器的输出电压、输入电压、输入电流与负载电流，建立输出电压与输入电压、输入电流与负载电流之间的关系模型；
S2：在满足约束条件的情况下，确定矩阵变换器的所有开关组合状态的数目；
S3：建立矩阵变换器负载端状态空间模型，获得负载电流在下一采样周期的预测值；
S4：建立矩阵变换器输入端滤波器状态空间模型，离散化后，获得下一采样周期输入电流与输入电压的预测值；
S5：确定评价函数；
S6：采用有限集模型预测控制策略，在每一采样周期，在所有开关组合状态中选择使评价函数值最小的开关状态，作为下一采样周期的开关状态；
S7：建立开关开路故障条件下的等效电路的状态空间模型，得到故障相负载电流在开关开路故障后的变化规律；
S8：对负载电流进行在线监测，当负载电流持续为零时，根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置。
以上步骤，通过实时监测三相负载电流，判断输出相电流保持为零(实际应用中，保持为零可以采用其持续时间是否超过阈值来判断)，检测矩阵变换器***是否有故障发生，并且根据当时预测控制器的开关状态定位故障开关，能及时诊断开路故障。
在实际应用中，在上述步骤的基础上，本发明另一优选实施例的基于预测控制的矩阵变换器开关开路故障诊断方法，增加容错运行步骤进行优化：
即在上述的步骤S8完成后，进行以下步骤：
S9：在识别故障开关的位置后，将所有开关组合状态中与故障开关的动作有关的开关状态去除，确定余下所有开关组合状态的数目；
S10：在余下的所有开关组合状态中，重新确定具体故障下的评价函数；
S11：采用有限集模型预测控制策略，在每一采样周期，在余下的所有开关组合状态中选择使评价函数值最小的开关组合状态，作为下一采样周期的开关状态。
这一优选实施例能进一步对故障进行容错处理，当矩阵变换器的一个开关发生开路故障时，利用剩余的正常开关组成新的开关状态组合，然后继续使用预测控制思想选择最合适的开关状态，使矩阵变换器***在单开关开路故障的情况下继续运行。本发明实施例可以避免因开路故障造成的器件损坏，进一步延长了器件的使用寿命。
下面将上述实施例应用到如图2所示的拓扑结构的矩阵变换器上，且以矩阵变换器的开关SAa发生开路故障为例，对本发明进行进一步分析：
矩阵变换器正常运行时，其主电路拓扑结构如图2所示；矩阵变换器有限集模型预测控制选择开关状态的过程如图3所示；实验参数如表1所示。
表1.实验参数表
本发明实施例的基于预测控制的矩阵变换器开关开路故障诊断方法，包括以下步骤：
S1：在矩阵变换器正常工作的情况下，依据采样周期采集矩阵变换器的输出电压、输入电压、输入电流与负载电流，建立输出电压与输入电压、输入电流与负载电流之间的关系模型如下：
其中，uoA、uoB和uoC分别为A、B和C三相的输出电压；uea,、ueb和uec分别为A、B和C三相的输入电压；iea、ieb和iec分别为A、B和C三相的输入电流；ioA、ioB和ioC分别为A、B和C三相的负载电流；S为开关状态矩阵；S中的元素SXY(X∈{A,B,C},Y∈{a,b,c})表示该开关的开关状态，SXY＝1表示开关SXY为闭合状态，SXY＝0表示开关SXY为断开状态。
S2：在满足约束条件的情况下，确定矩阵变换器的所有开关组合状态的数目。设立约束条件是为了确保矩阵变换器的安全运行。约束条件包括以下几个：
S201：输入端不短路，输出端不断路；
S202：在任意时刻，连接到同一相输出的三个双向开关中，有且仅有一个开关导通，另外两个开关关断；
S203：满足各个开关占空比的非负性；
在上述约束条件下，开关状态矩阵S中的每个元素满足：
可得矩阵变换器正常工作时的所有开关组合状态有3*3*3＝27种。
S3：建立矩阵变换器负载端状态空间模型，获得负载电流在下一采样周期的预测值。具体包括以下步骤：
S301：建立矩阵变换器负载端状态空间模型如下：
其中，io和uo分别为负载电流与输出电压；L和R分别表示负载的电感和电阻；
S302：定义一个开关周期Ts，根据向前欧拉公式得到负载电流在下一采样周期的预测值io k+1为：
其中，io k和uo k分别为当前采样时刻的负载电流值和输出电压值；系数F1和F2分别为：
S4：建立矩阵变换器输入端滤波器状态空间模型，离散化后，获得下一采样周期输入电流与输入电压的预测值。包括以下步骤：
S401：建立矩阵变换器输入端滤波器状态空间模型如下：
其中us、ue、is和ie分别表示电源电压、滤波电容电压、电源电流和输入电流；Ri、Li和Ci分别表示滤波电阻、电源和滤波电感以及滤波电容；
S402：结合(9)式和(10)式得到：
其中：
S403：将(11)式离散化，得到：
式中和分别表示在下一采样周期的滤波电容电压和电源电流值；和分别表示在当前采样时刻的电源电压、滤波电容电压、电源电流和输入电流值；且有：
S404：将(15)式代入(14)式，得到下一采样周期的输入电流以及滤波电容电压：
其中，E2、E3、D2和D3均为矩阵G的元素，E1、E4、D1和D4均为矩阵H的元素，由Ri、Li、Ci和采样周期Ts决定。
S5：确定评价函数，具体包括以下步骤：
S501：确定与矩阵变换器的输入电流相关的评价函数。
矩阵变换器输入电流的参考值为：
其中φ是输入电压的相角，为参考输入电流的幅值，并且有：
其中Usm为电源电压us的幅值；η为矩阵变换器的效率，由如下公式获得：
其中的Po和Pin分别表示输出功率和输入功率，分别由如下公式得到：
其中I* om表示参考负载电流i* o的幅值；期望输入电流is功率因数为1，所以式中θ为0；
S502：确定预测控制中关于输入电流is期望达到单位输入功率因数性能的评价函数g1为：
式中，△isα和△isβ分别表示两相静止坐标系下的输入电流在下一采样周期预测值is p与参考值is *偏差的α和β分量，由如下公式得到：
式中，Tabctoαβ是三相静止坐标系转换为两相旋转坐标系的变换矩阵，由下式得到：
S503：确定预测控制中与矩阵变换器负载电流相关的评价函数：
定义关于负载电流预测值期望达到跟踪参考负载电流性能的评价函数g2为：
式中，△ioα和△ioβ分别表示两相静止坐标系下的负载电流在下一采样周期预测值io p与参考值io *偏差的α和β分量，由如下公式得到：
S504：定义总评价函数为：
g＝λ*g1+g2 (28)
其中，λ是权重因子，表示输入电流的性能指标在整个矩阵变换器***的控制中的重要程度。
S6：在矩阵变换器运行过程中，采用有限集模型预测控制策略，在每个采样周期计算所有开关组合状态的总评价函数的g值，选择其中能使g值最小的一种开关状态，作为下一采样周期的开关状态，实现达到***期望的动态性能和稳态性能指标的运行。图4所示为矩阵变换器正常工作条件下的输入/输出电压以及输入/负载电流实验波形图。由图4可见，采用预测控制策略选择合适的开关状态可以使矩阵变换器的输入/负载电流正弦且很好地跟踪参考值，并且输入功率因数为1。
以下步骤为对矩阵变换器***进行实时的故障检测与诊断步骤：
S7：建立开关开路故障条件下的等效电路的状态空间模型，得到故障相负载电流在开关开路故障后的变化规律，具体包括以下步骤：
S701：假设如下条件S7011至S7013成立：
S7011：采样周期为Ts；
S7012：如图5所示，将输入电压周期分为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ和Ⅵ顺序排列的六个扇区，输入电压相角在扇区Ⅴ并且即将进入扇区Ⅵ；A输出相的参考电流值i* oA和实际值ioA都是正值；
S7013：SAa发生开路故障，在kth个采样周期Ts k的时间内，预测控制器为了缩小A输出相电流的实际值与参考值的偏差，始终给SAa导通信号。
则有，A输出相开路，A相负载电流ioA会流经箝位电路，电路如图6(a)所示。其电路状态空间模型描述如下：
式中，uN是负载中性点与电源接地点之间的电压。
输出相的电压表示如下：
其中，Ucp是箝位电路的电容电压，umax为输入三相电压的最大值。
由于A相负载电流流过箝位电路，所以Ucp远大于umax，表示为：
其中，umin表示输入三相电压的最小值。
根据(29)式和(30)式，得到A相负载电流ioA的变化率为：
由(31)可得：
将(33)代入(32)，得到：
当ioA>0时 (34)。
S703：构造能量函数v：
其变化率表示为：
于是，得到：
当ioA<0时 (37)
由(34)式和(37)式可知：当故障开关满足SAa＝1时，A相负载电流的幅值会越来越小，从而使其实际值与参考值之间的偏差越来越大；因为SAa在扇区VI和I连通A输出相对应的最大输入电压，预测控制器始终给开关SAa导通信号，然而，SAa因为故障无法导通，结果就是A相负载电流幅值减小并且连续几个采样周期内SAa＝1；在一定时间内，A相负载电流的幅值逐渐减小到零。
S704：当输入电压相角落在扇区III和IV时，A输出相的参考电流值i* oA和实际值ioA都是负值，等效电路如图6(b)所示，此时的能量函数变化率表示为：
在一定的时间内，A相负载电流的幅值也会减小到零。
步骤S8，包括以下步骤：
S801：当负载电流保持为零的持续时间Tzero达到故障诊断所需时间时，再根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置，故障情况与输入电压周期和负载电流基波周期To有关，将故障持续时间表示为：
其中，Tin和To分别表示输入电压基波周期和负载电流基波周期；
定义故障相的负载电流通过箝位电路变到零所需的时间Tdischarge，故障相电流保持为零的时间Tzero表示为：
Tzero＝Tfault-Tdischarge (39)。
为了区分故障相由于故障发生时电流保持为零和正常时电流过零点的情况，故障诊断时间需要满足：
kfTs>Tε (40)
其中，kf是故障诊断的敏感系数；Ts是正常运行情况下输出相的参考电流的过零(输出相的参考电流小于阈值ε)时间；阈值ε是为了减少故障诊断的误报率而设置的一个小的正值；
综合(38)～(40)式，故障诊断所需时间应该满足：
S802：根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关位置，假设SAa发生开路故障，故障开关定位的逻辑可以表示为：
如果几个连续的采样周期内，满足|ioA|≤ε且SAa＝1，那么SAa发生开路故障，即：
图7所示为以上矩阵变换器***进行实时的故障检测与诊断步骤的流程图，试验中故障诊断时间为1.4ms。
以下为在识别故障后，执行容错运行的步骤：
S9：在识别故障开关的位置后，将所有开关组合状态中与故障开关的动作有关的开关状态去除，确定余下所有开关组合状态的数目。例如当SAa开路时，矩阵变换器的运行的约束条件为：
从而余下的开关组合状态共有2*3*3＝18种；其余开关开路时，依此类推；
S10：在余下的所有开关组合状态中，不考虑与矩阵变换器输入电流相关的评价函数，即(28)式总评价函数中的λ取零值，重新确定具体故障下的评价函数为：
gF＝g2 (44)
S11：在矩阵变换器带故障运行过程中，采用有限集模型预测控制方法，参见图8，即在每个采样周期计算余下的所有开关组合状态下的gF值，选择其中使gF最小值的开关组合状态作为下一采样周期的开关状态；实现故障后的满足***动、稳态性能指标的容错运行。
图9所示为矩阵变换器发生开路故障的条件下加入诊断与容错措施前后的负载电流波形对比图。由图9(a)可见，矩阵变换器***前期正常运行时，负载电流基本正弦，谐波畸变率很低，并且很好地跟踪参考值，引入SAa开关开路故障后，负载电流ioA发生畸变并且迅速衰减为零，ioB和ioC也受到比较大的影响。如图9(b)所示，加入故障诊断措施，当引入SAa开关开路故障后，预测控制器在几个连续的采样周期内一直给SAa导通信号，然而，SAa由于故障无法导通，从而A输出相的电流会在几个连续的采样周期内保持为零，此时可诊断出SAa发生开路故障，该故障诊断过程可有图9(c)更明显地看出。在定位了故障开关后，加入容错措施，利用预测控制在剩余的开关组合状态中选择能使评价函数值最小的开关组合状态。由图可知，容错运行时，矩阵变换器的三相负载电流能基本平衡并且较好地跟踪参考正弦电流。
表2 是矩阵变换器正常、故障以及容错运行条件下的THD比较表。
表2
从表2中可以看出，开关SAa发生开路故障后，矩阵变换器A相负载电流的基波分量明显减少，谐波畸变率明显上升，B和C两相的负载电流的波形质量也受到了比较大的影响。在本发明实施例的容错运行状态下，A、B和C三相的负载电流的基波分量都有所增加，谐波畸变率都有所下降，波形质量得到提高，进一步证明了本发明的有效性。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于预测控制的矩阵变换器开关开路故障诊断方法，其特征在于，包括以下步骤：
S1：在矩阵变换器正常工作的情况下，依据采样周期采集矩阵变换器的输出电压、输入电压、输入电流与负载电流，建立输出电压与输入电压、输入电流与负载电流之间的关系模型；所述输出电压与输入电压、输入电流与负载电流之间的关系模型如下：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>e</mi> <mi>a</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>e</mi> <mi>b</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>e</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>S</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>a</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>b</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>c</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中，uoA、uoB和uoC分别为A、B和C三相的输出电压；uea,、ueb和uec分别为A、B和C三相的输入电压；iea、ieb和iec分别为A、B和C三相的输入电流；ioA、ioB和ioC分别为A、B和C三相的负载电流；S为开关状态矩阵；S中的元素SXY(X∈{A,B,C},Y∈{a,b,c})表示该开关的开关状态，SXY＝1表示开关SXY为闭合状态，SXY＝0表示开关SXY为断开状态；
S2：在满足约束条件的情况下，确定矩阵变换器的所有开关组合状态的数目；所述约束条件包括：
S201：输入端不短路，输出端不断路；
S202：在任意时刻，连接到同一相输出的三个双向开关中，有且仅有一个开关导通，另外两个开关关断；
S203：满足各个开关占空比的非负性；
在上述约束条件下，所述开关状态矩阵S中的每个元素满足：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>a</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>a</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>a</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
可得矩阵变换器正常工作时的所有开关组合状态有3*3*3＝27种；
S3：建立矩阵变换器负载端状态空间模型，获得负载电流在下一采样周期的预测值，包括以下步骤：
S301：建立矩阵变换器负载端状态空间模型如下：
<mrow> <mi>L</mi> <mfrac> <mrow> <msub> <mi>di</mi> <mi>o</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>u</mi> <mi>o</mi> </msub> <mo>-</mo> <msub> <mi>Ri</mi> <mi>o</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
其中，io和uo分别为负载电流与输出电压；L和R分别表示负载的电感和电阻；
S302：定义一个开关周期Ts，根据向前欧拉公式得到负载电流在下一采样周期的预测值io k+1为：
<mrow> <msubsup> <mi>i</mi> <mi>o</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msub> <mi>F</mi> <mn>1</mn> </msub> <msubsup> <mi>u</mi> <mi>o</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>F</mi> <mn>2</mn> </msub> <msubsup> <mi>i</mi> <mi>o</mi> <mi>k</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
其中，io k和uo k分别为当前采样时刻的负载电流值和输出电压值；系数F1和F2分别为：
<mrow> <msub> <mi>F</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>T</mi> <mi>s</mi> </msub> <mi>L</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>F</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mi>R</mi> </mrow> <mi>L</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
S4：建立矩阵变换器输入端滤波器状态空间模型，离散化后，获得下一采样周期输入电流与输入电压的预测值，包括以下步骤：
S401：建立矩阵变换器输入端滤波器状态空间模型如下：
<mrow> <msub> <mi>L</mi> <mi>i</mi> </msub> <mfrac> <mrow> <msub> <mi>di</mi> <mi>s</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>u</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>e</mi> </msub> <mo>-</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>C</mi> <mi>i</mi> </msub> <mfrac> <mrow> <msub> <mi>du</mi> <mi>e</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>i</mi> <mi>e</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中us、ue、is和ie分别表示电源电压、滤波电容电压、电源电流和输入电流；Ri、Li和Ci分别表示滤波电阻、电源和滤波电感以及滤波电容；
S402：结合(9)式和(10)式得到：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mrow> <msub> <mi>du</mi> <mi>e</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <msub> <mi>di</mi> <mi>s</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>A</mi> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>e</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mi>B</mi> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>s</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>e</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
其中：
<mrow> <mi>A</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mn>1</mn> <mo>/</mo> <msub> <mi>C</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>L</mi> <mi>i</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>B</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>C</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>/</mo> <msub> <mi>L</mi> <mi>i</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
S403：将(11)式离散化，得到：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>u</mi> <mi>e</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>i</mi> <mi>s</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>G</mi> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>u</mi> <mi>e</mi> <mi>k</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>i</mi> <mi>s</mi> <mi>k</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mi>H</mi> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>u</mi> <mi>s</mi> <mi>k</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>i</mi> <mi>e</mi> <mi>k</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
式中和分别表示在下一采样周期的滤波电容电压和电源电流值；和分别表示在当前采样时刻的电源电压、滤波电容电压、电源电流和输入电流值；且有：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>G</mi> <mo>=</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>AT</mi> <mi>s</mi> </msub> </mrow> </msup> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>H</mi> <mo>=</mo> <msup> <mi>A</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mi>G</mi> <mo>-</mo> <mi>I</mi> <mo>)</mo> </mrow> <mi>B</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
S404：将(15)式代入(14)式，得到下一采样周期的输入电流以及滤波电容电压：
<mrow> <msubsup> <mi>u</mi> <mi>e</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msub> <mi>E</mi> <mn>1</mn> </msub> <msubsup> <mi>u</mi> <mi>s</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>E</mi> <mn>2</mn> </msub> <msubsup> <mi>u</mi> <mi>e</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>E</mi> <mn>3</mn> </msub> <msubsup> <mi>i</mi> <mi>s</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>E</mi> <mn>4</mn> </msub> <msubsup> <mi>i</mi> <mi>e</mi> <mi>k</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>i</mi> <mi>s</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msub> <mi>D</mi> <mn>1</mn> </msub> <msubsup> <mi>u</mi> <mi>s</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>D</mi> <mn>2</mn> </msub> <msubsup> <mi>u</mi> <mi>e</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>D</mi> <mn>3</mn> </msub> <msubsup> <mi>i</mi> <mi>s</mi> <mi>k</mi> </msubsup> <mo>+</mo> <msub> <mi>D</mi> <mn>4</mn> </msub> <msubsup> <mi>i</mi> <mi>e</mi> <mi>k</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
其中，E2、E3、D2和D3均为矩阵G的元素，E1、E4、D1和D4均为矩阵H的元素，由Ri、Li、Ci和采样周期Ts决定；
S5：确定评价函数；包括以下步骤：
S501：确定与矩阵变换器的输入电流相关的评价函数；
矩阵变换器输入电流的参考值为：
<mrow> <msubsup> <mi>i</mi> <mi>s</mi> <mo>*</mo> </msubsup> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;phi;</mi> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mo>/</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>+</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mo>/</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
其中φ是输入电压的相角，为参考输入电流的幅值，并且有：
<mrow> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>&amp;eta;U</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> </msub> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;eta;U</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <msub> <mi>&amp;eta;R</mi> <mi>i</mi> </msub> <msup> <msubsup> <mi>I</mi> <mrow> <mi>o</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mn>2</mn> </msup> <mi>R</mi> </mrow> </msqrt> </mrow> <mrow> <mo>-</mo> <mn>2</mn> <msub> <mi>&amp;eta;R</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
其中Usm为电源电压us的幅值；η为矩阵变换器的效率，由如下公式获得：
<mrow> <mi>&amp;eta;</mi> <mo>=</mo> <mfrac> <msub> <mi>P</mi> <mi>o</mi> </msub> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
其中的Po和Pin分别表示输出功率和输入功率，分别由如下公式得到：
<mrow> <msub> <mi>P</mi> <mi>o</mi> </msub> <mo>=</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <msup> <msubsup> <mi>I</mi> <mrow> <mi>o</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mn>2</mn> </msup> <mi>R</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>U</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> </msub> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msup> <msubsup> <mi>I</mi> <mrow> <mi>s</mi> <mi>m</mi> </mrow> <mo>*</mo> </msubsup> <mn>2</mn> </msup> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
其中I* om表示参考负载电流i* o的幅值；期望输入电流is功率因数为1，所以式中θ为0；
S502：确定预测控制中关于输入电流is期望达到单位输入功率因数性能的评价函数g1为：
<mrow> <msub> <mi>g</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mi>&amp;Delta;i</mi> <mrow> <mi>s</mi> <mi>&amp;alpha;</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;Delta;i</mi> <mrow> <mi>s</mi> <mi>&amp;beta;</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
式中，Δisα和Δisβ分别表示两相静止坐标系下的输入电流在下一采样周期预测值is p与参考值is *偏差的α和β分量，由如下公式得到：
<mrow> <msub> <mi>&amp;Delta;i</mi> <mi>s</mi> </msub> <mo>=</mo> <msup> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;i</mi> <mrow> <mi>s</mi> <mi>&amp;alpha;</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&amp;Delta;i</mi> <mrow> <mi>s</mi> <mi>&amp;beta;</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>a</mi> <mi>b</mi> <mi>c</mi> <mi>t</mi> <mi>o</mi> <mi>&amp;alpha;</mi> <mi>&amp;beta;</mi> </mrow> </msub> <mrow> <mo>(</mo> <msubsup> <mi>i</mi> <mi>s</mi> <mi>P</mi> </msubsup> <mo>-</mo> <msubsup> <mi>i</mi> <mi>s</mi> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
式中，Tabctoαβ是三相静止坐标系转换为两相旋转坐标系的变换矩阵，由下式得到：
<mrow> <msub> <mi>T</mi> <mrow> <mi>a</mi> <mi>b</mi> <mi>c</mi> <mi>t</mi> <mi>o</mi> <mi>&amp;alpha;</mi> <mi>&amp;beta;</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> </msqrt> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>2</mn> </mfrac> </mtd> <mtd> <mrow> <mo>-</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>2</mn> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>
S503：确定预测控制中与矩阵变换器负载电流相关的评价函数：
定义关于负载电流预测值期望达到跟踪参考负载电流性能的评价函数g2为：
<mrow> <msub> <mi>g</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>&amp;Delta;i</mi> <mrow> <mi>o</mi> <mi>&amp;alpha;</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;Delta;i</mi> <mrow> <mi>o</mi> <mi>&amp;beta;</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>
式中，Δioα和Δioβ分别表示两相静止坐标系下的负载电流在下一采样周期预测值io p与参考值io *偏差的α和β分量，由如下公式得到：
<mrow> <msub> <mi>&amp;Delta;i</mi> <mi>o</mi> </msub> <mo>=</mo> <msup> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;i</mi> <mrow> <mi>o</mi> <mi>&amp;alpha;</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&amp;Delta;i</mi> <mrow> <mi>o</mi> <mi>&amp;beta;</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>a</mi> <mi>b</mi> <mi>c</mi> <mi>t</mi> <mi>o</mi> <mi>&amp;alpha;</mi> <mi>&amp;beta;</mi> </mrow> </msub> <mrow> <mo>(</mo> <msubsup> <mi>i</mi> <mi>o</mi> <mi>P</mi> </msubsup> <mo>-</mo> <msubsup> <mi>i</mi> <mi>o</mi> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
S504：定义总评价函数为：
g＝λ*g1+g2 (28)
其中，λ是权重因子，表示输入电流的性能指标在整个矩阵变换器***的控制中的重要程度；
S6：采用有限集模型预测控制策略，在每一采样周期，在所有开关组合状态中选择使评价函数值最小的开关状态，作为下一采样周期的开关状态；具体为：在矩阵变换器运行过程中，采用有限集模型预测控制策略，即在每个采样周期计算所有开关组合状态的总评价函数的g值，选择其中能使g值最小的一种开关状态，作为下一采样周期的开关状态，实现达到***期望的动态性能和稳态性能指标的运行；
S7：建立开关开路故障条件下的等效电路的状态空间模型，得到故障相负载电流在开关开路故障后的变化规律；
S8：对负载电流进行在线监测，当负载电流持续为零时，根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置。
2.根据权利要求1所述的故障诊断方法，其特征在于，所述步骤S8完成后，所述方法还包括：
S9：在识别所述故障开关的位置后，将所有开关组合状态中与故障开关的动作有关的开关状态去除，确定余下所有开关组合状态的数目；
S10：在余下的所有开关组合状态中，重新确定具体故障下的评价函数；
S11：采用有限集模型预测控制策略，在每一采样周期，在所述余下的所有开关组合状态中选择使评价函数值最小的开关组合状态，作为下一采样周期的开关状态。
3.根据权利要求2所述的故障诊断方法，其特征在于，所述步骤S7，包括以下步骤：
S701：假设如下条件S7011至S7013成立：
S7011：采样周期为Ts；
S7012：将输入电压周期分为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ和Ⅵ顺序排列的六个扇区，输入电压相角在扇区Ⅴ并且即将进入扇区Ⅵ；A输出相的参考电流值i* oA和实际值ioA都是正值；
S7013：SAa发生开路故障，在kth个采样周期Ts k的时间内，预测控制器为了缩小A输出相电流的实际值与参考值的偏差，始终给SAa导通信号；
则有，A输出相开路，A相负载电流ioA会流经箝位电路，其电路状态空间模型描述如下：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>Ri</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> <mo>-</mo> <mi>L</mi> <mfrac> <mrow> <msub> <mi>di</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>u</mi> <mi>N</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>Ri</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> <mo>-</mo> <mi>L</mi> <mfrac> <mrow> <msub> <mi>di</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>u</mi> <mi>N</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>Ri</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> <mo>-</mo> <mi>L</mi> <mfrac> <mrow> <msub> <mi>di</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>u</mi> <mi>N</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>
式中，uN是负载中性点与电源接地点之间的电压；
输出相的电压表示如下：
其中，Ucp是箝位电路的电容电压，umax为输入三相电压的最大值；
由于A相负载电流流过箝位电路，所以Ucp远大于umax，表示为：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>a</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>b</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mi>min</mi> </msub> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>a</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>b</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>p</mi> </mrow> </msub> <mo>&amp;GreaterEqual;</mo> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>min</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>
其中，umin表示输入三相电压的最小值；
根据(29)式和(30)式，得到A相负载电流ioA的变化率为：
<mrow> <mfrac> <mrow> <msub> <mi>di</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <mn>2</mn> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>p</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> <mo>-</mo> <mn>3</mn> <msub> <mi>Ri</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> </mrow> <mrow> <mn>3</mn> <mi>L</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>32</mn> <mo>)</mo> </mrow> </mrow>
由(31)可得：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>p</mi> </mrow> </msub> <mo>&amp;GreaterEqual;</mo> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>p</mi> </mrow> </msub> <mo>&amp;GreaterEqual;</mo> <msub> <mi>u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>33</mn> <mo>)</mo> </mrow> </mrow>
将(33)代入(32)，得到：
当ioA>0时 (34)
S703：构造能量函数v：
<mrow> <mi>v</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>Li</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>35</mn> <mo>)</mo> </mrow> </mrow>
其变化率表示为：
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>v</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>Li</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>36</mn> <mo>)</mo> </mrow> </mrow>
于是，得到：
当ioA<0时 (37)
由(34)式和(37)式可知：当故障开关满足SAa＝1时，A相负载电流的幅值会越来越小，从而使其实际值与参考值之间的偏差越来越大；因为SAa在扇区VI和I连通A输出相对应的最大输入电压，预测控制器始终给开关SAa导通信号，然而，SAa因为故障无法导通，结果就是A相负载电流幅值减小并且连续几个采样周期内SAa＝1；在一定时间内，A相负载电流的幅值逐渐减小到零；
S704：当输入电压相角落在扇区III和IV时，A输出相的参考电流值i* oA和实际值ioA都是负值，此时的能量函数变化率表示为：
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>v</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>i</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>u</mi> <mi>min</mi> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>p</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>B</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>o</mi> <mi>C</mi> </mrow> </msub> <mo>-</mo> <mn>3</mn> <msub> <mi>Ri</mi> <mrow> <mi>o</mi> <mi>A</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mn>3</mn> </mfrac> <mo>&lt;</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>37</mn> <mo>)</mo> </mrow> </mrow>
在一定的时间内，A相负载电流的幅值也会减小到零。
4.根据权利要求3所述的故障诊断方法，其特征在于，所述步骤S8，包括以下步骤：
S801：当所述负载电流保持为零的持续时间Tzero达到故障诊断所需时间时，再根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关的位置，故障情况与输入电压周期和负载电流基波周期To有关，将故障持续时间表示为：
<mrow> <msub> <mi>T</mi> <mrow> <mi>f</mi> <mi>a</mi> <mi>u</mi> <mi>l</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mo>{</mo> <mfrac> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mn>3</mn> </mfrac> <mo>,</mo> <mfrac> <msub> <mi>T</mi> <mi>o</mi> </msub> <mn>2</mn> </mfrac> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>38</mn> <mo>)</mo> </mrow> </mrow>
其中，Tin和To分别表示输入电压基波周期和负载电流基波周期；
定义故障相的负载电流通过箝位电路变到零所需的时间Tdischarge，故障相电流保持为零的时间Tzero表示为：
Tzero＝Tfault-Tdischarge (39)
故障诊断所需时间Tε需要满足：
kfTs>Tε (40)
其中，kf是故障诊断的敏感系数；Ts是正常运行情况下输出相的参考电流的过零时间；阈值ε是为了减少故障诊断的误报率而设置的一个小的正值；
综合(38)～(40)式，故障诊断所需时间Tε应该满足：
<mrow> <mi>min</mi> <mo>{</mo> <mfrac> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mn>3</mn> </mfrac> <mo>,</mo> <mfrac> <msub> <mi>T</mi> <mi>o</mi> </msub> <mn>2</mn> </mfrac> <mo>}</mo> <mo>-</mo> <msub> <mi>T</mi> <mrow> <mi>d</mi> <mi>i</mi> <mi>s</mi> <mi>c</mi> <mi>h</mi> <mi>arg</mi> <mi>e</mi> </mrow> </msub> <mo>&gt;</mo> <msub> <mi>k</mi> <mi>f</mi> </msub> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>&gt;</mo> <msub> <mi>T</mi> <mi>&amp;epsiv;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>41</mn> <mo>)</mo> </mrow> </mrow>
S802：根据有限集模型预测控制选择的开关状态进行故障诊断，识别故障开关位置，假设SAa发生开路故障，故障开关定位的逻辑可以表示为：
如果几个连续的采样周期内，满足|ioA|≤ε且SAa＝1，那么SAa发生开路故障，即：
5.根据权利要求4所述的故障诊断方法，其特征在于，所述步骤S9中当SAa开路时，矩阵变换器的运行的约束条件为：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>A</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>a</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>B</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>a</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>b</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>S</mi> <mrow> <mi>C</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>43</mn> <mo>)</mo> </mrow> </mrow>
从而余下的开关组合状态共有2*3*3＝18种；其余开关开路时，依此类推；
所述步骤S10中，在余下的所有开关组合状态中，不考虑与矩阵变换器输入电流相关的评价函数，即(28)式总评价函数中的λ取零值，重新确定具体故障下的评价函数为：
gF＝g2 (44)
所述步骤S11中，在矩阵变换器带故障运行过程中，采用有限集模型预测控制方法，即在每个采样周期计算余下的所有开关组合状态下的gF值，选择其中使gF最小值的开关组合状态作为下一采样周期的开关状态；实现故障后的满足***动、稳态性能指标的容错运行。
CN201510995957.9A 2015-12-28 2015-12-28 基于预测控制的矩阵变换器开关开路故障诊断方法 Active CN105548792B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201510995957.9A CN105548792B (zh) 2015-12-28 2015-12-28 基于预测控制的矩阵变换器开关开路故障诊断方法
PCT/CN2016/105673 WO2017114001A1 (zh) 2015-12-28 2016-11-14 基于预测控制的矩阵变换器开关开路故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510995957.9A CN105548792B (zh) 2015-12-28 2015-12-28 基于预测控制的矩阵变换器开关开路故障诊断方法

Publications (2)

Publication Number Publication Date
CN105548792A CN105548792A (zh) 2016-05-04
CN105548792B true CN105548792B (zh) 2017-12-15

Family

ID=55828113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510995957.9A Active CN105548792B (zh) 2015-12-28 2015-12-28 基于预测控制的矩阵变换器开关开路故障诊断方法

Country Status (2)

Country Link
CN (1) CN105548792B (zh)
WO (1) WO2017114001A1 (zh)

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105548792B (zh) * 2015-12-28 2017-12-15 中南大学 基于预测控制的矩阵变换器开关开路故障诊断方法
CN106093677B (zh) * 2016-06-08 2018-10-30 中国矿业大学 一种三电平有源滤波器igbt开路故障定位方法
CN106330034B (zh) * 2016-06-14 2019-03-29 北方工业大学 间接矩阵变换器简化采样的模型预测控制方法
CN106125012B (zh) * 2016-06-28 2019-03-19 南京航空航天大学 一种间接式矩阵变换器故障诊断方法
CN106208737B (zh) * 2016-08-24 2019-06-18 中南大学 基于三次谐波注入矩阵变换器的模型预测电流控制方法
GB2557294B (en) 2016-12-05 2022-03-30 Itt Mfg Enterprises Llc Matrix converter control method and system
CN106602906B (zh) * 2016-12-07 2019-03-05 中南大学 三电平整流器开路故障注入方法及故障注入器
CN107064707B (zh) * 2016-12-16 2019-05-31 南京航空航天大学 双级矩阵变换器在反向发电状态时的开路故障检测方法
CN107070269B (zh) * 2017-02-14 2019-11-26 中南大学 一种三相四线逆变器的控制方法及***
CN107219455B (zh) * 2017-05-22 2019-06-28 中国电子科技集团公司第四十一研究所 一种矩阵开关继电器故障检测装置及检测方法
CN107612446B (zh) * 2017-11-03 2020-06-02 天津工业大学 一种内置式永磁同步电机模型预测转矩控制方法
CN108132419B (zh) * 2017-12-06 2020-12-08 福建奥通迈胜电力科技有限公司 基于电流突变信息阵的简易配电网主动故障研判方法
CN110391670A (zh) * 2018-04-16 2019-10-29 广西师范大学 一种微电网***灾变预测方法及装置
US10838008B2 (en) * 2018-07-13 2020-11-17 GM Global Technology Operations LLC Battery switch testing system and method
CN109375029B (zh) * 2018-08-31 2020-12-01 中南大学 一种两电平变流器***开关器件开路故障诊断方法与***
CN109462251B (zh) * 2018-10-15 2022-07-12 上海电力学院 一种新型微网故障的限流方法
CN109815526B (zh) * 2018-12-07 2023-06-23 杭州电子科技大学 一种Boost型变换器大信号建模方法
CN109752648B (zh) * 2018-12-25 2020-11-03 河北科技师范学院 一种有载分接开关切换电弧能量在线监测方法及装置
CN109726491B (zh) * 2019-01-03 2022-09-27 陕西理工大学 基于宽负载范围Buck变换器输出纹波电压精确建模方法
CN110069741B (zh) * 2019-04-22 2023-03-17 济南大学 高压断路器机械特性模糊综合评价的校正方法及***
CN110532587B (zh) * 2019-07-11 2023-01-10 杭州电子科技大学 一种双重移相调制下双有源全桥变换器大信号建模方法
CN110543662B (zh) * 2019-07-18 2022-09-02 陕西理工大学 宽负载范围非最小相位开关Boost变换器参数优化方法
CN110880752B (zh) * 2019-10-24 2023-04-18 东北大学 一种upfc有限控制集模型预测控制方法
CN110737996B (zh) * 2019-10-28 2023-05-26 中国大唐集团科学技术研究院有限公司西北电力试验研究院 一种高压断路器分合闸线圈电流识别方法
CN110909518B (zh) * 2019-11-06 2023-04-11 国网福建省电力有限公司 基于改进平均值模型的mmc可靠性分析方法
CN110620405A (zh) * 2019-11-08 2019-12-27 南京航空航天大学 一种多电源供电方式下的矩阵式配电***及控制方法
CN111190066B (zh) * 2020-01-14 2021-01-05 中南大学 矩阵变换器电机驱动***的故障诊断方法
CN111539275B (zh) * 2020-04-14 2020-11-20 中南民族大学 一种基于负载特性可视化的电气负载分类方法和***
CN111463772B (zh) * 2020-04-22 2024-05-31 中国船舶重工集团公司第七0四研究所 可抑制微网母线电压波动的储能接口变换器控制方法
CN111983414B (zh) * 2020-08-12 2022-07-29 中南大学 针对轨道列车牵引变流器的开路故障诊断方法及***
CN112100946B (zh) * 2020-08-21 2021-05-25 北京科技大学 一种基于故障在线模拟的逆变器开路故障诊断方法及装置
CN111983451B (zh) * 2020-08-24 2023-01-31 宿迁学院 一种基于搅拌器的无刷直流电机开路故障诊断方法
CN112054656A (zh) * 2020-09-25 2020-12-08 安徽理工大学 一种固定开关频率的模型预测控制算法
CN112327206B (zh) * 2020-10-27 2023-06-27 河南理工大学 一种用于三电平逆变器的故障诊断方法
CN112217410B (zh) * 2020-11-17 2022-04-08 中南大学 三电平逆变器开路故障容错控制方法及***
CN112733320B (zh) * 2020-12-09 2022-03-22 淮阴工学院 一种基于delta算子的升压变换器执行器故障检测方法
CN112701676B (zh) * 2020-12-15 2023-07-07 深圳供电局有限公司 一种熄弧角预测与换相失败抑制方法
CN113030722B (zh) * 2021-03-05 2023-08-29 青岛大学 一种多相电机驱动***开路故障在线检测方法
CN113107975B (zh) * 2021-04-07 2022-04-26 华中科技大学 用于磁轴承的绕组控制器断路故障定位、容错方法和***
CN112994110B (zh) * 2021-04-25 2023-04-11 郑州轻工业大学 一种lc滤波型并网逆变器无参数预测电容电压控制方法
CN113504435B (zh) * 2021-07-20 2022-07-08 中南大学 三电平逆变器开路故障诊断方法及***
CN113625563B (zh) * 2021-08-06 2024-03-15 江南大学 一种直流电动机量化迭代学习容错控制方法
CN113794407B (zh) * 2021-08-24 2023-06-30 西北工业大学 三相两电平pmsm驱动逆变器单管开路故障的控制方法
CN114362548B (zh) * 2021-12-20 2023-10-13 中国矿业大学 一种双级式矩阵变换器最优开关序列模型预测控制算法
CN114859855B (zh) * 2022-04-22 2023-03-14 大连理工大学 基于参数依赖Lyapunov函数的汽车发动机LPV***故障诊断装置
CN115149806B (zh) * 2022-06-21 2024-03-12 西北工业大学 一种针对交错并联Boost变换器的自适应模型预测控制方法
CN115099345A (zh) * 2022-06-27 2022-09-23 东南大学 一种基于自适应聚类算法的mmc子模块故障监测方法
CN115267474B (zh) * 2022-07-21 2023-06-20 山东大学 一种基于特勒根定理的电力变流器故障诊断方法及***
CN115248128B (zh) * 2022-08-19 2023-07-14 中南大学 一种基于级联观测器的列车轮对空转故障诊断方法及***
TWI827218B (zh) * 2022-08-26 2023-12-21 台達電子工業股份有限公司 開關短路診斷方法
CN117318483B (zh) * 2023-11-30 2024-02-09 湖南昇辉新能源技术有限公司 一种tcm模式变换器的预测控制方法及***
CN117783929B (zh) * 2024-02-26 2024-05-03 西北工业大学 一种直流变换器的开关管故障诊断方法
CN118068767B (zh) * 2024-04-19 2024-07-19 中国民航大学 一种g1000航电***数据自动采集注入装置及方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008086094A (ja) * 2006-09-27 2008-04-10 Meidensha Corp マトリクスコンバータ
KR101207749B1 (ko) * 2011-02-25 2012-12-03 아주대학교산학협력단 간접 매트릭스 컨버터의 스위치 고장진단장치 및 방법
US10069430B2 (en) * 2013-11-07 2018-09-04 Regents Of The University Of Minnesota Modular converter with multilevel submodules
CN104682351B (zh) * 2015-02-13 2018-07-13 中南大学 基于数学构造调制的矩阵变换器开路故障容错控制方法
CN105548792B (zh) * 2015-12-28 2017-12-15 中南大学 基于预测控制的矩阵变换器开关开路故障诊断方法

Also Published As

Publication number Publication date
WO2017114001A1 (zh) 2017-07-06
CN105548792A (zh) 2016-05-04

Similar Documents

Publication Publication Date Title
CN105548792B (zh) 基于预测控制的矩阵变换器开关开路故障诊断方法
CN207819564U (zh) 带有电流检测功能的故障监控***
CN108872882B (zh) 一种三电平级联逆变器的故障诊断装置及其诊断方法
CN104682351B (zh) 基于数学构造调制的矩阵变换器开路故障容错控制方法
CN104300516B (zh) 一种基于Buck变换器的单向变压型高压直流断路器
CN107643447A (zh) 一种车辆绝缘检测电路及方法
CN104267296B (zh) 基于mmc的statcom故障诊断方法
CN103078400A (zh) 适用于大容量mmc柔性直流输电***的桥臂汇总保护***
CN107064707B (zh) 双级矩阵变换器在反向发电状态时的开路故障检测方法
CN102567783A (zh) 并联混合型电能质量调节器专家故障分析与诊断方法
CN106125012B (zh) 一种间接式矩阵变换器故障诊断方法
CN109541377A (zh) 一种诊断hbsm内部igbt开路故障的新方法
CN208520978U (zh) 一种三电平级联逆变器的故障诊断装置
CN101895155A (zh) 具有故障检测和自动重组功能的储能装置及重组方法
CN107359688A (zh) 供电设备的故障处理方法及装置
CN110333426A (zh) 一种模块化多电平储能***开路故障诊断装置及方法
CN112994429A (zh) 一种输入并联输出串联型boost变换器的故障容错控制方法
CN107092213A (zh) 一种具有故障电弧检测功能的无源滤波装置及方法
CN104901574A (zh) 桥式电路及其短路保护方法
CN112034392B (zh) 一种基于滞环控制的有源电力滤波器的故障诊断方法
CN112234576B (zh) 基于滑窗特征提取的模块化多电平换流器故障诊断方法
CN112710922B (zh) 一种双有源桥dc-dc变换器的开路故障诊断方法
Tang et al. Four-state active fault diagnosis method for cascaded H-bridge multilevel inverter
CN108267684A (zh) 一种变频器故障诊断方法
CN111130072A (zh) 直流***中一种新型主动式保护装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant