www.gusucode.com > 基于lingo求所以解，对潮流计算求出所有解 > matpower4.1/cdf2matp.m
function cdf2matp(cdf_file, matp_file, mpcver)
%CDF2MATP Converts data from IEEE Common Data Format to MATPOWER format.
% CDF2MATP : prompts for input & output file names
% CDF2MATP(CDF_FILE, MATP_FILE) : uses specified file names
% CDF2MATP(CDF_FILE, MATP_FILE, v) : and MATPOWER case version
% v = '1' or '2', default is '2'
%
% Optional arguments are the names of the input IEEE CDF file and the
% output MATPOWER case file. If these are not given you will be prompted
% to enter them.
%
% The IEEE CDF does not include some data need to run an optimal power
% flow. This script creates default values for some of this data as
% follows:
%
% Bus data:
% Vmin = 0.94 p.u.
% Vmax = 1.06 p.u.
% Gen data:
% Pmin = 0 MW
% Pmax = Pg + baseMVA
% Gen cost data:
% Quadratic costs with:
% c2 = 10 / Pg, c1 = 20, c0 = 0, if Pg is non-zero, and
% c2 = 0.01, c1 = 40, c0 = 0, if Pg is zero
% This should yield an OPF solution "close" to the
% existing solution (assuming it is a solved case)
% with lambdas near $40/MWh. See 'help caseformat'
% for details on the cost curve format.
%
% CDF2MATP may modify some of the data which are "infeasible" for
% running optimal power flow. If so, warning information will be
% printed out on screen.
%
% Note: Since our code can not handle transformers with variable tap,
% you may not expect to get exactly the same power flow solution
% using converted data. This is the case when we converted ieee300.cdf.
% MATPOWER
% $Id: cdf2matp.m,v 1.25 2011/12/01 19:06:47 cvs Exp $
% by Deqiang (David) Gan, PSERC Cornell & Zhejiang University
% Copyright (c) 1996-2010 by Power System Engineering Research Center (PSERC)
%
% This file is part of MATPOWER.
% See http://www.pserc.cornell.edu/matpower/ for more info.
%
% MATPOWER 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 3 of the License,
% or (at your option) any later version.
%
% MATPOWER is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY 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 MATPOWER. If not, see <http://www.gnu.org/licenses/>.
%
% Additional permission under GNU GPL version 3 section 7
%
% If you modify MATPOWER, or any covered work, to interface with
% other modules (such as MATLAB code and MEX-files) available in a
% MATLAB(R) or comparable environment containing parts covered
% under other licensing terms, the licensors of MATPOWER grant
% you additional permission to convey the resulting work.
%% define named indices into bus, gen, branch matrices
[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
[PW_LINEAR, POLYNOMIAL, MODEL, STARTUP, SHUTDOWN, NCOST, COST] = idx_cost;
rev = '$Revision: 1.25 $';
% ----- get the original load flow data -----
if nargin < 1
cdf_file = input('Please enter IEEE Common Data Format file name : ', 's');
end
if nargin < 2
matp_file = input('Please enter MATPOWER file name : ', 's');
end
if nargin < 3
mpcver = '2';
end
%% verify valid input filename
[cdf_path cdf_name cdf_ext] = fileparts(cdf_file);
if isempty(cdf_ext)
cdf_ext = '.cdf';
cdf_file = strcat(cdf_name , cdf_ext);
end
%% open input file
[fid, msg] = fopen(cdf_file, 'r');
if fid < 0
disp(msg);
error('cdf2matp: Can not read the input file: %s', cdf_file )
end
%% verify valid output filename
if isempty(matp_file)
matp_file = strcat(cdf_name, '.m');
rootname = cdf_name;
else
[matp_path matp_name matp_ext] = fileparts(matp_file);
rootname = matp_name;
matp_file = strcat(matp_name, matp_ext);
end
% set up some comments
comments = { sprintf('%s', upper(rootname)), ...
' Please see ''help caseformat'' for details on the case file format.', ...
' This data was converted from IEEE Common Data Format', ...
sprintf(' (%s) on %s by cdf2matp, rev. %s', cdf_file, date, rev(12:end-2)), ...
' See end of file for warnings generated during conversion.', ...
'' };
warnings = {};
% get baseMVA
title_cdf = fgetl(fid);
if isnumeric(str2num(title_cdf(32:37))) && not(isempty(str2num(title_cdf(32:37))))
baseMVA = str2num(title_cdf(32:37));
if length(findstr(title_cdf(2:9), '/')) == 2 %% date in the file
warnings{end+1} = sprintf('***** check the title format in the first line of the cdf file.');
end
else
error('cdf2matp:argChk','Error getting the Base MVA, check the title format in the first line of the file.')
end
% adding the cdf title
comments{end+1} = title_cdf;
% find string 'BUS DATA FOLLOWS'
while 1
line = fgetl(fid);
if line(1:16) == 'BUS DATA FOLLOWS', break, end
end
% ----- get bus data, feed them into matrix bus, gen, gencost
ibus = 0;
igen = 0;
iarea = 0;
while 1
line = fgetl(fid);
if line(1:4) == '-999', break, end
% feed bus data
ibus = ibus + 1;
bus(ibus, BUS_I) = str2num(line(1:4)); % bus number
busnames(ibus,:) = line(6:17); % bus names
bus(ibus, BUS_TYPE) = str2num(line(25:26));
if bus(ibus, BUS_TYPE) == 0 % bus type
bus(ibus, BUS_TYPE) = 1;
end
if (bus(ibus, BUS_TYPE) < 2) % Pd
bus(ibus, PD) = str2num(line(41:49)) - str2num(line(60:67));
elseif (bus(ibus, BUS_TYPE) >= 2)
bus(ibus, PD) = str2num(line(41:49));
end
bus(ibus, QD) = str2num(line(50:59)); % Qd
bus(ibus, GS) = baseMVA*str2num(line(107:114)); % Gs
bus(ibus, BS) = baseMVA*str2num(line(115:122)); % Bs
bus(ibus, BUS_AREA) = str2num(line(19:20)); % area
bus(ibus, VM) = str2num(line(28:33)); % Vm
bus(ibus, VA) = str2num(line(34:40)); % Va
bus(ibus, BASE_KV) = str2num(line(77:83)); % baseKV
bus(ibus, ZONE) = str2num(line(21:23)); % zone
bus(ibus, VMAX) = 1.06; % default voltage upper limit
bus(ibus, VMIN) = 0.94; % default voltage lower limit
% feed gen and gencost
Pg = str2num(line(60:67));
Qg = str2num(line(68:75));
Qmax = str2num(line(91:98));
Qmin = str2num(line(99:106));
if bus(ibus, BUS_TYPE) >= 2
igen = igen + 1;
if bus(ibus, BUS_TYPE) == 3, refgen = igen; end
gen(igen, GEN_BUS) = bus(ibus, BUS_I); % bus number
gen(igen, PG) = Pg; % Pg
if gen(igen, PG) < 0 % negative Pg is transformed as load
bus(ibus, PD) = bus(ibus, PD) - gen(igen, PG);
warnings{end+1} = sprintf('***** negative Pg at bus %g treated as Pd', bus(ibus, BUS_I));
fprintf('\n %s', warnings{end});
gen(igen, PG) = 0;
end
gen(igen, QG) = Qg; % Qg
gen(igen, QMAX) = Qmax; % Qmax
gen(igen, QMIN) = Qmin; % Qmin
if Qmax - Qmin < 0.01 % Qmax is modified
gen(igen, QMAX) = Qmin + 0.1 * baseMVA;
warnings{end+1} = sprintf('***** Qmax = Qmin at generator at bus %4i (Qmax set to Qmin + %g)', bus(ibus, BUS_I), baseMVA/10);
fprintf('\n %s', warnings{end});
end
gen(igen, VG) = str2num(line(85:90)); % specified voltage
gen(igen, MBASE) = baseMVA; % baseMVA
gen(igen, GEN_STATUS) = 1; % default status is 'on'
gen(igen, PMAX) = gen(igen, 2) + baseMVA; % Pmax
gen(igen, PMIN) = 0; % Pmin = 0 by default
gencost(igen, MODEL) = POLYNOMIAL; % by default, sets the model as polynomial
gencost(igen, STARTUP) = 0; % start up cost is zero by default
gencost(igen, SHUTDOWN) = 0; % shut down cost is zero by default
gencost(igen, NCOST) = 3; % number of coefficients in polynomial cost
% gencost(igen, COST) = 0.01; % default c2
% gencost(igen, COST+1) = 0.3; % default c1
% gencost(igen, COST+2) = 0.2; % default c0
end
end
totload = sum(bus(:, PD));
totgen = sum(gen(:, PG));
if totgen < 1.04 * totload
gen(refgen, PMAX) = gen(refgen, PG) + 1.1 * totload - totgen; % Pg at slack bus is modified
warnings{end+1} = sprintf('***** Insufficient generation, setting Pmax at slack bus (bus %d) to %g', gen(refgen, [GEN_BUS, PMAX]));
fprintf('\n %s', warnings{end});
end
% ----- set up the cost coefficients of generators
ng = size(gen, 1);
% gencost(:, COST) = zeros(ng, 1);
% gencost(:, COST+1) = 100*ones(ng, 1) ./ (gen(:, PG) + 10*ones(ng, 1));
% gencost(:, COST+2) = 100*ones(ng, 1) ./ (gen(:, PG) + 10*ones(ng, 1));
zg = find(gen(:, PG) == 0); %% for Pg = 0
gencost(zg, COST) = 0.01 * ones(size(zg));
gencost(zg, COST+1) = 40 * ones(size(zg));
nzg = find(gen(:, PG) ~= 0); %% Pg non-zero
gencost(nzg, COST) = 10 * ones(size(nzg)) ./ gen(nzg, PG);
gencost(nzg, COST+1) = 20 * ones(size(nzg));
gencost(:, COST+2) = zeros(ng, 1);
% find string 'BRANCH DATA FOLLOWS'
while 1
line = fgetl(fid);
if line(1:19) == 'BRANCH DATA FOLLOWS', break, end
end
% ----- get branch data, feed them into matrix branch
k = 0;
while 1
line = fgetl(fid);
if line(1:4) == '-999', break, end
k = k + 1;
branch(k, F_BUS) = str2num(line(1:4)); % fbus (also the tap bus)
branch(k, T_BUS) = str2num(line(6:9)); % tbus
branch(k, BR_R) = str2num(line(20:29)); % R
branch(k, BR_X) = str2num(line(30:40)); % X
branch(k, BR_B) = str2num(line(41:50)); % B
branch(k, RATE_A) = str2num(line(51:55)); % RATE A
if branch(k, RATE_A) < 0.000001
branch(k, RATE_A) = 99 * baseMVA; % RATE A is modified
warnings{end+1} = sprintf('***** MVA limit of branch %d - %d not given, set to %g', branch(k, [1, 2, 6]));
fprintf('\n %s', warnings{end});
end
branch(k, RATE_B) = str2num(line(57:61)); % RATE B
branch(k, RATE_C) = str2num(line(63:67)); % RATE C
branch(k, TAP) = str2num(line(77:82)); % transformer turns ratio
branch(k, SHIFT) = 0; % phase shifter can not be modelled
branch(k, BR_STATUS) = 1; % by default, branch is on
end
fprintf('\n');
fclose(fid);
% put in struct
mpc.baseMVA = baseMVA;
mpc.bus = bus;
mpc.branch = branch;
mpc.gen = gen;
mpc.gencost = gencost;
mpc = loadcase(mpc); %% convert to internal (e.g. v. '2') case format
% ----- write data in MATPOWER format -----
matp_file = savecase(matp_file, comments, mpc, mpcver);
% print the name buses
busnames = strtrim(busnames); %% remove leading and trailing blanks
if length(matp_file) < 4 || ~strcmp(matp_file(end-3:end), '.mat')
fid = fopen(matp_file, 'at'); %% open for append
fprintf(fid, '\n%% bus names\nmpc.busnames = [\n');
for i = 1:size(busnames, 1)
fprintf(fid, '\t''%s''\n', busnames(i, :));
end
fprintf(fid, '];\n');
fclose(fid);
end
% print conversion warnings as comments
if ~isempty(warnings)
if length(matp_file) < 4 || ~strcmp(matp_file(end-3:end), '.mat')
fid = fopen(matp_file, 'at'); % open for append
fprintf(fid, '\n%% Warnings from cdf2matp conversion:\n%%\n');
for i=1:length(warnings)
fprintf(fid, '%% %s\n', warnings{i});
end
fclose(fid);
end
end
% to do:
% PTI format can not be handled