第3章 Clark变换：从三相到两相，把复杂问题简单化

好，咱们进入FOC的第一个数学变换——Clark变换。说实话，我第一次接触这个变换时，觉得它就是一堆矩阵公式，没啥意思。直到我在一个项目里，用三相电流直接做PID控制，结果电机抖得像筛子一样……后来才明白，Clark变换不只是数学游戏，它是把三相耦合系统解耦的关键一步。

3.1 为什么要做Clark变换？

你想想看，电机里有三根相线，电流分别是i_a、i_b、i_c。这三相电流是相互耦合的——你动其中一相，另外两相也跟着变。直接控制它们，就像同时拽三根绳子，很难让电机平稳转动。

Clark变换的核心思想很简单：把三相坐标系（a-b-c）投影到两相静止坐标系（α-β）上。说白了，就是用两个正交的轴（α轴和β轴）来等效描述三相系统。这样一来，控制对象从三个变成了两个，而且这两个轴是正交解耦的，控制起来就舒服多了。

3.2 变换矩阵推导

推导过程其实不复杂。咱们从几何投影的角度来看：

1. 定义α轴：让α轴与a相轴线重合。
2. 定义β轴：β轴逆时针旋转90°，与α轴正交。
3. 投影计算：把b相和c相的磁动势分别投影到α轴和β轴上。

三相绕组的空间位置相差120°，所以投影关系是这样的：

• α轴分量：i_α = i_a + i_b·cos(120°) + i_c·cos(240°)
• β轴分量：i_β = 0 + i_b·sin(120°) + i_c·sin(240°)

代入cos(120°) = -1/2，cos(240°) = -1/2，sin(120°) = √3/2，sin(240°) = -√3/2，得到：

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

写成矩阵形式就是：

[ iα ]   [ 1    -1/2   -1/2 ] [ ia ]
[ iβ ] = [ 0    √3/2  -√3/2 ] [ ib ]
                              [ ic ]

这就是等幅值Clark变换矩阵。注意，这里没有乘以2/3系数，所以变换后的α-β分量幅值与三相幅值相同。

3.3 等幅值 vs 等功率变换

这里有个容易踩的坑。Clark变换有两种常见形式：

我曾经在一个项目里，用了某芯片厂商的库，它默认是等功率变换。我按等幅值去调PI参数，结果电流环怎么都稳不住。查了两天才发现是变换系数不匹配……从那以后，我每次用别人的库，第一件事就是确认Clark变换的系数。

3.4 仿真验证

光说不练假把式。咱们用Python做个简单的仿真，看看Clark变换到底干了什么。

import numpy as np
import matplotlib.pyplot as plt

# 生成三相正弦电流（120°相位差）
t = np.linspace(0, 0.02, 1000)
f = 50  # 频率50Hz
ia = 5 * np.sin(2 * np.pi * f * t)
ib = 5 * np.sin(2 * np.pi * f * t - 2*np.pi/3)
ic = 5 * np.sin(2 * np.pi * f * t - 4*np.pi/3)

# Clark变换（等幅值）
i_alpha = ia - 0.5*ib - 0.5*ic
i_beta = (np.sqrt(3)/2)*ib - (np.sqrt(3)/2)*ic

# 绘制结果
plt.figure(figsize=(10, 6))
plt.subplot(2,1,1)
plt.plot(t*1000, ia, label='ia')
plt.plot(t*1000, ib, label='ib')
plt.plot(t*1000, ic, label='ic')
plt.ylabel('三相电流 (A)')
plt.legend()
plt.grid(True)

plt.subplot(2,1,2)
plt.plot(t*1000, i_alpha, label='iα', linewidth=2)
plt.plot(t*1000, i_beta, label='iβ', linewidth=2)
plt.xlabel('时间 (ms)')
plt.ylabel('两相电流 (A)')
plt.legend()
plt.grid(True)
plt.tight_layout()
plt.show()

运行这段代码，你会看到：三相电流是三个相位差120°的正弦波，而变换后的i_α和i_β是两个相位差90°的正弦波。而且，i_α的幅值和i_a完全一样——这就是等幅值变换的效果。

3.5 知识体系总览

下面这张图把Clark变换在整个FOC中的位置说清楚了：

// 等幅值Clark变换（浮点版本）
void clark_transform(float ia, float ib, float ic, 
                     float *i_alpha, float *i_beta) {
    *i_alpha = ia - 0.5f * ib - 0.5f * ic;
    *i_beta  = 0.8660254f * ib - 0.8660254f * ic;
    // 0.8660254 = sqrt(3)/2
}

// 简化版本（仅用于三相平衡系统）
void clark_transform_balanced(float ia, float ib, 
                              float *i_alpha, float *i_beta) {
    *i_alpha = ia;
    *i_beta  = 0.57735027f * ia + 1.15470054f * ib;
    // 推导：利用 ic = -ia - ib 代入原公式
}