4. dq坐标系下的DFIG模型：解耦控制的理论基石
各位工程师朋友，咱们今天来啃一块硬骨头——两相旋转坐标系下的DFIG数学模型。
说实话，我刚接触DFIG控制那会儿，看到abc三相模型就头大。一堆时变系数，耦合得跟麻花似的。后来前辈告诉我：“你试试dq变换，世界就清净了”。嗯，确实如此。
这一节，我就带你一步步推导dq坐标系下的DFIG模型。别怕，咱们不搞纯数学推导，我会结合项目中的实际体会来讲。
4.1 为什么要做坐标变换？
先问个问题：三相静止坐标系下的DFIG模型，有什么毛病？
我直接说结论：电感矩阵是时变的。
你想想看，转子在转，定转子之间的互感就在不停变化。这导致电压方程里全是sin、cos函数，控制器设计起来简直噩梦。我在做第一个DFIG仿真项目时，直接用abc模型搭PI控制器，结果电流波形抖得像心电图——根本稳不住。
后来我明白了：必须把交流量变成直流量。这就是坐标变换的核心思想。
核心思路：
- abc三相静止 → αβ两相静止（Clark变换）
- αβ两相静止 → dq两相旋转（Park变换）
- 最终效果：基波分量变成直流，控制难度直线下降
4.2 坐标变换的数学基础
我个人习惯用等幅值变换，也就是变换前后电压电流幅值不变。这在工程中更直观。
Clark变换（abc → αβ）：
// 等幅值Clark变换矩阵
T_clark = (2/3) * [1, -1/2, -1/2;
0, √3/2, -√3/2;
1/2, 1/2, 1/2]Park变换（αβ → dq）：
// Park变换矩阵，θ是d轴与α轴的夹角
T_park = [cosθ, sinθ;
-sinθ, cosθ]这里有个关键点：θ怎么选？
在DFIG中，我们通常让d轴定向于定子磁链矢量。为什么？因为这样q轴分量就和有功功率直接挂钩了——这就是后面解耦控制的基础。
我的经验：刚开始做仿真时，我直接用电网电压定向，结果发现转子电流控制效果很差。后来改成定子磁链定向，动态响应明显改善。记住：DFIG的天然定向基准是定子磁链，不是电网电压。
4.3 dq坐标系下的DFIG电压方程
好，现在开始推导核心内容。咱们从三相静止坐标系下的电压方程出发：
定子电压方程（abc坐标系）：
u_sabc = R_s * i_sabc + d(ψ_sabc)/dt转子电压方程（abc坐标系）：
u_rabc = R_r * i_rabc + d(ψ_rabc)/dt经过Clark+Park变换后，得到dq坐标系下的方程：
定子电压（dq坐标系）：
u_sd = R_s * i_sd + d(ψ_sd)/dt - ω_s * ψ_sq
u_sq = R_s * i_sq + d(ψ_sq)/dt + ω_s * ψ_sd转子电压（dq坐标系）：
u_rd = R_r * i_rd + d(ψ_rd)/dt - ω_slip * ψ_rq
u_rq = R_r * i_rq + d(ψ_rq)/dt + ω_slip * ψ_rd其中：
- ω_s：同步角速度（定子磁场转速）
- ω_slip = ω_s - ω_r：转差角速度（这个量在DFIG控制中极其重要）
注意：转子方程里出现的是ω_slip，不是ω_s！我曾经在调试时写错了这个参数，结果转子电流控制完全发散。后来查了一整天，才发现是转差频率搞错了。这个坑，我替你们踩过了。
4.4 dq坐标系下的磁链方程
磁链方程在dq坐标系下变得非常简洁——电感矩阵不再是时变的了。
定子磁链：
ψ_sd = L_s * i_sd + L_m * i_rd
ψ_sq = L_s * i_sq + L_m * i_rq转子磁链：
ψ_rd = L_r * i_rd + L_m * i_sd
ψ_rq = L_r * i_rq + L_m * i_sq其中：
- L_s = L_ls + L_m：定子自感
- L_r = L_lr + L_m：转子自感
- L_m：激磁电感（定转子互感）
你看，所有电感都是常数了。这就是坐标变换的魔力——把时变系统变成了定常系统。
4.5 定子磁链定向下的简化模型
前面说了，我们通常让d轴定向于定子磁链矢量。这意味着：
ψ_sd = |ψ_s| (定子磁链幅值)
ψ_sq = 0 (q轴磁链分量为零)在这个条件下，模型可以大幅简化。我直接给出工程中常用的形式：
定子电压简化：
u_sd ≈ 0
u_sq ≈ ω_s * ψ_sd ≈ 电网电压幅值有功功率：
P_s = - (3/2) * u_sq * i_sq ≈ - (3/2) * U_s * i_sq无功功率：
Q_s = (3/2) * u_sq * i_sd ≈ (3/2) * U_s * i_sd关键结论：
- i_sd控制无功功率（励磁分量）
- i_sq控制有功功率（转矩分量）
- 定子电流的d、q分量实现了自然解耦
说白了，你只要控制i_sd和i_sq，就能独立控制DFIG的无功和有功。这就是功率解耦控制的理论基础。
4.6 转子侧控制的关键方程
对于转子侧变流器控制，我们需要知道转子电压和转子电流的关系。结合磁链方程和电压方程，可以得到：
转子电压-电流关系：
u_rd = R_r * i_rd + σ * L_r * d(i_rd)/dt - ω_slip * σ * L_r * i_rq
u_rq = R_r * i_rq + σ * L_r * d(i_rq)/dt + ω_slip * (L_m²/L_s * i_ms + σ * L_r * i_rd)其中：
- σ = 1 - L_m²/(L_s * L_r)：漏磁系数
- i_ms = ψ_sd / L_m：等效励磁电流
看到没？d轴和q轴之间仍然存在耦合项（带ω_slip的项）。这就是为什么我们需要前馈补偿来实现完全解耦。
我的做法：在实际控制器中，我会把耦合项作为前馈量加到PI输出上。比如控制u_rd时，把-ω_slip * σ * L_r * i_rq这个项直接加上去。这样PI控制器只需要处理纯积分环节，参数整定简单多了。
4.7 知识体系总览
为了让你对整个推导过程有个全局认识，我画了张图：
4.8 工程中的注意事项
推导完了，我再说几个实际项目中容易踩的坑：
- 坐标变换的符号约定要统一——不同文献可能用不同符号，我建议全程用发电机惯例（电流流出为正）。
- θ角的初始值要校准——编码器安装位置偏差会导致θ有固定偏移，需要在调试时补偿掉。
- 滤波器延迟会影响解耦效果——采样信号经过滤波器会有相位延迟，我一般用同步旋转坐标系下的低通滤波器来避免这个问题。
- 弱磁工况下模型会变——当转速超过额定值时，需要进入弱磁区，此时定子磁链定向不再适用，需要切换到电压定向。
曾经踩过的坑：我在做一台2MW双馈风机调试时，发现无功功率控制总是有静差。查了两天，最后发现是Park变换中的θ角用了电网电压相位，而不是定子磁链相位。两者在稳态时接近，但动态过程中差异很大。换成磁链观测器后，问题立刻解决。
4.9 小结
这一节我们完成了DFIG从abc到dq坐标系的完整推导。核心收获有三点：
- 坐标变换把时变系统变成了定常系统，控制器设计难度大幅降低
- 定子磁链定向实现了有功/无功的自然解耦，i_sd控无功，i_sq控有功
- 转子侧仍然存在耦合项，需要通过前馈补偿实现完全解耦
有了这个dq模型，下一节我们就可以开始设计具体的控制器了。嗯，不过那是下一章的事。今天先把模型吃透，后面控制器的参数整定才能有的放矢。