1. 电机控制里为什么要反复折腾坐标变换

做电机控制（尤其是永磁同步电机PMSM和异步电机FOC矢量控制）的工程师，对克拉克变换这个名字肯定不陌生。但说实话，我见过不少刚入行的同事，公式会背、代码能抄，问一句"为什么非得从三相ABC变到两相αβ"，就答不上来了。这篇文章我尽量把克拉克变换的来龙去脉讲透，包括数学推导、工程实现、以及那些教材里一般不写的坑。

先说清楚它解决的核心问题。电机里有三个绕组，通入三相正弦电流，每一相电流单独看都在随时间变化，三相加起来才能产生一个旋转的磁场。问题在于，三相系统里有三个变量（ia、ib、ic），它们之间还有约束关系（ia+ib+ic=0），控制器要同时盯着三个变量去调节，数学上麻烦，物理上也不直观。

克拉克变换做的事，就是把三相静止坐标系（ABC）投影到一个两相静止坐标系（αβ）上。α轴和A相轴对齐，β轴逆时针超前α轴90度。变换之后，原来三个耦合的交流量，变成两个正交的交流量。虽然电压电流仍然是随时间变化的正弦波，但至少从三个变量减到了两个，而且这两个变量在空间上正交。这一步做完，后续再接上Park变换（旋转变换），就能把交流量变成直流量，从而用经典的PI控制器做无静差调节。可以说，克拉克变换是整个FOC链路的第一块基石。

这篇文章适合谁看？如果你正在做电机驱动、逆变器控制、并网逆变器或者电能质量相关的项目，想彻底搞明白坐标变换的数学原理，还想知道实际写代码时有哪些细节容易出错，那么这篇内容应该能对你有帮助。基础部分我会从线性代数的角度推一遍，代码部分也会给可直接移植的C语言例程。

2. 等幅值与等功率：两种克拉克变换矩阵的来源

克拉克变换的公式，不同教材里长得不太一样，根源在于系数选择不同。主流有两种：等幅值变换（系数2/3）和等功率变换（系数√(2/3)）。很多初学者会在这里被绕晕，我建议先理解变换的几何意义，再记住系数差异。

2.1 从ABC到αβ的几何投影

三相绕组的轴线在空间上互差120度。假设A相轴线在0度，B相在120度，C相在240度。任何一个时刻，三相绕组产生的磁动势，可以看成三个空间矢量（分别沿A、B、C方向，模长对应三相电流）的矢量和。数学上，把这三个矢量投影到α轴和β轴上：

fα = fa·cos0° + fb·cos120° + fc·cos240° = fa - (1/2)fb - (1/2)fc

fβ = fa·sin0° + fb·sin120° + fc·sin240° = (√3/2)fb - (√3/2)fc

这是最原始的投影关系。注意，这个矩阵的行列式不是1（事实上计算一下，如果系数统一乘上k，矩阵的行列式是 k²·(3√3/2)），所以这个变换并不保持矢量的模长不变，它本质上是把一个三维向量投影到二维平面上。这里还没有确定前面的比例系数k，而k，正是两种变换流派的分歧点。

2.2 系数2/3是怎么来的：等幅值约束

等幅值变换要求变换前后合成矢量的幅值相等。怎么理解？假设三相电流是标准的正弦平衡系统：

ia = I·cos(ωt)

ib = I·cos(ωt - 120°)

ic = I·cos(ωt + 120°)

代入投影式（先不乘系数k）：

fα = I·cos(ωt) - (1/2)·I·cos(ωt-120°) - (1/2)·I·cos(ωt+120°)

用三角恒等式展开：cos(ωt-120°) + cos(ωt+120°) = 2·cos(ωt)·cos(120°) = -cos(ωt)

所以 fα = I·cos(ωt) - (1/2)·(-cos(ωt)) = (3/2)·I·cos(ωt)

同理，fβ = (√3/2)·I·cos(ωt-120°) - (√3/2)·I·cos(ωt+120°) = (3/2)·I·sin(ωt)

也就是说，直接投影的结果是（3/2）倍的原始幅值。要让它回到原幅值I，就得乘上2/3。这就是等幅值克拉克变换的由来：

| fα | | 1 -1/2 -1/2 | | fa | | | = 2/3 | | | | | fβ | | 0 √3/2 -√3/2 | | fb | | fc |

这里有个重要的工程含义：采用等幅值变换后，变换后的电压电流幅值，和你在示波器上看到的相电压、相电流幅值一致。做PI控制器参数调节时，可以直接用电流环给定值和采样值做比较，不用额外换算。这也是大多数电机控制项目（包括ST的MCSDK、TI的InstaSPIN）默认使用等幅值变换的原因。

2.3 系数√(2/3)是怎么来的：等功率约束

等功率变换的出发点不一样。如果想保持变换前后的功率不变，需要让变换矩阵保持正交性，即变换矩阵的逆等于其转置。这就要求变换矩阵的每一行是单位向量。算一下投影式每一行的模长：第一行 [1, -1/2, -1/2] 的模长是 √(1+1/4+1/4) = √(3/2)。要让它是单位长度，就得整体乘上 √(2/3)。同理可验证第二行。

| fα | | 1 -1/2 -1/2 | | fa | | | = √(2/3) | | | | | fβ | | 0 √3/2 -√3/2 | | fb | | fc |

等功率变换的好处是，变换前后功率表达式一致，做完变换后计算有功无功的时候不需要额外乘系数。这个性质在电网侧变流器（比如并网逆变器的功率计算）里很常用，因为直接关系到电能计量和功率因数校正。

两种变换本身没有谁对谁错。选择等幅值，PI调节器调起来直观；选择等功率，功率计算方便。我在实际项目中两种都用过，简单总结一下：

说起来还有个容易混淆的点。等幅值变换里，逆变换矩阵的系数不是2/3，而是1。因为正变换乘了2/3，逆变换要抵消回来。而等功率变换，正变换乘√(2/3)，逆变换由于正交性，可以直接用转置矩阵，也就不用再额外乘系数。有的教材逆变换写的是行矩阵，有的写的是列矩阵，本质上都是同一套东西，容易把人绕晕。我习惯在代码注释里写清楚用的是哪种系数，以及在变换和逆变换之间保持对应关系，否则很容易出现幅值对不上的问题。

2.4 第三行零序分量的真相

细心的读者会注意到，上面的变换矩阵只有两行，但输入是三相，这在数学上其实不完整。一个完备的三维到三维变换，还应该包含第三轴——零序轴。完整的克拉克变换会多一行：

f0 = k·(fa + fb + fc)

对于对称三相系统，fa + fb + fc = 0，所以零序分量为零，变换矩阵里的第三行不起作用。这也是为什么大多数应用里只取两行就可以的原因。

但我强烈建议在写代码的时候还是算一下零序分量，因为它的值能帮你诊断系统故障。比如在三相四线制系统里，如果三相电流不平衡，零序电流不为零，通过零序分量的幅值大小可以判断不平衡程度。在电机控制里，如果电流传感器采样有偏差或者某相断线，零序分量也会异常。我做过一个电机堵转检测方案，就是在克拉克变换的同时算零序电流，一旦发现零序分量超过阈值就报故障，比单纯的过流判断灵敏得多。

顺带说一句，在有些文献里，等功率变换的第三行系数是 1/√2；等幅值变换的第三行系数是 1/3（有的写1/2，取决于是否归一化），这个细节各家定义不一，参考时要注意上下文。

3. 代码实现之前先搞清楚的几个数学细节

明白了变换矩阵怎么来的，接下来要落地。但在写代码之前，有几个数学细节如果不注意，后面查bug会查到怀疑人生。

3.1 角度基准：α轴对齐A相是默认约定

不知道你有没有遇到过这样的场景——两个人联调，一个人用的是α轴与A相轴重合的坐标系，另一个人用的是α轴与AB线电压方向重合的坐标系，结果两边的d轴角度差了30度，电机怎么都调不起来。这其实不是算法问题，是坐标系基准不统一。

克拉克变换最基本的定义是α轴与A相绕组轴线重合。如果你的硬件上电角度传感器（编码器/霍尔）的零点不是对齐到A相轴，那就要在软件里做角度偏置补偿。这个偏置通常用“对齐上电”或者“角度自动辨识”的方式获取。

另一个相关的坑是：采样电阻上的电流方向定义。电流从逆变器流入电机绕组定义为正，还是从电机流入逆变器定义为正？不同硬件板子的定义可能不同。方向定义反了，变换出来的θ会偏180度，电机会反着转或者剧烈抖动。所以在做克拉克变换之前，先确认采样电流的正方向，在代码里统一处理。

3.2 三相电流采样与ia+ib+ic=0的使用

平衡三相系统满足基尔霍夫电流定律：ia + ib + ic = 0。所以实际采样时可以只采两相（比如ia和ib），第三相用公式推出来：ic = -(ia + ib)。这样做可以节省一个电流传感器，降低硬件成本。

但要注意，这个关系只在无零序通路时成立。电机三相绕组通常是Y型连接且中性点不引出，所以满足条件。如果是三相四线制或者电机中性点引出并接了零线，就不能这样省了。还有一个工程细节：如果采样噪声比较大，直接算出的ic可能和真实值有偏差。因为ia和ib各自有微小误差，用两路采样值推算第三路，误差不会累计出大问题，但如果ia和ib的噪声相位不一致，推算出来的ic噪声反而会被放大。我见过有人还是坚持采三路再平均，效果也未必更好。我的建议是：如果开环测试时电流波形毛刺很多，可以试试把推算的那一路改成滤波后再参与后续计算，或者干脆用三路采样取平均。这属于调优的经验范畴。

3.3 量纲问题：标幺值还是国际单位值

在DSP或者单片机里实现克拉克变换时，电流的单位可能是安培（A），也可能是标幺值（p.u.）。用标幺值的好处是所有量都归一到0~1（或±1附近），PI参数不会因为电流基数不同而变化，移植性强。TI或ST的库都是基于标幺值体系设计的。

如果直接用安培值算，需要注意数值范围。比如额定电流100A的驱动器，ADC采样出来的是电压（经霍尔或采样电阻转换成电压），先要标定成安培，再做变换。这里的标定系数如果出现量级错误，后果是电流环的PI参数怎么调都调不稳。

我在实际工作上更喜欢先归一化：定义Ibase = 额定电流峰值，id/d = i_adc / Ibase。这样克拉克变换前后的系数不随工况变化，后面Park变换、PI限幅、过流保护阈值都是基于标幺值去设计，系统里的所有环节都“眼见为实”。不过，标幺值有一个不那么直观的地方，等幅值变换下的标幺值电流，变换后幅值就是实际电流除以Ibase，如果Ibase取的是额定电流的有效值（RMS）而不是峰值，那尖峰工况下数值可能超过1，需要留意限幅逻辑。

4. 从公式到C代码的完整落地过程

理论清楚了，接下来看看怎么在嵌入式平台上写。下面给一个可直接用的角度变换模块，包含克拉克变换、Park变换和反变换，全部基于浮点实现。虽然很多量产代码会用定点或者查表法优化，但浮点版本逻辑最清晰，适合先跑通功能再优化。

4.1 克拉克变换的浮点实现

typedef struct { float alpha; float beta; float z0; // 零序分量，诊断用 } ClarkeOutput_t; typedef struct { float a; float b; float c; } AbcInput_t; // 等幅值克拉克变换 (2/3系数) ClarkeOutput_t Clarke_Transform(AbcInput_t* in) { ClarkeOutput_t out; float k = 2.0f / 3.0f; out.alpha = k * (in->a - 0.5f * in->b - 0.5f * in->c); out.beta = k * (0.8660254f * in->b - 0.8660254f * in->c); out.z0 = (in->a + in->b + in->c) / 3.0f; // 约等于零 return out; }

注意0.8660254就是√3/2的近似值。这里没有做输入校验，实际项目中如果ADC异常，NaN传进来会导致整个控制环失效，建议在调用这个函数之前对原始采样值做范围检查。

如果只采两相电流，可以少一个输入参数：

// 只采ia, ib，ic由基尔霍夫定律推算 ClarkeOutput_t Clarke_Transform_2Phase(float ia, float ib) { ClarkeOutput_t out; float ic = -(ia + ib); out.alpha = ia; // 代入原始公式化简后得到 out.beta = 0.5773503f * (ia + 2.0f * ib); // 1/√3 ≈ 0.5773503 out.z0 = 0.0f; return out; }

如果你推导一下就会发现，当把ic=-(ia+ib)代入标准克拉克变换公式后，α和β的表达式会大幅化简。但这里有个常见误解：化简后的公式只有在三相平衡（即ic严格等于-(ia+ib)）时才成立，实际系统中由于电机绕组不对称或采样时刻不同，直接用两路采样值推算第三路会引入等效误差。我建议如果控制精度要求高，还是采三路；如果硬件只支持两路采样，那上面这个化简式是标准做法，不用犹豫。

4.2 Park变换衔接：从αβ到dq

克拉克变换做完后，接着做Park变换，才能得到可用于PI控制的直流量。Park变换的输入是αβ坐标系下的量，以及电角度θ（转子位置角度，单位是电气角度）：

typedef struct { float d; float q; } ParkOutput_t; ParkOutput_t Park_Transform(float alpha, float beta, float theta_sin, float theta_cos) { ParkOutput_t out; out.d = theta_cos * alpha + theta_sin * beta; out.q = -theta_sin * alpha + theta_cos * beta; return out; }

为什么这里要传入sin和cos而不是直接传角度？因为在嵌入式实时控制里，算一次三角函数的开销相对较大（在Cortex-M4F上sin/cos约几十个周期，但如果有硬件FPU和三角函数加速指令会好很多），而且在中断里调用函数库函数容易出现不确定的延迟。更常见的做法是电角度计算模块同时输出θ的正余弦值（基于查表或CORDIC），直接传进来用，省去一次三角函数调用。

Park变换的系数这里是1，没有额外系数。也就是说等幅值克拉克变换出来的αβ量，直接做旋转变换就可以得到dq量，幅值保持。这也是FOC里默认采用等幅值框架的原因——电流环的给定值、反馈值、PI输出都在同一个幅值基准下，不需要乘来乘去。

4.3 逆变换：把dq电压拉回三相

电流环PI输出的Ud和Uq是在旋转坐标系下的直流量，要通过逆Park变换和逆克拉克变换变成三相占空比：

// 逆Park变换：dq -> αβ void InvPark_Transform(float ud, float uq, float theta_sin, float theta_cos, float* alpha, float* beta) { *alpha = theta_cos * ud - theta_sin * uq; *beta = theta_sin * ud + theta_cos * uq; } // 逆克拉克变换：αβ -> 三相 void InvClarke_Transform(float alpha, float beta, float* ua, float* ub, float* uc) { *ua = alpha; *ub = -0.5f * alpha + 0.8660254f * beta; *uc = -0.5f * alpha - 0.8660254f * beta; }

注意这里逆克拉克变换的系数判断：等幅值变换下，正变换乘2/3，逆变换不需要额外系数，因为αβ的幅值已经和三相幅值一致了。我见过有人在这里机械地套用"逆矩阵是转置乘以2/3"，算出来的占空比缩小了，电机带载能力明显不足。这其实还是对变换的物理含义理解不透。等幅值变换矩阵不是正交矩阵，它的逆矩阵不是简单的转置，一定要实际算一遍或者用验证过的代码。

逆变换的输入输出是电压，但FOC最终需要的是三相占空比。中间还要把ua、ub、uc换算成相对于直流母线中点的相电压，再除以母线电压的一半得到占空比，最后还要做SVPWM或者载波比较。有的集成库（比如ST的MC SDK）会把SVPWM和逆克拉克变换合在一起算，这时候你要搞清楚它的输出到底是"参考电压"还是"最终占空比"，不然容易在限幅环节出问题。

5. 工程实战中的高频坑与排查经验

理论、公式、代码都过了一遍，最后分享一些实际项目中容易被忽视的坑。这些坑往往不是算法本身的问题，而是算法与硬件、与系统其他部分交互时产生的问题。

5.1 母线电压波动对变换后电压的影响

逆克拉克变换输出的电压量，最终要调制到逆变器上。如果直流母线电压不稳，而你在算占空比时用的是额定母线电压而非实时采样值，那么实际加到电机上的电压就和控制算法期望的不一致。这在低压大电流系统里特别明显——电池电压随负载波动很大，瞬时压降可能超过10%。

我的做法是：在做逆变换的时候，直接把实时母线电压代入进行归一化。即先算出需要的相电压幅值，再除以(母线电压/2)得到占空比。这样母线波动不会影响电压环的增益，动态响应更均匀。另外，如果母线电压做了滤波，滤波器的时间常数不能太大，否则在突加负载的瞬间，占空比会瞬态饱和，表现为电流过冲。滤波时间常数我一般取几十到几百微秒的量级，具体看开关频率和控制周期。

5.2 采样时刻偏移导致的角度误差

FOC的电流采样有个老生常谈但永远有人踩的坑：必须在PWM的载波中心处采样。原因是逆变器开关动作会引起电流纹波，在载波顶点（或底点）附近，三相电流纹波处于对称位置，此时采样能得到接近于基波的平均值。如果在开关瞬间附近采样，采到的可能带有大的开关噪声，变换出来的αβ电流会有明显的毛刺。

但还有一个更隐蔽的问题：从采样时刻到真正执行克拉克变换的时刻之间存在延迟。如果在采完电流之后，读到的电角度不是采样瞬间的角度，而是当前时刻（比如晚了半个控制周期）的角度，那Park变换用的角度就和采样电流不匹配。这个角度误差在高速电机上会被放大，造成dq轴耦合、电流环振荡。

解决方法是记录采样时刻的角度值。比如使用ADC注入式采样（由PWM触发ADC，DMA将数据搬到内存），此时用定时器的计数器值去捕获电角度位置，确保αβ变换和角度用的是同一个时刻的"切片"。或者在没有硬件支持的情况下，用"上一次中断的角度+速度×时间"来外推采样时刻的角度。无论哪种方案，核心思想是：电流和角度必须严格同帧。

5.3 代码里的小数精度与定点化陷阱

浮点MCU上用float算克拉克变换绰绰有余。但如果你的平台是低成本的定点MCU（比如没有FPU的Cortex-M0+），就要考虑定点化。定点化的第一原则不是疯狂移位，而是先确定Q格式。比如电流采样值范围是±100A，ADC分辨率12位，那么电流变量取Q12格式可以获得比较好的精度，但要注意中途乘法运算的中间量可能会溢出，需要适时缩放到Q24或更高精度。

我只提醒几个容易出问题的地方：

• 0.8660254在定点里如果用整数近似，误差稍微大一点，变换结果在小电流时会失真。
• 乘加的累加顺序，先乘后加和先加后乘的舍入误差不同，工程上差异不太大，但一致性测试时最好固定。
• 如果用了volatile关键字或者某种编译优化选项，浮点运算的结果可能在不同优化级别下有微小差异。调试和发布版本用同一套优化选项，否则会遇到"调试一切正常，release后电流波形变了"这种玄学问题。

5.4 故障诊断：利用克拉克变换做健康监测

前面提到零序分量可以用来诊断系统故障，这里具体展开一下。

克拉克变换把三相量分解成正交的两相量和一个零序分量。在电机驱动系统中：

• 如果零序电流分量持续增大，多半是电流传感器零点漂移、三相绕组不对称或者某一相接触不良。
• 如果αβ量的轨迹（在αβ平面上画李萨育图）变成椭圆而不是圆，说明某一相增益异常，可能是采样电阻老化或者运放偏置问题。
• 如果αβ量的幅值突然跳动，可能是采样回路受到电磁干扰。

我在量产项目里就加过这样一个诊断逻辑：统计一个工频周期内零序电流的均方根，当它超过某阈值（比如额定电流的5%）且持续100ms，就报"电机电流不平衡"故障。这个方法简单可靠，能避免很多因为接线松动导致的售后问题。

还有一个技巧：用克拉克变换后的α、β电流算瞬时功率。在电机的无传感器算法里，观测器（比如滑模观测器）的输入就是αβ坐标系下的电压电流，这也是克拉克变换的一个重要延伸场景。如果你后面准备研究无感FOC，先把克拉克变换吃得透透的是非常必要的。

6. 关于系数选择与扩展思考

最后再聊一点我对两种变换系数的选择偏好，以及可能遇到的跨团队协作问题。

前面说过，等幅值和等功率没有绝对优劣。但工程项目的技术选型往往是"少数服从多数"且"跟随成熟方案"：电机控制生态（ST、TI、Infineon的参考设计）基本都用等幅值；电力系统/电网侧的功率计算更倾向等功率。如果你的系统里既有电机驱动模块和功率计算模块，比如车载主驱加DC-DC、或者是微网里的PCS储能变流器，那么建议全系统统一用同一种变换约定，避免接口处还要做系数换算。

结合我个人的项目经验，有两点特别想强调：

第一，变换矩阵本身只是几行几列的数学，决定系统性能的往往不是变换本身，而是它和采样、角度、PWM调制之间的配合。很多电机调不动或者调不稳，最后定位出的根因都不是克拉克变换，而是在变换前后某一环节系数错了或者数据不同步。排查这类问题时，我的办法是"分段校验"：人为施加一组已知的三相电流（比如在电流闭环开环状态下给一个直流偏置），用示波器核对αβ输出的幅值和相位，再逐个环节往后查。这样基本半小时内能找到问题。

第二，如果之后接触无感FOC、直接转矩控制（DTC）或者多相电机（比如五相、六相电机），你会发现坐标变换的思想仍然适用，只不过维度和变换矩阵扩展了。多相电机的Vector Space Decomposition（VSD），本质上就是克拉克变换的理念在高维空间的推广。所以把这一课学明白，后面很多高级控制算法都容易上手。

克拉克变换这个看似简单的线性变换，一头的物理解释是"把三绕组产生的旋转磁场描述成两轴分量"，另一头连着PWM调制、观测器、故障诊断等工程问题。无论从哪个方向切入，把矩阵系数、逆变换、采样时序这些细节都梳理清楚，对做电机控制或者电力电子的人来说都值得。