4. 坐标变换与Clarke/Park变换
永磁同步电机（PMSM）无传感器FOC控制的核心，在于将电机定子三相绕组中耦合的交流变量，通过数学变换解耦为相互正交的直流分量，从而能够像控制直流电机一样独立控制转矩和磁通。这一过程依赖两大数学工具：Clarke变换（三相静止abc → 两相静止αβ）与Park变换（两相静止αβ → 两相旋转dq）。本章将从物理概念出发，逐步推导变换矩阵的由来，并给出Simulink仿真模型搭建思路与嵌入式C代码实现，帮助读者彻底理解“坐标变换”这一FOC基石。
上图展示了从三相物理绕组电流到两相旋转坐标系直流分量的完整映射链路。Clarke变换将三相耦合系统降维为两相正交系统，而Park变换则通过转子位置角θe将静止坐标系“旋转”到与转子同步的坐标系中，使得基波交流量变为直流量。这一变换不仅是数学上的简化，更是FOC能够实现高性能转矩控制的前提。
4.1 核心原理：从三相耦合到两相解耦的数学本质
三相PMSM定子绕组在空间上对称分布，通入对称三相正弦电流后，合成旋转磁动势。若直接对三相电流进行PI控制，由于各相电流均为交流量，稳态误差无法消除。坐标变换的本质，是寻找一个合适的线性变换矩阵，将三相静止坐标系下的时变交流量映射到与转子同步旋转的坐标系中，从而将交流控制问题转化为直流控制问题。
Clarke变换（3s/2s）：将三相静止abc坐标系变换到两相静止αβ坐标系。其物理意义是：在功率不变或幅值不变的约束下，用一个正交的两相系统等效替代三相系统。以幅值不变为约束，变换矩阵为：
| iα | | 1 -1/2 -1/2 | | ia |
| iβ | = | 0 √3/2 -√3/2 | | ib |
(零序分量忽略，星形无中线) | ic |该矩阵将三相电流投影到αβ轴上，其中α轴与a相轴线重合，β轴超前α轴90°。对于无中线星形连接，ia+ib+ic=0，因此只需采样两相电流即可重构第三相。
Park变换（2s/2r）：将两相静止αβ坐标系变换到两相旋转dq坐标系。变换矩阵为：
| id | | cosθe sinθe | | iα |
| iq | = | -sinθe cosθe | | iβ |其中θe为转子电角度（由编码器或观测器提供）。经过Park变换后，id为励磁电流分量（对应转子磁链方向），iq为转矩电流分量（垂直磁链方向）。在表贴式PMSM中，通常控制id=0，此时电磁转矩Te = 1.5·p·ψf·iq，与iq呈线性关系。
4.2 判断方法与操作流程：如何选择变换系数与实现步骤
在实际工程中，变换系数的选择直接影响后续PI参数整定和标幺化处理。常见的约束有两种：等幅值变换与等功率变换。等幅值变换下，变换前后电流矢量的幅值不变，便于直观理解电流极限；等功率变换下，变换前后功率表达式不变，但幅值变为原来的√(2/3)倍。FOC控制中，绝大多数嵌入式实现采用等幅值变换，因为电流环PI输出直接对应相电压幅值，便于SVPWM过调制处理。
操作流程如下：
- 电流采样：通过采样电阻或霍尔传感器获取至少两相电流（如ia、ib），计算ic = -ia - ib。
- Clarke变换：将ia、ib、ic代入矩阵，得到iα、iβ。注意在数字实现中，需将浮点运算转换为Q格式定点运算以提升效率。
- 获取转子位置θe：无传感器方案中，θe由滑模观测器（SMO）或扩展卡尔曼滤波器（EKF）估算；有传感器方案中，由编码器读取。
- Park变换：利用θe的sin/cos值，将iα、iβ变换为id、iq。
- PI调节：对id、iq分别进行PI控制，输出vd、vq。
- 逆Park变换：将vd、vq变换回vα、vβ。
- SVPWM调制：由vα、vβ合成三相占空比，驱动逆变器。
| 步骤 | 输入 | 输出 | 关键参数 | 常见错误 |
|---|---|---|---|---|
| 1. 采样 | ADC值 | ia, ib | 采样时刻需避开PWM开关噪声 | 采样点不在中心对齐处 |
| 2. Clarke | ia, ib, ic | iα, iβ | 系数2/3或√(2/3) | 系数选择不一致导致幅值偏差 |
| 3. 位置估算 | 反电动势/电流 | θe | 观测器带宽与收敛速度 | θe滞后导致dq轴耦合 |
| 4. Park | iα, iβ, θe | id, iq | sin/cos查表或CORDIC | 角度单位未统一（弧度/度） |
| 5. 逆变换 | vd, vq, θe | vα, vβ | 与Park矩阵互为转置 | 矩阵求逆错误 |
4.3 案例拆解：Simulink仿真模型搭建与C代码实现
下面以一个典型的PMSM FOC仿真模型为例，演示坐标变换模块的搭建与验证。仿真环境为MATLAB/Simulink 2023a，电机参数：Rs=2.875Ω，Ld=Lq=8.5mH，ψf=0.175Wb，极对数p=4。
Simulink实现步骤：
- 从Simscape Electrical库中拖入PMSM模块，设置电机参数。
- 使用“Three-Phase V-I Measurement”模块测量三相电流。
- 编写MATLAB Function模块实现Clarke变换：
function [i_alpha, i_beta] = clarke(ia, ib, ic)
i_alpha = (2/3)*(ia - 0.5*ib - 0.5*ic);
i_beta = (2/3)*((sqrt(3)/2)*ib - (sqrt(3)/2)*ic);
end- 编写Park变换模块，输入为iα、iβ和电角度θe：
function [id, iq] = park(i_alpha, i_beta, theta)
id = i_alpha*cos(theta) + i_beta*sin(theta);
iq = -i_alpha*sin(theta) + i_beta*cos(theta);
end嵌入式C代码实现（基于ARM Cortex-M4，Q15格式）：
// 等幅值Clarke变换，输入Q15格式电流
void Clarke_Q15(int16_t ia, int16_t ib, int16_t ic, int16_t *i_alpha, int16_t *i_beta) {
int32_t tmp;
tmp = (int32_t)ia - ((int32_t)ib >> 1) - ((int32_t)ic >> 1); // 2/3系数后置处理
*i_alpha = (int16_t)((tmp * 21845) >> 15); // 21845 ≈ 2/3 * 32768
tmp = ((int32_t)ib - (int32_t)ic) * 28378; // √3/3 ≈ 0.57735 * 32768 ≈ 18918，此处用28378为√3/2*32768
*i_beta = (int16_t)((tmp * 21845) >> 15);
}
// Park变换，查表法获取sin/cos
void Park_Q15(int16_t i_alpha, int16_t i_beta, int16_t sin_theta, int16_t cos_theta, int16_t *id, int16_t *iq) {
*id = (int16_t)(((int32_t)i_alpha * cos_theta + (int32_t)i_beta * sin_theta) >> 15);
*iq = (int16_t)(((int32_t)i_beta * cos_theta - (int32_t)i_alpha * sin_theta) >> 15);
}在仿真中，给定电机转速1000rpm，负载转矩2N·m，观察id、iq波形。正确实现时，id应稳定在0附近，iq为恒定值约3.2A（对应转矩2N·m）。若id不为0，说明Park变换角度θe存在偏差或PI参数未整定好。
4.4 常见误区与总结复盘
误区一：混淆等幅值与等功率变换系数。若在等幅值变换下使用等功率的逆变换矩阵，会导致电压指令幅值偏差约15.5%，轻则电流环振荡，重则过流保护。务必在代码注释中明确标注变换约束。
误区二：忽略零序分量。对于三相四线制系统或非理想星形连接，零序电流通路存在，Clarke变换需保留零序分量。但在电机驱动中，通常为三相三线制，零序电流为零，可安全忽略。
误区三：Park变换角度未做归一化处理。θe通常由增量式编码器累加得到，若角度溢出未处理，会导致sin/cos计算跳变，引起电流尖峰。建议将θe归一化到[0, 2π)区间，或使用模运算处理。
误区四：采样时刻与PWM中心不对齐。电流采样应在PWM载波谷底或峰顶进行，此时开关噪声最小。若采样点偏移，会引入高频噪声，导致Clarke变换后的iα、iβ含有大量纹波。
总结复盘：坐标变换是FOC的“数学桥梁”，其正确性直接决定控制性能。建议读者在仿真中先验证开环变换（给定理想正弦电流，观察dq波形是否为直线），再接入闭环。在代码实现中，务必统一Q格式、明确变换系数、处理角度溢出，并通过示波器或DA输出观察id/iq波形，确保变换无误后再进行后续PI整定。