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.iqr,a.dat(:,8),a.dat(:,9))
fm_windup(a.idr,a.dat(:,10),a.dat(:,11))