ybus.m

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

function y = ybus(a,buslist)

global Bus

nb = Bus.n;
y = sparse(nb,nb);

if ~a.n, return, end

idx = [];
for i = 1:a.n
  jdx = find(buslist ~= a.bus(i));
  if ~isempty(jdx) & a.u(i), idx = [idx; i]; end
end

if isempty(idx), return, end

xdb = 0.5*(a.con(idx,10)+a.con(idx,15));
rdb = a.con(idx,7);
zdb = rdb + sqrt(-1)*xdb;
ydb = 1/zdb;

y = sparse(a.bus(idx),a.bus(idx),ydb,nb,nb);