MMC整流器混合FCS-MPC控制策略及Simulink复现指南

搞多电平变换器的兄弟应该都有过这种体验：一篇论文看完公式觉得懂了，真到自己动手复现的时候，光是Simulink里那个模块化多电平换流器（MMC）就够折腾一晚上。去年我复现一篇IEEE Trans. Power Electronics上关于混合有限集模型预测控制（FCS-MPC）的MMC整流论文，从拓扑搭到控制落地，中间踩坑踩到怀疑人生。网上能找到的模型要么是传统PI加调制器，要么是简化到只剩电流跟踪的MPC，"混合"到底混合了什么，几乎没人讲透。这篇文章就把整套Simulink仿真模型的搭建过程、控制策略的实现细节、参数配置和调参经验完整梳理一遍，算是一份可以直接对着做的复现笔记。项目核心就一句话：用混合FCS-MPC驱动MMC完成整流，在一个控制框架里协调好并网电流、电容电压均衡和桥臂环流抑制这三个问题。如果你正在做多电平变换器方向的课题、需要复现SCI论文，或者刚接触MPC想找一套能跑的Simulink模型，这篇内容应该能帮你省下不少时间。

1. 项目概述与整体设计思路

1.1 为什么选MMC和FCS-MPC这个组合

模块化多电平换流器（Modular Multilevel Converter, MMC）在中高压电力电子领域应用非常广，从柔性直流输电到中压电机驱动、从新能源并网到电能质量治理，都能看到它的身影。它和普通两电平变换器最大的区别是：用多个子模块级联取代了单个功率开关管，能输出接近正弦的阶梯波，谐波含量低，器件应力小，模块化设计也方便冗余。

选择MMC来做整流电路，目的很明确：验证FCS-MPC在复杂拓扑上的可行性和性能边界。因为MMC的控制难点恰恰是它的优势带来的——子模块数量多、状态变量多、每个桥臂又有环流问题，这给控制器的设计留足了发挥空间。

再说FCS-MPC，全称是Finite Control Set Model Predictive Control，有限控制集模型预测控制。它的思路比较直接：变换器可输出的开关状态是有限多个，在每个控制周期内把这所有状态全部枚举一遍，用系统模型预测下一时刻的变量值，然后带入代价函数，选出让代价函数最小的开关状态，直接驱动开关管。

相比传统PI加调制器的方案，FCS-MPC的最大亮点有三点。第一，动态响应快，因为它本质上是一个实时滚动优化，没有PI那种渐进逼近的滞后；第二，不需要线性调制器，输出的就是开关信号，省掉了SPWM、SVPWM这一层；第三，多目标控制很自然，代价函数里可以同时塞进电流误差、电容电压误差、环流误差，权重系数一调就能权衡。这些特性在MMC这种强非线性的多变量系统里特别合适。

1.2 复现SCI论文模型的整体思路

复现论文最忌讳一上来就把整个三相模型搬进Simulink，那只会让自己陷入无穷无尽的报错中。我的做法是先搭建单相MMC整流电路，把控制算法验证通过，再扩展成三相结构。

单相的好处是拓扑简单、状态变量少、波形现象直观，每个控制目标都能单独观测。论文里的公式在单相条件下已经能完整体现FCS-MPC的核心逻辑，等单相模型跑出正确的正弦电流、电容电压均衡、环流受控这三项结果后，再上三相就是水到渠成的事。

具体拆解下来，复现工作分成四层：

• 拓扑层：MMC功率电路，包括子模块、桥臂电感、交流侧电源和负载。
• 控制层：混合FCS-MPC算法，负责产生各桥臂的开关信号。
• 辅助层：电容电压均衡排序、环流抑制附加控制。
• 验证层：稳态波形、动态响应、谐波分析，与论文结果对照。

每一层都独立调试，能单独验证，最后再联调。这个分层思路也是我在这次复现中最大的收获，连跑不通的时候都能快速定位问题出在哪一层。

2. MMC整流电路的工作原理与控制目标

2.1 单相MMC拓扑结构与子模块状态

单相MMC整流器的拓扑结构可以这样理解：一条直流母线，分成上下两个桥臂，每个桥臂由N个子模块串联，再串一个桥臂电感Larm，上下桥臂在中点连接，中点接交流电源和负载。

子模块是MMC的最小单元，最常用的是半桥结构，由两个IGBT和一个直流电容构成。上下两个开关管的开关组合决定输出电压：

• 上管开通、下管关断：子模块投入，输出端电压等于电容电压Vc；
• 上管关断、下管开通：子模块切除，输出端电压约等于0；
• 两个都关断：闭锁状态，通常在启动或故障保护时使用。

所以每个桥臂能输出的电压台阶，取决于该桥臂当前投入了多少个子模块。理想情况下，上下桥臂投入的子模块数之和保持恒等于N，这样直流母线电压稳定在N倍的电容标称电压附近。

桥臂电流有两个关键分量：输出电流io=iu-il，以及环流ic=(iu+il)/2。环流是上下桥臂之间流动的电流，它不参与对外功率传输，却会给器件带来额外的电流应力和谐波问题。这部分在2.2节展开。

单相MMC整流电路工作在整流状态时，交流电源通过MMC的交流侧端口向直流负载供电。控制的目标就是要让输入电流呈现正弦波形且与电源电压同相位，实现单位功率因数整流，同时维持直流侧电压稳定。

2.2 三个核心控制目标的关系

混合FCS-MPC的"混合"之处，在于它要在一个控制框架中同时处理多个目标，而这些目标之间的耦合关系很微妙。

第一个目标是并网电流控制。整流器的最终任务是把交流能量转成直流，那么输入电流的质量直接决定谐波指标和功率因数。FCS-MPC对电流的控制核心是让实际电流在下一时刻尽量接近参考值。参考电流的幅值通常由直流电压外环给出，当直流侧负载变化或电压参考变化时，内环电流能快速跟随。

第二个目标是子模块电容电压均衡。每个子模块都有一个直流电容，如果某个电容电压偏高或偏低，分担不均会导致子模块间应力不同，严重时甚至损坏器件。MPC的代价函数里可以加入电容电压误差项，对偏离平均电压较大的子模块加大惩罚，从而在开关状态选择时优先平衡它们。

第三个目标是桥臂环流抑制。环流的产生源于上下桥臂瞬时电压不相等，它主要是二次谐波分量，在桥臂电感上产生额外损耗，还会让桥臂电流畸变、器件开关应力增大。传统做法是加独立的环流抑制器，而在混合FCS-MPC框架下，可以把环流预测误差作为一个代价项，或者用附加的环流闭环把二次谐波压下去。

这三个目标存在博弈：电流跟踪要的是快速响应，这往往带来较大的电压跳变；电压均衡要频繁切换子模块，这会影响开关损耗；环流抑制则要精细调节桥臂电压，又会和电流控制的效果叠加。混合控制的价值就是通过权重系数把这些矛盾的指标拉到一个可接受的平衡点上。

3. 混合FCS-MPC控制策略的核心实现

3.1 系统模型的离散化推导

FCS-MPC的基础是预测模型。模型不需要太复杂，关键是能准确描述控制周期内的电压电流关系。以单相MMC整流器为例，忽略开关损耗和线路电阻（或者用等效电阻Rs简化）后，可以写出两个核心微分方程。

输出电流的微分方程：

[ L_s \frac{di_o}{dt} + R_s i_o = v_s - v_o ]

其中vs是交流电源电压，vo是MMC交流侧输出电压，由上下桥臂输出电压之差决定：

[ v_o = \frac{v_l - v_u}{2} ]

环流的微分方程：

[ L_{arm} \frac{di_c}{dt} = \frac{V_{dc}}{2} - \frac{v_u + v_l}{2} ]

这两个方程描述了整流器最主要的动态特性。用前向欧拉法做离散化，采样周期Ts足够小（一般取几十微秒级别）时，可以把微分方程改写成一阶差分形式：

[ i_o(k+1) = i_o(k) + \frac{T_s}{L_s} \left[ v_s(k) - v_o(k) - R_s i_o(k) \right] ]

[ i_c(k+1) = i_c(k) + \frac{T_s}{L_{arm}} \left[ \frac{V_{dc}(k) - v_u(k) - v_l(k)}{2} \right] ]

这就是预测模型。给定当前测量值（io、ic、Vdc、vs）和候选开关状态（即候选的vu、vl组合），就能预测下一采样时刻的io和ic。

子模块电容电压的预测也类似。第j个子模块的电容电压变化率取决于它的开关状态sj和所在桥臂的电流方向：

[ v_{C,j}(k+1) = v_{C,j}(k) + \frac{T_s}{C} \cdot s_j(k) \cdot i_{arm}(k) ]

其中iarm是子模块所在桥臂的电流（上桥臂用iu，下桥臂用il）。C是子模块电容容值。

3.2 代价函数设计与权重系数的选择

代价函数是整个MPC的核心，它把所有控制目标压缩成一个标量。我在复现中用到的代价函数形式如下：

[ J = \lambda_1 [i_{o,ref}(k+1) - i_o(k+1)]^2

• \lambda_2 [i_{c,ref}(k+1) - i_c(k+1)]^2
• \lambda_3 \sum_{j=1}^{2N} [v_{C,avg}(k) - v_{C,j}(k+1)]^2 ]

第一项是输出电流误差项，保证整流电流跟踪参考值；第二项是环流误差项，参考值ic,ref通常设为零，因为理想情况下环流不存在；第三项是电容电压均衡项，让每个子模块电容电压向平均值靠拢。

权重系数λ1、λ2、λ3怎么调，这里有一句说一句，论文里一般只给出结果不告诉过程，实际调试时我一般遵守两个原则。

第一个原则是量纲统一。各项误差的数值数量级差别很大，比如输出电流可能是几十安培，电容电压误差可能只有几伏，如果不做标幺化或者归一化，某个项会直接压倒其他项。所以我在模型里先做基准值归一化，把电流、电压都折算到标幺值范围内，再调权重才有意义。

第二个原则是优先级递进。输出电流跟踪是主目标，权重必须最大；环流抑制次之；电容电压均衡项相对影响最小。我调试时先令λ1=1、λ2=0、λ3=0，把电流控制跑通、波形调准，再加入环流项，从λ2=0.1往上加，观察桥臂电流谐波变化；最后加λ3，从λ3=0.01往上调，看电容电压波动范围能不能压到5%以内。如果某两个目标互相冲突严重，就把冲突目标的权重降低而不是盲目加大。

权重系数调试的最终目标不是让某一项指标极致优秀，而是让三个指标同时满足工程可接受范围。比如电流THD小于5%，电容电压波动小于10%，环流分量小于额定电流的10%——这种平衡才是论文里所谓"混合控制优越性"的验证点。

3.3 有限控制集的构建与在线优化

FCS-MPC的"有限集"在MMC中体现为开关状态的组合数量。单个半桥子模块有2种有效开关状态，上下桥臂共2N个子模块，理论上有4的N次方种组合。

但实际控制中为了维持直流电压稳定，会让上下桥臂投入的子模块总数恒等于N，也就是满足约束：

[ n_{up} + n_{down} = N ]

这个约束下，有效组合数为排列组合C(2N, N)。我复现时用了N=4，也就是C(8,4)=70种组合。在每个控制周期里，MATLAB Function块内部分别计算这70种组合对应的预测值，并求出代价函数值，然后取最小值对应的开关模式输出。这个循环计算量对现代编译后的Simulink仿真来说完全扛得住，仿真速度也还能接受。

当然，如果子模块数很大，比如桥臂数N=10以上，组合数会达到184756种，单纯枚举就会非常吃力。这就是为什么论文中会引入排序法或分支定界法来压缩搜索空间。混合FCS-MPC在这里的优势是：可以把"代价函数选状态"和"排序法选子模块"结合起来，先把桥臂级需要投入的模块数量算出来，再通过排序法决定具体哪几个SM投入，计算量大幅下降。

3.4 环流抑制的实现方式

环流抑制这部分我单独拿出来说，是因为它最容易踩坑。在纯FCS-MPC方案里，环流抑制通常靠代价函数中的环流误差项实现，但实际效果往往不如预期，因为环流的动态时间尺度比输出电流慢，直接放在同一个代价函数里容易被快速电流项压制。

这也是"混合"概念的另一个体现：我最后采用的是FCS-MPC电流跟踪主环+独立环流抑制闭环的组合。环流闭环提取桥臂电流中的环流分量ic=(iu+il)/2，减去参考值ic_ref=0，误差送入比例谐振（PR）控制器，谐振频率设置在二倍频100Hz，输出一个附加的桥臂电压修正量，叠加到MPC输出的电压指令上。

如果不想用PR控制器，也可以用准比例谐振或者在代价函数中单独提高环流项的权重。但从工程复现角度看，独立的环流抑制闭环更容易调参、更容易和论文对标，也更容易解释波形改善的原因。

3.5 控制算法整体数据流

这个模型的控制流程可以概括为一个闭环逻辑：直流电压外环产生输出电流参考幅值，经过锁相同步后得到瞬时电流参考；FCS-MPC模块读取当前测量到的电流、电压、电容电压，枚举有限状态集合，预测未来值，计算代价函数，选出最优开关状态；环流抑制器独立计算附加电压修正量，叠加到桥臂电压指令上，最终生成驱动信号送进功率电路。整个闭环在每个采样周期（比如100微秒）内执行一次，滚动更新。

4. Simulink模型搭建与参数配置

4.1 功率电路搭建

功率电路用Simscape Electrical搭建，这是最能反映器件电平细节的方案。每个子模块用两个IGBT器件、两个二极管和一个直流电容按半桥结构连接，电容电压信号通过电压传感器引出来。

子模块数量少的时候可以直接复制粘贴，但数量多了就必须封装成子系统。我给子模块做成了带仿真端口的Mask封装，输入是开关脉冲，输出是电容电压，外部输入桥臂电流。这样桥臂上一串子模块通过总线连接，model结构清晰很多。

桥臂电感Larm直接用Simscape的串联R-L元件。交流侧电源用交流电压源，负载侧用电阻负载或电阻电感负载，方便测试动态响应。直流母线电压由MMC自身整流建立，不需要外部直流源。

需要注意电源连接极性。上下桥臂串联后，上桥臂首端接直流正端，下桥臂末端接直流负端，上下桥臂连接点接交流电源的一端，另一端接负载。接反了仿真一瞬间就会报错或者产生极大电流震荡。

4.2 控制系统搭建

控制系统全部放在Simulink普通域，用MATLAB Function块实现FCS-MPC算法。这个块的核心优势是支持循环和复杂逻辑，比用一堆Simulink基本模块搭MPC要高效得多。

MATLAB Function块内部结构我整理成四个步骤：

第一步，读取外部测量信号，包括交流电流io、上下桥臂电流iu和il、直流母线电压Vdc、交流电源电压vs、各子模块电容电压vc数组，以及直流电压外环给的电流参考幅值。

第二步，根据锁相环给出相位信息更新电流参考值：io_ref=Im_ref*sin(wt+theta)。锁相环用Simulink标准三相PLL或者单相PLL都可以。

第三步，枚举候选开关状态。先根据当前n_up=0到N生成上下桥臂投入数量组合，满足n_up+n_down=N；对每种组合计算对应的vu、vl预测值，再计算出下一时刻的io和ic预测值，同时根据桥臂电流方向预测电容电压变化。

第四步，计算代价函数，选出最优状态组合，输出相应的开关信号。开关信号是长度等于2N的向量，每一位对应一个子模块的开关管驱动。

这个流程里最难的一点是延迟补偿。实际数字控制里，从采样到执行开关状态穿过一个控制周期，直接施加k时刻算出的开关状态会有一步延迟。我的做法是用两步预测法：先用k时刻测量值预测k+1时刻状态，再利用预测出的k+1状态预测k+2时刻状态，用k+2的结果选最优开关状态。这样算出来导致的跟踪效果要好得多，尤其是采样频率不高时差异很明显。

4.3 典型参数配置清单

参数配置是复现论文的关键一环，给出我复现时的典型参数，你可以根据自己的论文目标参数进行调整。

直流电压外环是经典PI控制，输出限幅在电流参考幅值的合理范围（比如2倍额定电流内）。PI参数我用的是Kp=0.5，Ki=20，实际调整时观察直流电压阶跃响应再来微调。注意这里不是MPC的核心，用简单PI就能满足。

4.4 仿真求解器设置

Simulink求解器设置是决定仿真能否跑通的大前提。我踩过最狠的坑就是把求解器设成定步长离散而Simscape里又有连续积分环节，导致报错"Simscape does not support discrete solver"。

最终选择：功率电路（Simscape域）用变步长连续求解器ode23tb，控制部分用触发子系统以Ts为周期采样执行。这样既保证了Simscape器件的连续动态能正确仿真，又让MPC按固定的控制周期滚动计算。

控制信号的传递要在变步长仿真中实现固定采样，需要加Zero-Order Hold模块。把测量信号先经过ZOH（采样时间Ts）再进入MATLAB Function，保证每个控制周期读取的是同一时刻的测量值。开关信号输出端如果不加比例放大，Simscape IGBT模块的G端口可以直接接收逻辑电平，我的做法是加一个整数到逻辑转换。

步长精度方面，ode23tb的相对容差我设成1e-4，再小仿真速度下降明显，且结果差异不大。如果某一时刻波形突然发散，先把相对容差调到1e-5试试，往往就能稳住。

4.5 启动过程的处理

MMC整流器启动时所有电容电压为零，这时候如果直接跑MPC，桥臂电流瞬间冲击很大，容易出现仿真崩溃。我的处理方式是分两步启动：

第一步，预充电阶段。给直流侧接一个限流电阻和一个二极管整流桥，让电容先充到接近额定值。这个阶段不启用MPC，功率开关全部关断，靠二极管续流充电。

第二步，MPC切入。当直流电压达到额定值的80%以上后，切换到MPC控制，把外环的电流参考幅值从0慢慢爬升到额定值，斜升时间设成0.5秒，观察波形是否平顺。这个软启动操作在很多论文里不会写，但实际操作时少不了的，不加这个过程模型经常第一个仿真周期就崩掉。

5. 仿真结果分析与论文复现对比

5.1 稳态波形特征

模型跑通后，稳态波形的三条指标是最直观的验证。第一，交流侧电流io要呈现标准正弦波，且和电源电压同相位，功率因数接近1。第二，直流侧电压Vdc稳定在200V附近，纹波很小。第三，子模块电容电压Vc均围绕50V左右波动，波动范围控制在5%以内。桥臂电流中的环流分量被压到很小，波形上看不到明显的二倍频包络。

我把这些波形在Scope里存成了时间序列数据，然后导到MATLAB里做FFT分析。交流电流的THD在稳态下约为2.5%-4%，比传统PI控制下的5%-6%明显改善，这也符合FCS-MPC论文里的典型结论。

5.2 动态响应性能

动态响应测试我做了两组：负载突增和直流电压参考跳变。

负载突增测试：在t=1.0s时把负载电阻从20欧跳到10欧，观察直流电压跌落和恢复时间。FCS-MPC的响应很快，直流电压最多跌到约190V，0.2秒内恢复参考值，交流电流幅值立刻增大，整个过程没有明显超调震荡。

直流电压参考跳变测试：把Vdc_ref从200V阶跃到240V，观察电压跟随速度。因为外环是PI控制、内环MPC响应极快，总体电压恢复时间约0.1秒，电流波形在调节过程中畸变程度很小，动态响应指标优于传统调制方案。

5.3 环流与电容电压波动分析

环流抑制效果要看桥臂电流波形和频谱。未加环流抑制时，桥臂电流里有明显的100Hz分量，幅值可能达到基波幅值的一半。加入混合控制里的环流抑制闭环后，100Hz分量降到基波的5%以下，桥臂电流波形显著改善。

电容电压波动是另一个重点。N=4、C=4.7mF的参数下，稳态电容电压波动范围约±3%，比论文中可能给出的±5%略好，这是因为我把权重系数λ3调到了0.08左右，对电压均衡的惩罚力度更强。代价是开关频率略有上升，这也是权重博弈的直观体现。

5.4 与论文原文结果对照的注意事项

复现论文时最容易遇到的问题就是：参数完全按论文设置，结果却对不上。我总结了三个原因。

第一是标幺化基准不同。很多IEEE论文用标幺值系统，各变量基值定义不写明，导致实际数值差异很大。复现前先把论文的基值体系搞清楚，把标幺值换算回实际值再进Simulink。

第二是采样频率差异。论文往往只在图表注里写控制频率，不写清楚采样周期，直接套用参数很容易偏差。如果论文给了实验样机的DSP控制频率，就按那个来；如果只给了仿真，那就按常见的10kHz或20kHz尝试。

第三是负载工况不同。有的论文仿真用的是阻性负载，有的是阻感性负载，启动过程、稳态电容电压波动结论会差很多。对照论文结果时，先把负载条件和初始条件对齐。

复现的价值不在于逐像素还原论文的图，而在于把控制算法吃透验证。如果你的结果在趋势上、数量级上、动态特征上都与论文一致，那这个复现就是成功的。

6. 常见问题与排查技巧实录

6.1 问题速查表

下面这份速查表是我调试过程中整理出来的，几乎涵盖了常见故障：

6.2 权重系数的快速调节思路

权重系数调不出来是MPC最常见的卡点，这里给一个可操作的流程。先把输出电流项权重固定为1，环流和电容电压权重都设为零，跑通电流波形。这一步的目标是让整体控制闭环能稳定运行，对应一个所有子模块都参与输出的基准状态。

在这个基础上，加环流权重，从0.1开始，每次翻一倍观察桥臂电流的100Hz分量幅值，直到降到可接受范围。如果引入环流权重后电流波形出现明显畸变，说明环流项和电流项冲突，此时不继续加λ2，而是增大桥臂电感或者加入独立环流抑制闭环来分担任务。

最后加电容电压均衡权重，从0.01开始，观察电容电压波动范围。波动范围太大就增加λ3，但注意λ3过大会导致开关频率过高、开关损耗暴增。如果开关频率超过器件允许范围，说明均衡项惩罚太重，需要反过来降低，并通过增大电容容值解决波动问题。

6.3 仿真速度优化的三个技巧

MPC类模型跑起来慢是通病，我实际测试下来有三个有效优化手段。

第一个是减少输出诊断输出数据。Simscape默认记录大量中间节点波形，控制Scope一多了仿真速度急剧下降。我的做法是只记录关键信号，把不需要的To Workspace模块全部取消。

第二个是合理设定控制采样周期。Ts从50us增加到100us，计算量减半，控制性能下降不大。如果你的论文目标是验证控制架构而非极致波形质量，100us完全够用。

第三个是把MATLAB Function内部循环向量化。枚举70种状态的时候，用手写for循环嵌套在函数里会让仿真时间翻倍。改成先构建状态组合矩阵，再统一用矩阵运算算预测值，速度提升非常明显。这也是代码生成的推荐写法，对后续部署到控制器也有利。

6.4 复现过程中最值得注意的一个细节

最后说一个我自己吃了大亏的细节：开关信号和仿真采样时刻的对齐。Simulink里ZOH采样时刻和MATLAB Function输出更新时刻如果不在同一拍上，会产生一个额外的小延迟，这个延迟在高速动态下会造成电流跟踪抖动。解决办法是把ZOH的采样时间设置为Ts的整数倍，并且确保触发信号由同一个时钟源提供。我用的是一个周期为Ts的脉冲发生器作为整个控制系统的时基参考，所有ZOH、触发子系统都挂在同一个脉冲信号上。从此再也没有出现过诡异的高频抖动。

7. 写在最后：复现SCI论文的一点心得

论文复现这件事，看起来是照着公式搭模型，真正做起来才知道细节量有多大。我个人最大的感受是：不要试图一步到位，一定要为每个环节留出单独验证的时间，尤其是FCS-MPC这样的控制算法。预测模型、代价函数、枚举逻辑、延迟补偿、环流闭环，任何一个环节出问题，整个仿真都会表现得极不稳定，而且错误现象往往不直观，排查起来很耗时间。

顺着这次项目的经验，后续扩展方向其实很明确。把单相模型扩展成三相MMC整流器，验证三相不平衡工况下电流控制效果；或者在FCS-MPC框架里引入模型观测器，实现对非理想电源电压的自适应控制；再进一步，把仿真模型通过Embedded Coder自动生成C代码，部署到DSP或FPGA上做硬件在环验证，这才是这条技术路线真正落地的完整闭环。

如果你也在复现这类论文，记住这句话：先让模型跑起来，再去优化指标。一个能跑的粗糙模型，价值远大于一个永远调不通的完美设计。