我前阵子做完了基于PI电流控制器的PMSM矢量控制仿真，模型用MATLAB/Simulink 2018b搭建，最后还整理了一份说明报告。这是一套最经典、也最值得反复打磨的电机控制仿真组合：永磁同步电机做磁场定向控制，电流环用PI控制器，后面接SVPWM和逆变器。很多刚入门的同事都会从这套模型开始，但真正把它从“能动”调到“能看、能交、能复现”，还是有不少门道的。这篇文章我会从模型架构、PI参数整定、Simulink搭建过程、常见坑和报告写法五个方面展开，把我的实际经验交代清楚，希望能让正在做类似项目的人少踩几步。
1. 模型设计思路与整体方案选型
1.1 为什么核心选择PI电流控制器
永磁同步电机用矢量控制，本质上是想把交流电机复杂的三相耦合关系，转换成两个独立控制的直流分量。转子是永磁体，定子电流会产生旋转磁场，如果不做解耦，转矩和磁链相互牵连，动态性能会很差。矢量控制把定子电流分解到转子d-q坐标系里，d轴管励磁，q轴管转矩，这样电流控制逻辑就和直流电机几乎一样。
电流环是这个解耦思路的落脚点，目前工程上最常用的就是PI控制器。我选择PI而不是滑模、模型预测或滞环控制，原因很实在。第一，调试路径清晰，物理意义直观，带宽、比例增益、积分增益基本都能和电流环的数学模型对应上。第二，后续要往实际控制器移植时，底层代码大概率还是PI结构，仿真里用PI可以提前验证采样、限幅、抗积分饱和等工程问题。第三，对一份说明报告来说，PI参数整定过程有据可查，评审和接手的人都能快速理解，不需要额外解释复杂控制理论。
1.2 矢量控制整体框架与模块分工
这套模型走的是典型双闭环结构，外环是转速环，内环是d/q轴电流环，最后经SVPWM驱动机侧逆变器。整体思路并不复杂：转速误差经过转速PI生成q轴电流给定；d轴电流给定默认设0；两个电流PI分别输出d轴和q轴电压；再做Park逆变换得到静止坐标系下的电压分量；SVPWM模块把电压矢量转换成三相桥臂开关信号；逆变器再把这些信号变成实际加到电机端的三相电压。
我在搭建模型时特意把模块拆开，没有堆成一个巨大的子系统。这样好处很多：调试时能单独看每一级信号，报告里也好截图解释。下面这张表基本可以概括每个模块的作用：
| 模块 | 主要作用 | 关键实现点 |
|---|---|---|
| PMSM电机模块 | 提供三相电流、转子位置与转速 | 正确填入电机参数和机械负载 |
| 三相逆变器 | 将直流母线电压转换为三相交流 | 与SVPWM驱动信号严格对应 |
| Clarke/Park变换 | 三相电流转到旋转坐标系 | 电角度必须乘以极对数 |
| 电流PI控制器 | 调节d/q轴电流快速跟随给定 | 限幅和抗积分饱和 |
| 速度PI控制器 | 给出q轴电流指令 | 输出限幅要合理 |
| SVPWM模块 | 生成三相桥臂开关信号 | 扇区判断和过调制保护 |
如果你只是做电流环验证，不接转速环也能跑，但加上转速外环后，整个系统才真正成为一套可以对比工程曲线的驱动模型。说明报告里如果只有电流波形，说服力会弱很多，所以我还是建议按双闭环来搭。
1.3 为什么定在MATLAB/Simulink 2018b
我用的环境是2018b，理由很现实。Simulink自带的电力系统模块和PMSM电机模块，在这个版本里已经足够搭建带功率级闭环的仿真，不额外买工具箱也能跑起来。另一个重要原因是兼容性：2018b生成的文件能被后续版本自动升级，但反过来，新版文件拿到旧版里基本上打不开。很多工业项目里的验证电脑还停留在2018b，因为代码生成工具链都验证过，用老版本反而稳妥。
如果你现在拿到这套模型，打算在更高版本里打开，要注意几个名称差异：部分电气模块库路径重排过，个别模块参数弹窗布局也有变化，但核心数学关系完全一样。真遇到缺失模块，搜索库里对应的PMSM、Universal Bridge关键字，基本都能找到替代。模型本身不是版本绑定的，绑定的只是操作习惯。
2. PMSM数学模型与PI参数整定要点
2.1 PMSM在d-q坐标系下的数学模型
建模之前，先把PMSM的数学模型写清楚。这个模型是所有后续控制器设计的地基，不能只用一个黑盒子跑仿真。同步旋转坐标系下的定子电压方程是：
Vd = Rs·id + Ld·d(id)/dt - ωe·Lq·iq
Vq = Rs·iq + Lq·d(iq)/dt + ωe·(Ld·id + ψf)
电磁转矩方程为：
Te = 1.5·p·[ψf·iq + (Ld - Lq)·id·iq]
再配上机械运动方程：
J·dωm/dt = Te - TL - B·ωm
对于表贴式永磁同步电机，Ld等于Lq，转矩方程会变成Te = 1.5·p·ψf·iq。这意味着id对转矩没有贡献，采用id=0控制最省电流，效率也最好。很多模型默认电机是表贴式，所以d轴电流给定设0完全合理。如果被控对象是内置式PMSM，Ld和Lq不相等，最大转矩电流比控制会更复杂，但电流环PI的整定思路不变。
2.2 电流环PI参数的推导与计算
电流环PI参数是整个模型里最不能拍脑袋乱填的部分。把PI控制器到电机定子绕组的电流通路近似成一个一阶惯性环节，被控对象就是G(s) = 1/(Ls + R)，其中L是d轴或q轴电感，R是定子电阻。控制器写为：
H(s) = Kp + Ki/s
让PI的零点去对消被控对象的极点，也就是令Ki/Kp = R/L。代入之后，开环传递函数变成了Kp/(L·s)，相当于一个纯积分环节。闭环系统是一个一阶低通滤波，带宽为ωc = Kp/L。因此可以直接得到：
Kp = L·ωc
Ki = R·ωc
这里的ωc是角频率，单位是rad/s，不是Hz。如果期望电流环带宽fc = 500Hz，那么ωc = 2π·500 ≈ 3140rad/s。假设电机参数是Rs = 0.2Ω，L = 0.8mH，那么Kp = 0.0008 × 3140 ≈ 2.51，Ki = 0.2 × 3140 ≈ 628。这个起始值在绝大多数仿真里都不会发散，下一步再根据阶跃响应微调。
需要注意的是，实际数字控制系统还有采样延迟和PWM延迟，电流环带宽不能设计得太高，一般取开关频率的1/10到1/20。例如10kHz开关频率，电流环带宽取500Hz到1000Hz比较合理。太高会让数字延迟影响变大，太低则动态响应变慢。这个上限在实践中非常重要，仿真里如果发现电流震荡，第一步不是继续加带宽，而是先降下来。
2.3 坐标变换与SVPWM的关键算式
坐标变换是矢量控制的“翻译环节”。Clarke变换把三相电流从abc坐标系映射到静止αβ坐标系，公式可以直接写：
Iα = (2/3)·(Ia - 0.5·Ib - 0.5·Ic)
Iβ = (2/3)·(√3/2·Ib - √3/2·Ic)
Park变换再把αβ电流转到同步旋转坐标系：
Id = Iα·cosθe + Iβ·sinθe
Iq = -Iα·sinθe + Iβ·cosθe
这里θe是转子电角度，等于机械角度乘以极对数p。很多人第一次搭模型时，直接用PMSM模块输出的机械角度去做Park变换，结果电流解耦完全不对，d/q轴电流互相串扰。这是高频错误。
SVPWM方面，仿真里不用把每个细节都手写，但要清楚核心逻辑：根据参考电压矢量在αβ平面的位置判断所在扇区，计算相邻两个有效矢量各自的作用时间，再插入零矢量形成对称PWM。需要注意的是，当参考电压幅值超过线性调制范围时，T1+T2会大于PWM周期，必须做限幅处理，否则输出波形会畸变。最大相电压幅值约等于Vdc/√3，这个值在做电压限幅时非常有用。
2.4 速度环PI的简化设计
速度环作为外环，带宽要明显低于电流环，一般取电流环的1/5到1/10，避免内外环响应耦合引起振荡。从q轴电流到机械转速的传递函数，可以近似为一个积分环节：
Ω(s)/Iq(s) = (1.5·p·ψf)/(J·s)
如果希望速度环闭环带宽为ωn，且暂不考虑积分项，速度环比例增益初值可以取：
Kp_n = J·ωn / (1.5·p·ψf)
积分增益先给一个初值，例如Ki_n = Kp_n·ωn/5，然后看阶跃响应慢慢调。我举个例子，同样电机的J = 0.001kg·m²，p = 4，ψf = 0.05Wb，速度环带宽取50Hz，那么ωn = 2π·50 ≈ 314rad/s。Kp_n = 0.001 × 314 / (1.5 × 4 × 0.05) ≈ 1.05，Ki_n ≈ 65.8。这个初值跑出来一般不会疯狂超调，再根据转速响应细节微调。
3. Simulink仿真模型搭建实操过程
3.1 从零开始搭建模型的环境准备
我建议先把整个项目按工程目录来组织，而不是随手建一个slx文件到处放。我的习惯是建一个文件夹，里面放模型文件、初始化参数脚本和说明报告目录。等等，评分不满? 下面这些细节可以帮别人省时间。
| 文件 | 作用 |
|---|---|
FOC_PMSM.slx | Simulink仿真模型 |
init_params.m | 所有电机和控制参数初始化脚本 |
README.txt | 模型版本、依赖和运行说明 |
report | 说明报告存放目录 |
每次打开模型，第一步就是运行初始化脚本，把所有变量加载到工作区。这样可以避免两个最常见问题：一是Simulink模型里某个参数引用变量名，但工作区没有这个变量，导致仿真报错；二是改了电机参数后不知道哪些地方需要同步改。全部参数集中在一个脚本里，改一处就够了。
3.2 电机、逆变器与测量模块的配置
在Simulink库浏览器中搜索“Permanent Magnet Synchronous Machine”，拖入当前模型。双击模块后，把定子电阻、d/q轴电感、永磁体磁链、极对数和转动惯量等参数都填进去。最好通过初始化脚本赋值变量，而不是在参数框里直接填暗数字。PMSM模块的机械端口通常有两种输入方式，一种是输入负载转矩，一种是直接给转速；我在闭环模型里选的是转矩输入，把外部负载转矩接进去。如果是从零开始接触这个模块，注意看模块说明里的初始条件，别忽略转子初始角度设置。
逆变器部分我用的两电平三相桥，可以用“Universal Bridge”这类通用桥模块，桥臂数设3，开关器件根据场景选择理想开关或者带寄生二极管的IGBT模型。说明报告里建议保留死区时间设置，哪怕设一个很短的死区，比如2us，也能让SVPWM输出更贴近实际。直流侧需要一个大电容来稳定母线电压，母线电压由直流电源提供。电流和电压测量可以直接从PMSM模块的输出端口引出来，不一定非要额外加传感器模块。
3.3 电流环、速度环与SVPWM的建模
控制部分我全部用Simulink标准库模块手搭。两个电流环各自用一个PI控制器模块，我在参数设置里把积分输出限幅和整体输出限幅都填好了。速度环同样是PI控制器，输出作为q轴电流给定。为了保险，速度环输出要限幅，不然启动瞬间转速误差很大，q轴电流指令会冲到一个不合理的值，后面齿槽电流就会很难看。
SVPWM模块有两条路线：一是用MATLAB Function块写算法代码，二是纯Simulink模块搭。我更推荐MATLAB Function，因为代码形式便于逻辑检查和后续移植。函数输入是Uα、Uβ和母线电压，输出是三相占空比。函数内部按扇区判断、作用时间计算、限幅、矢量合成依次实现。最后把占空比和固定频率的三角载波比较，生成PWM信号送到逆变器。实际调试中，SVPWM模块最容易错的是扇区判断的符号处理，建议先在开环状态下用一组固定Uα、Uβ验证输出，六拍波形的开关顺序正确后再闭环。
3.4 仿真参数设置与初始化脚本
模型搭完后，仿真参数这一环也很关键。我把求解器设置为定步长模式，步长取1e-5，控制器的采样周期Ts保持1e-4，也就是10kHz。这样的时间尺度，能保证PWM载波和电流采样的时间序列清晰，也不至于让仿真跑得太慢。仿真时长先跑到0.2s，足够观察从零转速到额定转速的完整启动过程。
初始化脚本里除了电机参数，还会写入控制周期、母线电压、PI参数。模型里凡是需要填数字的地方，全部引用脚本里的变量。随后在Simulink的模型属性里，把初始化回调函数设置为执行init_params，这样别人拿到模型后直接点运行，也能自动加载参数，不用到处问“这个变量在哪里定义”。这一点对交付尤其重要。
3.5 结果观察与参数微调
第一次跑起来后，我习惯同时打开几个scope：转速值、三相电流、d/q轴电流、电磁转矩和直流母线电流。先看启动阶段：转速应该平滑上升到给定值，可能有少量超调；q轴电流在启动时冲到限幅值附近，随着转速接近给定逐渐回落；d轴电流基本维持在0附近的小波动；三相电流呈正弦，并且频率随转速上升而升高。
如果发现电流振荡剧烈，第一件事是把电流环带宽从500Hz往下调，而不是继续增加Kp。如果转速超调偏大，就降低速度环比例增益，或者加强速度环积分限幅。每调完一组参数，记录下波形变化，说明报告里就可以用几组对比图展示整定过程，这会比单纯给最终结果更有说服力。
4. 常见问题与排查技巧实录
4.1 电流发散和启动失败
最典型的现象是仿真一开始，电流瞬间冲到几千安培，转速乱跳。这个锅大半要甩给PI输出反号和角度初始值。最常见原因有两个：一是Park变换里的电角度忘了乘极对数，导致解耦完全错误；二是电流反馈的符号方向反了，PI为了把电流拉回给定，反而不断增大误差。排查时我会先断开转速环，给一个固定的q轴电流指令，例如2A，d轴给0，然后看电机的电流波形和转子运动方向是否正确。如果开环固定电流时电机状态平稳，再合上速度环闭环，就能确定问题出在速度环还是电流环。
4.2 角度错误导致的解耦失效
这个坑我在很多模型里都见过。症状是电机能转，但d轴电流不归零，q轴电流有奇怪的低频谐波，动态响应很差。检查后通常会发现，PMSM模块输出的转子角度是机械角度，而Park变换需要的是电角度。正确做法是把机械角度乘以极对数p。另一个隐含问题是角度单位，模块输出有的是弧度，有的是度，如果混用，角度跳变就会让坐标变换输出毛刺。仿真模型里如果能统一使用弧度，会少很多麻烦。
4.3 离散化与代数环问题
当控制部分采用离散PI，而电流反馈直接从连续时间电气网络里引出来时，Simulink很容易报告代数环，或者仿真速度奇慢。原因在于控制器输出计算依赖当前采样电流，而采样电流又依赖控制器刚输出的电压，形成了循环依赖。实际控制器里的采样保持过程在仿真里要用显式模块表示出来。简单处理办法是在电流反馈进控制器前加一个“Unit Delay”或“Memory”模块，模拟ZOH采样，同时把PI控制器的采样时间设置为Ts。这样既能消除代数环，也更贴近真实单片机控制行为。
4.4 SVPWM的临界与过调制
SVPWM进入过调制区时，表现在PWM占空比超过0到1的范围，逆变器输出会产生畸变，电流谐波明显增加。仿真里常常看不出严重问题，但如果后续要做代码生成和硬件测试，这个隐患会被放大。所以我在SVPWM函数内部对占空比做了硬限幅，不让它超出合理范围；同时在电流PI输出端加了电压限幅，参考电压幅值不超出Vdc/√3。这个组合能让模型在大多数工况下都运行在线性调制区，波形更干净，也方便和实际系统对标。
4.5 报告编写与模型交付的堵点
模型调通不算完成，报告写不出来等于没交付。我看过很多模型搭得不错，但报告里没有参数表，没有初始化脚本说明，别人拿到后根本跑不起来。堵点通常集中在三点。第一，截图没有对应运行条件，别人不知道是空载还是带载。第二，参数表里只有最终值，没有整定过程和选择依据。第三，版本依赖不写清楚，不知道在哪个MATLAB版本下验证过。解决方法是把模型文件、初始化脚本、说明报告放在同一个项目目录，README里写清楚打开步骤和运行顺序。
5. 实操心得与几条实用建议
5.1 把参数全部提到初始化脚本里
这是我最想强调的习惯。早期我在Simulink模块里直接填数字，后来换一版电机参数，就得修改十几个界面，还容易漏掉某个限幅值。把所有参数集中到一个init_params.m脚本里之后，模型的每个参数框只填变量名。要改电机或者调PI，只改脚本，仿真前运行一遍，所有模块自动更新。对别人接手也很友好，看脚本就能知道整套系统用了哪些关键参数。
5.2 从低频低压开始，逐级验证
搭建仿真时千万不要一口气把速度环、电流环、SVPWM全合上。我自己的顺序是：先用理想电压源直接开环带电机，确认电机模型参数正确；再接上SVPWM开环调试，给固定电压矢量，观察三相电流和转速；然后接电流环，用固定iq指令验证电流跟踪；最后才接速度环。每走一级都用示波器确认关键信号，尤其是SVPWM的相电压波形，要能看到六个扇区的切换痕迹。这样出了问题很快能定位，比全套接好后猜原因省太多时间。
5.3 说明报告到底写到什么程度
一份好的说明报告不是论文，而是能让一个后来人照着操作复现所有内容。我一般会写这些部分：项目背景和适用范围、模型整体框图、核心控制策略、参数整定过程、仿真结果与现象分析、常见问题与结论。重点不是堆公式，而是把“为什么这样设计”讲清楚。参数表一定要和初始化脚本保持一致，截图下面要标注工况和波形说明。报告里如果有仿真结果的异常现象，哪怕没有完全解决，也可以写出来，对后续接手的人很有参考价值。
最后说一个我真实踩过的坑：有一次我把电流环PI的Ki算错了，少了一个数量级，仿真结果看着还行，转速也稳，但负载一加上去，电流跟踪明显变慢，转速掉下去很久才回来。单纯看空载阶跃根本看不出问题。所以我现在每调完一组参数，一定会加一次满载突加试验，记下最大动态速降和恢复时间。这个习惯救了很多次急，也让我写报告时更有底气。
