4. PMSM数学模型
永磁同步电机（PMSM）的数学模型是整个FOC（磁场定向控制）算法体系的基石。无论是电流环PI参数整定、SVPWM调制策略，还是MTPA（最大转矩电流比）与弱磁控制，其本质都是在求解或利用电机的电压、磁链与转矩方程。本章将系统梳理PMSM在三相静止坐标系（ABC）、两相静止坐标系（αβ）以及两相旋转坐标系（dq）下的数学模型，重点剖析坐标变换的物理意义与方程形式的变化逻辑，为后续的控制器设计打下坚实的理论基础。
4.1 三相静止坐标系（ABC）下的PMSM方程
PMSM的定子绕组在空间上互差120°电角度，当通入三相对称交流电时，会产生旋转磁场。在三相静止坐标系下建立数学模型是最直接的方式，但也是最复杂的，因为定子绕组之间存在互感耦合，且转子位置角θ会随时间变化，导致方程呈现非线性、时变特性。
4.1.1 磁链方程
三相定子磁链由两部分组成：定子电流产生的电枢反应磁链和永磁体产生的励磁磁链。以A相为例，其磁链表达式为：
ψa = Laa·ia + Mab·ib + Mac·ic + ψfa(θe)其中，Laa为A相自感，Mab、Mac为A相与B、C相的互感，ψfa(θe)为永磁体在A相绕组中产生的磁链分量，其表达式为ψfa = ψf·cos(θe)。由于转子位置θe随时间变化，ψfa也是时变的。同理可写出B相和C相的磁链方程。将三相磁链方程写成矩阵形式，可表示为：
[ψabc] = [Ls]·[iabc] + ψf·[cos(θe), cos(θe-2π/3), cos(θe+2π/3)]ᵀ其中[Ls]为电感矩阵，包含自感和互感项。对于表贴式电机，Laa≈Lbb≈Lcc，互感也近似相等；对于内置式电机，电感随转子位置变化，矩阵更为复杂。
4.1.2 电压方程
根据基尔霍夫电压定律，三相定子绕组的电压方程可写为：
ua = Rs·ia + dψa/dt
ub = Rs·ib + dψb/dt
uc = Rs·ic + dψc/dt将磁链方程代入电压方程后，可以看到电压由三部分组成：电阻压降、电感压降（自感和互感）以及反电动势（永磁体旋转产生的运动电动势）。由于电感矩阵中存在转子位置θe的余弦项，电压方程是强耦合的时变微分方程组，直接用于控制器设计非常困难。
4.1.3 转矩方程
电磁转矩可以通过磁场储能对机械角位移求偏导得到，也可以从能量守恒角度推导。三相坐标系下的转矩方程为：
Te = (1/2)·p·(iabc)ᵀ·(∂[Ls]/∂θe)·(iabc) + p·ψf·(-ia·sinθe - ib·sin(θe-2π/3) - ic·sin(θe+2π/3))该方程包含两项：第一项为磁阻转矩（由电感随转子位置变化引起），第二项为永磁转矩（由永磁体与定子电流相互作用产生）。对于表贴式电机，第一项近似为零，转矩主要由第二项决定。
4.2 两相静止坐标系（αβ）下的PMSM方程
Clarke变换将三相ABC坐标系投影到两相正交的αβ坐标系中，其核心思想是用两个正交绕组等效替代三个互差120°的绕组。Clarke变换分为等幅值变换和等功率变换两种形式，工程中FOC算法通常采用等幅值变换，因为其物理意义更直观（变换后电流幅值与相电流幅值相等）。
4.2.1 Clarke变换矩阵
等幅值Clarke变换的矩阵形式为：
[iα] [ 1 -1/2 -1/2 ] [ia]
[iβ] = [ 0 √3/2 -√3/2 ] [ib]
[i0] [ 1/2 1/2 1/2 ] [ic]对于无中线系统，零序分量i0通常为零，可忽略。变换后，α轴与A相轴重合，β轴超前α轴90°电角度。
4.2.2 αβ坐标系下的磁链与电压方程
经过Clarke变换后，磁链方程简化为：
ψα = Ls·iα + ψf·cos(θe)
ψβ = Ls·iβ + ψf·sin(θe)其中Ls为等效两相绕组的电感（对于表贴式电机，Ls = Laa - Mab）。电压方程为：
uα = Rs·iα + Ls·(diα/dt) - ωe·ψf·sin(θe)
uβ = Rs·iβ + Ls·(diβ/dt) + ωe·ψf·cos(θe)可以看到，αβ坐标系下的方程仍然包含转子位置θe的正弦和余弦项，但电感矩阵已经对角化，不再存在互感耦合。反电动势项（-ωe·ψf·sinθe和ωe·ψf·cosθe）是转子位置的函数，这为后续的转子位置观测（如滑模观测器、龙伯格观测器）提供了数学基础。
4.2.3 αβ坐标系下的转矩方程
电磁转矩在αβ坐标系下可简洁地表示为：
Te = 1.5·p·(ψα·iβ - ψβ·iα)将磁链方程代入，可进一步展开为：
Te = 1.5·p·[Ls·(iα·iβ - iβ·iα) + ψf·(cosθe·iβ - sinθe·iα)] = 1.5·p·ψf·(cosθe·iβ - sinθe·iα)该方程表明，在αβ坐标系下，转矩与电流和转子位置的正余弦函数直接相关，仍然是一个交流量。但相比ABC坐标系，方程形式已大幅简化，且αβ轴电流的幅值和相位可以直接反映转矩的大小和方向。
4.3 两相旋转坐标系（dq）下的PMSM方程
Park变换将αβ坐标系进一步旋转到与转子同步旋转的dq坐标系中，其中d轴与转子永磁体磁链方向对齐，q轴超前d轴90°电角度。经过Park变换后，交流量变为直流量，这是FOC实现解耦控制的关键步骤。
4.3.1 Park变换矩阵
Park变换的矩阵形式为：
[id] [ cosθe sinθe ] [iα]
[iq] = [ -sinθe cosθe ] [iβ]其中θe为转子电角度，由编码器或霍尔传感器实时获取。变换后，id为励磁电流分量（与永磁体磁链方向一致），iq为转矩电流分量（与永磁体磁链正交）。
4.3.2 dq坐标系下的磁链方程
经过Park变换后，磁链方程变为：
ψd = Ld·id + ψf
ψq = Lq·iq其中Ld和Lq分别为d轴和q轴电感。对于表贴式电机（SPM），Ld≈Lq；对于内置式电机（IPM），Ld<Lq，这是产生磁阻转矩的根本原因。磁链方程在dq坐标系下完全解耦，不再包含转子位置信息，这是FOC能够实现线性控制的核心优势。
4.3.3 dq坐标系下的电压方程
电压方程是FOC电流环设计的基础，其标准形式为：
ud = Rs·id + Ld·(did/dt) - ωe·Lq·iq
uq = Rs·iq + Lq·(diq/dt) + ωe·(Ld·id + ψf)方程中出现了两项耦合项：-ωe·Lq·iq（d轴电压受q轴电流影响）和ωe·(Ld·id + ψf)（q轴电压受d轴电流和反电动势影响）。这两项是FOC中"前馈解耦"需要补偿的对象。当电机稳态运行时，did/dt和diq/dt为零，电压方程简化为：
ud = Rs·id - ωe·Lq·iq
uq = Rs·iq + ωe·(Ld·id + ψf)此时，uq中的ωe·ψf项即为反电动势，在高速运行时该值很大，需要足够的母线电压来克服。
4.3.4 dq坐标系下的转矩方程
电磁转矩在dq坐标系下的表达式为：
Te = 1.5·p·[ψf·iq + (Ld - Lq)·id·iq]该方程包含两项：第一项1.5·p·ψf·iq为永磁转矩，与iq成正比；第二项1.5·p·(Ld-Lq)·id·iq为磁阻转矩，与id和iq的乘积成正比。对于表贴式电机（Ld=Lq），磁阻转矩为零，转矩仅由iq控制；对于内置式电机（Ld<Lq），可以通过注入负的id电流来利用磁阻转矩，提高转矩密度，这就是MTPA控制的理论依据。
4.4 三种坐标系方程对比与工程选择
为了更直观地理解三种坐标系下PMSM数学模型的特点，下表从多个维度进行了对比：
| 对比维度 | ABC坐标系 | αβ坐标系 | dq坐标系 |
|---|---|---|---|
| 变量性质 | 交流时变量 | 交流时变量 | 直流量（稳态） |
| 电感矩阵 | 含自感和互感，随θe变化 | 对角化，常数 | 对角化，常数（Ld/Lq） |
| 转子位置依赖 | 强依赖（cos/sin项） | 依赖（反电动势项） | 仅通过ωe体现 |
| 耦合程度 | 强耦合 | 弱耦合（仅反电动势） | 存在交叉耦合项 |
| 转矩方程 | 复杂，含磁阻和永磁项 | Te=1.5p(ψα·iβ-ψβ·iα) | Te=1.5p[ψf·iq+(Ld-Lq)·id·iq] |
| 控制适用性 | 不适用（仅物理分析） | 适用于无传感器观测、DTC | 适用于FOC电流环/速度环 |
4.5 案例拆解：表贴式PMSM的FOC电流环方程应用
下面通过一个实际工程案例，展示dq坐标系方程在FOC电流环设计中的具体应用。
背景与约束：某伺服驱动器采用表贴式PMSM，额定参数为：Rs=0.5Ω，Ld=Lq=1.2mH，ψf=0.08Wb，极对数p=4，额定转速3000rpm（对应电角速度ωe=2π×3000/60×4=1256.6rad/s），母线电压310VDC。控制器采用id=0控制策略，目标iq=5A。
初始现象：在调试过程中发现，当电机转速超过2000rpm时，电流环出现饱和，实际iq无法跟踪指令值，且ud指令出现明显波动。
分析过程：根据dq电压方程，在id=0控制下，稳态电压指令为：
ud = -ωe·Lq·iq = -1256.6 × 0.0012 × 5 = -7.54V
uq = Rs·iq + ωe·ψf = 0.5×5 + 1256.6×0.08 = 2.5 + 100.5 = 103V合成电压幅值为√(ud²+uq²)≈103.3V。当转速升高到3000rpm时，uq中的反电动势项ωe·ψf达到100.5V，加上电阻压降后接近103V。考虑到SVPWM在母线电压310V下的最大线性输出相电压幅值为310/√3≈179V，理论上仍有裕量。但实际调试中发现，由于电流环PI参数未考虑前馈解耦，动态过程中ud和uq出现振荡，导致电压利用率下降。
解决方案：在电流环中加入前馈解耦项，将PI输出与解耦项叠加：
ud_ref = (Kp + Ki/s)·(id_ref - id) - ωe·Lq·iq
uq_ref = (Kp + Ki/s)·(iq_ref - iq) + ωe·(Ld·id + ψf)同时，将PI参数按照"带宽法"整定：设定电流环带宽为1000Hz，则Kp = Ls×带宽 = 0.0012×6283 = 7.54，Ki = Rs×带宽 = 0.5×6283 = 3141.5。修改后，ud指令在稳态时稳定在-7.5V左右，uq指令稳定在103V，电流环动态响应时间从原来的5ms缩短至1ms以内。
结果对比：加入前馈解耦后，在3000rpm满载运行时，电流环不再饱和，iq跟踪误差从±0.8A减小到±0.15A，系统效率提升约3%。
复盘经验：dq坐标系下的电压方程不仅是理论公式，更是工程调试的"地图"。在调试FOC电流环时，务必先根据电机参数和运行工况计算理论电压指令，再与示波器实测值对比，这样可以快速定位是PI参数问题、解耦补偿问题还是SVPWM调制问题。
4.6 总结与工程要点
PMSM的数学模型从ABC到αβ再到dq，是一个逐步简化和解耦的过程。ABC坐标系下的方程物理意义清晰但难以控制；αβ坐标系下的方程适合转子位置观测和直接转矩控制；dq坐标系下的方程实现了磁链和转矩的解耦，是FOC控制的核心。在实际工程中，需要重点掌握以下几点：
- 变换矩阵的基准：等幅值变换与等功率变换的系数差异，以及对应的转矩方程系数调整。
- 电感参数辨识：Ld和Lq的准确值直接影响解耦效果和PI参数整定，建议通过离线辨识（如注入高频信号）或在线辨识获取。
- 反电动势补偿：高速运行时，uq中的ωe·ψf项占主导，必须确保母线电压足够，否则需要弱磁控制。
- 耦合项处理：前馈解耦是标准做法，但在低成本MCU上可简化为查表补偿，以节省计算资源。