Power Systems Loss Coefficients (by Arango, Urrego, Rivera)

code to compute the B-coefficients (loss coefficiens using matpower)

AUTHORS: Dario Arango, Ricardo Urrego, Sergio Rivera

% Kron Coefficents (loss coefficients) for power systems using matpower (please note that matpower must be in you path)
% Code developed to calculate the kron coefficents (in order to get later robust
% loss coefficients) mentioned in the following paper
% 'Robust Loss Coefficients: Application to Power Systems with Solar and Wind Energy'
% by Dario Arango, Ricardo Urrego, Sergio Rivera ([e-mail])
% in International Journal of Power and Energy Conversion
%
% [B B0 B00 PL PLV]=kron(caso)
%
% caso must have consecutive bus numbering (i.e. internal bus numbering)
%
% LOSSES = (Pgen'*(B)*Pgen +B0*Pgen+B00)*Sb;
%
% Example:
% [B B0 B00 PL PLV]=KronCoeff(case30)
%
% results=runpf('case30');
% Sb=results.baseMVA;
% Pgen=results.gen(:,2)/Sb;
%
% LOSSES_KronCoefficients = (Pgen'*(B)*Pgen +B0*Pgen+B00)*Sb % It
% %must be equal to PL
%
% LOSSES_PowerFlow=sum(results.branch(:,14)+results.branch(:,16))
% %It mus be equal to PLV
%
% If you use the code please you can cite the paper in the following way:
% Arango, D.; Urrego, R.; Rivera, S. (in press). Robust Loss Coefficients:
% Application to Power Systems with Solar and Wind Energy. International
% Journal of Power and Energy Conversion. 2017