matpower2psat.m

电力系统分析计算程序

function check = matpower2psat(filename, pathname)
% MATPOWER2PSAT Matpower 2 PSAT filter
%
% CHECK = MATPOWER2PSAT(FILENAME,PATHNAME)
%       FILENAME name of the file to be converted
%       PATHNAME path of the file to be converted
%
%       CHECK = 1 conversion completed
%       CHECK = 0 problem encountered (no data file created)
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://www.uclm.es/area/gsee/Web/Federico
%
%*** Contributions ***
%
%Author:    Liulin
%Date:      13-Mar-2003
%Update:    01-May-2004
%E-mail:    [e-mail]
%
% Copyright (C) 2002-2007 Federico Milano

global Settings Path

check = 0;
pathname = [pathname,filesep];

fm_disp
fm_disp('Conversion from Matpower Data Format ...');
fm_disp(['Source data file "',pathname,filename,'"'])

baseMVA = [];
bus = [];
gen = [];
branch = [];
area = [];
gencost = [];

try
  cd(pathname)
  [baseMVA, bus, gen, branch, area, gencost] = feval(filename(1:end-2));
  cd(Path.local)
catch
  fm_disp(['Something wrong with the file ',pathname,filename])
  fm_disp('Conversion Process Interrupted',2)
  return
end
lasterr('')

try
  if baseMVA == 0, baseMVA = 100; end
  gen(find(gen(:,7)== 0),7) = baseMVA;
  bus(:,[3 4 5 6]) = bus(:,[3 4 5 6])/baseMVA;
  if sum(bus(:,10)) == 0, bus(:,10) = 1; end
  bus(:,9) = pi*bus(:,9)/180;
  gen(:,2) = gen(:,2)./gen(:,7);
  gen(:,3) = gen(:,3)./gen(:,7);
  gen(:,4) = gen(:,4)./gen(:,7);
  gen(:,5) = gen(:,5)./gen(:,7);
  gen(:,9) = gen(:,9)./gen(:,7);
  gen(:,10) = gen(:,10)./gen(:,7);
  branch(:,[6 7 8]) = branch(:,[6 7 8])/baseMVA;
  if bus,
    if length(bus(1,:)) ~= 13,
      error('Bus data is not in the standard format'),
    end,
  end
  if gen,
    if length(gen(1,:)) ~= 10,
      error('Generator data is not in the standard format'),
    end,
  end
  if branch,
    if length(branch(1,:)) ~= 11,
      error('Branch data is not in the standard format'),
    end,
  end
  if area,
    if length(area(1,:)) ~= 2,
      error('Area data is not in the standard format'),
    end,
  end
  if gencost,
    if length(gencost(1,:)) ~= max(gencost(:,4)),
      error('Generator cost data is not in the standard format'),
    end,
  end
catch
  fm_disp(['Error in data file ',pathname,filename])
  fm_disp(lasterr)
  fm_disp('Conversion process interrupted.',2)
  return
end

% definition of file name for PSAT data file
extension = findstr(filename,'.');
newfile = ['d_',filename(1:extension(end)-1),'.m'];
% open *.m file for writing data
fid = fopen([pathname,newfile], 'wt');
if fid == -1,
  fm_disp(['Can''t open file ',pathname,newfile],2),
  return,
end

% Bus data: Bus.con
% ----------------------------------------------------------------------

nbus = length(bus(:,1));
count = fprintf(fid,['%% ',datestr(now,2), ...
                    ' File data originated from Matpower data file\n']);
count = fprintf(fid,['%% \n\n']);
count = fprintf(fid, 'Bus.con = [ ...\n');
format = '   %4d %8.4g %8.4g %8.4g %4d %4d;\n';
count = fprintf(fid,format,bus(:,[1,10,8,9,7,11])');
count = fprintf(fid,'   ];\n\n');

% Swing Generator data: SW.con
% ----------------------------------------------------------------------

k = find(bus(:,2) == 3);
slackbus = k;
if ~isempty(k)
  n = length(k);
  h = [];
  for i = 1:n, h = [h; find(gen(:,1) == bus(k(i),1))]; end
  if length(h) ~= n,
    fm_disp('Slack bus not found. Conversion process interrupted.',2)
    return
  end
  swline = [bus(k,1),gen(h,7),bus(k,10),gen(h,6),bus(k,9),gen(h,[4,5]), ...
            bus(k,[12,13]),gen(h,2),(gen(h,8) > 0)];
  count = fprintf(fid, 'SW.con = [ ...\n');
  format = ['   %4d',repmat(' %8.4g',1,9),' 1 1 %2u;\n'];
  count = fprintf(fid,format,swline');
  count = fprintf(fid,'   ];\n\n');
end

% PV Generator data: PV.con
% ----------------------------------------------------------------------

%k = find(bus(:,2) == 2);
if ~isempty(gen)
  n = length(gen(:,1));
  h = [];
  for i = 1:n
    h = [h; find(bus(:,1) == gen(i,1))];
  end
  pvline = [bus(h,1),gen(:,7),bus(h,10),gen(:,[2,6,4,5]), ...
            bus(h,[12,13]),(gen(:,8) > 0)];
  k = find(h == slackbus);
  if k, pvline(k,:) = []; end
  count = fprintf(fid, 'PV.con = [ ...\n');
  format = ['%4d',repmat(' %8.4g',1,8),' 1 %2u;\n'];
  count = fprintf(fid,format,pvline');
  count = fprintf(fid,'   ];\n\n');
end

% Constant Power Load data: PQ.con
% ----------------------------------------------------------------------

k1 = find(bus(:,3) ~= 0 | bus(:,4) ~= 0);
% for buses without generators or loads, the Vmin/Vmax must be kept.
k2 = find(bus(:,3) == 0 & bus(:,4) == 0  & bus(:,2) == 1);
count = fprintf(fid, 'PQ.con = [ ...\n');
format = ['%4d ',num2str(baseMVA),repmat(' %8.4g',1,5),' 0 1;\n'];
count = fprintf(fid,format,bus([k1;k2],[1,10,3,4,12,13])');
count = fprintf(fid,'   ];\n\n');

% Shunt Impedance data: Shunt.con
% ----------------------------------------------------------------------

k = find(bus(:,5) | bus(:,6));
if ~isempty(k)
  count = fprintf(fid, 'Shunt.con = [ ...\n');
  format = ['%4d ',num2str(baseMVA),' %8.4g 60 %8.4g %8.4g 1;\n'];
  count = fprintf(fid,format,bus(k,[1 10 5 6])');
  count = fprintf(fid,'   ];\n\n');
end

% Branch data: Line.con
% ----------------------------------------------------------------------

if ~isempty(branch)
  busmax = max(bus(:,1));
  busint = zeros(busmax,1);
  for i = 1:nbus
    busint(round(bus(i,1))) = i;
  end
  nline = length(branch(:,1));
  Line_con = zeros(nline,13);
  Line_con(:,3) = bus(busint(branch(:,1)),10);
  fr = busint(branch(:,1));
  to = busint(branch(:,2));
  Line_con(:,4) = abs(sign(branch(:,9))).*bus(fr,10)./bus(to,10);
  Line_con(:,[1 2 5 6 7 8 9 10 11 12 13]) = ...
      branch(:,[1 2 3 4 5 9 10 6 7 8 11]);
  count = fprintf(fid, 'Line.con = [ ...\n');
  format = ['%4d %4d ',num2str(baseMVA),' %8.4g 60 0 ', ...
            repmat(' %8.4g',1,9), ' %2u;\n'];
  count = fprintf(fid,format,Line_con');
  count = fprintf(fid,'   ];\n\n');
end

% Area data
% ----------------------------------------------------------------------
if ~isempty(area),
  fm_disp(['Area data are not defined in PSAT for ', ...
           'OPF computations.'])
end

% Supply data
% ----------------------------------------------------------------------
if ~isempty(gencost)
  startup = find(gencost(:,2));
  shutdown = find(gencost(:,3));
  if startup
    fm_disp(['Generation startup not supported yet. ', ...
             'Startup costs will be ignored.'])
  end
  if shutdown
    fm_disp(['Generation shutdown not supported yet. ', ...
             'Shutdown costs will be ignored.'])
  end
  ngen = length(gen(:,1));
  ncost = length(gencost(:,1));
  coeff = zeros(ncost,3);
  h = find(gencost(:,1) == 1);
  if h,
    fm_disp(['Piecewise linear generator costs are ', ...
             'converted in polynomial approximations.'])
  end
  for i = 1:length(h)
    n = gencost(h(i),4);
    xidx = 4+(1:2:2*n);
    yidx = 4+(2:2:2*n);
    a = polyfit(gencost(h(i),xidx),gencost(h(i),yidx),min(2,n-1));
    if length(a) == 3,
      a = a(3:-1:1);
    elseif length(a) == 2,
      a = [a(2), a(1), 0];
    else
      a = [0 1 0];
    end
    coeff(h(i),:) = a;
  end
  h = find(gencost(:,1) == 2);
  if h & ~isempty(find(gencost(h,4) > 3)),
    fm_disp('Polynomial generator costs are reduced to 2nd order polynomials.'),
  end
  for i = 1:length(h)
    n = gencost(h(i),4);
    a = gencost(h(i),5:min(4+n,8));
    if length(a) == 3,
      a = a(3:-1:1);
    elseif length(a) == 2,
      a = [a(2), a(1), 0];
    else
      a = [0 1 0];
    end
    coeff(h(i),:) = a;
  end
  Supply_con = zeros(ngen,13);
  Supply_con(:,[1 2 4 5]) = gen(:,[1 7 9 10]);
  Supply_con(:,[7 8 9]) = coeff(1:ngen,:);
  if ncost == 2*ngen
    Supply_con(:,[10 11 12]) = coeff(ngen+1:ncost,:);
  end
  Supply_con(:,13) = gen(:,8);
  count = fprintf(fid, 'Supply.con = [ ...\n');
  format = ['%4d ',repmat(' %8.4g',1,11),' %2d 0 1;\n'];
  count = fprintf(fid,format,Supply_con');
  count = fprintf(fid,'   ];\n\n');
end

% Bus Names
% ----------------------------------------------------------------------

bus_name = cell(nbus,1);
for i = 1:nbus,
  bus_name{i,1} = deblank(['Bus',fvar(bus(i,1),5)]);
end
count = fprintf(fid, 'Bus.names = {...\n      ');
for i = 1:nbus-1
  count = fprintf(fid, ['''', deblank(bus_name{i}),'''; ']);
  if rem(i,5) == 0
    count = fprintf(fid,'\n      ');
  end
end
count = fprintf(fid, ['''', deblank(bus_name{end}),'''};\n\n']);

% end of operations
fm_disp(['Conversion into data file "',pathname,newfile,'" completed.'])
if Settings.beep
  beep
end
fclose(fid);
check = 1;