pqshunt.m

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

function y = pqshunt(a,idx)

global DAE Shunt Bus

g = a.con(idx,4)./(DAE.V(a.bus(idx)).^2);
b = -a.con(idx,5)./(DAE.V(a.bus(idx)).^2);
Shunt.g(a.bus(idx)) = Shunt.g(a.bus(idx)) + g;
Shunt.b(a.bus(idx)) = Shunt.b(a.bus(idx)) + b;
y = sparse(double(a.bus(idx)),double(a.bus(idx)),g+i*b,Bus.n,Bus.n);