pqdisplay.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 16 行

function idx = pqdisplay(a)

global PV SW

idx = 0;

for i = 1:a.n
  bpv = findbus(PV,a.bus(i));
  bsw = findbus(SW,a.bus(i));
  bpq = a.u(i)*a.con(i,4) > 0;
  if isempty(bpv) & isempty(bsw) & bpq
    idx = a.bus(i);
    break
  end
end