pqshunt.m

用于电力系统的一个很好的分析软件

function a = pqshunt(a)

global DAE Settings

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

if ~a.n, return, end

a.shunt = ~a.gen;
idx = find(a.shunt);
if isempty(idx), return, end
a.con(idx,7) = DAE.y(a.vbus(idx));
a.con(idx,8) = 0;