getxy.m

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

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

if ~a.n, return, end

h = find(ismember(a.bus,bus));

if ~isempty(h)
  x_temp = [a.slip(h); a.e1r(h); a.e1m(h); a.e2r(h); a.e2m(h)];
  idx = find(x_temp);
  x = [x; x_temp(idx)];
end