tiebreaks.m

来自「这是一个很适合研究和学习用的电力系统仿真软件」· M 代码 · 共 23 行

function KTBS = tiebreaks(a)

global OPF

if ~OPF.tiebreak,
  KTBS = zeros(a.n,1);
else
  % a Tiebreaks
  if length(a.con(1,:)) < 14,
    KTBS = zeros(a.n,1);
  else,
    KTBS = a.con(:,14);
  end
  idx = find(a.con(:,4) == 0);
  if ~isempty(idx)
    KTBS(idx) = 0;
  end
  idx = find(a.con(:,4));
  if ~isempty(idx)
    KTBS(idx) = KTBS(idx)./a.con(idx,4);
  end
end