Example_9bus test reduction code on 9-bus system (case9) Created by Yujia Zhu, yzhu54@asu.edu, Oct. 2014
0001 % Example_9bus 0002 % test reduction code on 9-bus system (case9) 0003 % 0004 % Created by Yujia Zhu, yzhu54@asu.edu, Oct. 2014 0005 0006 % MATPOWER 0007 % Copyright (c) 2014-2015 by Power System Engineering Research Center (PSERC) 0008 % by Yujia Zhu, PSERC ASU 0009 % 0010 % $Id: Example_9bus.m 2655 2015-03-18 16:40:32Z ray $ 0011 % 0012 % This file is part of MATPOWER. 0013 % Covered by the 3-clause BSD License (see LICENSE file for details). 0014 % See http://www.pserc.cornell.edu/matpower/ for more info. 0015 0016 %% load 9-bus case in MATPOWER format 0017 % If you have installed MATPOWER, you can use function "loadcase" instead 0018 % to read the case data 0019 mpc = loadcase('case9'); % input full model in MATPOWER case format 0020 %% give external bus indices 0021 ExBus=[1,5,8]'; % input list of external buses (buses to be eliminated) 0022 %% run reduction subroutine (MPReduction) 0023 %Input notes: 0024 % mpc: struct, the original full model in MATPOWER case format 0025 % ExBus: 1*n array, list of external buses 0026 [mpcreduced,Link,BCIRCr]=MPReduction(mpc,ExBus,0); % call reduction subroutine 0027 % Output notes: 0028 % 1. The output mpcreduced is the reduced 6 bus model 0029 % 2. All branch B shunts are converted to B shunts on buses 0030 % 3. There are 4 equivalent branches generated in the reduction process 0031 % branch between bus: 2-7, 2-9, 4-6, 7-9, all equivalent branches have 0032 % circuit number 99. 0033 % 4. Equivalent lines have no line ratings. 0034 % 5. Link gives generator bus mapping showing how generators are moved. 0035 % E.g: Generator on bus 1 is moved to bus 4, all other generators are 0036 % not moved since they were on retained buses 0037 % 6. The reduction process will generate equivalent lines with large 0038 % impedance. The software eliminate equivalent lines whose impedance is 0039 % larger than the threshold value. The threshold value in the software is 0040 % 10 times of maximum impedance (reactance) value in the original full 0041 % model.