getxy.m

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

function [x,y] = getxy(a,bus,x,y)

if ~a.n, return, end

h1 = ismember(a.bus1,bus);
h2 = ismember(a.bus2,bus);
h = find(h1+h2);

if ~isempty(h)
  x = [x; a.alpha(h); a.Pm(h)];
end