Y Bus Assembly (Using Singular Transformation technique)

Y_Bus_Assembly_Singular.m

function yBus = Y_Bus_Assembly_Singular()
    % Data from IEEE 5-bus Data
    % "http://shodhganga.inflibnet.ac.in/bitstream/10603/26549/14/14_appendix.pdf""
    lineData = csvread('Y_Bus_Data.csv', 1, 0);
    % Calculate the number of lines in the line data
    lines = length(lineData);
    % Calculate the number of buses in the line data
    buses = max(max(lineData(:,2)),max(lineData(:,3)));
    % Calculate the line admittance by taking element by element inverse
    lineAdmittance = lineData(:,4).^-1;
    % Transform the lineAdmittance array into a diagonal matrix
    y = diag(lineAdmittance);
    % Initialize the admittance matrix with dimensions (lines X buses)
    admittanceMatrix = zeros(lines,buses);

    % Fill the Admittance Matrix by reading the lineData
    for k = 1:lines
        sendingBus = lineData(k,2);
        receivingBus = lineData(k,3);
        % Fill in the sending nodes
        if sendingBus ~= 0               % condition used to remove reference node
            admittanceMatrix(k, sendingBus) = 1;
        end
        % Fill in the receiving nodes
        if receivingBus ~= 0              % condition used to remove reference node
            admittanceMatrix(k, receivingBus) = -1;
        end
    end

    % Calculate the Y Bus
    yBus = admittanceMatrix'*y*admittanceMatrix;
end

Y_Bus_Data.csv