MATLAB/Simulink永磁同步电机FOC仿真：数学模型、SVPWM与调试实战

前两年我接了一个电动工具控制器项目第一次真正把永磁同步电机PMSM驱动压在MATLAB里跑起来的时候我花了一个星期的晚上才把SVPWM那套查表逻辑弄清楚。后来带新人我发现大家卡住的点其实高度一致不是不会点鼠标而是对PMSM数学模型、坐标系、占空比到底是电压还是开关信号这些基础概念没打通。如果你正打算用MATLAB/Simulink搭一套永磁同步电机仿真模型或者已经搭了但波形总不对这篇文章正好能帮你把整条链路一次理顺。文章会从“为什么先搭仿真”讲起接着给出PMSM数学模型的完整梳理然后逐步拆解Simulink建模步骤附上可以直接抄走的初始化脚本和SVPWM核心函数代码最后把所有我踩过的坑、实测过的排查方法整理成速查思路。适合刚接触FOC的应届生、在调硬件但想先把算法跑通的工程师以及做MATLAB课程设计需要电机仿真模型的同学。全文不绕弯子直接按“理论→建模→代码→调参→排障”的顺序来。1. PMSM仿真模型整体设计与思路拆解1.1 为什么我先从仿真开始而不是直接上驱动板在你把控制板接到电机和逆变器之前先用仿真把控制算法吃透是我这几年做电机项目的一个基本原则。原因很简单仿真里你能看电流波形、转矩波形、转速阶跃响应、SVPWM调制波还能随意改参数看效果而硬件上随便一个接线错误、一个相序反接轻则烧驱动管重则把板子和电机一起送走。相比之下仿真模型的调试成本几乎为零一脚急刹车都花不了几毛钱。不过仿真也有自己的坑最大的坑就是“仿真能跑不代表硬件能跑”。如果仿真模型里的采样时间、逆变器死区、母线电压约束都被忽略那仿真波形再漂亮上板也会翻车。所以我搭仿真模型的时候会刻意保留几个硬件关键约束直流母线电压限制、PI输出限幅、PWM占空比饱和。这样仿真跑出来的结果才更有参考价值。另一个常见误解是“PMSM仿真必须用Simscape里的电机模块”。其实如果你只是想把控制算法调通用基本的数学方程建模也是完全可行的。但考虑到大多数读者要的是“能往上扩展、能接硬件”的模型我更推荐直接用Simscape Electrical里的Permanent Magnet Synchronous Machine模块它已经帮你把d/q轴方程、反电动势、机电耦合都实现好了只需要填参数。这样可以把精力集中在控制部分。1.2 仿真模型到底要搭成什么样FOC双闭环上图前先讲清楚控制结构。这套仿真模型我采用的是最主流的双闭环FOCField-Oriented Control磁场定向控制外环是速度环输入是转速指令和反馈转速输出是q轴电流指令 iq_ref内环是电流环d轴电流指令 id_ref 通常给0表贴式PMSM一般这样内置式PMSM以后你可以再扩展MTPA电流环PI输出的是 d/q 轴电压指令 u_d、u_qu_d、u_q 经过反Park变换得到 u_alpha、u_betau_alpha、u_beta 进入SVPWM调制得到逆变器开关管的占空比逆变器输出三相电压给电机电机产生转矩和转速从电机侧测回来的三相电流经过Clark变换、Park变换变成可反馈的 i_d、i_q。整个链路不复杂难点集中在两个地方坐标变换系数要一致SVPWM调制要正确。这两个点我后面都会重点展开。1.3 方案选型Simscape电机模块还是自建数学方程我见过不少教程用纯Simulink数学模块搭建PMSM本体好处是“模型完全透明”坏处是搭建费时且容易埋bug。我的建议是分场景如果你是初次接触PMSM仿真直接用Simscape Electrical的Permanent Magnet Synchronous Machine模块把参数填对就行如果你后面要接入控制器硬件在环HIL或者做参数辨识再考虑自建方程也不迟如果只是想验证某个控制算法比如滑模观测器、模型预测控制那自建一个简化的PMSM方程模块反而更灵活。我在下面给的示例模型就统一使用Simscape电机模块。版本上MATLAB 2019b及更新版本都能正常操作唯一要注意的是不同版本里Simscape电气库的模块路径可能略有不同建议直接用模块搜索框输入名字找。2. 数学基础PMSM数学模型与坐标变换2.1 为什么三相坐标系的电机方程很难直接用来做控制PMSM定子三相绕组产生的磁动势在空间上是静止的而转子在旋转所以定子电感、磁链等参数会随着转子位置不断变化。如果你坚持在a/b/c三相静止坐标系下写电压方程那自感、互感都和转子角度强耦合方程又长又难解。更麻烦的是三相电流是交流量控制器如果直接去跟踪交流正弦信号PI调节器很难做到无静差。所以工程上会把“三相静止坐标系”变换到“以转子磁极方向为基准的旋转坐标系”。在这个坐标系里正弦量变成了直流量PI控制器就能轻松做到无静差跟踪。这就是坐标变换的意义。从实现上PMSM建模必须用两条变换链先做Clark变换abc → αβ再做Park变换αβ → dq控制完成后再做反Parkdq → αβ最后通过SVPWM送到三相逆变器。2.2 dq旋转坐标系下的PMSM核心方程建模时我习惯把最常用的PMSM数学模型写出来方便和Simulink模块参数对照也方便排查问题的时候知道每个量代表什么。电压方程u_d Rs·i_d Ld·(d i_d/dt) - ω_e·Lq·i_qu_q Rs·i_q Lq·(d i_q/dt) ω_e·(Ld·i_d ψ_f)磁链方程常用于推导ψ_d Ld·i_d ψ_fψ_q Lq·i_q转矩方程T_e 1.5·pp·[ψ_f·i_q (Ld - Lq)·i_d·i_q]机械运动方程J·(d ω_m/dt) T_e - T_L - B·ω_m其中 ω_e pp·ω_mpp是极对数ω_m是机械角速度T_L是负载转矩B是粘滞摩擦系数ψ_f是永磁体磁链。注意这里的系数 1.5 是在“等幅值Park变换”前提下得到的如果你采用等功率变换转矩方程系数会变成 3/2 以外的值。这一点很多人容易踩坑。看到转矩方程里的 (Ld - Lq)·i_d·i_q 这一项了吗这就是内置式PMSMIPMSM特有的磁阻转矩。当Ld不等于Lq时如果只采用 id_ref0 控制磁阻转矩会被浪费掉想用MTPA的话需要根据电机具体参数计算最优i_d。表贴式PMSM一般Ld≈Lq所以id0控制够用。2.3 Clark/Park变换和那个最容易搞混的系数坐标变换公式网上很多我只强调最容易出错的点系数2/3和1。Clark变换等幅值i_alpha (2/3)·(i_a - 0.5·i_b - 0.5·i_c)i_beta (2/3)·(√3/2·i_b - √3/2·i_c)Park变换i_d i_alpha·cos(θ_e) i_beta·sin(θ_e)i_q -i_alpha·sin(θ_e) i_beta·cos(θ_e)反Park变换u_alpha u_d·cos(θ_e) - u_q·sin(θ_e)u_beta u_d·sin(θ_e) u_q·cos(θ_e)这里的θ_e是转子电角度不是机械角度计算时要乘以极对数θ_e pp·θ_m。很多仿真模型跑出来波形乱成麻花十有八九是忘记乘极对数。MATLAB自带的PMSM模块内部已经按等幅值约定处理了你只要保证外部的Clark/Park变换也采用等幅值形式两边才能对得上。另一点如果你在代码里看到教材上等功率变换的 √(2/3) 系数千万别直接照抄先确认你整个模型用的是哪套标准。我在实际项目中就见过仿真和硬件数据对不上的案例追了两天最后发现是变换系数不一致。3. Simulink模型搭建实操模块布局与参数设置3.1 你需要哪些模块搭建模型的思路是“先把动力部分搭好再搭控制部分最后连反馈”。下面是我常用的模块清单版本不同路径会变但模块名基本稳定。建议直接在Simulink搜索框里输入名字来找PowerguiSimscape Electrical Specialized Power SystemsPermanent Magnet Synchronous MachineSimscape Electrical Specialized Power Systems MachinesUniversal BridgeSimscape Electrical Specialized Power Systems Power ElectronicsDC Voltage Source直流母线电源Current Measurement电流测量接两相即可Machine Measurement Demux把PMSM的测量输出拆分出来PWM GeneratorSimscape Electrical Specialized Power Systems Control BlocksMATLAB Function用来写SVPWMDiscrete PI Controller离散PI控制器电流环2、速度环1各种常数、加法器、示波器这里的PWM Generator很关键。它的输入是调制信号输出是6路驱动脉冲正好可以直接推Universal Bridge。SVPWM算法在MATLAB Function里算出每相占空比再送给PWM Generator去和内部三角载波比较。3.2 动力部分怎么连线动力部分从直流侧开始DC Voltage Source正极接Universal Bridge的直流正端负极接直流负端。Universal Bridge的交流侧三根线分别接到PMSM定子三相绕组。中间可以串两个Current Measurement来采样电流分别放在A相和B相上就行C相电流可以由基尔霍夫定律算出来。PMSM模块的机械端口接一个常值负载转矩或阶跃负载用来模拟突加负载的情况。机械侧还有一个输出接口把转子位置、转速等信息通过Machine Measurement Demux引出到Simulink控制侧。这里有一个容易出问题的地方相序接反了转矩会是负的转速方向会和你预期相反。做仿真时很容易以为是自己控制器写错了其实只是把两相线换一下就行。我一般会把A/B/C三相标清楚如果出现转矩反向先交换任意两相试试。3.3 控制通路怎么连控制侧我把参考转速设为 Speed_ref例如150 rad/s约1432 rpm。速度环PI输出 iq_ref电流环d轴参考给0。三相电流反馈经过Clark/Park变换得到i_d、i_q与参考比较后进入两个电流环PI。PI输出u_d、u_q再经过反Park变换得到u_alpha、u_beta进入MATLAB Function的SVPWM模块。SVPWM算出来的占空比信号 d_a、d_b、d_c 送给PWM Generator的调制波输入端。PWM Generator内部会用开关频率的三角载波比较输出6路脉冲到Universal Bridge。这里要格外注意PWM Generator的调制波输入比例关系。我用的是0到1的占空比形式你可以参考帮助文档确认输入范围。最后别忘了一个细节电流反馈的符号。Simscape电流测量的方向如果接反反馈就变成了负反馈变正反馈电流一路发散。我建议在跑闭环前先用一个开环电压给电机确认电流测量和坐标变换的极性正确再做闭环控制。3.4 求解器和仿真步长设置我在仿真里用的求解器是Powergui Discrete采样周期设为开关周期的1/20左右。比如开关频率10kHz开关周期Ts_ctrl100us仿真步长Ts_sim可以取5us。设置太大会出现数值振荡看起来像电流纹波其实是数值误差设置太小则仿真速度慢调试时等着难受。如果只想快速验证控制算法可以把逆变器和PMSM一起换成平均模型SVPWM输出的占空比直接乘以直流母线电压当三相电压源用不接Universal Bridge和PWM Generator。这样仿真步长可以放宽到10us甚至20us速度会快很多。但代价是看不到开关纹波和逆变器非线性效应两种模型各有用途不要混为一谈。4. 完整代码初始化脚本与SVPWM函数4.1 初始化脚本变量集中管理在Simulink模型里直接敲数字很容易乱我习惯把所有参数集中在一个脚本 init_pmsm_sim.m 里每次运行前先跑一遍初始化脚本。这样改参数不用满模型找也方便以后做参数扫描。%% PMSM仿真模型初始化脚本 % 使用前先在此工作区运行再打开Simulink模型 %% 直流母线与开关频率 Vdc 300; % 直流母线电压单位V Fs 10e3; % 逆变器开关频率单位Hz Ts_ctrl 1/Fs; % 控制周期单位s Ts_sim Ts_ctrl/20; % 离散仿真步长 Tsim 0.5; % 仿真时长单位s %% 电机额定参数以某款3kW级PMSM为参考 pp 4; % 极对数 Rs 0.185; % 定子电阻单位Ohm Ld 0.835e-3; % d轴电感单位H Lq 1.95e-3; % q轴电感单位H flux 0.0725; % 永磁磁链单位Wb J 0.003; % 转动惯量单位kg*m^2 B 5e-4; % 粘滞摩擦系数单位N*m*s/rad %% 速度环参考 Speed_ref 150; % 目标转速单位rad/s约1432rpm %% 电流环PI参数计算 % 电流环带宽选1000 rad/sGPT惯用方法Kp L*bwKi Rs*bw bw_cc 1000; Kp_id bw_cc * Ld; Ki_id bw_cc * Rs; Kp_iq bw_cc * Lq; Ki_iq bw_cc * Rs; %% 速度环PI参数初值 % 速度环带宽选120 rad/sKp先用惯量法估Ki再根据阶跃响应微调 bw_sc 120; Kp_sp bw_sc * J / (1.5 * pp * flux); Ki_sp bw_sc * Kp_sp / 10; %% PI输出限幅 % SVPWM在线性区能输出的最大相电压幅值约为Vdc/sqrt(3) u_max Vdc / sqrt(3) * 0.95; iq_max 20; % 电流限幅单位A按电机额定或驱动器能力填写参数单位我写得很清楚复制到MATLAB里后可以直接用。需要特别提醒电机参数一定以实际电机型号手册为准不要直接拿我给的示例参数去对接真机否则可能偏差很大。我给的这套参数只是让模型能转起来方便你验证控制链路。4.2 SVPWM函数用零序注入方式实现说到SVPWM很多教程会先摆一大段扇区判断、作用时间计算、切换点顺序这也是我当年花了最多时间的地方。但后来我发现在Simulink仿真里有一个实现上更简洁、理解上也不绕的等价方法零序注入法又称载波比较型SVPWM。它将SVPWM等效为在正弦调制波上注入一个零序分量然后与三角载波比较得到的结果在谐波性能上与经典七段式SVPWM基本一致。而且代码极短不容易写错。下面是完整的MATLAB Function函数代码。你在Simulink里新建一个MATLAB Function模块把函数名和输入输出设置成和下面一致即可function [d_a, d_b, d_c] pmsm_svpwm(u_alpha, u_beta, Vdc) % 基于零序注入的SVPWM实现 % 输入u_alpha、u_beta 为相电压指令VVdc为母线电压V % 输出d_a、d_b、d_c 为占空比范围0~1 % 1. 由alpha/beta轴电压反推三相电压指令 ua u_alpha; ub -0.5 * u_alpha 0.8660254037844386 * u_beta; uc -0.5 * u_alpha - 0.8660254037844386 * u_beta; % 2. 计算零序分量并注入 umax max([ua, ub, uc]); umin min([ua, ub, uc]); offset -0.5 * (umax umin); va ua offset; vb ub offset; vc uc offset; % 3. 转为占空比 d_a (0.5 * Vdc va) / Vdc; d_b (0.5 * Vdc vb) / Vdc; d_c (0.5 * Vdc vc) / Vdc; % 4. 限幅防止过调制 d_a max(0, min(1, d_a)); d_b max(0, min(1, d_b)); d_c max(0, min(1, d_c)); end这段代码的核心逻辑就三条由αβ反推三相、算零序、转占空比。之所以它能等效SVPWM是因为SVPWM本质上就是在正弦调制波上叠加合适的零序分量使调制波变成鞍形波从而把直流母线电压利用率从SPWM的0.5倍提升到SVPWM的0.577倍左右。我在测试时会把d_a、d_b、d_c接到示波器上观察正常情况下它们应该是对称的鞍形波幅值在0到1之间频率与电流频率一致。如果看到平顶波说明正在过调制电流会变差。4.3 如何快速验证代码对不对很多人把模型搭完一运行看到转速飞上去就觉得代码没问题这是不够的。我建议按下面几步做快速验证先让速度环参考为0给一个很小的iq_ref比如1A看电流环是否能把i_q稳定在1A附近固定iq_ref给一步负载转矩看转速是否下降后又能稳住观察d_a/d_b/d_c波形正常应该是正弦加零序的鞍形波而不是纯正弦观察三相电流稳态时应该是对称正弦THD肉眼看起来没有明显毛刺把速度环闭合给一个速度阶跃观察超调量和调节时间。如果这几步都正常说明模型基本可靠可以放心往里面加更高级的控制算法。5. PI参数整定思路与调试流程5.1 电流环PI到底怎么算电流环是整个控制系统的内环响应要快。实际上电流环被控对象近似为一阶惯性环节G(s) 1/(Ls Rs)。用PI控制器零点可以抵消对象极点从而把电流环变成近似一阶低通特性。工程上常用下面这种初值估算方法Kp L · bw_ccKi Rs · bw_cc这里的 bw_cc 是期望电流环带宽单位是rad/s。我给的脚本里带宽取1000 rad/s对应时间常数约1ms对于10kHz开关频率来说留了足够的相位裕度。如果你发现q轴电流响应慢可以把带宽调高如果电流出现高频振荡就把带宽降下来。比如脚本里的Lq1.95mHRs0.185Ω那么Kp_iq 0.00195 × 1000 1.95Ki_iq 0.185 × 1000 185这个数值范围在仿真里很合适。注意实际电机参数里Lq通常比Ld大所以q轴Kp会比d轴Kp大一些这是正常现象。5.2 速度环PI初值怎么给速度环是外环带宽必须明显低于电流环。所以速度环带宽我取电流环的1/10左右也就是100~150 rad/s。速度环被控对象近似为积分环节电机轴系惯量所以Kp初值可以按下面公式给Kp_sp bw_sc · J / (1.5 · pp · flux)Ki_sp bw_sc · Kp_sp / 10Kp_sp算出来只是初值Ki先给一个比较小值然后在阶跃响应里看有没有稳态误差再看是否振荡一点点往上加。这是最实用的调法。速度环有一个常见问题如果不限制积分项速度阶跃时很大可能会超调上实物后很可能导致过流。所以速度环输出也就是iq_ref必须设置限幅。我脚本里限幅设的是±20A你可以根据电机额定电流填。5.3 从开环到双闭环的调试顺序我在带人的时候始终坚持“先开环、后电流环、再速度环”的顺序。千万别说你搭完模型就直接闭合速度环看波形那是拿运气当实力。第一步开环电压验证把PI都断开给定子施加一组低频三相正弦电压或者直接在dq轴给一个恒定电压看电机能不能转起来记录一下电流大小和波形是否正常。这能提前暴露测量模块、坐标变换极性、电机参数设置的问题。第二步电流环闭环固定iq_ref让id_ref0看id/iq能否快速跟踪。如果id和iq耦合严重或者跟踪不上先检查Park变换角度是否正确再看PI参数是否合理。第三步速度环闭环等电流环稳定后把速度环包上去。给一个速度阶跃看响应曲线是否平滑超调量是否可接受。如果转速发散优先检查反馈回来的是机械角度还是电角度以及有没有乘以极对数。我在实践中还发现一个规律很多人波形乱根子出在电流环还没稳定就开始调速度环到最后两个环都在震荡根本分不清是谁的锅。按顺序来能把问题范围缩到最小。6. 常见问题排查与经验6.1 典型现象速查表下面这张表是我这两年帮人看模型时遇到频率最高的问题按“现象→可能原因→排查方向”列出来现象可能原因排查方向仿真开始后数值直接NaN仿真步长太大把Powergui离散步长调小或降低电流环带宽电流发散、转速飞车电流反馈极性不对正反馈了检查Current Measurement方向交换相序或先开环验证电流环跟踪慢、有静差PI参数太小按KpL·bw、KiRs·bw重新计算必要时加大带宽转速环超调大速度环Ki太大或iq_ref限幅过大降Ki检查速度环输出限幅电流波形削顶、像梯形波SVPWM进入过调制区检查PI输出限幅u_max是否设置合理d轴电流不等于0Park变换角度错误检查是否漏乘极对数ppPWM模块无输出调制波极性/范围与PWM Generator不匹配确认占空比在0~1且PWM Generator输入模式选正确转矩始终为负相序接反交换A/B/C任意两相接线仿真发散但数值是NaN这个问题很多人第一反应是算法有bug。实际上多数原因是离散仿真步长比控制周期大导致PI连续积分时数值爆炸。先把Powergui采样时间调到Ts_ctrl/20甚至Ts_ctrl/50很多时候问题自己就消失了。6.2 不可忽视的细节与心得有几个细节我想单独拎出来说因为它们不起眼但影响极大。第一个是角度符号。Simscape电机模块输出的转角方向默认是逆时针为正如果你的坐标变换里sin/cos符号反了id/iq会颠倒或突变。如果你发现电流波形正常但id和iq完全交换了试试给Park变换里的角度加个负号通常能解决。第二个是母线电压在不同模型里的地位。Zero-order-hold、离散PI、PWM Generator都有各自的采样周期千万别把所有地方都填成仿真步长Ts_sim。控制周期是Ts_ctrlPWM载波周期也是Ts_ctrl而Ts_sim只是数值积分步长三者的关系要心里有数。否则会出现仿真波形看起来正确但硬件拿到代码后控制周期对不上的情况。第三个是给SVPWM函数输入的范围。u_alpha、u_beta是相电压指令的峰值不是线电压也不是标幺值。如果电流环PI输出没有限幅u_q可能会远超母线电压允许范围SVPWM过调制后电流波形失真。我建议在Simulink模型里把PI输出限幅设置成±(Vdc/sqrt(3))*0.95给SVPWM留出一点线性区余量。还有一点可能很多人没注意到MATLAB Function里不要写死Ts_ctrl最好作为参数传进来或者统一用基本工作区变量。我最初就是把开关频率硬编码在函数里后来改开关频率时总忘记同步模型跑出来的结果怎么都不对。现在我的习惯是函数里只用输入口所有参数都从初始化脚本传入尽量保持代码干净。最后一个经验是看波形。调PI的时候很多人都喜欢盯着转速波形看其实电流环调得好不好更该看的是d/q轴电流阶跃响应。我一般是先手动给iq_ref加一个方波看iq的上升时间、超调量、稳态误差然后用同样方法调id轴。等电流环在手感上已经合格了再去调速度环基本是一调一个准。写在调试台边上的话这篇文章写到这主体内容算是讲完了。我额外再分享一个个人习惯每次搭完一个新的PMSM仿真模型我都会先不接控制器直接给一个很小的开环三相正弦电压看电机能不能按预定方向缓慢旋转。这一步虽然看起来没什么技术含量却能在十分钟内确认动力链路、测量链路、坐标方向全都没问题之后再上双闭环就踏实很多。我那台电动工具项目的控制器最终上板前也是靠这层“笨办法”兜底才躲过了相序接反烧板子那类低级事故。希望这篇文章能帮你少踩一些我当年踩过的坑。