📄 test30_1.m
字号:
function [baseMVA, bus, gen, branch, areas, gencost] = test30_1(i)%CASE30 Power flow data for 30 bus, 6 generator case.% Please see 'help caseformat' for details on the case file format.%% Based on data from ...% Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",% IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,% 1974, pp. 745-751.% ... with branch parameters rounded to nearest 0.01, shunt values divided% by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out.% Generator locations, costs and limits and bus areas were taken from ...% Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis% in deregulated power systems using game theory", IEEE Transactions on% Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347.% Generator Q limits were derived from Alsac & Stott, using their Pmax% capacities. V limits and line |S| limits taken from Alsac & Stott.% MATPOWER% $Id: case30.m,v 1.7 2004/09/21 02:41:49 ray Exp $%%----- Power Flow Data -----%%%% system MVA basebaseMVA = 100;%% bus data% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vminbus =[ 28.0000 3.0000 0 0 0 0 1.0000 1.0000 0 135.0000 0 1.0500 0.9500 29.0000 2.0000 21.7000 12.7000 0 0 1.0000 1.0000 0 135.0000 0 1.1000 0.9500 25.0000 1.0000 2.4000 1.2000 0 0 1.0000 1.0000 0 135.0000 0 1.0500 0.9500 23.0000 1.0000 7.6000 1.6000 0 0 1.0000 1.0000 0 135.0000 0 1.0500 0.9500 30.0000 1.0000 0 0 0 0.1900 1.0000 1.0000 0 135.0000 0 1.0500 0.9500 2.0000 1.0000 0 0 0 0 1.0000 1.0000 0 135.0000 2.0000 1.0500 0.9500 26.0000 1.0000 22.8000 10.9000 0 0 1.0000 1.0000 0 135.0000 0 1.0500 0.9500 6.0000 1.0000 30.0000 30.0000 0 0 1.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 5.0000 1.0000 0 0 0 0 1.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 14.0000 1.0000 5.8000 2.0000 0 0 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 15.0000 1.0000 0 0 0 0 1.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 22.0000 1.0000 11.2000 7.5000 0 0 2.0000 1.0000 0 135.0000 0 1.0500 0.9500 27.0000 2.0000 0 0 0 0 2.0000 1.0000 0 135.0000 0 1.1000 0.9500 21.0000 1.0000 6.2000 1.6000 0 0 2.0000 1.0000 0 135.0000 0 1.0500 0.9500 1.0000 1.0000 8.2000 2.5000 0 0 2.0000 1.0000 0 135.0000 1.0000 1.0500 0.9500 24.0000 1.0000 3.5000 1.8000 0 0 2.0000 1.0000 0 135.0000 0 1.0500 0.9500 17.0000 1.0000 9.0000 5.8000 0 0 2.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 3.0000 1.0000 3.2000 0.9000 0 0 2.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 7.0000 1.0000 9.5000 3.4000 0 0 2.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 16.0000 1.0000 2.2000 0.7000 0 0 2.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 18.0000 1.0000 17.5000 11.2000 0 0 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 19.0000 2.0000 0 0 0 0 3.0000 1.0000 0 135.0000 3.0000 1.1000 0.9500 4.0000 2.0000 3.2000 1.6000 0 0 2.0000 1.0000 0 135.0000 3.0000 1.1000 0.9500 8.0000 1.0000 8.7000 6.7000 0 0.0400 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 9.0000 1.0000 0 0 0 0 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 10.0000 1.0000 3.5000 2.3000 0 0 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 11.0000 2.0000 0 0 0 0 3.0000 1.0000 0 135.0000 3.0000 1.1000 0.9500 20.0000 1.0000 0 0 0 0 1.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 12.0000 1.0000 2.4000 0.9000 0 0 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500 13.0000 1.0000 10.6000 1.9000 0 0 3.0000 1.0000 0 135.0000 3.0000 1.0500 0.9500];%% generator data% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmingen = [ 28 23.54 0 150 -20 1 100 1 80 0; 29 60.97 0 60 -20 1 100 1 80 0; 19 21.59 0 62.5 -15 1 100 1 50 0; 11 26.91 0 48.7 -15 1 100 1 55 0; 4 19.2 0 40 -10 1 100 1 30 0; 27 37 0 44.7 -15 1 100 1 40 0;];%% branch data% fbus tbus r x b rateA rateB rateC ratio angle statusbranch = [ 28 29 0.02 0.06 0.03 130 130 130 0 0 1; 28 25 0.05 0.19 0.02 130 130 130 0 0 1; 29 23 0.06 0.17 0.02 65 65 65 0 0 1; 25 23 0.01 0.04 0 130 130 130 0 0 1; 29 30 0.05 0.2 0.02 130 130 130 0 0 1; 29 2 0.06 0.18 0.02 65 65 65 0 0 1; 23 2 0.01 0.04 0 90 90 90 0 0 1; 30 26 0.05 0.12 0.01 70 70 70 0 0 1; 2 26 0.03 0.08 0.01 130 130 130 0 0 1; 2 6 0.01 0.04 0 32 32 32 0 0 1; 2 5 0 0.21 0 65 65 65 0 0 1; 2 14 0 0.56 0 32 32 32 0 0 1; 5 15 0 0.21 0 65 65 65 0 0 1; 5 14 0 0.11 0 65 65 65 0 0 1; 23 22 0 0.26 0 65 65 65 0 0 1; 22 27 0 0.14 0 65 65 65 0 0 1; 22 21 0.12 0.26 0 32 32 32 0 0 1; 22 1 0.07 0.13 0 32 32 32 0 0 1; 22 24 0.09 0.2 0 32 32 32 0 0 1; 21 1 0.22 0.2 0 16 16 16 0 0 1;% 24 17 0.08 0.19 0 16 16 16 0 0 1; 1 3 0.11 0.22 0 16 16 16 0 0 1; 3 7 0.06 0.13 0 16 16 16 0 0 1;% 19 16 0.03 0.07 0 32 32 32 0 0 1; 14 16 0.09 0.21 0 32 32 32 0 0 1; 14 17 0.03 0.08 0 32 32 32 0 0 1; 14 18 0.03 0.07 0 32 32 32 0 0 1; 14 19 0.07 0.15 0 32 32 32 0 0 1; 18 19 0.01 0.02 0 32 32 32 0 0 1; 1 4 0.1 0.2 0 16 16 16 0 0 1;% 19 8 0.12 0.18 0 16 16 16 0 0 1; 4 8 0.13 0.27 0 16 16 16 0 0 1; 8 9 0.19 0.33 0 16 16 16 0 0 1; 9 10 0.25 0.38 0 16 16 16 0 0 1; 9 11 0.11 0.21 0 16 16 16 0 0 1;% 20 11 0 0.4 0 65 65 65 0 0 1; 11 12 0.22 0.42 0 16 16 16 0 0 1; 11 13 0.32 0.6 0 16 16 16 0 0 1; 12 13 0.24 0.45 0 16 16 16 0 0 1; 6 20 0.06 0.2 0.02 32 32 32 0 0 1; 2 20 0.02 0.06 0.01 32 32 32 0 0 1;];%%----- OPF Data -----%%%% area dataareas = [ 1 8; 2 23; 3 26;];%% generator cost data% 1 startup shutdown n x0 y0 ... xn yn% 2 startup shutdown n c(n-1) ... c0gencost = [ 2 0 0 3 0.02 2 0; 2 0 0 3 0.0175 1.75 0; 2 0 0 3 0.0625 1 0; 2 0 0 3 0.00834 3.25 0; 2 0 0 3 0.025 3 0; 2 0 0 3 0.025 3 0;];return;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -