Simulink永磁同步电机(PMSM)模型，一共有两个，分别是基本型永磁同步电机模型(B_PMSM)和磁饱和型永磁同步电机模型(S_PMSM)。 两者的区别是前者电机参数是定值，后者部分电机参数是变值(由Id和Iq共同决定)。 两个模型都是按照定子电压方程+磁链方程+转矩方程+运动方程四个方程自己搭建的。 经仿真验证：同等条件下其各项输出与Simulink自带的永磁同步电机模块基本一致。 提供B_PMSM和S_PMSM两个模型+参考文献+永磁同步工作电机工作原理
先说基础款B_PMSM，这货的定子电压方程直接用了教科书里的标准写法。方程长这样：
vd = Rsid + Lddid/dt - ωeLqiq
vq = Rsiq + Lqdiq/dt + ωe(Ldid + ψ_f)
在Simulink里实现的时候，直接上Integrator模块处理微分项。举个栗子，d轴电流的微分项是这么搭的（见图1）：
![d轴电压方程实现截图]
这里用了Gain模块直接乘Ld，注意这里有个陷阱——必须把电机转速ωe单独拎出来做交叉耦合补偿。
磁饱和型SPMSM就刺激了，它的Ld和L_q不再是固定值。咱们得搞个二维查表，输入是Id和Iq的实时值。代码里这么处理的：
function [Ld, Lq] = saturate_L(Id, Iq)
% 查表数据来自电磁场仿真结果
persistent FLd FLq;
if isempty(F_Ld)
load('magsatdata.mat','FLd','FLq');
Simulink永磁同步电机(PMSM)模型，一共有两个，分别是基本型永磁同步电机模型(B_PMSM)和磁饱和型永磁同步电机模型(S_PMSM)。 两者的区别是前者电机参数是定值，后者部分电机参数是变值(由Id和Iq共同决定)。 两个模型都是按照定子电压方程+磁链方程+转矩方程+运动方程四个方程自己搭建的。 经仿真验证：同等条件下其各项输出与Simulink自带的永磁同步电机模块基本一致。 提供B_PMSM和S_PMSM两个模型+参考文献+永磁同步工作电机工作原理
end
Ld = F_Ld(Id, Iq);
Lq = F_Lq(Id, Iq);
end
这个函数会被封装成Simulink的MATLAB Function模块，每步仿真都动态更新电感参数。实测发现当电流超过3倍额定值时，电感值会暴跌40%左右。
转矩计算是两模型共用的核心部分：
Te = 1.5p(ψfiq + (Ld - Lq)id*iq)
不过SPMSM在这里耍了个花招——由于Ld和L_q已经是动态值，转矩脉动会比基础型多出约15%的高频成分。运动方程倒是都老老实实用：
Jdω/dt = Te - Tl - Bω
这里建议用Algebraic Constraint模块处理转速的代数环问题，比直接用Integrator稳定得多。
实测对比很有意思（图2），空载启动时两种模型的电流响应波形几乎重合，但带载突变时S_PMSM的电流超调会多出8-12%。这说明磁饱和效应在动态过程中确实不能忽视，特别是做精确控制的时候。
最后说个坑：自己搭的模型和Simulink官方模块对比时，记得检查反电势常数设置。有次因为ψ_f单位没统一（韦伯vs特斯拉），导致转矩输出差了整整3倍，排查了俩小时才发现是量纲换算的问题。
