gyreactive.m

来自「基于PSAT 软件的多目标最优潮流计算用于中小型电力系统的分析和管理」· M 代码 · 共 13 行

function Gyreactive(p)

global DAE

if ~p.n, return, end

DAE.J21(p.bus,:) = 0;
DAE.J22(p.bus,:) = 0;
DAE.J12(:,p.bus) = 0;
DAE.J22(:,p.bus) = 0;
DAE.J22(p.bus,p.bus) = DAE.J22(p.bus,p.bus) + speye(p.n);