pqshunt.m

这是一个很适合研究和学习用的电力系统仿真软件

function y = pqshunt(a)

global DAE Settings Shunt Bus

y = sparse(Bus.n,Bus.n);

if ~Settings.pq2z | Settings.init > 1, return, end

if ~a.n, return, end

idx = find(a.gen == 0);

if isempty(idx), return, end

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);