windup.m

电力系统分析计算程序

function windup(a)

if ~a.n, return, end

fm_windup(a.theta_p,999*ones(a.n,1),zeros(a.n,1))
fm_windup(a.iqs,a.con(:,21),a.con(:,22))
fm_windup(a.idc,-a.con(:,24),-a.con(:,23))