2. 坐标变换基础：Clark变换、Park变换、反Park变换的数学原理与实现

大家好，我是你们的老朋友。今天咱们来聊聊电机控制里最基础、也最绕不开的一块内容——坐标变换。

说实话，我刚入行那会儿，看到Clark、Park这些名字，第一反应是：这又是哪个数学家在折腾我们工程师？后来做项目做多了才明白，没有这些变换，PMSM根本玩不转。你想想看，电机里转的是三相正弦电流，可我们要控制的是转矩和磁通，这俩东西在静止坐标系下是耦合在一起的，解耦？不存在的。

所以，坐标变换说白了就是换个角度看问题。把复杂的、耦合的物理量，映射到一个我们好下手的坐标系里。今天我就带大家把这三个变换的数学原理和代码实现捋一遍。

2.1 为什么需要坐标变换？

先问个问题：PMSM的定子绕组上通的是三相交流电，对吧？这三相电流在空间上相差120度，时间上也相差120度。结果就是，它们合成一个旋转的磁动势。

但问题来了——这个旋转磁动势的大小和方向，跟转子位置、负载都有关系。如果我们直接在ABC三相坐标系下建模型，那方程里全是时变系数，解起来头大。

我当年在学校里推导PMSM的电压方程，推了满满两页纸，最后发现全是正弦余弦的乘积项。后来导师跟我说：你试试坐标变换。一试，嘿，方程简洁得像教科书里的例题。

所以，坐标变换的核心目的就两个：

• 解耦：把强耦合的物理量拆开，让转矩和磁通可以独立控制
• 简化：把时变系数变成常数，方便设计控制器

2.2 Clark变换：从ABC到αβ

Clark变换，也叫3s/2s变换。它把三相静止坐标系（ABC）映射到两相静止坐标系（αβ）。

数学上很简单，就是投影。假设三相电流分别是 i_a、i_b、i_c，那么：

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

等等，这里有个细节。上面这个公式是等幅值变换，也就是变换前后电流的幅值不变。还有一种叫等功率变换，系数会不一样。我个人习惯用等幅值变换，因为调试的时候看电流波形更直观。

完整的矩阵形式是这样的：

[ i_α ]   [ 1      -1/2    -1/2   ] [ i_a ]
[ i_β ] = [ 0      √3/2   -√3/2  ] [ i_b ]
[ i_0 ]   [ 1/√2   1/√2    1/√2  ] [ i_c ]

注意那个 i₀，它是零序分量。对于星形连接且中性点不接地的电机，i_a + i_b + i_c = 0，所以 i₀ 恒为0。实际代码里我们经常直接忽略它。

2.3 Park变换：从αβ到dq

Clark变换完了，我们得到了αβ坐标系下的两相正弦电流。但正弦信号还是不好控，对吧？PID控制器对直流信号的跟踪效果最好。

Park变换就是干这个的——它把静止的αβ坐标系，旋转到跟转子同步的dq坐标系。说白了，就是把交流变直流。

数学公式：

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

这里的 θ 是转子电角度，一般由编码器或霍尔传感器提供。

变换之后，i_d 对应励磁分量，i_q 对应转矩分量。你看，是不是清爽多了？控制 i_d 就能控制磁通，控制 i_q 就能控制转矩，互不干扰。

2.4 反Park变换：从dq回到αβ

控制器算出来的是 dq 坐标系下的电压指令 u_d* 和 u_q*，但最终要给到逆变器的是三相电压。所以我们需要反变换回去。

反Park变换就是Park变换的逆过程：

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

你看，跟Park变换的公式很像，只是符号位置换了一下。我记这个公式有个小窍门：正变换是「cos sin, -sin cos」，反变换是「cos -sin, sin cos」，矩阵是正交的，所以逆矩阵就是转置。

得到 u_α 和 u_β 之后，再经过反Clark变换（或者直接用SVPWM模块），就能生成三相PWM波了。

2.5 代码实现示例

下面给出一段C语言实现的代码，我在STM32F4上跑过，效果不错。注意用了浮点运算，如果资源紧张可以改成定点数。

// 坐标变换结构体
typedef struct {
    float alpha;
    float beta;
} Clarke_t;

typedef struct {
    float d;
    float q;
} Park_t;

// Clark变换
Clarke_t Clark_Transform(float ia, float ib, float ic) {
    Clarke_t out;
    out.alpha = ia;
    out.beta  = (ia + 2.0f * ib) / 1.7320508f;  // 1/√3
    return out;
}

// Park变换
Park_t Park_Transform(float alpha, float beta, float theta) {
    Park_t out;
    float sin_theta = sinf(theta);
    float cos_theta = cosf(theta);
    out.d =  alpha * cos_theta + beta * sin_theta;
    out.q = -alpha * sin_theta + beta * cos_theta;
    return out;
}

// 反Park变换
Clarke_t InvPark_Transform(float d, float q, float theta) {
    Clarke_t out;
    float sin_theta = sinf(theta);
    float cos_theta = cosf(theta);
    out.alpha = d * cos_theta - q * sin_theta;
    out.beta  = d * sin_theta + q * cos_theta;
    return out;
}

2.6 三种变换的关系图

为了让大家更直观地理解这三个变换之间的关系，我画了一张流程图：

_a, i_b, i_c 两相静止 i_α, i_β 两相旋转 i_d, i_q 反Park变换（电压） 转子电角度 θ 电流环中：采样→Clark→Park→PI调节→反Park→SVPWM

从这张图可以看得很清楚：电流从传感器出来，先经过Clark变换到αβ坐标系，再经过Park变换到dq坐标系。控制器在dq坐标系下算完电压指令后，先反Park变换回αβ坐标系，再通过SVPWM生成三相电压。整个过程，θ贯穿始终。

2.7 常见问题与避坑

最后，我把自己这些年踩过的坑整理一下，希望对大家有帮助：

问题	现象	原因	解决办法
角度延迟	高速时电流波形畸变	角度采样滞后于电流采样	使用角度预测或补偿
零漂误差	dq轴有直流偏置	电流传感器零漂未校准	上电时做零漂校准
变换系数错误	转矩输出偏小	等幅值/等功率混用	统一用一种，推荐等幅值
三角函数精度	低速时抖动	查表法分辨率不够	增加查表点数或插值

特别提醒：我曾经在一个项目里，因为用了不同版本的数学库，导致sin/cos的精度不一致，结果电机在低速时总是有轻微的转矩脉动。排查了整整两天才发现是库函数的问题。所以，请务必确保整个工程使用同一套数学库。

好了，关于坐标变换的基础知识就讲到这里。这些内容虽然看起来简单，但它们是整个电流环的基石。我建议大家动手写一遍代码，在示波器上看看变换前后的波形，感受一下「交流变直流」的神奇过程。

---

← 上一章 📖 返回目录 下一章 →