Clark变换与Park变换公式差异详解：从等幅值到等功率的FOC调试指南

1. 从一次调车翻车说起：为什么同一个Clark变换，两本书能写出两个样子

刚入行做电机控制那会儿，我照着某本教材把Clark变换敲进代码，电流环跑起来波形挺漂亮。后来换了个项目，参考另一份资料，同样的永磁同步电机，同样的采样电路，电流波形直接歪成麻花，转矩脉动大得能把桌子震响。当时第一反应是硬件坏了，查了半天采样电阻、运放偏置、ADC时序，全都没问题。最后把两份公式摆在一起逐行对比，才发现问题出在一个系数上——一个用了等幅值变换，一个用了等功率变换，而我的电流环PID参数是按等幅值整定的，套到等功率的坐标下，增益相当于被悄悄改了一截。

这件事让我意识到，Clark变换和Park变换虽然形式简单，但不同教材、不同厂商代码、不同仿真模型里的写法差异，足以让一个能跑的FOC系统直接趴窝。更麻烦的是，这些差异往往不会写在显眼的地方，而是藏在系数、符号、坐标系定义这些细节里。你如果只是复制粘贴，不去理解背后的约定，迟早会在某个项目上踩坑。

这篇内容就是把我这些年在这两个变换上踩过的坑、理过的逻辑、总结的对照方法，完整地摊开讲一遍。不管你是刚接触FOC的新手，还是已经能跑通电流环但被公式差异搞晕的工程师，都能从中找到可以直接用的判断方法和实操建议。核心关键词Clark变换、Park变换、电机控制、FOC、MATLAB/Simulink会贯穿始终，但我不打算写成教科书，而是按实际调试的顺序，把每个差异点拆开，告诉你它从哪来、影响什么、怎么统一。

2. 变换的本质：从三相到两相，到底在转什么

2.1 三相电流为什么需要变换

永磁同步电机的定子上有三个绕组，空间上互差120度电角度。通入三相对称正弦电流后，合成的是一个旋转的磁场。问题是，如果你直接在ABC坐标系下做控制，三个电流互相耦合，而且都是交流量，PID控制器对交流量的跟踪效果很差，稳态误差和相位滞后都很明显。

Clark变换的作用，是把这三个互差120度的物理量，投影到一个正交的静止坐标系上，也就是通常说的α-β坐标系。变换之后，三相的耦合关系被解开了，但α轴和β轴上的量仍然是交流的，因为合成磁场还在转。Park变换接着上场，把静止的α-β坐标系旋转到跟着转子一起转的d-q坐标系上。在d-q坐标系里，电流变成了直流量，PID控制器终于能发挥正常水平，这就是FOC的核心逻辑。

整个链条可以简单记成：三相交流 → Clark → 两相静止交流 → Park → 两相旋转直流 → PID控制 → 反Park → 反Clark → 三相输出。每一步都有明确的物理意义，但每一步的公式写法都可能因为约定不同而出现差异。

2.2 等幅值变换和等功率变换的分水岭

Clark变换最常见的两种形式，一种是等幅值变换，一种是等功率变换。等幅值变换的目标是让变换后的α轴电流幅值等于三相电流的幅值，系数里会出现2/3。等功率变换的目标是让变换前后的功率不变，系数里会出现根号下2/3。

为什么会有这两种？因为应用场景不同。做电流环控制的时候，你关心的是电流的幅值能不能被PID准确跟踪，所以等幅值变换更直观，PID参数整定起来也简单。做功率分析或者效率计算的时候，你关心的是功率守恒，等功率变换更合适。但问题在于，很多教材不会明确告诉你它用的是哪一种，你只能从系数上反推。

我个人的习惯是，在电机控制代码里统一用等幅值变换，因为电流环的参考值和反馈值都在同一个尺度下，PID参数不需要额外换算。如果仿真模型里用了等功率变换，我会在电流环入口和出口做一次系数补偿，避免把两套尺度混在一起。

2.3 系数差异对PID参数的实际影响

假设你从等幅值变换切换到等功率变换，α轴电流的幅值会变成原来的根号下3/2倍，大约是1.225倍。如果你没意识到这个变化，PID的比例增益相当于被放大了1.225倍，积分增益同理。对于电流环来说，这可能导致超调增大、振荡加剧，严重的时候直接发散。

反过来，从等功率切到等幅值，增益被缩小，响应变慢，转矩跟踪出现滞后。这些现象在调车的时候很容易被误判为电机参数不准或者采样有问题，实际上只是变换系数没对齐。

我的做法是，在代码里把变换系数定义成宏或者常量，注释清楚用的是哪一种，然后在电流环整定之前，先用示波器或者仿真确认α轴电流的幅值和三相电流的幅值是否一致。如果一致，说明是等幅值；如果α轴幅值是三相幅值的1.225倍，那就是等功率。这个检查只需要做一次，但能省掉后面很多莫名其妙的调试时间。

3. 不同教材公式差异的根源：约定、符号与坐标系

3.1 坐标系定义：α轴对齐A相还是超前90度

Clark变换的坐标系定义，不同教材有不同的选择。最常见的是α轴与A相绕组轴线重合，β轴超前α轴90度。但也有教材把α轴定义为超前A相90度，或者把β轴定义为滞后α轴90度。这些定义差异会直接改变变换矩阵里的符号。

比如，当α轴对齐A相时，A相电流在α轴上的投影是完整的，B相和C相各贡献一个余弦分量。如果α轴超前A相90度，那么A相电流在α轴上的投影就变成了零，全部投影到β轴上。这两种定义下，变换矩阵完全不同，但物理结果应该是一致的，只是坐标系的摆放方式不同。

问题在于，如果你从一份资料里抄了变换矩阵，又从另一份资料里抄了反变换矩阵，而两者的坐标系定义不一致，那电流环的反馈路径就错了。我见过一个项目，正变换用的是α轴对齐A相，反变换用的是α轴超前A相，结果电机能转，但电流波形完全不对，转矩脉动大得离谱。后来把两份矩阵统一到同一个定义下，问题立刻消失。

3.2 符号约定：β轴的正方向到底朝哪

β轴的正方向定义也是一个容易踩坑的地方。有的教材规定β轴超前α轴90度，有的规定滞后90度。这个符号差异会影响Park变换里sin和cos的符号，进而影响d轴和q轴电流的符号。

在电机控制里，q轴电流通常对应转矩分量，d轴电流对应励磁分量。如果β轴方向反了，q轴电流的符号也会反，PID控制器会朝着错误的方向调节，电机要么不转，要么飞车。这种问题在仿真里可能不明显，因为仿真模型可以自己定义正方向，但到了实际硬件上，编码器的计数方向和电流采样方向必须和变换的符号约定一致，否则就是灾难。

我的经验是，在项目初期就把所有符号约定写在一张纸上：α轴对齐A相，β轴超前α轴90度，d轴对齐转子磁链方向，q轴超前d轴90度。然后所有代码、仿真、文档都按这个约定来，任何人修改都要同步更新这张纸。这个习惯看起来笨，但能避免很多跨人协作时的低级错误。

3.3 三相到两相的矩阵形式：2/3系数放在哪

Clark变换的矩阵形式，2/3系数可以放在矩阵外面，也可以分配到矩阵里面的每个元素上。这两种写法在数学上等价，但在代码实现时会影响中间变量的范围，进而影响定点数的精度。

比如，对于定点DSP，如果2/3系数放在外面，中间变量可能会超出Q格式的表示范围，导致溢出。如果分配到矩阵里，每个元素都带一个小数系数，中间变量的范围更可控。我在用定点DSP做FOC的时候，倾向于把系数分配到矩阵里，并且根据电流的额定值选择合适的Q格式，确保中间变量不会溢出。

另外，有些教材会把Clark变换写成三相到两相的形式，有些会写成两相到三相的反变换形式，还有些会写成三相到零序加两相的形式。这些形式在数学上可以互相推导，但在代码里直接套用的时候，一定要确认输入输出的维度是否匹配。我见过有人把三相到两相的矩阵直接当成两相到三相用，结果电流环的反馈值完全错误，电机根本转不起来。

4. 手把手推导：从三相到d-q的完整链条

4.1 等幅值Clark变换的推导过程

假设三相电流为ia、ib、ic，满足ia+ib+ic=0。定义α轴与A相重合，β轴超前α轴90度。根据投影关系，α轴电流等于ia减去ib和ic在α轴上的投影。由于ib和ic在α轴上的投影都是负的，且幅值为1/2，所以：

iα = ia - (1/2)ib - (1/2)ic

因为ia+ib+ic=0，所以ib+ic=-ia，代入上式得到：

iα = ia - (1/2)(-ia) = (3/2)ia

为了保持幅值不变，需要乘以2/3，所以：

iα = (2/3)(ia - (1/2)ib - (1/2)ic)

β轴电流等于ib和ic在β轴上的投影之差。ib在β轴上的投影是(根号3/2)ib，ic在β轴上的投影是-(根号3/2)ic，所以：

iβ = (2/3)((根号3/2)ib - (根号3/2)ic) = (1/根号3)(ib - ic)

这就是等幅值Clark变换的经典形式。你可以看到，2/3系数被分配到了两个方程里，iα的系数是2/3，iβ的系数是1/根号3。

4.2 等功率Clark变换的系数从哪来

等功率变换的目标是让变换前后的功率相等。在三相系统中，功率等于三相电压和电流的点积。在α-β坐标系中，功率等于α轴和β轴的电压电流点积。为了让两者相等，变换矩阵需要满足正交条件，即矩阵的转置乘以矩阵等于单位矩阵。

对于三相到两相的变换，满足正交条件的矩阵系数是根号下2/3，而不是2/3。所以等功率Clark变换的形式是：

iα = 根号下(2/3) * (ia - (1/2)ib - (1/2)ic) iβ = 根号下(2/3) * ((根号3/2)ib - (根号3/2)ic) = (1/根号2)(ib - ic)

你可以验证一下，这个矩阵的转置乘以矩阵确实等于单位矩阵。这就是等功率变换的数学来源。

4.3 Park变换的旋转矩阵与角度定义

Park变换是把α-β坐标系旋转到d-q坐标系，旋转角度是转子电角度θ。标准的旋转矩阵是：

id = iα * cosθ + iβ * sinθ iq = -iα * sinθ + iβ * cosθ

这个形式对应的是d轴对齐转子磁链方向，q轴超前d轴90度。如果你看到的公式里sin和cos的符号不同，或者id和iq的位置互换，那说明角度定义或者坐标系定义不同。

比如，有些教材把q轴对齐转子磁链方向，d轴超前q轴90度，那么旋转矩阵就会变成：

id = iα * sinθ - iβ * cosθ iq = iα * cosθ + iβ * sinθ

这两种定义在物理上等价，但代码里的角度来源必须匹配。如果你用的编码器角度是d轴对齐转子磁链的，但Park变换用的是q轴对齐的公式，那d轴和q轴电流就互换了，电流环完全失控。

4.4 反变换的推导与代码实现要点

反Park变换是把d-q坐标系旋转回α-β坐标系，旋转角度是-θ。标准形式是：

iα = id * cosθ - iq * sinθ iβ = id * sinθ + iq * cosθ

反Clark变换是把α-β坐标系还原到三相坐标系。对于等幅值变换，反变换是：

ia = iα ib = -(1/2)iα + (根号3/2)iβ ic = -(1/2)iα - (根号3/2)iβ

在代码实现时，要注意反变换的输出是否满足ia+ib+ic=0。如果因为浮点误差导致和不完全为零，可以强制让ic = -ia - ib，避免零序分量引入额外的谐波。

另外，反Clark变换的输出通常要经过SVPWM或者SPWM调制，生成六路PWM波。如果反变换的系数和正变换不匹配，比如正变换用等幅值，反变换用等功率，那输出电压的幅值就会偏差1.225倍，电机的实际电流和参考值对不上，转矩控制精度直接崩掉。

5. MATLAB/Simulink仿真验证：用波形说话

5.1 搭建一个最小验证模型

在Simulink里验证Clark和Park变换，不需要完整的电机模型，只需要三个正弦波源、一个Clark变换模块、一个Park变换模块、一个角度源和一个示波器。三相正弦波的幅值设为1，频率设为50Hz，相位互差120度。角度源设为斜坡信号，斜率对应电频率。

Clark变换模块可以用Simulink的Fcn模块或者MATLAB Function模块实现。我习惯用MATLAB Function，因为可以直接写代码，方便对照。等幅值变换的代码如下：

function [ialpha, ibeta] = clark_amp(ia, ib, ic) ialpha = (2/3) * (ia - 0.5*ib - 0.5*ic); ibeta = (2/3) * (sqrt(3)/2*ib - sqrt(3)/2*ic); end

等功率变换的代码：

function [ialpha, ibeta] = clark_pow(ia, ib, ic) ialpha = sqrt(2/3) * (ia - 0.5*ib - 0.5*ic); ibeta = sqrt(2/3) * (sqrt(3)/2*ib - sqrt(3)/2*ic); end

Park变换的代码：

function [id, iq] = park(ialpha, ibeta, theta) id = ialpha * cos(theta) + ibeta * sin(theta); iq = -ialpha * sin(theta) + ibeta * cos(theta); end

把这三个模块串起来，用示波器观察ialpha、ibeta、id、iq的波形。如果一切正常，ialpha和ibeta应该是幅值为1的正弦波，相位差90度；id和iq应该是直流量，id等于0，iq等于1（假设三相电流的幅值为1，且初始相位对齐）。

5.2 观察等幅值和等功率的波形差异

把Clark变换模块从等幅值切换到等功率，重新运行仿真。你会看到ialpha和ibeta的幅值从1变成了1.225，id和iq的幅值也相应变化。如果Park变换的角度源没有变，id和iq的直流值会变成1.225和0。这个差异在仿真里一目了然，但在实际代码里，如果你不刻意去检查，很容易忽略。

我建议在仿真模型里加一个开关，可以在等幅值和等功率之间切换，然后对比两种模式下电流环的响应。你会发现，如果PID参数是按等幅值整定的，切到等功率后，电流环的阶跃响应会出现明显的超调，调节时间变长。反过来，如果PID参数是按等功率整定的，切到等幅值后，响应变慢，跟踪误差增大。

5.3 用FFT分析变换后的谐波含量

Clark和Park变换本身是线性变换，不会引入新的谐波。但如果三相电流不平衡，或者采样有偏置，变换后的α-β轴和d-q轴会出现二倍频或者直流分量。在Simulink里可以用FFT工具分析id和iq的频谱，如果看到明显的二倍频，说明三相电流不平衡；如果看到直流偏置，说明采样有零漂。

我在实际项目中遇到过一个问题：电机低速运行时，id和iq上出现了明显的二倍频波动，转矩脉动很大。查了半天，发现是三相电流采样电阻的温漂导致三相增益不一致。后来在Clark变换之前加了一个增益校准环节，问题才解决。这个经验说明，变换公式本身没问题，但输入信号的质量直接决定变换后的效果。

5.4 从仿真到DSP代码的移植检查清单

仿真跑通之后，把代码移植到DSP上，需要检查几个关键点。第一，角度来源是否一致。仿真里的角度是理想的斜坡信号，实际DSP里的角度来自编码器或者观测器，可能存在偏移和噪声。第二，定点数的Q格式是否匹配。仿真里用的是浮点，DSP里可能是定点，系数的精度和中间变量的范围需要重新评估。第三，采样和计算的时序是否对齐。仿真里没有延迟，实际DSP里ADC采样、Clark变换、Park变换、PID计算、SVPWM更新之间存在多个时钟周期的延迟，如果不对齐，电流环的相位裕度会下降。

我的做法是，在DSP代码里加一个调试变量，把Clark变换和Park变换的中间结果通过DAC或者串口输出，和仿真波形逐点对比。如果发现幅值或者相位有偏差，就逐级排查，直到两者一致。这个过程看起来繁琐，但能确保算法在硬件上的行为符合预期。

6. 常见问题与排查技巧实录

6.1 电流环振荡：先查变换系数再调PID

电流环振荡是最常见的问题，很多人第一反应是调PID参数，但有时候问题根本不在PID上。如果你从一份资料里抄了Clark变换，从另一份资料里抄了Park变换，两者的系数约定不一致，电流环的反馈值就是错的，PID怎么调都稳不了。

排查方法很简单：给电流环一个阶跃参考，观察实际电流的响应。如果实际电流的幅值和参考值对不上，比如参考值是1A，实际值是1.225A，那说明变换系数不匹配。这时候先把变换系数统一，再重新整定PID。

6.2 电机反转或飞车：检查β轴和q轴的符号

电机反转或者飞车，通常是符号约定出了问题。β轴的正方向反了，q轴电流的符号就会反，PID控制器会朝着错误的方向调节。编码器的计数方向反了，角度θ的符号就会反，Park变换的旋转方向也会反。

排查方法：用手转动电机，观察编码器的计数方向是否和Park变换的角度增加方向一致。如果不一致，要么改编码器的计数方向，要么改Park变换的符号。另外，给q轴一个正的小电流，观察电机是否朝着预期的方向转。如果反转，就把q轴电流的符号取反，或者把β轴的定义反过来。

6.3 低速时电流波形畸变：零序和死区的影响

低速时电流波形畸变，通常和零序分量、死区效应有关。Clark变换假设ia+ib+ic=0，但实际采样中，由于ADC偏置和死区，三相电流之和可能不为零。这个零序分量会在α-β轴上产生额外的波动，导致电流波形畸变。

解决方法：在Clark变换之前，先计算三相电流的平均值，然后从每相电流中减去这个平均值，强制ia+ib+ic=0。另外，死区补偿也可以改善低速时的波形，但死区补偿的参数需要根据具体的驱动电路和PWM频率来整定。

6.4 不同厂商代码的变换矩阵对照表

这张表不是绝对的，但可以帮你快速判断一份代码或者一份资料用的是哪种约定。拿到一份新代码，先看Clark系数是2/3还是根号下2/3，再看β轴的符号，最后看Park变换里sin和cos的位置。三个信息一对照，基本就能确定它的坐标系定义。

6.5 定点DSP上的溢出与精度问题

定点DSP做Clark和Park变换，最容易遇到溢出和精度损失。比如，三相电流的采样值是12位ADC，范围是0到4095，对应电流范围是-10A到10A。如果直接用整数做Clark变换，中间变量可能会超出16位或者32位的表示范围。

我的做法是，先把ADC采样值转换成标幺值，比如用Q15格式表示-1到1之间的电流。然后Clark变换的系数也用Q15表示，乘法之后右移15位，保持Q15格式。Park变换的sin和cos也用Q15查表，乘法之后右移15位。这样整个链条都在Q15格式下运行，不会溢出，精度也够用。

需要注意的是，Q15格式的乘法结果需要用32位累加器，否则会丢失精度。另外，查表法的角度分辨率要足够高，否则Park变换的谐波会变大。我一般用1024点的正弦表，角度分辨率是0.35度，对于大多数电机控制应用足够了。

7. 统一约定的工程实践：从代码规范到团队协作

7.1 在代码里显式定义变换系数

我在所有电机控制项目里，都会在头文件里定义变换系数，比如：

#define CLARK_COEFF (2.0f/3.0f) #define SQRT3_2 (0.86602540378f) #define SQRT3 (1.73205080757f)

然后在Clark变换函数里直接引用这些宏，而不是硬编码数字。这样做的好处是，如果以后需要从等幅值切换到等功率，只需要改一个宏定义，所有相关代码自动更新。另外，注释里要写清楚当前用的是哪种变换，以及坐标系定义。

7.2 仿真和代码使用同一套角度定义

仿真里的角度源和代码里的角度来源必须一致。我习惯在仿真里用编码器的实际角度，而不是理想的斜坡信号。这样可以在仿真阶段就发现角度偏移、噪声和延迟对电流环的影响。如果仿真里用的是理想角度，代码里用的是实际角度，两者的行为可能完全不同。

另外，仿真里的Park变换和代码里的Park变换要用同一个函数或者同一套公式。我通常会把变换函数写成独立的MATLAB Function和C函数，然后逐行对比，确保逻辑一致。这个习惯能避免很多“仿真能跑，代码跑不了”的问题。

7.3 团队协作中的文档与评审要点

在团队协作中，Clark和Park变换的约定必须写进设计文档，并且在代码评审时重点检查。我见过一个项目，两个工程师分别负责电流环和速度环，电流环用的是等幅值变换，速度环用的是等功率变换，结果速度环的输出经过电流环时，幅值偏差了1.225倍，速度响应一直有稳态误差。后来在评审时发现了这个问题，统一了变换系数，问题才解决。

评审要点包括：Clark系数是2/3还是根号下2/3，β轴方向是超前还是滞后，Park变换的sin和cos位置，角度来源是编码器还是观测器，定点数的Q格式是否一致。这些点看起来琐碎，但每一个都可能成为项目翻车的原因。

7.4 从教材到产品的最后一公里

教材里的公式是理想化的，产品里的代码要考虑采样延迟、死区、温漂、定点精度、电磁干扰。Clark和Park变换作为FOC的核心环节，必须经过仿真验证、硬件在环测试、实际电机调试三个阶段的反复确认。我个人的经验是，在实验室里跑通的代码，到了现场可能因为温度变化导致采样偏置漂移，电流波形畸变，转矩脉动增大。这时候需要重新校准采样电路，或者在软件里加温度补偿。

最后一公里没有捷径，只能靠细致的调试和完整的记录。每次调试都要把变换系数、角度定义、PID参数、采样校准值记录下来，形成项目自己的知识库。下次遇到类似问题，可以直接查记录，不用从头再来。

8. 写在最后：几个让我少走弯路的习惯

第一个习惯，拿到任何一份电机控制代码，先找Clark和Park变换的实现，确认系数和符号约定。这个动作只需要几分钟，但能避免后面几天的无效调试。

第二个习惯，在Simulink里搭一个最小验证模型，把变换公式和实际代码逐行对比。仿真波形和实际波形对不上的时候，不要急着改PID，先查变换。

第三个习惯，把坐标系定义和符号约定写进项目文档，并且在代码注释里重复一遍。人脑记不住这么多细节，但文档和注释可以。

第四个习惯，定点DSP上的变换系数用宏定义，不要硬编码。等幅值和等功率的切换只需要改一个宏，不用翻遍整个工程。

这些习惯看起来不起眼，但在我做过的十几个电机控制项目里，每一个都帮我省下了至少一天的调试时间。Clark和Park变换的公式差异，说到底不是数学问题，而是工程约定问题。把约定统一了，公式自然就一致了。