电机数学模型：三相坐标系、Clark变换、Park变换、dq轴模型、转矩方程

无刷电机FOC控制的本质，是将电机视为一个多变量、强耦合、非线性的被控对象，并通过坐标变换将其解耦为等效的直流电机模型。理解从三相静止坐标系到两相旋转坐标系的数学映射，是掌握磁场定向控制的理论基石。本节将系统梳理三相坐标系下的电机方程、Clark变换与Park变换的推导逻辑、dq轴等效电路模型以及转矩方程的物理意义，帮助学习者建立从物理电机到数学模型的完整认知链路。

3.1 三相静止坐标系下的电机方程

在自然坐标系（A-B-C三相）下，永磁同步电机（PMSM）的定子电压方程可以表示为：

u_a = Rs·i_a + dψ_a/dt
u_b = Rs·i_b + dψ_b/dt
u_c = Rs·i_c + dψ_c/dt

其中，ψ_a、ψ_b、ψ_c为三相定子磁链，包含定子电流产生的磁链和永磁体转子产生的耦合磁链。由于三相绕组在空间上互差120°电角度，且转子位置随时间变化，导致磁链方程中存在与转子位置θe相关的正弦、余弦项，使得方程呈现强耦合和非线性特征。直接在三相坐标系下设计控制器极为困难，因为每一相电压不仅受本相电流影响，还受到其他两相电流和转子位置的影响。

3.2 Clark变换：从三相静止到两相静止

Clark变换（又称3s/2s变换）将三相静止坐标系（A-B-C）映射到两相静止坐标系（α-β），其核心思想是保持磁动势等效。对于星形连接且无中性线引出的电机，三相电流满足i_a + i_b + i_c = 0，因此变换矩阵可以简化为：

[ iα ]   [ 1    -1/2   -1/2 ] [ ia ]
[ iβ ] = [ 0    √3/2  -√3/2 ] [ ib ]
                    [ ic ]

采用等幅值变换时，变换矩阵前乘以系数2/3。变换后，α轴与A轴重合，β轴超前α轴90°电角度。Clark变换的物理意义在于：将空间上互差120°的三个绕组等效为空间上正交的两个绕组，从而将三相变量简化为两相变量，但此时变量仍然是交流量，随转子位置变化。

3.3 Park变换：从两相静止到两相旋转

Park变换（又称2s/2r变换）将α-β静止坐标系映射到d-q旋转坐标系，旋转角速度与转子电角速度ωe同步。变换矩阵为：

[ id ]   [ cosθe   sinθe ] [ iα ]
[ iq ] = [ -sinθe  cosθe ] [ iβ ]

经过Park变换后，d轴与转子永磁体磁链方向对齐，q轴超前d轴90°电角度。在稳态运行时，id和iq变为直流量，从而可以通过PI控制器实现无静差调节。d轴电流分量id控制励磁磁场，q轴电流分量iq控制转矩输出，实现了磁场与转矩的解耦控制。

3.4 dq轴等效电路模型

在dq旋转坐标系下，PMSM的电压方程可以写为：

ud = Rs·id + Ld·did/dt - ωe·Lq·iq
uq = Rs·iq + Lq·diq/dt + ωe·Ld·id + ωe·ψf

其中，Ld和Lq分别为d轴和q轴电感，ψf为永磁体磁链。对于表贴式PMSM（SPMSM），Ld = Lq；对于内置式PMSM（IPMSM），Ld < Lq。方程中的ωe·Lq·iq和ωe·Ld·id为交叉耦合项，ωe·ψf为反电动势项。在FOC控制中，通常通过前馈补偿的方式消除这些耦合项，使d轴和q轴可以独立控制。

3.5 转矩方程与FOC控制策略

PMSM的电磁转矩方程为：

Te = 1.5·p·[ψf·iq + (Ld - Lq)·id·iq]

其中，p为极对数。第一项ψf·iq为永磁转矩，第二项(Ld-Lq)·id·iq为磁阻转矩。对于表贴式电机（Ld=Lq），磁阻转矩为零，转矩仅与iq成正比，因此采用id=0控制策略即可获得线性转矩特性。对于内置式电机（Ld<Lq），可以利用负的id产生磁阻转矩，在相同电流幅值下获得更大转矩，即最大转矩电流比（MTPA）控制。

3.6 坐标变换的离散化实现与工程注意事项

在数字控制器中，坐标变换以离散形式执行。典型FOC控制周期为10kHz-20kHz，每个周期内依次执行：电流采样→Clark变换→Park变换→PI调节→逆Park变换→SVPWM调制。需要注意以下几点：

• 角度对齐：Park变换使用的θe必须与当前采样时刻的转子位置一致，否则会产生角度延迟误差。高速运行时，建议在采样时刻补偿角度延迟：θe_comp = θe + ωe·T_sample。
• 电流采样同步：三相电流采样应在PWM载波的中心点进行，此时电流纹波最小，采样值最接近平均值。
• 标幺化处理：在定点DSP中，建议将电压、电流、磁链等物理量标幺化到[-1,1]区间，避免溢出并提高运算精度。
• 变换矩阵预计算：sinθe和cosθe的计算耗时较长，可查表或使用CORDIC算法加速。对于高开关频率应用，可预先计算并存储角度表。