双馈风力发电系统AC/DC/AC变流器的仿真模型 ，网侧变流器及转子侧变流器的控制方案，采用双PWM控制，网侧采用定子电压定向，转子侧采用定子磁链定向，实现超同步及次同步时功率的双向流动。 还有20+页的说明内容。

双馈风机变流器的控制就像在玩一场精密的速度游戏。网侧变流器和转子侧变流器这对搭档，一个盯着电网电压的节奏，一个追着磁链的方向跑，双PWM控制就是他们手里的指挥棒。咱们先扒开网侧变流器的代码看看——这货用的是定子电压定向，简单说就是把电网电压坐标系锁死在d轴上。

下面这段电压定向的Simulink函数是关键：

function [theta] = VoltageOrientation(U_alpha, U_beta) theta = atan2(U_beta, U_alpha); % 这玩意儿实时计算电压矢量角度 % 注意这里用了四象限反正切，防止角度跳变 end

电网电压的alpha-beta分量经过这个黑盒子，输出当前电压矢量的相位角。这个角度直接喂给Park变换，把交流量转成直流量控制。这时候的电流内环PI参数设置特别讲究，我一般先用临界比例法估算个大概：

Kp = 0.5 * L_filter / Ts; Ki = 0.5 * R_filter / Ts;

L_filter是网侧电抗，Ts是控制周期。实际调参时发现，当风速突变时，Kp值加个0.2能更好抑制直流母线电压波动。

转子侧那边玩的是定子磁链定向，代码里藏着玄机。磁链观测器的实现方式直接影响控制性能，这里用了电压模型和电流模型混合观测：

def flux_observer(u_alpha, u_beta, i_alpha, i_beta, R_s, L_s): psi_alpha = integrate(u_alpha - R_s*i_alpha) - L_s*i_alpha psi_beta = integrate(u_beta - R_s*i_beta) - L_s*i_beta return np.sqrt(psi_alpha**2 + psi_beta**2), np.arctan2(psi_beta, psi_alpha)

积分环节用梯形法代替欧拉法，防止低速时观测值漂移。仿真时发现，当转速低于同步速（次同步状态），转子电流相位要反过来跑，这时候前馈补偿项得加个负号，不然功率流动方向就乱套了。

超同步运行时的控制更有意思。这时候转子得往外送功率，网侧变流器的有功指令得从正值切到负值。代码里用了个平滑过渡函数：

float power_ramp(float target, float current, float rate) { float delta = target - current; if(fabs(delta) < rate) return target; return current + (delta > 0 ? rate : -rate); }

仿真模型里的直流母线电容取值是个坑，太小了电压波动大，太大了动态响应慢。试了N次发现，按能量守恒原则计算初始值再乘个1.2倍刚好合适。最后跑起来看波形，功率双向流动时网侧电流THD控制在3%以内，算是达标了。