三相整流器FCS-MPC模型预测控制Simulink仿真全解析

我在一个整流器预研项目上，最初还是老老实实用电压外环PI加电流内环PI、再配SVPWM调制的老路子。电流内环要解耦，电压外环带宽要和内环拉开，再加上电网电感L估计不准，整定参数花了两天，结果负载突变时电流还是超出了预期。后来把内环换成有限控制集模型预测控制（FCS-MPC），事情一下子简单了不少：不需要调制器，不需要内环解耦，每次控制周期直接选出最优开关状态输出，MATLAB/Simulink里一个自定义函数就能写完核心逻辑。

这篇博客就把我搭建的"三相整流器基于模型预测的仿真模型"完整拆开讲一遍。内容包括三相整流器在dq旋转坐标系下的离散数学模型、Simulink里的分层搭建方法、8个电压矢量的滚动优化代码、代价函数权重怎么调，以及我在实际调试过程中踩过的那些坑。希望给正在做三相PWM整流器仿真、或者想把MPC从理论落到Simulink模型里的朋友一个可以直接参考的完整路径。

1. 为什么我会用模型预测控制来做三相整流器，而不是继续堆PI参数

1.1 传统双闭环在"不调参"这件事上的痛苦

三相PWM整流器的常规做法是：外环稳住直流电压，内环控制dq轴电流，再用SVPWM或SPWM生成驱动脉冲。这套方案性能并不差，它在大量工业产品里经过验证，稳态精度和开关频率固定的特性都很好。但问题在于，每次换一套主电路参数，内环PI就得重新算一遍。

电流内环通常在dq坐标系下工作，d轴和q轴之间存在耦合项，理论教科书会告诉你用前馈解耦。但真正仿真的时候，电感L的标称值和实际值总有一点出入，解耦误差一变，PI的稳定裕度也跟着变。更麻烦的是，如果后面还要接不同短路比的电网，L从3mH变成8mH，原来那组PI参数就不再是最优，甚至可能出现振荡。

我在某供电系统预研项目里就是被这个反复调参的过程拖住了进度。仿真平台本身没问题，问题出在"每个新工况都要重新整定一套参数"这件事上。

1.2 有限控制集模型预测的基本套路

模型预测控制的核心想法很简单：先预测，再决策。

三相两电平整流桥只有6个功率管，上下桥臂不能同时导通，所以有效开关状态一共是8种组合，这8种组合就是所谓的"有限控制集"。控制算法在每个采样周期做这样几件事：

1. 读取当前采样值：电网电压、电感电流、直流电压。
2. 利用三相整流器的数学模型，对8个开关状态分别预测下一时刻的电流值。
3. 把预测电流和参考电流放进代价函数，误差最小的那个开关状态就是本周期要输出的结果。
4. 直接把对应的门极脉冲送给功率管，等下一个采样周期重复。

整个过程不需要三角载波，不需要调制器，也不需要电流内环PI。我最早看到这个思路的时候感觉像是"暴力枚举"，但跑起来之后发现，它对参数变化的适应能力反而比传统方法更直接，因为预测模型里的L和R是显式参与的，主电路一改，改两个常数就行，不需要重新凑PI参数。

1.3 它到底省掉了什么，又带来了什么

省掉的东西很明确：调制器、内环PI、电流解耦环节。

但代价也不是没有。模型预测控制对"模型精度"的依赖比对PI参数更显性，如果L参数估得不准，预测电流就会系统性偏差，结果可能表现为电流纹波偏大或者稳态误差去不掉。另外，它的开关频率不固定，不像SPWM那样由载波频率锁定，这点在工程上会带来滤波器设计上的麻烦。

所以，我的态度是：MPC不是来"取代"传统方法的，它更适合动态响应要求高、或者你不愿意花大量时间调内环参数的场景。工具选型这件事，先想清楚需求，再决定用哪套方案。

2. 从三相整流器数学方程到步长离散，仿真前必须想清楚的事

2.1 整流器在dq坐标系下的连续模型

仿真不是直接摆一堆Simulink模块就完事，先把数学模型写对，后面写代码才不会返工。

三相电网电压经过电感和等效电阻进入整流桥，整流桥交流侧可以看成由开关状态决定的受控电压源。直流侧是电容和负载并联。开关函数定义成：

• Sa = 1表示A相上桥臂导通，Sa = 0表示A相下桥臂导通；
• Sb、Sc同理。

整流桥交流侧相电压和直流电压Vdc的关系是：

va = Vdc * (2*Sa - Sb - Sc) / 3 vb = Vdc * (-Sa + 2*Sb - Sc) / 3 vc = Vdc * (-Sa - Sb + 2*Sc) / 3

这组式子非常关键，后面所有预测都是基于它。如果你不自己算，查表也行，8个开关状态对应8组va/vb/vc。

然后进入dq旋转坐标系。设d轴与电网电压矢量同方向，三相整流器的电压方程可以写成：

L * did/dt = ed - R*id - vd + ω*L*iq L * diq/dt = eq - R*iq - vq - ω*L*id

这里的ed、eq是电网电压的dq分量，vd、vq是整流桥交流侧电压的dq分量，ω是电网角频率。

2.2 一阶欧拉离散，直接把预测式写出来

模型预测控制在数字控制器里必须使用离散模型。工程上最常用的一阶欧拉法，形式简单，采样频率足够高时精度完全够用。

对电流微分做近似：

did/dt ≈ (id(k+1) - id(k)) / Ts diq/dt ≈ (iq(k+1) - iq(k)) / Ts

代入连续方程，整理出预测式：

id(k+1) = id(k) + (Ts/L) * [ed(k) - R*id(k) - vd(k) + ω*L*iq(k)] iq(k+1) = iq(k) + (Ts/L) * [eq(k) - R*iq(k) - vq(k) - ω*L*id(k)]

这组式子就是滚动优化的核心。每次采样后，把8组开关状态分别代入vd/vq，算出8个预测电流值，然后选代价函数最小的那一组。

这里有一个非常容易搞错的地方：vd/vq的正负号定义。如果你用Simulink自带的abc到dq0变换模块，变换方向可能和我上面写的相反，那会出现电流预测总是不对的现象。遇到这种情况不要硬改代码，先确认你的Park变换方向和电压方程方向是自洽的。

2.3 8个开关状态对应的交流侧电压表

为了减少仿真运行时反复做三角函数转换，我习惯把8个开关状态的vd/vq离线算好，写成常量表，在代码里直接查表。假设Vdc为当前直流电压，θ为PLL角度：

严格算完这一整张表需要按实际变换系数展开，我在仿真代码里是直接调用一个v_lookup(Vdc, theta, state_index)函数来算的。因为Vdc在启动过程会从0爬到600V，每次算一次也很快，8次循环对20kHz采样率来说开销完全可以接受。如果是想在更低成本的控制器上落地，再优化成离线表也不迟。

3. 按我的习惯，Simulink里分四层把MPC整流器搭起来

3.1 主电路和前级测量：不要把电压电流量测直接怼进MPC

Simulink里搭主电路我直接用Simscape Electrical的元件库，三相交流电源、三相RL阻抗支路、Universal Bridge、直流电容、直流负载，这些都有现成模块。三相两电平桥我选的是IGBT/Diode结构，因为整流工况下反并联二极管在启动阶段还要走不控整流，这个细节很关键。

搭建顺序：

1. 三相电源设置为线电压380V、频率50Hz。
2. 三相RL支路参数：L=5mH、R=0.05Ω，R把线路电阻和电感寄生电阻都算进去了。
3. Universal Bridge选择三相臂、IGBT/Diode，直流侧接电容2200μF。
4. 直流负载先用10Ω，后面做负载突变时可以换成阶跃投切电阻。

电压电流测量不要直接接进MATLAB Function，中间一定要通过信号转换和零阶保持器。Simscape物理信号和Simulink普通信号之间需要转换模块，这个看似不起眼，实际上很多人第一次跑模型报"数据类型不匹配"就是卡在这里。

我习惯这样处理：电流测量模块输出物理信号，加一个Simulink PS-Simulink Converter，转换成普通信号，再过零阶保持器，采样周期设成Ts=50μs，也就是20kHz的控制频率。这样做的好处是，MPC里面每一拍拿到的输入信号都是同一个采样时刻的快照，而不是连续斜坡信号。

3.2 坐标变换与PLL：锁不到相角后面全是白搭

MPC预测用的ed/eq和id/iq都在旋转坐标系下，所以必须先有准确的电网电压相角θ。我用的是Simscape电气库里的三相PLL模块，输入三相电压，输出角度θ。

有一点要注意：测量出来的三相电压是经过RL支路之后的整流桥输入电压，不是电源电压本身。如果用这个电压去锁相，得到的角度里会带一点电感压降造成的相位偏移。不过对于仿真入门来说，影响不大；如果做高精度控制，一般建议直接从电网母线侧取电压信号参与锁相。

θ出来后，用库里的abc_to_dq0变换模块，分别变换三相电网电压和三相输入电流，得到ed/eq和id/iq。dq变换系数选幅值不变形式，和我的预测模型保持一致。

3.3 MATLAB Function里写8个电压矢量的滚动优化

核心控制算法我放在一个MATLAB Function里，输入是ed、eq、id、iq、id_ref、iq_ref、Vdc、theta，输出是6路门极脉冲。

先给一个可以直接用的简化版本：

function gates = fcs_mpc_rectifier(ed, eq, id, iq, id_ref, iq_ref, Vdc, theta) % 三相两电平整流器 有限控制集模型预测电流控制 % 输入：电网电压dq分量、电流dq分量、电流参考值、直流电压、PLL角度 % 输出：6路门极脉冲 [Sa Sb Sc ~Sa ~Sb ~Sc] % 参数，可以在函数里直接定义，也可以作为持久变量从外部传入 L = 5e-3; R = 0.05; Ts = 50e-6; w = 2*pi*50; % 8个开关状态 states = [0 0 0; 0 0 1; 0 1 0; 0 1 1; 1 0 0; 1 0 1; 1 1 0; 1 1 1]; % 当前PLL角度 ct = cos(theta); st = sin(theta); % 初始化最优值 J_min = inf; best_idx = 1; for idx = 1:8 Sa = states(idx,1); Sb = states(idx,2); Sc = states(idx,3); % 整流桥交流侧相电压 va = Vdc * (2*Sa - Sb - Sc) / 3; vb = Vdc * (-Sa + 2*Sb - Sc) / 3; vc = Vdc * (-Sa - Sb + 2*Sc) / 3; % abc -> dq vd = 2/3 * (va*ct + vb*cos(theta-2*pi/3) + vc*cos(theta+2*pi/3)); vq = -2/3 * (va*st + vb*sin(theta-2*pi/3) + vc*sin(theta+2*pi/3)); % 预测下一拍电流 id_p = id + (Ts/L) * (ed - R*id - vd + w*L*iq); iq_p = iq + (Ts/L) * (eq - R*iq - vq - w*L*id); % 代价函数 J = (id_ref - id_p)^2 + (iq_ref - iq_p)^2; if J < J_min J_min = J; best_idx = idx; end end % 根据最优状态输出门极脉冲 Sa = states(best_idx,1); Sb = states(best_idx,2); Sc = states(best_idx,3); gates = [Sa; Sb; Sc; 1-Sa; 1-Sb; 1-Sc]; end

这段代码是用来跑通的版本，不是最终优化版本。后面提到的延迟补偿、开关惩罚项，都在这个框架上加。

3.4 电压外环PI的输出怎么限幅才能喂给内环

直流电压要稳住，必须有一个外环给出id_ref。这个外环我用的还是PI，因为MPC内部管的是电流动态，电压外环用PI完全没有问题。

电压PI的作用是：Vdc偏高就减小id_ref，Vdc偏低就增大id_ref。但如果id_ref不限幅，启动瞬间电容充电电流会非常吓人，可能烧管子或者让仿真发散。

我把PI输出限幅设在±40A。40A这个数值是结合启动电流峰值估计出来的，你可以按你的器件额定电流调整。积分环节也要限幅，不然限幅之后积分继续累加，会形成严重的超调。我用的初值是Kp=0.6、Ki=20，这个值不一定适用于所有系统，以我给的参数跑通之后再按阶跃响应微调即可。

iq_ref直接给0，意思是控制整流器工作在单位功率因数状态，也就是网侧电流只含功分量。如果你还想做无功补偿，iq_ref可以作为无功控制的输出通道。

4. 代价函数不是"公式抄完就完事"：权重系数与开关切换惩罚

4.1 最基础的代价函数长什么样

我上面给的代码里，代价函数是：

J = (id_ref - id_p)^2 + (iq_ref - iq_p)^2

两项都是电流误差平方，单位一致，所以默认权重都是1。这个形式作为入门完全够用，它让算法在8个开关状态里选一个让预测电流最接近参考值的状态。

但运行一段时间你会发现一个问题：算法为了把电流误差压到极小，会非常频繁地切换开关状态。因为相邻采样点参考电流变化一点，最优状态可能就从100跳到000，下一拍又跳回来。开关频率变高，电流纹波也没见得好多少，功率器件的开关损耗倒是上去了。

4.2 要不要加开关惩罚项，加多大合适

解决这个问题最直接的办法，是在代价函数里加一项开关切换惩罚项。

定义上一拍输出状态为Sa_old、Sb_old、Sc_old，当前候选状态为Sa、Sb、Sc，三者逐位比较，有几个位不同就记几次切换。用汉明距离来表示：

n_sw = abs(Sa - Sa_old) + abs(Sb - Sb_old) + abs(Sc - Sc_old)

代价函数变成：

J = (id_ref - id_p)^2 + (iq_ref - iq_p)^2 + λ_sw * n_sw

λ_sw这个系数就看你对开关频率的忍耐程度。它的物理含义可以理解为：多切换一个功率管，相当于允许引入多大的电流平方误差。假设你允许因为少切换一次而损失0.5A的跟踪精度，那电流误差项的差值是0.25 A²，λ_sw就可以从0.25附近开始试。

我自己的调整流程：

1. 先设λ_sw=0，跑一遍，记录平均开关频率和电流纹波。
2. 把λ_sw加到0.1，再看开关频率有没有下降。
3. 逐步增大到0.3左右，直到开关频率明显下降、但电流THD还在可接受范围。

每个系统的量纲一致，这个流程是通用的。

4.3 权重调整里容易被忽略的一个细节

有一点必须提醒，代价函数里的电流误差项是"预测值"而不是"实际值"的误差。模型偏差越大，代价函数本身的可信度越低。所以当你发现怎么调λ_sw都不对劲的时候，先回头确认离散模型参数有没有写错，而不是死磕权重。

另外，如果你后续想直接控制直流电压，代价函数里加Vdc项，电流项和电压项的量纲就不一致了。那时候不能简单把权重都设为1，而要先做归一化处理，或者用标幺值系统。这个我在后面扩展部分会再提一下。

5. 仿真结果跟预期一致吗：启动、稳态、动态和THD

5.1 启动过程和直流电压建立

用我上面给的参数启动仿真，直流电压不会一开始就受控，因为启动瞬间直流电容没电，MPC的输出在低电压下工作，桥臂基本走的是反并联二极管的不控整流路径。Vdc从0V先自然爬升到约510V左右，也就是380V线电压经不控整流得到的峰值电压，等Vdc超过这个值之后，Boost特性才起作用，MPC开始通过控制电流把电压拉向600V。

我仿真的结果是大约0.25s进入稳态，启动过程直流电压超调在5%左右，没有振荡。这个启动过程取决于电压外环PI参数和电流限幅值，如果你想更快建立电压，可以适当提高id_ref限幅，但代价是启动电流峰值变大。

5.2 负载突变时MPC的动态响应

整流器项目最关心的往往是动态性能。我把负载电阻从10Ω切换成5Ω，相当于输出功率从36kW变成72kW（按Vdc=600V估算）。

由于电压外环PI响应需要时间，直流电压会先跌落，然后id_ref增加，MPC在1到2个控制周期内就能把实际电流跟踪到新的参考值附近。整个恢复时间大概在20到30ms，具体取决于外环带宽。这个结果说明一个关键问题：MPC的内环动态很快，但系统的整体动态响应仍然受电压外环PI限制。

如果你希望直流电压跌得更少，可以加大电压环PI的Kp，但要注意不能太激进，否则启动超调和负载阶跃都会出现电压振铃。

5.3 网侧电流波形、单位功率因数和THD

稳态下看三相电流波形，接近正弦，但会比SPWM方案多一些“毛刺”，这是MPC开关频率不固定带来的高频分量。

把网侧A相电压和A相电流放到同一个示波器里，会发现电流相位和电压基本一致，功率因数接近1。对电流波形做FFT分析，我的仿真模型THD大约在3%到5%之间，主要谐波呈现为宽频分布而不是集中在某几次谐波。这是因为MPC的开关频率在频谱上较为分散。

5.4 和传统PI加SVPWM方案的对照

为了确认MPC到底值不值得用，我在同一套主电路上又跑了一版传统双闭环加SVPWM做对照，结果大体上：

我个人在项目里的感受是：如果你在做一个参数经常变化、动态要求又高的平台，MPC的优势非常明显；如果工况固定、追求极限稳态波形质量，那传统的调制方法仍然值得保留。

6. 这一路调试下来，我踩掉的那些坑

6.1 预测符号写反导致电流"倒灌"

最开始跑模型时，直流电压建立不起来，电流越来越大，像是整流器变成了逆变器，能量从直流侧倒灌回电网。查了半天，最后发现是dq变换正负号的问题。

这个问题特别隐蔽，因为不报错，波形看起来只是"方向不对"。我的排查方法是：在仿真里只给一组固定的开关状态，手动计算预测电流方向，和公式推导对照。一旦发现符号反了，电流偏差会越来越大，解决办法就是统一坐标变换方向。

拿到验证过的代码之后，我把这组基础模型封存起来，后面再搭维也纳整流器或者其它变换器，都以这版方向定义为基准。

6.2 Simulink数据类型和零阶保持器的细节

Simulink里MATLAB Function的输入输出默认是double类型。但Universal Bridge的门极信号有时候会被识别成Boolean，直接把0和1的double数组接进去，仿真会报类型不匹配。我的做法是在函数输出后面加一个DataType Conversion模块，再进Universal Bridge。

另外一个容易被忽略的问题是零阶保持器。如果不加零阶保持器，MATLAB Function的输入信号在仿真步进内部会被看成连续变化，导致同一个控制周期内多次执行函数，结果看似能用，但实际上控制时序和真实硬件对不上。所有进入控制函数的信号都过一遍采样器，采样时间设成Ts，模型行为才和真实控制器一致。

6.3 采样频率和求解器设置

MPC的预测模型使用一阶欧拉离散，Ts越小越准，但计算量越大。我在三相两电平整流器上用Ts=50μs跑20kHz，效果还不错。如果Ts大于100μs，一阶欧拉法的预测误差会明显增加，电流纹波变大，甚至可能不稳定。

Simulink求解器我建议用固定步长离散求解器，步长不超过10μs，也就是Ts的1/5。如果步长太大，PLL角度更新跟不上采样时序，开关状态和坐标变换会错位。我一开始用可变步长求解器，仿出了很多奇怪的高频毛刺，改成固定步长后问题消失。

6.4 延迟补偿：预测一拍比不预测强太多

实际控制系统从ADC采样到PWM输出，中间一定存在计算延迟。MPC在这个周期计算出的最优开关状态，最早也要到下一个周期才能施加到功率管上。如果你完全忽略这个延迟，在低采样频率下稳定性会受到很大影响。

补偿方法很简单：在MATLAB Function里做两步预测。先根据当前采样值预测出id(k+1)、iq(k+1)，再把id(k+1)、iq(k+1)当作当前值，代入同一个预测公式，算出id(k+2)、iq(k+2)。代价函数使用的就是k+2时刻的预测电流。

代码改动量不大，实际效果非常明显。做了两步预测之后，同样一组主电路参数下，允许的最大采样周期可以放宽不少，波形也更加干净。

7. 我最后留下的几句话和后续扩展方向

做三相整流器的MPC仿真，真正难的不是Simulink建模，也不是写那8行循环，而是把数学模型的每一个符号都验证到位。我在这套模型上踩过符号反号、零阶保持器、采样时序这些坑之后，最大的体悟是：模型预测控制把传统PI调参的痛苦转移到了"模型准确性"和"代价函数设计"上，你少写了一堆解耦代码，但必须比谁都清楚自己的离散方程。

如果你在这套模型基础上继续扩展，可以往这几个方向走。第一，把电流内环换成直接功率控制，代价函数里同时包含有功和无功误差，启动和动态响应会更快。第二，加入权重自适应，在线调整λ_sw，让稳态和动态都保持比较好表现。第三，加LCL滤波器改善网侧电流波形，然后做谐振抑制。第四，把三相两电平桥换成三电平NPC结构，MPC的控制集从8个变成27个，原理类似，但代价函数设计和开关状态表会复杂不少。

我自己的习惯是，每次改主电路参数，先跑一遍不控整流阶段，确认直流电压自然爬升正常，再切入MPC控制。这一步能帮你把主电路问题和算法问题分开，排查速度快很多。希望这篇博客能让你在三相整流器MPC仿真这条路上少走一点弯路。