getxy.m

电力系统分析计算程序

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.delta(h); a.omega(h); a.e1q(h); a.e1d(h); ...
       a.e2q(h); a.e2d(h); a.psiq(h); a.psid(h)]; 
  idx = find(x_temp);
  x = [x; x_temp(idx)];
  y = [y; a.vf(h); a.pm(h); a.p(h); a.q(h)];
end