calculo_Yb.m

%  <CÁLCULO DO FLUXO DE POTÊNCIA - LOAD FLOW APPLYING NEWTON-RAPHSON METHOD V1.0. 
%  This is the main source of this software that calculates the power flow of a power network described using an excel input data file >
%     Copyright (C) <2014>  <Sebastián de Jesús Manrique Machado>   <e-mail:sebajmanrique747@gmail.com>
% 
%     This program is free software: you can redistribute it and/or modify
%     it under the terms of the GNU General Public License as published by
%     the Free Software Foundation, either version 3 of the License, or
%     (at your option) any later version.
% 
%     This program is distributed in the hope that it will be useful,
%     but WITHOUT ANY WARRANTY; without even the implied warranty of
%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%     GNU General Public License for more details.
% 
%     You should have received a copy of the GNU General Public License
%     along with this program.  If not, see <http://www.gnu.org/licenses/>.

%Cálculo_Fluxo_de_Potência_ Newton Raphson
%   Sebastián de Jesús Manrique Machado
%   Estudante_Mestrado Em Engenharia Elétrica
%   UEL - 2014.
%Função para achar Y barras.

function [Y_barras, b_shunt_m] = calculo_Yb(Ys_linha, b_shunt, num_barras, num_linhas, nodo_i, nodo_j)

Y_barras = zeros(num_barras,num_barras);        %O tamanho desta matriz é igual ao número de barras
b_shunt_m = zeros(num_barras,num_barras);

for i = 1 : num_linhas
    % Elementos fora da diagonal
    Y_barras( nodo_i(i), nodo_j(i) ) = Y_barras( nodo_i(i), nodo_j(i) ) - Ys_linha(i);
    Y_barras( nodo_j(i), nodo_i(i) ) = Y_barras( nodo_i(i), nodo_j(i) );
    
    b_shunt_m( nodo_i(i), nodo_j(i) ) = b_shunt_m( nodo_i(i), nodo_j(i) ) + imag(b_shunt(i));
    b_shunt_m( nodo_j(i), nodo_i(i) ) = b_shunt_m( nodo_i(i), nodo_j(i) );
    
    % Elementos da diagonal
    Y_barras( nodo_i(i), nodo_i(i) ) = Y_barras( nodo_i(i), nodo_i(i) ) + Ys_linha(i) + b_shunt(i);
    Y_barras( nodo_j(i), nodo_j(i) ) = Y_barras( nodo_j(i), nodo_j(i) ) + Ys_linha(i) + b_shunt(i);
end



end