getxy.m

来自「电力系统分析计算程序」· M 代码 · 共 15 行

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

global Wind

if ~a.n, return, end

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

if ~isempty(h)
  vw = Wind.vw(a.wind(h));
  ws = Wind.ws(a.wind(h));
  x = [x; a.theta_p(h); a.omega_m(h); a.idr(h); a.iqr(h); vw];
  y = [y; ws];
end