We provide a MATLAB code for the analytical computation of nodal voltage sensitivity coefficients based on the proposed method in [1]. The numerical simulation presented in [1] can be ran by running the main.m (see Files for details). Other subfunctions are available and can work independently (see Files for details).
- [main.m] function that runs the simulation presented in [1]
- [SC_Voltage.m] function that computes complex, magnitude and angle nodal voltage sensitivity coefficients using the method presented in [1]
- [3ph Load Flow - IEEE34/Ymatrix_3ph.m] function that computes the compound admittance matrix. It takes as an input an excel file with a specific format.
- [3ph Load Flow - IEEE34/IEEE34Feeder.xlsx] Excel file formatted to fit the function Ymatrix_3ph.m. It contains the line parameters and base values values of the IEEE34 feeder [2] without the presence of the voltage regulators.
- [3ph Load Flow - IEEE34/NR_VoltageDepJacob.m] function that performs the Newton-Raphson (N-R) algorithm and outputs the network state (i.e. phase-to-ground nodal voltages at each node). It supports multiphases automitically depending on the inputted compound admittance matrix.
- [3ph Load Flow - IEEE34/NR_VoltageDepJacob_iteration.m] function that performs a sub-iteration of the Netwon-Raphson algorithm where the Jacobian of the k-th iteration is computed and the network state is updated using the inverted Jacobian and power mismatches of the k-th iteration.
- [LoadFlow input examples/E_star_IEEE34_example_noPQnodes.mat] setpoints for the magnitudes of nodal voltages for pv nodes for a simulation with no pv nodes. Input needed for N-R algorithm.
- [LoadFlow input examples/E_star_IEEE34_example_withPQnodes.mat] setpoints for the magnitudes of nodal voltages for pv nodes for a simulation with pv nodes. Input needed for N-R algorithm.
- [LoadFlow input examples/S_star_IEEE34_example.mat] setpoints for the apparent nodal power injections for pq nodes. Input needed for N-R algorithm.
- Running the function main.m should directly output all the results of the work in [1].
- All throughout the main.m function there are configurable inputs. These configurations can be set in the code line following the commented line: % (!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!) INPUT (!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!)
The following software is required to run the model:
- MATLAB ( in principle any version however tested on MATLAB version 2020a Update 1 (9.8.0.1359463) )
- [1] reference will be published soon. For more info contact the author at [email protected]
- [2] Kersting, W. H. (1991). Radial distribution test feeders. IEEE Transactions on Power Systems, 6(3), 975-985.