4. 永磁同步电机模型：PMSM的数学模型
聊到永磁同步电机（PMSM）的数学模型，很多初学者第一反应就是“dq变换好难”。其实说白了，这就是一个把复杂问题简单化的过程。我在刚入行时也被这些变换绕得晕头转向，直到有一次调试一个高速主轴电机，才真正理解了这些变换的价值。
4.1 为什么需要坐标变换？
先问一个问题：三相交流电机的电压方程长什么样？
在三相静止坐标系（ABC坐标系）下，PMSM的电压方程是一组时变系数的微分方程。因为转子在转，定子和转子之间的互感是随时间变化的。你想想看，要实时求解这样的方程，对控制器的算力要求太高了。
我早期做的一个项目，直接用三相模型做电流环，结果CPU算到冒烟，电流环带宽还是上不去。后来换了dq轴模型，问题迎刃而解。
核心思想：通过坐标变换，把时变系数的交流电机模型，变成常系数直流电机模型。这样控制起来就跟控制直流电机一样简单。
4.2 Clark变换：从ABC到αβ
Clark变换，也叫3/2变换。它的任务是把三相静止坐标系（A、B、C）变换到两相静止坐标系（α、β）。
数学表达式很简单：
// Clark变换（等幅值变换）
iα = iA
iβ = (iA + 2*iB) / √3
// 或者用矩阵形式
[iα] [ 1 -1/2 -1/2 ] [iA]
[iβ] = [ 0 √3/2 -√3/2 ] [iB]
[i0] [ 1/2 1/2 1/2 ] [iC]这里有个细节要注意：变换系数有等幅值和等功率两种选择。我个人习惯用等幅值变换，因为电流环的PI参数整定更直观。
小技巧：在实际代码中，如果电机是星形接法且中性点不接地，零序分量i0可以忽略。这样Clark变换只需要两个电流传感器就能完成。
4.3 Park变换：从αβ到dq
Park变换才是真正的“魔法”。它把两相静止坐标系（αβ）旋转到与转子同步旋转的坐标系（dq）。
变换公式：
// Park变换
id = iα * cos(θ) + iβ * sin(θ)
iq = -iα * sin(θ) + iβ * cos(θ)其中θ是转子位置角（电角度）。
为什么要这样做？因为变换之后，基波分量变成了直流量。你想想看，控制直流量和控制交流量，哪个更容易？
我曾经遇到过一个案例：现场工程师反馈电机低速运行时噪音大。检查后发现是Park变换中的角度补偿没做好，导致dq轴电流有纹波。把角度补偿加上后，问题就解决了。
注意：Park变换需要准确的转子位置信息。如果编码器分辨率不够或者安装有偏差，dq轴电流会存在耦合振荡。我建议在调试时先用示波器观察id和iq的波形，如果发现有正弦波动，先检查角度信号。
4.4 完整的坐标变换流程
把Clark和Park串起来，就是完整的坐标变换链路：
三相电流 → Clark变换 → 两相静止电流 → Park变换 → 两相旋转电流
iA,iB,iC iα,iβ id,iq反过来，从dq轴电压到三相电压的逆变换也一样：
dq轴电压 → 逆Park变换 → αβ轴电压 → 逆Clark变换 → 三相电压
ud,uq uα,uβ uA,uB,uC下面这张图展示了整个变换的流程：
4.5 PMSM的电压方程与磁链方程
经过dq变换后，PMSM的数学模型变得非常简洁。
磁链方程：
ψd = Ld * id + ψf
ψq = Lq * iq其中：
- ψd、ψq：d轴和q轴的磁链
- Ld、Lq：d轴和q轴的电感（对于表贴式PMSM，Ld = Lq）
- ψf：永磁体产生的磁链（常数）
电压方程：
ud = Rs * id + dψd/dt - ωe * ψq
uq = Rs * iq + dψq/dt + ωe * ψd其中：
- Rs：定子电阻
- ωe：电角速度（ωe = np * ωm，np是极对数，ωm是机械角速度）
把磁链方程代入电压方程，得到更常用的形式：
ud = Rs * id + Ld * did/dt - ωe * Lq * iq
uq = Rs * iq + Lq * diq/dt + ωe * (Ld * id + ψf)关键理解：从电压方程可以看出，d轴和q轴之间存在耦合项（-ωe*Lq*iq 和 ωe*Ld*id）。这就是为什么在高速运行时，需要做解耦控制。我在调试一个2kW的伺服驱动器时，没做解耦，结果高速时电流环直接发散。加上前馈解耦后，稳如泰山。
4.6 转矩方程
PMSM的电磁转矩方程是：
Te = 1.5 * np * [ψf * iq + (Ld - Lq) * id * iq]这个方程包含两部分：
| 转矩分量 | 表达式 | 说明 |
|---|---|---|
| 永磁转矩 | 1.5 * np * ψf * iq | 由永磁体与q轴电流相互作用产生，与iq成正比 |
| 磁阻转矩 | 1.5 * np * (Ld - Lq) * id * iq | 由d、q轴电感差异产生，仅在内嵌式PMSM中存在 |
对于表贴式PMSM（Ld = Lq），磁阻转矩为零，转矩只与iq有关：
Te = 1.5 * np * ψf * iq这就是为什么在表贴式电机中，我们通常让id=0来控制。而对于内嵌式电机，可以利用磁阻转矩来提高效率，这就是MTPA（最大转矩电流比）控制的基础。
经验之谈：我在做电动汽车驱动项目时，内嵌式电机用MTPA控制，同样的电流下转矩提升了约15%。但要注意，id为负值会削弱永磁体磁场，过大的id可能导致永磁体退磁。这个度要把握好。
4.7 机械运动方程
完整的PMSM模型还需要包含机械部分：
Te - TL = J * dωm/dt + B * ωm其中：
- TL：负载转矩
- J：转动惯量
- B：阻尼系数
- ωm：机械角速度
结合前面的电气方程，就构成了PMSM的完整数学模型。这个模型是矢量控制的基础，也是仿真建模的依据。
注意：实际电机参数会随工况变化。比如电阻Rs会随温度升高而增大，电感Ld、Lq会随电流增大而饱和。我在做参数辨识时发现，同一个电机在冷态和热态下，Rs能差30%以上。所以高性能伺服驱动都会做在线参数辨识和补偿。
4.8 本章小结
总结一下PMSM数学模型的核心要点：
- Clark变换：把三相静止变两相静止，减少变量数量
- Park变换：把两相静止变两相旋转，把交流量变直流量
- 磁链方程：ψd = Ld*id + ψf，ψq = Lq*iq
- 电压方程：包含电阻压降、电感压降和反电动势三项
- 转矩方程：永磁转矩 + 磁阻转矩