gams.m

这是一个很适合研究和学习用的电力系统仿真软件

function [n,idx,data] = gams(a,type)

global Bus GAMS Rmpg

n = int2str(a.n);
idx = sparse(a.bus,[1:a.n],1,Bus.n,a.n);

data = [];

if ~a.n, return, end

[Csa,Csb,Csc,Dsa,Dsb,Dsc] = costs(a);
[Psmax,Psmin] = plim(a);

if GAMS.loaddir
  Ps0 = a.con(:,3);
else
  Ps0 = a.con(:,6);
end

ksu = a.con(:,15);

if GAMS.method < 7
  data.val = [Ps0,Psmax,Psmin,Csa,Csb,Csc,Dsa,Dsb,Dsc,ksu];
else
  Rgup = a.con(:,18);
  Rgdw = a.con(:,19);
  data.val = [Ps0,Psmax,Psmin,Csa,Csb,Csc,Dsa,Dsb,Dsc,ksu,Rgup,Rgdw];
end

if type == 2 | type == 4
  data.val = [data.val, zeros(a.n,8)];
  data.val(Rmpg.sup,[5:9,11,12]) = Rmpg.con(:,[9,5,6,3,4,7,8]);
  data.val(:,4) = Csb;
  data.val(:,10) = a.con(:,13);
  data.labels = {cellstr(num2str([1:a.n]')), ...
              {'Ps0','Psmax','Psmin','Cs', ...
               'suc','mut','mdt','rut','rdt','u0','y0','z0'}};
else
  if GAMS.method < 7
    data.labels = {cellstr(num2str([1:a.n]')), ...
                {'Ps0','Psmax','Psmin','Csa','Csb', ...
                 'Csc','Dsa','Dsb','Dsc','ksu'}};
  else
    data.labels = {cellstr(num2str([1:a.n]')), ...
                {'Ps0','Psmax','Psmin','Csa','Csb', ...
                 'Csc','Dsa','Dsb','Dsc','ksu','RGup','RGdw'}};
  end
end

data.name = 'S';