在风力发电领域，双馈感应发电机(DFIG) 凭借“变速恒频、灵活无功调节”成为兆瓦级风电机组的主流机型（如维斯塔斯V52、金风科技GW 1.5MW）。其核心优势是：定子直接并网（工频），转子通过背靠背变频器馈电（转差功率双向流动），可在宽风速范围内实现最大功率捕获(MPPT)。

定子侧电压定向控制(VOC) 是DFIG的经典控制策略，核心思想是：将同步旋转坐标系的d轴与定子电压矢量重合，使定子电压q轴分量为零(usq=0)，从而实现有功功率Ps与无功功率Qs的解耦控制：

有功控制：通过d轴转子电流ird调节，跟踪MPPT给定功率；
无功控制：通过q轴转子电流irq调节，实现并网功率因数校正（如单位功率因数运行）。

核心价值：

解耦控制：有功/无功独立调节，动态响应快（电流环带宽>100Hz）；
并网性能优：输出电流谐波畸变率THD<3%（满足IEEE 1547标准）；
工程成熟度高：占DFIG风电机组控制的70%以上（尤其早期机型）。

挑战：

坐标变换复杂度：需掌握定子/转子坐标系的Clark/Park变换；
参数敏感性：转子漏感、互感漂移影响解耦精度（需参数辨识补偿）；
转子侧协调：定子侧VOC需与转子侧电流控制协同（双端变频协调）。

✅ 本文目标：从零搭建DFIG定子侧VOC控制系统，掌握“电压定向-功率解耦-双环控制”核心逻辑，验证其在风速波动下的MPPT跟踪与无功调节性能（有功误差<5%、无功误差<3%、THD<3%），为DFIG进阶控制（如矢量控制、直接功率控制）奠定基础。

二、核心原理：定子侧VOC的“三大支柱”

1. 控制思想：“电压定向→功率解耦→双环跟踪”

定子侧VOC通过“坐标定向+功率解耦”实现简化控制，具体分三步：

电压定向：将同步旋转坐标系(dq)的d轴与定子电压矢量us重合，使usd=Us（定子电压幅值）、usq=0（q轴分量为零），此时定子磁链ψs与d轴重合；
功率解耦：基于定向后的电压方程，推导出有功Ps仅与ird相关、无功Qs仅与irq相关（解耦表达式见下文）；
双环跟踪：外环（功率环）根据MPPT输出有功给定Ps∗、无功给定Qs∗，计算转子电流参考ird∗/irq∗；内环（电流环）跟踪转子电流，SVPWM驱转子侧变频器执行。

2. 数学原理：DFIG在定子电压定向下的模型

（1）DFIG基本结构与坐标系定义

定子侧：三相绕组（a-b-c）接电网（电压usa,usb,usc），频率fs=50Hz；
转子侧：三相绕组（a'-b'-c'）接背靠背变频器（转子电压ura,urb,urc），频率fr=sfs（s为转差率，s=(ns−nr)/ns，ns为同步转速，nr为转子转速）；
坐标系：同步旋转坐标系(dq)（与定子电压同步旋转）、转子旋转坐标系(d'r'q')（与转子同步旋转）。

（2）定子电压定向下的电压方程（关键推导）

将d轴与定子电压矢量重合（usq=0），忽略定子电阻Rs（工频下Rs≪Xs），定子电压方程简化为：

{usd=Us=ωsψsq≈0?修正：实际应为 usd=Us,usq=0ψsd=ψs(定子磁链幅值),ψsq=0

结合转子电压方程和磁链方程，最终推导出有功/无功解耦表达式（核心！）：

{Ps=23LsUsLmird(有功仅与转子d轴电流ird相关)Qs=23LsUs(Ls−Lm)irq(无功仅与转子q轴电流irq相关)

其中，Lm为定转子互感，Ls=Lls+Lm为定子总电感（Lls为定子漏感），Us为定子电压幅值。

（3）双环控制结构（外环功率+内环电流）

外环（功率环）：输入MPPT给定Ps∗、无功给定Qs∗（如Qs∗=0单位功率因数），通过PI控制器输出转子电流参考：
```
ird∗=3UsLm2LsPs∗,irq∗=3Us(Ls−Lm)2LsQs∗
```

内环（电流环）：跟踪ird∗/irq∗，输出转子电压参考urd∗/urq∗，经坐标变换后驱动转子侧SVPWM。

三、应用场景：兆瓦级DFIG风电机组变速恒频发电

场景设定

风机参数：叶轮半径r=35m，空气密度ρ=1.225kg/m3，最佳叶尖速比λopt=7，风能利用系数Cp_max=0.45；
DFIG参数：额定功率1.5MW，定子电压Us=690V（线电压），频率fs=50Hz，极对数p=2，同步转速ns=1500rpm，转子转速范围nr=1000∼1800rpm（转差率s=−0.33∼0.2），定子漏感Lls=0.1mH，转子漏感Llr=0.12mH，互感Lm=3mH，转动惯量J=100kg⋅m2；
控制需求：
- MPPT跟踪：风速vw=8m/s时捕获最大功率Pe_max≈800kW；
- 无功调节：并网功率因数cosφ=1（Qs=0）；
- 并网质量：输出电流THD<3%，转速跟踪误差<2%。

四、Simulink建模步骤（附核心代码与模块）

1. 主电路搭建（Simscape Electrical）

模块	作用	参数设置
`Wind Turbine Model`	风机气动模型（桨距角固定）	ρ=1.225kg/m3,r=35m,λopt=7,Cp_max=0.45
`Gearbox`	齿轮箱（增速比10:1）	增速比G=10（风机转速→DFIG转子转速）
`Doubly-Fed Induction Generator`	DFIG模型	p=2,Us=690V,fs=50Hz,Lls=0.1mH,Llr=0.12mH,Lm=3mH,J=100kg⋅m2
`Grid`	无穷大电网	线电压690V，频率50Hz，内阻0.01Ω
`Back-to-Back Converter`	转子侧变频器（IGBT）	开关频率5kHz，直流母线电压Vdc=1200V

2. 定子侧VOC核心模块实现（MATLAB Function代码）

（1）坐标变换模块（Clark+Park变换）

功能：将三相电压/电流转换为dq轴分量（定子侧定向用）。

function [usd, usq, isd, isq] = abc_to_dq_us_oriented(usa, usb, usc, isa, isb, isc, theta_s)
    % 输入：定子三相电压(usa,usb,usc)、电流(isa,isb,isc)，同步角度theta_s（定子电压定向角度）
    % 输出：dq轴电压(usd,usq)、电流(isd,isq)（d轴与us重合，usq=0）
    
    % Clark变换（3s→2s）
    us_alpha = (2/3)*(usa - 0.5*usb - 0.5*usc);
    us_beta = (2/3)*(0 + sqrt(3)/2*usb - sqrt(3)/2*usc);
    is_alpha = (2/3)*(isa - 0.5*isb - 0.5*isc);
    is_beta = (2/3)*(0 + sqrt(3)/2*isb - sqrt(3)/2*isc);
    
    % Park变换（2s→2r，d轴与us重合，theta_s为定子电压角度）
    usd = us_alpha*cos(theta_s) + us_beta*sin(theta_s);  % usd=Us（幅值），usq=0（定向保证）
    usq = -us_alpha*sin(theta_s) + us_beta*cos(theta_s);
    isd = is_alpha*cos(theta_s) + is_beta*sin(theta_s);
    isq = -is_alpha*sin(theta_s) + is_beta*cos(theta_s);
end

（2）功率计算与电流参考生成模块

功能：根据Ps∗/Qs∗计算转子电流参考ird∗/irq∗（解耦公式）。

function [ird_star, irq_star] = power_to_ird_irq(Ps_star, Qs_star, Us, Ls, Lm, Ts)
    % 输入：有功给定Ps_star(W)，无功给定Qs_star(var)，定子电压幅值Us(V)，定子总电感Ls(H)，互感Lm(H)
    % 输出：转子d/q轴电流参考ird_star, irq_star(A)
    Ls_total = Ls;  % Ls = Lls + Lm（定子总电感）
    ird_star = (2 * Ls_total * Ps_star) / (3 * Us * Lm);  % 有功→ird参考
    irq_star = (2 * Ls_total * Qs_star) / (3 * Us * (Ls_total - Lm));  % 无功→irq参考
    % 电流限幅（转子额定电流50A）
    ird_star = max(min(ird_star, 50), -50);
    irq_star = max(min(irq_star, 50), -50);
end

（3）转子电流环PI控制器（含解耦补偿）

功能：跟踪ird∗/irq∗，输出转子电压参考urd∗/urq∗。

function [urd_star, urq_star] = rotor_current_pi(ird_ref, irq_ref, ird_act, irq_act, wr, Llr, Lm, Rs, Rr, Ts)
    % 输入：转子电流参考/实际值，转子电角速度wr(rad/s)，转子漏感Llr，互感Lm，定转子电阻Rs/Rr
    % 输出：转子d/q轴电压参考urd_star, urq_star(V)
    
    % 转子d轴电流环PI（含解耦补偿）
    persistent int_ird; if isempty(int_ird), int_ird=0; end
    e_ird = ird_ref - ird_act;
    int_ird = int_ird + e_ird*Ts;
    urd_pi = Kp_r * e_ird + Ki_r * int_ird;  % PI输出
    urd_comp = wr * Llr * irq_act + wr * Lm * isq;  % 解耦补偿（抵消交叉耦合）
    urd_star = urd_pi + urd_comp;
    
    % 转子q轴电流环PI（含解耦补偿）
    persistent int_irq; if isempty(int_irq), int_irq=0; end
    e_irq = irq_ref - irq_act;
    int_irq = int_irq + e_irq*Ts;
    urq_pi = Kp_r * e_irq + Ki_r * int_irq;  % PI输出
    urq_comp = -wr * Llr * ird_act - wr * Lm * isd;  % 解耦补偿
    urq_star = urq_pi + urq_comp;
    
    % 电压限幅（≤直流母线电压1200V）
    urd_star = max(min(urd_star, 1200), -1200);
    urq_star = max(min(urq_star, 1200), -1200);
end

（4）SVPWM生成模块（转子侧变频器）

功能：将转子电压参考urd∗/urq∗转换为三相PWM信号。

function [ura, urb, urc] = svpwm_rotor(urd, urq, theta_r)
    % 输入：转子d/q轴电压urd, urq(V)，转子角度theta_r(rad)
    % 输出：转子三相电压ura, urb, urc(V)
    % 逆Park变换（d/q→α/β）
    ur_alpha = urd*cos(theta_r) - urq*sin(theta_r);
    ur_beta = urd*sin(theta_r) + urq*cos(theta_r);
    % 逆Clark变换（α/β→a/b/c）
    ura = ur_alpha;
    urb = -0.5*ur_alpha + sqrt(3)/2*ur_beta;
    urc = -0.5*ur_alpha - sqrt(3)/2*ur_beta;
end