gams.m

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

function [n,data,Li,Lj,Gh,Bh,Ghc,Bhc] = gams(a,method)

global Bus GAMS Settings

n = int2str(a.n);
nl = [1:a.n];
nb = Bus.n;

% Flows on transmission lines
tps = a.con(:,11).*exp(j*a.con(:,12)*pi/180);
r = a.con(:,8);
x = a.con(:,9);
chrg =  a.u.*a.con(:,10)/2;
z = r + j*x;
y = a.u./z;
g = real(y);
b = imag(y);
m = tps.*conj(tps);
y1 = tps.*y./m;
g1 = real(y1);
b1 = imag(y1);
g0 = g./m;

if GAMS.flow == 1 | GAMS.flow == 3
  b0 = chrg;
else
  b0 = chrg-b./m;
end

switch method
 case 1
  Li = [];
  Lj = [];  
 case 2
  Li = sparse(nl,a.fr,b,a.n,nb);
  Lj = sparse(nl,a.to,b,a.n,nb);
 otherwise
  Li = sparse(nl,a.fr,1,a.n,nb);
  Lj = sparse(nl,a.to,1,a.n,nb);
end

if method ~= 1
  Gh = real(a.Y);
  Bh = imag(a.Y);
  Gh(1,1) = Gh(1,1) + 1e-8;
  Bh(1,1) = Bh(1,1) + 1e-8;
else
  Gh = [];
  Bh = [];
end  

switch method
 case {4,6,7}
  if GAMS.line
    a.u(GAMS.line) = 0;
    a = build_y(a);
    Ghc = real(a.Y);
    Bhc = imag(a.Y);
  else
    Ghc = Gh;
    Bhc = Bh;
  end
 otherwise
  Ghc = [];
  Bhc = [];
end

if GAMS.flow
  Pijmax = getflowmax(a,GAMS.flow);
  Pjimax = Pijmax;
else
  Pijmax = 999*Settings.mva*ones(a.n,1);
  Pjimax = 999*Settings.mva*ones(a.n,1);
end

data.val = [g1,b1,g0,b0,Pijmax,Pjimax];
data.labels = {cellstr(num2str([1:a.n]')), ...
	    {'g','b','g0','b0','Pijmax','Pjimax'}};
data.name = 'N';