Simulink下永磁同步电机四种模型预测控制策略对比与实现

做永磁同步电机驱动控制仿真的同行应该都有过这种体验：传统双闭环PI加SVPWM调出来的波形成熟稳定，但要应对参数整定和动态响应之间的拉扯；一旦切换到有限集模型预测控制（FCS-MPC），思路就变了——直接把开关状态当成决策变量，用预测模型比较未来几个周期的电流轨迹。我这次在Simulink里搭了一套覆盖单矢量、占空比优化、双矢量和三矢量四种策略的PMSM模型预测控制仿真模型，同一台电机、同一套采样条件下把四种算法的性能差异全部拉出来对比，顺便把推导、代码和调参过程一起记录下来。

这套仿真模型适合正在研究永磁同步电机模型预测控制算法对比、准备写论文或者做课程设计的同行参考，也适合只听说过FCS-MPC但想亲手复现基础版本的人。全文不打算堆公式，而是讲清楚四种矢量策略到底差在哪里，以及如何在Simulink里把它们实现并得到可解释的仿真曲线。我会把每个模块的接线方式、核心MATLAB Function代码、采样周期怎么选、权重因子怎么调，都按我实际搭建时的顺序写出来，尽量让读者照着就能搭出一套可切换的PMSM-MPC对比平台。

1. 项目概览与控制问题设定

1.1 为什么把四种MPC放进同一个模型

我最初接触FCS-MPC时，觉得它和传统调制最大的区别，就是把逆变器的离散开关状态直接纳入优化。很多人上手先做的是单矢量方案：枚举几个开关组合，选一个代价最小的直接输出一整个控制周期。这个方案能跑通，但电流纹波和转矩脉动并不好看。于是后续的占空比优化、双矢量、三矢量，本质上都是在回答同一个问题：一个控制周期内，到底让哪些矢量作用、各自作用多久，才能让电流轨迹离参考值更近。

把它们放进同一个模型的原因很实际：不同策略的对比必须在完全相同的工作条件下才有说服力。如果电机参数、采样周期、坐标变换实现方式各不相同，对比出来的差异就容易混入实现误差。我在模型里留了一个算法切换端口，四种策略共用同一套电流采样、同一套坐标变换、同一套PMSM模型，切换时只改变控制器内部的矢量生成逻辑。

1.2 控制对象与模型参数

我选用的是表贴式永磁同步电机，额定参数如下表。这个电机模型比较常见，电磁时间常数也适合做50μs级的预测控制。

选择表贴式电机的好处是Ld=Lq，转矩方程里磁阻项为零，控制思路简单：d轴电流给0，q轴电流出转矩，这样四种MPC策略的差异最容易解释清楚。如果你手头是内置式电机，只需要把Ld和Lq改成不等值，把磁阻转矩项补上，算法代码里的通式已经保留了(Ld-Lq)项，直接替换参数即可。

1.3 控制目标与技术路线

整套系统控制目标明确：外环转速PI生成q轴电流参考，内环由MPC控制器直接输出功率管开关信号，d轴电流参考给0。仿真中测试两类工况：稳态匀速运行时看电流质量；转速阶跃或负载突变时看动态响应。

控制器周期取50μs，对应等效开关频率20kHz。逆变器用理想IGBT桥，并加入死区时间2μs。实际工程中如果要做代码生成，这个建模方式可以很自然地变成嵌入式软件里的MPC函数。顶层信号流向为：直流源→三相逆变桥→PMSM→电流/转速测量→坐标变换→速度PI→MPC控制器→门信号。

2. 四种矢量策略的原理与差异

2.1 从预测模型到代价函数

MPC的预测模型来自PMSM在dq同步旋转坐标系下的电压方程：

u_d = R_s * i_d + L_d * (di_d/dt) - ω_e * L_q * i_q u_q = R_s * i_q + L_q * (di_q/dt) + ω_e * (L_d * i_d + ψ_f)

用前向欧拉法离散化，取控制周期Ts，得到电流预测式：

i_d(k+1) = (1 - R_s*Ts/L_d) * i_d(k) + (Ts/L_d) * u_d(k) + ω_e*L_q/L_d * Ts * i_q(k) i_q(k+1) = (1 - R_s*Ts/L_q) * i_q(k) + (Ts/L_q) * u_q(k) - ω_e*Ts*(L_d/L_q * i_d(k) + ψ_f/L_q)

这套公式里假设了Ts远小于电气时间常数，50μs对大多数中小功率PMSM都成立。代价函数只取电流误差平方：

J = (i_d_ref - i_d(k+1))² + (i_q_ref - i_q(k+1))²

这套代价函数里没有任何权重因子。很多文章一上来就列各种加权多目标代价函数，反而掩盖了最核心的机制。先只用电流误差平方做对比，四种策略的差异一样能体现出来。

2.2 单矢量FCS-MPC的枚举逻辑

两电平逆变器有8种开关组合，对应7个空间电压矢量。单矢量策略在每一个控制周期枚举全部8个组合，对每种组合预测下一拍电流，把使J最小的那一组开关状态直接输出。

优点：实现最简单、计算量最小。缺点：一个周期只作用一个矢量，相当于只能用离散的七个方向去逼近任意参考电压，天然存在量化误差，所以电流纹波和转矩脉动偏大。打个比方，单矢量像朝目标方向卖力迈一大步，方向对了但落脚点很粗糙，剩下的偏差只能等下一拍再纠正。

2.3 占空比优化：给有效矢量补一段零矢量

既然单矢量是因为“一整周期都作用同一个方向”才产生过冲，那自然想到，让最优有效矢量只作用d*Ts的时间，剩余时间切换到零矢量。这里的d就是占空比。

设f0为零矢量作用下的电流变化率，f1为最优有效矢量作用下的电流变化率。终点电流为：

i(k+1) = i(k) + f0*Ts + d*(f1 - f0)*Ts

定义A为以零矢量为基准的误差向量，B为有效矢量相对零矢量的增量：

A_d = i_d_ref - i_d(k) - f0_d*Ts A_q = i_q_ref - i_q(k) - f0_q*Ts B_d = (f1_d - f0_d)*Ts B_q = (f1_q - f0_q)*Ts

代价函数J(d)=(A_d - dB_d)²+(A_q - dB_q)²，对d求导并令导数为零，得到解析解：

d = (A_d*B_d + A_q*B_q) / (B_d² + B_q²)

限幅在[0,1]之间即可。占空比方案不需要搜索d，直接算，计算量增加很小，但稳态电流质量提升明显。需要注意，零矢量在dq坐标下的电压为0，作用时的电流变化率由反电动势和电阻压降决定，不能简单写成0。

2.4 双矢量MPC：一个周期两个有效矢量

双矢量和占空比的区别在于，第二段不再用零矢量，而是用另一个有效矢量。通常第一个矢量仍按单矢量枚举选出，第二个矢量根据第一个矢量的相邻关系确定，两个矢量的作用时间同样用最小化电流误差的解析公式求，把上一节的f0替换成第二个矢量的变化率f2即可。

这种做法的实用场景：当工作点接近逆变器电压极限时，零矢量会让电流下跌、有效矢量让电流上升，单独占空比方案容易出现时间占比过大或过小的情况；两个有效矢量交替作用能更平滑地合成中间方向电压，动态性能更好。代价是开关切换次数比占空比方案高。

2.5 三矢量MPC：伏秒平衡求三段时间

三矢量方案直接换了一个思路：先算出本周期需要的参考电压矢量，再用两个相邻有效矢量加一个零矢量把它合成出来。参考电压可以用无差拍思想计算：

u_d_ref = L_d/Ts * (i_d_ref - i_d(k)) + R_s*i_d(k) - ω_e*L_q*i_q(k) u_q_ref = L_q/Ts * (i_q_ref - i_q(k)) + R_s*i_q(k) + ω_e*(L_d*i_d(k) + ψ_f)

再逆Park变换得到αβ轴参考电压：

V_α_ref = u_d_ref * cosθ_e - u_q_ref * sinθ_e V_β_ref = u_d_ref * sinθ_e + u_q_ref * cosθ_e

根据Vref所在扇区确定两个相邻有效矢量V1、V2，用伏秒平衡解时间：

Vref * Ts = V1*t1 + V2*t2 + V0*t0 t1 + t2 + t0 = Ts

三个矢量各作用对应时间，等效平均电压等于Vref，电流跟踪性能最好。三矢量本质上把FCS-MPC和SVPWM在瞬时电压合成上统一了。性能好，但计算量大，一个周期内开关切换次数也多。

2.6 四种策略横向对比

不同文献对双矢量和三矢量的定义略有不一致，有的把有效加零矢量也称作双矢量。我建议交流时先讲清楚矢量成分和时间分配方式，而不是只报一个名字。按本文的实现，四种策略的差异可以概括为：

3. Simulink模型架构与核心模块搭建

3.1 模型顶层结构与接线顺序

整个模型按这个顺序搭：直流源、Universal Bridge、PMSM、负载转矩、电流与速度测量、坐标变换、速度PI、MPC控制器、门信号生成。逆变器设置为IGBT/Diode，死区时间2μs。PMSM模块用Simscape Electrical里的标准模型，机械端口接负载转矩Step。

电气侧用固定步长求解器，步长1μs，求解器ode3。MPC控制器所在子系统单独用离散采样。速度PI的采样周期和控制周期保持一致，取50μs。坐标变换模块用MATLAB Function实现Clarke和Park变换，电角度θ_e由PMSM模块输出的机械角度乘以极对数得到，电角速度ω_e同样来自PMSM模块。

3.2 单矢量MPC的MATLAB Function实现

核心函数输入是电流参考、当前dq电流、电角速度和电角度，输出是三路开关信号。下面是完整的单矢量MPC函数代码：

function [Sa, Sb, Sc] = mpc_single(id_ref, iq_ref, id_k, iq_k, omega_e, theta_e) % 模型参数 Ts = 50e-6; Rs = 2.875; Ld = 8.5e-3; Lq = 8.5e-3; psi_f = 0.175; Vdc = 311; % 8种开关状态 [Sa Sb Sc] S = [0 0 0; 0 0 1; 0 1 0; 0 1 1; 1 0 0; 1 0 1; 1 1 0; 1 1 1]; % 预计算各状态对应的alpha/beta轴电压 Vab = zeros(8, 2); for idx = 1:8 sa = S(idx,1); sb = S(idx,2); sc = S(idx,3); ua = Vdc/3 * (2*sa - sb - sc); ub = Vdc/3 * (-sa + 2*sb - sc); uc = Vdc/3 * (-sa - sb + 2*sc); Vab(idx,1) = 2/3 * (ua - 0.5*ub - 0.5*uc); Vab(idx,2) = 2/3 * (sqrt(3)/2*ub - sqrt(3)/2*uc); end J_best = 1e10; best = 1; for idx = 1:8 % alpha/beta -> dq ud = Vab(idx,1) * cos(theta_e) + Vab(idx,2) * sin(theta_e); uq = -Vab(idx,1) * sin(theta_e) + Vab(idx,2) * cos(theta_e); % 前向欧拉预测 id_p = (1 - Ts*Rs/Ld) * id_k + Ts/Ld * ud + omega_e * Lq / Ld * Ts * iq_k; iq_p = (1 - Ts*Rs/Lq) * iq_k + Ts/Lq * uq - omega_e * Ts * (Ld/Lq * id_k + psi_f/Lq); J = (id_ref - id_p)^2 + (iq_ref - iq_p)^2; if J < J_best J_best = J; best = idx; end end Sa = S(best,1); Sb = S(best,2); Sc = S(best,3);

这段代码可以直接放进Simulink的MATLAB Function模块。需要注意MATLAB Function的采样时间必须设为Ts，保证控制器按50μs触发一次。id_k和iq_k来自Park变换，theta_e是转子电角度，omega_e是电角速度。

3.3 占空比、双矢量、三矢量的输出结构

占空比方案多了零矢量时间，MPC模块不能只输出一组开关状态，而是输出一串矢量序列和作用时间。我的做法是：MPC输出三个候选开关状态state1、state2、state3，以及对应的作用时间t1、t2、t0；后级用一个门信号生成子系统把时间展开成三相IGBT门极脉冲。

占空比核心代码是在单矢量找到最优有效矢量后，再计算一段duty：

% 最优有效矢量在dq轴的投影为ud_opt、uq_opt f1_d = (ud_opt - Rs*id_k + omega_e*Lq*iq_k) / Ld; f1_q = (uq_opt - Rs*iq_k - omega_e*Ld*id_k - omega_e*psi_f) / Lq; f0_d = (-Rs*id_k + omega_e*Lq*iq_k) / Ld; f0_q = (-Rs*iq_k - omega_e*Ld*id_k - omega_e*psi_f) / Lq; A_d = id_ref - id_k - f0_d*Ts; A_q = iq_ref - iq_k - f0_q*Ts; B_d = (f1_d - f0_d)*Ts; B_q = (f1_q - f0_q)*Ts; duty = (A_d*B_d + A_q*B_q) / (B_d^2 + B_q^2); duty = min(max(duty, 0), 1);

双矢量的计算结构与占空比基本一样，区别只是把f0换成第二个矢量的变化率，时间上限变成整个Ts而不需要补零矢量。三矢量则直接用矩阵法解t1、t2、t0，例如在扇区I中，V1和V2分别取100和110状态对应的空间矢量，联立伏秒平衡方程解出三段时间。

3.4 延迟补偿与两拍预测

真实数字控制器里，从采样到计算开关状态再到更新，至少有一拍延迟。如果不补偿，预测模型假设k时刻施加的矢量实际在k+1才生效，电流会出现一个采样周期的滞后，直观表现是d轴电流和q轴电流耦合加深。

解决办法是两步预测：在k时刻先用上一周期选择的矢量预测出当前电流，再以这个预测值为起点枚举所有候选矢量预测下一拍电流，代价函数比较下一拍电流与参考值。在Simulink中实现需要在函数外部把上一周期选择的矢量用Unit Delay模块保存回来，作为函数输入之一。不要试图在MATLAB Function里用persistent变量记录，因为Simulink对函数内状态的处理容易引入代数环或初始化问题。

3.5 求解器与仿真参数设置

固定步长建议取1~2μs，控制周期50μs，仿真时长0.5s。电气模型步长过大时，IGBT开关瞬间的电压尖峰容易造成数值噪声，甚至把MPC本身应该有的性能差异盖掉。MPC的采样时刻用离散触发，不要用过零检测，否则开关动作瞬间容易产生额外触发，导致控制周期不稳定。

4. 参数整定与仿真结果深度对比

4.1 工况设定

我设置的对比工况是：0~0.1s先以单矢量方式启动，稳定到额定转速；0.1s后通过切换端口换到目标策略；0.2s时突加负载转矩5N·m，记录稳态电流THD、转矩脉动峰值差、平均开关次数和动态调节时间。这里故意让电机先在相同策略下启动，是为了保证四种策略的起始条件完全一致。

4.2 稳态性能对比结果

按上述参数跑出来，相对趋势如下表。数值会随电机参数和采样周期变化，但排序基本稳定。

三矢量的稳态优势在低速轻载工况下格外明显，因为参考电压幅值小，单矢量容易出现相邻矢量来回切换的振荡。占空比方案虽然只比单矢量多了一段零矢量，THD改善却能接近一半，性价比很高。

4.3 为什么三矢量的电流波形最干净

三矢量方案在一个控制周期内把平均电压精确合成到参考值，等价于高频SVPWM加上预测控制的动态校正。单矢量方案只能选最接近的方向，相当于离散矢量量化误差直接进入电流回路，THD自然高。占空比方案补零矢量只是改善了幅值量化，方向仍只有7个；双矢量方案用两个有效矢量扩展了可选方向，进一步改善；三矢量则基本消除了方向受限问题。

4.4 调参过程中容易被忽略的细节

第一，代价函数里如果加入开关频率罚项，权重因子的数量级必须先归一化，否则罚项很容易把电流跟踪项淹没，导致动态响应迟钝。第二，Ts不能为了追求性能随意减小，真实控制器的计算时间是硬约束；仿真里可以作弊用很小的Ts，移到实验板就可能跑不完。第三，速度PI外环与MPC内环是非线性组合，速度环带宽不要设太高，否则iq_ref的抖动会直接放大成电流噪声。

5. 仿真中容易踩的坑与常见误区

5.1 把机械角度当成电角度送进Park变换

很多初版模型跑出来的id、iq波形出现周期性畸变，频率和转速相关，十有八九是这个问题。PMSM建模里的θ_e是电角度，是机械角度乘以极对数。如果直接接机械角度，电流波形的基频全错，PI输出也会发散。检查方法：把θ_e和u_α、u_β一起画出来，确认电压矢量和转子磁链的相对位置正确。

5.2 固定步长太大导致数值振荡

Simscape里的IGBT桥开关瞬间如果仿真步长太大，电流波形会出现锯齿状高频噪声。控制周期50μs时，仿真步长我建议不要超过5μs，最好1~2μs。否则对比THD时四条曲线差别不明显，甚至因为数值噪声把算法本身的性能优势盖掉。这个问题常被误认为算法问题，实际是求解器问题。

5.3 占空比方案里的零矢量选择影响开关次数

两电平逆变器有000和111两个零矢量，在理想模型里它们对电流作用完全等价，但门极切换次数完全不同。如果从000切到111，三相开关全翻转，开关次数立刻增加。仿真做算法对比时，统一用000；如果要估算实验中的实际损耗，建议按最小切换次数策略选择零矢量。这个细节直接决定了报告里的“平均开关频率”数值是否可信。

5.4 无差拍参考电压在过调制区的处理

三矢量方案中，参考电压Vref一旦超出电压六边形，伏秒平衡解出的t1或t2可能出现负值，门信号时间会乱。我采用的处理方式：判断|Vref|是否超过内接圆半径Vdc/√3，超出则做等比例限幅，缩回圆内。为什么缩到圆内而不是六边形内，因为圆内是线性调制区，不会产生负时间，代码简单且稳定。

5.5 算法切换对比的公平性问题

如果从仿真一开始就运行三矢量，单矢量方案从启动时刻起就在不同开关状态下工作，不能算同一个稳态工作点对比。我建议这么安排：0~0.1s先跑单矢量，让电机稳定到目标转速；0.1s后通过切换端口换到目标策略，再稳定0.1s后开始记录数据。这样四种策略的对比，就是同一台电机、同一转速、同一负载下纯粹由算法差异导致的结果。

6. 从这套模型还能继续往哪里走

6.1 把单矢量到三矢量的演进当成理解MPC的主线

我自己跑完这套模型的体会是，单矢量到三矢量不只是简单堆矢量数量，而是从“枚举离散方向”一步一步过渡到“连续合成参考电压”。理解这条主线之后，再去看各种变体，比如多矢量MPC、无差拍预测控制、模型辅助调制，就都能对号入座了。仿真阶段最重要的不是一开始就追求最复杂的策略，而是先把预测模型、坐标变换、延迟补偿这几个底层做扎实。

6.2 可以继续深挖的工程化方向

如果要做实际样机验证，下一步可以做几件事：把电机模型参数换成实测辨识值；把MATLAB Function用C代码生成，估算目标控制器的实际执行时间；叠加载磁观测器和参数在线辨识，弥补MPC对模型参数敏感的弱点。还可以把三矢量方案和SVPWM做统一，在设计硬件时直接复用传统调制器的外围接口，降低移植成本。

6.3 一点个人经验

如果你现在刚开始搭，我的建议是先只跑通单矢量版本，把电流波形、转速波形和坐标变换的对应关系看明白，再花半天把占空比加上。大多数模型跑出异常波形，问题都出在延迟补偿或者角度处理上，而不在MPC算法本身。希望这篇文章能帮你在Simulink里少走几步弯路，也欢迎你实际复现后把跑出来的差异反馈给我，后面我可以续一篇带参数摄动和扰动观测器的扩展对比。