3、坐标变换理论：Clark变换和Park变换的原理与实现

各位同学，今天我们来啃一块硬骨头——坐标变换。说实话，我刚入行那会儿，看到Clark和Park这两个名字，第一反应是“这又是哪个数学家在折腾我们”。但后来在风场调试时，我亲眼看着一台双馈风机因为坐标变换参数设错，直接报过流停机。那一刻我才明白：坐标变换不是纸上谈兵，它是矢量控制的命根子。

3.1 为什么非要从三相静止变到两相旋转？

你想想看，双馈电机的定子、转子都是三相绕组，通进去的是三相交流电。三相交流电有个麻烦——它随时间变化，而且三相之间还差120度。你要是直接拿三相电流去算转矩，那计算量能让你怀疑人生。

但如果我们换个思路：把三相静止坐标系下的交流量，变成两相旋转坐标系下的直流量。这样一来，控制问题就简化成了“调直流”的问题。PID控制器最喜欢直流信号了，稳得很。

我个人习惯把坐标变换比作“翻译官”：

• Clark变换：把三相静止（abc）翻译成两相静止（αβ）
• Park变换：把两相静止（αβ）翻译成两相旋转（dq）

这两步走完，你就能在dq坐标系下，像控制直流电机一样控制交流电机。说白了，这就是矢量控制的灵魂。

3.2 Clark变换：从abc到αβ

Clark变换的数学表达式其实很简单。假设三相电流为 i_a、i_b、i_c，变换到两相静止坐标系下的 i_α、i_β：

i_α = i_a
i_β = (i_a + 2*i_b) / √3

等等，这里有个坑。上面这个公式是等幅值变换，还有一种叫等功率变换，系数不一样。我在项目中遇到过有人把这两种搞混，结果算出来的转矩差了1.5倍，电机直接抖成筛子。

我个人习惯用等幅值变换，因为调试时看波形更直观。但如果你做的是功率计算，记得用等功率变换。两种变换的系数对比如下：

变换类型	α轴系数	β轴系数	适用场景
等幅值变换	1	1/√3	电流环控制、波形观察
等功率变换	√(2/3)	√(2/3)	功率计算、能量分析

3.3 Park变换：从αβ到dq

Clark变换做完，我们得到了αβ坐标系下的两相正弦波。但正弦波还是随时间变的，不好控。这时候Park变换上场了——它把静止的αβ坐标系，旋转起来，变成跟着磁场一起转的dq坐标系。

Park变换的公式：

i_d = i_α * cos(θ) + i_β * sin(θ)
i_q = -i_α * sin(θ) + i_β * cos(θ)

这里的θ是转子位置角，也就是d轴和α轴之间的夹角。在双馈风机里，这个θ通常来自编码器或者观测器。

你想想看，当坐标系以同步速旋转时，i_d和i_q就变成了直流量。d轴对应励磁分量，q轴对应转矩分量。这就是解耦控制的基础——调d轴控制无功，调q轴控制有功，互不干扰。

3.4 逆变换：从dq回到abc

控制算法算出来的是dq坐标系下的电压指令，但最终要给到逆变器的是三相电压。所以还得变回去：

• 逆Park变换：从dq变回αβ
• 逆Clark变换：从αβ变回abc

逆Park变换公式：

u_α = u_d * cos(θ) - u_q * sin(θ)
u_β = u_d * sin(θ) + u_q * cos(θ)

逆Clark变换公式（等幅值）：

u_a = u_α
u_b = (-u_α + √3 * u_β) / 2
u_c = (-u_α - √3 * u_β) / 2

嗯，这里要注意：逆变换的系数要和正变换保持一致。你正变换用了等幅值，逆变换也得用等幅值，否则算出来的电压幅值会不对。

3.5 坐标变换的完整流程

下面这张图展示了从三相电流采样到dq轴电流的完整流程：

3.6 代码实现示例

下面是我在DSP28335上实际跑过的代码片段，做了简化处理。你直接拿去用，注意改一下系数就行：

// Clark变换：三相电流 → αβ电流
void clark_transform(float ia, float ib, float ic, 
                     float *i_alpha, float *i_beta)
{
    *i_alpha = ia;                          // 等幅值变换
    *i_beta  = (ia + 2.0f * ib) / 1.732f;  // 1.732 ≈ √3
}

// Park变换：αβ电流 → dq电流
void park_transform(float i_alpha, float i_beta, float theta,
                    float *i_d, float *i_q)
{
    float cos_theta = cosf(theta);
    float sin_theta = sinf(theta);
    
    *i_d =  i_alpha * cos_theta + i_beta * sin_theta;
    *i_q = -i_alpha * sin_theta + i_beta * cos_theta;
}

// 逆Park变换：dq电压 → αβ电压
void inv_park_transform(float u_d, float u_q, float theta,
                        float *u_alpha, float *u_beta)
{
    float cos_theta = cosf(theta);
    float sin_theta = sinf(theta);
    
    *u_alpha = u_d * cos_theta - u_q * sin_theta;
    *u_beta  = u_d * sin_theta + u_q * cos_theta;
}

// 逆Clark变换：αβ电压 → 三相电压
void inv_clark_transform(float u_alpha, float u_beta,
                         float *u_a, float *u_b, float *u_c)
{
    *u_a = u_alpha;
    *u_b = (-u_alpha + 1.732f * u_beta) / 2.0f;
    *u_c = (-u_alpha - 1.732f * u_beta) / 2.0f;
}

3.7 常见问题与避坑

最后，我总结几个实际项目中容易踩的坑：

1. 角度对齐问题：Park变换用的θ必须是转子磁链位置，不是机械角度。在双馈风机里，这个角度需要根据转差频率和转子位置综合计算。
2. 变换系数一致性：正变换和逆变换的系数必须匹配。我曾经见过有人正变换用等幅值，逆变换用等功率，结果电压波形全是谐波。
3. 浮点精度：在定点DSP上做坐标变换，要注意三角函数查表法的精度。我建议用CORDIC算法或者查表+线性插值，精度能到0.1度以内。
4. 采样同步：三相电流采样必须在同一时刻完成，否则变换出来的dq分量会有纹波。这个坑我踩过，后来加了硬件同步采样才解决。

坐标变换这东西，看着是数学，其实是工程。你把它理解透了，矢量控制就入门了。剩下的就是调参数、调参数、再调参数——嗯，那是另一门学问了。

---

公众号:蓝海资料掘金营，微信deep3321