永磁同步电机数学模型：从坐标变换到控制基础
做电机控制这些年，我越来越觉得——数学模型就是控制算法的“地图”。没有这张地图，你连电机转不转得起来都说不清楚。今天咱们就把 PMSM 的数学底子彻底捋一遍。
2.1 为什么需要坐标变换？
三相永磁同步电机，定子上有三个绕组，通入三相对称电流。你想想看，三个变量互相耦合，控制起来多麻烦。我刚开始做项目时，直接在三相坐标系下写控制逻辑，结果调了两个月都没调顺。
后来才明白——坐标变换的本质，就是把三相静止坐标系下的交流量，变成两相旋转坐标系下的直流量。这样一来，控制就变成了调两个直流分量，简单太多了。
2.2 Clark 变换：从三相到两相
Clark 变换，也叫 3s/2s 变换。它把三相静止坐标系 (A, B, C) 映射到两相静止坐标系 (α, β)。
公式其实不复杂：
// 等幅值变换（最常用）
iα = ia
iβ = (ia + 2*ib) / √3
// 或者写成矩阵形式
[iα] [ 1 -1/2 -1/2 ] [ia]
[iβ] = [ 0 √3/2 -√3/2 ] [ib]
我的习惯：项目里我几乎只用等幅值变换。等功率变换虽然数学上更“干净”，但等幅值变换在调试时更直观——你给 1A 的电流指令，电机里就是 1A 的电流，不用换算。
Clark 变换的物理意义很简单：把三个绕组的磁动势，等效成两个垂直绕组产生的磁动势。说白了，就是“降维打击”。
2.3 Park 变换：从静止到旋转
Clark 变换做完，我们得到了 αβ 轴上的交流量。但交流量还是不好控——PID 调交流量，带宽永远追不上。
Park 变换（2s/2r 变换）就是来解决这个问题的。它把 αβ 静止坐标系，旋转到与转子磁极同步的 dq 旋转坐标系。
// Park 变换公式
id = iα * cos(θ) + iβ * sin(θ)
iq = -iα * sin(θ) + iβ * cos(θ)
// 其中 θ 是转子电角度
做完 Park 变换，id 和 iq 就变成了直流量。这时候你再上 PI 控制器，效果立竿见影。
我曾经踩过的坑：Park 变换用的角度，必须是电角度，不是机械角度。极对数不为 1 时，角度要乘以极对数。我第一次做 4 对极电机时忘了这茬，结果电流波形乱成一团，查了三天才找到原因。
2.4 电压方程：电机怎么“响应”你的指令？
在 dq 坐标系下，PMSM 的电压方程长这样：
ud = Rs * id + Ld * (did/dt) - ωe * Lq * iq
uq = Rs * iq + Lq * (diq/dt) + ωe * (Ld * id + ψf)
这里面：
- Rs：定子电阻，铜损的来源
- Ld, Lq：d 轴和 q 轴电感，凸极率的关键
- ωe：电角速度，ωe = 极对数 × 机械角速度
- ψf：永磁体磁链，电机的“底力”
你看这个方程，有个很有意思的地方——交叉耦合项。uq 方程里有个 ωe * Ld * id，ud 方程里有个 -ωe * Lq * iq。这意味着 d 轴和 q 轴不是独立的，转速越高，耦合越强。
关键理解：电压方程告诉我们——高速时，你给 uq 加电压，一部分会用来克服反电动势 (ωe * ψf)，一部分要抵消 d 轴耦合。这就是为什么高速时电压利用率会下降，也是弱磁控制的物理根源。
2.5 转矩方程：电机能出多少力？
转矩方程是 PMSM 控制的核心目标：
Te = 1.5 * p * [ ψf * iq + (Ld - Lq) * id * iq ]
拆开来看：
- ψf * iq：永磁转矩，主转矩分量。iq 越大，转矩越大。
- (Ld - Lq) * id * iq：磁阻转矩。只有 Ld ≠ Lq（凸极电机）时才有。
我刚开始做内置式 PMSM（IPMSM）时，一直不理解为什么还要用 id。后来才明白——对于 IPMSM，Ld < Lq，给 id 一个负值，磁阻转矩就变成正的。这就是 MTPA 的底层逻辑。
| 电机类型 | Ld vs Lq | 磁阻转矩 | 控制策略 |
|---|---|---|---|
| 表贴式 SPMSM | Ld ≈ Lq | 几乎为零 | id = 0 控制 |
| 内置式 IPMSM | Ld < Lq | 可利用 | MTPA + 弱磁 |
2.6 MTPA：最大转矩电流比
MTPA（Maximum Torque Per Ampere），说白了就是——用最小的电流，出最大的转矩。
对于 SPMSM，MTPA 就是 id = 0，很简单。但对于 IPMSM，需要找到 id 和 iq 的最优分配。
MTPA 的数学条件：
∂Te / ∂id = 0 （在电流幅值约束下）
推导结果：
id = -ψf / (2 * (Ld - Lq)) - sqrt( ψf² / (4 * (Ld - Lq)²) + iq² )
这个公式看着复杂，但实际工程中，我一般用查表法——离线算好 id_ref 和 iq_ref 的对应关系，在线直接查。
我的经验：MTPA 查表法比在线计算稳定得多。在线计算涉及开方和除法，定点 DSP 上跑容易出精度问题。我有个项目就是在线算 MTPA，结果低速时转矩抖动，换成查表后一次搞定。
2.7 弱磁控制：突破电压极限
电机转速升高，反电动势 (ωe * ψf) 越来越大。当它接近母线电压时，电流就控不住了。这时候就需要弱磁。
弱磁的本质：给 id 一个负值，抵消永磁体磁链，降低反电动势。
弱磁的电压约束：
ud² + uq² ≤ (Vdc / √3)² （SVPWM 最大不失真电压）
进入弱磁区后，id 会越来越负，iq 也要相应减小。这时候转矩会下降——这就是恒功率区的特性。
避坑指南：我曾经在一个项目中，弱磁深度太大，id 负到 -3 倍额定电流。结果永磁体发生了不可逆退磁，电机直接报废。弱磁不是无限度的，一定要留退磁安全裕量。
2.8 知识体系总览
下面这张图，把本章的核心逻辑串起来了：
从这张图可以看得很清楚：三相电流经过 Clark → Park 变换，变成 dq 轴直流量，然后基于电压方程和转矩方程，设计 id/iq 的指令分配策略。MTPA 和弱磁，本质上就是在不同转速下，如何分配 id 和 iq 的问题。
2.9 本章小结
嗯，内容不少，但核心就这几条：
- Clark 变换：三相变两相，降维
- Park 变换：交流变直流，好控
- 电压方程：告诉你电压和电流、转速的关系
- 转矩方程：告诉你电流怎么产生转矩
- MTPA：低速区，用最小电流出最大转矩
- 弱磁：高速区，牺牲转矩换转速
这些是 PMSM 控制的“内功心法”。后面讲三环整定时，你会反复用到这些概念。我个人建议你把电压方程和转矩方程背下来——调试时随手就能写出耦合补偿项，那种感觉真的很爽。