fm_matpower.m

电力系统的psat

function check = fm_matpower(filename, pathname, flag)
% FM_MATPOWER Conversion from Matpower Data format
%       into PSAT data format
%
% CHECK = FM_MATPOWER(FILENAME,PATHNAME,FLAG)
%       FILENAME name of the file to be converted
%       PATHNAME path of the file to be converted
%       FLAG = 1 append original file as a comment
%       FLAG = 0 do not append original file
%
%       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:  www.power.uwaterloo.ca
%
%*** Contributions ***
%
%Author:    Liulin
%Date:      13-Mar-2003
%Update:    01-May-2004
%E-mail:    [e-mail]
%
% Copyright (C) 2002-2005 Federico Milano
%
% This toolbox is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2.0 of the License, or
% (at your option) any later version.
%
% This toolbox is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANDABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
% General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this toolbox; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307,
% USA.

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

% load original file into cell string
if flag == 1,
  orig_file = textread([pathname,filename],'%s', ...
                       'delimiter','\n','whitespace','');
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
% ----------------------------------------------------------------------

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');
nbus = length(bus(:,1));
for i = 1:nbus-1
  count = fprintf(fid,'%4d %8.4g %8.4g %8.4g %4d %4d; ', ...
                  bus(i,[1,10,8,9,7,11]));
  if rem(i,3) == 0; count = fprintf(fid,'\n'); end
end
count = fprintf(fid, '%4d %8.4g %8.4g %8.4g %4d %4d];\n\n\n', ...
                bus(end,[1,10,8,9,7,11]));

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

k = find(bus(:,2) == 3);
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),zeros(n,1)];
  count = fprintf(fid, 'SW.con = [ ...\n');
  format = ['%4d',repmat(' %8.4g',1,9),' %2u'];
  count = fprintf(fid,[format,';\n'],swline(1:end-1,:)');
  count = fprintf(fid,[format,'];\n\n\n'],swline(end,:));
end

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

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

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

k1 = find(bus(:,3) ~= 0 | bus(:,4) ~= 0);
if ~isempty(k1)
  n = length(k1);
  pqline1 = [bus(k1,1),baseMVA*ones(n,1),bus(k1,[10,3,4,12,13]),zeros(n,1)];
end
% for buses without generators or loads, the Vmin/Vmax must be kept.
k2 = find(bus(:,3) == 0 & bus(:,4) == 0  & bus(:,2) == 1);
if isempty(k2)
  pqline2 = [];
else
  n = length(k2);
  pqline2 = [bus(k2,1),baseMVA*ones(n,1),bus(k2,[10,3,4,12,13]),zeros(n,1)];
end
pqline = [pqline1; pqline2];
n = size(pqline,1);
count = fprintf(fid, 'PQ.con = [ ...\n');
format = ['%4d',repmat(' %8.4g',1,6),' %2u'];
count = fprintf(fid,[format,';\n'],pqline(1:end-1,:)');
count = fprintf(fid,[format,'];\n\n\n'],pqline(end,:));

%k = find(bus(:,3) ~= 0 | bus(:,4) ~= 0);
%if ~isempty(k)
%  n = length(k);
%  pqline = [bus(k,1),baseMVA*ones(n,1),bus(k,[10,3,4,12,13]),zeros(n,1)];
%  count = fprintf(fid, 'PQ.con = [ ...\n');
%  format = ['%4d',repmat(' %8.4g',1,6),' %2u'];
%  count = fprintf(fid,[format,';\n'],pqline(1:end-1,:)');
%  count = fprintf(fid,[format,'];\n\n\n'],pqline(end,:));
%end

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

k = find(bus(:,5) | bus(:,6));
if ~isempty(k)
  count = fprintf(fid, 'Shunt.con = [ ...\n');
  n = length(k);
  format = ['%4d',repmat(' %8.4g',1,5)];
  shuntline = [bus(k,1),baseMVA*ones(n,1),bus(k,10),60*ones(n,1),bus(k,[5 6])];
  count = fprintf(fid,[format,';\n'],shuntline(1:end-1,:)');
  count = fprintf(fid,[format,'];\n\n\n'],shuntline(end,:));
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,15);
  Line_con(:,3) = baseMVA*ones(nline,1);
  Line_con(:,4) = bus(busint(branch(:,1)),10);
  Line_con(:,5) = 60*ones(nline,1);
  Line_con(:,7) = abs(sign(branch(:,9))).*bus(busint(branch(:,1)),10)./bus(busint(branch(:,2)),10);
  Line_con(:,[1 2 8 9 10 11 12 13 14 15]) = branch(:,[1 2 3 4 5 9 10 6 7 8]);
  count = fprintf(fid, 'Line.con = [ ...\n');
  format = ['%4d %4d',repmat(' %8.4g',1,13)];
  count = fprintf(fid,[format,';\n'],Line_con(1:end-1,:)');
  count = fprintf(fid,[format,'];\n\n\n'],Line_con(end,:));
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'];
  count = fprintf(fid,[format,';\n'],Supply_con(1:end-1,:)');
  count = fprintf(fid,[format,'];\n\n\n'],Supply_con(end,:));
end

% Bus Names: Varname.bus
% ----------------------------------------------------------------------

bus_name = cell(nbus,1);
for i = 1:nbus,
  bus_name{i,1} = deblank(['Bus',fvar(bus(i,1),5)]);
end
count = fprintf(fid, 'Varname.bus = {...\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']);

% append original data file
if flag == 1
  count = fprintf(fid, '\n\n%% ORIGINAL DATA FILE: \n\n');
  for i = 1:length(orig_file)
    count = fprintf(fid,'%% %s \n\n',orig_file{i,1});
  end
end

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