从零搭建PMSM FOC仿真模型：原理、步骤与避坑指南

做电机控制的，谁还没被FOC劝退过几次。上个月有个学弟跑来问我，说他在Simulink里搭的PMSM FOC模型调了一周，电流波形依旧不对，id、iq像两根缠在一起的弹簧，怎么解都解不开。我问他Park变换的角度用的是电角度还是机械角度，他愣了几秒，反问我：这不都是角度吗？问题就出在这——FOC表面上看是坐标变换加PI参数的事，但实际动手时，角度基准、模块标定、PI限幅、采样步长，任何一处小偏差都会让结果面目全非。

这篇我就把手把手从零搭PMSM FOC仿真模型的完整流程写清楚，包括电机本体的建模、Clark/Park变换、SVPWM发波、电流环和速度环的参数整定，以及我最想说的——那些文档里从来不会写、但注定让你在调试里耗掉几个通宵的坑。搭建环境是MATLAB R2021b，我调好的完整模型也打包好了，获取方式在第5章开头，需要的朋友可以直接拿去做对照。

1. 动手之前：FOC为什么值得花时间在Simulink里先过一遍

1.1 从"电流解耦"到"双闭环"，FOC的本质其实就一句话

FOC（磁场定向控制）做的事情，说穿了就是让三相交流电机的控制变得像直流电机一样简单。直流电机你只需要控制一个电枢电流，磁场和转矩天然正交，互不干扰。交流电机不行，三相电流在空间上旋转，磁场和转矩强耦合在一起，你踩了一脚油门，磁场也跟着被扰动，根本没法单独控制。

FOC的办法很直接：把三相静止坐标系下的电流先通过Clark变换投影到两相静止坐标系，再通过Park变换旋转到跟着转子磁场走的dq坐标系。到了dq坐标系里，id是励磁电流，iq是转矩电流，两者理论上完全解耦，可以各配一个PI控制器单独伺候。对表贴式PMSM来说，采用id=0控制策略时，电磁转矩近似只和iq成正比，这就意味着你只要控制一个电流量，就能线性地控制转矩。这也就是为什么FOC在永磁同步电机驱动里这么普及，本质上是把一个多变量、强耦合的非线性系统，一个变换就解耦成了两个独立的单输入单输出回路。

如果把FOC这套东西拆成模块，其实就是四件事：坐标变换、电流内环、速度外环、SVPWM发波。后面所有的工作，都是围绕这四个模块做细节。

1.2 三种搭建路线对比：为什么我推荐自己从零搭

搭一套FOC模型，面前通常有三条路，我先说结论：学习用，自己从零搭；项目赶进度，在官方例程上改；只想出个波形不想懂原理，直接用Motor Control Blockset拖出来。咱们这篇的目标是既懂原理又做出能跑的模型，所以按从零搭的路子走。

我的观点很明确：哪怕你最后还是要用官方库，第一遍也值得自己搭一次。Simulink的价值不在于“拖模块”，而在于它把数学公式变成了看得见信号的连线——搭一遍坐标变换，你对角度、符号、坐标定向的理解，比读十遍教材都有用。

2. 搭建前的准备工作：版本、工具箱和电机参数从哪来

2.1 软件版本与工具箱：够用就行，别追求最新

先说版本。我搭这套模型用的是MATLAB R2021b，这套流程在R2019b到R2023a之间都能复现，差异不大。很多新手喜欢装最新版，其实没必要，工具箱路径和模型兼容性反而容易踩坑。学生党建议直接用R2021b或R2022b，网上资料最多，遇到报错好搜索。

需要装齐的工具箱只有这几个：

• Simulink：不用多说了。
• Simscape Electrical：电机模型、逆变器桥臂、直流电源都在这个库里。老版本里它叫SimPowerSystems，R2019b之后统一改成Simscape Electrical，界面有变化，但模块路径基本一致。
• DSP System Toolbox：非必需，但如果要在电流反馈上加低通滤波，或者拿FFT分析谐波，建议装上。
• Control System Toolbox：非必需，但可以用sisotool辅助做PI整定，我个人觉得挺好用。

不建议为了省事直接加载Motor Control Blockset全套库。不是它不好，而是它封装程度太高，内部的PWM生成、ADC采样、坐标变换全是“黑盒”，波形一旦不对，你连从哪里入手查都不知道。自己搭一遍，基本概念通了，再回头用那些库，反而觉得轻松很多。

2.2 电机参数清单：仿真模型的"出厂数据"怎么凑出来

没有真实电机，仿真参数怎么定？很多教程直接丢一组参数，但换个电机就不知道怎么填。我这里给一个最常用的参数清单，顺便说下我对每个参数的工程理解。

没有手册的情况下，ψf可以通过一个很简单的实验估算：用外力把电机拖到某个转速，用示波器测任意两相之间的线反电动势峰值E_line_pk，同时记录机械转速n（rpm），那么电角速度ωe = 2π·p·n / 60，ψf ≈ E_line_pk / (√3·ωe)。前提是反电动势波形接近正弦，实测够用。

如果你手头没有真实电机，直接用下面这组教学参数，它来自MATLAB官方PMSM例子的变体，参数之间是自洽的：

• Rs = 0.958 Ω
• Ld = Lq = 5.25 mH（表贴式）
• ψf = 0.1827 Wb
• 极对数 p = 4
• J = 0.0006329 kg·m²
• 额定转速 3000 rpm
• 额定直流母线电压 300V

建议第一版仿真基于这组参数跑通，因为官方例子都拿类似参数做验证，网上遇到问题最好对照。

3. 核心搭建流程：从坐标变换到闭环控制的分步实操

3.1 第一步：电机本体——用Simscape PMSM模块还是自己搭数学模型

电机本体有两种做法。

做法A：直接用Simscape Electrical自带的Permanent Magnet Synchronous Motor模块。双击模块，设定极对数、定子电阻、电感和永磁磁链，它会自动处理电压方程和机械方程，还会提供一个m端口，直接引出三相电流、转子机械角度、机械转速和电磁转矩。我推荐第一次用这个，省事、不易错，而且它和逆变器桥臂的物理域接口是成套的。

做法B：完全不用电气库，用Simulink基本模块把PMSM电压方程和运动方程亲手搭一遍。好处是把公式亲手写一遍，理解会更深；坏处是后面接逆变器、接示波器都麻烦，单位换算也容易出bug。建议先用做法A把系统跑通，之后再回头用做法B验证自己的理解。

用PMSM模块时有几个细节需要注意：

• Back EMF waveform 选 Sinusoidal。
• Mechanical input 通常选 Torque TL，端口接常数0代表空载，或者接一个阶跃负载。
• 初始转子机械角度：注意和你坐标变换的角度基准对齐。默认设成0，Park变换用的角度初始值也必须对应0。
• 三相电枢绕组端子A/B/C接逆变器输出，Simscape的物理信号需要经过PS-Simulink Converter转成普通Simulink信号；反过来，控制信号进逆变器门极前，也要用Simulink-PS Converter回到物理域。千万不要把信号域和物理域搞混。

3.2 第二步：坐标变换模块的搭建与角度对齐

电机反馈出来的三相电流ia、ib、ic，先做Clark变换，把三相静止坐标变成两相静止坐标：

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必须是电角度。电机反馈回来的是机械角θm，必须先乘极对数：θe = p·θm。我见过太多新手把机械角直接接进Park变换，后果就是id、iq完全解耦不了，波形是两个带相位差的振荡信号，PI怎么调都压不出直线id=0。

判断角度方向对错有个土办法：让电机不带负载，把iq给定成一个小正值，如果id能保持在0附近、iq等于给定值，角度基本没问题；如果id和iq互相串扰，优先怀疑角度基准。在Simulink里，坐标变换可以用Fcn模块直接写表达式，但我更推荐用MATLAB Function，因为后续做C代码生成时这条路是通的。实现代码：

function [id, iq] = park_transform(ia, ib, ic, theta_e) % Clark变换 i_alpha = 2/3 * (ia - 0.5*ib - 0.5*ic); i_beta = 2/3 * (sqrt(3)/2 * ib - sqrt(3)/2 * ic); % Park变换 id = i_alpha * cos(theta_e) + i_beta * sin(theta_e); iq = -i_alpha * sin(theta_e) + i_beta * cos(theta_e); end

需要提醒一下，我这里用的是等幅变换，不带2/3归一的系数。等幅和等功率变换在系数上有差别，会影响PI参数的具体数值，但不会影响控制逻辑本身。你只要整个链路前后统一就行，别混用。

3.3 第三步：电流环PI调节器的参数初值计算

电流环PI参数是FOC仿真第一个硬骨头。网上公式很多，我给一个工程上够用的近似整定法。

假设PI形式是连续域的 Kp + Ki/s，电流环被控对象近似为一阶惯性加电阻电感环节，给定电流环带宽ωc（rad/s），经验公式：

Kp ≈ L * ωc Ki ≈ Rs * ωc

原理是让PI零点抵消被控对象的极点。电流环带宽一般取开关角频率的1/20到1/10。开关频率10kHz，开关角频率约62832 rad/s，电流环带宽取1000到3000 rad/s都合理。我们取ωc = 1000 rad/s，代入上面那组参数：

Kp = 0.00525 * 1000 = 5.25 V/A Ki = 0.958 * 1000 = 958 V/(A·s)

Simulink里用PID Controller模块，形式选PI，采样时间可以先设成连续域，跑通后再离散化。但有一件事必须现在就做——给PI输出加限幅。电流环输出是ud和uq，限幅要小于逆变器线性输出范围。母线电压300V时，SVPWM最大相电压幅值理论上能到Vdc/√3约173V，所以ud、uq限幅先设到±150V比较合适。不设限幅，启动瞬间PI可能把ud、uq积分到几千伏，波形直接崩。

实际调试顺序我强烈建议按这个来：

1. 先开环：固定给ud=0、uq=5V，看电机能不能转起来，确认电机和逆变器连接没问题。
2. 再电流环闭环：速度环断开，给定iq_ref=1A，id_ref=0，用上面的初值试，观察iq跟踪情况。如果振荡或发散，把Kp和Ki往下降。
3. 最后才上速度环：速度环带宽大概只有电流环的1/10，别一上来就把速度环调得很激进。

3.4 第四步：速度环与SVPWM逆变器的连接

速度环结构很简单：速度给定（rpm）减去反馈转速，进速度PI，输出就是iq_ref，id_ref直接给定0。速度PI参数起步建议：

• Kp先设小，比如0.1到0.5（具体跟转动惯量有关）；
• Ki比Kp低一个数量级；
• 先观察转速波形，再逐步加大参数；
• 波形持续振荡，通常是Kp太大或Ki积分太强；稳态误差明显，则需要加强Ki。

SVPWM是最后一块拼图，作用是把uα、uβ合成出三相逆变器桥臂的占空比。Simulink里实现有两条路：一是用Simscape Electrical自带的PWM发生器配合Universal Bridge，二是用MATLAB Function自己写SVPWM算法。我的建议是，哪怕后面用自带PWM发生器，也至少手写一遍SVPWM。矢量合成这个思想搞懂了，对理解调制极限、弱磁都会很有帮助。

以下是我验证过的一个SVPWM参考实现，输入uα、uβ、母线电压Vdc、PWM周期Tpwm，输出三相比较值Ta、Tb、Tc：

function [Ta, Tb, Tc] = svpwm_ref(u_alpha, u_beta, Vdc, Tpwm) % 扇区判断 V1 = u_beta; V2 = (sqrt(3)/2)*u_alpha - 0.5*u_beta; V3 = -(sqrt(3)/2)*u_alpha - 0.5*u_beta; N = (V1 > 0) + 2*(V2 > 0) + 4*(V3 > 0); % 矢量作用时间 X = sqrt(3)*u_beta*Tpwm/Vdc; Y = (1.5*u_alpha + sqrt(3)/2*u_beta)*Tpwm/Vdc; Z = (-1.5*u_alpha + sqrt(3)/2*u_beta)*Tpwm/Vdc; switch N case 1, T1 = Z; T2 = Y; case 2, T1 = Y; T2 = -X; case 3, T1 = -Z; T2 = X; case 4, T1 = -X; T2 = Z; case 5, T1 = X; T2 = -Y; otherwise, T1 = -Y; T2 = -Z; end % 限幅和比较值 T1 = min(max(T1, 0), Tpwm); T2 = min(max(T2, 0), Tpwm); t1 = (Tpwm - T1 - T2)/2; t2 = t1 + T1/2; t3 = t2 + T2/2; switch N case 1, Ta = t1; Tb = t2; Tc = t3; case 2, Ta = t2; Tb = t1; Tc = t3; case 3, Ta = t3; Tb = t1; Tc = t2; case 4, Ta = t3; Tb = t2; Tc = t1; case 5, Ta = t2; Tb = t3; Tc = t1; otherwise, Ta = t1; Tb = t3; Tc = t2; end end

生成的三路比较值送入Simulink-PS Converter，接Universal Bridge的门极，直流电压源接到桥臂正负端，桥臂输出接PMSM三相端子。需要留意的是，逆变器门极信号一旦是PWM，Simulink求解器就会变得很硬，第4章里说的采样步长问题至关重要。

4. 仿真调试中我踩过的六个坑，以及对应的波形特征

FOC仿真能不能跑起来，很大程度上是拿时间换来的。下面这六个坑是我和身边同事、学弟们踩过最多的，我按“症状—原因—处理”来写，方便你对号入座。

4.1 电机原地抖动不转：转子初始角度没对齐

症状很明显：一启动电流瞬间就到保护值，转子在零位附近来回抖，机械转速一直振荡，电磁转矩根本出不来。

最常见原因就是Park变换用的初始电角度和电机模块内部的转子初始位置不一致。比如电机模块里默认转子初始角是π/6，你这边从0开始做坐标变换，算法认为转子在0度，实际转子在30度，id、iq刚闭合就全是偏差，系统当然稳不住。

处理办法：在PMSM模块里把初始转子机械角度设成0，并且和你坐标变换的零位保持一致。如果后续接编码器接口，也要保证编码器零点对应A相轴线。真机调试时，这个问题的同义表达就是“转子初始位置检测”，这是所有无感FOC项目启动时最硬的骨头。仿真里先把起始位置对齐，确认算法没问题了，再回头研究预定位、高频注入这些内容。

4.2 电流环一闭环就炸：PI饱和与积分限幅

症状：电流环闭环瞬间，id或iq直接冲到几十上百安培，仿真数值发散，Scope里波形完全糊掉。

原因大概率不是公式错了，而是PI输出没有限幅，或者积分项没做处理。启动瞬间，给定1A，反馈0，误差1A，Kp=5时输出5V，这还好——但Ki积分在0.1s里能积出接近100V，加上Kp项，很容易就推高到几百伏。逆变器根本给不出这个电压，系统自然崩。

处理办法：

• PID Controller模块必须勾选输出饱和，ud、uq限幅先给±150V（对应母线300V）。
• 积分器最好勾选anti-windup，或者手动加积分限幅。
• 调试时不要给大幅值阶跃，先给定iq_ref=0.1A这类小信号，观察电流环响应再放大。

4.3 d轴电流收敛不到零：角度偏差或坐标变换符号问题

症状：id不是0，而是始终带一个直流偏置，比如稳定在0.5A附近；iq能勉强跟住给定，但波形上有同频振荡。

先查坐标变换符号。Park变换的iq公式，有的资料写成-iα·sinθ + iβ·cosθ，有的写成iα·sinθ - iβ·cosθ，两种对应不同旋转方向，必须和你θ的正方向匹配。再查θe方向，PMSM模块默认输出逆时针正转的机械角，如果电机在仿真里是反转的（转速显示负），那θe要取反。

处理办法：先用开环验证角度链路。开环给uq=2V，让电机转起来，用Scope同时记录θe、iα、iβ，看αβ坐标下的电流矢量方向是否和θe定义一致。这个方法在硬件调试里也通用，比单纯盯id、iq更快定位问题。

4.4 高速下转速上不去：反电动势饱和

症状：给定3000rpm，实际到2000rpm就上不去了，iq已经顶到限幅，uq长期处于饱和状态。

这是物理规律决定的。uq方程里有一项ωe·ψf，转速越高反电动势越高，逆变器能输出的电压上限固定，留给电流调节器的裕量会越来越小。到某个转速，dq轴电压矢量幅值达到SVPWM的六边形边界，电流环即使饱和也无法继续升压。

处理办法：仿真里先把直流母线电压从300V提高到350V，看最高转速有没有上升，一旦上升就可以确认是电压饱和问题。真实产品里，这时候就要上弱磁控制——让id变成负值，用去磁效应削弱气隙磁场，给电压挤出裕量。弱磁话题很大，这里先不展开，但你要知道：id=0不是永远最优，越过转折速度后必须弱磁。

4.5 开关频率与采样步长导致的电流毛刺

症状：电流波形在正弦基波上叠了一层高频锯齿，id、iq有周期性毛刺，PI输出跟着抖。

原因通常是PWM开关频率10kHz，但Simulink固定步长设成1e-4甚至更大，整个系统在PWM周期内只算了几个点，逆变器开关纹波根本分辨不出来。另一个模型层面的坑：Universal Bridge是理想开关模型，开关瞬间数值变化很陡，步长太大还会引起过零检测异常。

处理办法：

• 固定步长必须小于PWM周期的1/20。10kHz对应100μs周期，步长建议5e-6s甚至更小。
• 如果只关心控制性能、不关心开关纹波，把逆变器换成平均模型，步长可以放到1e-5，仿真速度快很多。
• 电流采样点尽可能对齐PWM中心。仿真里可以用PWM发生器产生的中心对齐信号触发采样，这和真实MCU里的ADC触发逻辑是一致的。

4.6 电流波形是"马鞍形"？别慌，那可能说明你搭对了

症状：相电流波形不是教科书里那种干净正弦，而是有点马鞍形，中间有隆起。

原因：SVPWM调制出的相电压本身就叠加了零序分量，所以相电压不是纯正弦，但线电压和电流经过电机电感滤波后依然接近正弦。如果相电流是干净正弦，反而可能是用了SPWM正弦调制；马鞍形恰恰是SVPWM的特征之一。

处理办法：不要急着改。用Scope同时看线电流和FFT，确认基波频率、幅值是否符合预期，谐波是否以开关频率附近的边带为主。如果出现明显的低次谐波，再回头查坐标变换和角度链路。

5. 拿到仿真文件之后：怎么改成你自己的电机参数，以及后续进阶

5.1 模型文件使用说明与参数修改清单

模型文件我打包好了，下载链接放在评论区置顶，拿不到的直接私信我。文件目录大概是：

• PMSM_FOC_Demo.slx：主模型，包含电机、逆变器、坐标变换、双闭环、SVPWM；
• init_pmsm_params.m：参数初始化脚本，所有电机参数和PI初值都在里面；
• README.txt：说明文件，记录版本、模块路径、常见问题。

拿到之后第一次打开，建议按这个顺序操作：

1. 先运行init_pmsm_params.m，把工作区参数刷出来。
2. 打开主模型，Ctrl+D更新图，确认没有模块缺失或版本兼容告警。
3. 在R2021b上直接Run，先复现我第3章说的波形，确认模型在自己手里是通的，再开始改参数。

要改你自己的电机时，只需在init_pmsm_params.m里改电机常量和母线电压，然后重新算Kp、Ki。我建议把计算逻辑直接写在脚本里，改完Rs、L自动算出初值，不要手算，既省时间又减少低级错误。模型里凡涉及PI限幅、逆变器母线电压的地方都要跟着参数一起检查——很多“改了参数就崩”的案例，其实就是限幅没跟着改。

5.2 从仿真到嵌入式代码的进阶路线

仿真模型跑通只是第一步。要做真机、做毕设或做产品，通常有两条路。

第一条路是自动代码生成。把控制器部分单独做成一个模型，用Simulink Coder或Embedded Coder生成C代码，再移植到STM32、TMS320F28379D这类MCU上。代码生成前必须把求解器配置为Fixed-step discrete，采样时间对应PWM中断周期，同时确保所用模块都支持代码生成。坐标变换、PI、SVPWM基本都是支持的，最容易翻车的是信号类型不一致，导致生成一堆低效的double转float。

第二条路是手动移植。手写C代码，把Simulink模型当算法参考。我接触的不少工程师更偏爱这种方式，因为MCU端底层驱动（ADC触发、PWM比较值刷新、编码器接口）本来就是自己写的，算法部分照抄模型里的公式就行。这条路要求你能从公式直接写出坐标变换、PI和SVPWM——所以还是要回到那句老话：自己从零搭一遍模型，后面写移植代码会轻松很多。

5.3 进阶玩法：MTPA、弱磁、无感FOC和无感启动

基础版id=0双闭环跑通后，可以加的东西非常多，按由易到难排：

• MTPA（最大转矩电流比）：内嵌式PMSM（Ld≠Lq）用id=0不是最优，需要根据电机参数算最优id、iq分配，让单位电流输出最大转矩。
• 弱磁控制：上一章提到的高速段升压问题，在电压饱和边界上引入id负向补偿，常用方案是电压外环法。
• 无感FOC：去掉编码器，用滑模观测器或龙伯格观测器从反电动势里估算转子位置和转速。这是个公认的大坑，尤其低转速和带载启动阶段。
• 无感启动策略：先预定位保证初始方向，再开环强拖到一定转速，观测器稳定后再切换回闭环。开环切闭环那个瞬间，经常能看到速度跳变和电流冲击，是整个调试过程里最磨人的部分。

每一条都能单独写一篇长文，这里先给大家指个方向。等你们基础模型都跑顺了，我们再一个一个拆开聊。

最后分享一点个人体会。我搭完第一版PMSM FOC模型时，最大的收获反而不是模型本身，而是重新理解了那些不起眼的细节——角度基准、PI限幅、采样时刻，每一个都在决定整套系统能不能稳定。Simulink的好处是，你改一个参数就能立刻看到波形变化，所以我特别建议你把模型参数故意改错几次：把θe不乘极对数，把Park变换符号反过来，或者把PI限幅去掉，亲眼看看波形怎么变成一锅粥，再一步一步改回来。这样折腾一晚上，你对FOC的理解会比直接照抄一份完美模型深得多。模型文件照旧在评论区或私信获取，有问题留言，我尽量回。