2、电机模型基础：永磁同步电机(PMSM)的数学模型、坐标变换(Clark/Park)、电压方程与磁链方程

无感FOC算法的核心在于通过精确的电机数学模型，在无位置传感器的条件下估算转子位置与速度。要理解滑模观测器(SMO)的设计原理，首先必须扎实掌握永磁同步电机(PMSM)在旋转磁场下的电气特性与数学描述。本章将系统梳理PMSM的物理结构、三相静止坐标系下的自然模型、Clark变换与Park变换的数学推导，以及最终在旋转dq坐标系下简化的电压方程与磁链方程。这些方程是后续构建滑模观测器、设计反电动势提取算法的基石。

2.1 永磁同步电机的物理结构与电磁特性

永磁同步电机（PMSM）的转子采用永磁体励磁，定子为三相对称绕组。根据永磁体在转子上的安装位置，PMSM可分为表贴式（SPMSM）与内置式（IPMSM）两种。表贴式电机的交轴电感Lq与直轴电感Ld近似相等（Ld≈Lq），而内置式电机由于磁路不对称，通常Lq>Ld，从而产生磁阻转矩。在无感FOC算法中，滑模观测器主要依赖反电动势（Back-EMF）信息，而反电动势的大小与转子磁链ψf和电角速度ωe成正比。对于表贴式电机，由于Ld=Lq，其数学模型更为简洁，反电动势在α-β轴系中呈现为正弦波形，这为SMO的设计提供了便利。

从电磁关系来看，定子绕组中通入三相对称正弦电流时，会产生旋转的定子磁场。该磁场与转子永磁体产生的恒定磁场相互作用，产生电磁转矩。当转子旋转时，定子绕组中会感应出反电动势，其相位与转子位置直接相关。因此，通过检测反电动势的相位，即可间接获得转子位置信息。这正是无感FOC中滑模观测器的核心思想——利用电压方程重构反电动势，进而提取转子位置与速度。

2.2 三相静止坐标系下的PMSM数学模型

在三相静止坐标系（A-B-C）下，PMSM的电压方程可以表示为：

u_A = R_s·i_A + dψ_A/dt
u_B = R_s·i_B + dψ_B/dt
u_C = R_s·i_C + dψ_C/dt

其中，R_s为定子相电阻，ψ_A、ψ_B、ψ_C为各相绕组的总磁链。磁链方程则更为复杂，它不仅包含定子电流产生的自感磁链和互感磁链，还包含转子永磁体产生的励磁磁链。以矩阵形式表示：

ψ_ABC = L_ABC·i_ABC + ψ_f·F(θ_e)

式中，L_ABC为3×3电感矩阵，其元素随转子位置θ_e变化，呈现周期性；ψ_f为永磁体磁链幅值；F(θ_e)为包含sin/cos函数的相位向量。这种强耦合、时变的数学模型在工程计算中极为不便，无法直接用于控制器设计。因此，必须引入坐标变换，将三相静止坐标系下的方程转换到两相旋转坐标系中，以实现解耦。

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

Clark变换（又称3s/2s变换）将三相静止坐标系（A-B-C）中的变量映射到两相静止坐标系（α-β）中。其数学表达式为：

[i_α]   [ 1    -1/2   -1/2 ] [i_A]
[i_β] = [ 0    √3/2  -√3/2 ] [i_B]
              [ 2/3 系数可选 ] [i_C]

这里存在两种常用的变换系数：等幅值变换（系数2/3）和等功率变换（系数√(2/3)）。在FOC控制中，通常采用等幅值变换，因为其物理意义直观——变换后α-β轴电流幅值与三相电流幅值相等。需要注意的是，Clark变换本身不涉及转子位置信息，它只是将三相绕组在空间上的120°分布关系映射为正交的α-β轴系。经过Clark变换后，原本耦合的三相变量被分解为两个正交分量，但这两个分量仍然是交流量，其频率与电频率相同。

在滑模观测器设计中，α-β轴系下的电压方程是核心。对于表贴式PMSM（Ld=Lq=Ls），α-β轴电压方程可简化为：

u_α = R_s·i_α + L_s·di_α/dt + e_α
u_β = R_s·i_β + L_s·di_β/dt + e_β

其中，e_α和e_β为反电动势分量，其表达式为：

e_α = -ψ_f·ω_e·sin(θ_e)
e_β =  ψ_f·ω_e·cos(θ_e)

从上述公式可以看出，反电动势的相位直接包含转子位置信息θ_e。SMO的核心任务就是通过滑模面设计，从电压方程中准确提取e_α和e_β，进而通过反正切函数或锁相环（PLL）计算出转子位置与速度。

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

Park变换（又称2s/2r变换）将两相静止坐标系（α-β）中的变量转换到随转子同步旋转的d-q坐标系中。变换矩阵为：

[i_d]   [ cos(θ_e)  sin(θ_e) ] [i_α]
[i_q] = [-sin(θ_e)  cos(θ_e) ] [i_β]

经过Park变换后，交流量变为直流量，这为PI控制器提供了无静差控制的基础。在d-q坐标系下，PMSM的电压方程和磁链方程呈现出简洁的解耦形式：

u_d = R_s·i_d + L_d·di_d/dt - ω_e·L_q·i_q
u_q = R_s·i_q + L_q·di_q/dt + ω_e·(L_d·i_d + ψ_f)

磁链方程为：

ψ_d = L_d·i_d + ψ_f
ψ_q = L_q·i_q

其中，ω_e为电角速度，与机械转速的关系为ω_e = n_p·ω_m（n_p为极对数）。在FOC控制中，通常令i_d=0（表贴式电机）或采用最大转矩电流比（MTPA）控制（内置式电机），以实现单位电流下的最大转矩输出。值得注意的是，Park变换需要准确的转子位置角θ_e，而在无感FOC中，这个角度正是需要估算的量。因此，SMO估算的θ_e质量直接决定了Park变换的准确性，进而影响整个控制系统的性能。

2.5 坐标变换在无感FOC中的工程应用与案例拆解

在实际的无感FOC系统中，坐标变换的流程是：首先通过采样电路获取三相电流i_A、i_B、i_C（通常只需采样两相，第三相通过基尔霍夫定律计算），然后依次进行Clark变换和Park变换，得到d-q轴电流i_d、i_q，送入电流环PI控制器。同时，电压指令u_d、u_q经过反Park变换得到u_α、u_β，再通过SVPWM调制生成PWM波驱动逆变器。

下面以一个典型的表贴式PMSM无感FOC驱动系统为例，说明坐标变换的具体应用。假设电机参数为：额定功率400W，极对数n_p=4，定子电阻R_s=2.5Ω，定子电感L_s=6.5mH，永磁体磁链ψ_f=0.085Wb。在某一运行工况下，采样得到三相电流：i_A=2.0A，i_B=-1.2A，i_C=-0.8A（三者和为零）。

第一步：Clark变换（采用等幅值变换，系数2/3）

i_α = (2/3)·(2.0 - 0.5·(-1.2) - 0.5·(-0.8)) = (2/3)·(2.0 + 0.6 + 0.4) = 2.0A
i_β = (2/3)·(0 + (√3/2)·(-1.2) - (√3/2)·(-0.8)) = (2/3)·(-1.039 + 0.693) = -0.231A

第二步：假设SMO估算的转子位置为θ_e=30°，进行Park变换：

i_d = 2.0·cos(30°) + (-0.231)·sin(30°) = 1.732 - 0.1155 = 1.616A
i_q = -2.0·sin(30°) + (-0.231)·cos(30°) = -1.0 - 0.2 = -1.2A

从计算结果可以看出，如果θ_e估算准确，i_d和i_q将呈现直流量特性。若θ_e存在偏差，则i_d和i_q中会混入交流分量，导致电流环控制性能下降。因此，SMO的估算精度至关重要。

在实际调试中，我们曾遇到一个典型问题：电机在低速运行时（<200rpm），SMO估算的角度存在约15°的滞后，导致i_d无法稳定在0A附近，电机转矩波动明显。通过分析发现，问题根源在于低速时反电动势幅值较小（e_α = -ψ_f·ω_e·sin(θ_e)，ω_e较小时e_α很小），信噪比降低，滑模观测器的抖振现象加剧。解决方案是：在SMO中引入自适应增益调节，即根据电角速度ω_e动态调整滑模增益k，同时在反电动势提取后增加低通滤波器（截止频率设为500Hz），并采用锁相环（PLL）替代反正切函数来提取角度，从而有效抑制了低速抖动。

2.6 常见误区与设计要点总结

下表总结了三种坐标系下的关键方程与适用场景，便于对照理解：

坐标系	电压方程形式	磁链方程	适用场景
三相静止（A-B-C）	u_ABC = R_s·i_ABC + dψ_ABC/dt	ψ_ABC = L_ABC·i_ABC + ψ_f·F(θ_e)	理论推导，不用于控制
两相静止（α-β）	u_αβ = R_s·i_αβ + L_s·di_αβ/dt + e_αβ	ψ_α = L_s·i_α + ψ_f·cos(θ_e)	SMO设计、反电动势观测
两相旋转（d-q）	u_dq = R_s·i_dq + L_dq·di_dq/dt + ω_e·J·ψ_dq	ψ_d = L_d·i_d + ψ_f，ψ_q = L_q·i_q	电流环PI控制、转矩控制

在无感FOC的工程实现中，坐标变换的数值精度直接影响控制性能。建议在嵌入式代码中使用浮点运算或Q格式定点运算，并注意角度θ_e的周期性处理（0~2π或-π~π）。此外，电流采样与坐标变换的时序必须严格同步，通常在PWM载波周期的中心点进行采样，以确保电流值的准确性。